Search NASA⌕ Search

SEARCH · Search NASA

Results for “barrier function”

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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Quantal diffusion approach for multinucleon transfer processes in the 58,64 Ni + 208 Pb reactions: Toward the production of unknown neutron-rich nuclei

In recent years, substantial efforts have been made for the study of multinucleon transfer reactions at energies around the Coulomb barrier both experimentally and theoretically, aiming at the production of unknown neutron-rich heavy nuclei. It is crucial to provide reliable theoretical predictions based on microscopic theories with sufficient predictive power. Purpose: This article aims to clarify the applicability of the quantal diffusion approach based on the stochastic mean-field (SMF) theory for multinucleon transfer processes. Isotope production cross sections are evaluated for the reactions of 64 Ni + 208 Pb at E c.m. = 268 MeV and 58 Ni + 208 Pb at E c.m. = 270 MeV and are compared with available experimental data. Methods: Three-dimensional time-dependent Hartree-Fock (TDHF) calculations are carried out for a range of initial orbital angular momenta with Skyrme SLy4d functional. Quantal diffusion equations, derived based on the SMF theory, for variances and covariance of neutron and proton numbers of reaction products are solved, with microscopic drift and diffusion coefficients obtained from time evolution of occupied single-particle orbitals in TDHF. Secondary de-excitation processes, both particle evaporation and fission, are simulated by a statistical compound-nucleus de-excitation model, gemini++. Dynamics of a fast isospin equilibration process followed by a slow drift toward the mass symmetry are commonly observed, as expected. Various reaction outcomes are evaluated, including average mass and charge numbers of reaction products, total kinetic energy loss (TKEL), scattering angle, contact time, and production cross sections for primary and secondary products. By comparing with the experimental data, we find that SMF and TDHF quantitatively reproduce experimental data for few-nucleon-transfer channels around the average values. In contrast, for many-nucleon-transfer channels, we find that the SMF approach provides much better description of the experimentally measured isotopic distributions. The results underline the importance of beyond-mean-field effects, especially one-body (mean-field) fluctuations and correlations, in describing multinucleon transfer processes. Moreover, through a combined analysis of SMF with a statistical model, gemini++, we find a significant contribution of transfer-induced fission, which is consistent with the experimental observation. In some cases, the SMF approach overestimates the isotopic width, requiring further improvements of the theoretical description. Possible ways to improve the description are discussed. The SMF approach is designed to describe the quantum many-body problem according to an ensemble of mean-field trajectories, taking into account part of many-body correlations in the description. As it requires feasible computational costs comparable to the ordinary TDHF approach, together with further model improvements, it will be a promising tool in the search for optimal reaction conditions to produce yet-unknown neutron-rich heavy nuclei through the multinucleon transfer reaction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Atomistic and mesoscale simulations to determine effective diffusion coefficient of fission products in SiC

The silicon carbide (SiC) layer in tristructural isotropic (TRISO) particles serves as the barrier to prevent escape of fission products produced in the fuel kernel. Knowing the diffusion coefficient of fission products through SiC is critical to determining whether fission gas can escape from the particle. It has been observed in experiments that Ag accumulated in grain boundaries and triple junctions in SiC. It is hypothesized that grain boundary diffusion is the primary pathway by which fission products penetrate the SiC layer. In this report, the effective diffusion coefficient of the fission product Ag through the grain boundary network is calculated using a combination of atomistic and phase-field methods. The grain boundary diffusion coefficient is calculated using molecular dynamics simulations. The bulk diffusion coefficient is determined using a combination of density functional theory and nudged elastic band methods. An effective diffusion coefficient is calculated, accounting for the grain structure using a phase-field method. The effective diffusion coefficient will be incorporated into Bison and fission product release calculations are compared to available experimental data.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Mortar: An Open Testbed for Portable Building Analytics

Access to large amounts of real-world data has long been a barrier to the development and evaluation of analytics applications for the built environment. Open datasets exist, but they are limited in their span (how much data is available) and context (what kind of data is available and how it is described). Evaluation of such analytics is also limited by how the analytics themselves are implemented, often using hard-coded names of building components, points and locations, or unique input data formats. To advance the methodology for how such analytics are implemented and evaluated, we present Mortar: an open testbed for portable building analytics, currently spanning 90 buildings and containing over 9.1 billion data points. All buildings in the testbed are described using Brick, a recently developed metadata schema, providing rich functional descriptions of building assets and subsystems. We also propose a simple architecture for writing portable analytics applications that are robust to the diversity of buildings and can configure themselves based on context. We demonstrate the utility of Mortar by implementing 11 applications from the literature.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Mechanistic Insights on the Low-Temperature Oxidation of CO Catalyzed by Isolated Co Ions in N-Doped Carbon

