Search NASA⌕ Search

SEARCH · Search NASA

Results for “first principles simulations”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

First-principles equation of state database for warm dense matter computation

We put together a first-principles equation of state (FPEOS) database for matter at extreme conditions by combining results from path integral Monte Carlo and density functional molecular dynamics simulations of the elements H, He, B, C, N, O, Ne, Na, Mg, Al, and Si as well as the compounds LiF , B 4 C , BN , CH 4 , CH 2 , C 2 H 3 , CH , C 2 H , MgO , and MgSiO 3 . For all these materials, we provide the pressure and internal energy over a density-temperature range from ~ 0.5 to 50 g cm - 3 and from ~ 10 4 to 10 9 K, which are based on ~ 5000 different first-principles simulations. We compute isobars, adiabats, and shock Hugoniot curves in the regime of L - and K -shell ionization. Invoking the linear mixing approximation, we study the properties of mixtures at high density and temperature. Furthermore, we derive the Hugoniot curves for water and alumina as well as for carbon-oxygen, helium-neon, and CH-silicon mixtures. We predict the maximal shock compression ratios of H 2 O , H 2 O 2 , Al 2 O 3 , CO , and CO 2 to be 4.61, 4.64, 4.64, 4.89, and 4.83, respectively. Finally we use the FPEOS database to determine the points of maximum shock compression for all available binary mixtures. We identify mixtures that reach higher shock compression ratios than their end members. We discuss trends common to all mixtures in pressure-temperature and particle-shock velocity spaces. In the Supplemental Material, we provide all FPEOS tables as well as computer codes for interpolation, Hugoniot calculations, and plots of various thermodynamic functions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Synergic Effects of Surface Chemistry and Applied Potentials on the Kinetics of the Electrocatalytic Oxygen Evolution Reaction in IrO 2

Improving the efficacy of the oxygen evolution reaction (OER) through water oxidation is critical for advancing photoelectrochemical water splitting. Among the catalysts, IrO 2 exhibits a high OER catalytic activity and stability under acidic conditions. The OER mechanism in this system has been a topic of intense research; however, many mechanistic understandings are lacking. In this work, we carried out first-principles calculations with an implicit solvation model at a constant potential to study the OER behavior on the IrO 2 (110) surface. We find that the surface hydrogen coverage has a significant effect on OER kinetics and transition states. We then develop a microkinetics model that accounts for the continuous evolution of both OER activation energy and hydrogen coverage as a function of an applied potential. We show that this inclusion leads to significant improvement in the simulated Tafel plot compared to available experiments. In conclusion, our results point to a complex interplay between surface chemistry and the applied potential on OER kinetics.

08 HYDROGEN↗
DOE OSTI2021

Probing the electronic properties of the electrified silicon/water interface by combining simulations and experiments.

Silicon (Si) is broadly used in electrochemical and photoelectrochemical devices, where the capacitive and Faradaic reactions at the Si/water interfaces are critical for signal transduction or noise generation. However, probing the electrified Si/water interface at the microscopic level remains a challenging task. Here we focus on hydrogenated Si surfaces in contact with water, relevant to transient electronics and photoelectrochemical modulation of biological cells and tissues. We show that by carrying out first principles molecular dynamics simulations of the Si(100)/water interface in the presence of an electric field we can realistically correlate the computed flat-band potential and tunneling current images at the interface with experimentally measured capacitive and Faradaic currents. Specifically, we validate our simulations in the presence of bias by performing pulsed chronoamperometry measurements on Si wafers in solution. Consistent with prior experiments, our measurements and simulations indicate the presence of voltage-dependent capacitive currents at the interface. We also find that Faradaic currents are weakly dependent on the applied bias, which we relate to surface defects present in newly prepared samples.

Ye, Zifan↗
DOE OSTI2026

Machine Learning Interatomic Potentials for Modeling Framework Flexibility and Water Uptake in NbOFFIVE-1-Ni Metal–Organic Framework

