Dynamics of Coupled Electron–Boson Systems with the Multiple Davydov D1 Ansatz and the Generalized Coherent State

Optical-cavity manipulation strategies of singlet fission systems mediated by conical intersections: Insights from fully quantum simulations

We offer a theoretical perspective on simulation and engineering of polaritonic conical-intersection-driven singlet-fission (SF) materials. We begin by examining fundamental models, including Tavis-Cummings and Holstein-Tavis-Cummings Hamiltonians, exploring how disorder, non-Hermitian effects, and finite temperature conditions impact their dynamics, setting the stage for studying conical intersections and their crucial role in SF. Using rubrene as an example and applying the numerically accurate Davydov Ansatz methodology, we derive dynamic and spectroscopic responses of the system and demonstrate key mechanisms capable of SF manipulation, viz. cavity-induced enhancement/weakening/suppression of SF, population localization on the singlet state via engineering cavity-mode excitation, polaron/polariton decoupling, and collective enhancement of SF. We outline unsolved problems and challenges in the field and share our views on the development of the future lines of research. We emphasize the significance of careful modeling of cascades of polaritonic conical intersections in high excitation manifolds and envisage that collective geometric phase effects may remarkably affect the SF dynamics and yield. We argue that the microscopic interpretation of the main regulatory mechanisms of polaritonic conical-intersection-driven SF can substantially deepen our understanding of this process, thereby providing novel ideas and solutions for improving conversion efficiency in photovoltaics.

Coherent Transient Localization Mechanism of Interchain Exciton Transport in Regioregular P3HT: A Quantum-Dynamical Study

Transient localization has been proposed as a transport mechanism in organic materials, for both charge carriers and excitons. Here, we characterize a quantum coherent transient localization mechanism using full quantum simulations of an H-aggregated model system representative of regioregular poly(3-hexylthiophene) (rrP3HT). A Frenkel-Holstein Hamiltonian parametrized from first principles is considered, including local high-frequency modes and anharmonic, site-correlated interchain modes. Quantum-dynamical calculations are carried out using the Multi-Layer Multi-Configuration Time-Dependent Hartree (ML-MCTDH) method for a 13-site system with 195 vibrational modes, under periodic boundary conditions. It is shown that temporary localization of exciton polarons alternates with resonant transfer driven by interchain modes. While the transport process is mainly determined by exciton-polarons at the low-energy band edge, persistent coupling with the excitonic manifold is observed, giving rise to a nonadiabatic excitonic flux. This elementary transport mechanism remains preserved for limited static disorder and gives way to Anderson localization when the static disorder becomes dominant.

Simulating optical linear absorption for mesoscale molecular aggregates: An adaptive hierarchy of pure states approach

In this paper, we present dyadic adaptive HOPS (DadHOPS), a new method for calculating linear absorption spectra for large molecular aggregates. This method combines the adaptive HOPS (adHOPS) framework, which uses locality to improve computational scaling, with the dyadic HOPS method previously developed to calculate linear and nonlinear spectroscopic signals. To construct a local representation of dyadic HOPS, we introduce an initial state decomposition that reconstructs the linear absorption spectra from a sum over locally excited initial conditions. We demonstrate the sum over initial conditions can be efficiently Monte Carlo sampled and that the corresponding calculations achieve size-invariant [i.e., O(1)] scaling for sufficiently large aggregates while trivially incorporating static disorder in the Hamiltonian. We present calculations on the photosystem I core complex to explore the behavior of the initial state decomposition in complex molecular aggregates as well as proof-of-concept DadHOPS calculations on an artificial molecular aggregate inspired by perylene bis-imide to demonstrate the size-invariance of the method.

Formally Exact Simulations of Mesoscale Exciton Diffusion in a Light-Harvesting 2 Antenna Nanoarray

