Lipeng Chen

Variational dynamics of the sub-Ohmic spin-boson model on the basis of multiple Davydov D1 states

Dynamics of the sub-Ohmic spin-boson model is investigated by employing a multitude of the Davydov D1 trial states, also known as the multi-D1 Ansatz. Accuracy in dynamics simulations is improved significantly over the single D1 Ansatz, especially in the weak system-bath coupling regime. The reliability of the multi-D1 Ansatz for various coupling strengths and initial conditions is also systematically examined, with results compared closely with those of the hierarchy equations of motion and the path integral Monte Carlo approaches. In addition, a coherent-incoherent phase crossover in the nonequilibrium dynamics is studied through the multi-D1 Ansatz. The phase diagram is obtained with a critical point sc = 0.4. For sc < s < 1, the coherent-to-incoherent crossover occurs at a certain coupling strength, while the coherent state recurs at a much larger coupling strength. For s < sc, only the coherent phase exists.

Effect of high-frequency modes on singlet fission dynamics

Singlet fission is a spin-allowed energy conversion process whereby a singlet excitation splits into two spin-correlated triplet excitations residing on adjacent molecules and has a potential to dramatically increase the efficiency of organic photovoltaics. Recent time-resolved nonlinear spectra of pentacene derivatives have shown the importance of high frequency vibrational modes in efficient fission. In this work, we explore impacts of vibration-induced fluctuations on fission dynamics through quantum dynamics calculations with parameters from fitting measured linear and nonlinear spectra. We demonstrate that fission dynamics strongly depends on the frequency of the intramolecular vibrational mode. Furthermore, we examine the effect of two vibrational modes on fission dynamics. Inclusion of a second vibrational mode creates an additional fission channel even when its Huang-Rhys factor is relatively small. Addition of more vibrational modes may not enhance the fission per se, but can dramatically affect the interplay between fission dynamics and the dominant vibrational mode.

Dynamics of a One-Dimensional Holstein Polaron with the Hierarchical Equations of Motion Approach

Dynamics of a one-dimensional Holstein molecular crystal model is investigated by making use of the hierarchical equations of motion (HEOM) introduced by Tanimura and Kubo [J. Phys. Soc. Jpn. 1989, 104, 101]. Our extended, numerically exact HEOM approach is capable of treating exciton-phonon coupling in a nonperturbative manner and is applicable to any temperature. It is revealed that strong exciton phonon coupling leads to excitonic localization, while a large exciton transfer integral facilitates exciton transport. Temperature effects on excitonic scattering have also been examined. A proof of concept, our work also serves as a benchmark for future comparisons with other numerical approaches to Holstein polaron dynamics.

Theory of femtosecond coherent double-pump single-molecule spectroscopy: Application to light harvesting complexes.

We develop a first principles theoretical description of femtosecond double-pump single-molecule signals of molecular aggregates. We incorporate all singly excited electronic states and vibrational modes with significant exciton-phonon coupling into a system Hamiltonian and treat the ensuing system dynamics within the Davydov D1 Ansatz. The remaining intra- and inter-molecular vibrational modes are treated as a heat bath and their effect is accounted for through lineshape functions. We apply our theory to simulate single-molecule signals of the light harvesting complex II. The calculated signals exhibit pronounced oscillations of mixed electron-vibrational (vibronic) origin. Their periods decrease with decreasing exciton-phonon coupling.

Optimal Energy Transfer in Light-Harvesting Systems.

Photosynthesis is one of the most essential biological processes in which specialized pigment-protein complexes absorb solar photons, and with a remarkably high efficiency, guide the photo-induced excitation energy toward the reaction center to subsequently trigger its conversion to chemical energy. In this work, we review the principles of optimal energy transfer in various natural and artificial light harvesting systems. We begin by presenting the guiding principles for optimizing the energy transfer efficiency in systems connected to dissipative environments, with particular attention paid to the potential role of quantum coherence in light harvesting systems. We will comment briefly on photo-protective mechanisms in natural systems that ensure optimal functionality under varying ambient conditions. For completeness, we will also present an overview of the charge separation and electron transfer pathways in reaction centers. Finally, recent theoretical and experimental progress on excitation energy transfer, charge separation, and charge transport in artificial light harvesting systems is delineated, with organic solar cells taken as prime examples.

Ground-state properties of sub-Ohmic spin-boson model with simultaneous diagonal and off-diagonal coupling

By employing a variational approach, the density matrix renormalization group (DMRG), the exact diagonalization, and symmetry and mean-field analyses, the ground-state properties of the two-bath spin-boson model with simultaneous diagonal and off-diagonal coupling are systematically studied in the sub-Ohmic regime. A quantum phase transition from a doubly degenerate “localized phase” to the other doubly degenerate “delocalized phase” is uncovered. Via the multi-D1 Ansatz as the variational wave function, transition points are determined accurately, consistent with the results from DMRG and exact diagonalization. An effective spatial dimension deff = 2.37(6) is then estimated, which is found to be compatible with the mean-field prediction. Furthermore, the quantum phase transition is inferred to be of first order for the baths described by a continuous spectral density function. In the single-mode case, however, the transition is softened.

