Search NASA⌕ Search

SEARCH · Search NASA

Results for “DFT 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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Simulation of multi-shell fullerenes using Machine-Learning Gaussian Approximation Potential

Multi-shell fullerenes ”buckyonions ” were simulated, starting from initially random configurations, using a density-functional-theory (DFT)-trained machine-learning carbon potential within the Gaussian Approximation Potential (GAP) Framework [Volker L. Deringer and Gábor Csányi, Phys. Rev. B 95, 094203 (2017)]. Fullerenes formed from seven different system sizes, ranging from 60 ~ 3774 atoms, were considered. The buckyonions are formed by clustering and layering starting from the outermost shell and proceeding inward. Inter-shell cohesion is partly due to interaction between delocalized π electrons protruding into the gallery. The energies of the models were validated ex post facto using density functional codes, VASP and SIESTA , revealing an energy difference within the range of 0.02 - 0.08 eV/atom after conjugate gradient energy convergence of the models was achieved with both methods.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Additively manufactured Ni-20Cr to V functionally graded material: Computational predictions and experimental verification of phase formations

Throughout this work, a database for the Cr-Ni-V system was constructed by modeling the binary Cr-V and ternary Cr-Ni-V systems using the CALPHAD approach aided by density functional theory (DFT)-based first-principles calculations and ab initio molecular dynamics (AIMD) simulations. To validate this new database, a functionally graded material (FGM) using Ni-20Cr and elemental V was fabricated using directed energy deposition additive manufacturing (DED AM) and experimentally characterized. The deposited Ni-20Cr was pure fcc phase, while increasing the amount of V across the gradient resulted in the formation of sigma phase, followed by the bcc phase. The experimentally measured phase data was compared with computational predictions made using a Cr-Ni-V thermodynamic database from the literature as well as the database developed in the present work. The newly developed database was shown to better predict the experimentally observed phases due to its accurate modeling of binary systems within the database and the ternary liquid phase, which is critical for accurate Scheil calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Helium interaction with solutes and impurities in neutron-irradiated nanostructured ferritic alloys: A first principles study

Density functional theory calculations are performed to explore the binding between He and alloying solutes, impurities, and transmutation products expected in neutron irradiated nanostructured ferritic alloys (NFAs), here 14YWT is taken as an example. Elements that exhibit significant binding (attraction) with an interstitial He are Y (binding energy = 0.46 eV), Mg (0.32), O (0.33), Ti (0.16), and C (0.15). Those that provide significant binding to a substitutional He are O (1.44), Y (1.24), N (0.73), H (0.56), Mg (0.52), Ti (0.34), Si (0.34), C (0.33), Al (0.32), Ni (0.26), Ta (0.23), and Mn (0.16). The presence of these elements in Fe matrix could reduce the transport of He towards oxide particles, dislocations, and internal boundaries, and could promote He bubble nucleation in the matrix. For convenience, we compile existing binding energy data of He with He n and He n V (He-vacancy) clusters. Dissociation pathway analysis reveals that, in general, the most likely dissociation of a He n V cluster is by a sequential emission of individual He atoms. Furthermore, larger bubbles are more prone to dissociation than smaller ones. In addition, we estimate the binding energy (segregation energy) of He in bulk Y 2 Ti 2 O 7 (YTO) single crystal, YTO/Fe interface, and YTO particle embedded in Fe, with respect to interstitial He in Fe, from existing formation energies of He in these structures. We also compile available data of He binding with Fe self-interstitial atom (SIA), SIA clusters, and edge and screw dislocations. Note that given the absence of DFT data, the binding with SIA clusters and dislocations are gathered from simulations with empirical potentials. Finally, the data presented in this paper is important to inform multiscale simulations of He bubble accumulation.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Modelling single atom catalysts for water splitting and fuel cells: A tutorial review