The photosynthetic apparatus of plants and bacteria combine atomically precise pigment-protein complexes with dynamic membrane architectures to control energy transfer on the 10-100 nm length scales. Recently, synthetic materials have integrated photosynthetic antenna proteins to enhance exciton transport, though the influence of artificial packing on the excited-state dynamics in these biohybrid materials is not fully understood. Here, we use the adaptive hierarchy of pure states (adHOPS) to perform a formally exact simulation of excitation energy transfer within artificial aggregates of light-harvesting complex 2 (LH2) with a range of packing densities. We find that LH2 aggregates support a remarkable exciton diffusion length ranging from 100 nm at a biological packing density to 300 nm at the densest packing previously suggested in an artificial aggregate. The unprecedented scale of these formally exact calculations also underscores the efficiency with which adHOPS simulates excited-state processes in molecular materials.

Dynamics of the spin-boson model: The effect of bath initial conditions

The dynamics of the (sub-)Ohmic spin-boson model under various bath initial conditions is investigated by employing the Dirac-Frenkel time-dependent variational approach with the multiple Davydov D₁ Ansatz in the interaction picture. The validity of our approach is carefully checked by comparing the results with those of the hierarchy equations of motion method. By analyzing the features of nonequilibrium dynamics, we identify the phase diagrams for different bath initial conditions. We find that for the spectral exponent s < s_c, there exists a transition from coherent to quasicoherent dynamics with increasing coupling strengths. For s_c < s ≤ 1, the coherent to incoherent crossover occurs at a certain coupling strength and the quasicoherent dynamics emerges at much larger couplings. The initial preparation of the bath has a considerable influence on the dynamics.

The hierarchy of Davydov's Ansätze: From guesswork to numerically "exact" many-body wave functions

This Perspective presents an overview of the development of the hierarchy of Davydov's Ansätze and a few of their applications in many-body problems in computational chemical physics. Davydov's solitons originated in the investigation of vibrational energy transport in proteins in the 1970s. Momentum-space projection of these solitary waves turned up to be accurate variational ground-state wave functions for the extended Holstein molecular crystal model, lending unambiguous evidence to the absence of formal quantum phase transitions in Holstein systems. The multiple Davydov Ansätze have been proposed, with increasing Ansatz multiplicity, as incremental improvements of their single-Ansatz parents. For a given Hamiltonian, the time-dependent variational formalism is utilized to extract accurate dynamic and spectroscopic properties using Davydov's Ansätze as its trial states. A quantity proven to disappear for large multiplicities, the Ansatz relative deviation is introduced to quantify how closely the Schrödinger equation is obeyed. Three finite-temperature extensions to the time-dependent variation scheme are elaborated, i.e., the Monte Carlo importance sampling, the method of thermofield dynamics, and the method of displaced number states. To demonstrate the versatility of the methodology, this Perspective provides applications of Davydov's Ansätze to the generalized Holstein Hamiltonian, variants of the spin-boson model, and systems of cavity-assisted singlet fission, where accurate dynamic and spectroscopic properties of the many-body systems are given by the Davydov trial states.

Equation-of-Motion Methods for the Calculation of Femtosecond Time-Resolved 4-Wave-Mixing and N-Wave-Mixing Signals

Femtosecond nonlinear spectroscopy is the main tool for the time-resolved detection of photophysical and photochemical processes. Since most systems of chemical interest are rather complex, theoretical support is indispensable for the extraction of the intrinsic system dynamics from the detected spectroscopic responses. There exist two alternative theoretical formalisms for the calculation of spectroscopic signals, the nonlinear response-function (NRF) approach and the spectroscopic equation-of-motion (EOM) approach. In the NRF formalism, the system-field interaction is assumed to be sufficiently weak and is treated in lowest-order perturbation theory for each laser pulse interacting with the sample. The conceptual alternative to the NRF method is the extraction of the spectroscopic signals from the solutions of quantum mechanical, semiclassical, or quasiclassical EOMs which govern the time evolution of the material system interacting with the radiation field of the laser pulses. The NRF formalism and its applications to a broad range of material systems and spectroscopic signals have been comprehensively reviewed in the literature. This article provides a detailed review of the suite of EOM methods, including applications to 4-wave-mixing and N-wave-mixing signals detected with weak or strong fields. Under certain circumstances, the spectroscopic EOM methods may be more efficient than the NRF method for the computation of various nonlinear spectroscopic signals.

Quantum dynamical study of inter-chain exciton transport in a regioregular P3HT model system at finite temperature: HJ vs. H-aggregate models

