SEARCH · Search NASA

Results for “BORN APPROXIMATION”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Long-Range Interaction between Molecules with Electronic Degeneracy: Beyond the Born–Oppenheimer Approximation

Long-range intermolecular interaction plays an essential role in (ultra)cold collisions. For monomers with degenerate electronic states, an accurate description of the long-range interaction can be challenging due to the potential breakdown of the Born–Oppenheimer approximation. In this work, we developed a general method to construct long-range interaction potential energy surfaces (IPESs) between arbitrary molecules using the degenerate perturbation theory of Schrödinger. To further simplify the application of this method, a group-theoretical method is introduced to identify nonzero components of electrostatic properties of monomers, especially for the nonlinear molecules. To demonstrate this approach, we report the long-range IPESs of the O( 3 P)–OH( 2 Π) and OH( 2 Π)–OH( 2 Π) systems. The excellent agreement with ab initio calculations in large separations, even near electronic degeneracies, confirms the accuracy of our approach. The method is expected to lay the groundwork for the accurate investigation of (ultra)cold collisional dynamics involving open-shell species.

ab initio calculations↗

DOE OSTI2024

Deriving Effective Coupling Strength with Born-Oppenheimer Approximation

The development of high-fidelity quantum gates is paramount to the scalability of quantum processors. Tunable couplers have played a key role in the realization of these low error rate quantum gates in superconducting circuits. However, the derivation of the effective coupling strength between the qubits is usually enabled by a Schrieffer-Wolff transformation, which is complex and requires prior-knowledge of the proper generator for the transformation. We propose a simpler method of obtaining this effective coupling using an approach similar to the Born-Oppenheimer approximation.

Wichmann, Conrad↗

DOE OSTI2025

Tetraquarks in the Born-Oppenheimer approximation

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.

Effective Field Theories of QCD↗

DOE OSTI2025

Microscopic calculations with noniterative finite amplitude methods and the application to neutron radiative captures and inelastic scatterings

We derive the fully self-consistent quasiparticle random-phase approximation (QRPA) equations with noniterative finite amplitude methods and calculate the transition strengths of giant resonances. Then, we apply the QRPA results to both neutron radiative capture calculations based on the statistical Hauser-Feshbach theory and inelastic scattering calculations based on distorted-wave Born approximation (DWBA). We compare the calculated results with available experimental data and demonstrate how our approach can reproduce giant resonances and various nuclear reactions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

coh3

CoH3 (CoH ver.3) is an optical model, exciton pre-equilibrium, and Hauser-Feshbach statistical model code, which calculates nuclear reaction cross sections for medium to heavy targets in the keV to MeV energy region. This program is written in standard C++, divided into approximately 200 source and header files. CoH solves the Schroedinger equation for optical potentials defined in the code, and calculates differential elastic scattering, reaction, and total cross sections, for neutron, proton, deuteron, triton, 3He, and alpha-particle. Deformed optical potentials are solved with the coupled-channels method, in which the ground state rotational band members, or vibrational phonon states are coupled. The optical model gives particle transmission coefficients that are fed into the statistical model calculations. CoH includes the pre-equilibrium model (exciton model), the direct/semidirect capture model, and the multi-stage Hauser-Feshbach statistical decay with width fluctuation correction based on the Gaussian orthogonal ensemble. For weakly coupled levels, the DWBA (distorted wave Born approximation) method is used to calculate the direct inelastic scattering process to the excited states.

Kawano, Toshihiko↗

DOE OSTI2025

Classical eikonal from Magnus expansion

