Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Interatomic potentials”

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

Inverse methods for design of soft materials

Functional soft materials, comprising colloidal and molecular building blocks that self-organize into complex structures as a result of their tunable interactions, enable a wide array of technological applications. Inverse methods provide a systematic means for navigating their inherently high-dimensional design spaces to create materials with targeted properties. Furthermore, while multiple physically motivated inverse strategies have been successfully implemented in silico, their translation to guiding experimental materials discovery has thus far been limited to a handful of proof-of-concept studies. In this perspective, we discuss recent advances in inverse methods for design of soft materials that address two challenges: (1) methodological limitations that prevent such approaches from satisfying design constraints and (2) computational challenges that limit the size and complexity of systems that can be addressed. Strategies that leverage machine learning have proven particularly effective, including methods to discover order parameters that characterize complex structural motifs and schemes to efficiently compute macroscopic properties from the underlying structure. We also highlight promising opportunities to improve the experimental realizability of materials designed computationally, including discovery of materials with functionality at multiple thermodynamic states, design of externally directed assembly protocols that are simple to implement in experiments, and strategies to improve the accuracy and computational efficiency of experimentally relevant models.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Self-assembly of a triply periodic continuous mesophase with Fddd symmetry in simple one-component liquids

Triply periodic continuous morphologies (networks) arising as a result of the microphase separation in block copolymer melts have so far never been observed self-assembled in systems of particles with spherically symmetric interaction. Here, we report a molecular dynamics simulation where two simple one-component liquids form upon cooling an equilibrium network with the Fddd space group symmetry. This complexity reduction in the liquid network formation in terms of the particle geometry and the number of components evidences the generic nature of this class of phase transition suggesting opportunities for producing these structures in a variety of new systems.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Simulating sticky particles: A Monte Carlo method to sample a stratification

Many problems in materials science and biology involve particles interacting with strong, short-ranged bonds that can break and form on experimental timescales. Treating such bonds as constraints can significantly speed up sampling their equilibrium distribution, and there are several methods to sample probability distributions subject to fixed constraints. We introduce a Monte Carlo method to handle the case when constraints can break and form. More generally, the method samples a probability distribution on a stratification: a collection of manifolds of different dimensions, where the lower-dimensional manifolds lie on the boundaries of the higher-dimensional manifolds. We show several applications of the method in polymer physics, self-assembly of colloids, and volume calculation in high dimensions.

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

Xenon mobility in γ-uranium and uranium–molybdenum alloys

Diffusion in bcc uranium and U–Mo alloys is of great interest because fission gas and other fission products impact the performance of nuclear fuels. We investigate the mobility of xenon and molybdenum in bcc uranium (⁠$γ$-U) and metallic U–Mo alloys by calculating the migration energies of xenon and molybdenum for various local compositions using density functional theory. We also calculate the solute–vacancy binding energies of different solutes to vacancies in bcc uranium. We find that the solute–vacancy binding energy in bcc uranium is significantly higher than it is in other bcc metals (e.g., Fe and W). We also find that the migration energy of molybdenum is substantially higher than the migration energy of xenon, indicating that xenon is much more mobile than molybdenum in bcc uranium. The presence of molybdenum in the nearest-neighbor shell around a xenon atom typically increases the migration energy of xenon, which indicates a reduction of xenon mobility in U–Mo alloys compared to pure bcc uranium.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Activated relaxation in supercooled monodisperse atomic and polymeric WCA fluids: Simulation and ECNLE theory

Here, we combine simulation and Elastically Collective Nonlinear Langevin Equation (ECNLE) theory to study the activated relaxation in monodisperse atomic and polymeric Weeks–Chandler–Andersen (WCA) liquids over a wide range of temperatures and densities in the supercooled regime under isochoric conditions. By employing novel crystal-avoiding simulations, metastable equilibrium dynamics is probed in the absence of complications associated with size polydispersity. Based on a highly accurate structural input from integral equation theory, ECNLE theory is found to describe well the simulated density and temperature dependences of the alpha relaxation time of atomic fluids using a single system-specific parameter, a c , that reflects the nonuniversal relative importance of local cage and collective elastic barriers. For polymer fluids, the explicit dynamical effect of local chain connectivity is modeled at the fundamental dynamic free energy trajectory level based on a different parameter, N c , that quantifies the degree of intramolecular correlation of bonded segment activated barrier hopping. For the flexible chain model studied, a physically intuitive value of N c ≈ 2 results in good agreement between simulation and theory. A direct comparison between atomic and polymeric systems reveals that chain connectivity can speed up activated segmental relaxation due to weakening of equilibrium packing correlations but can slow down relaxation due to local bonding constraints. The empirical thermodynamic scaling idea for the alpha time is found to work well at high densities or temperatures but fails when both density and temperature are low. The rich and subtle behaviors revealed from simulation for atomic and polymeric WCA fluids are all well captured by ECNLE theory.

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