Isolated cobalt ions on nitrogen-doped carbon (Co-N-C) catalyze CO oxidation at temperatures as low as 196 K, but the active site and mechanism for this reaction remain elusive. In this work, steady-state CO oxidation around 273 K over Co-N-C revealed nearly first order behavior in both CO and O 2 as well as a negative apparent activation energy. Isotopic transient analysis of the reaction confirmed a rapid turnover frequency and low surface coverage of adsorbed intermediates leading to CO 2 (< 10% of the Co). Results from kinetics experiments combined with quantum chemical calculations and molecular dynamics simulations are consistent with a reaction path involving weak adsorption of CO onto Co ions followed by a low barrier for the CO-assisted activation of weakly adsorbed O 2 . Furthermore, this proposed mechanism for dioxygen activation does not involve a redox cycle with the transition metal ion and may be important in other low temperature catalytic reactions involving O 2 .

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

Ultrathin Magnesium-Based Coating as an Efficient Oxygen Barrier for Superconducting Circuit Materials

Scaling up superconducting quantum circuits based on transmon qubits necessitates substantial enhancements in qubit coherence time. Over recent years, tantalum (Ta) has emerged as a promising candidate for transmon qubits, surpassing conventional counterparts in terms of coherence time. However, amorphous surface Ta oxide layer may introduce dielectric loss, ultimately placing a limit on the coherence time. In this study, a novel approach for suppressing the formation of tantalum oxide using an ultrathin magnesium (Mg) capping layer is presented. Synchrotron-based X-ray photoelectron spectroscopy studies demonstrate that oxide is confined to an extremely thin region directly beneath the Mg/Ta interface. Additionally, it is demonstrated that the superconducting properties of thin Ta films are improved following the Mg capping, exhibiting sharper and higher-temperature transitions to superconductive and magnetically ordered states. Moreover, an atomic-scale mechanistic understanding of the role of the capping layer in protecting Ta from oxidation is established based on computational modeling. Further, this work provides valuable insights into the formation mechanism and functionality of surface tantalum oxide, as well as a new materials design principle with the potential to reduce dielectric loss in superconducting quantum materials. Ultimately, the findings pave the way for the realization of large-scale, high-performance quantum computing systems.

36 MATERIALS SCIENCE↗
DOE OSTI2026

First principles band structure of interacting phosphorus and boron/aluminum δ -doped layers in silicon

Silicon can be heavily doped with phosphorus in a single atomic layer (a δ layer), significantly altering the electronic structure of the conduction bands within the material. Recent progress has also made it possible to further dope silicon with acceptor-based δ layers using either boron or aluminum, making it feasible to create devices with interacting δ layers with opposite polarity. Using density functional theory, we calculate the electronic structure of a phosphorus-based δ layer interacting with a boron or aluminum δ layer, varying the distances between the δ layers. At separations 1 nm and smaller, the dopant potentials overlap and largely cancel each other out, leading to an electronic structure closely mimicking intrinsic silicon. At separations greater than 1 nm, the two δ layers behave independently of one another, with an equivalent electronic structure to a p–n diode with an intrinsic layer taking the place of the depletion region. One mechanism for charge transfer between δ layers at larger distances could be tunneling, where we see a tunneling probability exceeding what would be seen for a standard silicon 1.1 eV triangular barrier, indicating that the interaction between delta layers may enhance tunneling compared to a traditional junction.

delta layers↗
DOE OSTI2022

Development of a web-based screening tool for ground source heat pump applications

Ground source heat pump (GSHP) technology has great potential to help the nation meet its energy and decarbonization goals, but several barriers hinder the wide application of GSHP. Important barriers include the lack of a coherent toolset for analyzing the technical feasibility and economic viability of the GSHP application. The current design and analysis methods are ineffective and require significant expertise to apply. Although building energy modeling is increasingly important in designing buildings, the tools for GSHP modeling and simulation are lacking. A web-based free-to-use tool is being developed for quick techno-economic analysis of GSHP applications in nearly any building in the United States. This tool is enabled by improvements in the calculation methodology to allow rapid sizing of borehole configurations that provide significant cost savings. The screening tool currently uses US Department of Energy (DOE) prototype building models and an extended g-function library to size ground heat exchangers and simulate the performance of GSHP systems. The team is integrating with DOE's Oak Ridge National Laboratory's AutoBEM program to automatically create a building model based on user inputs. This paper introduces the structure, components, features, and results of the web-based screening tool for GSHP applications. Future directions for further developing the tool are also discussed.

