Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic simulation”

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 1,459 records · Page 81
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Unveiling mechanisms and onset threshold of humping in high-speed laser welding

The fabrication of fuel cells relies on a rapid laser welding process. However, challenges arise with the occurrence of humping when the welding speed surpasses a critical threshold, which poses difficulties in achieving a smooth surface finish and a consistent weld strength. This study aims to elucidate the humping mechanisms by analyzing the morphology of molten pool and the characteristics of melt flow at varying welding speeds via in situ synchrotron high-speed X-ray imaging and computational fluid dynamics simulations. Our findings indicate that the short keyhole rear wall, the high backward melt velocity, and the prolonged tail of molten pool are the primary factors contributing to the onset of humping. Furthermore, a dimensionless humping index (π h ) was introduced, which successfully captured the onset threshold of humping across different literatures. This index not only provides a quantitative description of the humping formation tendency but also serves as a valuable tool for optimizing the laser welding process.

42 ENGINEERING↗
DOE OSTI2022

On the mechanisms of ion adsorption to aqueous interfaces: air-water vs. oil-water

The adsorption of ions to water-hydrophobe interfaces influences a wide range of phenomena, including chemical reaction rates, ion transport across biological membranes, and electrochemical and many catalytic processes; hence, developing a detailed understanding of the behavior of ions at water-hydrophobe interfaces is of central interest. Here, we characterize the adsorption of the chaotropic thiocyanate anion (SCN − ) to two prototypical liquid hydrophobic surfaces, water-toluene and water-decane, by surface-sensitive nonlinear spectroscopy and compare the results against our previous studies of SCN − adsorption to the air-water interface. For these systems, we observe no spectral shift in the charge transfer to solvent spectrum of SCN − , and the Gibb’s free energies of adsorption for these three different interfaces all agree within error. We employed molecular dynamics simulations to develop a molecular-level understanding of the adsorption mechanism and found that the adsorption for SCN − to both water-toluene and water-decane interfaces is driven by an increase in entropy, with very little enthalpic contribution. This is a qualitatively different mechanism than reported for SCN − adsorption to the air-water and graphene-water interfaces, wherein a favorable enthalpy change was the main driving force, against an unfavorable entropy change.

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

Comparing ab initio and quantum-kinetic approaches to electron transport in warm dense matter

Accurate knowledge of the electronic transport properties of warm dense matter is one of the main concerns of research in high-energy-density physics. Three modern approaches with vastly different levels of fidelity are reviewed and compared: the Kubo–Greenwood (KG) approach based on density-functional-theory molecular dynamics simulations (QMD), quantum kinetic theory based on average-atom models, and time-dependent density functional theory. Throughout, emphasis is placed on the connection between static properties of the electrons (e.g., density of states) and transport properties. Overall, it is found that whenever the conduction electrons can be modeled as being nearly free, fair to excellent agreement is found between QMD and kinetic theory approaches. Such a circumstance is required for modeling warm dense matter as a plasma of ions and free electrons, which is assumed in most kinetic theory approaches. The sensitivity of transport properties to the electronic structure is further highlighted by comparing different exchange–correlation approximations in QMD and KG calculations. It is found that the inclusion of exact exchange via thermal hybrid functionals can make a pronounced impact on electrical and thermal conduction in warm dense matter. We also investigate dynamic screening physics via kinetic theory and time-dependent density functional theory calculations of the mean free path of an electron in a hot dense plasma. In sum, we identify three axes along which to make progress in predicting electron transport in warm dense matter.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Negative capacitance regime in ferroelectrics demystified from nonequilibrium molecular dynamics

Negative capacitance regime in ferroelectrics is associated with the negative slope in the dependence of the polarization on electric field and is prohibited by the laws of equilibrium thermodynamics. Very recently, it developed into a controversial but very productive area with a promise for next generation low-power transistors. At least, in part, the controversies arise from the lack of direct experimental insight into the mechanism for such phenomena. Although the experimental insight is, indeed, challenging at present, we probe the effect in ferroelectric PbTiO 3 bulk and a thin film with realistic first-principles-based nonequilibrium molecular dynamics simulations which model most closely experimental conditions. Here, the following insights emerge from the simulations: (i) the ferroelectric never enters the thermodynamically forbidden region even under extreme rates of electric-field application; (ii) the negative slope naturally arises when the compensating charge is either scarce or slowed down by the external circuitry; (iii) the negative slope is well controlled by the interplay between the applied step pulse field and the amount of the screening charge. The latter can be used to design the negative capacitance regime.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