Elucidating size effects on the yield strength of single-crystal Cu via the Richtmyer–Meshkov instability

Capturing the dynamic response of a material under high strain-rate deformation often demands challenging and time consuming experimental effort. While shock hydrodynamic simulation methods can aid in this area, a priori characterizations of the material strength under shock loading and spall failure are needed in order to parameterize constitutive models needed for these computational tools. Moreover, parameterizations of strain-rate-dependent strength models are needed to capture the full suite of Richtmyer–Meshkov instability (RMI) behavior of shock compressed metals, creating an unrealistic demand for these training data solely on experiments. Herein, we sweep a large range of geometric, crystallographic, and shock conditions within molecular dynamics (MD) simulations and demonstrate the breadth of RMI in Cu that can be captured from the atomic scale. In this work, yield strength measurements from jetted and arrested material from a sinusoidal surface perturbation were quantified as Y RMI = 0.787 ± 0.374 GPa, higher than strain-rate-independent models used in experimentally matched hydrodynamic simulations. Defect-free, single-crystal Cu samples used in MD will overestimate Y RMI , but the drastic scale difference between experiment and MD is highlighted by high confidence neighborhood clustering predictions of RMI characterizations, yielding incorrect classifications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Chemical order transitions within extended interfacial segregation zones in NbMoTaW

Interfacial segregation and chemical short-range ordering influence the behavior of grain boundaries in complex concentrated alloys. In this study, we use atomistic modeling of a NbMoTaW refractory complex concentrated alloy to provide insight into the interplay between these two phenomena. Hybrid Monte Carlo and molecular dynamics simulations are performed on columnar grain models to identify equilibrium grain boundary structures. Here our results reveal extended near-boundary segregation zones that are much larger than traditional segregation regions, which also exhibit chemical patterning that bridges the interfacial and grain interior regions. Furthermore, structural transitions pertaining to an A2-to-B2 transformation are observed within these extended segregation zones. Both grain size and temperature are found to significantly alter the widths of these regions. An analysis of chemical short-range order indicates that not all pairwise elemental interactions are affected by the presence of a grain boundary equally, as only a subset of elemental clustering types are more likely to reside near certain boundaries. The results emphasize the increased chemical complexity that is associated with near-boundary segregation zones and demonstrate the unique nature of interfacial segregation in complex concentrated alloys.

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

Combined molecular and spin dynamics simulation of BCC iron with vacancy defects

Utilizing an atomistic computational model, which handles both translational and spin degrees of freedom, combined molecular and spin dynamics simulations have been performed to investigate the effect of vacancy defects on spin wave excitations in ferromagnetic iron. Fourier transforms of space- and time-displaced correlation functions yield the dynamic structure factor, providing characteristic frequencies and lifetimes of the spin wave modes. A comparison of the system with a 5% vacancy concentration with pure lattice data shows a decrease in frequency and a decrease in lifetime for all transverse spin wave excitations observed. In addition, the clearly defined transverse spin wave excitations are distorted with the introduction of vacancy defects, and we observe reduced excitation lifetimes due to increased magnon–magnon scattering. We observe further evidence of increased magnon–magnon scattering, as the peaks in the longitudinal spin wave spectrum become less distinct. Finally, similar impacts are observed in the vibrational subsystem, with a decrease in characteristic phonon frequency and flattening of lattice excitation signals due to vacancy defects.

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

Medium-range order and compositional correlation in metallic glasses

The compositional atomic ordering in metallic glasses was studied by simulation focusing on the medium-range order (MRO). Many metallic alloy liquids and glasses show MRO characterized by the oscillations in the atomic pair-distribution function (PDF) beyond the first peak, which decay exponentially with distance. To study the effects of the local chemical order on MRO, we examine the compositionally resolved PDF and its MRO for models of various binary metallic alloy glasses. We show that compositional ordering is limited mostly to the nearest-neighbor atoms and the MRO is largely independent of the compositional order. For some elements that strongly repel each other in the alloy, a second MRO periodicity is observed owing to the distinct correlations among them. These results are discussed in light of the idea that the MRO oscillations in the PDF describe the correlations in the atomic density fluctuations, rather than the detailed local atomic structure.

Atomic structure↗
DOE OSTI2022

Modification of a force field for molecular dynamics simulations of silicon etching by chlorine atoms

Here, a modified classical molecular dynamics (MD) force field is presented for the interaction between silicon (Si) and chlorine (Cl). The original version of the force field is shown to significantly overestimate the probability of Si etching by thermal Cl atoms. However, the modified force field corrects this problem and results in generally good agreement with experimental data. Further, it is shown that while the modification of the force field improves the prediction of Si spontaneous etching with Cl atoms, it does not degrade predictions of atomic-layer etching of Si with Cl[Formula: see text] molecules.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Kinetic Monte Carlo simulations of aging in δ -Pu

