Driven by the impetus to simulate quantum dynamics in photosynthetic complexes or even larger molecular aggregates, we have established a hierarchy of forward-backward stochastic Schrödinger equation in the light of stochastic unravelling of the symmetric part of the influence functional in the path-integral formalism of reduced density operator. The method is numerically exact and is suited for Debye-Drude spectral density, Ohmic spectral density with an algebraic or exponential cutoff, as well as discrete vibrational modes. The power of this method is verified by performing the calculations of time-dependent population differences in the valuable spin-boson model from zero to high temperatures. By simulating excitation energy transfer dynamics of the realistic full FMO trimer, some important features are revealed.

A time-dependent wavepacket diffusion method is used to investigate the effects of charge transfer (CT) states, singlet exciton and multiexciton migrations on singlet fission (SF) dynamics in organic aggregates. The results reveal that the incorporation of CT states can result in a different SF dynamics from the direct interaction between singlet exciton and multiexciton, and an obvious SF interference is also observed between the direct channel and the indirect channel mediated by CT states. In the case of direct interaction, although the fast population transfer of singlet exciton in monomers, by the increase of exciton-exciton interaction, can accelerate the SF process, the spatial coherence alternatively has a counter-productive effect, and their competition leads to an optimal exciton-exciton interaction at which SF has a maximal rate. This trade-off relationship in SF dynamics is further analysed from different perspectives, specifically in spatial and energy representations, is also confirmed through the indication that static energy disorders can speed up SF process by destructing the coherence.Meanwhile, it is found that the couplings among multiexciton states decrease SF rates by the multiexciton coherence and backward conversion from multiexciton to singlet exciton states.

The vibrationally resolved absorption spectra of zinc phthalocyanine (ZnPc) aggregates (up to 70 monomers) are explored using the non-Markovian stochastic Schrödinger equation. Various types of local excitations, charge-transfer (CT) excitations, and exciton–phonon couplings are explicitly included in a comprehensive model Hamiltonian, which is parameterized by first-principles calculations. The absorption spectral simulations clarify that the two absorption bands in the Q-band region observed in experiments can be assigned to the contribution from the CT-mediated interactions, rather than the mixtures of different-type aggregates, as prevailingly assumed. Furthermore, the relative intensities of the two bands are found to be closely related to the intermolecular distance and molecular number in a ZnPc aggregate. From the investigation of the decoherence process after optical excitation, it is found that CT states can induce coherence regeneration as the time scale of charge separation is much faster than that of the vibration-induced decoherence. However, they would instead boost the decoherence process as the two time scales become comparable. The two different effects of CT states may suggest a novel way to regulate the decoherence process in excitation energy relaxation.

No abstract

A number of non‐Markovian stochastic Schrödinger equations, ranging from the numerically exact hierarchical form toward a series of perturbative expressions sequentially presented in an ascending degrees of approximations are revisited in this short review, aiming at providing a systematic framework which is capable to connect different kinds of the wavefunction‐based approaches for an open system coupled to the harmonic bath. One can optimistically expect the extensive future applications of those non‐Markovian stochastic Schrödinger equations in large‐scale realistic complex systems, benefiting from their favorable scaling with respect to the system size, the stochastic nature which is extremely suitable for parallel computing, and many other distinctive advantages. In addition, we have presented a few examples showing the excitation energy transfer in the Fenna‐Matthews‐Olson complex, a quantitative measure of decoherence timescale of hot exciton, and the study of quantum interference effects upon the singlet fission processes in organic materials, since a deep understanding of both mechanisms is very important to explore the underlying microscopic processes and to provide novel design principles for highly efficient organic photovoltaics. This article is categorized under: Theoretical and Physical Chemistry > Reaction Dynamics and Kinetics Structure and Mechanism > Computational Materials Science

The unified coherent-to-diffusive energy relaxation of hot exciton in organic aggregates or polymers, which still remains largely unclear and is also a great challenge theoretically, is investigated from a time-dependent wavepacket diffusive approach. The results demonstrate that in the multiple time scale energy relaxation dynamics, the fast relaxation time essentially corresponds to the dephasing time of excitonic coherence motion, whereas the slow time is related to a hopping migration, and a suggested kinetic model successfully connects these two processes. The dependencies of those times on the initial energy and delocalization of exciton wavepacket as well as exciton-phonon interactions are further analyzed. The proposed method together with quantum chemistry calculations has explained an experimental observation of hot exciton energy relaxation in the low-bandgap copolymer PBDTTPD.

Based on the stochastic unravelling of the reduced density operator in the Feynman path integral formalism for an open quantum system in touch with harmonic environments, a new non-Markovian stochastic Schrödinger equation (NMSSE) has been established that allows for the systematic perturbation expansion in the system-bath coupling to arbitrary order. This NMSSE can be transformed in a facile manner into the other two NMSSEs, i.e., non-Markovian quantum state diffusion and time-dependent wavepacket diffusion method. Benchmarked by numerically exact results, we have conducted a comparative study of the proposed method in its lowest order approximation, with perturbative quantum master equations in the symmetric spin-boson model and the realistic Fenna-Matthews-Olson complex. It is found that our method outperforms the second-order time-convolutionless quantum master equation in the whole parameter regime and even far better than the fourth-order in the slow bath and high temperature cases. Besides, the method is applicable on an equal footing for any kind of spectral density function and is expected to be a powerful tool to explore the quantum dynamics of large-scale systems, benefiting from the wavefunction framework and the time-local appearance within a single stochastic trajectory.

Understanding current-induced bond rupture in single-molecule junctions is both of fundamental interest and a prerequisite for the design of molecular junctions, which are stable at higher-bias voltages. In this work, we use a fully quantum mechanical method based on the hierarchical quantum master equation approach to analyze the dissociation mechanisms in molecular junctions. Considering a wide range of transport regimes, from off-resonant to resonant, non-adiabatic to adiabatic transport, and weak to strong vibronic coupling, our systematic study identifies three dissociation mechanisms. In the weak and intermediate vibronic coupling regime, the dominant dissociation mechanism is stepwise vibrational ladder climbing. For strong vibronic coupling, dissociation is induced via multi-quantum vibrational excitations triggered either by a single electronic transition at high bias voltages or by multiple electronic transitions at low biases. Furthermore, the influence of vibrational relaxation on the dissociation dynamics is analyzed and strategies for improving the stability of molecular junctions are discussed.