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

Abstract

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 light harvesting complexes of FMO and 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. * Electronic address: jianshu@mit.edu 2

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

“…26,[32][33][34][35][36] Using a Hubbard-Stratonovich transformation, 26,32,34,37 it was shown that the influence functional may be unraveled by an auxiliary stochastic field. The ensuing imaginary time evolution may then be interpreted as one governed by a time-dependent…”

Efficient energy transfer in light-harvesting systems: Quantum-classical comparison, flux network, and robustness analysis

“…Without trapping, hopping kinetics always imposes P n (t → ∞) ∝ exp(−βε n ). Instead, the quantum steady-state distribution is evaluated by ρ(t → ∞) ∝ Tr B {exp(−βH tot )}, which arises from the rigorous quantum Boltzmann distribution with the consideration of both system and bath [53]. The steady-state population of the lowest energy trap site can be decreased, implying a lower probability of energy being trapped (energy transfer efficiency).…”

A continued fraction resummation form of bath relaxation effect in the spin-boson model

“…Consequently, the equilibrium population predicted by the RQKE rate can deviate from the classical Boltzmann distribution and approach to the exact quantum Boltzmann distribution, P eq;n ∝ [Tr B {exp(− βH tot )}] nn . 39,40 As a verification, we extend our previous study at a high temperature T = 300 K to lower temperatures. Since the equilibrium population is always one half in the unbiased system, we only consider the biased system, ε 12 = 100 cm −1 , with J = 100 cm −1 and λ = 100 cm −1 .…”

“…The hierarchy equation with the Matsubara frequency summation is used to obtain the exact equilibrium population, which is numerically the same as the result of the stochastic path integral. 39,40 Our numerical calculation shows that each correction term δ 2 j(<k−1) is different for the forward and backward rates, and the deviation increases as temperature decreases. As shown in Fig.…”

“…V, the temperature-dependent equilibrium population is calculated using the RQKE rates, which is also compared with the exact stochastic path integral result. 39,40 In Sec. VI, we summarize our studies.…”