RELAP-7 Application and Enhancement for FLEX Strategies and ATF Behavior under Extended Loss of AC Power Conditions

This report summarizes the results of a three-year research project sponsored by the U.S. Department of Energy (DOE) Nuclear Energy University Program (NEUP) to enhance and apply the RELAP-7 code by adding and improving several important components (e.g., a mechanistic Reactor Core Isolation Cooling (RCIC) system model) for thermal hydraulic studies of LWRs under ELAP conditions and evaluating the time available for transition to portable FLEX equipment. The project team included University of Massachusetts–Lowell, The Ohio State University, Texas A&M University, Idaho National Laboratory and Oak Ridge National Laboratory. In the Fukushima accident, it was found that the RCIC system played a crucial role in delaying core meltdown by almost three days in Fukushima Daiichi Unit 2, because of self-regulated operation of the steam driven RCIC turbine-pump injection system. Steam flow in the convergent-divergent nozzles of the RCIC Terry turbine is two-phase non-equilibrium transonic flow with homogenous nucleation condensation. To more accurately predict the dynamic process and behavior of the transonic compressible steam flow, a one-dimensional transient two-phase analytical model is presented. A simplified four-fluid model was employed in the present work with the consideration of four separate fluid fields: vapor, liquid film, entrained droplets and condensed droplets. The mass, momentum and energy interactions between the fluids were considered and modeled. An extended seven-equation non-equilibrium critical flow model was developed to obtain the critical pressure and velocities of each phase at the nozzle throat. To predict the wetness in the divergent section, a mechanistic nucleation condensation model was integrated in the nozzle analysis model, considering the generation and consequent growth of droplets. The governing differential equations on a staggered grid were discretized using the second-order Lax-Wendroff scheme with a flux limiter, and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE) algorithm was employed to solve the discrete linear system. To demonstrate the predictability and reliability of the physical models and the numerical method proposed in the present work, three representative nozzles were modeled and simulated. The results show good agreement with the available experimental data, even for condensation shock. Then, the 1D nozzle model was employed to obtain nozzle flow tables of the Terry turbine nozzle for different working pressures which can cover the operation pressure range of the RCIC system. A mechanistic RCIC turbine-pump system model was developed and implemented in the system code TRACE to simulate dynamic responses of the RCIC system under Beyond Design Basis Accident (BDBA) conditions. The turbine-pump governing equations are based on the control volume approach of the angular momentum balance. The physics based mechanistic RCIC model was developed using the TRACE control system components (i.e., signal variables, control blocks, and tables), and incorporated into a TRACE boiling water reactor (BWR) model. The TRACE model in this report has a detailed nodalization of the reactor pressure vessel (RPV), and all of the major flow paths and system components, including the safety relief valves (SRVs) and the containment suppression pool and drywell. Based on the nozzle flow tables generated from the 1D nozzle model developed, the turbine drive torque can be calculated from table lookup. Since the detailed specifications of the RCIC pump are unavailable, the homologous curves for a Bingham pump were used in the current pump component. A station black-out (SBO) accident test problem was selected to demonstrate the TRACE RCIC model. The short-term SBO simulations were performed for two cladding materials: Zircaloy and FeCrAl, to demonstrate the effect of the accident tolerant fuel cladding on fuel heat-up under BDBA accident conditions. The wetwell plays a vital safety role in SBO and other BWR accident scenarios in that it can reduce containment pressure and supply additional core make-up water. The suppression pool temperature distribution has a very large impact on both RPV and containment pressure. Thus, another novel contribution of the project comes mainly from an improved, systems-level wetwell model which can capture buoyancy-induced thermal stratification effects due to steam injection and condensation. A two-zone stratified wetwell model has been implemented in RELAP-7 and some results from that model are presented. This wetwell model is capable of simulating thermal stratification due to a low steam mass injection rate. With a low mass flow rate, the model assumes that all the steam condenses within the pipe and the resulting plume can be approximated with a purely buoyant, heat-source driven model. The wetwell model developed with these assumptions is adequate to simulate slow transients such as extended SBO transients.

