Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ab initio data”

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 1,009 records · Page 56

Neutral weak-current two-body contributions in inclusive scattering from {sup 12}C

An {\it ab initio} calculation of the sum rules of the neutral weak response functions in $^{12}$C is reported, based on a realistic Hamiltonian, including two- and three-nucleon potentials, and on realistic currents, consisting of one- and two-body terms. We find that the sum rules of the response functions associated with the longitudinal and transverse components of the (space-like) neutral current are largest and that a significant portion ($\simeq 30$\%) of the calculated strength is due to two-body terms. This fact may have implications for the MiniBooNE and other neutrino quasi-elastic scattering data on nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

LANL Contributions to pyiron: An integrated development environment for computational materials science

To support and accelerate the development of simulation protocols in atomistic modelling, we introduce an integrated development environment (IDE) for computational materials science called pyiron (http://pyiron.org). The pyiron IDE combines a web based source code editor, a job management system for build automation, and a hierarchical data management solution. The core components of the pyiron IDE are pyiron objects based on an abstract class, which links application structures such as atomistic structures, projects, jobs, simulation protocols and computing resources with persistent storage and an interactive user environment. The simulation protocols within the pyiron IDE are constructed using the Python programming language. To highlight key concepts of this tool as well as to demonstrate its ability to simplify the implementation and testing of simulation protocols we discuss two applications. In these examples we show how pyiron supports the whole life cycle of a typical simulation, seamlessly combines ab initio with empirical potential calculations, and how complex feedback loops can be implemented. While originally developed with focus on ab initio thermodynamics simulations, the concepts and implementation of pyiron are general thus allowing to employ it for a wide range of simulation topics.

Janssen, Jan↗

Ab initio Molecular Dynamics Assessment of Thermodynamic and Transport Properties in (K,Li)Cl and (K, Na)Cl Molten Salt Mixtures

Molten salt mixtures are integral part of highly important technological applications such as nuclear reactors. However, due to inherent difficulties associated with experiment at high temperatures and the intrinsic complexity of liquid-phase multi-component systems, understanding their properties at a molecular level remains a challenge. Here, we report on an ab initio molecular dynamics investigation on structural, electronic, transport, and thermal properties of two common molten salt mixtures, (K, Li)Cl and (K,Na)Cl, at five different compositions and three temperatures. Most of the properties were found to depend on both composition and temperature. While properties, like atomic charges, show additive behaviors, other properties, such as electrical conductivity, show considerable deviations from additivity. We shall show that the mixing of the molten salt mixtures is mainly driven by entropy, and that the KCl and LiCl mix better than KCl and NaCl. Our computational results are in general consistent with available experimental data. Comparison with available theoretical data is also provided.

Nguyen, Manh Thuong↗

Ab initio intermolecular potential energy surface and thermophysical properties of nitrous oxide

We present an analytical intermolecular potential energy surface (PES) for two rigid nitrous oxide (N{sub 2}O) molecules derived from high-level quantum-chemical ab initio calculations. Interaction energies for 2018 N{sub 2}O–N{sub 2}O configurations were computed utilizing the counterpoise-corrected supermolecular approach at the CCSD(T) level of theory using basis sets up to aug-cc-pVQZ supplemented with bond functions. A site-site potential function with seven sites per N{sub 2}O molecule was fitted to the pair interaction energies. We validated our PES by computing the second virial coefficient as well as shear viscosity and thermal conductivity in the dilute-gas limit. The values of these properties are substantiated by the best experimental data.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Relativistic corrections for lepton-nucleus scattering in the short-time approximation

Here, we present an approach for including relativistic corrections in lepton-nucleus scattering calculations within the short-time approximation (STA). Previous ab initio studies employed electromagnetic currents expanded in powers of 𝑞/𝑚, where 𝑞 is the momentum transfer and 𝑚 is the nucleon mass, restricting their validity to low-𝑞 kinematics. We adopt an expansion scheme that treats the initial nucleon momentum perturbatively while allowing for arbitrary momentum transfer, thereby extending the applicability of the STA to high-𝑞 regimes. Additionally, we incorporate a relativistic treatment of the two-nucleon final-state energies. Calculations for 3 He and 4 He inclusive electron-scattering cross sections show a substantial improvement over previous results, achieving good agreement with experimental data in the quasi-elastic region for both low- and high-momentum transfer.

Andreoli, Lorenzo [Old Dominion Univ., Norfolk, VA↗

Theoretical and experimental studies relevant to interpretation of auroral emissions

The accomplishments achieved over the past year are detailed with emphasis on the interpretation or auroral emissions and studies of potential spacecraft-induced contamination effects. Accordingly, the research was divided into two tasks. The first task is designed to add to the understanding of space vehicle induced external contamination. An experimental facility for simulation of the external environment for a spacecraft in low earth orbit was developed. The facility was used to make laboratory measurements of important phenomena required for improving the understanding of the space vehicle induced external environment and its effect on measurement of auroral emissions from space-based platforms. A workshop was sponsored to provide a forum for presentation of the latest research by nationally recognized experts on space vehicle contamination and to discuss the impact of this research on future missions involving space-based platforms. The second task is to add an ab initio auroral calculation to the extant ionospheric/thermospheric global modeling capabilities. Once the addition of the code was complete, the combined model was to be used to compare the relative intensities and behavior of various emission sources (dayglow, aurora, etc.). Such studies are essential to an understanding of the types of vacuum ultraviolet (VUV) auroral images which are expected to be available within two years with the successful deployment of the Ultraviolet Imager (UVI) on the ISTP POLAR spacecraft. In anticipation of this, the second task includes support for meetings of the science working group for the UVI to discuss operational and data analysis needs. Taken together, the proposed tasks outline a course of study designed to make significant contributions to the field of space-based auroral imaging.

Keffer, Charles E.↗

Probabilistic Modeling of Aircraft Trajectories for Dynamic Separation Volumes

With a proliferation of new and unconventional vehicles and operations expected in the future, the ab initio airspace design will require new approaches to trajectory prediction for separation assurance and other air traffic management functions. This paper presents an approach to probabilistic modeling of the trajectory of an aircraft when its intent is unknown. The approach uses a set of feature functions to constrain a maximum entropy probability distribution based on a set of observed aircraft trajectories. This model can be used to sample new aircraft trajectories to form an ensemble reflecting the variability in an aircraft's intent. The model learning process ensures that the variability in this ensemble reflects the behavior observed in the original data set. Computational examples are presented.

Lewis, Timothy A.↗

A theoretical kinetic study of ĊH 3 + ṄH 2 : From electronic structure to NH 3 /CH 4 combustion modelling implications

Carbon–nitrogen interaction reactions play an important role in governing the reactivity of ammonia blended fuels. However, there remains uncertainties regarding their detailed reaction pathways and rate constants, hampering the development of high-fidelity chemical kinetic models. In this study, the kinetics of ĊH 3 + ṄH 2 , a key C–N interaction reaction in ammonia/methane blend combustion have been investigated. The potential energy surface has been explored using the high-level ANL0F method, yielding highly accurate stationary point energies that agree with ATcT values within 0.1 kcal mol –1 . Variable reaction coordinate transition state theory is used to treat the barrierless association and decomposition reaction channels, based on directly sampled radical-radical interaction energies at the CASPT2-F12(2e,2o)/cc-pVTZ-F12 level of theory. The minimum transitional mode numbers of states obtained are then coupled with the RRKM/master equation to calculate temperature- and pressure-dependent rate constants. Our a priori calculations capture available experimental measurements from the literature very well. The calculated rate constants have been incorporated into an NH 3 /CH 4 chemical kinetic model currently under development at the University of Galway. The effect of the updated kinetic data for ĊH 3 + ṄH 2 on model predicted NH 3 /CH 4 fuel reactivity is elucidated.

ab initio↗

Modelling the local atomic structure of molybdenum in nuclear waste glasses with ab initio molecular dynamics simulations

The nature of chemical bonding of molybdenum in high level nuclear waste glasses has been elucidated by ab initio molecular dynamics simulations. Two compositions, (SiO2)57.5 – (B2O3)10 – (Na2O)15 – (CaO)15 – (MoO3)2.5 and (SiO2)57.3 – (B2O3)20 – (Na2O)6.8 – (Li2O)13.4 – (MoO3)2.5 , were considered in order to investigate the effect of ionic and covalent components on the glass structure and the formation of the crystallisation precursors (Na2MoO4 and CaMoO4). The coordination environments of Mo cations and the corresponding bond lengths calculated from our model are in excellent agreement with experimental observations. The analysis of the first coordination shell reveals two different types of molybdenum host matrix bonds in the lithium sodium borosilicate glass. Based on the structural data and the bond valence model, we demonstrate that the Mo cation can be found in a redox state and the molybdate tetrahedron can be connected with the borosilicate network in a way that inhibits the formation of crystalline molybdates. These results significantly extend our understanding of bonding in Mo-containing nuclear waste glasses and demonstrate that tailoring the glass composition to specific heavy metal constituents can facilitate incorporation of heavy metals at high concentrations. K.K. was supported through the Impact Studentship scheme at UCL co-funded by the IHI Corporation and UCL. P.V.S. thanks the Royal Society, which supported preliminary work on this project, and the Laboratory Directed Research and Development program at PNNL, a multiprogram national laboratory operated by Battelle for the U.S. Department of Energy. Via our membership of the UK's HEC Materials Chemistry Consortium, which is funded by EPSRC (EP/L000202), this work used the ARCHER UK National Supercomputing Service (http://www.archer.ac.uk).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Microscopic scattering approach to in-gap states: Cr adatoms on superconducting 𝛽−Bi 2 ⁢Pd

We develop a microscopic scattering formalism to describe Yu-Shiba-Rusinov (YSR) states due to a single Cr adatom on the Bi-terminated surface of 𝛽−Bi 2 ⁢Pd, by combining ab initio Wannier functions with a real-space Green's function approach in the Bogoliubov–de Gennes formalism. Our framework reproduces key scanning tunneling spectroscopy features, including a single particle-hole asymmetric YSR peak and isotropic 𝑑⁢𝐼/𝑑⁢𝑉 maps around the impurity. Decomposing the YSR states reveals contributions from four nearly degenerate 𝐶 4⁢𝑣 representations, with energy broadening masking their individual signatures. Spin-orbit coupling induces partial spin polarization, while the spatial asymmetry between particle and hole components arises from Cr 𝑑−Bi 𝑝 hybridization. These results highlight the importance of realistic band structures and microscopic modeling for interpreting STM data and provide a foundation for studying impurity chains hosting topological excitations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Ca 2 MnO 4 structural path: Following the negative thermal expansion at the local scale

The oxygen octahedral rotations in Ca 2 MnO 4 , the first member of the CaO(CaMnO 3 )n Ruddlesden-Popper family, is probed through a set of complementary techniques, including temperature-dependent neutron and x-ray diffraction, combined with local probe studies and ab initio calculations. Here we demonstrate the enhancement of the uniaxial negative thermal expansion coefficient from -1.26 ± 0.25 to -21 ± 1.8 ppm/K at the second order I4 1 /acd to I4/mmm structural phase transition, providing direct evidence for the corkscrew atomic mechanism. We establish, also, that the predicted I4/mmm high symmetry is attained around 1050 K. At lower temperatures, within the 10–1000 K temperature range, our first-principles calculations and detailed analysis of the Ca local environment reveals that the reported Aba2 structural phase, coexisting with the I4 1 /acd one, cannot describe correctly this compound. On the other hand, our data allow for the coexistence of the locally identical I4 1 /acd and Acam structural phases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Calibration and validation of the foundation for a multiphase strength model for tin

In this work, the Common Model of Multi-phase Strength and Equation of State (CMMP) model was applied to tin. Specifically, calibrations of the strength-specific elements of the CMMP foundation were developed with a combination of experiments and theory, and then the model was validated experimentally. The first element of the foundation is a 10 multi-phase analytic treatment of the melt temperature and the shear modulus for the solid phases. These models were parameterized for each phase based on ab initio calculations using the software VASP (Vienna Ab initio Simulations Package) based on density functional theory (DFT). The shear modulus model for the ambient β phase was validated with ultrasonic sound speed measurements as a function of pressure and temperature. The second element of the foundation is a viscoplastic strength model for the β phase, upon which strength for inaccessible higher-pressure phases can be scaled as necessary. The stress-strain response of tin was measured at strain rates of 10 -3 to 3 x 10 3 s -1 and temperatures ranging from 87 to 373 K. The Preston-Tonks-Wallace (PTW) strength model was fit to that data using Bayesian model calibration. For validation, six forward and two reverse Taylor impact experiments were performed at different velocities to measure large plastic deformation of tin at strain rates up to 10 5 s -1 . The PTW model accurately predicted the deformed shapes of the cylinders, with modest discrepancies attributed to the inability 20 of PTW to capture the effects of twinning and dynamic recrystallization. Some material in the simulations of higher velocity Taylor cylinders reached the melting temperature, thus testing the multiphase model because of the presence of a second phase, the liquid. In simulations using a traditional modeling approach, the abrupt reduction of strength upon melt resulted in poor predictions of the deformed shape and non-physical temperatures. With CMMP, the most deformed material points evolved gradually to a mixed solid-liquid but never fully liquid state, never fully lost strength, 25 stayed at the melt temperature as the latent heat of fusion was absorbed, and predicted the deformed shape well.

36 MATERIALS SCIENCE↗

Microscopic Scattering Approach to In-Gap States

We develop a microscopic scattering formalism to describe Yu-Shiba-Rusinov (YSR) states due to a single Cr adatom on the Bi-terminated surface of beta Bi2Pd, by combining ab initio Wannier functions with a real-space Green's function approach in the Bogoliubov-de Gennes formalism[1]. Our framework reproduces key scanning tunneling spectroscopy features, including a single particle-hole asymmetric YSR peak and isotropic dIdV maps around the impurity. Decomposing the YSR states reveals contributions from four nearly degenerate C4v representations, with energy broadening masking their individual signatures. Spin-orbit coupling induces partial spin polarization, while the spatial asymmetry between particle and hole components arises from Cr d-Bi p hybridization. These results highlight the importance of realistic band structures and microscopic modeling for interpreting STM data for magnetic in-gap states on superconductors. Further advances examining layered 2D material surfaces, such as NbSe2, will be described[2]. For this system the superconducting properties are obtained from a full anisotropic Eliashberg calculation of the superconducting order parameter along with the charge density wave gap. Additional features associated with proposals to measure the dynamics of these individual YSR states will be presented. [1] arXiv:2507.08740 [2] arXiv:2507.11856

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

BH{sub 2} revisited: New, extensive measurements of laser-induced fluorescence transitions and ab initio calculations of near-spectroscopic accuracy

The spectroscopy of gas phase BH{sub 2} has not been explored experimentally since the pioneering study of Herzberg and Johns in 1967. In the present work, laser-induced fluorescence (LIF) spectra of the A{sup ~} {sup 2}B{sub 1}(Π{sub u})−X{sup ~2}A{sub 1} band system of {sup 11}BH{sub 2}, {sup 10}BH{sub 2}, {sup 11}BD{sub 2}, and {sup 10}BD{sub 2} have been observed for the first time. The free radicals were “synthesized” by an electric discharge through a precursor mixture of 0.5% diborane (B{sub 2}H{sub 6} or B{sub 2}D{sub 6}) in high pressure argon at the exit of a pulsed valve. A total of 67 LIF bands have been measured and rotationally analyzed, 62 of them previously unobserved. These include transitions to a wide variety of excited state bending levels, to several stretch-bend combination levels, and to three ground state levels which gain intensity through Renner-Teller coupling to nearby excited state levels. As an aid to vibronic assignment of the spectra, very high level hybrid ab initio potential energy surfaces were built starting from the coupled cluster singles and doubles with perturbative triples (CCSD(T))/aug-cc-pV5Z level of theory for this seven-electron system. In an effort to obtain the highest possible accuracy, the potentials were corrected for core correlation, extrapolation to the complete basis set limit, electron correlation beyond CCSD(T), and diagonal Born-Oppenheimer effects. The spin-rovibronic states of the various isotopologues of BH{sub 2} were calculated for energies up to 22 000 cm{sup −1} above the X{sup ~} (000) level without any empirical adjustment of the potentials or fitting to experimental data. The agreement with the new LIF data is excellent, approaching near-spectroscopic accuracy (a few cm{sup −1}) and has allowed us to understand the complicated spin-rovibronic energy level structure even in the region of strong Renner-Teller resonances.

74 ATOMIC AND MOLECULAR PHYSICS↗

Reuniting crystallography with real space: Ab initio structure elucidation with 4D-STEM

Structure elucidation via single-crystal methods has historically lacked experimental access to real-space information, instead relying exclusively on diffraction-space measurements of Bragg reflections. Here we exploit the dual-space imaging power of 4D scanning transmission electron microscopy to meaningfully integrate real-space information into the crystallographic workflow. We show that virtual apertures assembled by segmentation of high-angle annular dark-field images enable i) pixel-by-pixel separation of coherent Bragg signal from clusters of closely spaced nanocrystals and ii) selective extraction of integrated intensities from thinner subregions of individual specimens, facilitating retroactive tuning of multiple scattering artifacts. This strategy empowers us to simply pick and choose whichever nanoscale regions of interest generate the highest-quality diffraction patterns, allowing us to solve several independent structures of the metal-organic framework UiO-66 from specimens whose agglomerated morphology proved intractable for conventional microcrystal electron diffraction. Our method is compatible with both rotational and serial approaches to data processing, ultimately divulging the first scanning nanobeam electron diffraction structures determined by direct methods at subangstrom resolution.

