Corrections to the Born-oppenheimer Approximation, Part II
Adiabatic corrections to long-range Born- Oppenheimer potentials
ADIABAT↗
SEARCH · Engineering Papers
Results for “Adiabatic 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.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Development of fixed-nuclei approximations of electron-molecule scattering and the adiabatic theory of rotational excitation
Fixed nuclei approximations of electron-molecule scattering and adiabatic theory of rotational excitation
Temkin, A.↗
Modeling the Electronic Absorption Spectra of the Indocarbocyanine Cy3
Accurate modeling of optical spectra requires careful treatment of the molecular structures and vibronic, environmental, and thermal contributions. The accuracy of the computational methods used to simulate absorption spectra is limited by their ability to account for all the factors that affect the spectral shapes and energetics. The ensemble-based approaches are widely used to model the absorption spectra of molecules in the condensed-phase, and their performance is system dependent. The Franck–Condon approach is suitable for simulating high resolution spectra of rigid systems, and its accuracy is limited mainly by the harmonic approximation. In this work, the absorption spectrum of the widely used cyanine Cy3 is simulated using the ensemble approach via classical and quantum sampling, as well as, the Franck–Condon approach. The factors limiting the ensemble approaches, including the sampling and force field effects, are tested, while the vertical and adiabatic harmonic approximations of the Franck–Condon approach are also systematically examined. Our results show that all the vertical methods, including the ensemble approach, are not suitable to model the absorption spectrum of Cy3, and recommend the adiabatic methods as suitable approaches for the modeling of spectra with strong vibronic contributions. We find that the thermal effects, the low frequency modes, and the simultaneous vibrational excitations have prominent contributions to the Cy3 spectrum. The inclusion of the solvent stabilizes the energetics significantly, while its negligible effect on the spectral shapes aligns well with the experimental observations.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Biased degenerate ground-state sampling of small Ising models with converged quantum approximate optimization algorithm
The quantum alternating operator ansatz, a generalization of the quantum approximate optimization algorithm (QAOA), is a quantum algorithm used for approximately solving combinatorial optimization problems. QAOA typically uses the transverse field mixer as the driving Hamiltonian. One of the interesting properties of the transverse field driving Hamiltonian is that it results in nonuniform sampling of degenerate ground states of optimization problems. In this study, we numerically examine the fair sampling properties of the transverse field mixer QAOA, and Grover mixer QAOA (GM-QAOA), which provides theoretical guarantees of fair sampling of degenerate optimal solutions, up to a large enough p such that the mean expectation value converges to an optimal approximation ratio of 1. This comparison is performed with high-quality heuristically computed, but not necessarily optimal, QAOA angles, which give strictly monotonically improving solution quality as p increases. These angles are computed using the Julia based numerical simulation software JuliQAOA. Fair sampling of degenerate ground states is quantified using the Shannon entropy of the ground-state amplitudes distribution. The fair sampling properties are reported on several quantum signature Hamiltonians from previous quantum annealing fair sampling studies. Small random fully connected spin glasses are shown, which exhibit exponential suppression of some degenerate ground states with transverse field mixer QAOA. The transverse field mixer QAOA simulations show that some problem instances clearly saturate the Shannon entropy of 0 with a maximally biased distribution that occurs when the learning converges to an approximation ratio of 1 while other problem instances never deviate from a maximum Shannon entropy (uniform distribution) at any p step. Published by the American Physical Society 2025
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Vibrationally resolved cross sections for proton collisions with hydrogen molecules
The semi-classical molecular convergent close-coupling approach has been applied to proton collisions with molecular hydrogen in the ground electronic and vibrational state. Cross sections for electron loss and electron capture agree well with the experimental data in the intermediate energy range where coupling between electronic reaction channels is strong and previously available calculations produce inconsistent results. Results for elastic scattering and total excitation are compared with effective one-electron coupled-channel calculations, showing that accurate target structure makes an important difference to these cross sections at intermediate energies. The excitation cross sections are calculated using both the fixed-nuclei and adiabatic-nuclei approximations for transitions from the $X^1Σ^+_{\textrm{g}}$ state to the $B^1Σ^+_{\textrm{u}}$, $B'$$^1Σ^+_{\textrm{u}}$, $B''$$^1Σ^+_{\textrm{u}}$, $EF$$^1Σ^+_{\textrm{g}}$, $GK$$^1Σ^+_{\textrm{g}}$, $H$$^1Σ^+_{\textrm{g}}$, $C$$^1Π_{\textrm{u}}$, $D$$^1Π_{\textrm{u}}$, $D'$$^1Π_{\textrm{u}}$, $I$$^1Π_{\textrm{u}}$, and $J$$^1Δ_{\textrm{g}}$ states for proton energies from to keV. We find significant differences between the present ab initio calculations and the equivelocity scaled electron-scattering data that is currently used for collisional-radiative modelling. Furthermore, the adiabatic-nuclei calculations enable us to resolve the final vibrational level after excitation, producing a complete set of cross sections for excitation of ground-state molecular hydrogen for all electronic states up to $n = 3$, where $n$ is the united-atoms-limit principle quantum number.
74 ATOMIC AND MOLECULAR PHYSICS↗
Stochastic and mixed density functional theory within the projector augmented wave formalism for simulation of warm dense matter
Stochastic density functional theory (DFT) and mixed stochastic-deterministic DFT are burgeoning approaches for the calculation of the equation of state and transport properties in materials under extreme conditions. In the intermediate warm dense matter regime, a state between correlated condensed matter and kinetic plasma, electrons can range from being highly localized around nuclei to delocalized over the whole simulation cell. The plane-wave basis pseudopotential approach is thus the typical tool of choice for modeling such systems at the DFT level. Unfortunately, stochastic DFT methods scale as the square of the maximum plane-wave energy in this basis. To reduce the effect of this scaling and improve the overall description of the electrons within the pseudopotential approximation, we present stochastic and mixed DFT approaches developed and implemented within the projector augmented wave formalism. In conclusion, we compare results between the different DFT approaches for both single-point and molecular dynamics trajectories and present calculations of self-diffusion coefficients of solid density carbon from 1 to 50 eV.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Scattering of electrons from neon atoms
Scattering of electrons from neon atoms is investigated by the polarized-orbital method. The perturbed orbitals calculated with use of the Sternheimer approximation lead to the polarizability 2.803 a(0)-cube in fairly good agreement with the experimental value 2.66 a(0)-cube. Phase shifts for various partial waves are calculated in the exchange, exchange-adiabatic, and polarized-orbital approximations. They are compared with the previous results. The calculated elastic differential, total, and momentum-transfer cross sections are compared with the experimental results. The polarized-orbital approximation yields results which show general improvement over the exchange-adiabatic approximation.
Dasgupta, A.↗
Signature of 0 + excited state and shape coexistence in 94 Kr through 93 Kr(d,p) 94 Kr reaction
A measurement of the excitation spectrum of 94 Kr via one-neutron transfer to the ground state of 93 Kr using the 93 Kr(d,p) 94 Kr reaction at 8A MeV, observing the outgoing protons, performed with the IRIS facility at TRIUMF is reported. Two states in 94 Kr, at 1.50±0.14 MeV and 2.20±0.14 MeV, were observed. An adiabatic-wave approximation analysis of the differential cross sections leads us to identify the lower energy state as being populated with neutron transfer to the 3 s1/2 orbital. This leads to the first observation of the lowest 0 + excited state in 94 Kr, hence signaling shape co-existence. Theoretical calculations performed within the in-medium similarity renormalization group framework are presented that also show the existence of a low-energy 0 + state, aligning qualitatively with the observation.
74 ATOMIC AND MOLECULAR PHYSICS↗
Importance of Electron Correlation on the Geometry and Electronic Structure of [2Fe–2S] Systems: A Benchmark Study of the [Fe 2 S 2 (SCH 3 ) 4 ] 2–,3–,4– , [Fe 2 S 2 (SCys) 4 ] 2– , [Fe 2 S 2 (S- p -tol) 4 ] 2– , and [Fe 2 S 2 (S- o -xyl) 4 ] 2– Complexes
Iron–sulfur clusters are crucial for biological electron transport and catalysis. Obtaining accurate geometries, energetics, manifolds of their excited electronic states, and reduction energies is important to understand their role in these processes. Here, using a [2Fe–2S] model complex with Fe II and Fe III oxidation states, which leads to different charges, i.e., [Fe 2 S 2 (SMe) 4 ] 2–,3–,4– , we benchmarked a variety of computational methodologies ranging from density functional theory (DFT) to post-Hartree–Fock methods, including complete active space self-consistent field (CASSCF), multireference configuration interaction, the second-order N-electron valence state perturbation theory (NEVPT2), and the linearized integrand approximation of adiabatic connection (AC0) approaches. Additionally, we studied three experimentally well-characterized complexes, [Fe 2 S 2 (SCys) 4 ] 2– , [Fe 2 S 2 (S-o-tol) 4 ] 2– , and [Fe 2 S 2 (S-o-xyl) 4 ] 2– , via DFT methods. We conclude that the dynamic electron correlation is important for accurately predicting the geometry of these complexes. Broken symmetry (BS) DFT correctly predicts experimental geometries of low-spin multiplicity, while CASSCF does not. However, BS-DFT significantly overestimates the difference between the low- and high-spin electronic states for a given oxidation state. At the same time, CASSCF underestimates it but provides relative energies closer to the reference NEVPT2 results. Finally, AC0 provides energetics of NEVPT2 quality with the additional advantage of being able to use large CASSCF sizes. NEVPT2 gives the best estimates of the Fe III /Fe III → Fe II /Fe III (4.27 eV) and Fe II /F III → Fe II /F II (7.72 eV) reduction energies. The results provide insight into the electronic structure of these complexes and assist in the understanding of their physical properties.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
All-electron molecular tunnel ionization based on the weak-field asymptotic theory in the integral representation
Tunnel ionization (TI) underlies many important ultrafast processes, such as high-harmonic generation andstrong-field ionization. Among the existing theories for TI, many-electron weak-field asymptotic theory (ME-WFAT) is by design capable of accurately treating many-electron effects in TI. An earlier version of ME-WFATrelied on an accurate representation of the asymptotic tail of the orbitals, which hindered its implementation inGaussian-basis-set-based quantum chemistry programs. In this work, we reformulate ME-WFAT in the integralrepresentation, which makes the quality of the asymptotic tail much less critical, hence greatly facilitating itsimplementation in standard quantum chemistry packages. The integral reformulation introduced here is thereforemuch more robust when applied to molecules with arbitrary geometry. Here, we present several case studies, amongwhich is the CO molecule where some earlier theories disagree with experiments. Here we find that ME-WFATproduces the largest ionization probability when the field points from C to O, as experiments suggest. Anattractive feature of ME-WFAT is that it can be used with various types of multielectron methods whether ofdensity functional or multiconfiguration types, this inturn facilitates tunnel ionization calculation in systems exhibiting a strong multireference character.
74 ATOMIC AND MOLECULAR PHYSICS↗
Asymptotic errors in adiabatic evolution
The adiabatic theorem in quantum mechanics implies that if a system is in a discrete eigenstate of a Hamiltonian and the Hamiltonian evolves in time arbitrarily slowly, the system will remain in the corresponding eigenstate of the evolved Hamiltonian. Understanding corrections to the adiabatic result that arise when the evolution of the Hamiltonian is slow—but not arbitrarily slow—has become increasingly important, especially since adiabatic evolution has been proposed as a method of state preparation in quantum computing. Here, this paper identifies two regimes, an adiabatic regime in which corrections are generically small and can depend on details of the evolution throughout the path, and a hyperadiabatic regime in which the error is given by a form similar to an asymptotic expansion in the inverse of the evolution time with the coefficients depending principally on the behavior at the endpoints. However, the error in this hyperadiabatic regime is neither given by a true asymptotic series nor solely dependent on the endpoints: the coefficients combine the contributions from both endpoints, with relative phase factors that depend on the average spectral gaps along the trajectory, multiplied by the evolution time. The central result of this paper is to identify a quantity, referred to as the typical error, which is obtained by appropriately averaging the error over evolution times that are small compared to the evolution time itself. This typical error is characterized by an asymptotic series and depends solely on the endpoints of the evolution, remaining independent of the details of the intermediate evolution.
adiabatic approximation↗
Rotational-vibrational coupling in the theory of electron-molecule scattering
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.
Temkin, A.↗
Elastic scattering and rotational excitation of a polyatomic molecule by electron impact - Acetylene
Differential, integral, momentum transfer, and partial cross sections have been calculated for elastic scattering and rotational excitation of C2H2 by 10-eV electrons. The effective potential includes static, exchange, and polarization interactions calculated by the INDOX/1s method and the semiclassical exchange approximation with adiabatic polarization at large electron-molecule distances. The scattering is treated by well converged rotational close coupling using the centrifugal dominant scheme to select the channels included and including up to 32 coupled channels for a given total angular momentum. The calculated integral cross sections for pure elastic scattering and rotation excitation are 54.5 and 41.4 a(0)squared, respectively. These are much larger than the values (34.4 and 18.6 a(0)squared) previously (Onda and Truhlar, 1979) calculated for the isoelectronic molecule N2, at this energy. This illustrates how the greater spatial extent of C2H2 greatly increases the cross sections for pure elastic and rotationally inelastic scattering.
Thirumalai, D.↗
Onion-shell model of cosmic ray acceleration in supernova remnants
A method is devised to approximate the spatially averaged momentum distribution function for the accelerated particles at the end of the active lifetime of a supernova remnant. The analysis is confined to the test particle approximation and adiabatic losses are oversimplified, but unsteady shock motion, evolving shock strength, and non-uniform gas flow effects on the accelerated particle spectrum are included. Monoenergetic protons are injected at the shock front. It is found that the dominant effect on the resultant accelerated particle spectrum is a changing spectral index with shock strength. High energy particles are produced in early phases, and the resultant distribution function is a slowly varying power law over several orders of magnitude, independent of the specific details of the supernova remnant.
Bogdan, T. J.↗
The Compressible Flow Past Various Plane Profiles Near Sonic Velocity
In an earlier report UM No.1117 by Gothert,the single-source method was applied to the compressible flow around circles, ellipses, lunes, and around an elongated body of revolution at different Mach numbers and the results compared as far as possible with the calculations by Lamla ad Busemann. Essentially, it was found that with favorable source arrangement the single-source method is in good agreement with the calculations of the same degree of approximation by.Lamla and Busemann. Near sonic velocity the number of steps must be increased considerably in order to sufficiently approximate the adiabatic curve. After exceeding a certain Mach number where local supersonic fields occur already, it was no longer possible, in spite of the substantially increased number of steps, to obtain a systematic solution because the calculation diverged. This result,was interpreted to mean that above this point of divergence the symmetrical type of flow ceases to exist and changes into the unsymmetrical type characterized by compressibility shocks.
Goethert, B.↗
Approximate bounds and temperature dependence of adiabatic connection integrands for the uniform electron gas
Thermal density functional theory is commonly used in simulations of warm dense matter, a highly energetic phase characterized by substantial thermal effects and by correlated electrons demanding quantum mechanical treatment. Methods that account for temperature dependence, such as Mermin–Kohn–Sham finite-temperature density functional theory and free energy density functional theory, are now employed with more regularity and available in many standard code packages. However, approximations from zero-temperature density functional theory are still often used in temperature-dependent simulations using thermally weighted electronic densities as an input to exchange–correlation functional approximations, a practice known to miss temperature-dependent effects in the exchange–correlation free energy of these systems. Here, in this work, the temperature-dependent adiabatic connection is demonstrated and analyzed using a well-known parameterization of the uniform electron gas free energy. Useful tools based on this formalism for analyzing and constraining approximations of the exchange–correlation at zero temperature are leveraged for the finite-temperature case. Inspired by the Lieb–Oxford inequality, which provides a lower bound for the ground-state exchange–correlation energy, bounds for the exchange–correlation at finite temperatures are approximated for various degrees of electronic correlation.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Approximations to the two-state impact- parameter treatment of heavy-particle collisions.
Adiabatic and distortion approximations to two state treatment of ground state H atoms collisions at low energies
Flannery, M. R.↗