Corrections to the Born-oppenheimer Approximation, Part I
Corrections to Born-Oppenheimer approximation for diatomic systems
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Corrections to Born-Oppenheimer approximation for diatomic systems
The conventional loosely bound molecule interpretation of the X(3872) is not compatible with the recent LHCb experimental measurement of the ratio of branching fractions $\mathcal{R}$ = Br(X → ψ′γ)/Br(X → ψγ). We systematically determine the entire tetraquark spectrum for J = 0, 1, 2 and refine the calculation of $\mathcal{R}$ in an improved Born-Oppenheimer description of the X(3872) compact tetraquark. This refinement yields a significantly better agreement with experimental data on $\mathcal{R}$ and on the spectroscopy of the states themselves. Extending the diquark-antidiquark paradigm to encompass tetraquarks that are linear super-positions of open charm singlets and color octets, we discover that these exotic resonances manifest as compact shallow bound states of quarks in color force potentials.
Adiabatic corrections to long-range Born- Oppenheimer potentials
The adiabatic-nuclei approximation of vibrational-rotational excitation of homonuclear diatomic molecules can be simply augmented to describe the vibrational-rotational coupling by including the dependence of the vibrational wave function on j. Appropriate formulas are given, and the theory, is applied to e-H2 excitation, whereby it is shown that deviations from the simple Born-Oppenheimer approximation measured by Wong and Schultz can be explained. More important, it can be seen that the inclusion of the j-dependent centrifugal term is essential for transitions involving high-rotational quantum numbers.
Isotopic compositions of alkanes are typically assumed to be kinetically controlled, but recently is has been proposed that alkanes can isotopically equilibrate for both C and H isotopes during natural gas generation. Evaluation of this requires knowledge of the isotopic equilibrium between alkanes and other common hydrogen and carbon bearing species. Here, in this study, we calculate isotopic equilibria within and between gaseous dihydrogen (H 2 ), water (H 2 O), methane (CH 4 ), ethane (C 2 H 6 ) and propane (C 3 H 8 ), including isotope fractionation among molecules, clumped isotope effects, as well as among sites of propane (i.e., the site-specific isotope effects) from 0°C to 500°C using a path-integral method paired with high-level descriptions of molecular potentials and the diagonal correction to the Born Oppenheimer approximation. While path-integral calculations with high- level CCSD(T) potentials are available for the isotopic equilibria involving methane, the path-integral calculations for ethane and propane have only been performed based on lower-level descriptions of the molecular potentials. We analyze the relative importance of various approximations that are commonly employed when isotopic equilibria are evaluated. We find that clumped isotope effects can be calculated to the same accuracy using computationally inexpensive combination of the Bigeleisen-Mayer-Urey model with the molecular potential from density functional theory. In contrast, fractionation and site preferences of both deuterium and carbon-13 benefit from the use of the higher level CCSD(T) potentials and accounting for anharmonic effects. Additionally, for fractionation and site preference of deuterium corrections to Born-Oppenheimer approximation can also be important.
We propose a superconducting qubit based on engineering the first and second harmonics of the Josephson energy and phase relation 𝐸 𝐽1 cos 𝜑 and 𝐸 𝐽2 cos 2𝜑. By constructing a circuit such that 𝐸 𝐽2 is negative and |𝐸 𝐽1 | ≪ |𝐸 𝐽2 |, we create a periodic potential with two nondegenerate minima. The qubit, which we dub “harmonium,” is formed from the lowest-energy states of each minimum. Bit-flip protection of the qubit arises due to the localization of each qubit state to their respective minima, while phase-flip protection can be understood by considering the circuit within the Born-Oppenheimer approximation. We demonstrate with time-domain simulations that single- and two-qubit gates can be performed in approximately 100 ns. Finally, we compute the qubit coherence times using numerical diagonalization of the complete circuit in conjunction with state-of-the-art noise models. We estimate out-of-manifold heating times on the order of milliseconds, which can be treated as erasure errors using conventional dispersive readout. We estimate pure-dephasing times on the order of many tens of milliseconds, and bit-flip times on the order of seconds.
Current methods for modeling non-adiabatic molecular dynamics face fundamental limitations when treating geometric phase effects: quantum mechanical phenomena where nuclear wavepackets acquire phase shifts when encircling conical intersections. Existing approaches either neglect these effects entirely or rely on potential energy surfaces arising from the Born-Oppenheimer approximation, which introduce artificial singularities and can overestimate geometric phase contributions. This project developed a new theoretical framework based on exact factorization (XF) methods to overcome these limitations. We derived mathematical formulations for hybrid quantum-classical XF dynamics that selectively treat critical nuclear degrees of freedom quantum mechanically while propagating others classically. This approach addresses the computational intractability that has previously limited exact methods to toy systems. Key innovations include a new approach to systematically identifying nuclear coordinates requiring quantum treatment, as well as novel implementation strategies that interface with existing quantum chemistry codes. The project also developed a proof-of-concept code for treating Jahn-Teller systems and creation of educational materials on non-adiabatic dynamics geared at the graduate level. The theoretical framework developed will enable future systematically improvable calculations of nuclear quantum effects in realistic molecular systems, filling a critical gap in non-adiabatic dynamics methods. This foundation supports future development of predictive tools for designing energy-relevant photochemical processes where quantum coherence effects may be exploited to control reaction outcomes.
The motion of atoms in programmable optical tweezer arrays offers many new opportunities for neutral atom quantum science. These include inter- and intra-site atom motion for resource-efficient implementations of fermionic and bosonic modes, respectively, as well as tweezer transport for efficient compilation of arbitrary circuits. However, the exploitation of atomic motion for all three purposes and others is limited by the inertia of the atoms. We present a comprehensive architectural blueprint for the use of fermionic metastable helium-3 (He) atoms -- the lightest trappable atomic species -- in programmable optical tweezer arrays. This includes a concrete analysis of atomic structure considerations as well as Rydberg-mediated interactions. We show that inter-tweezer hopping of He atoms can be faster than previous demonstrations with lithium-6. We also demonstrate a new toolbox for encoding and manipulating qubits directly in the tweezer trap potential, uniquely enabled by the light mass of He. Finally, we provide several examples of new opportunities for fermionic quantum simulation and computation that leverage the transport and inter-tweezer hopping of He atom arrays. These tools present new methods to improve the resource efficiency of neutral atom quantum science that may also enable quantum simulations of lattice gauge theories and quantum chemistry outside the Born-Oppenheimer approximation
An ultraviolet spectral probe for a hydrogen-rich planetary atmosphere, such as that of Jupiter, is suggested, utilizing discrete lines in the H2(+) 2p pi u - /s sigma g electronic transition. For the Jovian atmosphere, the dominant mechanism for exciting H2(+) to its 2p pi u state appears to be photoexcitation, principally through absorption of the solar Lyman-alpha line. The critical role of corrections to the Born-Oppenheimer approximation in the use of an H2(+) probe is discussed.
The paper examines the basic principles underlying the molecular transitions responsible for interstellar molecular spectra. The energy levels of molecules are discussed in detail with special attention given to the Born-Oppenheimer approximation, the electronic Hamiltonian, and the parameters of vibrational and rotational energy. The probabilities for radiative and collisional transitions are calculated. A brief review of techniques for molecular spectroscopy is presented along with methods used to determine collision cross sections on both an experimental and a theoretical basis.
Huang and Hughes have discussed the hyperfine splitting delta mu of muonic helium using a variational approach. In this paper, the Born-Oppenheimer approximation is used to simplify the evaluation of delta mu in the nonrelativistic limit. The first-order perturbed wave function of the electron is obtained in closed form by modifying the method used by Dalgarno and Lynn. The result is close to those of Huang and Hughes which required a large Hylleraas expansion and considerable extrapolation.
It is pointed out that with increased computer power and improved computational techniques, such as the gradients developed in recent years, it is becoming practical to compute spectra ab initio, from the fundamental constants of nature, for systems of increasing complexity. The present investigation has the objective to explore several possible ab initio approaches to spectra, giving particular attention to infrared and nonresonance Raman. Two approaches are discussed. The sequential approach, in which first the electronic part and then later the nuclear part of the Born-Oppenheimer approximation is solved, is appropriate for small systems. The simultaneous approach, in which the electronic and nuclear parts are solved at the same time, is more appropriate for many-atom systems. A review of the newer quantum gradient techniques is provided, and the infrared and Raman spectral band contours for the water molecule are computed.
We have calculated static properties of the ground state of the HD(+) ion and its lowest-lying P-state without making use of the Born-Oppenheimer approximation, as was done in the case of H2(+) and D2(+) [Phys. Rev. A 58, 2787 (1998)]. The ion is treated as a three-body system whose ground state is spherically symmetric. The wavefunction is of generalized Hylleraas type, but it is necessary to include high powers of the internuclear distance to localize the nuclear motion. We obtain good values of the energies of the ground S-state and lowest P-state and compare them with earlier calculations. Expectation values are obtained for various operators, the Fermi contact parameters, and the permanent quadrupole moment. The cusp conditions are also calculated. The polarizability was then calculated using second-order perturbation theory with intermediate P pseudostates. Since the nuclei in HD(+) are not of equal mass there is dipole coupling between the lowest two rotational states, which are almost degenerate. This situation is carefully analyzed, and the Stark shift is calculated variationally as a function of the applied electric field.
Wavefunctions of Hylleraas type were used earlier to calculate energy levels of muonic systems. Recently, we found in the case of the molecular ions H2+, D2+ and HD+ that it was necessary to include high powers of the internuclear distance in the Hylleraas functions to localize the nuclear motion when treating the ions as three-body systems without invoking the Born-Oppenheimer approximation. We try the same approach in a muonic system, td(mu-). Improved convergence is obtained for J = 0 and 1 states for shorter expansions when we use this type of generalized Hylleraas function, but as the expansion length increases the high powers are no longer useful. We obtain good energy values for the two lowest J = 0 and J = 1 states and compare them with the best earlier calculations. Expectation values are obtained for various operators, the Fermi contact parameters, and the permanent quadrupole moment. The cusp conditions are also calculated. The polarizability of the ground state is then calculated using second-order perturbation theory with intermediate J = 1 pseudostates. It should be possible to measure the polarizability by observing Rydberg states of atoms with td(mu-) acting as the nucleus.
Here, we perform electronic structure calculations for a set of molecules with degenerate spin-dependent ground states ( 3 CH 2 , 2 CH$^{•}_{3}$, 3 O 2 ) going beyond the Born–Oppenheimer approximation and accounting for nuclear motion. According to a phase space approach that parameterizes electronic states (|Φ⟩) and electronic energies (E) by nuclear position and momentum [i.e., |Φ(R, P)⟩ and E(R, P)], we find that the presence of degenerate spin degrees of freedom leads to broken symmetry ground states. More precisely, rather than a single degenerate minimum at (R, P) = (R min , 0), the ground state energy has two minima at (R,P)=(R' min ,±P min ) (where R' min is close to R min ), dramatically contradicting the notion that the total energy of the system can be written in separable form as E = $\frac{P^2}{2M}$ + V el . Although we find that the broken symmetry solutions have small barriers between them for the small molecules, we hypothesize that the barriers should be macroscopically large for metallic solids, thus offering up a new phase-space potential energy surface for simulating the Einstein–de Haas effect.
The Born-Oppenheimer potentials for QCD with light quarks include adjoint-hadron potentials that are repulsive at short distances and heavy-hadron-pair potentials that approach thresholds at large distances. The adjoint-hadron potentials must connect smoothly to the heavy-hadron-pair potentials at intermediate distances. We identify exotic hidden-heavy hadrons as bound states and resonances in adjoint-hadron potentials that cross below a heavy-hadron-pair threshold before approaching it. This explains why many exotic hidden-charm and hidden-bottom hadrons have energies near heavy-hadron-pair thresholds. The remarkable properties of some exotic hidden-heavy mesons can be explained by fine tunings of adjoint-meson energies in QCD.
In order to model experimental non-resonant Raman optical activity, chemists must compute a host of second-order response tensors (e.g., the electric-dipole–magnetic-dipole polarizability) and their nuclear derivatives along a set of vibrational modes. While these response functions are almost always computed within a Born–Oppenheimer (BO) framework, here we provide a natural interpretation of the electric-dipole–magnetic-dipole polarizability within phase space electronic structure theory, a beyond-BO model whereby the electronic structure depends on nuclear momentum (P) in addition to nuclear position (R). By coupling to nuclear momentum, phase space electronic structure theory is able to capture the asymmetric response of the electronic properties to an external field, in so far as for a vibrating (non-stationary) molecule, $\frac{∂μ}{∂B}$≠$\frac{∂m}{∂F}$, where μ and m are the electrical linear and magnetic dipoles, and F and B are electric and magnetic fields. As an example, for a prototypical methyloxirane molecule, we show that phase space electronic structure theory is able to deliver a reasonably good match with experimental results in a manner that is formally invariant to gauge origin G 0 —provided that one uses a complete basis or, alternatively, gauge invariant atomic orbitals.
Coupled Schroedinger equations for diatomic molecules internal motion solved using best adiabatic approximation to obtain energy correction to Born-Oppenheimer method