Search NASA⌕ Search

SEARCH · Search NASA

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Ab-initio informed cluster dynamics simulation of self- and Xe diffusivity in uranium mononitride under irradiation

Uranium mononitride (UN) is one of the ceramic nuclear fuel alternatives to oxide fuels considered for light water reactor and advanced reactor designs, as it presents significant advantages such as high uranium density (better economics) and high thermal conductivity and melting point (increased safety). Self- and fission gas diffusivities need to be better understood, given that they influence key fuel performance phenomena like swelling and fission gas release. Recently, radiation enhanced diffusivity was investigated in UN by means of cluster dynamics simulations relying on empirical potential-based parameterizations, the reliability of which highly depends on the interatomic potential accuracy. Here, in this work, we refine this approach by determining, using ab-initio calculations, the properties of defect clusters containing vacancies, self-interstitials and Xe impurities. We also consider larger clusters than previous studies. The obtained dataset (formation enthalpies, entropies, and migration barriers) is used to parameterize a cluster dynamics model of mobile clusters, and to calculate the defect cluster concentrations under irradiation. This gives us access to the radiation enhanced self- and fission gas diffusivities. Although the resulting diffusivities are close to the values reported in the literature, we find important qualitative differences in the diffusion mechanisms. Capturing the correct mechanisms is crucial to properly describe the chemistry and fission rate dependence of the model.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Real-space visualization of short-range antiferromagnetic correlations in a magnetically enhanced thermoelectric

Short-range magnetic correlations can significantly increase the thermopower of magnetic semiconductors, representing a noteworthy development in the decades-long effort to develop high-performance thermoelectric materials. Here, we reveal the nature of the thermopower-enhancing magnetic correlations in the antiferromagnetic semiconductor MnTe. Using magnetic pair distribution function analysis of neutron scattering data, we obtain a detailed, real-space view of robust, nanometer-scale, antiferromagnetic correlations that persist into the paramagnetic phase above the Neel temperature $T_N$ = 307 K. In this work, the magnetic correlation length in the paramagnetic state is significantly longer along the crystallographic c axis than within the ab plane, pointing to anisotropic magnetic interactions. Ab initio calculations of the spin-spin correlations using density functional theory in the disordered local moment approach reproduce this result with quantitative accuracy. These findings constitute the first real-space picture of short-range spin correlations in a magnetically enhanced thermoelectric and inform future efforts to optimize thermoelectric performance by magnetic means.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Predicting thermophysical properties of molten salts in the MgCl 2 -NaCl-KCl-LiCl system with a shell-model potential

Ternary eutectic salts composed of MgCl 2 , NaCl, and KCl, referred to as MNK salts, have recently emerged as promising candidates as high-temperature heat transfer fluids and thermal energy storage media. Here we performed classical molecular dynamics (MD) simulations to predict the densities, specific heat capacities, viscosities, and ionic self-diffusivities for MNK salts over a wide temperature range. The impact of LiCl additive on their thermophysical properties was also investigated. To capture the electronic polarization of Cl anions by neighboring cations, we developed a novel shell-model potential using the force-matching method and a dataset of ab initio calculated interatomic forces. Our extensive MD simulations predict structure and properties for pure salts and binary/ternary salt mixtures in the MgCl 2 -NaCl-KCl-LiCl system in overall good agreement with available experimental and theoretical data, which corroborates the accuracy and reliability of our developed potential.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Layered topological semimetal GaGeTe: New polytype with non-centrosymmetric structure

GaGeTe is a layered van der Waals material composed of germanene and GaTe sublayers that has been recently predicted to be a basic Z 2 topological semimetal. To date, only one polytype of GaGeTe is known with trigonal centrosymmetric structure (α phase, space group R-3m, No. 166). Here we show that as-grown samples of GaGeTe show traces of at least another polytype with hexagonal non-centrosymmetric structure (β phase, space group P6 3 mc, No. 186). Moreover, we suggest that another bulk hexagonal polytype (γ phase, space group P-3m1, No. 164) could also be found near room conditions. Both α and β polytypes have been identified and characterized by means of X-ray diffraction and Raman scattering measurements with the support of ab initio calculations. We provide the vibrational properties of both polytypes and show that the Raman spectrum reported for GaGeTe almost forty years ago and attributed to the α phase, was, in fact, that of the secondary β phase. Additionally, we show that a Fermi resonance occurs in α-GaGeTe under non-resonant excitation conditions, but not under resonant excitation conditions. Theoretical calculations show that bulk β-GaGeTe is a non-centrosymmetric weak topological semimetal with even smaller lattice thermal conductivity than centrosymmetric bulk α-GaGeTe. In perspective, our work paves the way for the control and engineering of GaGeTe polytypes to design and implement complex van der Waals heterostructures formed by a combination of centrosymmetric and non-centrosymmetric layers of up to three different polytypes in a single material, suitable for a number of fundamental studies and technological applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Anisotropic topological surface states in thin-film monoclinic Ag 2 Se

