Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ab initio calculations”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

The Heat of Formation of HNO

The HNO molecule is of interest in both combustion and atmospheric chemistry. For example, Guadagnini et al. have recently presented ab initio potential energy surfaces for the three lowest lying electronic states of HNO and then used these in examining several chemical reactions that take place in the combustion of nitrogen containing fuels and in the oxidation of atmospheric nitrogen. We have previously studied the ground state potential energy surface (i.e., stationary points along the HNO reversible reaction HON path), vibrational spectrum (using an accurate quartic force field), zero-point energy, and bonding of HNO using coupled-cluster ab initio methods. HNO is also very interesting because of the unique nature of its bonding characteristics. That is, the potential energy surface is very flat along the H-N bonding coordinate thereby giving unusual harmonic and fundamental vibrational frequencies, and the H-N bond energy is rather weak in comparison to other H-N bond energies. In fact, using experimental heats of formation for HNO, H, and NO, the H- bond energy is computed to be only 49.9 kcal/ mol (298 K). However, ab initio calculations of isodesmic reaction energies involving HNO, FNO, ClNO, and several other molecules have shown that there is an inconsistency in the experimental heats of formation of the XNO (X double bond H, F, and Cl) species. Hence the motivation for this study was to determine a very accurate (Delta)H(sub f, sup o) value for HNO using state-of-the-art ab initio methods. Based on many recent studies it is evident that the singles and doubles coupled-cluster method that includes a perturbational estimate of the effects of connected triple excitations, denoted CCSD(T), in conjunction with large one-particle basis sets should be reliable to better than +/- 0.8 kcal/mol for this quantity. The computational methodology is described in the next section followed by our results and discussion. Conclusions are presented in the final section.

Lee, Timothy J.↗
NASA NTRS1995

Electronic transition moment for the B(2)II-X(2)II system of NO

The upsilon' = 0-3 and 7 vibrational levels of the NO B(2)II state have been selectively excited by laser radiation. The fluorescence spectra together with calculated Franck-Condon factors and r-centroids have been used to evaluate the electronic transition moment. The results for upsilon' = 0-3 are in very good agreement with recent chemiluminescence measurements and ab initio calculations. Furthermore, the data from upsilon' = 7 have been used to extend the empirically determined moment to limits 1.23 and 1.78 A, improving agreement with experimentally determined lifetimes.

Luque, Jorge↗
DOE OSTI2025

Precision Mass Measurements Reveal Low Neutron Pairing in Tin beyond 𝑁=82 and Its Impact on Stellar Nucleosynthesis

We present a study on neutron-rich tin (𝑍 =50) isotopes beyond the doubly closed shell of 𝑁 = 82 through high-precision mass measurements, including the first-ever measurements of the masses of 136 Sn, 137 Sn, and 138 Sn isotopes. These measurements enhance our understanding of the nuclear structure and astrophysical nucleosynthesis in this previously unexplored region. The new mass data are used for evaluation of the final abundances of mass numbers 𝐴 =135 and 137 in 𝑟-process network calculations. Our findings reveal a notable change in the empirical pairing gap for tin isotopes beyond the 𝑁 = 82 closed shell and a shift in the two-neutron-separation energy slope compared to heavier elements above the shell closure. A new set of ab initio calculations effectively describes these observed trends.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Electronic Structure and Bonding of US, SUO, and US 2

Anion photoelectron spectra of US – and US 2 – were recorded using the third (355 nm) and fourth (266 nm) harmonics of an Nd:YAG laser, which yielded vertical detachment energies (VDEs) of 1.71 and 2.02 eV, respectively. The experimental results are supported by extensive relativistic ab initio calculations, primarily at the coupled cluster level of theory, with systematic sequences of correlation consistent basis sets. Calculations include the closely related SUO and SUO – molecules, as well as the oxide congeners UO/UO – and UO 2 /UO 2 – which are well-known experimentally and provide benchmark systems for the sulfide calculations. Adiabatic electron detachment energies (ADEs) are computed for UO – , UO 2 – , US – , SUO – , and US 2 – using the Feller–Peterson–Dixon (FPD) composite approach. Additionally, ADEs are determined for UO – and US – using a spinor-based coupled cluster approach where spin–orbit coupling is included at the orbital level. VDEs are derived from the ab initio results from Franck–Condon simulations of the photoelectron spectra.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1990

