SEARCH · Engineering Papers

Results for “ab-initio”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

DOE OSTI2024

A Theoretical Analysis on the Dynamic Behavior and Electromechanical Response of PZT

Piezoelectric materials are used as a power source in stockpile components for safety and reliability assurances. The objective of this project is to gain insights into the fundamental pressure-induced behavior and electromechanical response of lead zirconate titanate (PZT). Specifically, to establish the basis for an accurate, physics-based material model for robust model-based design to rapidly optimize PZT-based materials and components for performance studies. Established ab-initio methods are used to interrogate and understand the dynamic behavior of PZT as a function of composition (50/50, 65/35, 80/20, 95/5) and dopants (La, Nb) to overcome the costly and time-consuming experimental methodologies. New cold curves for pure and doped single-crystal PZT are obtained as the reference equation-of-state (EOS). A negligible change in the pressure responses was observed for the systems and strains studied. Dielectric and piezoelectric responses of pure and doped single-crystal PZT were also calculated as a function of pressure. For undoped PZT, there is a clear and orderly pressure and composition dependence. For doped PZT, there is a significant increase in the responses, but the behavior is very disordered and inconclusive.

36 MATERIALS SCIENCE↗

DOE OSTI2024

FY24 Progress Report on Viscosity and Thermal Conductivity Measurements of Nuclear Industry Relevant Chloride Salts: An Experimental and Computational Study

As presented in this report, experimental and computational techniques were performed to assess the viscosity and thermal conductivity of key alkali and actinide chloride mixtures for molten salt reactor developers. These mixtures were pure LiCl, NaCl-KCl, LiCl-NaCl, LiCl-KCl, LiCl-NaCl-KCl, and NaCl-UCl 3 . Experimental measurements of viscosity were performed with a rolling ball viscometer, whereas experimental measurements of thermal conductivity were performed with a variable gap apparatus. Additional benchmarking work was performed using both property measurement systems to prepare for x-ray radiography in stainless-steel crucibles for viscosity and to ensure that calibration methods were accurate for thermal conductivity before assessing the NaCl-UCl 3 system. Validation data for the NaCl-UCl 3 in literature are minimal. Details on the calibration methods, salt measurement processes, and sources of error and uncertainty are discussed in detail for both property measurements. The computational methods described herein involved ab-initio molecular dynamics (AIMD) calculations using CP2K. The calculations were performed for the LiCl-KCl-NaCl and NaCl-UCl 3 systems. These calculations not only provided thermophysical property estimations for comparison to experimental data, but they also allowed for the determination of diffusion coefficients, coordination numbers, and radial distribution functions to provide insight into ion mobility and local coordination environments, which is linked to macroscopic property trends.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2024

Atomistic Simulations for Thermophysical Properties of Uranium-Containing Halide Molten Salts

Characterizing the thermophysical properties in both fuel and coolant salts are critical in modeling, developing, process optimizing and utilizing molten salt reactors (MSRs), as these properties directly relate to operation metrics and can inform on the selection of candidate salts. The demand for consistent, accurate and publicly available thermophysical property data has become more apparent in recent years as interests have increased from molten salt reactor developers. There are a number of challenges in experimentally measuring properties such as thermal conductivity, viscosity, density and heat capacity , which have led to sparse and often times conflicting data points or molten salts in general. Additionally, there are a number of hazards to consider when synthesizing, storing, using, treating and disposing of molten salts. With the advances in computational capabilities over the last 10 years, the use of atomistic simulations can be implemented to support these efforts. The primary objective of this work is characterize the thermophysical transport properties in a number of molten chloride salts, and in particular NaCl-UCl 3 using ab-initio molecular dynamic (AIMD) simulations. In this binary salt the UCl 3 acts as the primary fissile material and NaCl acts as a carrier salt due with its’ high solubility for actinides A number of studies on the thermophysical properties of NaCl-UCl 3 have been published but there is not a vast amount of viscosity data for this system. In 1975, Desyatnik, et al published a study reporting dynamic viscosities that were calculated from kinematic viscosity measurements, and using the coefficients provided the viscosity in a 70:30 NaCl:UCl 3 mixture is 2.29 cP and 2.88 for a 60:40 mixture. Termini et al. recently reported viscosities in the range of 2.75 – 3 cP for the 63:37 NaCl-UCl 3 mixture in the same temperature range using rolling ball viscosity measurements. Computational viscosity of a similar mixture (64:36) can be obtained from the work Andersson et al. using the reported diffusion coefficients, and the hydrodynamic radius from the pair-radial distribution functions (RDFs). Using Eq (1) (vida infra), the viscosity would be 2.50 cP at 1100K. This is not to say that these values are incorrect due to the varying reported values, but aims to highlight the necessity of this work. The data reported in this ongoing work are computations on a 64:36 mixture of NaCl-UCl 3 at 987K. This work is likely to be expanded into varying concentrations of this mixture along with the inclusion of other salt candidate mixtures.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