Silver chalcogenide systems have recently attracted significant attention due to their promising topological insulating properties. Here we conducted systematic low-temperature local scanning tunneling microscopy/spectroscopy and first-principles studies on the surface states of monoclinic Ag 2 Se thin films grown using molecular beam epitaxy. Through the use of quasi-particle interference patterns, we have observed evidence for topological surface states on the selenium-terminated surfaces with different types of defect densities. The results of ab initio calculations confirm the existence of nontrivial topological surface states in the monoclinic Ag 2 Se structure, for which such properties have not been previously reported. The energy dispersion determined using voltage-dependent standing wave patterns suggests that these topological states have an anisotropic Dirac cone structure. In conclusion, this discovery may lead to new applications for monoclinic Ag 2 Se in the rapidly growing fields of nanoelectronics and spintronics.

36 MATERIALS SCIENCE↗
DOE OSTI2023

The spectroscopic quadrupole moment of the 2 1 + state of 12 C: A benchmark of theoretical models

The spectroscopic quadrupole moment of the first 2 + state of 12 C has been measured employing the Coulomb-excitation re-orientation technique. Our result of $Q_{s} (2^{+}_{1}) = +9.3^{+3.5}_{-3.8}$ efm 2 suggests a larger oblate deformation than previously reported. Combining this with the consistently re-analyzed adopted value, we present the most precise value to date of $Q_{s} (2^{+}_{1})$ = +9.5(18) efm 2 , which is consistent with a geometrical rotor description. This simple outcome is compared to state-of-the-art shell-model, mean-field, ab initio calculations, cluster-based and geometrical-like theories, which show varying degrees of emergent quadrupole collectivity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Modeling inclusive electron-nucleus scattering with Bayesian artificial neural networks

We introduce a Bayesian protocol based on artificial neural networks that is suitable for modeling inclusive electron-nucleus scattering on a variety of nuclear targets with quantified uncertainties. Unlike previous applications in the field, which directly parameterize the cross sections, our approach employs artificial neural networks to represent the longitudinal and transverse response functions. In contrast to cross sections, which depend on the incoming energy, scattering angle, and energy transfer, the response functions are determined solely by the energy and momentum transfer to the system, allowing the angular component to be treated analytically. We assess the accuracy and predictive power of our framework against the extensive data in the quasielastic inclusive electron-scattering database. Additionally, we present novel extractions of the longitudinal and transverse response functions and compare them with previous experimental analysis and nuclear ab-initio calculations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Ab-initio nucleon-nucleon correlations and their impact on high energy 16 O+ 16 O collisions

Investigating nucleon-nucleon correlations inherent to the strong nuclear force is one of the core goals in nuclear physics research. We showcase the unique opportunities offered by collisions of 16 O nuclei at high-energy facilities to reveal detailed many-body properties of the nuclear ground state. We interface existing knowledge about the geometry of 16 O coming from ab-initio calculations of nuclear structure with transport simulations of high-energy 16 O+ 16 O collisions. Bulk observables in these processes, such as the elliptic flow or the fluctuations of the mean transverse momentum, are found to depend significantly on the input nuclear model and to be sensitive to realistic clustering and short-range repulsive correlations, effectively opening a new avenue to probe these features experimentally. This finding demonstrates collisions of oxygen nuclei as a tool to elucidate initial conditions of small collision systems while fostering connections with effective field theories of nuclei rooted in quantum chromodynamics (QCD).