We report on quantum dynamical simulations of inter-chain exciton transport in a model of regioregular poly(3-hexylthiophene), rr-P3HT, at finite temperature, using the Multi-Layer Multi-Configuration Time-Dependent Hartree (ML-MCTDH) method for a system of up to 63 electronic states and 180 vibrational modes. A Frenkel Hamiltonian of HJ aggregate type is used, along with a reduced H-aggregate representation; electron-phonon coupling includes local high-frequency modes as well as anharmonic intermolecular modes. The latter are operative in mediating inter-chain transport, by a mechanism of transient localization type. Strikingly, this mechanism is found to be of quantum coherent character and involves non-adiabatic effects. Using periodic boundary conditions, a normal diffusion regime is identified from the exciton mean-squared displacement, apart from early-time transients. Diffusion coefficients are found to be of the order of 3 x 10-3 cm2/s, showing a non-monotonous increase with temperature.

Simulation of absorption spectra of molecular aggregates: A hierarchy of stochastic pure state approach

Simulation of spectroscopic observables for molecular aggregates with strong and structured coupling of electronic excitation to vibrational degrees of freedom is an important but challenging task. The Hierarchy of Pure States (HOPS) provides a formally exact solution based on local, stochastic trajectories. Exploiting the localization of HOPS for the simulation of absorption spectra in large aggregates requires a formulation in terms of normalized trajectories. Here, we provide a normalized dyadic equation where the ket- and bra-states are propagated in different electronic Hilbert spaces. This work opens the door to applying adaptive HOPS methods for the simulation of absorption spectra.

Variational Squeezed Davydov Ansatz for Realistic Chemical Systems with Nonlinear Vibronic Coupling

Chemical systems normally possess strong nonlinear vibronic couplings at both zero and finite temperature. For the lowest-order quadratic couplings, here, we introduce a squeezing operator into a variational coherent-state-based method, Davydov ansatz, to simulate the quantum dynamics and the respective spectroscopy. Two molecular systems, pyrazine and the 2-pyridone dimer, are taken as calculated model systems, both of which involve nontrivial quadratic vibronic couplings in high- and low-frequency regions, respectively. Upon a comparison with the benchmarks, the method manifests its advantage for nonlinear couplings. The squeezed bases are also proven to be applicable for the finite temperature by adapting with the thermofield dynamics.

Dynamics of disordered Tavis-Cummings and Holstein-Tavis-Cummings models

By employing the time-dependent variational principle and the versatile multi-D₂ Davydov trial states, in combination with the Green's function method, we study the dynamics of the Tavis-Cummings model and the Holstein-Tavis-Cummings model in the presence of diagonal disorder and cavity-qubit coupling disorder. For the Tavis-Cummings model, time evolution of the photon population, the optical absorption spectra, and the hetero-entanglement between the qubits and the cavity mode are calculated by using the Green's function method to corroborate numerically exact results of Davydov's Ansätze. For the Holstein-Tavis-Cummings model, only the latter is utilized to simulate effects of disorder on the photon population dynamics and the absorption spectra. We have demonstrated that the complementary employment of analytical and numerical methods permits uncovering a fairly comprehensive picture of a variety of complex behaviors in disordered multidimensional polaritonic cavity quantum electrodynamics systems.

Signatures of coherent vibronic exciton dynamics and conformational control in two-dimensional electronic spectroscopy of conjugated polymers

Two-dimensional electronic spectroscopy (2DES) signals for homo-oligomer J-aggregates are computed, with a focus on the role of structural change induced by low-frequency torsional modes along with quasi-stationary trapping effects induced...

Simulation of Time- and Frequency-Resolved Four-Wave-Mixing Signals at Finite Temperatures: A Thermo-Field Dynamics Approach

