Jeremy M. Moix

A hybrid stochastic hierarchy equations of motion approach to treat the low temperature dynamics of non-Markovian open quantum systems

The hierarchical equations of motion technique has found widespread success as a tool to generate the numerically exact dynamics of non-Markovian open quantum systems. However, its application to low temperature environments remains a serious challenge due to the need for a deep hierarchy that arises from the Matsubara expansion of the bath correlation function. Here we present a hybrid stochastic hierarchical equation of motion (sHEOM) approach that alleviates this bottleneck and leads to a numerical cost that is nearly independent of temperature. Additionally, the sHEOM method generally converges with fewer hierarchy tiers allowing for the treatment of larger systems. Benchmark calculations are presented on the dynamics of two level systems at both high and low temperatures to demonstrate the efficacy of the approach. Then the hybrid method is used to generate the exact dynamics of systems that are nearly impossible to treat by the standard hierarchy. First, exact energy transfer rates are calculated across a broad range of temperatures revealing the deviations from the Förster rates. This is followed by computations of the entanglement dynamics in a system of two qubits at low temperature spanning the weak to strong system-bath coupling regimes.

Equilibrium-reduced density matrix formulation : influence of noise, disorder, and temperature on localization in excitonic systems

An exact method to compute the entire equilibrium-reduced density matrix for systems characterized by a system-bath Hamiltonian is presented. The approach is based upon a stochastic unraveling of the influence functional that appears in the imaginary time path integral formalism of quantum statistical mechanics. This method is then applied to study the effects of thermal noise, static disorder, and temperature on the coherence length in excitonic systems. As representative examples of biased and unbiased systems, attention is focused on the well-characterized complexes of the Fenna-Matthews-Olson (FMO) protein and the light harvesting complex of purple bacteria, LH2, respectively. Due to the bias, FMO is completely localized in the site basis at low temperatures, whereas LH2 is completely delocalized. In the latter, the presence of static disorder leads to a plateau in the coherence length at low temperature that becomes increasingly pronounced with increasing strength of the disorder. The introduction of noise, however, precludes this effect. In biased systems, it is shown that the environment may increase the coherence length, but only decrease that of unbiased systems. Finally it is emphasized that for typical values of the environmental parameters in light harvesting systems, the system and bath are entangled at equilibrium in the single excitation manifold. That is, the density matrix cannot be described as a product state as is often assumed, even at room temperature. The reduced density matrix of LH2 is shown to be in precise agreement with the steady state limit of previous exact quantum dynamics calculations.

Accuracy of second order perturbation theory in the polaron and variational polaron frames

In the study of open quantum systems, the polaron transformation has recently attracted a renewed interest as it offers the possibility to explore the strong system-bath coupling regime. Despite this interest, a clear and unambiguous analysis of the regimes of validity of the polaron transformation is still lacking. Here we provide such a benchmark, comparing second order perturbation theory results in the original untransformed frame, the polaron frame, and the variational extension with numerically exact path integral calculations of the equilibrium reduced density matrix. Equilibrium quantities allow a direct comparison of the three methods without invoking any further approximations as is usually required in deriving master equations. It is found that the second order results in the original frame are accurate for weak system-bath coupling; the results deteriorate when the bath cut-off frequency decreases. The full polaron results are accurate for the entire range of coupling for a fast bath but only in the strong coupling regime for a slow bath. The variational method is capable of interpolating between these two methods and is valid over a much broader range of parameters.

Identifying reactive trajectories using a moving transition state.

A time-dependent no-recrossing dividing surface is shown to lead to a new criterion for identifying reactive trajectories well before they are evolved to infinite time. Numerical dynamics simulations of a dissipative anharmonic two-dimensional system confirm the efficiency of this approach. The results are compared to the standard fixed transition state dividing surface that is well-known to suffer from recrossings and therefore requires trajectories to be evolved over a long time interval before they can reliably be classified as reactive or nonreactive. The moving dividing surface can be used to identify reactive trajectories in harmonic or moderately anharmonic systems with considerably lower numerical effort or even without any simulation at all.

Coherent quantum transport in disordered systems: I. The influence of dephasing on the transport properties and absorption spectra on one-dimensional systems

Excitonic transport in static-disordered one dimensional systems is studied in the presence of thermal fluctuations that are described by the Haken–Strobl–Reineker model. For short times, non-diffusive behavior is observed that can be characterized as the free-particle dynamics on the length-scale bounded by the Anderson localized system. Over longer time scales, the environment-induced dephasing is sufficient to overcome the Anderson localization caused by the disorder and allow for transport to occur which is always seen to be diffusive. In the limiting regimes of weak and strong dephasing quantum master equations are developed, and their respective scaling relations imply the existence of a maximum in the diffusion constant as a function of the dephasing rate that is confirmed numerically. In the weak dephasing regime, it is demonstrated that the diffusion constant is proportional to the square of the localization length which leads to a significant enhancement of the transport rate over the classical prediction. Finally, the influence of noise and disorder on the absorption spectrum is presented and its relationship to the transport properties is discussed.

