Modeling many-body quantum systems with strong interactions is one of the core challenges of modern physics. A range of methods has been developed to approach this task, each with its own idiosyncrasies, approximations, and realm of applicability. However, there remain many problems that are intractable for existing methods. In particular, many approaches face a huge computational barrier when modeling large numbers of coupled electrons and ions at finite temperature. Here, we address this shortfall with a new approach to modeling many-body quantum systems. On the basis of the Bohmian trajectory formalism, our new method treats the full particle dynamics with a considerable increase in computational speed. As a result, we are able to perform large-scale simulations of coupled electron-ion systems without using the adiabatic Born-Oppenheimer approximation.
Modeling many-body quantum systems with strong interactions is one of the core challenges of modern physics. A range of methods has been developed to approach this task, each with its own idiosyncrasies, approximations, and realm of applicability. However, there remain many problems that are intractable for existing methods. In particular, many approaches face a huge computational barrier when modeling large numbers of coupled electrons and ions at finite temperature. Here, we address this shortfall with a new approach to modeling many-body quantum systems. On the basis of the Bohmian trajectory formalism, our new method treats the full particle dynamics with a considerable increase in computational speed. As a result, we are able to perform large-scale simulations of coupled electron-ion systems without using the adiabatic Born-Oppenheimer approximation.
The theoretical investigation of nonadiabatic processes is hampered by the complexity of the coupled electron-nuclear dynamics beyond the Born-Oppenheimer approximation. Classically, the simulation of such reactions is limited by the unfavorable scaling of the computational resources as a function of the system size. While quantum computing exhibits proven quantum advantage for the simulation of real-time dynamics, the study of quantum algorithms for the description of nonadiabatic phenomena is still unexplored. In this Letter, we propose a quantum algorithm for the simulation of fast nonadiabatic chemical processes together with an initialization scheme for quantum hardware calculations. In particular, we introduce a first-quantization method for the time evolution of a wave packet on two coupled harmonic potential energy surfaces (Marcus model). In our approach, the computational resources scale polynomially in the system dimensions, opening up new avenues for the study of photophysical processes that are classically intractable.
The method of many-body Green's functions is developed for arbitrary systems of electrons and nuclei starting from the full (beyond Born-Oppenheimer) Hamiltonian of Coulomb interactions and kinetic energies. The theory presented here resolves the problems arising from the translational and rotational invariance of this Hamiltonian that afflict the existing many-body Green's function theories. We derive a coupled set of exact equations for the electronic and nuclear Green's functions and provide a systematic way to approximately compute the properties of arbitrary many-body systems of electrons and nuclei beyond the Born-Oppenheimer approximation. The case of crystalline solids is discussed in detail.
In this work we present a numerical analysis of the so-called ‘interacting quantum trajectory method’ (IQTM), which is based on a formal theory developed by Poirier [Bohmian mechanics without pilot waves, Chem. Phys. 370 (2010), pp. 4–14.], and later extended by Schiff and Poirier [Communication: Quantum mechanics without wavefunctions, J. Chem. Phys. 136 (2012), p. 031102.], for representing the quantum mechanical state of a particle in terms of an ensemble of Bohmian quantum trajectories. This has led to a number of numerical algorithms over the years, mostly for time-independent applications. A numerically robust, exact time-dependent wavepacket propagation scheme would have to await the refinements of Dupuy et al. [L. Dupuy, F. Talotta, F. Agostini, D. Lauvergnat, B. Poirier, and Y. Scribano, Adiabatic and nonadiabatic dynamics with interacting quantum trajectories, J. Chem. Theory Comput. 18 (2022), pp. 6447–6462.], which serves as the starting point for our continued development in the present study. Here we will discuss the numerical protocols comprising IQTM, and present detailed error analysis for several 1D model problems involving simple free particles and oscillators. These studies reveal some key advantages of IQTM; however, they also highlight problems for which the algorithms need further work.
ConspectusSimulating molecular dynamics (MD) within a comprehensive quantum framework has been a long-standing challenge in computational chemistry. An exponential scaling of computational cost renders solving the time dependent Schrödinger equation (TDSE) of a molecular Hamiltonian, including both electronic and nuclear degrees of freedom (DOFs), as well as their couplings, infeasible for more than a few DOFs. In the Born-Oppenheimer (BO), or adiabatic, picture, electronic and nuclear parts of the wave function are decoupled and treated separately. Within this framework, the nuclear wave function evolves along potential energy surfaces (PESs) computed as solutions to the electronic Schrödinger equation parametrized in the nuclear DOFs. This approximation, together with increasingly elaborate numerical approaches to solve the nuclear time dependent Schrödinger equation (TDSE), enabled the treatment of up to a few dozens of degrees of freedom (DOFs). However, for particular applications, such as photochemistry, the BO approximation breaks down. In this regime of non-adiabatic dynamics, solving the full molecular problem including electron-nuclear couplings becomes essential, further increasing the complexity of the numerical solution. Although valuable methods such as multiconfigurational time-dependent Hartree (MCTDH) have been proposed for the solution of the coupled electron-nuclear dynamics, they remain hampered by an exponential scaling in the number of nuclear DOFs and by the difficulty of finding universal variational forms.In this Account, we present a perspective on novel quantum computational algorithms, aiming to alleviate the exponential scaling inherent to the simulation of many-body quantum dynamics. In particular, we focus on the derivation and application of quantum algorithms for adiabatic and non-adiabatic quantum dynamics, which include efficient approaches for the calculation of the BO potential energy surfaces (PESs). Thereafter, we study the time-evolution of a model system consisting of two coupled PESs in first and second quantization. In a first application, we discuss a recently introduced quantum algorithm for the evolution of a wavepacket in first quantization and exploit the potential quantum advantage of mapping its spatial grid representation to logarithmically many qubits. For the second demonstration, we move to the second quantization framework and review the scaling properties of two alternative time-evolution algorithms, namely, a variational quantum algorithm (VQA) (based on the McLachlan variational principle) and conventional Trotter-type evolution (based on a Lie-Trotter-Suzuki formula). Both methods clearly demonstrate the potential of quantum algorithms and their favorable scaling compared to the available classical approaches. However, a clear demonstration of quantum advantage in the context of molecular quantum dynamics may require the implementation of these algorithms in fault-tolerant quantum computers, while their application in near-term, noisy quantum devices is still unclear and deserves further investigation.
Quantum computing introduces a revolutionary approach to solving problems that are computationally beyond the reach of classical systems. Among its most promising applications is the simulation of many-body systems, where strong particle interactions lead to complex quantum behaviors. Traditional computational techniques often fail due to exponential resource requirements, but quantum simulators can efficiently model these interactions through entanglement and superposition. Recent advances in superconducting qubits, trapped ions, and photonic quantum devices have enabled small- scale experimental demonstrations of spin models, Hubbard systems, and correlated electron dynamics. These developments offer new insights into condensed matter physics and quantum phase transitions. As quantum technologies continue to evolve, they hold the potential to transform the study of materials, superconductivity, and other many-body phenomena with remarkable precision and scalability.
We present a new approach to calculate real-time quantum dynamics in complex systems. The formalism is based on the partitioning of a system's environment into "core" and "reservoir" modes with the former to be treated quantum mechanically and the latter classically. The presented method only requires the calculation of the system's reduced density matrix averaged over the quantum core degrees of freedom which is then coupled to a classically evolved reservoir to treat the remaining modes. We demonstrate our approach by applying it to the spin-boson problem using the noninteracting blip approximation to treat the system and core, and Ehrenfest dynamics to treat the reservoir. The resulting hybrid methodology is accurate for both fast and slow baths, since it naturally reduces to its composite methods in their respective regimes of validity. In addition, our combined method is shown to yield good results in intermediate regimes where neither approximation alone is accurate and to perform equally well for both strong and weak system-bath coupling. Our approach therefore provides an accurate and efficient methodology for calculating quantum dynamics in complex systems.