Here, in this tutorial review, we report the state-of-the-art of the modeling approaches of Single Atom Catalysts (SAC) for water splitting and fuel cells reactions. The discussion applies for Hydrogen Evolution Reaction (HER), Oxygen Reduction Reaction (OER), Hydrogen Oxidation Reaction (HOR), and Oxygen Reduction Reaction (ORR). The main scope of this work is to underline the relevant aspects of SACs modelling. On the one hand, the review could help computational chemists aiming to start the study of SACs. On the other hand, experimentalists could find the critical analysis of DFT results of interest to understand better the strengths and weaknesses of simulations, and how to interpret computational results. After an introductory section of SACs, we start by briefly presenting the state-of-the-art methodologies. Then, we analyze the critical aspects for a reliable prediction of the electronic properties, and we discuss the robustness of the structural models and ways to validate them. Furthermore, we discuss the main approaches to predict catalytic activity and selectivity, which is the final goal of the computational catalysis. We conclude this review with a critical analysis of the current challenges in the field, and the main limitations of the modeling approaches that are described.

25 ENERGY STORAGE↗
DOE OSTI2020

A Bayesian analysis of nuclear deformation properties with Skyrme energy functionals

In spite of numerous scientific and practical applications, there is still no comprehensive theoretical description of the nuclear fission process based solely on protons, neutrons and their interactions. Furthermore, the most advanced simulations of fission are currently carried out within nuclear density functional theory (DFT). In spite of being fully quantum-mechanical and rooted in the theory of nuclear forces, DFT still depends on a dozen or so parameters characterizing the energy functional. Calibrating these parameters on experimental data results in uncertainties that must be quantified for applications. This task is very challenging because of the high computational cost of DFT calculations for fission. In this paper, we use Gaussian processes to build emulators of DFT models in order to quantify and propagate statistical uncertainties of theoretical predictions for a range of nuclear deformations relevant to describing the fission process.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Towards a Deeper Fundamental Understanding of (Al,Sc)N Ferroelectric Nitrides

Density functional theory (DFT) calculations, within the virtual crystal alloy approximation, are performed, along with the development of a Landau-type model employing a symmetry-allowed analytical expression of the internal energy and having parameters determined from first principles, to investigate properties and energetics of Al1-xScxN ferroelectric nitrides in their hexagonal forms. These DFT computations and this model predict the existence of two different types of minima, namely, the fourfold-coordinated wurtzite (WZ) polar structure and a five-fold coordinated paraelectric hexagonal phase (denoted as H5), for any Sc composition up to 40%. The H5 minimum progressively becomes the lowest-energy state within hexagonal symmetry as the Sc concentration increases from 0 to 0.4. Furthermore, the model points to several key findings. Examples include the crucial role of the coupling between polarization and strains to create the WZ minimum, in addition to polar and elastic energies, and that the origin of the H5 state overcoming the WZ phase as the global minimum within hexagonal symmetry when increasing the Sc composition mostly lies in the compositional dependency of only two parameters-one linked to the polarization and another one being purely elastic in nature. Other examples are that forcing Al1-xScxN systems to have no or a weak change in lattice parameters when heating them allows us to reproduce their finite-temperature polar properties well and that a value of the axial ratio close to that of the ideal WZ structure implies a large polarization at low temperatures but not necessarily at high temperatures because of the ordered-disordered character of the temperature-induced formation of the WZ state. Such findings should allow for a better fundamental understanding of (Al,Sc)N ferroelectric nitrides, which may be used to design efficient devices having, e.g., low operating voltages.

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

Surface Wettability Analysis from Adsorption Energy and Surface Electrical Charge

Surface wettability is determined by the attraction of a liquid phase to a solid surface. It is typically quantified by using contact angle measurements at mineral surfaces in the case of the flotation of mineral particles. Contact angle research to describe wettability has been investigated at different scales by sessile drop measurements, molecular dynamic simulation, and atomic force microscopy. In this study, the density functional theory (DFT) was employed for predicting the surface free energy and contact angles of a well-known hydrophobic phyllosilicate mineral talc and a well-known hydrophilic phyllosilicate mineral muscovite based on the calculated interfacial energy and surface charge. The results revealed that the predicted contact angle at the atomic scale was larger than the experimental value, and identified two interactions: electrostatic interaction and hydrogen bonding, between the hydrophilic muscovite surface and the water layer, while a water-exclusion zone of 3.346 Å was found between the hydrophobic talc surface and the first water layer. This investigation gives a new perspective for wettability determination at the atomic scale.

Geochemistry & Geophysics↗
DOE OSTI2023

Kohn–Sham accuracy from orbital-free density functional theory via Δ-machine learning