Saha, Ambarneil↗

Anomalous thermoelectricity at the two-dimensional structural transition of SnSe monolayers

The thermoelectric figure of merit ZT comprises electronic and vibrational contributions that drastically change across phase transitions, and the most common theoretical ab initio approach to thermoelectricity fails to describe the evolution of ZT across finite-temperature structural transitions in its entirety. Furthermore, while the thermoelectric behavior of bulk SnSe has been extensively studied, SnSe monolayers have been experimentally realized only recently, and the existent prediction of thermoelectricity on this two-dimensional material is unreliable because it misses its structural transition altogether. SnSe monolayers (and similar GeS, GeSe, GeTe, SnS, and SnTe monolayers) experience a temperature-induced two-dimensional Pnm2 1 → P4/nmm structural transition precipitated by the softening of vibrational modes, and we describe their thermoelectric properties across the phase transition, using molecular dynamics data to inform both electronic and vibrational coefficients directly and within the same footing. Similar to recent experimental observations pointing to an overestimated ZT past the transition temperature in bulk SnSe, we find a smaller ZT on SnSe monolayers when compared to its value predicted by the standard paradigm, due to the dramatic changes in the electrical conductivity and lattice thermal conductivity as the structural transition ensues. Here, the process described here lends a strong focus to both the vibrational and electronic evolutions throughout the structural transition, and it applies to thermoelectric materials undergoing thermally driven solid-to-solid structural phase transitions in one, two, and three dimensions.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