42 ENGINEERING↗
DOE OSTI2022

Correlations for the specific heat capacity of ( U x Pu 1 - x ) 1 - y Gd y O 2 - z derived from molecular dynamics

We report UO 2 is the primary conventional fuel used in most nuclear reactors with Gd 2 O 3 commonly added as a burnable absorber to produce a more level power distribution in the reactor core at the beginning of operation. It can also be mixed with other actinide oxides to produce mixed oxide (MOx) fuel. In this study, molecular dynamics simulations were used to predict the specific heat capacity of Gd-doped PuO 2 , UO 2 and (U, Pu)O 2 MOx accommodating Gd 3+ substituted at cation sites via two charge compensation mechanisms - oxygen vacancy formation and the oxidation of U 4+ to U 5+ . The specific heat capacity values for PuO 2 and UO 2 are in good agreement with other studies showing a distinct peak at high temperatures - above 1800 K. As Gd 3+ is added, the peak height reduces for each composition considered. An analytical fit was applied to the data where Gd 3+ was fully charge compensated by either oxygen vacancies or U 5+ . The expression was then validated by predicting the specific heat capacity for three compositions of (Ux Pu 1-x ) 1-y Gd y O 2-z containing both oxygen vacancies and U 5+ , and compared to molecular dynamics data.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Structural Characterization and Dynamics of AdhE Ultrastructures from Clostridium thermocellum Show a Containment Strategy for Toxic Intermediates

Clostridium thermocellum, a cellulolytic thermophilic anaerobe, is considered by many to be a prime candidate for the realization of consolidated bioprocessing (CBP) and is known as an industry standard for biofuel production. C. thermocellum is among the best biomass degraders identified to date in nature and produces ethanol as one of its main products. Many studies have helped increase ethanol titers in this microbe; however, ethanol production using C. thermocellum is still not economically viable. Therefore, a better understanding of its ethanol synthesis pathway is required. The main pathway for ethanol production in C. thermocellum involves the bifunctional aldehyde-alcohol dehydrogenase (AdhE). To better understand the function of the C. thermocellum AdhE, we used cryo-electron microscopy (cryo-EM) to obtain a 3.28 A structure of the AdhE complex. This high-resolution structure, in combination with molecular dynamics simulations, provides insight into the substrate channeling of the toxic intermediate acetaldehyde, indicates the potential role of C. thermocellum AdhE to regulate activity and cofactor pools, and establishes a basis for future engineering studies. The containment strategy found in this enzyme offers a template that could be replicated in other systems where toxic intermediates need to be sequestered to increase the production of valuable biochemicals.

09 BIOMASS FUELS↗
DOE OSTI2026

Molecular dynamics based study on the effects of cation size on the local structure and diffusion in polymerized ionic liquids

We have used coarse-grained molecular dynamics simulations to understand the effect of cation size on ion diffusion in polymerized ionic liquids at temperatures well above the glass transition temperature of the polymers. We investigated dependencies of the diffusion constant on the cation radius and static dielectric constant and interpreted these results in terms of underlying structural changes, and decoupling phenomenon. We have found non-monotonic effects of the cation radius on the diffusion constant, with a maximum at an intermediate radius resulting from two different size dependent effects. Changes in the radial distribution functions characterizing spatial distribution of cations with respect to anions result in a monotonic increase in the coordination number. This yields an increase in the effective interaction energy as a function of cation radius, which can explain the decrease in the diffusion constant. This increase is counteracted at small cation radii by a corresponding increase in decoupling of cation and anion dynamics. In addition, the diffusion constant is found to increase on increasing the static dielectric constant for all radii of the cations. These results highlight the significant role of decoupling in designing polymerized ionic liquids with an enhanced diffusion constant of cations.