We propose a new approach to simulate four-wave-mixing signals of molecular systems at finite temperatures by combining the multiconfigurational Ehrenfest method with the thermo-field dynamics theory. In our approach, the four-time correlation functions at finite temperatures are mapped onto those at zero temperature in an enlarged Hilbert space with twice the vibrational degrees of freedom. As an illustration, we have simulated three multidimensional spectroscopic signals, time- and frequency-resolved fluorescence spectra, transient-absorption pump-probe spectra, and electronic two-dimensional (2D) spectra at finite temperatures, for a conical intersection-mediated singlet fission model of a rubrene crystal. It is shown that a detailed dynamical picture of the singlet fission process can be extracted from the three spectroscopic signals. An increasing temperature leads to lower intensities of the signals and broadened vibrational peaks, which can be attributed to faster singlet-triplet population transfer and stronger bath-induced electronic dephasing at higher temperatures.

Quantum Dynamics of Exciton Transport and Dissociation in Multichromophoric Systems

Due to the subtle interplay of site-to-site electronic couplings, exciton delocalization, nonadiabatic effects, and vibronic couplings, quantum dynamical studies are needed to elucidate the details of ultrafast photoinduced energy and charge transfer events in organic multichromophoric systems. In this vein, we review an approach that combines first-principles parameterized lattice Hamiltonians with accurate quantum dynamical simulations using advanced multiconfigurational methods. Focusing on the elementary transfer steps in organic functional materials, we address coherent exciton migration and creation of charge transfer excitons in homopolymers, notably representative of the poly(3-hexylthiophene) material, as well as exciton dissociation at polymer:fullerene heterojunctions. We emphasize the role of coherent transfer, trapping effects due to high-frequency phonon modes, and thermal activation due to low-frequency soft modes that drive a diffusive dynamics.

Efficient simulation of time- and frequency-resolved four-wave-mixing signals with a multiconfigurational Ehrenfest approach

We have extended the multiconfigurational Ehrenfest approach to the simulation of four-wave-mixing signals of systems involving multiple electronic and vibrational degrees of freedom. As an illustration, we calculate signals of three widely used spectroscopic techniques, time- and frequency-resolved fluorescence spectroscopy, transient absorption spectroscopy, and two-dimensional (2D) electronic spectroscopy, for a two-electronic-state, twenty-four vibrational-mode conical intersection model. It has been shown that all these three spectroscopic signals characterize fast population transfer from the higher excited electronic state to the lower excited electronic state. While the time- and frequency-resolved spectrum maps the wave packet propagation exclusively on the electronically excited states, the transient absorption and 2D electronic spectra reflect the wave packet dynamics on both electronically excited states and the electronic ground state. Combining trajectory-guided Gaussian basis functions and the nonlinear response function formalism, the present approach provides a promising general technique for the applications of various Gaussian basis methods to the calculations of four-wave-mixing spectra of polyatomic molecules.

Photon-assisted Landau-Zener transitions in a periodically driven Rabi dimer coupled to a dissipative mode

We investigate multiple photon-assisted Landau-Zener (LZ) transitions in a hybrid circuit quantum electrodynamics device in which each of two interacting transmission-line resonators is coupled to a qubit, and the qubits are driven by periodic driving fields and also coupled to a common phonon mode. The quantum state of the entire composite system is modeled using the multi-D₂Ansatz in combination with the time-dependent Dirac-Frenkel variational principle. Applying a sinusoidal driving field to one of the qubits, this device is an ideal platform to study the photon-assisted LZ transitions by comparing the dynamics of the two qubits. A series of interfering photon-assisted LZ transitions takes place if the photon frequency is much smaller than the driving amplitude. Once the two energy scales are comparable, independent LZ transitions arise and a transition pathway is revealed using an energy diagram. It is found that both adiabatic and nonadiabatic transitions are involved in the dynamics. Used to model environmental effects on the LZ transitions, the common phonon mode coupled to the qubits allows for more available states to facilitate the LZ transitions. An analytical formula is obtained to estimate the short time phonon population and produces results in reasonable agreement with numerical calculations. Equipped with the knowledge of the photon-assisted LZ transitions in the system, we can precisely manipulate the qubit state and successfully generate the qubit dynamics with a square-wave pattern by applying driving fields to both qubits, opening up new venues to manipulate the states of qubits and photons in quantum information devices and quantum computers.

Quantum dynamical simulations of intra-chain exciton diffusion in an oligo (para-phenylene vinylene) chain at finite temperature

