Nonperturbative spin–boson and spin–spin dynamics and nonlinear Fano interferences: A unified dissipaton theory based study

Generalized system-bath entanglement theorem for Gaussian environments

The entanglement between system and bath often plays a pivotal role in complex systems spanning multiple orders of magnitude. A system-bath entanglement theorem was previously established for Gaussian environments in J. Chem. Phys. 152, 034102 (2020) regarding linear response functions. This theorem connects the entangled responses to the local system and bare bath properties. In this work, we generalize it to correlation functions. Key steps in derivations involve using the generalized Langevin dynamics for hybridizing bath modes and the Bogoliubov transformation that maps the original finite-temperature reservoir to an effective zero-temperature vacuum by employing an auxiliary bath. The generalized theorem allows us to evaluate the system-bath entangled correlations and the bath mode correlations in the total composite space, as long as we know the bare-bath statistical properties and obtain the reduced system correlations. To demonstrate the cross-scale entanglements, we utilize the generalized theorem to calculate the solvation free energy of an electron transfer system with intramolecular vibrational modes.

Predicting rate kernels via dynamic mode decomposition

Simulating dynamics of open quantum systems is sometimes a significant challenge, despite the availability of various exact or approximate methods. Particularly when dealing with complex systems, the huge computational cost will largely limit the applicability of these methods. In this work, we investigate the usage of dynamic mode decomposition (DMD) to evaluate the rate kernels in quantum rate processes. DMD is a data-driven model reduction technique that characterizes the rate kernels using snapshots collected from a small time window, allowing us to predict the long-term behaviors with only a limited number of samples. Our investigations show that whether the external field is involved or not, the DMD can give accurate prediction of the result compared with the traditional propagations, and simultaneously reduce the required computational cost.

Open quantum systems with nonlinear environmental backactions: Extended dissipaton theory vs core-system hierarchy construction

In this paper, we present a comprehensive account of quantum dissipation theories with the quadratic environment couplings. The theoretical development includes the Brownian solvation mode embedded hierarchical quantum master equations, a core-system hierarchy construction that verifies the extended dissipaton equation of motion (DEOM) formalism [R. X. Xu et al., J. Chem. Phys. 148, 114103 (2018)]. Developed are also the quadratic imaginary-time DEOM for equilibrium and the λ(t)-DEOM for nonequilibrium thermodynamics problems. Both the celebrated Jarzynski equality and Crooks relation are accurately reproduced, which, in turn, confirms the rigorousness of the extended DEOM theories. While the extended DEOM is more numerically efficient, the core-system hierarchy quantum master equation is favorable for "visualizing" the correlated solvation dynamics.

Coherent excitation energy transfer in model photosynthetic reaction center: Effects of non-Markovian quantum environment

Excitation energy transfer (EET) and electron transfer (ET) are crucially involved in photosynthetic processes. In reality, the photosynthetic reaction center constitutes an open quantum system of EET and ET, which manifests an interplay of pigments, solar light and phonon baths. So far theoretical studies have been mainly based on master equation approaches in the Markovian condition. The non-Markovian environmental effect, which may play a crucial role, has not been sufficiently considered. In this work, we propose a mixed dynamic approach to investigate this open system. The influence of phonon bath is treated via the exact dissipaton equation of motion (DEOM) while that of photon bath is via the Lindblad master equation. Specifically, we explore the effect of non-Markovian quantum phonon bath on the coherent transfer dynamics and its manipulation on the current--voltage behavior. Distinguished from the results of completely Markovian Lindblad equation and those adopting classical environment description, the mixed DEOM--Lindblad simulations exhibittransfer coherence up to a few hundreds femtosecondsand the related environmental manipulation effect on current.These non-Markovian quantum coherent effects may be extended tomore complex and realistic systems and be helpful to thedesign of organic photovoltaic devices.

Electron Transfer under the Floquet Modulation in Donor-Bridge-Acceptor Systems

Electron transfer (ET) processes are of broad interest in modern chemistry. With the advancements of experimental techniques, one may modulate the ET via such events as light-matter interactions. In this work, we study the ET under a Floquet modulation occurring in the donor-bridge-acceptor systems, with the rate kernels projected out from the exact dissipaton equation of motion formalism. This together with the Floquet theorem enables us to investigate the interplay between the intrinsic non-Markovianity and the driving periodicity. The observed rate kernel exhibits a Herzberg-Teller-like mechanism induced by the bridge fluctuation subject to effective modulation.

Application of the imaginary time hierarchical equations of motion method to calculate real time correlation functions

We investigate the application of the imaginary time hierarchical equations of motion method to calculate real time quantum correlation functions. By starting from the path integral expression for the correlated system–bath equilibrium state, we first derive a new set of equations that decouple the imaginary time propagation and the calculation of auxiliary density operators. The new equations, thus, greatly simplify the calculation of the equilibrium correlated initial state that is subsequently used in the real time propagation to obtain the quantum correlation functions. It is also shown that a periodic decomposition of the bath imaginary time correlation function is no longer necessary in the new equations such that different decomposition schemes can be explored. The applicability of the new method is demonstrated in several numerical examples, including the spin-Boson model, the Holstein model, and the double-well model for proton transfer reaction.