Liu, Xiaobing↗
DOE OSTI2024

Circumventing Kinetic Barriers to Metal Hydride Formation with Metal–Ligand Cooperativity

We report the two-electron, one-proton mechanism of cobalt hydride formation for the conversion of [Co III Cp(P Ph 2 N Bn 2 )(CH 3 CN)] 2+ to [HCo III Cp(P Ph 2 NBn 2 )] + . This complex catalytically converts CO 2 to formate under CO 2 reduction conditions, with hydride formation as a key elementary step. Through a combination of electrochemical measurements, digital simulations, theoretical calculations, and additional mechanistic and thermochemical studies, we outline the explicit role of the P Ph 2 N Bn 2 ligand in the proton-coupled electron transfer (PCET) reactivity that leads to hydride formation. We reveal three unique PCET mechanisms, and we show that the amine on the P Ph 2 N Bn 2 ligand serves as a kinetically accessible protonation site en route to the thermodynamically favored cobalt hydride. Cyclic voltammograms recorded with proton sources that span a wide range of pK a values show four distinct regimes where the mechanism changes as a function of acid strength, acid concentration, and timescale between electrochemical steps. Peak shift analysis was used to determine proton transfer rate constants where applicable. Furthermore, this work highlights the astute choices that must be made when designing catalytic systems, including the basicity and kinetic accessibility of protonation sites, acid strength, acid concentration, and timescale between electron transfer steps, to maximize catalyst stability and efficiency.

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

Mechanistic Insights into the Catalytic Condensation of Methyl Ketones on MgO Surfaces

Ketone coupling via aldol condensation is one of the promising routes to produce cyclic and value-added precursors for renewable hydrocarbon biofuels. A first-principles-based microkinetic modeling is performed to evaluate the surface-mediated reaction mechanisms and the role of water molecules in the observed activities for 2-pentanone and 3-pentanone aldol condensation on dehydroxylated MgO(111) surface and hydroxylated terminated surface[OH-MgO(111)]. Here, we have identified the enhancement of the surface OH group to MgO(111) surface catalytic activity by destabilizing the binding strength of reaction intermediates and reducing the energy barriers of rate-determining steps(proton transfer and dehydration steps). The 2-pentanone has one elementary step less in the complete reaction mechanism of aldol condensation and preferable energy barrier for proton transfer and dehydration steps, revealing 2-pentanone as terminal ketone is more reactive than 3-pentanone as central ketone. The water molecules dominated the OH-MgO( 111) surface after further addition of water, leading to the reduction of turnover frequency of the aldol condensation dimer product as the loss of aldol condensation reaction intermediates in competitive adsorption with water molecules.

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

From Coherence to Function: Exploring the Connection in Chemical Systems