Solar flare gamma-ray line shapes

A computer code has been developed which is used to calculate ab initio the laboratory shapes and energy shifts of gamma-ray lines from (C-12)(p, gamma/4.438/)p-prime(C-12) and (O-16)(p, gamma/6.129/)p-prime(O-16) reactions and to calculate the expected shapes of these lines from solar flares. The sensitivity of observable solar flare gamma-ray line shapes to the directionality of the incident particles is investigated for several projectile angular distributions. Shapes of the carbon and oxygen lines are calculated assuming realistic proton energy spectra for particles in circular orbits at the mirror points of magnetic loops, for particle beams directed downward into the photosphere, and for isotropic particle distributions. Line shapes for flare sites near the center of the sun and on the limb are shown for both thin-target and thick-target interaction models.

Werntz, C.↗
NASA NTRS1982

Rotational and vibrational spectra of ethynol from quantum-mechanical calculations

It is noted that ethynol (HCCOH), despite the theoretical prediction of its stability to tautomerization to ketene, has thus far not been observed. It is shown here that the identification of this unknown molecule, both in space and in the laboratory, can be aided by an ab initio calculation of spectroscopic parameters. At the HF/3-21G level, harmonic vibrational frequencies are computed by way of analytic second differentiation of the Hartee-Fock (HF) energy with respect to the nuclear coordinates. After applying an empirical scale factor, the resultant frequencies are (per cm) 473, 517, 773, 841, 1003, 1217, 2206, 3285, and 3418. The computed dipole moment at the CISD/DZ+P level is 1.79 D. At the CISD+Q/DZ+P level, the molecule's rotational constants are determined. After scaling by empirical correction factors, they are used in deriving the 4(04) - 3(03) frequency of 76.81 + or - 0.3 GHz with a triplet splitting of 0.30 + or - 0.01 GHz. The triplet splitting involves 4(14) - 3(13) and 4(13) - 3(12) relative to the 4(04) - 3(03) transition as the central line.

Defrees, D. J.↗
NASA NTRS1991

Calculations concerning the reaction C + H3(+) - CH(+) + H2

The ion-molecule reaction of C + H3(+) - CH(+) + H2 has not been studied in the laboratory but is thought to be important in the gas phase synthesis of organic molecules in dense interstellar clouds. This reaction was studied theoretically by performing quantum chemical ab initio calculations on the potential surface. It is found that there is no activation barrier to the reaction and that it proceeds smoothly to the first excited electronic state of CH(+). The rate coefficient as a function of temperature can then be estimated using the proper long-range potentials. The rate coefficient at 10 K is calculated to be 2.9 x 10 to the -9th cu cm/sec.

Talbi, D.↗
NASA NTRS1990

Non-LTE model calculations for SN 1987A and the extragalactic distance scale

This paper presents model atmospheres for the first week of SN 1987A, based on the luminosity and density/velocity structure from hydrodynamic models of Woosley (1988). The models account for line blanketing, expansion, sphericity, and departures from LTE in hydrogen and helium and differ from previously published efforts because they represent ab initio calculations, i.e., they contain essentially no free parameters. The formation of the UV spectrum is dominated by the effects of line blanketing. In the absorption troughs, the Balmer line profiles were fit well by these models, but the observed emissions are significantly stronger than predicted, perhaps due to clumping. The generally good agreement between the present synthetic spectra and observations provides independent support for the overall accuracy of the hydrodynamic models of Woosley. The question of the accuracy of the Baade-Wesselink method is addressed in a detailed discussion of its approximations. While the application of the standard method produces a distance within an uncertainty of 20 percent in the case of SN 1987A, systematic errors up to a factor of 2 are possible, particularly if the precursor was a red supergiant.

Schmutz, W.↗
NASA NTRS1989

Theoretical study of the structures and electron affinities of the dimers and trimers of the group IB metals (Cu, Ag, and Au)

