Search NASA⌕ Search

SEARCH · Search NASA

Results for “Density functional theory development”

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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Experimental trends and theoretical descriptors for electrochemical reduction of carbon dioxide to formate over Sn-based bimetallic catalysts

The electrochemical carbon dioxide reduction reaction (CO 2 RR) using renewable energy sources is a promising solution for mitigating CO 2 emissions. In particular, CO 2 RR to formate represents a commercially profitable target. However, a comprehensive understanding of the catalytic mechanisms of Sn-based catalysts under reaction conditions, including the real-time structural evolution of catalysts and the role of all key reaction intermediates in influencing the CO 2 RR selectivity, is still lacking. The current study reports a framework to study the selectivity preference of Sn-based bimetallic catalysts using a combination of electrochemical measurements, in situ characterization, and density functional theory (DFT) calculations. The addition of a second metal (Co, Ni, Ag, Zn, Ga, Bi) was found to play a vital role in affecting the CO 2 RR performance. In situ X-ray absorption near edge structure (XANES) measurements revealed a dynamic evolution in the Sn valence state induced by different secondary metals. A multidimensional descriptor involving all the key reaction intermediates was developed to assess formate selectivity using a 2-dimensional volcano plot. Finally, this research offers an effective framework for understanding CO 2 RR catalytic selectivity by considering both the real-time structural evolution of catalysts and all the key intermediates involved.

58 GEOSCIENCES↗
DOE OSTI2022

Thermal Neutron Scattering Law for Beryllium Hydride and Critical Mass Calculations

The thermal neutron scattering law (TSL) for crystalline beryllium hydride (BeH 2 ) is developed from first-principles ab initio lattice dynamics calculations and the impact of neutron thermalization in this material on critical mass is estimated. BeH 2 has a body-centered orthorhombic crystal structure with 12 molecules per unit cell and a theoretical density of 0.755 g/cm 3 . The vibrational (phonon) densities of states for H and Be bound in BeH 2 are determined using VASP density functional theory and PHONON lattice dynamics calculations. The TSLs for H bound in BeH 2 , H(BeH 2 ), and Be bound in BeH 2 , Be(BeH 2 ), are then evaluated in the incoherent approximation from the calculated H and Be partial phonon density of states using FLASSH. Finally, critical mass as a function of 235 U loading density for bare and reflected BeH 2 moderated spheres is predicted from MC21 Monte Carlo neutron transport calculations using ENDF/B-VIII.0 cross sections and the H(BeH 2 ) and Be(BeH 2 ) TSL evaluations. Comparisons are made to water (H 2 O), polyethylene (CH 2 ), and beryllium oxide (BeO) as moderators. These critical mass predictions are a refinement upon the prior work by Rao and Srinivasan that neglected thermal neutron scattering effects. The minimum critical mass of a BeH 2 moderated assembly is estimated to be 0.207 kg 235 U for a 0.20 m thick BeO reflected sphere and 0.178 kg 235 U for a 0.40 m thick BeO reflected sphere.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Incorporating Coverage-Dependent Reaction Barriers into First-Principles-Based Microkinetic Models: Approaches and Challenges

Mean-field microkinetic models (MKMs) are appealing for their relatively facile construction, computational tractability, and high-throughput catalyst screening capabilities. As such, they will continue to be a valuable tool for materials design in heterogeneous catalysis even as the field aims to describe more complex systems. Numerous prior reports have provided the groundwork for constructing first-principles-based MKMs, including the analysis of strategies for incorporating lateral interactions into thermodynamic parameters (e.g., adsorption energies). Yet, there remains a need for concerted dialogue on methods for calculating and incorporating coverage-dependent kinetic parameters into MKMs. In this Perspective, we assess strategies for doing so, including the corresponding key physical implications and computational challenges. Here, we emphasize that decoupling thermodynamic and kinetic parameters within MKMs can violate thermodynamic consistency and risk unphysical solutions. For some reactions and catalyst materials, scaling relationships can predict coverage-dependent activation energies, but there are several exceptions evident in the literature, indicating that this approach is not universally applicable and that the field could benefit from research aimed at elucidating the limitations. Conducting high-coverage transition state searches is a rigorous but computationally costly strategy, and the effects of various methods for mitigating this cost on resulting energetics have yet to be broadly explored and validated. The goal of this Perspective is to generate discussion on and inspire focused research into the physical relevance of approaches for describing coverage-dependent reaction barriers in MKMs, including the development of computationally tractable methodologies, to advance the applicability of MKMs across diverse reaction chemistries and conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Thermal Activation of Zirconium(IV) Acetylacetonate Catalysts to Enhance Polyurethane Synthesis and Reprocessing