The role of quantum mechanical coherences or coherent superposition states in excited state processes has received considerable attention in the last two decades largely due to advancements in ultrafast laser spectroscopy. These coherence effects hold promise for enhancing the efficiency and robustness of functionally relevant processes, even when confronted with energy disorder and environmental fluctuations. Understanding coherence deeply drives us to unravel mechanisms and dynamics controlled by order and synchronization at a quantum mechanical level, envisioning optical control of coherence to enhance functions or create new ones in molecular and material systems. In this frontier, the interplay between electronic and vibrational dynamics, specifically the influence of vibrations in directing electronic dynamics, has emerged as the leading principle. Here, two energetically disparate quantum degrees of freedom work in-sync to dictate the trajectory of an excited state reaction. Moreover, with the vibrational degree being directly related to the structural composition of molecular or material systems, new molecular designs could be inspired by tailoring certain structural elements. In the realm of chemical kinetics, our understanding of the dynamics of chemical transformations is underpinned by fundamental theories, such as transition state theory, activated rate theory, and Marcus theory. These theories elucidate reaction rates by considering the energy barriers that must be overcome for reactants to transform into products. Those barriers are surmounted by the stochastic nature of energy gap fluctuations within reacting systems, emphasizing that the reaction coordinate, the pathway from reactants to products, is not rigidly defined by a specific vibrational motion but encompasses a diverse array of molecular motions. While less is known about the involvement of specific intramolecular vibrational modes, their significance in certain cases cannot be overlooked. In this Account, we summarize key experimental findings that offer deeper insights into the complex electronic–vibrational trajectories encompassing excited states afforded from state-of-the-art ultrafast laser spectroscopy in three exemplary processes: photoinduced electron transfer, singlet–triplet intersystem crossing, and intramolecular vibrational energy flow in molecular systems. We delve into the rapid decoherence, or loss of phase and amplitude correlations, of vibrational coherences along promoter vibrations during subpicosecond intersystem crossing dynamics in a series of binuclear platinum complexes. This rapid decoherence illustrates the vibration-driven reactive pathways from the Franck–Condon state to the curve crossing region. We also explore the generation of new vibrational coherences induced by impulsive reaction dynamics rather than by the laser pulse in these systems, which sheds light on specific energy dissipation pathways and thereby on the progression of the reaction trajectory in the vicinity of the curve crossing on the product side. Another property of vibrational coherences, amplitude, reveals how energy can flow from one vibration to another in the electronic excited state of a terpyridine–molybdenum complex hosting a nonreactive dinitrogen substrate. In conclusion, a slight change in vibrational energy triggers a quasi-resonant interaction, leading to constructive wavepacket interference and ultimately intramolecular vibrational redistribution from a Franck–Condon active terpyridine vibration to a dinitrogen stretching vibration, energizing the dinitrogen bond.

Electrical energy↗
DOE OSTI2021

Towards Efficient Uncertainty estimation in deep learning for robust energy prediction in crystal materials

In material science, recent studies have started to explore the potential of using deep learning to improve property prediction from high-fidelity simulations, e.g, density functional theory (DFT). However, the design spaces are sometimes too large and intractable to sample completely. This results in a critical question that is how to evaluate the confidence and robustness of the prediction. In this paper, we propose an efficient approach to estimate uncertainty in deep learning using a single forward pass and then apply it for robust prediction of the total energy in crystal lattice structures. Our approach is built upon the deep kernel learning (DKL) that originally introduces to leverage the expressiveness of deep neural networks as input with a probabilistic prediction of Gaussian processes (GPs) as output. Existing DKL methods have difficulties in the accuracy of predictive uncertainty, training stability, and scaling to large datasets, which lead to significant barriers in real-world applications. We propose to address these challenges by using an inducing point approximate GP in feature space combined with spectral normalization as a regularization. We finally demonstrate our robust performance on an artificial example and a real-world application from materials chemistry.

Bi, Sirui↗
DOE OSTI2024

Self-compensation of group-V acceptors in CdTe

Cadmium Telluride is at the core of an important thin-film technology for photovoltaics that is already commercially available, yet the CdTe-based solar cell efficiency remains limited at 22%, well below the theoretical limit of ~30%. Increasing the hole concentration is crucial for higher efficiency, and group-V elements such as As, P, and Sb are potential dopants as they are shallow acceptors. Nevertheless, group-V doped p-type CdTe often exhibits low doping activation, and the compensation source remains debated. Here, we performed hybrid density functional calculations with spin-orbit coupling to investigate possible sources of hole compensation in group-V doped CdTe. First, regarding possible self-compensation of the group-V dopants, we find that the formation of AX centers is unlikely since they are found to be unstable relative to the shallow acceptor forms. However, if the group-V dopants come in during growth (such as dimer molecules As2, P2, and Sb2), we find that the impurity atoms, which would occupy nearest neighbor sites, maintain the V-V bonds, limiting the hole density. For the native defects, our study reveals that Cd interstitial is the lowest energy donor defect in p-type CdTe. Still, it has a small migration barrier of 0.5 eV, making it highly mobile and unstable at room temperature. The Te vacancy is the next lowest formation energy donor. The migration barrier of 1.4 eV indicates that the Te vacancy is stable at room temperature. The antisite CdTe is also a donor, with low formation energy and stable at room temperature, potentially limiting the hole concentration. Our results, therefore, shed light on possible compensation centers and some guidance on how to avoid them.

14 SOLAR ENERGY↗
DOE OSTI2021

Valuation of Surface Coatings in High-Energy Density Lithium-ion Battery Cathode Materials