Gillespie, Colin [ORNL] (ORCID:0000000320287080)↗
DOE OSTI2023

Molecular dynamics study of grain boundary and radiation effects on tritium population and diffusion in zirconium

Tritium population thermodynamics and transport kinetics critically define the tritium storage performance of zirconium tritides that can be used for a variety of nuclear applications including tritium-producing burnable absorber rods. Both thermodynamic and kinetic properties can be sensitive to grain sizes of materials and can be significantly altered by irradiated defects during operation under the reactor environments. A thorough experimental characterization of how these properties evolve under different reactor conditions and different initial grain structures is extremely challenging. Here molecular dynamics simulations are used to investigate tritium population and diffusion in zirconium with and without different planar symmetric and asymmetric tilt grain boundaries and irradiated defects. Here, we found that in addition to trapping tritium, the most significant effect of planar grain boundaries is to increase tritium diffusivity on the boundary plane. Furthermore, fine grain structures are found to mitigate the change of tritium diffusivity due to irradiated point defects as these point defects are likely to migrate to and sink at grain boundaries.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Zinc binding alters the conformational dynamics and drives the transport cycle of the cation diffusion facilitator YiiP

YiiP is a secondary transporter that couples Zn 2+ transport to the proton motive force. Structural studies of YiiP from prokaryotes and Znt8 from humans have revealed three different Zn 2+ sites and a conserved homodimeric architecture. These structures define the inward-facing and outward-facing states that characterize the archetypal alternating access mechanism of transport. To study the effects of Zn 2+ binding on the conformational transition, we use cryo-EM together with molecular dynamics simulation to compare structures of YiiP from Shewanella oneidensis in the presence and absence of Zn 2+ . To enable single-particle cryo-EM, we used a phage-display library to develop a Fab antibody fragment with high affinity for YiiP, thus producing a YiiP/Fab complex. To perform MD simulations, we developed a nonbonded dummy model for Zn 2+ and validated its performance with known Zn 2+ -binding proteins. Using these tools, we find that, in the presence of Zn 2+ , YiiP adopts an inward-facing conformation consistent with that previously seen in tubular crystals. After removal of Zn 2+ with high-affinity chelators, YiiP exhibits enhanced flexibility and adopts a novel conformation that appears to be intermediate between inward-facing and outward-facing states. This conformation involves closure of a hydrophobic gate that has been postulated to control access to the primary transport site. Comparison of several independent cryo-EM maps suggests that the transition from the inward-facing state is controlled by occupancy of a secondary Zn 2+ site at the cytoplasmic membrane interface. This work enhances our understanding of individual Zn 2+ binding sites and their role in the conformational dynamics that govern the transport cycle.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Immobilizing Pertechnetate in Ettringite via Sulfate Substitution

Until recently, a nuclear waste form did not exist for the immobilization and long-term storage of the anionic form of radionuclide technetium-99 (99Tc), pertechnetate (TcO 4 - ). Across international nuclear waste sites, the inability to capture TcO 4 - poses a threat to the local environment due to the anion’s long half-life and environmental mobility under environmentally oxic conditions. To combat this challenge, most waste forms would capture 99Tc using materials that reduce TcO 4 - to Tc(IV) species, e.g. 99Tc oxides or 99Tc sulfides. However, this is a short-lived solution unless reducing capacities can be regenerated to last over the long life expectancy of some nuclear waste forms (~10,000 years). In light of recent experimental observations gained from the development of cementitious waste forms and unconfirmed hypotheses from the literature, this work explores how ettringite, a common mineral formed in cementitious materials with demonstrated success in incorporating oxyanions, may be used to immobilize TcO 4 - directly and over the life-time of the waste form. Using batch precipitation experiments, solid phase characterization techniques, and ab initio molecular dynamics simulations, we demonstrate herein successful incorporation of TcO 4 - into ettringite via substitution for SO 4 2- within the ettringite crystal structure. Here, these results are first of their kind and confirm the hypothesis that ettringite acts as an immobilizing mineral within cementitious waste forms for the highly mobile TcO 4 - oxyanion, enhancing the overall performance of these waste forms over time.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Polyamorphism in a solute-lean Al–Ce metallic glass

