Search NASA⌕ Search

SEARCH · Search NASA

Results for “materials and molecular modeling”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Model for Humidity-Mediated Diffusion on Aluminum Surfaces and Its Role in Accelerating Atmospheric Aluminum Corrosion

Bare aluminum metal surfaces are highly reactive, which leads to the spontaneous formation of a protective oxide surface layer. Because many subsequent corrosive processes are mediated by water, the structure and dynamics of water at the oxide interface are anticipated to influence corrosion kinetics. Here, using molecular dynamics simulations with a reactive force field, we model the behavior of aqueous aluminum metal ions in water adsorbed onto aluminum oxide surfaces across a range of ion concentrations and water film thicknesses corresponding to increasing relative humidity. We find that the structure and diffusivity of both the water and the metal ions depend strongly on the humidity of the environment and the relative height within the adsorbed water film. Aqueous aluminum ion diffusion rates in water films corresponding to a typical indoor relative humidity of 30% are found to be more than 2 orders of magnitude slower than self-diffusion of water in the bulk limit. Connections between metal ion diffusivity and corrosion reaction kinetics are assessed parametrically with a reductionist model based on a 1D continuum reaction–diffusion equation. Our results highlight the importance of incorporating the properties specific to interfacial water in predictive models of aluminum corrosion.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Porous Covalent Organic Polymers for Efficient Fluorocarbon‐Based Adsorption Cooling

Abstract Adsorption‐based cooling is an energy‐efficient renewable‐energy technology that can be driven using low‐grade industrial waste heat and/or solar heat. Here, we report the first exploration of fluorocarbon adsorption using porous covalent organic polymers (COPs) for this cooling application. High fluorocarbon R134a equilibrium capacities and unique overall linear‐shaped isotherms are revealed for the materials, namely COP‐2 and COP‐3. The key role of mesoporous defects on this unusual adsorption behavior was demonstrated by molecular simulations based on atomistic defect‐containing models built for both porous COPs. Analysis of simulated R134a adsorption isotherms for various defect‐containing atomistic models of the COPs shows a direct correlation between higher fluorocarbon adsorption capacities and increasing pore volumes induced by defects. Combined with their high porosities, excellent reversibility, fast kinetics, and large operating window, these defect‐containing porous COPs are promising for adsorption‐based cooling applications.

Zheng, Jian↗
DOE OSTI2020

Hardening in Au-Ag nanoboxes from stacking fault-dislocation interactions

Porous, nano-architected metals with dimensions down to ~10 nm are predicted to have extraordinarily high strength and stiffness per weight, but have been challenging to fabricate and test experimentally. Here, we use colloidal synthesis to make ~140 nm length and ~15 nm wall thickness hollow Au-Ag nanoboxes with smooth and rough surfaces. In situ scanning electron microscope and transmission electron microscope testing of the smooth and rough nanoboxes show them to yield at 130 ± 45 MPa and 96 ± 31 MPa respectively, with significant strain hardening. A higher strain hardening rate is seen in rough nanoboxes than smooth nanoboxes. Finite element modeling is used to show that the structure of the nanoboxes is not responsible for the hardening behavior suggesting that material mechanisms are the source of observed hardening. Molecular dynamics simulations indicate that hardening is a result of interactions between dislocations and the associated increase in dislocation density.

42 ENGINEERING↗
NASA NTRS1996

Organics In Meteorites

The variety of classes of organic compounds that occur in carbonaceous meteorites suggests a rich pre-planetary chemistry with possible connections to interstellar, solar nebular and parent body processes. Structural diversity prevails within all classes examined in detail. Among amino acids for instance, all possible isomers are found up to species containing 4-6 carbon atoms, with abundances decreasing with increasing molecular weight. Such diversity seems limited to those carbonaceous meteorites which show evidence of having been exposed to liquid water; meteorites lacking such evidence also show much lower abundances and less structural diversity in their organic contents. This apparent dependency on water suggests a role for cometary ices in the chemical evolution of organic compounds on parent bodies. Measurements of the stable isotope compositions of C, H, N and S in classes of compounds and at the individual compound level show strong deviations from average chondritic values. These deviations are difficult to explain by solar system or parent body processes, and precedents for some of these isotopic anomalies exist in interstellar (e.g., high D/H ratios) and circumstellar chemistry. Therefore, presolar origins for much if not all of the meteoritic organic compounds (or their precursors) is a distinct possibility. In contrast, evidence of solar nebular origins is either lacking or suspect. Results from molecular and isotopic analyses of meteoritic organics, from laboratory simulations and from a model of interstellar grain reactions will be used to flesh out the hypothesis that this material originated with interstellar chemistry, was distributed within the early solar system as cometary ices, and was subsequently altered on meteorite parent bodies to yield the observed compounds.

