Engineering PapersSearch

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 37 records · Page 2

Ab initio calculation of the X 1 Sigma + state of CsH

Stevens et al. (1981) considered CsH as a two electron problem, but simulated the Cs core electrons by an empirical pseudopotential following the work by Bardsley (1970). These potentials, since they are derived empirically, include the relativistic effects and the atomic core-valence interactions. However, molecular core-core interactions are not accounted for. Stevens et al. obtained an R(e) value which is too small. This result was attributed to uncorrected core (proton)-core interactions in the molecule. The present investigation is concerned with a resolution of the discrepancies between the results of earlier studies. The X 1 Sigma + ground state of CsH is recalculated. The calculation employs a nine valence electron relativistic effective core potential (RECP) for Cs. The bonding in CsH is found to involve a Cs(6s)-H(1s) bond but with a significant ionic (Cs+H-) component.

Laskowski, B. C.

Zeeman effect of P II.

An investigation was conducted of the Zeeman effect of PII because it provides unique J-value assignments together with a sensitive check on the accuracy of wave functions obtained from least-squares fitting of energy levels. About forty new lines were added to the spectrum, and g values were obtained for 76 levels of PII. By incorporating this new experimental material in a systematic program of least-squares calculations, supplemented by ab initio calculations, it was possible to complete the description of the 3p5p and 3p5d configurations.

Li, H.

Edge dislocation mediated anomalous charge transfer in face centered cubic high entropy alloys

The charge transfer in alloys dictates their structural and functional properties. The presence of the line defect can influence the charge transfer characteristics in alloys. The edge dislocation-volume misfit interaction in high entropy alloys is a critical area of interest concerning the development of accurate solid solution strengthening models for these materials. The evolution of the volume misfit of atoms is either studied empirically using Vegard's law or mechanistically using \emph{ab initio} calculations, without considering the effect of edge dislocation on charge transfer explicitly. In this work, we show the influence of edge dislocation on the charge transfer in CoCrFeMnNi, CoCrNi and CoNi using large-scale \emph{ab initio} calculations. The anomalous charge transfer characteristics, in terms of deviations from electronegativity trends, are demonstrated, which are dictated by electronegativity equalisation rather than pair-atom interactions. The deviation of atomic volume in the compressive and tensile stress regions of the edge dislocation is rationalised in terms of anomalous magneto-volume fluctuations in such alloys.

Ghosh, Swarnava [ORNL] (ORCID:0000000338005264)

Reversible hydrogen storage in multilayer graphane: Lattice dynamics, compressibility, and heat capacity studies

Multilayer graphane (hydride of graphite) is a crystalline hydrocarbon of composition CH, which can be synthesized from graphite and molecular hydrogen at pressures above 2GPa [V.E. Antonov et al. Carbon 100 (2016) 465]. Using X-ray diffraction, this compound was tentatively identified as the “graphane II” phase of 3D-graphane predicted by ab initio calculations [X.-D. Wen et al. PNAS 108 (2011) 6833] and consisting of layers of 2D-graphane in the “chair” conformation. When heated in a vacuum, the compound does not form any intermediate hydrocarbons and reversibly decomposes back into graphite and hydrogen at 770–920 K. In the present work, almost single-phase samples of graphite hydride and deuteride were synthesized at 7.4 GPa and 870 K. Their investigation by inelastic neutron scattering supplemented by ab initio calculations gave spectra g(E) of the phonon density of states with a gap of about 15 meV at approx. 100 meV, which is a unique identifier for the chair form of graphane. The equation of state V(P) of the hydride was studied at room temperature and hydrogen pressures up to 53 GPa by synchrotron X-ray diffraction in a diamond anvil cell. Further, the graphane II phase did not react with the surrounding hydrogen and did not undergo any phase transformations upon the compression and after heating to 1500 K at 53GPa. The high thermal and pressure stability of this exotic phase makes it an important part of the C–H system. The obtained g(E) spectra of graphite hydride and deuteride were used to calculate temperature dependences of their heat capacity. Measurements of the heat capacity at temperatures 120–673 K confirmed the good accuracy of these calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

From Electronic Structure to Ion Transport: Photoelectron Spectroscopy and Molecular Dynamics Simulations Reveal the Role of Anions in Lithium Battery Electrolytes