Zhang, Chunjian [Fudan University, Shanghai (China↗
DOE OSTI2021

Probing the antiknock effect of anisole through an ignition, speciation and modeling study of its blends with isooctane

In order to unravel the reaction pathways relevant to anisole co-oxidation within a fuel blend, a detailed study of isooctane/anisole blends was performed with the ULille RCM. Ignition delays as well as mole fraction profiles were measured during a two-stage ignition delay using sampling and GC techniques. These results are used to validate a kinetic model developed from ab initio calculations for the most relevant rate constants which included H-atom abstraction reactions from anisole, and reactions on the potential energy surfaces of methoxyphenyl + O 2 and anisyl + O 2 . Pressure dependent rate constants were computed for the methoxyphenyl + O 2 and anisyl + O 2 reactive systems using master equation code analysis. The new kinetic model shows good agreement with the experimental data. Dual brute-force sensitivity analysis was performed, on both first- and second-stages of ignition, allowing the identification of the most important reactions in the prediction of both ignition delays. It was observed that while pure anisole does not show NTC behavior, a 60/40 isooctane/anisole blend displays such behavior, as well as two-stage ignition. This suggests anisole addition may not be as beneficial to knock resistance as expected from its high octane number. The kinetic modeling results demonstrate the importance of H-abstraction reactions both from the methoxy group and from the aryl ring in ortho-position and the addition of the resultant radicals to O2, mostly leading to the formation of polar or non-aromatic products.

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

Metastable orientation relationships in thin film Cu-Cr bilayers

Metastable orientation relationships (ORs) between Cu and Cr were kinetically stabilized via epitaxial thin film deposition on MgO(001). Both Cu(001) and Cr(001) grow epitaxially on MgO(001). The Bain OR was observed by x-ray diffraction and scanning transmission electron microscopy for Cr(001) / Cu(001) / MgO(001). In contrast, three Cr/Cu ORs were found for Cr deposition on Cu(001) / MgO(001): the Pitsch OR, and two previously unreported ORs related to the Bain and Pitsch ORs, respectively. Ab initio calculations predict the energetics of these ORs, and reveal that the deformation resistance of Cr leads to the three observed ORs on Cu(001).

Kaspar, Tiffany C.↗
DOE OSTI2021

Stability and mobility of tungsten clusters on tungsten (110) surface: Ab initio and atomistic simulations

Quantifying the surface transport properties of tungsten (W) is of prime importance to understand the formation of nano-fuzz in fusion plasma-facing conditions. The stability and mobility of W adatom clusters (W n , n = 2-9) on the W(110) surface has been investigated by computer simulations, including ab initio calculations using density functional theory (DFT) and molecular statics (MS) simulations with multiple W interatomic potentials. The DFT results demonstrate that the sequential binding energy generally increases with number of W adatoms, except for the 5 th and 7 th W adatoms. The most common elemental migration steps of W n (n>2) clusters are observed to consist of monomer and dimer hops, while larger W n clusters can also diffuse by dissociation and recombination of smaller clusters. The threshold migration energy of W 9 is the highest, then followed by W 8 , W 4 , and W 6 , while W 3 , W 5 , and W 7 have similar migration energies. Compared to DFT, each interatomic potential evaluated overestimates the binding energies of W n clusters. Additionally, the embedded-atom potential developed by Juslin and Wirth adequately predicts the threshold migration energy of W n (n>2) clusters, although it predicts different underlying migration mechanisms. The results show that interaction mechanism between W adatoms controls the stability and mobility of W n clusters on the W(110) surface.

36 MATERIALS SCIENCE↗
DOE OSTI2023

The P( 4 S) + NH( 3 Σ – ) and N( 4 S) + PH( 3 Σ – )reactions as sources of interstellar phosphorus nitride

Phosphorus nitride (PN) is believed to be one of the major reservoirs of phosphorus in the interstellar medium (ISM). For this reason, understanding which reactions produce PN in space and predicting their rate coefficients is important for modelling the relative abundances of P-bearing species and clarifying the role of phosphorus in astrochemistry. Here, in this work, we explore the potential energy surfaces of the $\textrm{P}(^4\textrm{S}) + \textrm{NH}(^3\Sigma^-)$ and $\textrm{N}(^4\textrm{S}) + \textrm{PH}(^3\Sigma^-)$ reactions and the formation of $\textrm{H}(^2\textrm{S}) + \textrm{PN}(^1\Sigma^+)$ through high accuracy ab initio calculations and the variable reaction coordinate transition state theory (VRC-TST). We found that both reactions proceed without an activation barrier and with similar rate coefficients that can be described by a modified Arrhenius equation ( $k(T)=\alpha\!\left( T/300 \right)^{\beta} \exp\!{(\!-\!\gamma/T)})$ with $\alpha=0.93\times 10^{-10}\rm cm^3\,s^{-1}$ , $\beta=-0.18$ and $\gamma=0.24\, \rm K$ for the $\textrm{P} + \textrm{NH} \longrightarrow \textrm{H} + \textrm{PN}$ reaction and $\alpha=0.88\times 10^{-10}\rm cm^3\,s^{-1}$ , $\beta=-0.18$ and $\gamma=1.01\, \rm K$ for the $\textrm{N} + \textrm{PH} \longrightarrow \textrm{H} + \textrm{PN}$ one. Both reactions are expected to be relevant for modelling PN abundances even in the cold environments of the ISM. Given the abundance of hydrogen in space, we have also predicted rate coefficients for the destruction of PN via H + PN collisions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Predicting and Understanding Non-Covalent Interactions Using Novel Forms of Symmetry-Adapted Perturbation Theory

Although sometimes derided as “weak” interactions, non-covalent forces play a critical role in ligand binding and crystal packing and in determining the conformational landscape of flexible molecules. Symmetry-adapted perturbation theory (SAPT) provides a framework for accurate ab initio calculation of intermolecular interactions and furnishes a natural decomposition of the interaction energy into physically meaningful components: semiclassical electrostatics (rigorously obtained from monomer charge densities), Pauli or steric repulsion, induction (including both polarization and charge transfer), and dispersion. This decomposition helps to foster deeper understanding of non-covalent interactions and can be used to construct transferable, physics-based force fields. Separability of the SAPT interaction energy also provides the flexibility to construct composite methods, a feature that we exploit to improve the description of dispersion interactions. Furthermore, these are challenging to describe accurately because they arise from nonlocal electron correlation effects that appear for the first time at second order in perturbation theory but are not quantitatively described at that level.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Oxygen Affinity: The Missing Link Enabling Prediction of Proton Conductivities in Doped Barium Zirconates

Proton-conducting oxides, specifically doped barium zirconates, have garnered much attention as electrolytes for solid-state electrochemical devices operable at intermediate temperatures (400–600 °C). In chemical terms, hydration energy, E hyd , and proton–dopant association energy, E as , are two key parameters that determine whether an oxide exhibits fast proton conduction, but to date ab initio studies have for the most part studied each parameter separately, with no clear correlation with proton conductivity identified in either case. Here, we demonstrate that the oxygen affinity, E O.dopant , defined as the energy released when an oxide ion enters an oxygen vacancy close to a dopant atom, is the missing link between these two parameters and correlates well with experimental proton conductivities in doped barium zirconates. Ab initio calculations of point defects and their complexes in Sc-, In-, Lu-, Er-, Y-, Gd-, and Eu-doped barium zirconates are used to determine E hyd , E as , E O.dopant , and the hydrogen affinity, E H.host , of each system. These four energy terms are related by E hyd = E O.dopant + 2E H.host + 2E as . Complementary impedance spectroscopy measurements reveal that the stronger the calculated oxygen affinity of a system, the higher the proton conductivity at 350 °C. Although the proton trapping site is also an important factor, the results show that oxygen affinity is an excellent predictor of proton conductivity in these materials.

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

Competitive Adsorption of NH 3 and H 2 O in Metal–Organic Framework Materials: MOF-74

Elucidating the interaction between coadsorbed H 2 O and NH 3 in metal–organic frameworks (MOFs) is of paramount importance to uncover mechanistic details of their competitive coadsorption behavior as well as to guide the design of new materials for enhanced NH 3 adsorption in humid environments. Nevertheless, molecular competition between NH 3 and H 2 O within the confined nanopores of MOFs was rarely explored and is poorly understood due to challenges in characterization. Here, we combine in situ infrared spectroscopy with ab initio calculations to unveil the competition of NH 3 and H 2 O for occupying active adsorption sites in the representative MOF-74 material by analyzing the kinetics and energetics of the molecular exchange process. We find that at a high NH 3 /H 2 O ratio, the incoming NH 3 is capable of displacing metal-bound H 2 O and moving it to secondary adsorption sites due to the stronger binding of NH 3 compared with H 2 O. Interestingly, the reverse process of H 2 O displacing metal-bound NH 3 is also possible upon increasing water concentration. Our calculations show that H 2 O exchanging the preabsorbed NH 3 at the metal site is driven not only by a reduced kinetic barrier but also by a favorable energetical state resulting from the formation of water clusters at metal sites and intermolecular H-bonding between the metal-coordinated H 2 O and displaced NH 3 . Our finding emphasizes that the description of molecular occupation in MOFs at equilibrium cannot simply be established by comparing molecules’ binding energies at their strongest binding sites derived by single-component measurements; rather, intermolecular interactions can greatly affect molecular distribution at equilibrium. Furthermore, we show that vibrational modes of adsorbed NH 3 are markedly perturbed upon contact with water molecules, accompanied by a large frequency shift (>30 cm –1 ) and considerable intensity decrease, which arises from the freezing of NH 3 vibrations by coadsorbed H 2 O. As a result, the mechanistic insight obtained through our study sheds light on molecular coadsorption processes in MOFs and helps to assess NH 3 removal efficiency of MOFs containing open-metal sites under realistic conditions, particularly in the presence of humidity.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Synthesis, Stability, and Magnetic Properties of Antiperovskite Co 3 PdN

Experimental synthesis and characterization of theoretically predicted compounds are important steps in the materials discovery pipeline. Here, we report on the synthesis of Co 3 PdN, which was recently predicted to be a stable magnetic antiperovskite. The Co 3 PdN thin films were grown by reactive sputtering and were confirmed to form in an antiperovskite crystal structure. The thermal stability of the compound is demonstrated up to 600 K by in situ X-ray diffraction, though the phase persists at slightly higher temperatures (700 K) in an air-free magnetometer. Both ab initio calculations and magnetization measurements find Co 3 PdN to be ferromagnetic with an experimentally determined Curie temperature of T C = 560 ± 5 K. The saturation magnetization of 1.2 μ B /Co found in the experiment is slightly lower than the 1.7 μ B /Co value expected by theory. A narrow magnetic hysteresis loop with a coercive field of 100 Oe at low temperature suggests that Co 3 PdN might be useful in electronic applications requiring fast switching of the magnetization vector. While prior prediction of Co 3 PdN showed a gapped electronic band structure for each spin channel, we show that this was due to incomplete sampling of Brillouin zone paths and that band crossings exist along R-X|M and X|M-R paths. The metallic nature of Co 3 PdN is further confirmed by temperature-dependent transport measurements, which also show a considerable anomalous Hall effect. Altogether, this work represents an appreciable step toward understanding the synthesis, structure, stability, and properties of a new magnetic material.

14 SOLAR ENERGY↗
DOE OSTI2020

Giant Magnetoelectric Coupling and Magnetic-Field-Induced Permanent Switching in a Spin Crossover Mn(III) Complex

We investigate giant magnetoelectric coupling at a Mn 3+ spin crossover in [Mn III L]BPh 4 (L = (3,5-diBr-sal) 2 323) with a field-induced permanent switching of the structural, electric, and magnetic properties. An applied magnetic field induces a first-order phase transition from a high spin/low spin (HS-LS) ordered phase to a HS-only phase at 87.5 K that remains after the field is removed. We observe this unusual effect for DC magnetic fields as low as 8.7 T. The spin-state switching driven by the magnetic field in the bistable molecular material is accompanied by a change in electric polarization amplitude and direction due to a symmetry-breaking phase transition between polar space groups. The magnetoelectric coupling occurs due to a γη 2 coupling between the order parameter γ related to the spin-state bistability and the symmetry-breaking order parameter η responsible for the change of symmetry between polar structural phases. We also observe conductivity occurring during the spin crossover and evaluate the possibility that it results from conducting phase boundaries. We perform ab initio calculations to understand the origin of the electric polarization change as well as the conductivity during the spin crossover. Thus, in this work, we demonstrate a giant magnetoelectric effect with a field-induced electric polarization change that is 1/10 of the record for any material.

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

Revisiting Competitive Adsorption of Small Molecules in the Metal–Organic Framework Ni-MOF-74

To precisely evaluate the potential of metal-organic frameworks (MOFs) for gas separation and purification applications, it is crucial to understand how various molecules competitively adsorb inside MOFs. In this paper, we combine in situ infrared spectroscopy with ab initio calculations to investigate the mechanisms associated with co-adsorption of a number of small molecules including CO, NO, and CO 2 inside the prototypical framework Ni-MOF-74. Surprisingly, we find that the displacement of CO bound inside Ni-MOF-74 (binding energy of 53 kJ/mol) is readily driven by CO 2 exposure, even though CO 2 has a noticeably weaker binding energy of only 41 kJ/mol; meanwhile, the significantly more strongly binding NO molecule (90 kJ/mol) is not able to easily displace bound CO inside Ni-MOF74. These results show that single-phase binding energies of a molecule inside the MOF cannot completely describe their interaction with the MOF in the presence of other guest molecules. Here, we unveil a number of crucial factors such as the kinetic barrier, partial pressure, secondary binding sites, and guest-host/lateral interactions that control the co-adsorption process and combined with the binding energy are better descriptors of the behavior and adsorption of gas mixtures inside MOFs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