Engineering Papers⌕ Search

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.

At least 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

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↗
DOE OSTI2018

Uncertainty quantification due to optical potentials in models for ( d , p ) reactions

Recent work has studied the uncertainty in predictions for $A(d,p)B$ reactions using the distorted-wave Born approximation (DWBA), coming from the parametrization of the effective $dA$ interactions [Lovell et al. Phys. Rev. C 95, 024611 (2017)]. There are different levels of sophistication in reaction theories for one-nucleon transfer reactions, including the adiabatic wave approximation (ADWA) which takes deuteron breakup into account to all orders. In this paper, we quantify the uncertainties associated with the ADWA method that come from the parametrization of the NA interactions and compare ADWA with DWBA. Similar to Lovell et al., we use nucleon elastic-scattering data on a wide variety of targets at the appropriate incoming and outgoing energies to constrain the optical potential input to the ADWA theory. Pulling from the χ 2 distribution, we obtain 95% confidence bands for the elastic distributions. From the resulting parameters, we predict 95% confidence bands for the ( d , p ) transfer cross sections. Results obtained with the standard uncorrelated χ 2 are compared to those using the correlated χ 2 of Lovell et al. We also repeat the DWBA calculations for the exact same reactions for comparison purposes. We find that NA elastic-scattering data provide a significant constraint to the interactions, and, when the uncertainties are propagated to the transfer reactions using ADWA, predictions are consistent with the transfer data. The angular distributions for ADWA differ from those predicted by DWBA, particularly at small angles. As in Lovell et al., confidence bands obtained using the uncorrelated χ 2 function are unrealistically narrow and become much wider when the correlated χ 2 function is considered. For most cases, the uncertainty bands obtained in ADWA are narrower than DWBA when using elastic data of similar quality and range. However, given the large uncertainties predicted from the correlated χ 2 function, at this point, the transfer data cannot discriminate between these two methods.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Validity of the Born-Oppenheimer approximation in the indirect-dissociative-recombination process

An alternative method is introduced to solve a simple two-dimensional model describing vibrational excitation and dissociation processes during the electron-molecule collisions. The model works with one electronic and one nuclear degree of freedom. The two-dimensional R matrix can be constructed simultaneously on the electronic and nuclear surfaces using all three forms developed previously for electron-atom and electron-molecule collisions. These are the eigenchannel R-matrix form, inversion technique of Nesbet and Robicheaux, and the Wigner-Eisenbud-type form using expansion over the poles of the symmetrized Hamiltonian. The 2D R-matrix method is employed to solve a simple model tailored to describe the dissociative recombination and the vibrational excitation of H 2 + cation in the singlet ungerade symmetry 1 Σ u . These results then serve as a (near-exact) benchmark for the following calculation in which the R-matrix states are replaced by their Born-Oppenheimer approximations. Furthermore, the accuracy of this approach and its correction with the first-order nonadiabatic couplings are discussed.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

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↗
DOE OSTI2026

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↗
DOE OSTI2023

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↗
NASA NTRS1984

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.↗
DOE OSTI2014

Braiding fluxes in Pauli Hamiltonian

Aharonov and Casher showed that Pauli Hamiltonians in two dimensions have gapless zero modes. We study the adiabatic evolution of these modes under the slow motion of N fluxons with fluxes Φ{sub a}∈R. The positions, r{sub a}∈R{sup 2}, of the fluxons are viewed as controls. We are interested in the holonomies associated with closed paths in the space of controls. The holonomies can sometimes be abelian, but in general are not. They can sometimes be topological, but in general are not. We analyse some of the special cases and some of the general ones. Our most interesting results concern the cases where holonomy turns out to be topological which is the case when all the fluxons are subcritical, Φ{sub a}<1, and the number of zero modes is D=N−1. If N≥3 it is also non-abelian. In the special case that the fluxons carry identical fluxes the resulting anyons satisfy the Burau representations of the braid group.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

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↗
DOE OSTI2024

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↗
DOE OSTI2014