Here, we present a Δ-machine learning model for obtaining Kohn–Sham accuracy from orbital-free density functional theory (DFT) calculations. In particular, we employ a machine-learned force field (MLFF) scheme based on the kernel method to capture the difference between Kohn–Sham and orbital-free DFT energies/forces. We implement this model in the context of on-the-fly molecular dynamics simulations and study its accuracy, performance, and sensitivity to parameters for representative systems. We find that the formalism not only improves the accuracy of Thomas–Fermi–von Weizsäcker orbital-free energies and forces by more than two orders of magnitude but is also more accurate than MLFFs based solely on Kohn–Sham DFT while being more efficient and less sensitive to model parameters. We apply the framework to study the structure of molten Al 0.88 Si 0.12 , the results suggesting no aggregation of Si atoms, in agreement with a previous Kohn–Sham study performed at an order of magnitude smaller length and time scales.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Accurate machine-learning predictions of coercivity in high-performance permanent magnets

Increased demand for high-performance permanent magnets in the electric vehicle and wind-turbine industries has prompted the search for cost-effective alternatives. Discovering magnetic materials with the desired intrinsic and extrinsic permanent magnet properties presents a significant challenge to researchers because of issues with the global supply of rare-earth elements, material stability, and a low maximum magnetic energy product BH max . While first-principles density functional theory (DFT) predicts materials’ magnetic moments, magnetocrystalline anisotropy constants, and exchange interactions, it cannot compute extrinsic properties such as coercivity (H c ). Although it is possible to calculate H c theoretically with micromagnetic simulations, the predicted value is larger than the experiment by almost an order of magnitude due to the Brown paradox. To circumvent these issues, we employ machine-learning (ML) methods on an extensive database obtained from experiments, DFT calculations, and micromagnetic modeling. The use of a large experimental dataset enables realistic H c predictions for materials such as Ce-doped Nd 2 ⁢Fe 14 ⁢B, comparing favorably against micromagnetically simulated coercivities. Remarkably, our ML model accurately identifies uniaxial magneto-crystalline anisotropy as the primary contributor to H c . With DFT calculations, we predict the Nd-site-dependent magnetic anisotropy behavior in Nd 2 ⁢Fe 14 ⁢B, confirming that Nd 4⁢g sites mainly contribute to uniaxial magnetocrystalline anisotropy, and also calculate the Curie temperature (T c ). Finally, both calculated results are in good agreement with the experiments. The coupled experimental dataset and ML modeling with DFT input predict H c with far greater accuracy and speed than was previously possible using micromagnetic modeling. Further, we reverse engineer the grain-boundary and intergrain exchange coupling with micromagnetic simulations by employing the ML predictions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Developing reliable machine learning interatomic potential for Fe–Cr–Ni austenitic alloys

Gaining atomistic understanding of mechanical behavior of heat-resistant structural materials such as Fe–Cr–Ni-based alloys requires an approach with an accuracy close to density functional theory (DFT) that considers the intrinsic properties of the bulk lattice and important defects such as stacking faults, grain boundaries, and surfaces. This work aims to develop reliable machine learning interatomic potential (MLIAP) at cross-scale for Fe–Cr–Ni ternary alloys with a focus on the face-centered-cubic (fcc) solid solution structure. Leveraging the advantages of moment tensor potentials, which typically necessitate a relatively small training dataset and enable rapid calculations using the large-scale atomic/molecular massively parallel simulator package, we ensure the stability and accuracy of the trained potentials. Important defects such as stacking faults, grain boundaries, and surfaces for wide-range compositions are investigated. Structural, thermal, elastic, and defect properties are determined from molecular dynamics simulations comprising several thousand atoms, generated via canonical Monte Carlo simulations guided by the trained potential. The trained potential allows efficient atomic simulations of structural, thermal, and mechanical properties of fcc Fe–Cr–Ni solid solution alloys as a function of composition and temperature. Therefore, the MLIAP approach represents a major advancement from DFT calculations that are limited to small simulation sizes and traditional molecular dynamics simulations using relatively low accuracy potentials. Furthermore, this work outlines a practical foundation for further investigating the structural evolution and mechanical behavior of austenitic stainless steel and nickel-based alloys in a wide array of applications in extreme environments.

Crystal structure↗
DOE OSTI2022

Final Technical Report for "Cloud-based Low-Scaling Quantum Chemistry Simulations for Materials"