Quantum dissipation with nonlinear environment couplings: Stochastic fields dressed dissipaton equation of motion approach

Accurate and efficient simulation on quantum dissipation with nonlinear environment couplings remains a challenging task nowadays. In this work, we propose to incorporate the stochastic fields, which resolve just the nonlinear environment coupling terms, into the dissipaton-equation-of-motion (DEOM) construction. The stochastic fields are introduced via the Hubbard-Stratonovich transformation. After the transformation, the resulted stochastic-fields-dressed (SFD) total Hamiltonian contains only linear environment coupling terms. On the basis of that, SFD-DEOM can then be constructed. The resultant SFD-DEOM, together with the ensemble average over the stochastic fields, constitutes an exact and nonperturbative approach to quantum dissipation under nonlinear environment couplings. It is also of relatively high efficiency and stability due to the fact that only nonlinear environment coupling terms are dealt with stochastic fields, while linear couplings are still treated as the usual DEOM. Numerical performance and demonstrations are presented with a two-state model system.

Role of Pigment-Protein Coupling in the Energy Transport Dynamics in the Fenna-Matthews-Olson Complex

The role of pigment-protein coupling in the dynamics of photosynthetic energy transport in chromophoric complexes has not been fully understood. The excitation energy transfer in the photosynthetic system is tremendously efficient. In particular, we investigate the excitation energy transport in the Fenna-Matthews-Olson (FMO) complex. The exciton dynamics and excitation energy transfer (EET) depend on the interaction between the excited chromophores and their environment. Most theoretical models believe that all bacteriochlorophyll-a (BChla) sites are surrounded by the same local protein environment, which is contradicted by the structural analysis of the FMO complex. Based on different values of pigment-protein coupling for different sites, measured in the adiabatic limit, we have theoretically investigated the effect of the heterogeneous local protein environment on the EET process. By the realistic and site-dependent model of the system-bath couplings, the results show that this interaction may have a critical value for the coherent energy-transfer process. Furthermore, we verify that the two transport pathways are coherent and stable to the important parameter reorganization energy of environmental interactions. The quantum dynamical simulations show that the correlation fluctuation keeps the oscillation of the coherent excitation on a long timescale. In addition, due to the inhomogeneous pigment-protein coupling, different BChl sites have asymmetric excitation oscillation timescales.

Theoretical Study of the Spectral Differences of the Fenna-Matthews-Olson Protein from Different Species and Their Mutants

The structural basis for the spectral differences between the Fenna-Matthews-Olson (FMO) proteins from Chlorobaculum tepidum (C. tepidum) and Prosthecochloris aestuarii 2K (P. aestuarii) is yet to be fully understood. Mutation-induced perturbation to the exciton structure and the optical spectra of the complex provide a suitable means to investigate the critical role played by the protein scaffold. In this work, we have performed quantum-mechanics/molecular-mechanics calculations over the molecular dynamics simulation trajectories with the polarized protein-specific charge scheme for both wild-type FMOs and two mutants. Our result reveals that a single-point mutation in the vicinity of BChl 6, namely, W183F of C. tepidum, significantly affects the absorption spectrum, resulting in a switch of the absorption spectrum from type 2, for which the 806 nm band is more pronounced than the 815 nm band, to type 1, for which the 815 nm band is pronounced. Our observations agree with the single-point mutation experiments reported by Saer et al. (Biochim. Biophys. Acta, Bioenerg. 2017, 1858, 288-296) and Khmelnitskiy et al. (J. Phys. Chem. Lett. 2018, 9, 3378-3386). In contrast, the absorption spectrum of the P. aestuarii experiences the opposite transition (from type 1 to type 2) upon the same mutation. Furthermore, by comparing the contributions of individual pigments to the spectra in the wild type and its mutant, we find that a single-point mutation near BChl 6 not only induces changes in excitation energy of BChl 6 per se but also affects the excitonic structures of the neighboring BChls 5 and 7 through strong interpigment electronic couplings, resulting in a significant change in the absorption spectra.

Correlated vibration-solvent effects on the non-Condon exciton spectroscopy

Excitation energy transfer is crucially involved in a variety of systems. During the process, the non-Condon vibronic coupling and the surrounding solvent interaction may synergetically play important roles. In this work, we study the correlated vibration-solvent influences on the non-Condon exciton spectroscopy. Statistical analysis is elaborated for the overall vibration-plus-solvent environmental effects. Analytic solutions are derived for the linear absorption of monomer systems. General simulations are accurately carried out via the dissipaton-equation-of-motion approach. The resulted spectra in either the linear absorption or strong field regime clearly demonstrate the coherence enhancement due to the synergetic vibration-solvent correlation.