Chang, Sherwood↗
DOE OSTI2020

Electrokinetic Analysis of Poorly Conductive Electrocatalytic Materials

While many factors can influence the activity of the electrochemical interface at atomic and molecular levels, the electron and ion transport have a big impact on the overall electrochemical performance. Here, we investigate how electrical resistance (conductivity) of the electrode affects the measured electrocatalytic activity using the hydrogen evolution reaction as a model reaction. We clearly demonstrate the effect of the poor electronic conductivity of MoS 2 on the electrochemical performance, as both values of the exchange current density as well as the Tafel slope depend on the level of utilization of the active material and the ohmic losses. A quantitative analysis using the composite electrode model is provided, showing the importance of the experimental conditions and good control of electrode/electrolyte system for the proper evaluation and comparison of electrocatalysts’ kinetic parameters.

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

Quantifying multipoint ordering in alloys

A central problem in multicomponent lattice systems is to systematically quantify multipoint ordering. Ordering in such systems is often described in terms of pairs, even though this is not sufficient when three-point and higher-order interactions are included in the Hamiltonian. Current models and parameters for multipoint ordering are often only applicable for very specific cases or require approximating a subset of correlated occupational variables on a lattice as being uncorrelated. In this paper, cluster order parameters are introduced to systematically quantify arbitrary multipoint ordering motifs in substitutional systems through direct calculations of normalized cluster probabilities. These parameters can describe multipoint chemical ordering in crystal systems with multiple sublattices, multiple components, and systems with reduced symmetry. These are defined in this paper and applied to quantify four-point chemical ordering motifs in platinum/palladium alloy nanoparticles that are of practical interest to the synthesis of catalytic nanocages. Impacts of chemical ordering on nanocage stability are discussed. It is demonstrated that approximating four-point probabilities from superpositions of lower-order pair probabilities is not sufficient in cases where three- and four-body terms are included in the energy expression. Furthermore, conclusions about the formation mechanisms of nanocages may change significantly when using common pair approximations.

36 MATERIALS SCIENCE↗
NASA NTRS2019

Recent Advancements in Modeling and Simulation of Entry Systems at NASA

This paper describes recent development of modeling and simulation technologies for entry systems in support of NASA's exploration missions. Mission-tailored research and development in modeling of entry systems occurs across the Agency (e.g., within the Orion and Mars 2020 Programs), however the aim of this paper is to discuss the broad, cross-mission research conducted by NASA's Entry Systems Modeling (ESM) Project, which serves as the Agency's only concerted effort toward advancing entry systems across a range of technical disciplines. Technology development in ESM is organized and prioritized from a system-level perspective, resulting in four broad technical areas of investment: (1) Predictive material modeling, (2) Shock layer kinetics and radiation, (3) Computational and experimental aerosciences, and (4) Guidance, navigation, and control. Investments in thermal protection material modeling are geared toward high-fidelity, predictive models capable of handling complex structures, with an eye toward optimizing design performance and quantifying thermal protection system reliability. New computational tools have been developed to characterize material properties and behavior at the microstructural level, and experimental techniques (molecular beam scattering, micro-computed tomography, among others) have been developed to measure material kinetics, morphology, and other parameters needed to inform and validate detailed simulations. Advancements have also been made in macrostructural simulation capability to enable 3-D system-scale calculations of material response with complex topological features, including differential recession of tile gaps. Research and development in the area of shock layer kinetics has focused on air and CO2-based atmospheres. Capacity and capability of the NASA Ames Electric Arc Shock Tube (EAST) have been expanded in recent years and analysis of resulting data has led to several improvements in kinetic models, while simultaneously reducing uncertainties associated with radiative heat transfer predictions. First-principles calculations of fundamental kinetic, thermodynamic, and transport data, along with state-specific models for non-equilibrium flow regimes, have also yielded new insights and have the potential to vastly improve model fidelity. Aerosciences is a very broad area of interest in entry systems, yet a number of important challenges are being addressed: Coupled fluid-structure simulations of parachute inflation and dynamics; Experimental and computational studies of vehicle dynamics; Multi-phase flow with dust particles to simulate entry environments at Mars during dust storms; Studies of roughness-induced heating augmentation relevant to tiled and woven thermal protection systems; and Advanced numerical methods to optimize computational analyses for desired accuracy versus cost. Guidance and control in the context of entry systems has focused on development of methods for multi-axis control (i.e. pitch and yaw, rather than bank angle alone) of spacecraft during entry and descent. With precision landing requirements driven by Mars human exploration goals, recent efforts have yielded 6-DOF models of multi-axis control with propulsive descent of both inflatable and rigid ellipsled-like architectures.