The objective of phase I was to develop the basic infrastructure for performing efficient hybrid DFT calculations on solid materials with a turnaround time amenable to the use in large-scale data-driven approaches. To fulfill the goal we have developed a novel algorithm that significantly reduces the cost of exchange matrix formation. The algorithm does so by making use of three ingredients (a) interpolative decomposition of the electron integrals (b) robust pseudospectral method and (c) occ-RI approach. Our published work demonstrates that the algorithm is orders of magnitude faster than any other hybrid-DFT method and is on the order of only 3-4 times slower than pure DFT. All steps of the algorithm developed during phase I can be accelerated to the point that, for large enough systems, the diagonalization of the Fock matrix will become the most expensive step. To treat such large systems we have created an interface to the ASCR-funded PEXSI library. The PEXSI method is used to compute the density matrix directly from the Fock matrix in a manner that preserves the sparsity of the local representation.

Shiozaki, Toru↗
DOE OSTI2023

Probing the self-ionization of liquid water with ab initio deep potential molecular dynamics

The chemical equilibrium between self-ionized and molecular water dictates the acid–base chemistry in aqueous solutions, yet understanding the microscopic mechanisms of water self-ionization remains experimentally and computationally challenging. Herein, Density Functional Theory (DFT)–based deep neural network (DNN) potentials are combined with enhanced sampling techniques and a global acid–base collective variable to perform extensive atomistic simulations of water self-ionization for model systems of increasing size. The explicit inclusion of long-range electrostatic interactions in the DNN potential is found to be crucial to accurately reproduce the DFT free energy profile of solvated water ion pairs in small (64 and 128 H 2 O) cells. The reversible work to separate the hydroxide and hydronium to a distance S is found to converge for simulation cells containing more than 500 H 2 O, and a distance of ∼ 8 Å is the threshold beyond which the work to further separate the two ions becomes approximately zero. The slow convergence of the potential of mean force with system size is related to a restructuring of water and an increase of the local order around the water ions. Calculation of the dissociation equilibrium constant illustrates the key role of long-range electrostatics and entropic effects in the water autoionization process.

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

Time-dependent density-functional theory study on nonlocal electron stopping for inertial confinement fusion

Understanding laser–target coupling is of the utmost importance for achieving high performance in laser-direct-drive (LDD) inertial confinement fusion (ICF) experiments. Thus, accurate modeling of electron transport and deposition through ICF-relevant materials and conditions is necessary to quantify the total thermal conduction and ablation. The stopping range is a key transport quantity used in thermal conduction models; in this work, we review the overall role that the electron mean free path (MFP) plays in thermal conduction and hydrodynamic simulations. The currently used modified Lee–More model employs various physics approximations. We discuss a recent model that uses time-dependent density functional theory (TD-DFT) to eliminate these approximations in both the calculation of the electron stopping power and corresponding MFP in conduction zone polystyrene (CH) plasma. In general, the TD-DFT calculations showed a larger MFP (lower stopping power) than the standard modified Lee–More model. Using the TD-DFT results, an analytical model for the electron deposition range, λTD−DFT(ρ,T,K), was devised for CH plasmas between ρ=[0.05−1.05] g/cm3, kBT=[100−1000] eV. We implemented this model into LILAC, for simulations of a National Ignition Facility-scale LDD implosion and compared key physics quantities to ones obtained by simulations using the standard model. The implications of the obtained results and the path moving forward to calculate this same quantity in conduction-zone deuterium–tritium plasmas are further discussed, to hopefully close the understanding gap for laser target coupling in LDD-ICF simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Identifying factors controlling the selective ethane dehydrogenation on Pt-based catalysts from DFT based micro-kinetic modeling

Herein, a detailed EDH reaction mechanism on Pt(111), Pt(211), Pt 3 Sn(111) and Pt 3 Sn(211) surfaces has been systematically computed with DFT. A mean-field micro-kinetic model has been established to simulate the product distribution and selectivity on each surface at experimentally relevant conditions. We have uncovered a number of different factors controlling the catalytic performance of Pt-based EDH catalysts, thus leading to a better understanding of experimentally detected phenomenon. In this study, the (111) and (211) model surfaces of Pt and Pt 3 Sn catalysts were chosen to represent the terrace and step-edge sites when simulating the EDH reaction mechanisms.

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

Additively manufactured refractory high-entropy alloys with superior radiation resistance