Equilibrium and transient thermodynamics: A unified dissipaton-space approach

This work presents a unified dissipaton-equation-of-motion (DEOM) theory and its evaluations on the Helmholtz free energy change due to the isotherm mixing of two isolated subsystems. One is a local impurity, and the other is a nonlocal Gaussian bath. DEOM constitutes a fundamental theory for such open quantum mixtures. To complete the theory, we also construct the imaginary-time DEOM formalism via an analytical continuation of dissipaton algebra, which would be limited to equilibrium thermodynamics. On the other hand, the real-time DEOM deals with both equilibrium structural and nonequilibrium dynamic properties. Its combination with the thermodynamic integral formalism would be a viable and accurate means to both equilibrium and transient thermodynamics. As illustrations, we report the numerical results on a spin-boson system, with elaborations on the underlying anharmonic features, the thermodynamic entropy vs the von Neumann entropy, and an indication of "solvent-cage" formation. Beside the required asymptotic equilibrium properties, the proposed transient thermodynamics also supports the basic spontaneity criterion.

Theoretical Study on the Effect of Environment on Excitation Energy Transfer in Photosynthetic Light-Harvesting Systems

In recent years, two-dimensional (2D) electronic spectroscopy experiments prove that the excitation energy transfer (EET) in photosynthetic light-harvesting systems presents long-lived electronic quantum beating signals. After being discovered in the light-harvesting system, the quantum coherence effect has aroused widespread discussion. To illustrate the EET process in the Fenna-Matthews-Olson (FMO) and phycocyanin 645 (PC645) complex, the local protein environment is often thought to be the same; however, this is ambivalent to the practical structural analysis of the light-harvesting complex. By adopting the dissipaton equation of motion theory, we present the effect of a heterogeneous protein environment on the energy transfer process with accurate numerical results. We demonstrate that the energy transfer process relies on the local heterogeneous environment for the FMO complex. A similar good agreement is found for the PC645 complex. Furthermore, we discuss the optimal value of different chromophores in the excitation energy transfer process by controlling the environmental characteristics.

A unified stochastic formulation of dissipative quantum dynamics. II. Beyond linear response of spin baths

We use the "generalized hierarchical equation of motion" proposed in Paper I [C.-Y. Hsieh and J. Cao, J. Chem. Phys. 148, 014103 (2018)] to study decoherence in a system coupled to a spin bath. The present methodology allows a systematic incorporation of higher-order anharmonic effects of the bath in dynamical calculations. We investigate the leading order corrections to the linear response approximations for spin bath models. Two kinds of spin-based environments are considered: (1) a bath of spins discretized from a continuous spectral density and (2) a bath of localized nuclear or electron spins. The main difference resides with how the bath frequency and the system-bath coupling parameters are distributed in an environment. When discretized from a continuous spectral density, the system-bath coupling typically scales as ∼1/N_B where N_B is the number of bath spins. This scaling suppresses the non-Gaussian characteristics of the spin bath and justifies the linear response approximations in the thermodynamic limit. For the nuclear/electron spin bath models, system-bath couplings are directly deduced from spin-spin interactions and do not necessarily obey the 1/N_B scaling. It is not always possible to justify the linear response approximations in this case. Furthermore, if the spin-spin Hamiltonian is highly symmetrical, there exist additional constraints that generate highly non-Markovian and persistent dynamics that is beyond the linear response treatments.

Theories of quantum dissipation and nonlinear coupling bath descriptors

The quest of an exact and nonperturbative treatment of quantum dissipation in nonlinear coupling environments remains in general an intractable task. In this work, we address the key issues toward the solutions to the lowest nonlinear environment, a harmonic bath coupled both linearly and quadratically with an arbitrary system. To determine the bath coupling descriptors, we propose a physical mapping scheme, together with the prescription reference invariance requirement. We then adopt a recently developed dissipaton equation of motion theory [R. X. Xu et al., Chin. J. Chem. Phys. 30, 395 (2017)], with the underlying statistical quasi-particle ("dissipaton") algebra being extended to the quadratic bath coupling. We report the numerical results on a two-level system dynamics and absorption and emission line shapes.

Onsets of hierarchy truncation and self-consistent Born approximation with quantum mechanics prescriptions invariance

The issue of efficient hierarchy truncation is related to many approximate theories. In this paper, we revisit this issue from both the numerical efficiency and quantum mechanics prescription invariance aspects. The latter requires that the truncation approximation made in Schrödinger picture, such as the quantum master equations and their self-consistent-Born-approximation improvements, should be transferable to their Heisenberg-picture correspondences, without further approximations. We address this issue with the dissipaton equation of motion (DEOM), which is a unique theory for the dynamics of not only reduced systems but also hybrid bath environments. We also highlight the DEOM theory is not only about how its dynamical variables evolve in time, but also the underlying dissipaton algebra. We demonstrate this unique feature of DEOM with model systems and report some intriguing nonlinear Fano interferences characteristics that are experimentally measurable.