In a classical scattering problem, the classical eikonal is defined as the generator of the canonical transformation that maps in-states to out-states. It can be regarded as the classical limit of the log of the quantum S-matrix. In a classical analog of the Born approximation in quantum mechanics, the classical eikonal admits an expansion in oriented tree graphs, where oriented edges denote retarded/advanced worldline propagators. The Magnus expansion, which takes the log of a time-ordered exponential integral, offers an efficient method to compute the coefficients of the tree graphs to all orders. We exploit a Hopf algebra structure behind the Magnus expansion to develop a fast algorithm which can compute the tree coefficients up to the 12th order (over half a million trees) in less than an hour. In a relativistic setting, our methods can be applied to the post-Minkowskian (PM) expansion for gravitational binaries in the worldline formalism. We demonstrate the methods by computing the 3PM eikonal and find agreement with previous results based on amplitude methods. Importantly, the Magnus expansion yields a finite eikonal, while the naïve eikonal based on the time-symmetric propagator is infrared-divergent from 3PM on.

Black Holes↗

DOE OSTI2025

Observation of three-proton-decaying resonant states in 21 Al

The structure of the proton-unbound nucleus 21 Al has been studied via proton transfer from a 9 Be target to a 20 Mg beam at E/A = 103 MeV. Three resonant states in 21 Al are identified by the invariant-mass reconstruction of 18 Ne + 3p events. The decay-energy spectra for the two subsystems of 21 Al , i.e. 18 Ne+ 1 p and 18 Ne+ 2 p , indicate that the observed resonant state with the lowest energy has decay branches consisting of three sequential steps of one-proton decay. Based on the ab initio and the distorted-wave Born approximation (DWBA) calculations, this resonant state at Ε ⁎ = 6.09(11) MeV likely has a spin-parity of (1/2 + - 7/2 + ). Meanwhile, the possibility that the observed resonant peak represents an admixture of several states cannot be ruled out. The good isospin-symmetry of the low-lying states with d 5/2 , d 3/2 configurations in the 21 Al- 21 O pair and the ab initio calculations indicate weak configuration mixing between the s - and d -orbital components. The systematics of experimental energy differences between the 5/2 + ground state and lowest-lying 3/2 + excited state for isotopes (isotones) of Z(N) = 13 suggest a likely semimagic nature of 21 Al.

21Al↗

DOE OSTI2025

Valence 1⁢𝑠−0⁢𝑑 proton vacancy of the 32 Si ground state

The 32 Si( 3 He,d) 33 P reaction was studied in inverse kinematics at 6.3 MeV/u. States in 33 P corresponding to the proton 1s-0d single-particle orbitals were identified up to ≈ 4.5 MeV in excitation energy. The ( 3 He,d) spectroscopic factors were determined from distorted-wave Born approximation calculations. When combined with complementary neutron-adding data, the 1s-0d proton vacancies in the 32 Si ground state were extracted. In conjunction with a reanalysis of data from previous single-particle measurements, the trends in proton and neutron vacancy were explored across the 28,30,32,34 Si isotopes. Both proton and neutron vacancy data show gradual changes in their occupancies. The proton 1s 1/2 orbitals in 32 Si and 34 Si are both consistent with being empty. In conclusion, the ground-state nucleon distributions are described by shell-model calculations constrained to the 1s-0d model space.

