We study the dynamics of a molecule’s nuclear wave function near an avoided crossing of two electronic energy levels for one nuclear degree of freedom. We derive the general form of the Schrödinger equation in the n th superadiabatic representation for all . Using these results, we obtain closed formulas for the time development of the component of the wave function in an initially unoccupied energy subspace when a wave packet travels through the transition region. In the optimal superadiabatic representation, which we define, this component builds up monotonically. Finally, we give an explicit formula for the transition wave function away from the avoided crossing, which is in excellent agreement with high-precision numerical calculations.
We study the dynamics of a molecule’s nuclear wave function near an avoided crossing of two electronic energy levels for one nuclear degree of freedom. We derive the general form of the Schrödinger equation in the n th superadiabatic representation for all . Using these results, we obtain closed formulas for the time development of the component of the wave function in an initially unoccupied energy subspace when a wave packet travels through the transition region. In the optimal superadiabatic representation, which we define, this component builds up monotonically. Finally, we give an explicit formula for the transition wave function away from the avoided crossing, which is in excellent agreement with high-precision numerical calculations.
Abstract Avoided crossings of level pairs with opposite slopes can form potential-energy minima for the external degree of freedom of quantum particles, giving rise to metastable states on the avoided crossings (MSACs). Nonadiabatic decay of MSACs is studied by solving the two-component Schrödinger equation in diabatic and adiabatic representations. Non-perturbative lifetime values are found by evaluating wave function flux and scattering phases of time-independent solutions, as well as wave-function decay of time-dependent solutions. The values from these methods generally agree well, validating the utilized approaches. As the adiabaticity parameter, V, of the system is increased by about a factor of ten across the mixed diabatic/adiabatic regime, the MSAC character transitions from marginally to highly stable, with the lifetimes increasing by about ten orders of magnitude. The dependence of MSAC lifetime on the vibrational quantum number, ν, is discussed for several regimes of V. Time-dependent perturbation theory yields lifetimes that deviate by ≲30% from non-perturbative results, over the range of V and ν studied, while a semi-classical model based on Landau–Zener tunneling is up to a factor of twenty off. The results are relevant to numerous atomic and molecular systems with metastable states on intersecting, coupled potential energy curves.
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.
They have analytically tested several quantum transitions in molecular dynamics. In addition, they have introduced applications of theoretical formalism in physics to economics. First, they have completed the semiclassical theory of nonadiabatic transitions in the exponential potential model. They treated a two-state quantum system, which commonly appears as a subsystem in many-state problems. They were concerned with one-dimensional Schroedinger equation, although the type or semiclassical approach which they used can be applied to higher-dimension potential surfaces by means of Fourier transform. In particular, they solved the case of repulsive exponentially growing potentials, that is the forward-back scattering. They re-expressed the nonadiabatic transition matrix, which describes the transition between the two states for one-direction passage, in terms of contour integrals in the complex coordinate plane. The integrands are adiabatic momenta derived from the eigenvalues of potential matrix. Such type of concept has the advantage that it can be applied to other models than just the exponential potential model. Based on this, they have obtained an analytical tool for describing one-transition point class of potential curve systems. Examples are Rosen-Zener-Demkov model and Landau-Zener-Stueckelberg model. The formalism, which they contributed to, is useful in describing physical processes such as the charge transfer, and it is suitable for shaping external fields or modulating their frequencies in order to achieve a desired state or a product with the theoretical 100% probability. Second, they have derived the first exact analytical solution for a diabatically avoided crossing model. The physical system is described with one-dimensional coupled potential well and barrier which do not intersect in the whole coordinate range. On the basis of special mathematical functions they derived the complete quantum mechanical solution for the wave functions and the transmission coefficient. They also analyzed this model semiclassically and compared it with the exact results. This comparison yields interesting implications on the validity and accuracy of semiclassical approximation. In particular, they solved the case of dramatically avoided crossing with exponentially steep potential curves, and obtained the probability of one-way passage through the interaction region. The overall behavior corresponds to the probability of tunneling through the single potential barrier. In addition, there are complete reflection resonances at discrete energies, resulting from the trapping of a wave in the potential well. The formalism of diabatically avoided crossing model is useful in laser control of molecular reactions. This is because the potential energy curves in realistic molecular systems can be shifted by the photon energy; the curves are coupled by the field intensity via the dipole moment. That is why the new topology of potential energy curves, which was so far broadly omitted, becomes practically interesting due to the development in lasers techniques. Third, they treated several two-state quantum systems common in molecular physics. The semiclassical conditions for the complete transmission and the complete reflection in one-dimensional scattering were given analytically. They discovered the possibility of these phenomena in energy ranges in which the complete reflection or the complete transmission are not commonly expected to occur. In particular, they have made use of the broadly shared assumption that the transitions between the states are localized phenomena. This enables them to separate the whole coordinate region into distinct subsets and to analyze the physics in each subset separately. Applying this approach, they obtained a fully algebraic problem, which can be treated e.g. by diagrammatic techniques. Let them also mention that the complete transmission or complete reflection conditions simplify the wave function, and hence the mathematical solution. The formalism which they used while deriving the semiclassical conditions for the resonant phenomena mentioned above can be further generalized. This will probably yield a whole spectrum of tools in the laser control of molecular processes. Fourth, they developed a powerful technique for solving integral equations with singular kernel. One of the physical formulations of scattering is the time-independent approach based on the Lippmann-Schwinger equation for the Green operator. They have developed a detailed technique how to treat the singularity in this equation. In particular, any scattering problem can be divided into two parts. The first one is the geometry of the system and the static interaction of its components. The second one is more dynamical: the singular equation for the Green operator including the variable energy of the scattering. They fix the most costly calculations in the static part. Then they solve the singularity analytically and thus obtain the cross section as a function of energy. The so far common approach is quite opposite: to fix the singularity first and then to recalculate heavily all the static part for any new value of energy. Obviously, their approach saves a considerable amount of computational time; therefore it is a good standpoint for treating large systems. The formalism above was proved to be extremely efficient taking the example of dissociative recombination H2+ + e- → H + H* because this is a diatomic system, quantum chemical data are easily available in the literature and they could demonstrate the superiority of their approach to the other approaches. Fifth, they suggested applications of basic physical concepts to economics. They introduced aggregate quantities in financial markets and interpreted these quantities as potential energy surfaces. They also formalized interactions between cyclic processes in terms of state to state transitions. Summarizing, they treated the exponential potential model of nonadiabatic transitions, introduced a new exactly solvable model of diabatically avoided crossing, discussed the semiclassical conditions for complete transmission and complete reflection from the viewpoint of laser control, contributed to scattering calculations, and introduced some interdisciplinary applications.
Dynamics based on quantum trajectories with approximate quantum potential is generalized to nonadiabatic systems and its semiclassical properties are discussed. The formulation uses the mixed polar-coordinate space representation of a wave function. The polar part describes the overall time evolution of the wave-function components semiclassically using the single-surface approximate quantum potential. The coordinate part represents a complex``population'' amplitude, which in case of localized coupling can be solved for quantum mechanically in an efficient manner. In the high-energy regime this is accomplished by using a small basis determined by the coupling between surfaces. An illustration is given for a typical curve-crossing problem. The energy-resolved probabilities obtained from the time evolution of two wave packets for a wide range of energies are in excellent agreement with exact results for energies above the threshold of the diabatic reaction, including the case of total nonadiabatic transition.
Preserving quantum coherence with the increase of a system's size and complexity is a major challenge. Molecules, with their diverse sizes and complexities and many degrees of freedom, are an excellent platform for studying the transition from quantum to classical behavior. While most quantum-control studies of molecules focus on vibrations and rotations, we focus here on creating a quantum superposition between two nuclear-spin isomers of the same molecule. We present a scheme that exploits an avoided crossing in the spectrum to create strong coupling between two uncoupled nuclear-spin-isomer states, hence creating an isomeric qubit. We model our scheme using a four-level Hamiltonian and explore the coherent dynamics in the different regimes and parameters of our system. Our four-level model and approach can be applied to other systems with a similar energy-level structure.
We derive an exact solution of an explicitly time-dependent multichannel model of quantum mechanical nonadiabatic transitions. Our model corresponds to the case of a single linear diabatic energy level interacting with a band of an arbitrary N states, for which the diabatic energies decay with time according to the Coulomb law. We show that the time-dependent Schrödinger equation for this system can be solved in terms of Meijer functions whose asymptotics at a large time can be compactly written in terms of elementary functions that depend on the roots of an N th order characteristic polynomial. Our model can be considered a generalization of the Demkov–Osherov model. In comparison to the latter, our model allows one to explore the role of curvature of the band levels and diabatic avoided crossings.