Artificial barriers, usually with either electrochemically active or inactive coating materials, are deployed on cathode material surfaces to mitigate detrimental side reactions by suppressing direct contact of cathode and electrolyte called surface coatings. These surface coatings are commonly known to increase the wettability of liquid electrolyte and reduce the interfacial charge transfer resistance. An important caveat is the selection of appropriate coating material with appropriate thickness for achieving enhanced electrochemical performance. As modern battery materials are increasingly developed with some type of surface coating, a careful and thorough examination of their role in mitigating the cycle life issues of cathode materials is paramount. Here, this comprehensive review article extensively covers the selection criteria of coating materials based on their chemical and physical properties and electrochemical functionalities. Additionally, the article discusses the concept of critical coating thickness and methods of achieving homogeneous coating architectures that deliver desired performance benefits. Furthermore, this comprehensive article summarizes the recent advancements, effectiveness, necessity of cathode surface coatings and identifies the key aspect of structure-property correlation between coating type/thickness and lithium-ion diffusion through coating layers as the linchpin that validates surface coating approaches especially for high capacity nickel-rich cathodes.

25 ENERGY STORAGE↗
DOE OSTI2022

Structure, Stability, and Electronic Properties of Boron Suboxide: A Density Functional Theory Study

Boron suboxide (B 6 O) is a boron-rich compound derived from the α-rhombohedral boron lattice with extreme hardness and unusual semiconducting properties. In this work, density functional theory (DFT) was used to show that unit cell volume, mechanical strength, band gaps, and thermodynamic stabilities of B 6 O were influenced by the interstitial elements and point defects at the icosahedral sites. While the hexagonal unit cell volume (HUCV) varies with interstitial occupancy, it is the icosahedral defect that weakens the intrinsic bulk modulus of B 6 O. Using the hybrid HSE functional, we confirmed that the perfect B 6 O bulk is a p-type semiconductor with a direct band gap of 2.8 eV. Furthermore, by screening α-boron compounds systematically, we found that a simple octet rule may offer a consistent explanation for the variations in the computed electronic structures. The formation free energies calculated over a wide range of temperatures (0–2500 K) and pressures (0–80 GPa) predict that formations of interstitial defects become favorable only at higher temperatures (ca. 1800 K) in bulk B 6 O lattices. The nudged elastic band (NEB) method was employed to identify the minimum energy pathways for the diffusions of dislocated B and O atoms. The diffusion of icosahedral B atoms has an energy barrier of 0.16 eV. More complex B diffusion paths involving the reorganization of icosahedral boron atoms incur higher barriers (>1 eV). In contrast, the diffusion of interstitial O atoms is facile with a barrier of 0.4 eV. Lastly, successive O insertions into the α-B lattice were performed using DFT to generate a basic understanding of the oxidation process. Furthermore, these calculations provide fundamental atomistic insights into the growth of B 6 O crystals and control of their point defects.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Acceleration of Pd-V intermetallic diffusion by hydrogen

Vanadium-based membranes have great potential for hydrogen purification due to their perfect selectivity, high permeability, and relatively low cost. With appropriate surface cleaning, V efficiently permeates hydrogen at elevated temperature, but performance declines due to its affinity to absorb impurities. Here, the application of palladium thin films maintains a clean surface that catalyzes hydrogen dissociation and recombination. Hydrogen permeation in Pd-V-Pd membranes initially reach theoretical permeability, but declines due to Pd-V interdiffusion. The objective of this work was to quantify the intermetallic diffusion process as a function of temperature and ambient. Pd-V composites were subjected to various annealing treatments and characterized using Auger electron spectroscopy, X-ray diffraction, and energy dispersive X-ray spectroscopy, as well as correlated to measurements of membrane permeability. In an inert environment Pd-V interdiffusion was observable as low as 300 °C, and the diffusion coefficient had an activation energy of 44 kJ/mol. Furthermore, the presence of hydrogen at partial pressures > 10 kPa accelerated interdiffusion six-fold at T = 400 °C. Membrane performance degraded with an activation energy 75 kJ/mol, suggesting that intermetallic diffusion leads to both a loss of catalytic activity and as well as degradation of bulk permeability. These findings provide a baseline for evaluating hydrogen permeable interdiffusion barriers to overcome these challenges.

08 HYDROGEN↗
DOE OSTI2022

Energy-Efficient Microbial Electrochemical Lignin and Alkaline Hydroxide Recovery from DMR Black Liquor

