The accurate solution of dissipative quantum dynamics plays an important role in the simulation of open quantum systems. Here, we propose a machine learning-based universal solver for the hierarchical equations of motion, one of the most widely used approaches which takes into account non-Markovian effects and nonperturbative system-environment interactions in a numerically exact manner. We develop a neural quantum propagator model by utilizing the neural network architecture, which avoids time-consuming iterations and can be used to evolve any initial quantum state for arbitrarily long times. To demonstrate the efficacy of our model, we apply it to the simulation of population dynamics and linear and two-dimensional spectra of the Fenna-Matthews-Olson complex.
The accurate solution of dissipative quantum dynamics plays an important role in the simulation of open quantum systems. Here, we propose a machine learning-based universal solver for the hierarchical equations of motion, one of the most widely used approaches which takes into account non-Markovian effects and nonperturbative system-environment interactions in a numerically exact manner. We develop a neural quantum propagator model by utilizing the neural network architecture, which avoids time-consuming iterations and can be used to evolve any initial quantum state for arbitrarily long times. To demonstrate the efficacy of our model, we apply it to the simulation of population dynamics and linear and two-dimensional spectra of the Fenna-Matthews-Olson complex.
The accurate (or even approximate) solution of the equations that govern the dynamics of dissipative quantum systems remains a challenging task in quantum science. While several algorithms have been designed to solve those equations with different degrees of flexibility, they rely mainly on highly expensive iterative schemes. Most recently, deep neural networks have been used for quantum dynamics, but current architectures are highly dependent on the physics of the particular system and usually limited to population dynamics. Here we introduce an artificial-intelligence-based surrogate model that solves dissipative quantum dynamics by parametrizing quantum propagators as Fourier neural operators, which we train using both data set and physics-informed loss functions. Compared with conventional algorithms, our quantum neural propagator avoids time-consuming iterations and provides a universal superoperator that can be used to evolve any initial quantum state for arbitrarily long times. To illustrate the wide applicability of the approach, we employ our quantum neural propagator to compute the population dynamics and time-correlation functions of the Fenna-Matthews-Olson complex.
Steady states play pivotal roles in many equilibrium and non-equilibrium open system studies. Their accurate evaluations call for exact theories with rigorous treatment of system-bath interactions. Therein, the hierarchical equations-of-motion (HEOM) formalism is a nonperturbative and non-Markovian quantum dissipation theory, which can faithfully describe the dissipative dynamics and nonlinear response of open systems. Nevertheless, solving the steady states of open quantum systems via HEOM is often a challenging task, due to the vast number of dynamical quantities involved. In this work, we propose a self-consistent iteration approach that quickly solves the HEOM steady states. We demonstrate its high efficiency with accurate and fast evaluations of low-temperature thermal equilibrium of a model Fenna-Matthews-Olson pigment-protein complex. Numerically exact evaluation of thermal equilibrium Rényi entropies and stationary emission line shapes is presented with detailed discussion.
Quantum computing offers promising new avenues for tackling the long-standing challenge of simulating the quantum dynamics of complex chemical systems, particularly open quantum systems coupled to external baths. However, simulating such nonunitary dynamics on quantum computers is challenging since quantum circuits are specifically designed to carry out unitary transformations. Furthermore, chemical systems are often strongly coupled to the surrounding environment, rendering the dynamics non-Markovian and beyond the scope of Markovian quantum master equations like Lindblad or Redfield. In this work, we introduce a quantum algorithm designed to simulate non-Markovian dynamics of open quantum systems. Our approach enables the implementation of arbitrary quantum master equations on noisy intermediate-scale quantum (NISQ) computers. We illustrate the method as applied in conjunction with the numerically exact hierarchical equations of motion (HEOM) method. The effectiveness of the resulting quantum HEOM algorithm is demonstrated as applied to simulations of the non-Lindbladian electronic energy and charge transfer dynamics in models of the carotenoid-porphyrin-C 60 molecular triad dissolved in tetrahydrofuran and the Fenna–Matthews–Olson complex.
The Non-Markovian Stochastic Schrodinger Equation (NMSSE) offers a promising approach for open quantum simulations, especially in large systems, owing to its low scaling complexity and suitability for parallel computing. However, its application at low temperatures faces significant convergence challenges. While short-time evolution converges quickly, long-time evolution requires a much larger number of stochastic trajectories, leading to high computational costs. To this end,we propose a scheme that combines neural network techniques with simulations of the non-Markovian stochastic Schrodinger equation. By integrating convolutional neural networks (CNNs) and long short-term memory recurrent neural networks (LSTMs),along with the iterative attentional feature fusion (iAFF) technique, this approach significantly reduces the number of trajectories required for long-time simulations, particularly at low temperatures, thereby substantially lowering computational costs and improving convergence. To demonstrate our approach, we investigated the dynamics of the spin-boson model and the Fenna-Matthews-Olson (FMO) complex across a range of parameter variations.
Describing the dynamics of strong-laser driven open quantum systems is a very challenging task that requires the solution of highly involved equations of motion. While machine learning techniques are being applied with some success to simulate the time evolution of individual quantum states, their use to approximate time-dependent operators (that can evolve various states) remains largely unexplored. In this work, we develop driven (NQP), a universal neural network framework that solves driven-dissipative quantum dynamics by approximating propagators rather than wave functions or density matrices. NQP can handle arbitrary initial quantum states, adapt to various external fields, and simulate long-time dynamics, even when trained on far shorter time windows. Furthermore, by appropriately configuring the external fields, our trained NQP can be transferred to systems governed by different Hamiltonians. We demonstrate the effectiveness of our approach by studying the spin-boson and the three-state transition Gamma models.
The non-Markovian stochastic Schrödinger equation (NMSSE) offers a promising approach for open quantum simulations owing to its low scaling complexity and suitability for parallel computing. However, its application at low temperatures faces significant convergence challenges. While short-time evolution converges quickly, long-time evolution requires a much larger number of stochastic trajectories, leading to high computational costs. To this end, we propose a scheme that utilizes neural networks to extract effective information from long-time simulations of the non-Markovian stochastic Schrödinger equation, fine-tunes the networks using converged short-time data, and predicts long-term system evolution, eliminating the influence of stochastic oscillations. By integrating convolutional neural networks (CNNs) and long short-term memory recurrent neural networks (LSTMs), along with the iterative attentional feature fusion (iAFF) technique, this approach significantly reduces the number of trajectories required for long-time simulations, particularly at low temperatures, thus substantially lowering computational costs and improving convergence. To demonstrate our approach, we investigated the dynamics of the spin-boson model and the Fenna-Matthews-Olson (FMO) complex across a range of parameter variations.