Refractory high-entropy alloys (RHEAs) are promising candidates for next-generation nuclear and high-temperature applications. Among many approaches to manufacture RHEAs, additive manufacturing (AM) represents the most recent and advanced metal manufacturing method which allows near-net-shape manufacturing to reduce material waste and post-processing time. However, performance of AM RHEAs under complex irradiation conditions remains largely unexplored. Here, in this study, we demonstrate for the first time the response of directed energy deposition (DED) AM quaternary RHEAs (HfTaVW, CrTaVW) subjected to sequential dual-beam ion irradiation, consisting of helium pre-implantation followed by high-dose heavy ion bombardment. Compositions of DED AM RHEAs were selected using Monte Carlo (MC) simulations based on a cluster expansion (CE) Hamiltonian parameterized by density functional theory (DFT). Post-irradiation microstructural characterization revealed that the AM RHEA maintained remarkable stability, with suppressed helium bubble growth and reduced defect accumulation compared to conventional alloys. Even at high doses (∼100 dpa), the alloy exhibited no void swelling, a low density of dislocation loops, and no evidence of severe degradation. These results highlight the intrinsic ability of AM-derived microstructures and multicomponent chemistry to synergistically mitigate irradiation effects. Our findings establish AM RHEAs as a class of materials with superior resistance to radiation damage under conditions relevant to advanced fusion and fission environments and demonstrate the importance of sequential ion beam studies in evaluating their long-term performance.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Incorporating Electronic Information into Machine Learning Potential Energy Surfaces via Approaching the Ground-State Electronic Energy as a Function of Atom-Based Electronic Populations

Machine learning (ML) approximations to density functional theory (DFT) potential energy surfaces (PESs) are showing great promise for reducing the computational cost of accurate molecular simulations, but at present, they are not applicable to varying electronic states, and in particular, they are not well suited for molecular systems in which the local electronic structure is sensitive to the medium to long-range electronic environment. With this issue as the focal point, we present a new machine learning approach called "BpopNN" for obtaining efficient approximations to DFT PESs. Conceptually, the methodology is based on approaching the true DFT energy as a function of electron populations on atoms; in practice, this is realized with available density functionals and constrained DFT (CDFT). Additionally, the new approach creates approximations to this function with neural networks. These approximations thereby incorporate electronic information naturally into a ML approach, and optimizing the model energy with respect to populations allows the electronic terms to self-consistently adapt to the environment, as in DFT. We confirm the effectiveness of this approach with a variety of calculations on LinHn clusters.

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

Development of a Transferable Density-Functional Tight-Binding Model for Organic Molecules at the Water/Platinum Interface

A computationally efficient and transferable approach for modeling reactions at metal/water interfaces could significantly accelerate our understanding and ultimately the development of new catalytic transformations, particularly in the context of the emerging field of biomass conversion. Here, we present a parametrization of Pt–X (X = H, O, C) density-functional tight-binding (DFTB) for addressing this need. We first constructed Pt–H, Pt–O, and Pt–C repulsive potential splines. These pairwise parameters were then augmented to include many-body interactions using the Chebyshev Interaction Model for Efficient Simulation (ChIMES). We compare the geometrical and energetic performances of both DFTB and DFTB/ChIMES methods with DFT reference data across a variety of organic molecules at the platinum surface from nanoparticles to single-crystal surfaces. DFTB shows limited transferability between extended crystal surfaces and small nanoparticles. This transferability is significantly improved through the introduction of three-body interactions with Pt in DFTB/ChIMES, which provides consistent results across various systems, with reductions in the RMSD from around 30 kcal/mol in DFTB to around 10 kcal/mol. We demonstrate the stability and reliability of the obtained parameters by performing metadynamic simulations for the adsorption of phenol on Pt(111). We observe that DFTB itself is undersolvating the surface, leading to only one or two chemisorbed water molecules in a c(4 × 6) unit cell. In contrast, DFTB/ChIMES leads to a coverage of about 0.5 ML and successfully captures the chemisorbed mode of phenol at both the solid/liquid and the solid/gas interfaces. Furthermore, in agreement with experimental measurements, the adsorption at the solid/liquid interface is significantly weaker than that at the solid/gas interface. As a result, we highlight that even with DFTB, where we can accumulate dynamics for more than 1 ns for a given system, the simulations are not fully converged.

Adsorption↗