Black liquor is a complex waste stream generated from the fast-expanding biorefinery industry, and the high management cost has been known as a major development barrier. In this work, we demonstrated a rationally designed 4-chamber microbial electrolysis process that was capable of recovering the value-added lignin, NaOH, and H2 products while removing waste organics from the Deacetylation and Mechanical Refining (DMR) black liquor, which could potentially reduce the disposal cost by up to 40%. The lab-scale MEC achieved 52% organic removal from a highly concentrated feedstock (28,960 ± 212 mg/L COD) and recovered 83.2 mM NaOH under a small voltage application (1.2 V). Lignin recovery was observed to be as high as 13.98 +/- 0.28 g/L, and 0.35 L/L/day H2 was produced with a faradaic efficiency of 93.5%. This work provided a good example on how complex liquid waste stream can be valorized via multi-functional microbial electrochemical technology, and how the traditionally unfavored pH gradient in an electrochemical reactor can be manipulated for lignin precipitation and alkaline recovery at the same time. The versatility and good performance in concurrent contaminants removal, lignin precipitation, H2 production and NaOH recovery demonstrate a good potential for microbial electrolysis in black liquor treatment.

BIOMASS FUELS,INORGANIC, ORGANIC, PHYSICAL, AND AN↗
DOE OSTI2020

Subsurface Nitrogen Dissociation Kinetics in Lithium Metal from Metadynamics

The dissociation of molecular nitrogen in lithium is of interest for several promising technologies, such as the catalytic synthesis of ammonia in ambient or mild conditions. In this work we simulate nitrogen dissociation in the lithium BCC (110) surface at ambient and elevated temperatures using density functional theory (DFT) metadynamics simulations. The rate constants at temperatures of 300, 400, and 500 K are calculated by statistical analysis of the reaction time distributions from the accelerated simulations. This approach finds and estimates rate constants for transition pathways out of the initial state; the required input is the stable initial state and a reasonable choice of collective variable. A single collective variable is used in this case: the N–N distance. The results are robust to changes in metadynamics parameters, and the reaction time distributions follow the expected exponential distribution. We show that the metadynamics-derived rate constants are in agreement with results from the conventional harmonic approximation approach using a climbing image nudged elastic band (NEB) transition state search. The reaction barriers from metadynamics and the NEB/harmonic approximation agree to within 0.02–0.04 eV at all temperatures studied. This paper demonstrates that the harmonic approximation provides an accurate description of the rate constants for nitrogen dissociation in lithium metal, even at temperatures near or above the melting point of lithium, lending credence to previous and future theoretical studies using this approximation. Moreover, this work demonstrates a step toward the automated exploration and discovery of reaction mechanisms and associated rate constants for elementary surface-catalyzed reactions using DFT-based metadynamics.

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

Accurate Fe–He machine learning potential for studying He effects in BCC-Fe

Nanostructured Ferritic Alloys (NFAs) containing oxide nanoparticles are prospective advanced structural materials in future fusion reactors. Lacking fusion neutron sources with an adequate flux, modeling is critical to understand radiation damage in the NFAs, including helium effects. Machine learning interatomic potentials (MLPs) have been demonstrated to be capable of describing atomic interactions in chemically complex systems with comparable accuracy to density functional theory (DFT) but with much less computational costs. Hence, MLPs are chosen to study He bubble accumulation in chemically complex NFAs. As a preliminary step, we develop a Fe–He potential, based on ~10,000 atomic configurations. The developed MLP accurately predicts the bulk properties of BCC-Fe compared to DFT, including the lattice constant (2.832 Å), elastic constants (c 11 = 271 GPa, c 12 = 141, c 44 = 93), and phonon frequencies (maximum error < 3.5%). The MLP predicts He binding energies in He n V bubbles and He n clusters with a mean absolute error (MAE) of only ~70 meV, which is ~3 to ~5 times smaller than the MAE of binding energies estimated using empirical potentials. The MLP also accurately predicts the migration barrier of interstitial He (64 meV). Furthermore, the MLP correctly predicts the maximum number of He atoms (critical size) at which a He n V bubble (He n cluster) can kick-out a self-interstitial atom (SIA), forming a bigger He n V 2 bubble (He n V bubble). Subsequently, we use the MLP to perform molecular dynamics simulations to investigate the effect of temperature on the critical size of He n V bubbles and He n clusters. We find that the critical size decreases with increasing temperature — indicating that smaller bubbles and clusters can kick-out an SIA at higher temperatures, thereby escalating the He bubble growth. We report the current Fe–He MLP can be further developed to include all chemical interactions in NFAs.

36 MATERIALS SCIENCE↗