Polyamorphism discovered in lanthanide-rich metallic glasses (MGs) has been attributed to the electronic transition of the lanthanide element as a solvent element. In this work, we report that pressure-induced polyamorphism still exists in a Ce-poor Al 93 Ce 7 binary MG where the 4f electron element serves as a solute and solute–solute avoidance is expected. The polyamorphic transition, observed by in situ high-pressure synchrotron x-ray diffraction, is accompanied by a volume collapse of ~0.78% and occurs over a narrow pressure range from ~0.8 to ~1.8 GPa. Further synchrotron Ce L 3 -edge x-ray absorption spectroscopy measurements reveal that pressure-induced 4f electron delocalization underlies the polyamorphic transition. Molecular dynamics simulations confirm that the Ce atoms in the MG are completely isolated by the solvent Al atoms. This result demonstrates that 4f element-bearing alloys with extremely dilute concentrations can also exhibit polyamorphic states originating from electronic transitions, extending the compositional space of polyamorphism of MGs into very dilute regions. Our work suggests that tunable properties under compressive stress could be achieved in MGs by even minor doping of elements prone to electronic transitions.

36 MATERIALS SCIENCE↗
DOE OSTI2023

MD Simulation of Water Using a Rigid Body Description Requires a Small Time Step to Ensure Equipartition

In simulations of aqueous systems, it is common to freeze the bond vibration and angle bending modes in water to allow for a longer time step δt for integrating the equations of motion. Thus, δt = 2 fs is often used in simulating rigid models of water. We simulate the SPC/E model of water using δt from 0.5 to 3.0 fs and up to 4 fs using hydrogen mass repartitioning. In these simulations, we find that for all but δt = 0.5 fs, equipartition is not obtained between translational and rotational modes, with the rotational modes exhibiting a lower temperature than the translation modes. To probe the reasons for the lack of equipartition, we study the autocorrelation of the translational velocity of the center of mass and the angular velocity of the rigid water molecule, respectively. Here, we find that the rotational relaxation occurs on a timescale comparable to vibrational periods, calling into question the original motivations for freezing the vibrations. Furthermore, a time step with δt ≥ 1 fs is not able to capture accurately the fast rotational relaxation, which reveals its impact as an effective slowing-down of rotational relaxation. The fluctuation–dissipation relation then leads to the conclusion that the rotational temperature should be cooler for δt greater than the reference value of 0.5 fs. Consideration of fluctuation–dissipation in equilibrium molecular dynamics simulations also emphasizes the need to capture the temporal evolution of fluctuations with fidelity and the role of δt in this regard. The time step also influences the solution thermodynamic properties: both the mean system potential energies and the excess entropy of hydration of a soft repulsive cavity are sensitive to δt.

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

Demystify radiation-enhanced hydrogen isotope diffusion in Fe-Ni-Cr austenitic stainless steels

Understanding and containing hydrogen isotope diffusion is crucial for many nuclear applications. In situ experiments have consistently shown that radiation significantly enhances isotope diffusion in austenitic stainless steels. Despite extensive research, the mechanism behind this phenomenon remains elusive, as most radiation-induced defects (e.g., vacancies, dislocations, and grain boundaries) typically trap hydrogen, thereby slowing diffusion. While grain boundaries may increase in-plane diffusivity and interstitials may enhance diffusion due to material swelling, these effects are relatively minor. Utilizing an Fe-Ni-Cr-H interatomic potential for stainless steels, we conducted extensive molecular dynamics simulations to investigate the origins of radiation-enhanced diffusion. Here, our findings reveal that when a system is resolidified, mimicking defects created by radiation displacements, the resulting structure contains a mixture of phases, boundaries, and dislocation networks. This defective structure significantly increases hydrogen diffusivity, enhancing it by approximately 1.7 times at 900 K. These results suggest that the complex defect structures formed during radiation displacements are the primary drivers of the observed diffusion enhancement, providing valuable insights into the mechanisms underlying radiation-enhanced diffusion in nuclear materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Enhanced Activity from Coordinatively Unsaturated and Dynamic Zeolite-Bound Organoscandium Species