Carbamate formation and exchange catalysts enable efficient polyurethane (PU) manufacturing, as well as emerging recycling and reprocessing methods for PU thermosets. Zirconium β-diketonate complexes, such as Zr acetylacetonate [Zr(acac) 4 ], are effective alternatives to toxic organotin catalysts that have been used for PU reprocessing. Here, we report that Zr(acac) 4 undergoes a thermally activated process in the PU network during reprocessing that transforms it into a more active carbamate exchange catalyst. This process is associated with the irreversible loss of acetylacetonate ligands and is not observed for the more sterically hindered Zr 2,2,6,6-tetramethyl-3,5-heptanedione [Zr(tmhd) 4 ] complex. Crossover experiments between PU thermoplastics indicated enhanced carbamate exchange after the thermal activation of Zr(acac) 4 in the presence of one of the PUs, whereas a sample of Zr(acac) 4 activated in the absence of the PU had no catalytic activity. Thermal gravimetric analysis suggested that this process is associated with the loss of one protonated acac ligand. Stress relaxation analysis of PU thermosets indicated a distinct change in the characteristic relaxation time associated with the thermal activation of Zr(acac) 4 at temperatures above 140 °C; no such change was observed for samples reprocessed using Zr(tmhd) 4 . Density functional theory and molecular experiments suggest that irreversible ligand exchange of acac with alkoxide or carbamate reduces the activation energy for urethane formation and reversion. Furthermore, the Zr(acac) 4 catalyst activated in the presence of a PU’s polyol precursor provided more porous and less dense PU foams compared to those made using the unactivated Zr(acac) 4 catalyst. Furthermore, these findings are important for developing improved PU synthesis and recycling processes. Thermally activating a catalyst during reprocessing may provide more nuanced control of the in-use and reprocessing characteristics of PU thermosets.

Alcohols↗
DOE OSTI2024

Degradation of performance in ICF implosions due to Rayleigh–Taylor instabilities: A Hamiltonian perspective

The Rayleigh–Taylor instability (RTI) is an ubiquitous phenomenon that occurs in inertial-confinement-fusion (ICF) implosions and is recognized as an important limiting factor of ICF performance. To analytically understand the RTI dynamics and its impact on ICF capsule implosions, we develop a first-principle variational theory that describes an imploding spherical shell undergoing RTI. The model is based on a thin-shell approximation and includes the dynamical coupling between the imploding spherical shell and an adiabatically compressed fluid within its interior. Using a quasilinear analysis, we study the degradation trends of key ICF performance metrics (e.g., stagnation pressure, residual kinetic energy, and areal density) as functions of initial RTI parameters (e.g., the initial amplitude and Legendre mode), as well as the 1D implosion characteristics (e.g., the convergence ratio). We compare analytical results from the theory against nonlinear results obtained by numerically integrating the governing equations of this reduced model. Our findings emphasize the need to incorporate polar flows in the calculation of residual kinetic energy and demonstrate that higher convergence ratios in ICF implosions lead to significantly greater degradation of key performance metrics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Thermal Neutron Scattering Research and Development at Oak Ridge National Laboratory [Slides]

This presentation provides an overview of the project, specifically describing how neutron cross sections can be broken up into several regions. It also discusses how thermal scattering law (TSL) $\textit{S(α, β)}$ is related to dynamic structure factor $\textit{S(q, E)}$ and how $\textit{S(q, E)}$ can be broken up into 2 primary components, incoherent and coherent, and how it can be calculated in several ways, including density functional theory (DFT) and molecular dynamics (MT). Additionally, it states that ENDF currently stores coherent elastic (Bragg edges), incoherent elastic and incoherent inelastic. Also discussed are the experimental capabilities, evaluations, and uncertainty quantification of the project.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Improving stability of ALD ZrN thin film coatings over U-Mo dispersion fuel