Watwood, N. [Argonne National Laboratory (ANL), Ar↗

DOE OSTI2024

Single-neutron adding on $^{34}$S

For this work, single-neutron adding data was collected in order to determine the distribution of the single-neutron strength of the 0ƒ 7/2 , 1p 3/2 , 1p 1/2 and 0ƒ 5/2 orbitals outside of Z = 16, N = 18, 34 S. The 34 S(d,p) 35 S reaction has been measured at 8 MeV/u to investigate cross sections to excited states in 35 S. Outgoing proton yields and momenta were analyzed by the Super-Enge Split-Pole Spectrograph in conjunction with the CeBrA demonstrator located at the John D. Fox Laboratory at Florida State University. Angular distributions were compared with Distorted Wave Born Approximation calculations in order to extract single-neutron spectroscopic overlaps. Spectroscopic overlaps and strengths were determined for states in S up through 6 MeV in excitation energy. Each orbital was observed to have fragmented strength where a single level carried the majority. The single-neutron centroids of the 0ƒ 7/2 , 1p 3/2 , 1p 1/2 and 0ƒ 5/2 orbitals were determined to be 2360$^{+90}_{-40}$ keV, 3280$^{+80}_{-50}$ keV, 4780$^{+60}_{-40}$ keV, and ≳7500 keV, respectively. A previous discrepancy in the literature with respect to the distribution of the neutron 1p 1/2 strength was resolved. The integration of the normalized spectroscopic strengths, up to 5.1 MeV in excitation energy, revealed fully-vacant occupancies for the 0ƒ 7/2 , 1p 3/2 , and 1p 1/2 orbitals, as expected. The spacing in the single-neutron energies highlighted a reduction in the traditional N = 28 shell-gap, relative to both the 1p spin-orbit energy difference (N = 32) and the lower limit on the N = 34 shell spacing.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

Stable isotope equilibria in the dihydrogen-water-methane-ethane-propane system. Part 1: Path-integral calculations with CCSD(T) quality potentials

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.

03 NATURAL GAS↗

DOE OSTI2025

Symmetry breaking as predicted by a phase space Hamiltonian with a spin Coriolis potential

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.

Berry connection↗

DOE OSTI2026

Shadow molecular dynamics for flexible multipole models

Shadow molecular dynamics provide an efficient and stable atomistic simulation framework for flexible charge models with long-range electrostatic interactions. Shadow molecular dynamics simulations are driven by approximate “shadow” Born–Oppenheimer potentials for which the exact charges and forces are directly accessible without relying on costly (and approximate) iterative solvers. While previous implementations have been limited to atomic monopole charge distributions, we extend this approach to flexible multipole models. We derive detailed expressions for the shadow energy functions, potentials, and force terms, explicitly incorporating monopole–monopole, dipole–monopole, and dipole–dipole interactions. In our formulation, both atomic monopoles and atomic dipoles are treated as extended dynamical variables alongside the propagation of the nuclear degrees of freedom. We demonstrate that introducing the additional dipole degrees of freedom preserves the stability and accuracy previously seen in monopole-only shadow molecular dynamics simulations. In addition, we present a shadow molecular dynamics scheme where the monopole charges are held fixed while the dipoles remain flexible. Our extended shadow dynamics provide a framework for stable, computationally efficient, and versatile molecular dynamics simulations involving long-range interactions between flexible multipoles. This is of particular current interest in combination with machine-learned interatomic potentials, including long-range electrostatic interactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Attosecond response of molecules to impulsive ionization

When matter interacts with energetic radiation it can undergo sudden, or impulsive, ionization. This process can drive chemical change and occurs widely in space and planetary atmospheres, yet its comprehensive description challenges our current theoretical and computational capabilities as it requires advanced treatment of electron correlation and nonadiabatic dynamics beyond the Born–Oppenheimer approximation. Here, in this study, we measure the response of the para-aminophenol molecule to sudden ionization. Using attosecond X-ray absorption spectroscopy, we resolve the ultrafast dynamics of the ionized molecule with atomic precision. A subfemtosecond decay corresponds to states undergoing non-radiative decay, whereas few-femtosecond oscillatory signatures are associated with electronic wavepacket motion in stable cation states that later couple to nuclear motion. We compare our measurement with state-of-the-art computational modelling, qualitatively reproducing the observed response across multiple timescales. These results provide a benchmark for computational models of sudden ionization and ultrafast charge motion in matter.

Driver, Taran [SLAC National Accelerator Laborator↗

DOE OSTI2025

SchrödingerNet: A Universal Neural Network Solver for the Schrödinger Equation

Recent advances in machine learning have facilitated numerically accurate solution of the electronic Schrödinger equation (SE) by integrating various neural network (NN)-based wave function ansatzes with variational Monte Carlo methods. Nevertheless, such NN-based methods are all based on the Born–Oppenheimer approximation (BOA) and require computationally expensive training for each nuclear configuration. In this work, we propose a novel NN architecture, SchrödingerNet, to solve the full electronic-nuclear SE by defining a loss function designed to equalize local energies across the system. This approach is based on a translationally, rotationally and permutationally symmetry-adapted total wave function ansatz that includes both nuclear and electronic coordinates. Furthermore, this strategy not only allows for an efficient and accurate generation of a continuous potential energy surface at any geometry within the well-sampled nuclear configuration space, but also incorporates non-BOA corrections, through a single training process. Comparison with benchmarks of atomic and small molecular systems demonstrates its accuracy and efficiency.

Chemical calculations↗

DOE OSTI2025

Direct Observation of the ππ* to nπ* Transition in 2-Thiouracil via Time-Resolved NEXAFS Spectroscopy

The photophysics of nucleobases has been the subject of both theoretical and experimental studies over the past decades due to the challenges posed by resolving the steps of their radiationless relaxation dynamics, which cannot be described in the framework of the Born–Oppenheimer approximation (BOA). In this context, the ultrafast dynamics of 2-thiouracil has been investigated with a time-resolved NEXAFS study at the Free Electron Laser FLASH. Near Edge X-ray Absorption Fine Structure spectroscopy (NEXAFS) can be used to observe electronic transitions in ultrafast molecular relaxation. We performed time-resolved UV-pump/X-ray probe absorption measurements at the sulfur 2s (L1) and 2p (L2/3) edges. We are able to identify absorption features corresponding to the S2 (ππ*) and S1 (nπ*) electronic states. We observe a delay of 102 ± 11 fs in the population of the nπ* state with respect to the initial optical excitation and interpret the delay as the time scale for the S2 → S1 internal conversion. We furthermore identify oscillations in the absorption signal that match a similar observation in our previous X-ray photoelectron spectroscopy study on the same molecule.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Non-Markovian relaxation spectroscopy of fluxonium qubits

Recent studies have shown that parasitic two-level systems (TLS) in superconducting qubits, which are a leading source of decoherence, can have relaxation times longer than the qubits themselves. However, the standard techniques used to characterize qubit relaxation is only valid for measuring T 1 under the Born-Markov approximation and could mask environmental memory effects in practice. Here, we introduce two-timescale relaxometry, a technique to probe the qubit and environment relaxation simultaneously and efficiently. We apply it to high-coherence fluxonium qubits over a frequency range of 0.1-0.4 GHz, and reveal a discrete spectrum of TLS with millisecond lifetimes. Our analysis of the spectrum is consistent with a random distribution of TLS in the aluminum oxide tunnel barrier of the Josephson junction chain of the fluxonium, with a spectral and volumetric density and average electric dipole similar to previous TLS studies at much higher frequencies. Our study suggests that investigating and mitigating TLS in the junction chain is crucial to the development of various types of noise-protected qubits in circuit QED.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Constrained nuclear–electronic orbital method for periodic density functional theory: Application to H 2 chemisorption on Si(001) surfaces

The nuclear–electronic orbital (NEO) method provides a powerful computational framework for incorporating nuclear quantum effects (NQE) in electronic structure calculations beyond the Born–Oppenheimer approximation. By incorporating additional constraints to the position operator on quantum particles like protons, the NEO method enables calculation of effective potential that accounts for NQE. Here, in this work, we present a new constrained NEO (cNEO) formulation for density functional theory (cNEO-DFT) calculations in the context of extended periodic systems. Using the nudged elastic band method, we discuss an application of the cNEO-DFT approach to studying the adsorption of a hydrogen molecule on the Si(001) surfaces. The calculation shows how NQE impacts the reaction energetics. The proton density changes are computed along the reaction pathways. This work demonstrates the capability of the new cNEO-DFT method to study a wide range of chemical processes, such as surface reactions where the quantum nature of light atoms like protons is non-negligible.

Chemical processes↗

DOE OSTI2025

Nondegenerate Noise-Resilient Superconducting Qubit

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.

quantum error correction↗