The simulation of open quantum dynamics is a critical tool for understanding how the non-classical properties of matter might be functionalised in future devices. However, unlocking the enormous potential of molecular quantum processes is highly challenging due to the very strong and non-Markovian coupling of 'environmental' molecular vibrations to the electronic 'system' degrees of freedom. Here, we present an advanced but general computational strategy that allows tensor network methods to effectively compute the non-perturbative, real-time dynamics of exponentially large vibronic wave functions of real molecules. We demonstrate how ab initio modelling, machine learning and entanglement analysis can enable simulations which provide real-time insight and direct visualisation of dissipative photophysics, and illustrate this with an example based on the ultrafast process known as singlet fission.
The simulation of open quantum dynamics is a critical tool for understanding how the non-classical properties of matter might be functionalised in future devices. However, unlocking the enormous potential of molecular quantum processes is highly challenging due to the very strong and non-Markovian coupling of 'environmental' molecular vibrations to the electronic 'system' degrees of freedom. Here, we present an advanced but general computational strategy that allows tensor network methods to effectively compute the non-perturbative, real-time dynamics of exponentially large vibronic wave functions of real molecules. We demonstrate how ab initio modelling, machine learning and entanglement analysis can enable simulations which provide real-time insight and direct visualisation of dissipative photophysics, and illustrate this with an example based on the ultrafast process known as singlet fission.
We develop a powerful tree tensor network states method that is capable of simulating exciton-phonon quantum dynamics of larger molecular complexes and open quantum systems with multiple bosonic environments. We interface this method with ab initio density functional theory to study singlet exciton fission (SF) in a pentacene dimer. With access to the full vibronic many-body wave function, we track and assign the contributions of different symmetry classes of vibrations to SF and derive energy surfaces, enabling us to dissect, understand, and describe the strongly coupled electronic and vibrational dynamics, relaxation, and reduced state cooling. This directly exposes the rich possibilities of exploiting the functional interplay of molecular symmetry, electronic structure and vibrational dynamics in SF material design. The described method can be similarly applied to other complex (bio-) molecular systems, characterised by a rich manifold of electronic states and vibronic coupling driving non-adiabatic dynamics.
Understanding decoherence and dissipation in quantum systems is essential for exploiting quantum effects in emerging technological applications. A proper understanding requires accurate modelling of the system-environment interaction which can only be achieved through treating the system and environment on an equal footing. Simulating the entire system+environment wave-function is a challenging task due to the exponential scaling which is characteristic of quantum mechanics. By using many-body methods based on the powerful concept of a tensor network, one can, nonetheless, obtain numerically exact quantum dynamics for systems with hundreds of degrees of freedom. In addition, a recent innovation has shown that this many-body approach can be extended to finite temperatures without a significant increase in computational cost or algorithmic complexity. This has allowed us to directly compare quantum dynamical simulations with experimental results. Specifically, simulations of the linear vibronic coupling Hamiltonian, parameterised using abinitio electronic structure methods, were used to compute the absorption spectrum of the dye Methylene blue in water. This study has elucidated certain unexplained features of the absorption and shown the how non-adiabatic dynamics resulting from a conical intersection can lead to intensity borrowing between bright and dark excited states. Several methodological advances are also made, including an improvement to the underlying time evolution method TDVP which brings with it an order of magnitude speed-up compared to existing methods.
Light-harvesting in nature is known to work differently than conventional man-made solar cells. Recent studies found electronic excitations, delocalised over several chromophores, and a soft, vibrating structural environment to be key schemes that might protect and direct energy transfer yielding increased harvest efficiencies even under adversary conditions. Unfortunately, testing realistic models of noise assisted transport at the quantum level is challenging due to the intractable size of the environmental wave function. I developed a powerful tree tensor network states (TTNS) method that finds an optimally compressed explicit representation of the combined electronic and vibrational quantum state. With TTNS it is possible to simulate exciton-phonon quantum dynamics from small molecules to larger complexes, modelled as an open quantum system with multiple bosonic environments. After benchmarking the method on the minimal spin-boson model by reproducing ground state properties and dynamics that have been reported using other methods, the vibrational quantum state is harnessed to investigate environmental dynamics and its correlation with the spin system. To enable simulations of realistic non-Born-Oppenheimer molecular quantum dynamics, a clustering algorithm and novel entanglement renormalisation tensors are employed to interface TTNS with ab initio density functional theory (DFT). A thereby generated model of a pentacene dimer containing 252 vibrational normal modes was simulated with TTNS reproducing exciton dynamics in agreement with experimental results. Based on the environmental state, the (potential) energy surfaces, underlying the observed singlet fission dynamics, were calculated yielding unprecedented insight into the super-exchange mediated avoided crossing mechanism that produces ultrafast and high yield singlet fission. This combination of DFT and TTNS is a step towards large scale material exploration that can accurately predict excited states properties and dynamics. Furthermore, application to biomolecular systems, such as photosynthetic complexes, may give valuable insights into novel environmental engineering principles for the design of artificial light-harvesting systems.
This paper explores the application of tensor network methods for simulating molecular dynamics on quantum computers. Traditional molecular dynamics simulations are computationally expensive, particularly for complex systems exhibiting strong quantum correlations. Tensor networks provide a powerful framework for representing and manipulating high-dimensional quantum states, potentially offering a pathway to efficiently simulate molecular dynamics by directly leveraging quantum computational resources. This work details the theoretical framework, outlining how tensor network techniques, specifically Matrix Product States (MPS) and Tree Tensor Networks (TTN), can be adapted for this purpose. We present a simplified model of a diatomic molecule to illustrate the core concept: the construction of a tensor network representation of the molecular wavefunction, followed by the application of time evolution operators to simulate dynamics. The focus is on demonstrating the potential of this approach rather than a fully optimized simulation, highlighting key challenges and future research directions. The efficiency gains compared to classical methods are discussed, acknowledging the significant technological hurdles involved in realizing this approach on near-term quantum devices. We investigate the scalability of the method, considering the limitations imposed by the bond dimension of the tensor network. The ultimate goal is to demonstrate how tensor networks can provide a viable strategy for tackling computationally intractable molecular dynamics problems, paving the way for a new era of quantum-enhanced molecular simulations.
The difficulty to simulate the dynamics of open quantum systems resides in their coupling to many-body reservoirs with exponentially large Hilbert space. Applying a tensor network approach in the time domain, we demonstrate that effective small reservoirs can be defined and used for modeling open quantum dynamics. The key element of our technique is the timeline reservoir network (TRN), which contains all the information on the reservoir's characteristics, in particular, the memory effects timescale. The TRN has a one-dimensional tensor network structure, which can be effectively approximated in full analogy with the matrix product approximation of spin-chain states. We derive the sufficient bond dimension in the approximated TRN with a reduced set of physical parameters: coupling strength, reservoir correlation time, minimal timescale, and the system's number of degrees of freedom interacting with the environment. The bond dimension can be viewed as a measure of the open dynamics complexity. Simulation is based on the semigroup dynamics of the system and effective reservoir of finite dimension. We provide an illustrative example showing the scope for new numerical and machine learning-based methods for open quantum systems.
Non-Markovian open quantum networks are among the most challenging problems in quantum theory, as they suffer from both, the exponentially scaling many-body state space as well as long memory times when the sites are strongly coupled to structured environments. Recent efforts to address them by combining tensor networks with real-time path integrals, such as TEMPO and process tensors, were however limited to systems with rather weak system-environment couplings because strong couplings lead to large bond dimensions. Here, we devise a tensor network scheme (CTEMPO) resulting in particularly small bond dimensions for single-site problems. We then extend this to a practical scheme to simulate non-Markovian quantum networks by time-evolving a corresponding comb tensor network. Benchmarking the single-site algorithm on the spin-boson model we find orders-of-magnitude improvement over established methods. The multi-site algorithm is shown to solve the notoriously difficult problem of the 7-site Fenna-Matthews-Olson photosynthetic complex, whose spectral density contains 62 peaks. We then apply our approach to explore the scaling behavior of the superradiance of up to N=20 quantum dots individually coupled to super-Ohmic phonon baths. Our approach provides a turnkey solution to non-Markovian quantum networks and is compatible with more general tensor network topologies.