Atomic layer deposition (ALD) of ZrN is a candidate technology for coating U-Mo dispersion fuel as a diffusion barrier coating. During the early development stages of the coating the ZrN deposition showed to be mechanically unstable and ultimately resulted in spalling. Based on experimental outcomes it was found that instabillity can be eliminated through the introduction of a thin amorphous (a) Al 2 O 3 interlayer coating deposited in between the ~1-μm thick ZrN and the U-Mo. substrate. To elucidate these findings simulations were performed with density functional theory (DFT) to measure work of adhesions at different interfaces while finite element modeling (FEM) was performed to measure the residual stress distribution. DFT indicated that the ZrN coating deposited over (a)-Al 2 O 3 is ~2.5 times stronger when compared to direct deposition over U or UO 2 substrate. While calculations from FEM recognized; (1) large stress concentrations can originate from the irregularly distributed native surface oxides (UO 2 ) over U-Mo and (2) stress concentrations can be reduced if those surface UO 2 can be modified to a uniform layer (achieved after application of a 8 nm (a)-Al 2 O 3 interlayer). These results provided explanations and confirmed the role played by interlayer in eliminating the original mechanical instability of the ZrN.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Role of Undercoordinated Sites for the Catalysis in Confined Spaces Formed by Two-Dimensional Material Overlayers

Adding a two-dimensional (2D) overlayer on a metal surface is a promising route for activating reactants confined in the interfacial space. However, an atomistic understanding of the role played by undercoordinated sites of the 2D overlayer in the activation of molecules in this nanoscaled confined space is yet to be developed. In this paper, we study CO dissociation as a prototypical reaction to investigate CO activation in the confined space enclosed by Rh(111) and a monolayer of hexagonal boron nitride (h-BN). The effect of the space size (i.e., the distance between h-BN and the metal surface), the type of undercoordinated sites, and the size of the defect are explicitly studied by density functional theory with dispersion correction. The following temperature-programmed X-ray photoelectron spectroscopy measurement suggests that a small portion of the CO dissociated during the desorption, leaving the residual atomic oxygen incorporated into the h-BN lattice, which validates the theoretical prediction. Finally, the combination of theory and experiment calls for further attention to be paid to the role of undercoordinated sites in the 2D overlayers in confined systems forming potential new catalytic environments.

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

Nanoprecipitates to Enhance Radiation Tolerance in High-Entropy Alloys

The growth of advanced energy technologies for power generation is enabled by the design, development, and integration of structural materials that can withstand extreme environments, such as high temperatures, radiation damage, and corrosion. High-entropy alloys (HEAs) are a class of structural materials in which suitable chemical elements in four or more numbers are mixed to typically produce single-phase concentrated solid solution alloys (CSAs). Many of these alloys exhibit good radiation tolerance like limited void swelling and hardening up to relatively medium radiation doses (tens of displacements per atom (dpa)); however, at higher radiation damage levels (>50 dpa), some HEAs suffer from considerable void swelling limiting their near-term acceptance for advanced nuclear reactor concepts. In this study, we developed a HEA containing a high density of Cu-rich nanoprecipitates distributed in the HEA matrix. The Cu-added HEA, NiCoFeCrCu 0.12 , shows excellent void swelling resistance and negligible radiation-induced hardening upon irradiation up to high radiation doses (i.e., higher than 100 dpa). The void swelling resistance of the alloy is measured to be significantly better than NiCoFeCr CSA and austenitic stainless steels. Density functional theory simulations predict lower vacancy and interstitial formation energies at the coherent interfaces between Cu-rich nanoprecipitates and the HEA matrix. The alloy maintained a high sink strength achieved via nanoprecipitates and the coherent interface with the matrix at a high radiation dose (~50 dpa). From our experiments and simulations, the effective recombination of radiation-produced vacancies and interstitials at the coherent interfaces of the nanoprecipitates is suggested to be the critical mechanism responsible for the radiation tolerance of the alloy. Furthermore, the materials design strategy based on incorporating a high density of interfaces can be applied to high-entropy alloy systems to improve their radiation tolerance.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Molecular surface programming of rectifying junctions between InAs colloidal quantum dot solids