We have developed a first-passage kinetic Monte Carlo approach for materials aging to investigate the sensitivity of void swelling to model parameters, including helium bubble density and size distribution. In addition to explicitly accounting for the spatial distribution of individual point defects, bubbles, and voids, our approach can simulate total doses equivalent to 100 years of natural aging on statistically representative volumes of materials. This technique enables us to study the effects on swelling and radiation damage evolution due to temperature and dose rate (as altered in artificially aged experiments), differences in effective interaction radii between vacancies and interstitials, and varying defect diffusion activation energies, while providing more detailed information than previous rate-equation based approaches. In conclusion, our results indicate that spatial effects that are not modeled in mean-field rate theories could play a significant role in void swelling initiation and growth for certain regimes of model parameters.

Actinides↗
DOE OSTI2021

Implementing a neural network interatomic model with performance portability for emerging exascale architectures

The two main thrusts of computational science are increasingly accurate predictions and faster calculations; to this end, the zeitgeist in molecular dynamics (MD) simulations is pursuing machine learned and data driven interatomic models, e.g. neural network potentials, and novel hardware architectures, e.g. GPUs. Current implementations of neural network potentials are orders of magnitude slower than traditional interatomic models and while looming exascale computing offers the ability to run large, accurate simulations with these models, achieving portable performance for MD with new and varied exascale hardware requires rethinking traditional algorithms, using novel data structures, and library solutions. We re-implement a neural network interatomic model in CabanaMD, an MD proxy application, built on libraries developed for performance portability. Our implementation shows significantly improved thread scaling in this complex kernel as compared to a current LAMMPS implementation, across both strong and weak scaling. Our single-source solution enables simulations up to 20 million atoms on a single CPU node and 4 million atoms with improved performance on a single GPU. Furthermore, we also explore parallelism and data layout choices (using flexible data structures called AoSoAs) and their effect on performance, seeing up to ~50% and ~5% improvements in performance on a GPU by choosing the right level of parallelism and data layout respectively.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Multi-state Catalysts Modulated by Mechanical Force (Final Report)

The development of more efficient catalytic processes and new approaches to control catalytic activity and selectivity are central to the realization of more selective, atom economic, and energy efficient routes to value added chemicals and polymers. The reactivity and selectivity of a transition metal catalyst is intimately related to the ligand-sphere geometry and, in many cases, the ideal ligand geometry for one step of a catalytic cycle is poorly matched to the ideal ligand geometry for another, resulting in sub-optimal efficiency. Macroscopic mechanical forces are both large, potentially much larger than interatomic forces, and are directional and localized to an extent that differentiates them from other forms of energy input such as heat or light. As such, mechanical force represents a heretofore untapped approach to modulate catalyst geometry, with the potential to reversibly modulate catalyst geometry on the timescale of catalytic turnover or monomer enchainment. This project has addressed the fundamental challenges in material-to-molecule strain coupling associated with the development of a new class of mechanically responsive catalysts (mechanocatalysts) in which active organotransition metal catalysts are strategically embedded in a flexible polymer network such that application of external mechanical force (stretching or deformation) leads to modulation of catalyst geometry, and hence reactivity and selectivity. Our efforts during the tenure of this grant were directed toward the elucidation of force-reactivity relationships of elementary transformations that occur within the first coordination sphere of a transition metal complex employing stiff stilbene photoswitches tethered to a flexible bidentate phosphine ligand derived from MeOBiphep as molecular force probes which provide a range of compressive and extension forces to the coupled transition metal complex depending on the geometry of the stiff stilbene and length of the tethering chains. During the tenure of this grant, we have quantified the rate of C(sp 2 )-C(sp 2 ) reductive elimination from platinum(II) diaryl complexes containing bis(phosphine) force probe ligands as a function of mechanical force; compressive forces decreased the rate of reductive elimination whereas extension forces increased the rate relative to the strain-free MeOBiphep complex with a 3.4-fold change in rate over a ~290 pN range of restoring forces. In a similar manner, we have quantified the rate of oxidative addition of bromobenzene to low-ligated palladium(0) complexes containing force probe ligands as a function of mechanical force; compressive forces increase the rate of oxidative addition, whereas tensile forces decrease the rate with a ~6 fold change in rate across ~340 pN of force applied to the complexes. In both cases, experimental and computational analyses argue strongly against any significant force-induced perturbation of ground state geometry within the first coordination sphere of the reactant complexes. Rather, the force/rate behavior observed for these transformations across these ranges of forces is attributed to the coupling of force to the nuclear motion comprising the reaction coordinates for reductive elimination and oxidative addition. These results together inform the development of catalysts whose activity can be tuned by an external force that is adjusted within a catalytic cycle and suggest opportunities to experimentally map geometry changes associated with reactions in transition metal complexes and potential strategies for force-modulated catalysis.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2023