Förster resonance energy transfer, absorption and emission spectra in multichromophoric systems. III. Exact stochastic path integral evaluation

A numerically exact path integral treatment of the absorption and emission spectra of open quantum systems is presented that requires only the straightforward solution of a stochastic differential equation. The approach converges rapidly enabling the calculation of spectra of large excitonic systems across the complete range of system parameters and for arbitrary bath spectral densities. With the numerically exact absorption and emission operators, one can also immediately compute energy transfer rates using the multi-chromophoric Förster resonant energy transfer formalism. Benchmark calculations on the emission spectra of two level systems are presented demonstrating the efficacy of the stochastic approach. This is followed by calculations of the energy transfer rates between two weakly coupled dimer systems as a function of temperature and system-bath coupling strength. It is shown that the recently developed hybrid cumulant expansion (see Paper II) is the only perturbative method capable of generating uniformly reliable energy transfer rates and emission spectra across a broad range of system parameters.

Semiclassical initial value series representation in the continuum limit: Application to vibrational relaxation

A recently formulated continuum limit semiclassical initial value series representation (SCIVR) of the quantum dynamics of dissipative systems is applied to the study of vibrational relaxation of model harmonic and anharmonic oscillator systems. As is well known, the classical dynamics of dissipative systems may be described in terms of a generalized Langevin equation. The continuum limit SCIVR uses the Langevin trajectories as input, albeit with a quantum noise rather than a classical noise. Combining this development with the forward-backward form of the prefactor-free propagator leads to a tractable scheme for computing quantum thermal correlation functions. Here we present the first implementation of this continuum limit SCIVR series method to study two model problems of vibrational relaxation. Simulations of the dissipative harmonic oscillator system over a wide range of parameters demonstrate that at most only the first two terms in the SCIVR series are needed for convergence of the correlation function. The methodology is then applied to the vibrational relaxation of a dissipative Morse oscillator. Here, too, the SCIVR series converges rapidly as the first two terms are sufficient to provide the quantum mechanical relaxation with an estimated accuracy on the order of a few percent. The results in this case are compared with computations obtained using the classical Wigner approximation for the relaxation dynamics.

Förster resonance energy transfer, absorption and emission spectra in multichromophoric systems. II. Hybrid cumulant expansion

We develop a hybrid cumulant expansion method to account for the system-bath entanglement in the emission spectrum in the multi-chromophoric Förster transfer rate. In traditional perturbative treatments, the emission spectrum is usually expanded with respect to the system-bath coupling term in both real and imaginary time. This perturbative treatment gives a reliable absorption spectrum, where the bath is Gaussian and only the real-time expansion is involved. For the emission spectrum, the initial state is an entangled state of the system plus bath. Traditional perturbative methods are problematic when the excitations are delocalized and the energy gap is larger than the thermal energy, since the second-order expansion cannot predict the displacement of the bath. In the present method, the real-time dynamics is carried out by using the 2nd-order cumulant expansion method, while the displacement of the bath is treated more accurately by utilizing the exact reduced density matrix of the system. In a sense, the hybrid cumulant expansion is based on a generalized version of linear response theory with entangled initial states.

Coherent quantum transport in disordered systems: A unified polaron treatment of hopping and band-like transport

Quantum transport in disordered systems is studied using a polaron-based master equation. The polaron approach is capable of bridging the results from the coherent band-like transport regime governed by the Redfield equation to incoherent hopping transport in the classical regime. A non-monotonic dependence of the diffusion coefficient is observed both as a function of temperature and system-phonon coupling strength. In the band-like transport regime, the diffusion coefficient is shown to be linearly proportional to the system-phonon coupling strength and vanishes at zero coupling due to Anderson localization. In the opposite classical hopping regime, we correctly recover the dynamics described by the Fermis Golden Rule and establish that the scaling of the diffusion coefficient depends on the phonon bath relaxation time. In both the hopping and band-like transport regimes, it is demonstrated that at low temperature, the zero-point fluctuations of the bath lead to non-zero transport rates and hence a finite diffusion constant. Application to rubrene and other organic semiconductor materials shows a good agreement with experimental mobility data.

Quantum Diffusion on Molecular Tubes: Universal Scaling of the 1D to 2D Transition

The transport properties of disordered systems are known to depend critically on dimensionality. We study the diffusion coefficient of a quantum particle confined to a lattice on the surface of a tube, where it scales between the 1D and 2D limits. It is found that the scaling relation is universal and independent of the temperature, disorder, and noise parameters, and the essential order parameter is the ratio between the localization length in 2D and the circumference of the tube. Phenomenological and quantitative expressions for transport properties as functions of disorder and noise are obtained and applied to real systems: In the natural chlorosomes found in light-harvesting bacteria the exciton transfer dynamics is predicted to be in the 2D limit, whereas a family of synthetic molecular aggregates is found to be in the homogeneous limit and is independent of dimensionality.