Finite-temperature time-dependent variation with multiple Davydov states

The Dirac-Frenkel time-dependent variational approach with Davydov Ansätze is a sophisticated, yet efficient technique to obtain an accurate solution to many-body Schrödinger equations for energy and charge transfer dynamics in molecular aggregates and light-harvesting complexes. We extend this variational approach to finite temperature dynamics of the spin-boson model by adopting a Monte Carlo importance sampling method. In order to demonstrate the applicability of this approach, we compare calculated real-time quantum dynamics of the spin-boson model with that from numerically exact iterative quasiadiabatic propagator path integral (QUAPI) technique. The comparison shows that our variational approach with the single Davydov Ansätze is in excellent agreement with the QUAPI method at high temperatures, while the two differ at low temperatures. Accuracy in dynamics calculations employing a multitude of Davydov trial states is found to improve substantially over the single Davydov Ansatz, especially at low temperatures. At a moderate computational cost, our variational approach with the multiple Davydov Ansatz is shown to provide accurate spin-boson dynamics over a wide range of temperatures and bath spectral densities.

Polaron dynamics with off-diagonal coupling: beyond the Ehrenfest approximation

Treated traditionally by the Ehrenfest approximation, the dynamics of a one-dimensional molecular crystal model with off-diagonal exciton-phonon coupling is investigated in this work using the Dirac-Frenkel time-dependent variational principle with the multi-D2Ansatz. It is shown that the Ehrenfest method is equivalent to our variational method with the single D2Ansatz, and with the multi-D2Ansatz, the accuracy of our simulated dynamics is significantly enhanced in comparison with the semi-classical Ehrenfest dynamics. The multi-D2Ansatz is able to capture numerically accurate exciton momentum probability and help clarify the relation between the exciton momentum redistribution and the exciton energy relaxation. The results demonstrate that the exciton momentum distributions in the steady state are determined by a combination of the transfer integral and the off-diagonal coupling strength, independent of the excitonic initial conditions. We also probe the effect of the transfer integral and the off-diagonal coupling on exciton transport in both real and reciprocal space representations. Finally, the variational method with importance sampling is employed to investigate temperature effects on exciton transport using the multi-D2Ansatz, and it is demonstrated that the variational approach is valid in both low and high temperature regimes.

Transient dynamics of a one-dimensional Holstein polaron under the influence of an external electric field

Following the Dirac-Frenkel time-dependent variational principle, transient dynamics of a one-dimensional Holstein polaron with diagonal and off-diagonal exciton-phonon coupling in an external electric field is studied by employing the multi-D2 Ansatz, also known as a superposition of the usual Davydov D2 trial states. Resultant polaron dynamics has significantly enhanced accuracy, and is in perfect agreement with that derived from the hierarchy equations of motion method. Starting from an initial broad wave packet, the exciton undergoes typical Bloch oscillations. Adding weak exciton-phonon coupling leads to a broadened exciton wave packet and a reduced current amplitude. Using a narrow wave packet as the initial state, the bare exciton oscillates in a symmetric breathing mode, but the symmetry is easily broken by weak coupling to phonons, resulting in a non-zero exciton current. For both scenarios, temporal periodicity is unchanged by exciton-phonon coupling. In particular, at variance with the case of an infinite linear chain, no steady state is found in a finite-sized ring within the anti-adiabatic regime. For strong diagonal coupling, the multi-D2 Anstaz is found to be highly accurate, and the phonon confinement gives rise to exciton localization and decay of the Bloch oscillations.

Dynamics of the two-spin spin-boson model with a common bath

Dynamics of the two-spin spin-boson model in the presence of Ohmic and sub-Ohmic baths is investigated by employing a multitude of the Davydov D1 trial states, also known as the multi-D1 Ansatz. Its accuracy in dynamics simulations of the two-spin SBM is improved significantly over the single D1 Ansatz, especially in the weak to moderately strong coupling regime. Validity of the multi-D1 Ansatz for various coupling strengths is also systematically examined by making use of the deviation vector which quantifies how faithfully the trial state obeys the Schrödinger equation. The time evolution of population difference and entanglement has been studied for various initial conditions and coupling strengths. Careful comparisons are carried out between our approach and three other methods, i.e., the time-dependent numerical renormalization group (TD-NRG) approach, the Bloch-Redfield theory, and a method based on a variational master equation. For strong coupling, the multi-D1 trial state yields consistent results as the TD-NRG approach in the Ohmic regime while the two disagree in the sub-Ohmic regime, where the multi-D1 trial state is shown to be more accurate. For weak coupling, the multi-D1 trial state agrees with the two master-equation methods in the presence of both Ohmic and sub-Ohmic baths, but shows considerable differences with the TD-NRG approach in the presence of a sub-Ohmic bath, calling into question the validity of the TD-NRG results at long times in the literature.