Coupled trajectory mixed quantum classical (CTMQC) dynamics is a rigorous approach to trajectory-based non-adiabatic dynamics, which has recently seen an improvement to energy conservation via the introduction of the CTMQC-E algorithm. Despite this, the method's two key quantities distinguishing it from Ehrenfest dynamics, the modified Born-Oppenheimer momentum and the quantum momentum, require regularisation procedures in certain circumstances. Such procedures in the latter can cause instabilities leading to undesirable effects such as energy drift and spurious population transfer, which is expected to become increasingly prevalent the larger the system as such events would happen more frequently. We propose a further modification to CTMQC-E which includes a redefinition of the quantum momentum, CTMQC-EDI (Double Intercept), such that it has no formal divergences. We then show for Tully models I-IV that the algorithm has greatly improved total energy conservation and negligible spurious population transfer at all times, in particular in regions of strong non-adiabatic coupling. CTMQC-EDI therefore shows promise as a numerically robust non-adiabatic dynamics technique that accounts for decoherence from first principles and that is scalable to large molecular systems and materials.
Coupled trajectory mixed quantum classical (CTMQC) dynamics is a rigorous approach to trajectory-based non-adiabatic dynamics, which has recently seen an improvement to energy conservation via the introduction of the CTMQC-E algorithm. Despite this, the method's two key quantities distinguishing it from Ehrenfest dynamics, the modified Born-Oppenheimer momentum and the quantum momentum, require regularisation procedures in certain circumstances. Such procedures in the latter can cause instabilities leading to undesirable effects such as energy drift and spurious population transfer, which is expected to become increasingly prevalent the larger the system as such events would happen more frequently. We propose a further modification to CTMQC-E which includes a redefinition of the quantum momentum, CTMQC-EDI (Double Intercept), such that it has no formal divergences. We then show for Tully models I-IV that the algorithm has greatly improved total energy conservation and negligible spurious population transfer at all times, in particular in regions of strong non-adiabatic coupling. CTMQC-EDI therefore shows promise as a numerically robust non-adiabatic dynamics technique that accounts for decoherence from first principles and that is scalable to large molecular systems and materials.
We analyze a mixed quantum-classical algorithm recently derived from the exact factorization equations [Min, Agostini, Gross, PRL 115, 073001 (2015)] to show the role of the different terms in the algorithm in bringing about decoherence and wavepacket branching. The algorithm has the structure of Ehrenfest equations plus a "coupled-trajectory" term for both the electronic and nuclear equations, and we analyze the relative roles played by the different nonadiabatic terms in these equations, including how they are computed in practice. In particular, we show that while the coupled-trajectory term in the electronic equation is essential in yielding accurate dynamics, that in the nuclear equation has a much smaller effect. A decoherence time is extracted from the electronic equations and compared with that of augmented fewest-switches surface-hopping. We revisit a series of nonadiabatic Tully model systems to illustrate our analysis.
We analyze a mixed quantum-classical algorithm recently derived from the exact factorization equations [Min, Agostini, Gross, PRL 115, 073001 (2015)] to show the role of the different terms in the algorithm in bringing about decoherence and wavepacket branching. The algorithm has the structure of Ehrenfest equations plus a "coupled-trajectory" term for both the electronic and nuclear equations, and we analyze the relative roles played by the different nonadiabatic terms in these equations, including how they are computed in practice. In particular, we show that while the coupled-trajectory term in the electronic equation is essential in yielding accurate dynamics, that in the nuclear equation has a much smaller effect. A decoherence time is extracted from the electronic equations and compared with that of augmented fewest-switches surface-hopping. We revisit a series of nonadiabatic Tully model systems to illustrate our analysis.
The method, described in the paper, incorporates ideas of the Bohmian formulation of quantum mechanics into a trajectory approach designed for systems comprised of light (quantum) and heavy (nearly classical) particles. The goal is to treat all degrees of freedom on the same footing but, for practical reasons, to include the quantum-mechanical effects into dynamics for just the light particles. The total wavefunction, taken as a product of slow and fast components representing the heavy and light particles respectively is represented in terms of trajectories. The fast component depends on the trajectory-guided configurations of the heavy particles. The initial conditions for the trajectories are defined by the full-dimensional initial wavefunction according to Bohm's prescription of associating the trajectory momentum with the gradient of the wavefunction phase. The time-dependent Schrödinger equation is solved for the quantum degrees of freedom using an ensemble of the quantum trajectories (or other time-propagation method) for each guiding trajectory of the classical subspace. Neglect of the quantum force acting on the heavy particles allows decoupling of the quantum calculations associated with different guiding trajectories. The approach allows reconstruction of the full-dimensional wavefunction and conserves the total wavefunction energy as illustrated for the vibrationally non-adiabatic model.
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.
The dynamics of two coupled oscillators can become quite complex if anharmonic potential energy functions are employed. This type of system therefore represents a good model for an investigation of the performance of mixed quantum-classical methods. In this work, the motion of two coupled particles with a mass ratio of one to ten is studied with three different mixed quantum-classical methods in the presence of anharmonic potential terms for a comparison with exact quantum mechanical calculations. The mixed quantum-classical approaches include the multitrajectory Ehrenfest, the mixed quantum-classical Bohmian (MQCB), and the so-called coupled Schrödinger equations (CSE) formalisms. The analysis shows that while the description of a weakly anharmonic system by the Ehrenfest and MQCB schemes is accurate if proper sampling techniques are applied, both approximations break down rapidly if the anharmonic terms are increased. The performance of the simple CSE prescription, which corresponds to a reduction of the full two-dimensional wave function to two one-dimensional wave functions representing two quantum oscillators coupled via the potential energy in a classical fashion, decreases if the width of the initial wave packet is enlarged. The dependence of the CSE method on the diffuseness of the initial wave packet is therefore opposite to that of the MQCB method, which is more accurate for wide wave packets. Overall, the multitrajectory Ehrenfest ansatz is found to be most successful in reproducing the exact quantum results.
We present a novel quantum-classical approach to nonadiabatic dynamics, deduced from the coupled electronic and nuclear equations in the framework of the exact factorization of the electron-nuclear wave function. The method is based on the quasiclassical interpretation of the nuclear wave function, whose phase is related to the classical momentum and whose density is represented in terms of classical trajectories. In this approximation, electronic decoherence is naturally induced as an effect of the coupling to the nuclei and correctly reproduces the expected quantum behavior. Moreover, the splitting of the nuclear wave packet is captured as a consequence of the correct approximation of the time-dependent potential of the theory. This new approach offers a clear improvement over Ehrenfest-like dynamics. The theoretical derivation presented in this Letter is supported by numerical results that are compared to quantum mechanical calculations.