Theoretical study of nonadiabatic hydrogen atom scattering dynamics on metal surfaces using the hierarchical equations of motion method

Simulating Non-Markovian Quantum Dynamics on NISQ Computers Using the Hierarchical Equations of Motion

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-C60 molecular triad dissolved in tetrahydrofuran and the Fenna-Matthews-Olson complex.

Nonadiabatic Quantum Dynamics of Molecules Scattering from Metal Surfaces

Nonadiabatic coupling between electrons and molecular motion at metal surfaces leads to energy dissipation and dynamic steering effects during chemical surface dynamics. We present a theoretical approach to the scattering of molecules from metal surfaces that incorporates all nonadiabatic and quantum nuclear effects due to the coupling of the molecular degrees of freedom to the electrons in the metal. This is achieved with the hierarchical equations of motion (HEOM) approach, combined with a matrix product state representation in twin space. The method is applied to the scattering of nitric oxide from Au(111), for which strongly nonadiabatic energy loss during scattering has been experimentally observed, thus presenting a significant theoretical challenge. Since the HEOM approach treats the molecule-surface coupling exactly, it captures the interplay between nonadiabatic and quantum nuclear effects. Finally, the data obtained by the HEOM approach are used as a rigorous benchmark to assess various mixed quantum-classical methods, from which we derive insights into the mechanisms of energy dissipation and the suitable working regimes of each method.

Integrating Self-Initialized Local Thermalizing Lindblad Operators for Variational Quantum Algorithm with Quantum Jump: Implementation and Performance

Quantum computing holds great potential for simulating quantum systems, such as molecular systems, due to its inherent ability to represent and manipulate quantum states. However, simulating nonunitary dynamics of open quantum systems on a quantum computer is still challenging. Here, we present a scheme denoted as VQS-QJ-LTLME, which adopts the trajectory average of quantum jump Monte Carlo wave function variational evolution to evolve the system's density matrix and local thermalizing Lindblad operators to describe the system-environment interactions. This combination allows the quantum circuit to be initialized after a period of evolution, thereby eliminating the accumulation of errors. The VQS-QJ-LTLME algorithm requires only log2(n) qubits for system size n and holds a time complexity of O ( T n 3 ) , making it particularly suitable for operation on noisy intermediate-scale quantum devices. To accelerate the Monte Carlo sampling process in VQS-QJ-LTLME, we introduce an efficient sampling method, named No-evolution sampling (NES). The VQS-QJ-LTLME aided by NES requires simulating only n3 + 1 trajectories on quantum computers, and then makes 106 to 107 samplings on classical computers in a few seconds. To demonstrate the performance of the VQS-QJ-LTLME algorithm, we simulate the dynamics of a two-level spin-boson model and a four-level reduced Fenna-Matthews-Olson system with real superconducting quantum computers and classical simulators. It is shown that the simulations produced by VQS-QJ-LTLME algorithm align closely with those obtained from the purely classical methods.

Hierarchical Equations of Motion for Quantum Chemical Dynamics: Recent Methodology Developments and Applications

