Exact quantum master equation for a molecular aggregate coupled to a harmonic bath

Reconciling nonlinear dissipation with the bilinear model of two Brownian particles

The Brownian motion of a single particle is a paradigmatic model of the nonequilibrium dynamics of dissipative systems. In the system-plus-reservoir approach, one can derive the particle's equations of motion from the reversible dynamics of the system coupled to a bath of oscillators representing its thermal environment. However, extending the system-plus-reservoir approach to multiple particles in a collective environment is not straightforward, and conflicting models have been proposed to that end. Here, we set out to reconcile some aspects of the nonlinear and the bilinear models of two Brownian particles. We show how the nonlinear dissipation originally derived from exponential system-reservoir couplings can alternatively be obtained from the bilinear Lagrangian, with a modified spectral function that explicitly depends on the distance between the particles. We discuss applications to the contexts of anomalous diffusion and of hydrodynamic interactions. Our results thus broaden the applicability of the bilinear model.

On simulating the dynamics of electronic populations and coherences via quantum master equations based on treating off-diagonal electronic coupling terms as a small perturbation

Quantum master equations provide a general framework for describing the dynamics of electronic observables within a complex molecular system. One particular family of such equations is based on treating the off-diagonal coupling terms between electronic states as a small perturbation within the framework of second-order perturbation theory. In this paper, we show how different choices of projection operators, as well as whether one starts out with the time-convolution or the time-convolutionless forms of the generalized quantum master equation, give rise to four different types of such off-diagonal quantum master equations (OD-QMEs), namely, time-convolution and time-convolutionless versions of a Pauli-type OD-QME for only the electronic populations and an OD-QME for the full electronic density matrix (including both electronic populations and coherences). The fact that those OD-QMEs are given in terms of the interaction picture makes it non-trivial to obtain Schrödinger picture electronic coherences from them. To address this, we also extend a procedure for extracting Schrödinger picture electronic coherences from interaction picture populations recently introduced by Trushechkin in the context of time-convolutionless Pauli-type OD-QME to the other three types of OD-QMEs. The performance of the aforementioned four types of OD-QMEs is explored in the context of the Garg-Onuchic-Ambegaokar benchmark model for charge transfer in the condensed phase across a relatively wide parameter range. The results show that time-convolution OD-QMEs can be significantly more accurate than their time-convolutionless counterparts, particularly in the case of Pauli-type OD-QMEs, and that rather accurate Schrödinger picture coherences can be obtained from interaction picture electronic inputs.

Hierarchical Equations-of-Motion Method for Momentum System-Bath Coupling

For a broad class of quantum models of practical interest, we demonstrate that the Hamiltonian of the system nonlinearly coupled to a harmonic bath through the system and bath coordinates can be equivalently mapped into the Hamiltonian of the system bilinearly coupled to the bath through the system and bath momenta. We show that the Hamiltonian with bilinear system-bath momentum coupling can be treated by the hierarchical equations-of-motion (HEOM) method and present the corresponding proof-of-principle simulations. The developed methodology creates the opportunity to scrutinize a new family of nonlinear quantum systems by the numerically accurate HEOM method.

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.

Origin of Unexpectedly Simple Oscillatory Responses in the Excited-State Dynamics of Disordered Molecular Aggregates

Unraveling the many facets of coherent and incoherent exciton motion in an ensemble of chromophores is an inherently complex quantum mechanical problem that has triggered a vivid debate on the role of quantum effects in molecular materials and biophysical systems. Here the dynamics of a statistical ensemble of molecular aggregates consisting of identical chromophores is investigated within a new theoretical framework. Taking account of intrinsic properties of the system, the Hamiltonian of the aggregate is partitioned into two mutually commuting vibrational and vibronic operators. This representation paves the way for an analysis that reveals the role of static disorder in ensembles of aggregates. Using analytical methods, it is demonstrated that after a critical time τ_D ≃ 2π/σ (σ being the dispersion of the disorder) any dynamic variable of the aggregate exhibits purely vibrational dynamics. This result is confirmed by exact numerical calculations of the time-dependent site populations of the aggregate. These findings may be useful for the interpretation of recent femtosecond spectroscopic experiments on molecular aggregates.