P SI4 1.4: Open-source software for high-throughput quantum chemistry

PSI4 is a free and open-source ab initio electronic structure program providing implementations of Hartree–Fock, density functional theory, many-body perturbation theory, configuration interaction, density cumulant theory, symmetry-adapted perturbation theory, and coupled-cluster theory. Most of the methods are quite efficient, thanks to density fitting and multi-core parallelism. The program is a hybrid of C++ and Python, and calculations may be run with very simple text files or using the Python API, facilitating post-processing and complex workflows; method developers also have access to most of PSI4 's core functionalities via Python. Job specification may be passed using The Molecular Sciences Software Institute (MolSSI) QCSCHEMA data format, facilitating interoperability. We announce a rewrite of our top-level computation driver, and concomitant adoption of the MolSSI QCARCHIVE INFRASTRUCTURE project, makes the latest version of PSI4 well suited to distributed computation of large numbers of independent tasks. The project has fostered the development of independent software components that may be reused in other quantum chemistry programs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ground-state properties of the Hubbard model in one and two dimensions from the Gutzwiller conjugate gradient minimization theory

We introduce Gutzwiller conjugate gradient minimization (GCGM) theory, an ab initio quantum many-body theory for computing the ground-state properties of infinite systems. GCGM uses the Gutzwiller wave function but does not use the commonly adopted Gutzwiller approximation (GA), which is a major source of inaccuracy. Instead, the theory uses an approximation that is based on the occupation probability of the on-site configurations, rather than approximations that decouple the site-site correlations as used in the GA. We test the theory in the one-dimensional and two-dimensional Hubbard models at various electron densities and find that GCGM reproduces energies and double occupancies in reasonable agreement with benchmark data at a very small computational cost.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