The molecular structure of both the neutral and negatively charged diatomic and triatomic systems containing the Cu, Ag, and Au metals are determined from ab initio calculations. For the neutral triatomic systems, the lowest energy structure is found to be triangular. The relative stability of the 2A1 and 2B2 structures can be predicted simply by knowing the constituent diatomic bond distances and atomic electron affinities (EAs). The lowest energy structure is linear for all of the negative ions. For anionic clusters containing Au, the Au atom(s) preferentially occupy the terminal position(s). The EAs of the heteronuclear systems can be predicted relatively accurately from a weighted average of the corresponding homonuclear systems. Although the theoretical EAs are systematically too small, accurate predictions for the EAs of the triatomics are obtained by uniformly scaling the ab initio results using the accurate experimental EA values available for the atoms and homonuclear diatomics.

Bauschlicher, Charles W., Jr.↗
NASA NTRSDate not supplied

Elucidation of Electrolyte Decomposition Pathways in Li--O2 Batteries using Ab Initio Computations

Li-O2 batteries offer exciting possibilities for energy-storage applications requiring higher energy content than current Li-ion battery technologies. However, development of commercial Li-O2 batteries faces many challenges including electrolyte decomposition in order to realize the full potential of this emerging technology. In this work, we present ab initio calculations that probe the role of reactive oxygen species including singlet oxygen, superoxide, and peroxides on electrolyte decomposition for potential organic electrolytes for use in Li-O2 batteries. We investigate the influence that the harsh electrochemical environment near the cathode plays in promoting electrolyte decomposition by using increasingly detailed computational models. Our results highlight existing problems with organic electrolytes for Li-O2 battery development and suggest potential design strategies for reducing electrolyte decomposition.

Li-O2 Batteries↗
NASA NTRS2010

Can Vertical Profiles of Tropospheric Methane on Titan Be Derived from Radio-Occultation Soundings?

The intensity of the received signal at Earth in the radio occultations of Titan is attenuated both by refractive defocusing and pressure-induced absorption from N2-N2 and CH4-N2 pairs. Because the absorption strength is different for the two sets of pairs, matching the retrieved absorptivity profile can in principle yield the vertical variation in gaseous methane in the troposphere. There are two factors that make this difficult. The first is the propagation of noise in the phase and amplitude of the received signal in the absorption retrieval. The phase data is first inverted to retrieve vertical profiles of refractivity, from which the refractive defocusing is calculated. This is then subtracted from the observed. intensity attenuation of the received signal to generate a profile of atmospheric absorption. The second problem is the uncertainty in the pressure-induced absorption coefficients. Laboratory data at radio wavelengths is only available near room temperature (see, e.g., [1] for N2-N2), and the extrapolation to the low temperatures in Titan's troposphere is not well established. Ab initio calculations by Borysow et al. [2, 3] provide absorption coefficients at low temperatures and long wavelengths, but their accuracy has come into question. We present examples from Cassini radio occultations of Titan to illustrate the difficulties. For methane mole fractions in the lower troposphere comparable to that inferred from the Huygens probe (approximately 0.05), it will be difficult to separate the contributions of N2-N2 collisions from those of N2-CH4, collisions to the retrieved absorption. However, higher concentrations of CH4 and/or a higher signal-to-noise ratio from a future uplink experiment could result in a successful separation of the two components. However, key to this are highly accurate estimates of the absorption from a combination of laboratory measurements at love temperatures and long wavelengths, and possibly improved theoretical calculations.

Flasar, F. Michael↗
DOE OSTI2025

Probing autoionization decay lifetimes of the 4d−16ℓ core-excited states in xenon using attosecond noncollinear four-wave-mixing spectroscopy