Barnhardt, Michael D.↗
NASA NTRS2021

Recent Advancements in Modeling and Simulation of Entry Systems at NASA

This paper describes recent development of modeling and simulation technologies for entry systems in support of NASA’s exploration missions. Mission-tailored research and development in modeling of entry systems occurs across the Agency (e.g., within the Orion and Mars 2020 Programs), however the aim of this paper is to discuss the broad, cross-mission research conducted by NASA’s Entry Systems Modeling (ESM) Project, which serves as the Agency’s only concerted effort toward advancing entry systems across a range of technical disciplines. Technology development in ESM is organized and prioritized from a system-level perspective, resulting in four broad technical areas of investment: (1) Predictive material modeling, (2) Shock layer kinetics and radiation, (3) Computational and experimental aerosciences, and (4) Guidance, navigation, and control. Investments in thermal protection material modeling are geared toward high-fidelity, predictive models capable of handling complex structures, with an eye toward optimizing design performance and quantifying thermal protection system reliability. New computational tools have been developed to characterize material properties and behavior at the microstructural level, and experimental techniques (molecular beam scattering, micro-computed tomography, among others) have been developed to measure material kinetics, morphology, and other parameters needed to inform and validate detailed simulations. Advancements have also been made in macrostructural simulation capability to enable 3-D system-scale calculations of material response with complex topological features, including differential recession of tile gaps. Research and development in the area of shock layer kinetics has focused on air and CO2-based atmospheres. Capacity and capability of the NASA Ames Electric Arc Shock Tube (EAST) have been expanded in recent years and analysis of resulting data has led to several improvements in kinetic models, while simultaneously reducing uncertainties associated with radiative heat transfer predictions. First-principles calculations of fundamental kinetic, thermodynamic, and transport data, along with state-specificmodels for non-equilibrium flow regimes, have also yielded new insights and have the potential to vastly improve model fidelity. Aerosciences is a very broad area of interest in entry systems, yet a number of important challenges are being addressed: Coupled fluid-structure simulations of parachute inflation and dynamics; Experimental and computational studies of vehicle dynamics; Multi-phase flow with dust particles to simulate entry environments at Mars during dust storms; Studies of roughness-induced heating augmentation relevant to tiled and woven thermal protection systems; and Advanced numerical methods to optimize computational analyses for desired accuracy versus cost. Guidance and control in the context of entry systems has focused on development of methods for multi-axis control (i.e. pitch and yaw, rather than bank angle alone) of spacecraft during entry and descent. With precision landing requirements driven by Mars human exploration goals, recent efforts have yielded 6-DOF models of multi-axis control with propulsive descent of both inflatable and rigid ellipsled-like architectures. Results for both configurations have demonstrated the ability to land within the 50-meter precision requirement demanded by Mars human exploration missions, while also reducing propellant requirements by enabling more efficient control through entry and descent. Ongoing research in GN&C is developing mechanical specifications for the systems and establishing engineering feasibility.

aerothermodynamics↗
DOE OSTI2021

Modeling of Peptides with Classical and Novel Machine Learning Force Fields: A Comparison

The replacement of classical force fields (FFs) with novel neural-network-based frameworks is an emergent topic in molecular dynamics (MD) simulations. In contrast to classical FFs, which have proven their capability to provide insights into complex soft matter systems at an atomistic resolution, the machine learning (ML) potentials have yet to demonstrate their applicability for soft materials. However, the underlying philosophy, which is learning the energy of an atom in its surrounding chemical environment, makes this approach a promising tool. In particular for the exploration of novel chemical compounds, which have not been considered in the original parametrization of classical FFs. In this article, we study the performance of the ANI-2x ML model and compare the results with those of two classical FFs, namely, CHARMM27 and the GROMOS96 43a1 FF. We explore the performance of these FFs for bulk water and two model peptides, trialanine and a 9-mer of the α-aminoisobutyric acid, in vacuum and water. The results for water describe a highly ordered water structure, with a structure similar to those using ab initio molecular dynamics simulations. The energy landscape of the peptides described by Ramachandran maps show secondary structure basins similar to those of the classical FFs but differ in the position and relative stability of the basins. Details of the sampled structures show a divergent performance of the different models, which can be related either to the short-ranged nature of the ML potentials or to shortcomings of the underlying data set used for training. These findings highlight the current state of the applicability of ANI-2x ML potential for MD simulations of soft matter systems. Simultaneously, they provide insights for future improvements of current ML potentials.

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