Shadow energy functionals and potentials in Born–Oppenheimer molecular dynamics

In Born–Oppenheimer molecular dynamics (BOMD) simulations based on the density functional theory (DFT), the potential energy and the interatomic forces are calculated from an electronic ground state density that is determined by an iterative self-consistent field optimization procedure, which, in practice, never is fully converged. The calculated energies and forces are, therefore, only approximate, which may lead to an unphysical energy drift and instabilities. Here, we discuss an alternative shadow BOMD approach that is based on backward error analysis. Instead of calculating approximate solutions for an underlying exact regular Born–Oppenheimer potential, we do the opposite. Instead, we calculate the exact electron density, energies, and forces, but for an underlying approximate shadow Born–Oppenheimer potential energy surface. In this way, the calculated forces are conservative with respect to the approximate shadow potential and generate accurate molecular trajectories with long-term energy stabilities. We show how such shadow Born–Oppenheimer potentials can be constructed at different levels of accuracy as a function of the integration time step, δt, from the constrained minimization of a sequence of systematically improvable, but approximate, shadow energy density functionals. For each energy functional, there is a corresponding ground state Born–Oppenheimer potential. These pairs of shadow energy functionals and potentials are higher-level generalizations of the original “zeroth-level” shadow energy functionals and potentials used in extended Lagrangian BOMD. The proposed shadow energy functionals and potentials are useful only within this extended dynamical framework, where also the electronic degrees of freedom are propagated as dynamical field variables together with the atomic positions and velocities. The theory is quite general and can be applied to MD simulations using approximate DFT, Hartree–Fock, or semi-empirical methods, as well as to coarse-grained flexible charge models.

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

Temperature-dependent rotationally inelastic collisions of OH – and He

We have studied the fundamental rotational relaxation and excitation collision of OH – J = 0 ↔ 1 with helium at different collision energies. Using state-selected photodetachment in a cryogenic ion trap, the collisional excitation of the first excited rotational state of OH – has been investigated and absolute inelastic collision rate coefficients have been extracted for collision temperatures between 20 and 35 K. Furthermore, the rates are compared with accurate quantum scattering calculations for three different potential-energy surfaces. Good agreement is found within the experimental accuracy, but the experimental trend of increasing collision rates with temperature is only in part reflected in the calculations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

van der Waals interaction as the starting point for an effective field theory

We consider the system of three 4 He atoms to assess whether a pure van der Waals potential can be used as a starting point for an effective field theory to describe three-body processes in ultracold atomic systems. Using a long-range van der Waals interaction in combination with short-distance two-body counterterms, we analyze the dependence of two- and three-body observables on the short-distance regulator that is required due to the singular nature of the van der Waals interaction. We benchmark our approach with results obtained with the realistic 4 He- 4 He LM2M2 potential and find good agreement. We furthermore show that in this effective field theory approach no three-body force is required at leading order and that universal van der Waals physics leads to a universal correlation between three-body observables in the absence of an Efimov three-body parameter.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Composition-transferable machine learning potential for LiCl-KCl molten salts validated by high-energy x-ray diffraction

Unraveling the liquid structure of multicomponent molten salts is challenging due to the difficulty in conducting and interpreting high-temperature diffraction experiments. Here, motivated by this challenge, we developed composition-transferable Gaussian approximation potential (GAP) for molten LiCl-KCl. A DFT-SCAN accurate GAP is active-learned from only ~1100 training configurations drawn from 10 unique mixture compositions enriched with metadynamics. The GAP-computed structures show strong agreement across high-energy x-ray diffraction experiments, including for a eutectic not explicitly included in model training, thereby opening the possibility of composition discovery.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Phase Diagram of a Deep Potential Water Model

Using the Deep Potential methodology, we construct a model that reproduces accurately the potential energy surface of the SCAN approximation of density functional theory for water, from low temperature and pressure to about 2400 K and 50 GPa, excluding the vapor stability region. The computational efficiency of the model makes it possible to predict its phase diagram using molecular dynamics. Satisfactory overall agreement with experimental results is obtained. Here, the fluid phases, molecular and ionic, and all the stable ice polymorphs, ordered and disordered, are predicted correctly, with the exception of ice III and XV that are stable in experiments, but metastable in the model. The evolution of the atomic dynamics upon heating, as ice VII transforms first into ice VII" and then into an ionic fluid, reveals that molecular dissociation and breaking of the ice rules coexist with strong covalent fluctuations, explaining why only partial ionization was inferred in experiments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