Metal–organic frameworks (MOFs), with their distinctive porous structures and tunable chemical properties, have shown immense promise in the separation and storage of gases. Currently, the accurate simulation of their adsorptive properties remains challenging, especially for systems where the molecules fit very tightly into the pores. Traditional simulation methods often approximate the frameworks as rigid and do not account for the framework flexibility seen in materials such as NbOFFIVE-1-Ni. First-principles molecular dynamics (FPMD) simulations offer the desired accuracy in modeling this flexibility but are limited by their extensive computational demands, rendering them impractical for long simulations. Conversely, classical force field-based simulations offer computational efficiency but lack the necessary accuracy. Here, to break this accuracy-efficiency trade-off, we have developed machine learning interatomic potentials trained on energies and forces from FPMD to model the framework flexibility of NbOFFIVE-1-Ni in the presence of water over nanosecond time scales. Furthermore, by integrating MLIP-driven molecular dynamics (MLIP-MD) with grand canonical Monte Carlo (GCMC) simulations, we further incorporated framework flexibility into adsorption predictions, yielding water adsorption isotherms that better align with experimental data compared to those of conventional GCMC simulations. These advances offer new opportunities for the design and optimization of MOFs in gas storage and separation applications.

adsorption↗
DOE OSTI2025

Artificial Intelligence/Deep Learning FRNN Software for Prediction & Real-Time Control of DIII-D Plasma Control System (PCS)

This collaborative project integrated an improved version of the Artificial Intelligence/Deep Learning FRNN prediction and control software into the real-time DIII-D PCS (plasma control system). A key AI/DL software challenge is to build a modern high-performance computing (HPC) enabled “synthetic plasma simulator” capable of carrying out HPC-driven real-time plasma control applications. This involves development of a deep learning framework to train the surrogate model for a first-principles-based instability analysis simulator (“SGTC”) derived from the global gyrokinetic code GTC. The role of SGTC is to provide accurate and detailed plasma instability information from a real-time AI-based simulator capability to complement the deep learning prediction and control from experimentally-measured signals, such as ECE Imaging, supplemented by synthetic SGTC-ECEI.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Thermodynamics of Liquid Uranium from Atomistic and Ab Initio Modeling

We present thermodynamic properties for liquid uranium obtained from classical molecular dynamics (MD) simulations and the first-principles theory. The coexisting phases method incorporated within MD modeling defines the melting temperature of uranium in good agreement with the experiment. The calculated melting enthalpy is in agreement with the experimental range. Classical MD simulations show that ionic contribution to the total specific heat of uranium does not depend on temperature. The density of states at the Fermi level, which is a crucial parameter in the determination of the electronic contribution to the total specific heat of liquid uranium, is calculated by ab initio all electron density functional theory (DFT) formalism applied to the atomic configurations generated by classical MD. The calculated specific heat of liquid uranium is compared with the previously calculated specific heat of solid γ-uranium at high temperatures. The liquid uranium cannot be supercooled below T sc ≈ 800 K or approximately about 645 K below the calculated melting point, although, the self-diffusion coefficient approaches zero at T D ≈ 700 K. Uranium metal can be supercooled about 1.5 times more than it can be overheated. The features of the temperature hysteresis are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Changes in nitrate binding with lanthanides in BLPhen complexes