Scandium borohydride grafted into the micropores of the faujasite zeolite HY30 catalyzes the C–H borylation of benzene, whereas silica-grafted species are inactive. This catalytic activity may originate from grafting at a Brønsted acid site leading to a more electron-deficient rare earth center. Herein, we apply multinuclear double-resonance nuclear magnetic resonance (NMR) experiments to probe the structure and dynamics of zeolite- and silica-bound scandium borohydride complexes. The experiments reveal that scandium centers located within the zeolite micropores, in proximity to Al-created Brønsted sites, are more dynamic than rigid scandium sites grafted on silanols. Through a combination of NMR and molecular dynamics simulations, we show that the coordination of the scandium in the zeolite is labile, with the metal exchanging between two binding sites. As a result, the weak electron donation from the support that enables the movement of the Sc center leads to the formation of an undercoordinated metal center that cannot exist on silica, ultimately leading to the new catalytic activity of the species.

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

Effects of carbon concentration on the local atomic structure of amorphous GST

Ge-Sb-Te (GST) alloys are leading phase-change materials for data storage due to the fast phase transition between amorphous and crystalline states. Ongoing research aims at improving the stability of the amorphous phase to improve retention. This can be accomplished by the introduction of carbon as a dopant to Ge 2 Sb 2 Te 5 , which is known to alter the short- and mid-range structure of the amorphous phase and form covalently bonded C clusters, both of which hinder crystallization. The relative importance of these processes as a function of C concentration is not known. Here we used molecular dynamics simulation based on density functional theory to study how carbon doping affects the atomic structure of GST-C. Carbon doping results in an increase in tetrahedral coordination, especially of Ge atoms, and this is known to stabilize the amorphous phase. We observe an unexpected, non-monotonous trend in the number of tetrahedral bonded Ge with the amount of carbon doping. Our simulations show an increase in the number of tetrahedral bonded Ge up to 5 at.% C, after which the number saturates and begins to decrease above 14 at.% C. The carbon atoms aggregate into clusters, mostly in the form of chains and graphene flakes, leaving less carbon to disrupt the GST matrix at higher carbon concentrations. Different degrees of carbon clustering can explain divergent experimental results for recrystallization temperature for carbon doped GST.

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

Water transport through monolayer fullerene membrane

Water transport through nanoporous materials is important in water treatment, desalination, and nanofiltration. Two-dimensional (2D) membranes such as porous graphene have been explored for high-permeance water transport. However, water transport through a new class of 2D membranes based on two-dimensional covalently linked fullerene monolayers has not been fully explored. Here we use classical molecular dynamics simulations to investigate both vapor and liquid water transport through a monolayer fullerene membrane. We find that a quasi-tetragonal phase fullerene membrane possesses the right pore size and geometry that allows fast water vapor transport (∼ 50 g m −2 day −1 Pa −1 ) and water liquid transport (∼ 2.0 g m −2 day −1 Pa −1 ). Furthermore, simulation of sea water transport through the fullerene membrane shows 100 % salt rejection. The much faster vapor transport rate is attributed to the funnel-shaped pore and the optimal size that allows free rotation of water molecules permeating through, while the slower liquid transport is due to the need to desolvate a water molecule to break its hydrogen-bond network across the hydrophobic pore. This work shows the great potential of using monolayer fullerene membranes as 2D membranes for fast and selective water transport.

Desalination↗
DOE OSTI2020

Neural network potential from bispectrum components: A case study on crystalline silicon

In this article, we present a systematic study on developing machine learning force fields (MLFFs) for crystalline silicon. While the main-stream approach of fitting a MLFF is to use a small and localized training set from molecular dynamics simulations, it is unlikely to cover the global features of the potential energy surface. Additionally, to remedy this issue, we used randomly generated symmetrical crystal structures to train a more general Si-MLFF. Furthermore, we performed substantial benchmarks among different choices of material descriptors and regression techniques on two different sets of silicon data. Our results show that neural network potential fitting with bispectrum coefficients as descriptors is a feasible method for obtaining accurate and transferable MLFFs.

36 MATERIALS SCIENCE↗