Molecular dynamics study on interface formation and bond strength of impact-welded Mg-steel joints

It was recently demonstrated that the vaporizing foil actuator welding (VFAW) method can directly join immiscible magnesium and steel alloys without coating or a third chemical element based intermetallic compound layer. The VFAW Mg/steel joint exhibits a mixed interface layer of up to 200μm thickness consisting of Mg matrix and Fe particles. Computer simulations have suggested the formation of the interlayer is from the high-velocity frictional shearing between the Mg/steel substrates during the oblique impact in the VFAW process. This paper investigates the formation of Mg-Fe interlayer under VFAW condition with different shearing velocities using molecular dynamics (MD) model, and studies the bonding strength under different scenarios. Finally, the results elucidate the critical role of shearing velocity and surface roughness in achieving Mg/Fe joint.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Applying machine learning and quantum chemistry to predict the glass transition temperatures of polymers

Glass transition temperature (T g ) is important for understanding the physical and mechanical properties of a polymer material because it relates to the thermal energy required to transition between a hard glassy state and a soft rubbery one. Over the years, various models have been developed for predicting this thermal property from molecular structure to aid in designing novel polymers in selected classes. This work builds on those efforts by utilizing both machine learning (ML) and quantum chemistry (QC) techniques to develop models that can predict T g values from the molecular structure under different data availability scenarios and for a wide variety of polymer types. For the ML model, a graph convolutional network (GCN) was used to map topological polymer features; this model was trained against a dataset of more than 7500 T g values and resulted in a root mean square error (RMSE) of 38.1 °C. The QC-based regression model was trained on 83 T g values and produced an RMSE of 34.5 °C. In conclusion, this work demonstrated that while both model techniques produce accurate predictions and are suitable for different data availability scenarios, the QC-based regression model offered a more interpretable model framework with significantly less training data.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Controlling morphology in hybrid isotropic/patchy particle assemblies

Brownian dynamics is used to study self-assembly in a hybrid system of isotropic particles (IPs), combined with anisotropic building blocks that represent special “designer particles.” Those are modeled as spherical patchy particles (PPs) with binding only allowed between their patches and IPs. In this study, two types of PPs are considered: Octahedral PPs (Oh-PPs) and Square PPs (Sq-PPs), with octahedral and square arrangements of patches, respectively. The self-assembly is additionally facilitated by the simulated annealing procedure. Here, the resultant structures are characterized by a combination of local correlations in cubatic ordering and a symmetry-specific variation of bond orientation order parameters (SymBOPs). By varying the PP/IP size ratio, we detected a sharp crossover between two distinct morphologies in both types of systems. High symmetry phases, NaCl crystal for Oh-PP and square lattice for Sq-PP, are observed for larger size ratios. For the smaller ones, the dominant morphologies are significantly different, e.g., Oh-PPs form a compact amorphous structure with predominantly face-to-face orientation of neighboring PPs. Unusually, for a morphology without a long-range order, it is still possible to identify well organized coherent clusters of this structure, thanks to the adoption of our SymBOP-based characterization.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2026

Molecular dynamics simulations of reflection and sputtering behavior of boron under deuterium ion irradiation