Electrolyte anions are pivotal for lithium battery performance, yet their fundamental electronic structural properties are not well understood. In this work, we employ a combination of negative-ion photoelectron spectroscopy (NIPES), ab initio calculations, and molecular dynamics (MD) simulations to investigate the electronic structures of three representative electrolyte anions. This multiscale approach enables us to elucidate how their intrinsic electronic properties govern anion–solvent interactions in gas-phase clusters, as well as lithium-ion (Li + ) solvation structures and ion transport behavior in the condensed phase. NIPES reveals that difluoro(oxalato)borate (DFOB – ), bis(fluorosulfonyl)imide (FSI – ), and bis(oxalato)borate (BOB – ) all exhibit high electron binding energies, with vertical/adiabatic detachment energies increasing from DFOB – (6.09/5.70 eV) to FSI – (6.80/6.10 eV) to BOB – (6.82/6.40 eV), correlating with enhanced oxidation stability. Ab initio calculations reveal that DFOB – /FSI – –solvent complexes bind Li + ∼ 10 kcal/mol stronger than BOB – series, aligning with the strength of a Li + –anion model. DFOB – exhibits pronounced charge localization on both oxygen and fluorine atoms, enabling their involvement in Li + coordination. In contrast, fluorine atoms in FSI – are largely electron-depleted and remain excluded from direct Li + binding. MD simulations further demonstrate that LiDFOB and LiFSI systems exhibit Li + diffusion coefficients three and five times higher than those of LiBOB across four common solvents. Notably, LiFSI salt in acetonitrile (AN) exhibits the fastest Li + diffusion among 12 electrolyte systems, highlighting the synergistic effect of FSI – and AN in promoting ion mobility. In conclusion, these findings provide a molecular-level understanding of the critical roles of anion and its microsolvation in optimizing Li + diffusion dynamics, once again emphasizing the positioning of FSI – and DFOB – as prime candidates for next-generation electrolytes.

25 ENERGY STORAGE

The rotational spectra of HOCO/plus/, HOCN, HN3, and HNCO from quantum mechanical calculations

Ab initio molecular orbital theory has been used to determine the equilibrium geometries, rotational constants, and rotational spectra of four isoelectronic molecules. Two of these, HOCO(plus) and HOCN, are candidate interstellar molecules. The other two, HNCO and HN3, have known rotational constants. Theoretical rotational constants and spectra for the two unknown species were corrected with the mean experimental to theoretical ratios from the two known species. This procedure resulted in predicted frequencies of 83.75 plus or minus 0.2 GHz for the 4(04) to 3(03) transition in HOCN and 85.08 plus or minus 0.2 GHz for the same transition in HOCO(plus). These are the central lines of triplets whose other members are the 4(14) to 3(13) and 4(13) to 3(12) transitions. The triplet splittings were predicted to be 0.36 plus or minus 0.01 GHz for HOCN and 0.33 plus or minus 0.01 GHz for HOCO(plus). These results indicate that HOCO(plus) is a better candidate for the source of a series of lines reported by Thaddeus, Guelin, and Linke than is HOCN.

Defrees, D. J.

Interpretation of experimental differential elastic scattering cross section for H/+/ + Ne.

The interatomic V(r) found by potential-model calculations is compared to an ab initio calculation of the intermolecular potential for NeH+ due to Peyerimhoff (1965). The reason for a rather significant difference in the two results is discussed, and a new method due to Remler (1971) involving Regge poles is applied. The method can be used to calculate the differential cross section, where the starting point for such calculations is based on intuition in light of the classical deflection function.

Bobbio, S. M.

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

Electric dipole polarizability of 58 Ni

The electric dipole strength distribution in 58 Ni between 6 and 20 MeV has been determined from proton inelastic scattering experiments at very forward angles at RCNP, Osaka. The experimental data are rather well reproduced by quasiparticle random-phase approximation calculations including vibration coupling, despite a mild dependence on the adopted Skyrme interaction. They allow an estimate of the experimentally inaccessible high-energy contribution above 20 MeV, leading to an electric dipole polarizability α D ⁡ ( 58 Ni) = 3.48 (31)⁢ fm 3 . This serves as a test case for recent extensions of coupled-cluster calculations with chiral effective field theory interactions to nuclei with two nucleons on top of a closed-shell system.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS

Understanding Adsorption and Reactions at Aqueous Oxide Interfaces with Neural Network Potential Molecular Dynamics

Chemical processes at metal oxide−water interfaces are of central importance in geochemistry, biology, and energy technologies. A better understanding of these processes would allow us to make a significant step toward optimizing and controlling them, which could in turn lead to broader impacts. Computational modeling is indispensable to accomplishing this task because complexity and disorder often make it difficult to extract atomistic information from experiments. Balancing computational cost and accuracy, simulation schemes based on efficient machine learning representations of the potential energy surface (PES) predicted by ab initio calculations have become increasingly popular over the past decade. In particular, several studies have demonstrated the ability of machine learning models to accurately reproduce the complex ab initio PESs of aqueous oxide interfaces, allowing simulations of systems and processes that are not accessible with ab initio methods. In this Account, we review our recent efforts to understand adsorption processes and reactions at aqueous oxide interfaces using deep potential molecular dynamics (DPMD), a simulation scheme employing deep neural networks (DNNs), which has proven to be quite successful in accurately describing many different systems in the condensed phase. After summarizing the DPMD methodology, we first review our work on the acid−base chemistry of oxide surfaces in contact with water, a fundamental characteristic that controls proton transfer and surface charge at the interface. We focus on the aqueous interface of rutile IrO 2 , an oxide material thus far considered the best catalyst for the oxygen evolution reaction (OER). We show that this interface is characterized by a large fraction of dissociated water and a strong Brønsted acidity of the surface sites, in good agreement with the experimentally measured value of the point of zero proton charge. In our second example, we investigate how the adsorption of organic species from ambient air or water affects the structure and wettability of the aqueous interfaces of TiO 2 , a prototypical photocatalytic material. This is a question that is relevant to understanding the UV-induced hydrophilicity of TiO 2 surfaces, a property at the basis of self-cleaning windows and related applications. Specifically focusing on formic and acetic acids, the two most common atmospheric organic acids, our simulations reveal that these acids control the wettability of TiO 2 largely through acid−base chemistry at the interface rather than chemisorption on the oxide surface, a finding that could help improve the design of self-cleaning surfaces and photocatalytic devices. Finally, we review our recent study of methanol at TiO 2 −water interfaces, a system whose interest is largely motivated by the role of methanol in enhancing photocatalytic hydrogen evolution on TiO 2 . Our simulations provide mechanistic insights into the coupled roles of the organic adsorbate and water at the TiO 2 interface, with implications for how methanol enhances the activity of H 2 evolution.