Heavy-metal-free III–V colloidal quantum dots (CQDs) show promise in optoelectronics: Recent advancements in the synthesis of large-diameter indium arsenide (InAs) CQDs provide access to short-wave infrared (IR) wavelengths for three-dimensional ranging and imaging. In early studies, however, we were unable to achieve a rectifying photodiode using CQDs and molybdenum oxide/polymer hole transport layers, as the shallow valence bandedge (5.0 eV) was misaligned with the ionization potentials of the widely used transport layers. This occurred when increasing CQD diameter to decrease the bandgap below 1.1 eV. Here, we develop a rectifying junction among InAs CQD layers, where we use molecular surface modifiers to tune the energy levels of InAs CQDs electrostatically. Previously developed bifunctional dithiol ligands, established for II-VI and IV-VI CQDs, exhibit slow reaction kinetics with III-V surfaces, causing the exchange to fail. Here, we study carboxylate and thiolate binding groups, united with electron-donating free end groups, that shift upward the valence bandedge of InAs CQDs, producing valence band energies as shallow as 4.8 eV. Photophysical studies combined with density functional theory show that carboxylate-based passivants participate in strong bidentate bridging with both In and As on the CQD surface. The tuned CQD layer incorporated into a photodiode structure achieves improved performance with EQE (external quantum efficiency) of 35% (>1 μm) and dark current density < 400 nA cm -2 , a >25% increase in EQE and >90% reduced dark current density compared to the reference device. This work represents an advance over previous III-V CQD short-wavelength IR photodetectors (EQE < 5%, dark current > 10,000 nA cm -2 ).

36 MATERIALS SCIENCE↗
DOE OSTI2020

Inducing metallicity in graphene nanoribbons via zero-mode superlattices

The design and fabrication of robust metallic states in graphene nanoribbons (GNRs) are challenging because lateral quantum confinement and many-electron interactions induce electronic band gaps when graphene is patterned at nanometer length scales. Recent developments in bottom-up synthesis have enabled the design and characterization of atomically precise GNRs, but strategies for realizing GNR metallicity have been elusive. Here we demonstrate a general technique for inducing metallicity in GNRs by inserting a symmetric superlattice of zero-energy modes into otherwise semiconducting GNRs. We verify the resulting metallicity using scanning tunneling spectroscopy as well as first-principles density-functional theory and tight-binding calculations. Our results reveal that the metallic bandwidth in GNRs can be tuned over a wide range by controlling the overlap of zero-mode wave functions through intentional sublattice symmetry breaking.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Methanol tolerance of atomically dispersed single metal site catalysts: mechanistic understanding and high-performance direct methanol fuel cells

Proton-exchange membrane fuel cells (PEMFCs) and direct methanol fuel cells (DMFCs) are promising power sources from portable electronic devices to vehicles. The high-cost issue of these low-temperature fuel cells can be primarily addressed by using platinum-group metal (PGM)-free oxygen reduction reaction (ORR) catalysts, in particular atomically dispersed metal–nitrogen–carbon (M–N–C, M = Fe, Co, Mn). Furthermore, a significant advantage of M–N–C catalysts is their superior methanol tolerance over Pt, which can mitigate the methanol cross-over effect and offer great potential of using a higher concentration of methanol in DMFCs. In this work, we investigated the ORR catalytic properties of M–N–C catalysts in methanol-containing acidic electrolytes via experiments and density functional theory (DFT) calculations. FeN 4 sites demonstrated the highest methanol tolerance ability when compared to metal-free pyridinic N, CoN 4 , and MnN 4 active sites. The methanol adsorption on MN 4 sites is even strengthened when electrode potentials are applied during the ORR. The negative influence of methanol adsorption becomes significant for methanol concentrations higher than 2.0 M. However, the methanol adsorption does not affect the 4e - ORR pathway or chemically destroy the FeN 4 sites. The understanding of the methanol-induced ORR activity loss guides the design of promising M–N–C cathode catalyst in DMFCs. Accordingly, we developed a dual-metal site Fe/Co–N–C catalyst through a combined chemical-doping and adsorption strategy. Instead of generating a possible synergistic effect, the introduced Co atoms in the first doping step act as “scissors” for Zn removal in metal–organic frameworks (MOFs), which is crucial for modifying the porosity of the catalyst and providing more defects for stabilizing the active FeN 4 sites generated in the second adsorption step. The Fe/Co–N–C catalyst significantly improved the ORR catalytic activity and delivered remarkably enhanced peak power densities (i.e., 502 and 135 mW cm -2 ) under H 2 –air and methanol–air conditions, respectively, representing the best performance for both types of fuel cells. Notably, the fundamental understanding of methanol tolerance, along with the encouraging DMFC performance, will open an avenue for the potential application of atomically dispersed M–N–C catalysts in other direct alcohol or ammonia fuel cells.