Preorganized ligands such as bis-lactam-1,10-phenanthroline (BLPhen) offer unique selectivity trends in solvent extraction of rare-earth ions, particularly from aqueous nitrate solutions. However, due to the ligand’s lipophilicity it is experimentally challenging to obtain crystal structures of the corresponding lanthanide (Ln):ligand complexes to describe the local structure of the first coordination shell and the role of the nitrate ions in the organic phase. In this work, we report first principles molecular dynamics (FPMD) simulations of the complexation of trivalent lanthanide ions (La, Nd, and Eu) with the BLPhen ligand in the presence of nitrate ions in the dichloroethane (DCE) solvent. We find that two nitrate anions are present in the first solvation shell and the third nitrate anion is far from the first solvation shell in all three [Ln(BLPhen) 2 ] 3+ complexes examined. Moving along the lanthanide series from La to Nd to Eu, the binding pocket formed by the two BLPhen ligands shrinks in size; as a result, the two nitrate configuration changes from one-monodentate/one-bidentate binding mode for La and Nd to all monodentate for Eu. Meanwhile, the total coordination number drops from 11 for La and Nd to 10 for Eu. More interestingly, the bidentate binding in [La(BLPhen) 2 ] 3+ and [Nd(BLPhen) 2 ] 3+ complexes is highly dynamic, frequently switching to monodentate and back. The FPMD insights into the differing modes and dynamics of nitrates in the first coordination shell will be useful for further atomistic understanding of the complex structure, formation, and stability in the organic phase.

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

Pressure-induced transformation of Nb 2 O 5 under shock compression from first principles

Ab initio molecular dynamics (AIMD) simulations were carried out to investigate the equation of state of Nb 2 O 5 and its pressure-density relationship under shock conditions. Here, the focus of this study is on the monoclinic B-Nb 2 O 5 (C2/c) polymorph. Enthalpy calculations from AIMD trajectories at 300 K show that the pressure-induced transformation between the thermody-namically most stable crystalline monoclinic parent phase H-Nb 2 O 5 (P2/m) and B-Nb 2 O 5 occurs at ~1.9 GPa. This H→B transition is energetically more favorable than the H→L(Pmm2) pressure-induced transition recently observed at ~5.9-9.0 GPa. The predicted shock properties of Nb 2 O 5 polymorphs are also compared to their Nb and NbO 2 counterparts to assess the impact of niobium oxidation on shock response.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

FINETUNA: fine-tuning accelerated molecular simulations

Abstract Progress towards the energy breakthroughs needed to combat climate change can be significantly accelerated through the efficient simulation of atomistic systems. However, simulation techniques based on first principles, such as density functional theory (DFT), are limited in their practical use due to their high computational expense. Machine learning approaches have the potential to approximate DFT in a computationally efficient manner, which could dramatically increase the impact of computational simulations on real-world problems. However, they are limited by their accuracy and the cost of generating labeled data. Here, we present an online active learning framework for accelerating the simulation of atomic systems efficiently and accurately by incorporating prior physical information learned by large-scale pre-trained graph neural network models from the Open Catalyst Project. Accelerating these simulations enables useful data to be generated more cheaply, allowing better models to be trained and more atomistic systems to be screened. We also present a method of comparing local optimization techniques on the basis of both their speed and accuracy. Experiments on 30 benchmark adsorbate-catalyst systems show that our method of transfer learning to incorporate prior information from pre-trained models accelerates simulations by reducing the number of DFT calculations by 91%, while meeting an accuracy threshold of 0.02 eV 93% of the time. Finally, we demonstrate a technique for leveraging the interactive functionality built in to Vienna ab initio Simulation Package (VASP) to efficiently compute single point calculations within our online active learning framework without the significant startup costs. This allows VASP to work in tandem with our framework while requiring 75% fewer self-consistent cycles than conventional single point calculations. The online active learning implementation, and examples using the VASP interactive code, are available in the open source FINETUNA package on Github.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Machine-learned quantum molecular dynamics calculations of warm dense equation of state and ionic transport coefficients of deuterated water