adsorption

Hyperfine-to-rotational energy transfer in ultracold atom–molecule collisions of Rb and KRb

Energy transfer between different mechanical degrees of freedom in atom–molecule collisions has been studied and largely understood. However, systems involving spins remain less explored. Here, in this study, we directly observed energy transfer from atomic hyperfine to molecular rotation in the 87 Rb ($| {F}_{a},{M}_{{F}_{a}}\rangle =| 2,2\rangle$) + 40 K 87 Rb (X 1 Σ + , rotational state N = 0) ⟶ Rb ($| 1,1\rangle$) + KRb ( N = 0, 1, 2) collision with state-to-state precision. We also performed quantum scattering calculations that rigorously included the coupling between spin and rotational degrees of freedom at short range under the assumption of rigid-rotor KRb monomers moving along a single potential energy surface. The calculated product rotational state distribution deviates from the observations even after extensive tuning of the atom–molecule potential energy surface. In addition, our ab initio calculations indicate that spin–rotation coupling is enhanced close to a conical intersection that is energetically accessible at short range. This, together with the deviation, suggests that vibrational degrees of freedom and conical intersections play an important part in the coupling. Our observations confirm that spin is coupled to mechanical rotation at short range and establish a benchmark for future theoretical studies.

chemical physics

RIKEN TRIP Magnets Database

This dataset contains ab-initio calculation results for the temperature-dependent anomalous Hall conductivity, the anomalous Nernst effect, and the Seebeck coefficient. All calculations are based on ab-inito Quantum Espresso (PWSCF v.6.3) + Wannier90 (v.3.0.0). The dependence on carrier doping is also calculated. For all calculations a ferromagnetic order has been assumed, which might not correspond to the true ground state of the system. Tabulated values for the magnetic moments and essential input files for Quantum Espresso are available for download as attachments. This project has been supported by the RIKEN Transformative Research Innovation Platform (TRIP), Use Case: Many-body Electron Systems.

36 MATERIALS SCIENCE

GPR_calculator: An on-the-fly surrogate model to accelerate massive nudged elastic band calculations

We present GPR_calculator, a package based on Python and C++ programming languages to build an on-the-fly surrogate model using Gaussian Process Regression (GPR) to approximate computationally expensive electronic structure calculations. The key idea is to dynamically train a GPR model during the simulation that can accurately predict energies and forces with uncertainty quantification. When the uncertainty is high, the costly electronic structure calculation is performed to obtain the ground truth data, which is then used to update the GPR model. To illustrate the effectiveness of GPR_calculator, we demonstrate its application in Nudged Elastic Band (NEB) simulations of surface diffusion and reactions, achieving 3-10 times acceleration compared to pure ab initio calculations. The source code is available at https://github.com/MaterSim/GPR_calculator.

Gaussian process regression

Probing the 𝑗 dependence of angular distributions and 𝑁=20 shell rigidity via the 36 S (𝑝,𝑑) 35 S reaction

An investigation of the N=20 36 S nucleus has been performed through a detailed study of the 3 6S (p,d)3 5S neutron-removal reaction, employing a 66-MeV proton beam at iThemba Laboratory for Accelerator Based Sciences and an innovative target design. Using the high-resolution K = 600 magnetic spectrometer, 98 states in 3 5S were identified up to 16-MeV excitation energy, including 47 previously unobserved states. Angular distributions and spectroscopic factors, including isobaric analog-state contributions, were extracted for 81 levels. A pronounced j dependence in the angular distributions of ℓ=2 states provides refined insights into the spin-orbit splitting. Finite-range adiabatic distorted-wave approximation calculations qualitatively reproduce the observed j dependence. Comparisons of the measured 1d 5/2 spectroscopic strength distribution with large-scale shell-model and ab initio calculations show good agreement overall, and the robustness of the N=20 shell closure in 36 S is confirmed when comparing the relatively low fp orbital occupancies in 40 Ca and 3 6S across the Fermi surface. This study underscores the utility of neutron-removal reactions in probing nuclear structure and the Fermi surface of sd nuclei and beyond. The findings advance our understanding of shell evolution and offer constraining data for theoretical models.

20 ≤ A ≤ 38

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

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

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.