ConspectusQuantum effects are critical to understanding many chemical dynamical processes in condensed phases, where interactions between molecules and their environment are usually strong and non-Markovian. In this Account, we review recent progress from our group in development and application of the hierarchical equations of motion (HEOM) method, highlighting its ability to address some challenging problems in quantum chemical dynamics.In the HEOM method, the bath degrees of freedom are represented using effective modes from exponential decomposition of the bath correlation function. Complex spectral densities and low temperature simulations often require a larger number of modes, making the simulations very expensive. Recent advances, such as the barycentric spectral decomposition (BSD) technique, can significantly reduce the number of effective modes, allowing to handle complex spectral densities and enabling simulations at very low temperatures, including near-zero temperature dynamics.Another key improvement in the computational efficiency is the use of tensor network methods like matrix product states and hierarchical tensor networks. These techniques allow for efficient HEOM propagation with thousands of effective modes, crucial for simulating large molecular systems interacting with multiple baths. This combination enables simulations of excitation energy transfer (EET) in systems like the Fenna-Matthews-Olson (FMO) complex and even larger systems with experimentally determined spectral densities.The versatility of the HEOM method is demonstrated through applications to a wide range of chemical dynamics problems. Simulations of EET and related ultrafast spectroscopy are first briefly covered. Applications of the HEOM to quantum tunneling effects in proton transfer reactions are then presented. Early works have studied the non-Kramers dependence of the rate constant as a function of bath friction due to deep tunneling and revealed vibrationally nonadiabatic dynamics within the so-called nontraditional view of proton transfer reactions. A recent work on the large kinetic isotope effects in soybean lipoxygenase also indicated that many quantum correction approximations to classical transition-state theory may fall short in describing deep tunneling effects.Charge transport and separation dynamics in organic semiconductors are another area where the HEOM method has been instrumental. We first demonstrate that the HEOM provides a unified description of both band-like and thermally assisted charge carrier transport in organic materials. The effect of non-nearest neighbor transitions is then investigated by combining generalized master equations with exact memory kernels. The HEOM method also enables simulation of charge separation in organic photovoltaics (OPVs) and reveals how factors such as external electric fields, entropy, and charge delocalization influence the charge separation barrier and dynamics.Moreover, HEOM has been applied to investigate hydrogen atom scattering on the Au(111) surface and vibrational energy relaxation at molecule-metal interfaces. These studies provide deeper insights into how electron-hole pair excitations and temporary charge transfer states influence the nuclear motion, offering a new framework for simulating nonadiabatic dynamics on metal surfaces.In summary, the HEOM method has developed into a robust tool for simulating quantum effects in condensed phases. Future developments in algorithm efficiency and computational power will likely expand its applicability to even more complex systems.

Simulating Chemistry on Bosonic Quantum Devices

Bosonic quantum devices offer a novel approach to realize quantum computations, where the quantum two-level system (qubit) is replaced with the quantum (an)harmonic oscillator (qumode) as the fundamental building block of the quantum simulator. The simulation of chemical structure and dynamics can then be achieved by representing or mapping the system Hamiltonians in terms of bosonic operators. In this Perspective, we review recent progress and future potential of using bosonic quantum devices for addressing a wide range of challenging chemical problems, including the calculation of molecular vibronic spectra, the simulation of gas-phase and solution-phase adiabatic and nonadiabatic chemical dynamics, the efficient solution of molecular graph theory problems, and the calculations of electronic structure.

Artificial-Intelligence-Based Surrogate Solution of Dissipative Quantum Dynamics: Physics-Informed Reconstruction of the Universal Propagator

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.

Electronic friction near metal surface: Incorporating nuclear quantum effect with ring polymer molecular dynamics

The molecular dynamics with electronic friction (MDEF) approach can accurately describe nonadiabatic effects at metal surfaces in the weakly nonadiabatic limit. That being said, the MDEF approach treats nuclear motion classically such that the nuclear quantum effects are completely missing in the approach. To address this limitation, we combine Electronic Friction with Ring Polymer Molecular Dynamics (EF-RPMD). In particular, we apply the averaged electronic friction from the metal surface to the centroid mode of the ring polymer. We benchmark our approach against quantum dynamics to show that EF-RPMD can accurately capture zero-point energy as well as transition dynamics. In addition, we show that EF-RPMD can correctly predict the electronic transfer rate near metal surfaces in the tunneling limit as well as the barrier crossing limit. We expect that our approach will be very useful to study nonadiabatic dynamics near metal surfaces when nuclear quantum effects become essential.

A Linearized Semiclassical Dynamics Study of the Multiquantum Vibrational Relaxation of NO Scattering from a Au(111) Surface

The vibrational relaxation of NO molecules scattering from a Au(111) surface has served as the focus of efforts to understand nonadiabatic energy transfer at metal-molecule interfaces. Experimental measurements and previous theoretical efforts suggest that multiquantal NO vibrational energy relaxation occurs via electron-hole pair excitations in the metal. Here, using a linearized semiclassical approach, we accurately predict the vibrational relaxation of NO from the νi = 3 state for different incident translational energies. We also accurately capture the central role of transient electron transfer from the metal to the molecule in mediating the vibrational relaxation process but fall short of quantitatively predicting the full extent of multiquantum relaxation for high incident vibrational excitations (νi = 16).