Simulation of Quantum Dynamics of Excitonic Systems at Finite Temperature: an efficient method based on Thermo Field Dynamics

Quantum electron-vibrational dynamics in molecular systems at finite temperature is described using an approach based on Thermo Field Dynamics theory. This formulation treats temperature effects in the Hilbert space without introducing the Liouville space. The solution of Thermo Field Dynamics equations with a novel technique for the propagation of Tensor Trains (Matrix Product States) is implemented and discussed. The methodology is applied to the study of the exciton dynamics in the Fenna-Mathews-Olsen complex using a realistic structured spectral density to model the electron-phonon interaction. The results of the simulations highlight the effect of specific vibrational modes on the exciton dynamics and energy transfer process, as well as call for careful modeling of electron-phonon couplings.

Exact vs. asymptotic spectral densities in the Garg-Onuchic-Ambegaokar charge transfer model and its effect on Fermi's golden rule rate constants

The Garg-Onuchic-Ambegaokar model [J. Chem. Phys. 83, 4491 (1985)] has been used extensively for benchmarking methods aimed at calculating charge transfer rates. Within this model, the donor and acceptor diabats are described as shifted parabolas along a single primary mode, which is bilinearly coupled to a harmonic bath consisting of secondary modes, characterized by an Ohmic spectral density with exponential cutoff. Rate calculations for this model are often performed in the normal mode representation, with the corresponding effective spectral density given by an asymptotic expression derived at the limit where the Ohmic bath cutoff frequency is much larger than the primary mode frequency. We compare Fermi's golden rule rate constants obtained with the asymptotic and exact effective spectral densities. We find significant deviations between rate constants obtained from the asymptotic spectral density and those obtained from the exact one in the deep inverted region. Within the range of primary mode frequencies commonly employed, we find that the discrepancies increase with decreasing temperature and with decreasing primary mode frequency.

Coherent Exciton Dynamics in the Presence of Underdamped Vibrations

Recent ultrafast optical experiments show that excitons in large biological light-harvesting complexes are coupled to molecular vibration modes. These high-frequency vibrations will not only affect the optical response, but also drive the exciton transport. Here, using a model dimer system, the frequency of the underdamped vibration is shown to have a strong effect on the exciton dynamics such that quantum coherent oscillations in the system can be present even in the case of strong noise. Two mechanisms are identified to be responsible for the enhanced transport efficiency: critical damping due to the tunable effective strength of the coupling to the bath, and resonance coupling where the vibrational frequency coincides with the energy gap in the system. The interplay of these two mechanisms determines parameters responsible for the most efficient transport, and these optimal control parameters are comparable to those in realistic light-harvesting complexes. Interestingly, oscillations in the excitonic coherence at resonance are suppressed in comparison to the case of an off-resonant vibration.

Polaron dynamics in two-dimensional photon-echo spectroscopy of molecular rings

We have developed a new approach to the computation of third-order spectroscopic signals of molecular rings, by incorporating the Davydov soliton theory into the nonlinear response function formalism. The Davydov D1 and D Ansätze have been employed to treat the interactions between the excitons and the primary phonons, allowing for a full description of arbitrary exciton-phonon coupling strengths. As an illustration, we have simulated a series of optical 2D spectra for two models of molecular rings.

Optically promoted bipartite atomic entanglement in hybrid metallic carbon nanotube systems

We study theoretically a pair of spatially separated extrinsic atomic type species (extrinsic atoms, ions, molecules, or semiconductor quantum dots) near a metallic carbon nanotube, that are coupled both directly via the inter-atomic dipole-dipole interactions and indirectly by means of the virtual exchange by resonance plasmon excitations on the nanotube surface. We analyze how the optical preparation of the system by using strong laser pulses affects the formation and evolution of the bipartite atomic entanglement. Despite a large number of possible excitation regimes and evolution pathways, we find a few generic scenarios for the bipartite entanglement evolution and formulate practical recommendations on how to optimize and control the robust bipartite atomic entanglement in hybrid carbon nanotube systems.