The decay of core-excited states is a sensitive probe of autoionization dynamics and correlation effects in many-electron systems, occurring on ultrafast timescales. Xenon, with its dense manifold of autoionizing resonances that can be coupled with near-infrared light, provides a platform to investigate these processes. In this work, the autoionization decay lifetimes of 4d-16ℓ (ℓ = s, p, d, …) core-excited states in xenon atoms are probed with extreme ultraviolet (XUV) attosecond noncollinear four-wave-mixing (FWM) spectroscopy. The 4d{5/2,3/2}-16p XUV-bright states (optically dipole-allowed) exhibit decay lifetimes of ∼6 fs, which is consistent with spectator-type decay. In contrast, the 4d{5/2,3/2}-16s and 4d{5/2,3/2}-16d XUV-dark states (optically dipole forbidden) show longer decay lifetimes of ∼20 fs. Photoionization calculations confirm that all core-hole states with 4d character should decay via spectator channels in ≤6 fs, suggesting that the apparently longer dark-state decay times arise from an alternative mechanism. A few-level simulation of the FWM process shows that the inclusion of a nearby, longer-lived dark state can mimic the experimental FWM signal, suggesting population cycling with a second electronic state with non-4d character. Ab initio calculations support the presence of such multielectron excited states in the 60-70 eV range. These results demonstrate that FWM signals can encode coupled-state dynamics when probing complex systems, highlighting the importance of combining theoretical and experimental approaches to disentangle accurate core-level decay pathways and lifetimes.

Puskar, Nicolette G↗
DOE OSTI2025

Defect Complexes in CrSBr Revealed Through Electron Microscopy and Deep Learning

Atomic defects underpin the properties of van der Waals materials, and their understanding is essential for advancing quantum and energy technologies. Scanning transmission electron microscopy is a powerful tool for defect identification in atomically thin materials, and extending it to multilayer and beam-sensitive materials would accelerate their exploration. Here, we establish a comprehensive defect library in a bilayer of the magnetic quasi-1D semiconductor CrSBr by combining atomic-resolution imaging, deep learning, and calculations. We apply a custom-developed machine learning work flow to detect, classify, and average point vacancy defects. This classification enables us to uncover several distinct Cr interstitial defect complexes, combined Cr and Br vacancy defect complexes, and lines of vacancy defects that extend over many unit cells. We show that their occurrence is in agreement with our computed structures and binding energy densities, reflecting the intriguing layer interlocked crystal structure of CrSBr. Our ab initio calculations show that the interstitial defect complexes give rise to highly localized electronic states. These states are of particular interest due to the reduced electronic dimensionality and magnetic properties of CrSBr and are, furthermore, predicted to be optically active. Our results broaden the scope of defect studies in challenging materials and reveal new defect types in bilayer CrSBr that can be extrapolated to the bulk and to over 20 materials belonging to the same FeOCl structural family.

deep learning↗
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 OSTI2025

Unconventional Anomalous Hall Effect in Hexagonal Polar Magnet Y 3 Co 8 Sn 4

A rare realization of unconventional anomalous Hall effect (UAHE) both below and above the magnetic transition temperature (T C ) in a hexagonal noncentrosymmetric magnet Y 3 Co 8 Sn 4 , using a combined experimental and ab initio calculations. Occurrence of such UAHE is mainly attributed to the reciprocal ($\mathscr{KS}$) topology (i.e., the presence of topological Weyl points at/near the Fermi level), along with some contribution from the topological magnetic texture (at low temperatures), as inferred from the measured field-dependent ac susceptibility. The effect of UAHE on the measured transport behavior however evolves differently with temperature above and below T C = 53 K, suggesting different physical mechanism responsible in the two phases. A unique planar ferrimagnetic ordering is found to be the most stable state with ab-plane as the easy plane below T C . The simulated net magnetization and the moment per Co atom agrees fairly well with the experimentally measured values. A reasonably large anomalous Hall conductivity (AHC) is also observed in both the phases (above and below T C ) of the present compound, which is again not so ubiquitous. In conclusion, the results underscore the family of R 3 Co 8 Sn 4 (R = rare earth) polar magnets as a compelling backdrop for exploring the synergy of magnetism and topological non-trivial electronic states, pivotal for spintronic applications.

Density functional theory calculations↗
DOE OSTI2024

A modular and extensible CHARMM-compatible model for all-atom simulation of polypeptoids

