Over the last few decades, Born-Oppenheimer (BO) and Born-Oppenheimer-Huang (BOH) treatments for electron-nuclear equation of motion have been considered as the cornerstone for qualitative as well as quantitative exploration and understanding of chemical and physical processes. The BO approximation due to the large mass difference between electrons and nuclei as well as the energy separation among the electronic states could able to predict some of the processes (chemical and physical) satisfactorily at lower energy regime. Such approximation (BO) quite often fails for those cases with the involvement of excited electronic states, in particular, due to the coupling among the states. On the other hand, the BO treatment provides two important quantities, namely, the adiabatic potential energy surfaces (PESs) and the nonadiabatic coupling terms (NACTs) determining the forces, which control the motion of the nuclei in a molecule. The NACTs play a major role in photoexcitation spectra, charge transfer chemical reactions, ferroelectric property and other processes, where significant anomalies are observed between single and multiple PES(s) results. Since the NACTs often exhibit singularity at the close vicinity of degenerate states forming conical intersections (CIs), such coupling terms could break the symmetry of the molecular processes changing the observables qualitatively.
AI Summary
This book-length review addresses how electron–nuclear coupling modifies molecular spectra, reactive scattering, solid-state distortions and molecule–surface collisions when a single-surface Born–Oppenheimer description is insufficient. Its central subject is constructing diabatic Hamiltonians from adiabatic potential-energy surfaces and nonadiabatic derivative couplings, particularly near conical intersections.
The treatment develops the Born–Huang coupled equations, geometric phases, Jahn–Teller symmetry-breaking models and adiabatic-to-diabatic transformations for multiple electronic states. It reviews symbolic and numerical implementations, including the OpenMP-parallelized ADT package, and applications to molecular photoelectron spectra, three-state $\mathrm{H}_3^+$ and $\mathrm{F}+\mathrm{H}_2$ scattering, four-state $\mathrm{HeH}_2^+$ scattering, and transition-metal complexes. These applications use electronic-structure calculations followed by transformed-Hamiltonian nuclear dynamics, with comparisons to available experimental and theoretical spectra, cross-sections and rate constants. The presentation consolidates earlier studies rather than identifying one new quantitative benchmark.
A further component treats finite-temperature molecule–surface scattering through effective mean-field potentials for phonons and electron–hole excitations. Six molecular degrees of freedom plus a dissipation coordinate produce a $(6+1)$-dimensional scattering formulation.
The stated conditions and approximations are important: exact removal of derivative couplings within a chosen electronic subspace requires the curl condition; diabatic Hamiltonians are nonunique; several Jahn–Teller spectral formulas employ adiabatic, harmonic or radial–angular separation approximations; and surface dissipation is treated at mean-field level. Incorporating anharmonic molecule–surface and phonon effects is described as ongoing work.
No sign of AI-generated textOverclaimsUnclear noveltyUnclear writingPossible error
Error found: The abstract over-attributes symmetry breaking to singular derivative couplings, whereas the displayed Jahn–Teller models establish vibronic distortion and splitting without showing that derivative-coupling singularity itself is the cause.
The NO3 discussion incorrectly includes the antisymmetric A2′ component in the symmetrized squares of E′ and E″, thereby admitting an incorrect a2′ linear Jahn–Teller coupling channel.
The introduction incorrectly claims that singular nonadiabatic derivative couplings intrinsically cause energy nonconservation in the adiabatic representation, rather than potentially causing numerical difficulties.
Written by gpt-6.1-sol from the paper's full text. Automated: it can be wrong.
Over the last few decades, Born-Oppenheimer (BO) and Born-Oppenheimer-Huang (BOH) treatments for electron-nuclear equation of motion have been considered as the cornerstone for qualitative as well as quantitative exploration and understanding of chemical and physical processes. The BO approximation due to the large mass difference between electrons and nuclei as well as the energy separation among the electronic states could able to predict some of the processes (chemical and physical) satisfactorily at lower energy regime. Such approximation (BO) quite often fails for those cases with the involvement of excited electronic states, in particular, due to the coupling among the states. On the other hand, the BO treatment provides two important quantities, namely, the adiabatic potential energy surfaces (PESs) and the nonadiabatic coupling terms (NACTs) determining the forces, which control the motion of the nuclei in a molecule. The NACTs play a major role in photoexcitation spectra, charge transfer chemical reactions, ferroelectric property and other processes, where significant anomalies are observed between single and multiple PES(s) results. Since the NACTs often exhibit singularity at the close vicinity of degenerate states forming conical intersections (CIs), such coupling terms could break the symmetry of the molecular processes changing the observables qualitatively.
This book-length review addresses how electron–nuclear coupling modifies molecular spectra, reactive scattering, solid-state distortions and molecule–surface collisions when a single-surface Born–Oppenheimer description is insufficient. Its central subject is constructing diabatic Hamiltonians from adiabatic potential-energy surfaces and nonadiabatic derivative couplings, particularly near conical intersections.
The treatment develops the Born–Huang coupled equations, geometric phases, Jahn–Teller symmetry-breaking models and adiabatic-to-diabatic transformations for multiple electronic states. It reviews symbolic and numerical implementations, including the OpenMP-parallelized ADT package, and applications to molecular photoelectron spectra, three-state H3+ and F+H2 scattering, four-state HeH2+ scattering, and transition-metal complexes. These applications use electronic-structure calculations followed by transformed-Hamiltonian nuclear dynamics, with comparisons to available experimental and theoretical spectra, cross-sections and rate constants. The presentation consolidates earlier studies rather than identifying one new quantitative benchmark.
A further component treats finite-temperature molecule–surface scattering through effective mean-field potentials for phonons and electron–hole excitations. Six molecular degrees of freedom plus a dissipation coordinate produce a (6+1)-dimensional scattering formulation.
The stated conditions and approximations are important: exact removal of derivative couplings within a chosen electronic subspace requires the curl condition; diabatic Hamiltonians are nonunique; several Jahn–Teller spectral formulas employ adiabatic, harmonic or radial–angular separation approximations; and surface dissipation is treated at mean-field level. Incorporating anharmonic molecule–surface and phonon effects is described as ongoing work.
We analyze the entanglement between electronic and nuclear motions in molecular wave functions widely used by theoretical chemists, namely, (i) Born-Oppenheimer factorization in the adiabatic picture, (ii) the transformation into a diabatic picture, (iii) the use of a Born-Huang expansion, and (iv) the eigenfunction of the full molecular Hamiltonian. Our showcase is based on two one-electron one-dimensional molecular Hamiltonians (H2+ and the Shin-Metiu model). We find that within the Born-Oppenheimer approximation, any molecular state (although factorizable) is always entangled, and its entanglement content may be assessed by the variation of the electronic wave function along the different nuclear geometries, with the nuclear wave function indeed playing the role of a tester. The presence of avoided crossings among the adiabatic potential energy curves brings about dramatic changes in the entanglement content of the wave function: sharp avoided crossings favor a diabatic picture (real crossings between potential energy curves), while in broad avoided crossings, the adiabatic picture prevails. The total eigenfunction of the molecular Shin-Metiu Hamiltonian indicates that nuclear densities accommodate well within the diabatic curves for strong adiabatic couplings but within adiabatic curves for weak ones. Consequently, we find that the electron-nuclei entanglement content is a valid witness to unveil strong or weak nonadiabatic couplings in molecules. In terms of entanglement, we also find that the Born-Huang expansion, based on Born-Oppenheimer adiabatic electronic states, does not provide a correct trend of entanglement compared with that of the total molecular eigenfunction, thus indicating a very slow convergence of this expansion.
The role of electron-nuclear entanglement on the validity of the Born-Oppenheimer (BO) approximation is investigated. While nonadiabatic couplings generally lead to entanglement and to a failure of the BO approximation, surprisingly the degree of electron-nuclear entanglement is found to be uncorrelated with the degree of validity of the BO approximation. This is because while the degree of entanglement of BO states is determined by their deviation from the corresponding states in the crude BO approximation, the accuracy of the BO approximation is dictated, instead, by the deviation of the BO states from the exact electron-nuclear states. In fact, in the context of a minimal avoided crossing model, extreme cases are identified where an adequate BO state is seen to be maximally entangled, and where the BO approximation fails but the associated BO state remains approximately unentangled. Further, the BO states are found to not preserve the entanglement properties of the exact electron-nuclear eigenstates, and to be completely unentangled only in the limit in which the BO approximation becomes exact.
The Born-Oppenheimer (BO) and single-surface Born-Huang (BH) approximations replace the full molecular problem by nuclear motion associated with one clamped electronic state, whereas the full BH expansion and exact factorization (EF) are exact representations. We use two analytically transparent models to identify what controls one-surface error. For Fernández's bilinearly coupled oscillators, the Brattsev ground-state bracket was already known; we derive closed-form differences proving strictness for every positive nuclear mass and nonzero stable coupling. Using Hunter's conditional-amplitude construction, we place the closed-form ground-state EF potential alongside BO and single-surface BH. Exact excited-state calculations show that the bracket does not extend through the spectrum. The diagonal BH correction is the quantum-metric coefficient of the retained electronic state divided by twice the nuclear mass. Fernández's sixth-order expansion shows that the next omitted correction contains the same derivative couplings, an additional inverse electronic-energy separation, and a factor depending on the nuclear state. We therefore add a two-state model with a position-dependent gap. Decreasing its minimum separation makes the metric higher and narrower while leaving the total Fubini-Study length equal to pi/2. The exact EF, BO, and BH potentials can differ visibly even when their nuclear ground-state amplitudes are nearly identical. The two models show that one-surface accuracy depends on nuclear mass, electronic-state variation, energy separation from omitted states, and the nuclear wavefunction.
The role of electron-nuclear entanglement on the validity of the Born-Oppenheimer (BO) approximation is investigated. Although nonadiabatic couplings generally lead to entanglement and to a failure of the BO approximation, surprisingly the degree of electron-nuclear entanglement is found to be uncorrelated with the degree of validity of the BO approximation. This is because while the degree of entanglement of BO states is determined by their deviation from the corresponding states in the crude BO approximation, the accuracy of the BO approximation is dictated, instead, by the deviation of the BO states from the exact electron-nuclear states. In fact, in the context of a minimal avoided crossing model, extreme cases are identified where an adequate BO state is seen to be maximally entangled, and where the BO approximation fails but the associated BO state remains approximately unentangled. Further, the BO states are found to not preserve the entanglement properties of the exact electron-nuclear eigenstates, and to be completely unentangled only in the limit in which the BO approximation becomes exact.
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 trajectory calculations for electrons are a useful tool in the field of molecular dynamics, e.g. to understand processes in ultrafast spectroscopy. They have, however, two limitation: On the one hand, such calculations are typically based on the Born-Oppenheimer approximation (BOA) and the electron dynamics for stationary nuclei is considered, thus neglecting quantum effects of the nuclei. On the other hand, even if the quantum nuclear motion would be taken into account, a BOA dynamics on a single potential energy surface would not provide any electron trajectories because the electronic part is treated as a stationary problem. By using the exact factorization method, we overcome these limitations and generalize the theory of electronic quantum trajectories to a fully quantum-mechanical treatment of the nuclei. After reviewing the time-dependent theory of quantum hydrodynamics and quantum trajectories, we show that the nuclei can be viewed as a quantum clock for the electronic motion and we develop a fully quantum-mechanical clock-dependent version of quantum hydrodynamics. This theory is used to obtain electronic trajectories for quantum nuclei, as is exemplified for a model system of a proton-coupled electron transfer dynamics. Our work generalizes the concept of quantum trajectories and lays the foundations for the development of trajectory-based simulation methods of electron dynamics beyond the BOA.