Boronization is a commonly used method of wall conditioning in fusion reactors. The application of boron films to the plasma-facing materials results in enhanced plasma performance due to the reduction of intrinsic impurities. This is primarily driven by a reduction in oxygen content that is chemically trapped in the boron film. The reactive nature of these boron films also raises questions concerning interactions with hydrogen isotopes. In this work, boron-deuterium interactions were studied using molecular dynamics (MD). Reactive force field potentials were used to model the chemical interactions between B and D. An amorphous boron substrate was irradiated by D atoms at varying incident energies, 10 eV < E i < 150 eV and angles, 0° < α < 85°. The reflection probability was calculated and compared to results from the commonly-used binary collision approximation (BCA) method. This comparison found that the BCA underestimated the reflection probability at E i < 35 eV and α > 45 ∘ . The source of this discrepancy was found to be the surface binding energy model. The BCA calculation with an isotropic surface binding energy model was more closely aligned to the MD result. This, in combination with a correction function based on the MD results allows for corrections to the reflection probability of deuterium impinging on boron surfaces. The sputtering of the substrate material was also studied. While this study did not contain sufficient events to quantitatively describe the sputtering behavior, some qualitative results emerged: namely, chemical sputtering of B and D-containing molecules (BD, BD 2 , BD 3 ) at low ( < 20 eV) incident deuterium energies. This result suggests that chemical sputtering could be a significant factor in limiting boron coating lifetime when exposed to lower ion energies, such as those in detached plasmas. The results show that chemical interactions should be taken into account when modeling the interactions between D ions and B surfaces.

Boron↗
DOE OSTI2021

Quantum Machine-Learning for Eigenstate Filtration in Two-Dimensional Materials

Quantum machine-learning algorithms have emerged to be a promising alternative to their classical counterparts as they leverage the power of quantum computers. Such algorithms have been developed to solve problems like electronic structure calculations of molecular systems and spin models in magnetic systems. However, the discussion in all these recipes focuses specifically on targeting the ground state. Herein we demonstrate a quantum algorithm that can filter any energy eigenstate of the system based on either symmetry properties or a predefined choice of the user. The workhorse of our technique is a shallow neural network encoding the desired state of the system with the amplitude computed by sampling the Gibbs–Boltzmann distribution using a quantum circuit and the phase information obtained classically from the nonlinear activation of a separate set of neurons. We show that the resource requirements of our algorithm are strictly quadratic. To demonstrate its efficacy, we use state filtration in monolayer transition metal dichalcogenides which are hitherto unexplored in any flavor of quantum simulations. We implement our algorithm not only on quantum simulators but also on actual IBM-Q quantum devices and show good agreement with the results procured from conventional electronic structure calculations. We thus expect our protocol to provide a new alternative in exploring the band structures of exquisite materials to usual electronic structure methods or machine-learning techniques that are implementable solely on a classical computer.

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

Capacitance of Carbon Nanotube/Graphene Composite Electrodes with [BMIM + ][BF$_4^–$]/Acetonitrile: Fixed Voltage Molecular Dynamics Simulations

The capacitance of nanoporous carbon electrode materials is dictated by the physical interactions with electrolyte ions and molecules at the accessible interior electrode surface. While significant progress has been made in designing and synthesizing carbon materials with well-defined and relatively homogeneous nanoporosity, the majority of materials remain heterogeneous, and the resulting properties reflect an average over this distribution. In this regard, computer simulations can be a valuable tool to predict or validate structure/property relationships by systematically investigating well-defined, model electrode morphologies. In this work, we utilize fixed-voltage molecular dynamics simulations to predict structure/capacitance relationships for five different model morphologies of carbon nanotube/graphene (CNT/G) composite electrodes with 1-butyl-3-methylimidazolium tetrafluoroborate/acetonitrile electrolyte. The CNT/G electrode models are inspired by experimental “layer-by-layer” syntheses and strike a balance between realism and computational tractability. Comparison between different model electrode architectures elucidates important structure/property relationships. We find that CNT/graphene contact points serve as “hot spots” with significantly enhanced charge separation relative to the rest of the electrode. Furthermore, we demonstrate a specific nanoconfinement motif that provides substantial 3–4$\times$ enhancement of local capacitance, resulting in a ~40% increase of the total electrode differential capacitance. Because the accessible surface area of the model CNT/G electrodes is precisely determined, a comparison of per-area capacitance across systems is unambiguous. Our results thus complement a prior computational demonstration of capacitance enhancement in nanoconfinement while elucidating additional interaction motifs at CNT/G electrochemical interfaces.

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

High Energy Density Physics of Inertial Confinement Fusion Ablator Materials