DOE OSTI2025

Development of a High-Rate Lithium-Air Battery Using a Gaseous CO 2 Reactant

Li-air batteries are considered a potential alternative to Li-ion batteries for transportation applications due to their high theoretical specific energy. Most works in this area focus on use of O 2 as the reactant. However, newer concepts for using gaseous reactants (such as CO 2 , which has a theoretical specific energy density of 1,876 Wh/kg) provide opportunities for further exploration. The main objective of this project was the development of a novel strategy that enables operation of Li-CO 2 batteries at high-capacity and high-rate, with a long-cycle-life. The team was able to: (1) Synthesize two novel transition metal chalcogenide (TMC) catalysts that work in synergy with ionic liquid-based electrolytes to enhance the efficiency of reactions during discharge and charge processes; (2) Fabricate high-porosity cathode electrodes with 3D printing to increase electrode surface area and gas permeability; (3) Develop a multiscale modeling framework that integrates Density Functional Theory (DFT), Ab-Initio Molecular Dynamics (AIMD), classical Molecular Dynamics (MD), and Finite Element Analysis (FEA) to investigate atomic and cell-level properties of Li-CO 2 batteries; (4) Assemble a stackable Li-CO 2 pouch-cell able to deliver a capacity of >200 mAh. These achievements were realized through an integrated approach based on materials synthesis, testing, characterization, analysis, and computation. This project produced a thorough understanding of key chemical, electronic, and kinetic parameters that govern the operation of Li- CO 2 batteries in realistic conditions. The methodologies employed, and the insight generated, will be valuable beyond advancing the field of Li-CO 2 batteries

25 ENERGY STORAGE↗

DOE OSTI2026

Planar Systems for Quantum Information

This project aims to develop two‐dimensional (2D) moiré materials as a quantum simulator to implement model Hamiltonians and their phase diagrams. Progress in quantum information science (QIS) requires the development of advanced quantum materials systems. The rich family of layered van der Waals materials and their heterostructures present opportunities to create previously unrealized types of applications for QIS. Specifically, when two layers of van der Waals materials are overlaid with a small twist angle or/and lattice mismatch, a moiré superlattice with a period of about ten nanometers is formed. This provides a periodic trapping potential for electrons. Electrons can tunnel between the traps and repel each other by their mutual Coulomb interactions. The platform of 2D moiré materials provides many attractive features, including tunability of length and energy scales, charge density, and even lattice symmetry. It presents new possibilities for realizing quantum simulation of the many-body physics in a solid-state platform. This integrated team of six investigators seeks to develop relevant theoretical treatments to link ab-initio studies of 2D moiré materials to model Hamiltonians and to evaluate correlated phases predicted by these model Hamiltonians in the relevant regimes. On the experimental side, the team aims to develop methods to realize a homogeneous and highly controlled potential landscape for the electrons and to initiate, protect, and measure their quantum many-body states.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Studying electroweak few-body observables in chiral effective field theory

The use of nuclei to study electroweak probes is becoming increasingly relevant experimentally. The success of dark matter and neutrino experiments strongly depends on the ability to control nuclear effects in order to extract the fundamental parameters associated with external probes. Therefore, reliable theoretical calculations of nuclear structure and reactions, with well-controlled errors, are crucial for the success of experimental efforts. Currently, chiral effective field theory ($\chi$EFT) coupled with {\it ab-initio} methods represents one of the best approaches that fulfills these requirements. To use this approach as a tool for studying fundamental physics, it is essential to validate it against experimental data for which the calculations are well under control, such as the elastic scattering of electrons on nuclei. In this proceeding, I will present recent developments in the fitting of electromagnetic currents derived using $\chi$EFT and the calculation of electromagnetic form factors of light nuclei. The results of these calculations demonstrate the strength of the theory in describing the interaction of nuclei with electromagnetic probes over a broad range of momentum transfers and highlight the robustness of $\chi$EFT for analyzing future experimental data aimed at extracting fundamental parameters.