White dwarf models require accurate equations of state and ionic transport coefficients in the warm dense matter regime, where kinetic theory models and tabulated equations of state are often inaccurate. In this work, spectral-partitioned density functional theory and machine-learned interatomic potentials are combined to perform large-scale, first-principles quantum molecular dynamics simulations of deuterated water (D 2 O) near the principal Hugoniot. This approach retains Kohn-Sham accuracy while achieving orders-of-magnitude speedup, yielding converged equation of state and transport properties over a broad pressure and temperature range. The results reveal the thermodynamic conditions under which ionic transport models for interdiffusivity and shear viscosity converge and identify those in closest agreement with density functional theory benchmarks at temperatures in the warm dense matter regime. The present framework extends first-principles transport calculations to higher temperatures than previously achieved, and provides an efficient, scalable, and general approach for studying transport properties in complex multicomponent mixtures.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Ion Solvation and Transport in Narrow Carbon Nanotubes: Effects of Polarizability, Cation-π Interaction, and Confinement

Understanding ion solvation and transport under confinement is critical for a wide range of emerging technologies, including water desalination and energy storage. While molecular dynamics (MD) simulations have been widely used to study the behavior of confined ions, considerable deviations between simulation results depending on the specific treatment of intermolecular interactions remain. In the following, we present a systematic investigation of the structure and dynamics of two representative solutions, that is, KCl and LiCl, confined in narrow carbon nanotubes (CNTs) with a diameter of 1.1 and 1.5 nm, using a combination of first-principles and classical MD simulations. Our simulations show that the inclusion of both polarization and cation-π interactions is essential for the description of ion solvation under confinement, particularly for large ions with weak hydration energies. Beyond the variation in ion solvation, we find that cation-π interactions can significantly influence the transport properties of ions in CNTs, particularly for KCl, where our simulations point to a strong correlation between ion dehydration and diffusion. Finally, our study highlights the complex interplay between nanoconfinement and specific intermolecular interactions that strongly control the solvation and transport properties of ions.

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

Homogeneous ice nucleation in an ab initio machine-learning model of water

Molecular simulations have provided valuable insight into the microscopic mechanisms underlying homogeneous ice nucleation. While empirical models have been used extensively to study this phenomenon, simulations based on first-principles calculations have so far proven prohibitively expensive. Here, we circumvent this difficulty by using an efficient machine-learning model trained on density-functional theory energies and forces. We compute nucleation rates at atmospheric pressure, over a broad range of supercoolings, using the seeding technique and systems of up to hundreds of thousands of atoms simulated with ab initio accuracy. The key quantity provided by the seeding technique is the size of the critical cluster (i.e., a size such that the cluster has equal probabilities of growing or melting at the given supersaturation), which is used together with the equations of classical nucleation theory to compute nucleation rates. We find that nucleation rates for our model at moderate supercoolings are in good agreement with experimental measurements within the error of our calculation. We also study the impact of properties such as the thermodynamic driving force, interfacial free energy, and stacking disorder on the calculated rates.

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

First-principles-integrated Study of Blazar Synchrotron Radiation and Polarization Signatures from Magnetic Turbulence

Blazar emission is dominated by nonthermal radiation processes that are highly variable across the entire electromagnetic spectrum. Turbulence, which can be a major source of nonthermal particle acceleration, can widely exist in the blazar emission region. The Turbulent Extreme Multi-Zone (TEMZ) model has been used to describe turbulent radiation signatures. Recent particle-in-cell (PIC) simulations have also revealed the stochastic nature of the turbulent emission region and particle acceleration therein. However, radiation signatures have not been systematically studied via first-principles-integrated simulations. In this paper, we perform combined PIC and polarized radiative transfer simulations to study synchrotron emission from magnetic turbulence in the blazar emission region. We find that the multiwavelength flux and polarization are generally characterized by stochastic patterns. Specifically, the variability timescale and average polarization degree (PD) are governed by the correlation length of the turbulence. Interestingly, magnetic turbulence can result in polarization angle swings with arbitrary amplitudes and duration, in either direction, that are not associated with changes in flux or PD. Surprisingly, these swings, which are stochastic in nature, can appear either bumpy or smooth, although large-amplitude swings (>180°) are very rare, as expected. Our radiation and polarization signatures from first-principles-integrated simulations are consistent with the TEMZ model, except that in the latter, there is a weak correlation, with zero lag, between flux and degree of polarization.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Estimating ionization states and continuum lowering from ab initio path integral Monte Carlo simulations for warm dense hydrogen