Peptoids (N-substituted glycines) are a class of sequence-defined synthetic peptidomimetic polymers with applications including drug delivery, catalysis, and biomimicry. Classical molecular simulations have been used to predict and understand the conformational dynamics of single chains and their self-assembly into morphologies including sheets, tubes, spheres, and fibrils. The CGenFF-NTOID model based on the CHARMM General Force Field has demonstrated success in accurate all-atom molecular modeling of peptoid structure and thermodynamics. Extension of this force field to new peptoid side chains has historically required reparameterization of side chain bonded interactions against ab initio data. This fitting protocol improves the accuracy of the force field but is also burdensome and precludes modular extensibility of the model to arbitrary peptoid sequences. In this work, we develop and demonstrate a Modular Side Chain CGenFF-NTOID (MoSiC-CGenFF-NTOID) as an extension of CGenFF-NTOID employing a modular decomposition of the peptoid backbone and side chain parameterizations, wherein arbitrary side chains within the large family of substituted methyl groups (i.e., –CH 3 , –CH 2 R, –CHRR', and –CRR'R") are directly ported from CGenFF. We validate this approach against ab initio calculations and experimental data to develop a MoSiC-CGenFF-NTOID model for all 20 natural amino acid side chains along with 13 commonly used synthetic side chains and present an extensible paradigm to efficiently determine whether a novel side chain can be directly incorporated into the model or whether refitting of the CGenFF parameters is warranted. We make the model freely available to the community along with a tool to perform automated initial structure generation.

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

The Role of Surface Hydroxyls in Dehydration and Dehydrogenation of Formic Acid on Fe 3 O 4 (001)

Understanding the role of surface structure and hydroxylation in catalytic reactions on metal oxide surfaces is important for developing a mechanistic insight into the complex interface processes. Here, we investigate the reactivity of formic acid on reconstructed Fe 3 O 4 (001) using a combination of X-ray photoelectron spectroscopy, infrared reflection absorption spectroscopy, temperature-programmed reaction spectroscopy, low energy electron diffraction, and electronic structure calculations. We find that formic acid initially dissociates at low temperatures (< 80 K) into bidentate formate and a hydroxyl up to an initial dosed coverage of two HCOOH per Fe 3 O 4 (001) unit cell. At higher temperatures (> 450 K), formate largely decomposes along the dehydration pathway, producing CO and H 2 O, with dehydrogenation to CO 2 being a minority side reaction. As a first step, water formation leads to surface oxygen extraction via the Mars-van Krevelen mechanism. Computational studies reveal formate embedded in oxygen vacancies as a key intermediate in the CO formation mechanism. CO formation proceeds via two reaction pathways with desorption that peaks at 530 K on the hydroxyl-rich surface and 560 K on the hydroxyl-deficient surface. Atomic hydrogen coadsorption experiments and ab initio calculations reveal that the presence of surface hydroxyls reduces the CO formation barrier. Furthermore, these results highlight the complex interactions between substrate and intermediate species occurring during reactions on metal oxide surfaces.

Adsorption↗
NASA NTRS1985

A measurement of the vibrational band strength for the upsilon sub 3 band of the HO2 radical

The HO2 radicals generated in a discharge-flow system were observed with tunable diode laser absorption in the P-branch of the nu(3) vibrationall band at 1080/cm. The observed line positions agree with those calculated from the molecular constants for the nu(3) bland obtained from a previous study using laser magnetic resonance spectroscopy. The band strength was determined by observing line center absoptions when HO2 is produced in the reaction F + H2O2 yields HO2 + HF (k1) with a measured concentration of atomic fluorine and excess hydrogen peroxide. F-atom concentrations are measured by diode laser absorption of the spin-orbit transition at 404/cm. The analysis accounts for HO2 losses due to the reactions of HO2 + HO2 yields H2O2 + O2 (k3) and F + HO2 yields HF + O2 (k4). The line strength for the 6(15) 7(16) F(1) transition is 2.9 x 10 to the 21st power sq cm/molecule/cm which corresponds to a nu(3) band strength of 34 +/- 9 sq/cm(STP atm). This value is a factor of 6 lower than previous ab initio calculations. These results will be useful in assessing the feasibility of atmospheric measurements of HO2 using infrared absorption techniques.

Zahniser, M. S.↗