Gnech, Alex [Old Dominion Univ., Norfolk, VA (Unit↗

DOE OSTI2026

Computing neutrino cross sections from Euclidian responses

Energy integrated neutrino cross sections are integrals of nuclear responses weighted with kinematic prefactors. We decompose the prefactors into a limited set of functions of energy transfer and show the relevant integrals are the moments of the responses, and integrals weighted with $1/(a+ω)^n$ with $n\leq 2$. These can be directly obtained from the Euclidean response, avoiding the need for inversion of the Laplace transform. As a proof of concept we study the procedure with toy-model responses for the quasielastic peak. We show that the different contributions can be straightforwardly organized in terms of relative importance, and how flux-averaged cross sections can be obtained. Using a realistic model for the response and numerical uncertainty we show that it is feasible to obtain the required integrals from the Euclidean response, with large uncertainties only for the third moment. Due to kinematic restrictions, the integrals contain contributions from the unphysical region for neutrino scattering, coming from high-momentum nucleons. We show that (in the absence of two-body currents) robust corrections for this contamination are obtained from the single-nucleon momentum distribution. These results present an opportunity to compute certain neutrino cross sections with ab-initio methods with controlled uncertainties.

Nikolakopoulos, A. [Washington U., Seattle]↗

DOE OSTI2026

Computing neutrino cross sections from Euclidian responses

Nikolakopoulos, A. [Washington U., Seattle]↗

DOE OSTI2025

Static Subspace Approximation for Random Phase Approximation Correlation Energies: Applications to Materials for Catalysis and Electrochemistry

Modeling complex materials using high-fidelity, ab initio methods at low cost is a fundamental goal for quantum chemical software packages. The GW approximation and random phase approximation (RPA) provide a unified description of both electronic structure and total energies using the same physics in a many-body perturbative approach that can be more accurate than generalized-gradient density functional theory (DFT) methods. However, GW/RPA implementations have historically been limited to either specific materials classes or application toward small chemical systems. Here, the static subspace approximation allows for reduced cost full-frequency GW/RPA calculations and has previously been benchmarked thoroughly for GW calculations. Here, we describe our approach to including partial occupations of electronic orbitals in full-frequency GW and RPA calculations for the study of electrocatalysts. We benchmarked RPA total energy calculations using the subspace approximation across a diverse test suite of materials for a variety of computational parameters. The benchmarking quantifies the impact of different extrapolation procedures for representing the static polarizability at infinite screened cutoff, and shows that using screened cutoffs above 20-25 Ryd result in diminishing accuracy returns for predicting RPA total energies. Additionally, for moderately sized electrocatalytic models, 2-3 times fewer computational resources are used to compute RPA total energies by representing the static polarizability with 20-30% of the static subspace basis, with an error of approximately 0.01 eV or better in RPA adsorption energy calculations. Finally, we show that for these electrochemical models RPA can shift DFT adsorption energy shifts by up to 0.5 eV and that GW can frequently shift DFT eigenvalues of surface and adsorbate states by approximately 0.5-1 eV.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Observation of halogen-like behavior of gold in fluorinated bimetallic CoAuF 1-2 − and CuAuF 1-2 − clusters: Anion photoelectron spectroscopy and density functional theory

Using size-selected anion photoelectron spectroscopy and density functional theory, we investigated the structures and properties of fluorinated bimetallic clusters CoAuF 1-2 − and CuAuF 1-2 − and their neutrals. Both experimental and theoretical results show that in these cluster anions, Au behaves like a halogen atom. For example, the measured vertical detachment energies (VDEs) of CoAuF − (2.00 ± 0.08 eV) and CuAuF − (3.8 ± 0.1 eV) are close to those of CoF 2 − (2.12 ± 0.08 eV) and CuF 2 − (3.58 ± 0.08 eV), respectively. The theoretical results show that the geometries and electronic structures of CoAuF − and CuAuF − are similar to those of CoF 2 − and CuF 2 − . The natural population analysis and natural electron configuration analyses further confirm that the electronic properties of Au in MAuF − (M = Co, Cu) mimic those of MF 2 − . In addition, the electron localization function analyses show that the M-Au chemical bonds are similar to the corresponding M-F chemical bonds, providing evidence for the ionic nature of the interactions. When a second F atom is attached to the CoAuF − and CuAuF − clusters, the VDEs of the resulting CoAuF 2 − and CuAuF 2 − are 4.38 ± 0.08 eV and 3.71 ± 0.08 eV, respectively, indicating their superhalogen character as these values are higher than those of halogen anions. The results may be useful for understanding the properties of gold at the nanoscale that play an important role in catalysis and nanotechnology.

Ab-initio methods↗

DOE OSTI2024

Anion photoelectron spectroscopy and chemical bonding of ThS 2 – and ThSO –

Anion photoelectron spectra of ThSO – and ThS 2 – were recorded using the third (355 nm) harmonic of an Nd-YAG laser; these provided the measured vertical detachment energies of each anion. The experiments are supported by extensive coupled cluster calculations on ThSO, ThSO – , ThS 2 , and ThS 2 – , as well as the oxygen congeners ThO 2 and ThO 2 – . The ab initio calculations, which included complete basis set extrapolations, spin–orbit effects using four-component coupled cluster, and higher-order correlation contributions through CCSDT(Q), yielded an adiabatic electron affinity for ThO 2 – that was within 0.02 eV of the previously determined experimental value. The singly occupied molecular orbital (SOMO) in all three anions corresponds primarily to the 7s orbital on Th. Successive substitution of S for each O in ThO 2 leads to larger electron affinities and smaller bond angles in the neutral molecules, but larger angles in the anions. As demonstrated by Franck–Condon simulations of the spectra using the CCSD(T) spectroscopic constants, substitution of O by S significantly complicates the resulting detachment spectra due to the lower vibrational frequencies in the sulfur species. Altogether the calculated vertical detachment energies are in very good agreement with the experiment. In addition to the adiabatic electron affinities of each species, atomization energies and heats of formation have also been determined via the FPD approach with expected uncertainties of 1–2 kcal/mol.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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

Model-free Rayleigh weight from x-ray Thomson scattering measurements

X-ray Thomson scattering (XRTS) has emerged as a powerful tool for the diagnostics of matter under extreme conditions. In principle, it gives one access to important system parameters such as the temperature, density, and ionization state, but the interpretation of the measured XRTS intensity usually relies on theoretical models and approximations. In this context, a key property is given by the Rayleigh weight that describes the electronic localization around the ions. Here, we show that it is possible to extract the Rayleigh weight directly from the experimental data without the need for any model calculations or simulations. As a practical application, we consider an experimental measurement of strongly compressed Be at the National Ignition Facility [Döppner et al., Nature 618, 270–275 (2023)]. We demonstrate that experimental results for the Rayleigh weight open up new avenues for the interpretation of XRTS experiments by matching the measurement with ab initio simulations such as density functional theory or path integral Monte Carlo. Interestingly, this new procedure leads to significantly lower density compared to previously used chemical models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Excited electronic states of Na 2 and K 2 : The potential for long-lived “reservoir” states leading to collision induced population inversions

Potential energy curves (PECs) for the spin-free (ΛS) and spin–orbit (Ω) states associated with the four lowest-lying dissociation channels of Na 2 and K 2 were calculated at the SA-CASSCF/SO-CASPT2/aug-cc-pwCVQZ-DK level. The PECs of Na 2 were consistent with the experimental data and with the FS-CCSD (2,0) calculations, reproducing the double-well and the “shelf” character for some of the potentials of the excited states. For K 2 , the PECs behaved in a similar way and the spectroscopic parameters for the ground and the excited states are in good agreement with the available experimental values. The dissociation energy of K 2 was predicted to be D e = 4454 cm −1 , within an agreement of 5 cm −1 with the experiments. For Na 2 , D e = 5789 cm −1 compared to the experimental value of 6022 cm −1 . The inclusion of spin–orbit coupling effects resulted in avoided crossings, which affect the PECs. Spin–orbit changes the predicted curves for some excited Ω states arising from ΛS states that overlap each other, affecting their associated vibrational frequencies and bond distances. Here, the current studies of the low-lying states in K 2 reveal a similar structure to those of Na 2 , which suggests the accessibility of long-lived energy storing reservoir states and possible population inversions in K 2 following prior experimental work on the reaction of halogen atoms with Na 3 to produce excited states of Na 2 .

Ab-initio methods↗

DOE OSTI2025

Influence of surface chemistry on Li nucleation energetics on graphene-based surfaces

Lithium metal is a promising high-capacity anode material for solid-state batteries, but it typically suffers from poor cyclability. Carbon scaffold hosts have the potential to improve this performance due to their high electronic conductivity and large surface area, which facilitates lithium-ion adsorption and desorption. Scaffold surface chemistry is known to significantly influence performance outcomes, but the details of these interactions are not fully understood. Here, this study employs first-principles simulations to explore lithium transport and nucleation on graphene anodes with various surface chemistries. Using enhanced sampling techniques, ab initio molecular dynamics, and density functional theory calculations, we find that although surface chemistry has a minimal impact on lithium interfacial transport, it influences surface nucleation significantly. Both heteroatom dopants and intrinsic defects lower the nucleation barrier, creating a more favorable environment for lithium nucleation compared to pristine graphene. In addition, our results reveal a complex interplay between surface lithium concentration, lithium transport, and nucleation kinetics. These findings highlight the potential of surface modifications to precisely control nucleation processes on carbon-based anodes and provide design guidance for reducing dendrite formation and improving the cycle life of solid-state batteries.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Accurate and efficient parameterization of an atomic cluster expansion (ACE) potential for ammonia under extreme conditions

We present a machine learning interatomic potential for ammonia designed to capture its complex multiphase behavior, including both molecular and superionic phases. The potential is based on the atomic cluster expansion (ACE) formulation and has been parameterized to facilitate high-fidelity molecular dynamics simulations of ammonia under extreme conditions, for pressures up to 100 GPa and for temperatures above 500 K and up to 6000 K. A diverse range of configurations was generated through high-quality ab initio molecular dynamics simulations, covering insulating and superionic ice phases, liquid ammonia, molecular nitrogen (N 2 ) and hydrogen (H 2 ), and metastable compounds that form upon dissociation, including $NH^{+}_{4}$, $H^{+}_{3}$, N 2 H 4 , and N 3 H. We demonstrate that the ammonia ACE potential accurately reproduces experimental and density functional theory predicted isotherms and Hugoniots. Crucially, the potential is able to capture the intricate phase behavior of ammonia, including the transition from insulating molecular fluid to the superionic phase. This work provides a robust interatomic potential that can be used for large-scale, accurate simulations of ammonia under extreme thermodynamic conditions, offering a powerful tool for investigating its behavior in various phases and applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

How the choice of exchange–correlation functional affects DFT-based simulations of the hydrated electron

Hydrated electrons are anionic species that are formed when an excess electron is introduced into liquid water. Building an understanding of how hydrated electrons behave in solution has been a long-standing effort of simulation methods, of which density functional theory (DFT) has come to the fore in recent years. The ability of DFT to model the reactive chemistry of hydrated electrons is an attractive advantage over semi-classical methodologies; however, relatively few density functional approximations (DFAs) have been used for the hydrated electron simulations presented in the literature. Here, we simulate hydrated electron systems using a series of exchange–correlation (XC) functionals spanning Jacob’s ladder. We calculate a variety of experimental and other observables of the hydrated electron and compare the XC functional dependence for each quantity. We find that the formation of a stable localized hydrated electron is not necessarily limited to hybrid XC functionals and that some hybrid functionals produce delocalized hydrated electrons or electrons that react with the surrounding water at an unphysically fast rate. Here, we further characterize how different DFAs impact the solvent structure and predicted spectroscopy of the hydrated electron, considering several methods for calculating the hydrated electron’s absorption spectrum for the best comparison between structures generated using different density functionals. None of the dozen or so DFAs that we investigated are able to correctly predict the hydrated electron’s spectroscopy, vertical detachment energy, or molar solvation volume.

Ab-initio molecular dynamics↗

DOE OSTI2025

GPU acceleration of hybrid functional calculations in the SPARC electronic structure code

We present a Graphics Processing Unit (GPU)-accelerated version of the real-space SPARC electronic structure code for performing hybrid functional calculations in generalized Kohn–Sham density functional theory. In particular, we develop a batch variant of the recently formulated Kronecker product-based linear solver for the simultaneous solution of multiple linear systems. We then develop a modular, math kernel based implementation for hybrid functionals on NVIDIA architectures, where computationally intensive operations are offloaded to the GPUs, while the remaining workload is handled by the central processing units (CPUs). Considering bulk and slab examples, we demonstrate that GPUs enable up to 8× speedup in node-hours and 80× in core-hours compared to CPU-only execution, reducing the time to solution on V100 GPUs to around 300 s for a metallic system with over 6000 electrons, and significantly reducing the computational resources required for a given wall time.

Kohn-Sham density functional theory↗