Warm dense matter (WDM) is an active field of research, with applications ranging from astrophysics to inertial confinement fusion. Ionization degree and continuum lowering are important quantities to understand how materials behave under these conditions, but can be difficult to diagnose since experimental campaigns are limited and often require model-dependent analysis. This is especially true for hydrogen, which has a comparably low scattering cross section, making high-quality data particularly difficult to obtain. Consequently, building equation of state tables often relies on simulations in combination with untested approximations to extract properties from experiments. Here, we investigate an approach for extracting the ionization potential depression and ionization degree—quantities which are otherwise not directly accessible from the physical model—from first-principles path integral Monte Carlo (PIMC) simulations utilizing a chemical model. In contrast to experimental measurements, where noise and nonequilibrium effects add to the uncertainty of the inferred parameters, PIMC simulations provide a clean signal with well-defined thermodynamic conditions. Comparisons against commonly used models show a qualitative agreement, but we find deviations primarily for the high-density and high-temperature cases. We also demonstrate the decreasing sensitivity of the dynamic structure factor with respect to both ionization and continuum lowering for increasing scattering angles in x-ray Thomson scattering experiments. Our work has important implications for the design of future experiments, but also offers qualitative understanding of structure factors and the imaginary-time correlation function obtained from first-principles quantum Monte Carlo simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Understanding dynamic properties of materials using neutron spectroscopy and atomistic simulation

Recent developments in first-principles lattice dynamics and classical force field based molecular dynamics are revolutionising the field of neutron spectroscopy. Herein we present a short review of these methods, their critical role in the supporting of cutting-edge experiments, and how they are improved by matching experimental data. We begin with a brief overview of how lattice dynamics calculations can be compared to inelastic neutron scattering (INS) and molecular dynamics simulations to both INS and quasi-elastic neutron scattering (QENS). We then provide a series of exemplar applications where lattice dynamics and molecular dynamics have been used in conjunction with neutron spectroscopy to bring significant understanding to topical areas of materials science namely: (i) lattice dynamics and INS for the study of hybrid organic-inorganic perovskites (ii) lattice dynamics and INS for the study of flexible porous solids and (iii) molecular dynamics and QENS for probing molecular behaviour in zeolite catalysis. In all three cases, the understanding gained through the synergy of experiment and computation would have been significantly reduced using either in isolation. Finally, we consider the current state of the art, describing outstanding challenges and suggesting future directions in this exciting and fertile area of physical science.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

Jet Noise Prediction for Chevron Nozzle Flows with Wall-Modeled Large-Eddy Simulation

This paper presents results from ongoing research on jet noise prediction with wall-modeled large eddy simulations (WMLES) performed with the LAVA computational framework. In particular the focus of this study is on mixing enhancements from a single stream chevron nozzle at Reynolds number of 1×10(exp 6). Although the concept of chevron nozzles to reduce jet noise is not new, our understanding of its impact on the overall noise is still not well understood. As a first step towards predicting noise reduction due to mixing enhancement concepts from first principles with WMLES, we simulate the noise generated by a single stream chevron nozzle SMC001 as well as its equivalent axisymmetric round jet nozzle SMC000. Detailed comparisons are made with a dedicated experiment conducted at NASA’s Glenn Research Center and good agreement was achieved. Two different approaches to introduce a turbulent boundary layer were compared but show no major impact on the results. This is especially important given future work were multi-stream nozzles are considered and extended costs of resolving the internal BL would have a bigger cost impact. A permeable Ffowcs Williams Hawkings (FWH) surface enclosing the jet is used to predict far-field noise from the simulated flow-field and excellent comparison to microphone array measurements is achieved within the resolved frequency band. Sensitivity of far-field noise predictions to grid resolution is systematically documented. Near-field comparisons to PIV data shows great agreement for both velocity and normal stresses, however a systematic TKE overshoot at the nozzle exit is seen in the shear-layer. The paper shows a cost reduction of an order of magnitude compared to an earlier study of this configuration due to algorithmic and software improvements and demonstrates that WMLES can be used as a cost-competitive approach for jet noise predictions.

