Accelerating quantum dynamical simulations with quantum computing has received considerable attention but remains a significant challenge. In variational quantum algorithms for quantum dynamics, designing an expressive and shallow-depth parametrized quantum circuit (PQC) is a key difficulty. Here, we propose a multiset variational quantum dynamics algorithm (MS-VQD) tailored for nonadiabatic dynamics involving multiple electronic states. The MS-VQD employs multiple PQCs to represent the electronic-nuclear coupled wave function, with each circuit adapting to the motion of the nuclear wavepacket on a specific potential energy surface. By simulating excitation energy transfer dynamics in molecular aggregates described by the Frenkel-Holstein model, we demonstrate that the MS-VQD achieves the same accuracy as the traditional VQD while requiring significantly shallower PQCs. Notably, its advantage increases with the number of electronic states, making it suitable for simulating nonadiabatic quantum dynamics in complex molecular systems.
Accelerating quantum dynamical simulations with quantum computing has received considerable attention but remains a significant challenge. In variational quantum algorithms for quantum dynamics, designing an expressive and shallow-depth parametrized quantum circuit (PQC) is a key difficulty. Here, we propose a multiset variational quantum dynamics algorithm (MS-VQD) tailored for nonadiabatic dynamics involving multiple electronic states. The MS-VQD employs multiple PQCs to represent the electronic-nuclear coupled wave function, with each circuit adapting to the motion of the nuclear wavepacket on a specific potential energy surface. By simulating excitation energy transfer dynamics in molecular aggregates described by the Frenkel-Holstein model, we demonstrate that the MS-VQD achieves the same accuracy as the traditional VQD while requiring significantly shallower PQCs. Notably, its advantage increases with the number of electronic states, making it suitable for simulating nonadiabatic quantum dynamics in complex molecular systems.
Dynamical simulation of quantum many-body systems is a central task in quantum chemistry and requires efficient wavefunction representations. Although tensor-network and neural-network quantum states have achieved considerable success in ground-state calculations, entanglement growth hinders their application to quantum dynamics. Quantum computing may offer a route to quantum advantage, but algorithms such as Trotterization generally require fault-tolerant quantum computers, whereas near-term hardware supports only limited circuit sizes. Variational quantum simulation (VQS) instead represents the evolving state with a single shallow, fixed-size parameterized quantum circuit (PQC). Here, we introduce nonorthogonal variational quantum simulation (NOVQS), which applies linear combinations of parameterized quantum states to real- and imaginary-time evolutions. We design a shallow, hardware-friendly ansatz tailored to second-quantized electronic-structure Hamiltonians, together with resource-efficient protocols for measuring the matrices and vectors in the parameter equations of motion. Error analysis and resource estimation are also provided. Numerical simulations of hydrogen chains and the nitrogen molecule demonstrate that a collection of shallow, or even single-layer, PQCs can match or outperform a much deeper PQC in VQS. In particular, we identify a trade-off between circuit number and depth. NOVQS therefore provides a flexible route to enhancing wavefunction expressivity under circuit-depth constraints, making it a promising framework for quantum dynamics simulations on near-term quantum processors.
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.
Investigating ultrafast electron-nuclear dynamics in molecules requires robust theoretical methods that offer predictive capabilities beyond interpretation purposes. However, the vast number of electronic and nuclear degrees of freedom often prevents an accurate description of chemical systems on classical computers. To address these challenges, hybrid approaches combining quantum and classical computing can leverage the strengths of both to enhance the computational efficiency and fully exploit hardware capabilities. In this work, we perform excited-state molecular dynamics simulations of the formaldimine molecule by combining the surface hopping nonadiabatic molecular dynamics technique with quantum algorithms. We replace conventional electronic energy and overlap calculations with the variational quantum deflation (VQD) algorithm and the Hadamard test, respectively, employing a configuration state function (CSF) spin-restricted ansatz. Finally, we evaluate the accuracy and feasibility of three different gradient calculation approaches integrated with the quantum computing components, demonstrating the advantages of an analytical gradient method over finite difference approaches for hybrid quantum-classical integration.
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.
Abstract In this Perspective, we discuss how quantum computers may advance the simulation of nonadiabatic molecular dynamics, a framework central to describing excited-state processes in photochemistry, biology, and materials science. Classical approaches span from exponentially scaling full quantum dynamics to more approximate mixed quantum–classical techniques such as surface hopping and Ehrenfest dynamics. Hybrid quantum–classical algorithms – particularly those based on the variational quantum eigensolver – offer a transformative alternative by providing access to the key electronic properties needed to drive nonadiabatic molecular dynamics simulations, including energies, gradients, and nonadiabatic couplings. We examine recent proof-of-principle quantum simulations of reduced model systems which, despite being restricted to small molecules and limited active spaces due to constraints of qubit number and device noise, already showcase the potential of quantum devices to capture phenomena such as conical intersections and ultrafast relaxation. Although practical applications are not yet feasible in the present noisy intermediate-scale quantum era, these efforts underline the conceptual and methodological advances of quantum algorithms paving the way for large-scale quantum simulations of nonadiabatic processes. Framed within the 2025 International Year of Quantum Science and Technology, such progress exemplifies how quantum computing may open new horizons for chemistry and beyond.
Abstract Variational Quantum Eigensolver (VQE) and its extension, Variational Quantum Deflation (VQD), have emerged as promising algorithms for computing ground and excited state energy eigenvalues of Hamiltonians, particularly in quantum physics, chemistry, and materials science. Recently, the VQE and VQD algorithms were implemented to calculate the electronic band structure described by the tight-binding Hamiltonian. Despite their success for simple models, scalability remains a challenge due to the need for a large number of circuit executions, limited qubit coherence times, and low gate fidelities. In this work, we implemented the VQE and VQD algorithms for tight-binding models of diamond and zincblende crystal structures (Sn, C, Si, Ge, AlP, AlAs, AlSb, and GaP) using a quantum computer simulator. We investigate the scalability challenges, focusing on the significant overhead caused by a large number of quantum circuit executions required to find the energy eigenvalues. Our analysis highlights the substantial amount of both measurement and optimization overheads in VQD, particularly as the number of excited states increases. Based on our findings, we critically discuss the implications of implementing VQD on current noisy intermediate-scale quantum (NISQ) devices.