A full-dimensional multilayer multiconfiguration time-dependent Hartree study on the ultraviolet absorption spectrum of formaldehyde oxide

Employing the multilayer multiconfiguration time-dependent Hartree (ML-MCTDH) method in conjunction with the multistate multimode vibronic coupling Hamiltonian (MMVCH) model, we perform a full dimensional (9D) quantum dynamical study on the simplest Criegee intermediate, formaldehyde oxide, in five lower-lying singlet electronic states. The ultraviolet (UV) spectrum is then simulated by a Fourier transform of the auto-correlation function. The MMVCH model is built based on extensive MRCI(8e,8o)/aug-cc-pVTZ calculations. To ensure a fast convergence of the final calculations, a large number of ML-MCTDH test calculations is performed to find an appropriate multilayer separations (ML-trees) of the ML-MCTDH nuclear wave functions, and the dynamical calculations are carefully checked to ensure that the calculations are well converged. To compare the computational efficiency, standard MCTDH simulations using the same Hamiltonian are also performed. A comparison of the MCTDH and ML-MCTDH calculations shows that even for the present not-too-large system (9D here) the ML-MCTDH calculations can save a considerable amount of computational resources while producing identical spectra as the MCTDH calculations. Furthermore, the present theoretical B{sup ~} {sup 1}A{sup ′}←X{sup ~} {sup 1}A{sup ′} UV spectral band and the corresponding experimental measurements [J. M. Beames, F. Liu, L. Lu, and M. I. Lester, J. Am. Chem. Soc. 134, 20045–20048 (2012); L. Sheps, J. Phys. Chem. Lett. 4, 4201–4205 (2013); W.-L. Ting, Y.-H. Chen, W. Chao, M. C. Smith, and J. J.-M. Lin, Phys. Chem. Chem. Phys. 16, 10438–10443 (2014)] are discussed. To the best of our knowledge, this is the first theoretical UV spectrum simulated for this molecule including nuclear motion beyond an adiabatic harmonic approximation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Benchmark calculations of electron impact electronic excitation of the hydrogen molecule

We present benchmark integrated and differential cross-sections for electron collisions with H 2 using two different theoretical approaches, namely, the R-matrix and molecular convergent close-coupling. This is similar to comparative studies conducted on electron–atom collisions for H, He and Mg. Electron impact excitation to the ${\mathrm{b}}^{3}{{\Sigma}}_{\mathrm{u}}^{+}$, ${\mathrm{a}}^{3}{{\Sigma}}_{\mathrm{g}}^{+}$, ${\mathrm{B}}^{1}{{\Sigma}}_{\mathrm{u}}^{+}$, c3Πu, ${\mathrm{E}\mathrm{F}}^{1}{{\Sigma}}_{\mathrm{g}}^{+}$, C1Πu, ${\mathrm{e}}^{3}{{\Sigma}}_{\mathrm{u}}^{+}$, ${\mathrm{h}}^{3}{{\Sigma}}_{\mathrm{g}}^{+}$, ${\mathrm{B}}^{\prime 1}{{\Sigma}}_{\mathrm{u}}^{+}$ and d3Πu excited electronic states are considered. Calculations are presented in both the fixed nuclei and adiabatic nuclei approximations, where the latter is shown only for the ${\mathrm{b}}^{3}{{\Sigma}}_{\mathrm{u}}^{+}$ state. Good agreement is found for all transitions presented. Where available, we compare with existing experimental and recommended data.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2015

CMB hemispherical asymmetry from non-linear isocurvature perturbations