The historic December 5, 2022 experiment at Lawrence Livermore National Lab’s (LLNL) National Ignition Facility (NIF) reached fusion energy ignition for the first time. This is the most important scientific breakthrough of the 21st century paves the way to future clean inertial fusion energy (IFE). The diamond (high density carbon (HDC)) ablator material used in this experiment displays detrimental effects due to the development of hydrodynamic instabilities at the diamond/fuel interface under shock compression. New alternatives to diamond ablators are required to step up the energy yield in ICF experiments. The unique combination of mechanical strength (approaching that of diamond), the ability to accommodate high-Z dopants (in contrast to diamond), and the tunability of the properties (through synthesis material with varying sp 3 content) make amorphous carbon (a-C) a promising material for next-generation IFE ablative capsules. However, despite its critical importance to the IFE program, the behavior of a-C carbon at extreme temperatures and pressures remains largely unexplored. The primary goals of this project were to perform groundbreaking dynamic compression experiments and predictive simulations to uncover the fundamental high-energy-density physics of amorphous carbon. Our goals were (1) to uncover the metastability range of amorphous carbon and probe phase transitions to diamond or metastable supercooled liquid carbon; (2) to acquire high-quality equation of state (EOS) data and develop an experimentally validated EOS from machine-learning MD simulations of the complex states of carbon; and (3) to uncover the complex behavior of carbon liquid in both thermodynamically stable and metastable supercooled states by accessing large areas of carbon phase diagram with amorphous samples with variable sp 3 content. Our proposed experimental program included measurements of equation of state and diffraction measurements using the Omega EP laser at the Laboratory of Laser Energetics at the University of Rochester. The theoretical/simulation program involved the development of machine-learning models of the complex response of amorphous carbon under dynamic compression by performing molecular dynamics simulations at experimental time and length scales using leadership class DOE supercomputers. Simulations guided experiments to observe predicted phenomena and acquire critical experimental data in specific pressure-temperature domains to validate theoretical models. This research delivered fundamental properties of novel amorphous carbon IFE ablator material, including phase diagram and EOS. These results will aid in IFE target design and implosion experiments. A unique combination of predictive simulations and dynamic and static experiments provided a highly inspirational intellectual environment for graduate students and postdocs involved in this project.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

When more data hurts: Optimizing data coverage while mitigating diversity-induced underfitting in an ultrafast machine-learned potential

Machine-learned interatomic potentials (MLIPs) are becoming an essential tool in materials modeling. However, optimizing the generation of training data used to parametrize the MLIPs remains a significant challenge. This is because MLIPs can fail when encountering local environments too different from those present in the training data. The difficulty of determining a priori the environments that will be encountered during molecular dynamics simulation necessitates diverse, high-quality training data. Here, this study investigates how training data diversity affects the performance of MLIPs using the Ultra-Fast force field (UF 3 ) to model amorphous silicon nitride. We employ expert and autonomously generated data to create the training data and fit four force field variants to subsets of the data. Our findings reveal a critical balance in training data diversity: insufficient diversity hinders generalization, while excessive diversity can exceed the MLIP's learning capacity, reducing simulation accuracy. Specifically, we found that the UF 3 variant trained on a subset of the training data, in which nitrogen-rich structures were removed, offered vastly better prediction and simulation accuracy than any other variant. By comparing these UF 3 variants, we highlight the nuanced requirements for creating accurate MLIPs, emphasizing the importance of application-specific training data to achieve optimal performance in modeling complex material behaviors.

ab initio molecular dynamics↗
DOE OSTI2024

Changes in dislocation punching behavior due to hydrogen-seeded helium bubble growth in tungsten

The accumulation of gas atoms in tungsten is a topic of long-standing interest to the plasma-facing materials community due the metal's use as a divertor material in some tokamak fusion reactors. The nucleation and growth of He/H gas bubbles (along with their isotopes) can result from impinging fluxes of these gases which give rise to damage at the W divertor surface. The inclusion of He or H in W has been studied extensively by the community, finding that He bubbles modify the surface through periodic dislocation punching and bursting mechanisms while H bubbles impact the metal through plastic-strain induced material failure. However, the mechanisms which are present during the combined flux of both He and H is not well-studied atomistically. Motivated by this, an atomistic modeling study is conducted using molecular dynamics to assess the behavior of mixed concentration He:H bubbles in W. Here we find that the introduction of H into a growing He bubble results in a dramatic change in the nature and presence of dislocation loops which are typically generated via dislocation punching in over-pressurized He bubbles. Most notably, at high H concentrations, there is a switchover in energetic favorability from glissile 1/2<111> dislocations to sessile <100> dislocations. This thermodynamic crossover could imply a significant reduction in W surface morphology changes than with pure He bubbles and, additionally, we show this to have implications on the trapping of H in the bubbles and their associated dislocations.

36 MATERIALS SCIENCE↗