We report on quantum dynamical simulations of exciton diffusion in an oligo(para-phenylene vinylene) chain segment with 20 repeat units (OPV-20) at finite temperature, complementary to our recent study of the same system at T = 0 K [R. Binder and I. Burghardt, J. Chem. Phys. 152, 204120 (2020)]. Accurate quantum dynamical simulations are performed using the multi-layer multi-configuration time-dependent Hartree method as applied to a site-based Hamiltonian comprising 20 electronic states of Frenkel type and 460 vibrational modes, including site-local quinoid-distortion modes along with site-correlated bond-length alternation (BLA) modes, ring torsional modes, and an explicit harmonic-oscillator bath. A first-principles parameterized Frenkel-Holstein type Hamiltonian is employed, which accounts for correlations between the ring torsional modes and the anharmonically coupled BLA coordinates located at the same junction. Thermally induced fluctuations of the torsional modes are described by a stochastic mean-field approach, and their impact on the excitonic motion is characterized in terms of the exciton mean-squared displacement. A normal diffusion regime is observed under periodic boundary conditions, apart from transient localization features. Even though the polaronic exciton species are comparatively weakly bound, exciton diffusion is found to be a coherent-rather than hopping type-process, driven by the fluctuations of the soft torsional modes. Similar to the previous observations for oligothiophenes, the evolution for the most part exhibits a near-adiabatic dynamics of local exciton ground states (LEGSs) that adjust to the local conformational dynamics. However, a second mechanism, involving resonant transitions between neighboring LEGSs, gains importance at higher temperatures.

Polaron Diffusion in Pentathienoacene Crystals

Molecular crystals have been used as prototypes for studying the energetic and dynamic properties of charge carriers in organic electronics. The growing interest in oligoacenes and fused-ring oligothiophenes in the last two decades is due, in particular, to the success achieved in conceiving pentacene-based organic photovoltaic devices. In the present work, a one-dimensional Holstein-Peierls model is designed to study the temperature-dependent polaron transport in pentathienoacene (PTA) lattices. The tight-binding Hamiltonian employed here takes into account intra and intermolecular electron-lattice interactions. Results reveal that polarons in PTAs can be stable structures even at high temperatures, about 400 K. During the dynamical process, these charge carriers present a typical 1D random walk diffusive motion with a low activation energy of 13 meV and a room temperature diffusivity constant of 1.07 × 10⁻³ cm² s⁻¹. Importantly, these critical values for the polaron diffusion and activation energy are related to the choice of model parameters, which are adopted to describe pristine lattices.

Density matrix dynamics in twin-formulation: An efficient methodology based on tensor-train representation of reduced equations of motion

The twin-formulation of quantum statistical mechanics is employed to describe a new methodology for the solution of the equations of motion of the reduced density matrix in their hierarchical formulation. It is shown that the introduction of tilde operators and of their algebra in the dual space greatly simplifies the application of numerical techniques for the propagation of the density matrix. The application of tensor-train representation of a vector to solve complex quantum dynamical problems within the framework of the twin-formulation is discussed. Next, applications of the hierarchical equations of motion to a dissipative polaron model are presented showing the validity and accuracy of the new approach.

Dissipative dynamics in a tunable Rabi dimer with periodic harmonic driving

Recent progress on qubit manipulation allows application of periodic driving signals on qubits. In this study, a harmonic driving field is added to a Rabi dimer to engineer photon and qubit dynamics in a circuit quantum electrodynamics device. To model environmental effects, qubits in the Rabi dimer are coupled to a phonon bath with a sub-Ohmic spectral density. A nonperturbative treatment, the Dirac-Frenkel time-dependent variational principle together with the multiple Davydov D₂ ansatz, is employed to explore the dynamical behavior of the tunable Rabi dimer. In the absence of the phonon bath, the amplitude damping of the photon number oscillation is greatly suppressed by the driving field, and photons can be created, thanks to the resonance between the periodic driving field and the photon frequency. In the presence of the phonon bath, one can still change the photon numbers in two resonators and indirectly alter the photon imbalance in the Rabi dimer by directly varying the driving signal in one qubit. It is shown that qubit states can be manipulated directly by the harmonic driving. The environment is found to strengthen the interqubit asymmetry induced by the external driving, opening up a new venue to engineer the qubit states.