We investigate whether non-adiabatic perturbations from inflation could produce an asymmetric distribution of temperature anisotropies on large angular scales in the cosmic microwave background (CMB). We use a generalised non-linear δ N formalism to calculate the non-Gaussianity of the primordial density and isocurvature perturbations due to the presence of non-adiabatic, but approximately scale-invariant field fluctuations during multi-field inflation. This local-type non-Gaussianity leads to a correlation between very long wavelength inhomogeneities, larger than our observable horizon, and smaller scale fluctuations in the radiation and matter density. Matter isocurvature perturbations contribute primarily to low CMB multipoles and hence can lead to a hemispherical asymmetry on large angular scales, with negligible asymmetry on smaller scales. In curvaton models, where the matter isocurvature perturbation is partly correlated with the primordial density perturbation, we are unable to obtain a significant asymmetry on large angular scales while respecting current observational constraints on the observed quadrupole. However in the axion model, where the matter isocurvature and primordial density perturbations are uncorrelated, we find it may be possible to obtain a significant asymmetry due to isocurvature modes on large angular scales. Such an isocurvature origin for the hemispherical asymmetry would naturally give rise to a distinctive asymmetry in the CMB polarisation on large scales.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Adiabatic regularisation of power spectra in k-inflation

We look at the question posed by Parker et al. about the effect of UV regularisation on the power spectrum for inflation. Focusing on the slow-roll k-inflation, we show that up to second order in the Hubble and sound flow parameters, the adiabatic regularisation of such model leads to no difference in the power spectrum apart from certain cases that violate near scale-invariant power spectra. Furthermore, extending to non-minimal k-inflation, we establish the equivalence of the subtraction terms in the adiabatic regularisation of the power spectrum in Jordan and Einstein frames.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Angle-dependent strong-field ionization and fragmentation of carbon dioxide measured using rotational wave packets

In this work, we experimentally study the angle-dependent single ionization of carbon dioxide (CO 2 ) by linearly and circularly polarized pulses. The angle dependence of the ionization probability by linearly polarized pulses extracted from time-domain measurements on an impulsively excited rotational wave packet is compared with data obtained from a direct angle-scan measurement. The results from the measurement with linear and circular polarization are consistent with the adiabatic ionization approximation. Further, we extend the time-domain method to extract the dependence of the asymptotic momentum distribution of fragment ions on the orientation of the molecular axis, and apply it to investigate dissociative double ionization of CO 2 . We show that such measurements can directly test the validity of the axial recoil approximation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Temporal coherent control of resonant two-photon double ionization of the hydrogen molecule via doubly excited states

Here, we use time-delayed, counter-rotating, circularly polarized few-cycle attosecond nonoverlapping pulses to study the temporal coherent control of the resonant process of two-photon double ionization (TPDI) of hydrogen molecule via doubly excited states for pulse propagation direction along $\hat k$ either parallel or perpendicular to the molecular axis $\hat R$. For $\hat k$∥ $\hat R$ and a pulse carrier frequency of 36 eV resonantly populating the Q 2 1 Π$^+_u$ (1) doubly excited state as well as other 1 Π$^+_u$ doubly excited states, we find that the indirect ionization pathway through these doubly excited states changes the character of the kinematical vortex-shaped momentum distribution produced by the two direct ionization pathways from fourfold to twofold rotational symmetry. This result is similar to what found in TPDI of the He atom involving 1 P$^ o_ {±1}$ doubly excited states; however, angular distributions exhibiting a quantum beat effect between the ground state and a doubly excited state seen for the He atom are observed here for its molecular counterpart with an anomaly in shape and magnitude, not in frequency. The sixfold differential probability integrated over the azimuthal angle of the photoelectron pair shows that this anomaly is due to autoionization decays and quantum beats between doubly excited states. For $\hat k$⟂$\hat R$ and a broadband pulse carrier frequency of 30 eV populating the Q 1 1 Π$^+_u$(1), Q 1 1 Σ$^+_u$(1), Q 2 1 Π$^+_u$(1) , and Q 1 1 Σ$^+_u$ (2) doubly excited states, the momentum distribution is shown to exhibit dynamical electron vortices with four spiral arms, which originates from the interplay between the 1 Δ$^+_ g$ , 1 Π$^+_ g$, and 1 Σ$^+_ g$ dynamical ionization amplitudes. Our treatment within either the adiabatic-nuclei approximation or fixed-nuclei approximation shows that the latter provides a very good account for this correlated process.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

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↗
DOE OSTI2025

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↗