25 ENERGY STORAGE↗
DOE OSTI2022

Machine learning prediction and experimental verification of Pt-modified nitride catalysts for ethanol reforming with reduced precious metal loading

Ethanol is the smallest molecule containing C–O, C–C, C–H, and O–H bonds present in biomass-derived oxygenates. The development of inexpensive and selective catalysts for ethanol reforming is important towards the renewable generation of hydrogen from biomass. Transition metal nitrides (TMN) are interesting catalyst support materials that can effectively reduce precious metal loading for the catalysis of ethanol and other oxygenates. Herein theoretical and experimental methods were used to probe platinum-modified molybdenum nitride (Pt/Mo 2 N) surfaces for ethanol reforming. Computations using density-functional theory and machine learning predicted monolayer Pt/Mo 2 N to be highly active and selective for ethanol reforming. Temperature-programmed desorption (TPD) experiments verified that ethanol primarily underwent decomposition on Mo 2 N, and the reaction pathway shifted to reforming on Pt/Mo 2 N surfaces. Additionally, high-resolution electron energy loss spectroscopy (HREELS) results further indicated that while Mo2N decomposed the ethoxy intermediate by cleaving C–C, C–O, and C–H bonds, Pt-modification preserved the C–O bond, resulting in ethanol reforming.

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

Stabilization of point-defect spin qubits by quantum wells

Defect-based quantum systems in wide bandgap semiconductors are strong candidates for scalable quantum-information technologies. However, these systems are often complicated by charge-state instabilities and interference by phonons, which can diminish spin-initialization fidelities and limit room-temperature operation. Here, we identify a pathway around these drawbacks by showing that an engineered quantum well can stabilize the charge state of a qubit. Using density-functional theory and experimental synchrotron X-ray diffraction studies, we construct a model for previously unattributed point defect centers in silicon carbide as a near-stacking fault axial divacancy and show how this model explains these defects' robustness against photoionization and room temperature stability. These results provide a materials-based solution to the optical instability of color centers in semiconductors, paving the way for the development of robust single-photon sources and spin qubits.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Non-covalent planarizing interactions yield highly ordered and thermotropic liquid crystalline conjugated polymers

Controlling the multi-level assembly and morphological properties of conjugated polymers through structural manipulation has contributed significantly to the advancement of organic electronics. In this work, a redox active conjugated polymer, TPT–TT, composed of alternating 1,4-(2-thienyl)-2,5-dialkoxyphenylene (TPT) and thienothiophene (TT) units is reported with non-covalent intramolecular S⋯O and S⋯H–C interactions that induce controlled main-chain planarity and solid-state order. As confirmed by density functional theory (DFT) calculations, these intramolecular interactions influence the main chain conformation, promoting backbone planarization, while still allowing dihedral rotations at higher kinetic energies (higher temperature), and give rise to temperature-dependent aggregation properties. Thermotropic liquid crystalline (LC) behavior is confirmed by cross-polarized optical microscopy (CPOM) and closely correlated with multiple thermal transitions observed by differential scanning calorimetry (DSC). This LC behavior allows us to develop and utilize a thermal annealing treatment that results in thin films with notable long-range order, as shown by grazing-incidence X-ray diffraction (GIXD). Specifically, we identified a first LC phase, ranging from 218 °C to 107 °C, as a nematic phase featuring preferential face-on π–π stacking and edge-on lamellar stacking exhibiting a large extent of disorder and broad orientation distribution. A second LC phase is observed from 107 °C to 48 °C, as a smectic A phase featuring sharp, highly ordered out-of-plane lamellar stacking features and sharp tilted backbone stacking peaks, while the structure of a third LC phase with a transition at 48 °C remains unclear, but resembles that of the solid state at ambient temperature. Furthermore, the significance of thermal annealing is evident in the ∼3-fold enhancement of the electrical conductivity of ferric tosylate-doped annealed films reaching 55 S cm −1 . More importantly, thermally annealed TPT–TT films exhibit both a narrow distribution of charge-carrier mobilities (1.4 ± 0.1) × 10 −2 cm 2 V −1 s −1 along with a remarkable device yield of 100% in an organic field-effect transistor (OFET) configuration. This molecular design approach to obtain highly ordered conjugated polymers in the solid state affords a deeper understanding of how intramolecular interactions and repeat-unit symmetry impact liquid crystallinity, solution aggregation, solution to solid-state transformation, solid-state morphology, and ultimately device applications.