CST↗
NASA NTRSDate not supplied

Advanced Materials for the Lunar Surface: Multiscale Computational Design of Refractory Alloys and Carbides

Emerging operational environments, such as the lunar surface, present novel challenges for NASA and drive the need for advanced materials in applications like fission surface power systems. To address these demands, computational materials science is rapidly evolving to augment or replace costly and hazardous empirical testing. Although materials selection at NASA remains predominantly experimentally driven, advanced simulation methodologies are being steadily integrated into the engineering lifecycle. This work details the application of multiscale simulation techniques—including first-principles calculations, CALPHAD, dislocation dynamics, and molecular dynamics—at NASA's Ames Research Center to evaluate advanced materials for extreme environments. First, we present contributions to the Space Nuclear Propulsion Project. Be-cause propellant channel coatings in nuclear thermal rockets must withstand high-pressure, high-temperature hydro-gen, optimizing these materials is critical. First-principles calculations were employed to establish a rigorous quantitative and qualitative understanding of the behavior of the refractory carbides ZrC, NbC, and their mixtures in high-enthalpy hydrogen environments. This necessitated the generation of high-fidelity thermodynamic models for both stoichiometric and carbon-depleted carbides, both with and without the presence of hydrogen. Furthermore, we highlight efforts under the Refractory Alloy Additive Manufacturing Build Optimization (RAAMBO) project, where existing and novel alloy compositions were assessed for additive manufacturing printability and subsequent performance in applications such as heat pipes and rocket nozzle extensions. This was accomplished through a comprehensive multiscale simulation framework that bridged the gap from the nanometer to the millimeter scale. Across both initiatives, rigorous validation against empirical data was prioritized. By systematically employing a verified and validated computational frame-work, we demonstrate how simulation effectively supports multidisciplinary engineering efforts, builds project-wide confidence, and drives critical materials development.

computational materials↗
DOE OSTI2021

Ab initio ultrafast spin dynamics in solids

Spin relaxation and decoherence is at the heart of spintronics and spin-based quantum information science. Currently, theoretical approaches that can accurately predict spin relaxation of general solids including necessary scattering pathways and are capable of nanosecond to millisecond simulation time are urgently needed. We present a first-principles real-time density-matrix approach based on Lindblad dynamics to simulate ultrafast spin dynamics for general solid-state systems. Through the complete first-principles descriptions of pump, probe, and scattering processes including electron-phonon, electron-impurity, and electron-electron scatterings with self-consistent electronic spin-orbit couplings, our method can directly simulate the ultrafast pump-probe measurements for coupled spin and electron dynamics over nanoseconds at any temperatures and doping levels. We first apply this method to a prototypical system GaAs and obtain excellent agreement with experiments. We found that the relative contributions of different scattering mechanisms and phonon modes differ considerably between spin and carrier relaxation processes. In sharp contrast to previous work based on model Hamiltonians, we point out that the electron-electron scattering is negligible at room temperature but becomes dominant at low temperatures for spin relaxation in n-type GaAs. We further examine ultrafast dynamics in novel spin-valleytronic materials: monolayer and bilayer WSe 2 with realistic defects. We find that spin relaxation is highly sensitive to local symmetry and chemical bonds around defects. For the bilayer WSe 2 , we identify the scattering pathways in ultrafast dynamics and determine relevant dynamical properties, essential to its utilization of unique spin-valley-layer locking effects. In conclusion, our work provides a predictive computational platform for spin dynamics in solids, which has potential for designing new materials ideal for spintronics and quantum information technology.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