Chemistry↗
DOE OSTI2020

Electronic structure and optical properties of quantum crystals from first principles calculations in the Born–Oppenheimer approximation

In this work, we develop a formalism to accurately account for the renormalization of the electronic structure due to quantum and thermal nuclear motions within the Born–Oppenheimer approximation. We focus on the fundamental energy gap obtained from electronic addition and removal energies from quantum Monte Carlo calculations in either the canonical or grand-canonical ensembles. The formalism applies as well to effective single electron theories such as those based on density functional theory. We show that the electronic (Bloch) crystal momentum can be restored by marginalizing the total electron–ion wave function with respect to the nuclear equilibrium distribution, and we describe an explicit procedure to establish the band structure of electronic excitations for quantum crystals within the Born–Oppenheimer approximation. Based on the Kubo–Greenwood equation, we discuss the effects of nuclear motion on optical conductivity. Our methodology applies to the low temperature regime where nuclear motion is quantized and, in general, differs from the semi-classical approximation. We apply our method to study the electronic structure of C2/c-24 crystalline hydrogen at 200 K and 250 GPa and discuss the optical absorption profile of hydrogen crystals at 200 K and carbon diamond at 297 K.

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

Effective pair potential for random fcc CoCrFeMnNi alloys

The single-phase equiatomic CoCrFeMnNi alloy is a random solid solution of five elements on the face-centered cubic lattice, whose pure constituents crystallize in very different structures and exhibit diverse magnetic properties. Due to the randomness of the alloy, 80% of nearest neighbor bonds are between unlike elements and thus the details of bonding in pure structures are less important. The elastic moduli of this alloy give rise to small Cauchy pressure C 12 - C 44 , which suggests that the dominant part of bonding may be described by a simple pair potential. We test this hypothesis by developing a long-range Lennard-Jones potential in which the equilibrium crystal structures of pure constituents are taken as reference. The standard mixing rules for regular solutions are then adopted to obtain parameters for bonds between unlike elements in the quinary system. The transferability of this potential to quaternary CoCrFeNi, ternary CoCrNi, and binary FeNi alloys is investigated and the predictions compared with experiments and density functional theory calculations. By sampling over a large number of random configurations, we investigate the effect of compositional randomness on misfit volumes, energies of point defects and stacking faults, and the dislocation friction stresses experienced by moving edge and screw dislocations.

Materials Science↗
DOE OSTI2025

Developing Machine Learning Interatomic Potential for Fe-Cr-Ni Alloys

Accurate prediction of creep and fatigue behavior of stainless steel at elevated temperatures in hydrogen environment requires fundamental understanding of alloy-hydrogen interaction at cross-scale including bulk lattice and key defects such as vacancies, grain boundaries, surfaces, stacking faults, dislocations, and precipitates. This project aims to predict creep behavior of 347H stainless steel with H using machine learning interatomic potentials based on first-principles density functional theory simulations. The Moment Tensor Potentials platform is adopted for this work since it demonstrates a fine balance between model accuracy and computational efficiency. The potential is well trained based on large amount of high-fidelity density functional theory calculations. The validation is carried out by comparing various important properties including short range order, coefficient of thermal expansion, elastic properties, stacking fault energy, grain boundary energy, and surface energy. This work lays the foundation for reliable atomistic simulation of high temperature hydrogen attack of stainless steel.

density functional theory (DFT)↗