Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Molecular 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,405 records · Page 78
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Molecular-level solvation and selectivity behavior of Na + , K + , and Li + within glycerol-derived solvents

Ether-based acyclic solvents have a tendency to confine positively charged atoms, forming structures called podands. Due to the stability of these structures, podands have been identified as promising candidates for several different electrochemical applications. In this work, we use molecular dynamics simulations to investigate the solvation and selectivity of Na + , K + , and Li + by different glycerol-derived solvents. The structural properties suggest that the solvent-cation association is stable over a very wide temperature range. Further, a strong correlation is observed between the more electronegative oxygen sites and preferential interaction with the cations, which suggests that the electrostatic oxygen-ion interactions dominate the coordination between these species. Significantly more intense interactions are observed between the solvents and Li + , versus Na + or K + . These molecular-level thermodynamic insights provide fundamental support for further exploring the performance of glycerol-derived podand complexes in Li + capture applications and/or electrolyte formulations.

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

Machine learning approaches for structural and thermodynamic properties of a Lennard-Jones fluid

Predicting the functional properties of many molecular systems relies on understanding how atomistic interactions give rise to macroscale observables. However, current attempts to develop predictive models for the structural and thermodynamic properties of condensed-phase systems often rely on extensive parameter fitting to empirically selected functional forms whose effectiveness is limited to a narrow range of physical conditions. Here, we illustrate how these traditional fitting paradigms can be superseded using machine learning. Specifically, we use the results of molecular dynamics simulations to train machine learning protocols that are able to produce the radial distribution function, pressure, and internal energy of a Lennard-Jones fluid with increased accuracy in comparison to previous theoretical methods. The radial distribution function is determined using a variant of the segmented linear regression with the multivariate function decomposition approach developed by Craven et al. [J. Phys. Chem. Lett. 11, 4372 (2020)]. The pressure and internal energy are determined using expressions containing the learned radial distribution function and also a kernel ridge regression process that is trained directly on thermodynamic properties measured in simulation. The presented results suggest that the structural and thermodynamic properties of fluids may be determined more accurately through machine learning than through human-guided functional forms.

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

Effects of Coulomb coupling on friction in strongly magnetized plasmas

The friction force on a test particle traveling through a plasma that is both strongly coupled and strongly magnetized is studied using molecular dynamics simulations. In addition to the usual stopping power component aligned antiparallel to the velocity, a transverse component that is perpendicular to both the velocity and Lorentz force is observed. This component, which was recently discovered in weakly coupled plasmas, is found to increase in both absolute and relative magnitude in the strongly coupled regime. Strong coupling is also observed to induce a third component of the friction force in the direction of the Lorentz force. These first-principles simulations reveal novel physics associated with collisions in strongly coupled, strongly magnetized, plasmas that are not predicted by existing kinetic theories. The effect is expected to influence macroscopic transport in a number of laboratory experiments and astrophysical plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

The pervasive impact of critical fluctuations in liquid–liquid extraction organic phases

Liquid-liquid extraction is an essential chemical separation technique where polar solutes are extracted from an aqueous phase into a nonpolar organic solvent by amphiphilic extractant molecules. A fundamental limitation to the efficiency of this important technology is third phase formation, wherein the organic phase splits upon sufficient loading of polar solutes. The nanoscale drivers of phase splitting are challenging to understand in the complex hierarchically structured organic phases. In this study, we demonstrate that the organic phase structure and phase behavior are fundamentally connected in a way than can be understood with critical phenomena theory. For a series of binary mixtures of trialkyl phosphate extractants with linear alkane diluents, we combine small angle x-ray scattering and molecular dynamics simulations to demonstrate how the organic phase mesostructure over a wide range of compositions is dominated by critical concentration fluctuations associated with the critical point of the third phase formation phase transition. These findings reconcile many longstanding inconsistencies in the literature where small angle scattering features, also consistent with such critical fluctuations, were interpreted as reverse micellar-like particles. Altogether, this study shows how the organic phase mesostructure and phase behavior are intrinsically linked, deepening our understanding of both and providing a new framework for using molecular structure and thermodynamic variables to control mesostructure and phase behavior in liquid-liquid extraction.

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

DeePMD-kit v2: A software package for deep potential models

DeePMD-kit is a powerful open-source software package that facilitates molecular dynamics simulations using machine learning potentials known as Deep Potential (DP) models. This package, which was released in 2017, has been widely used in the fields of physics, chemistry, biology, and material science for studying atomistic systems. The current version of DeePMD-kit offers numerous advanced features, such as DeepPot-SE, attention-based and hybrid descriptors, the ability to fit tensile properties, type embedding, model deviation, DP-range correction, DP long range, graphics processing unit support for customized operators, model compression, non-von Neumann molecular dynamics, and improved usability, including documentation, compiled binary packages, graphical user interfaces, and application programming interfaces. This article presents an overview of the current major version of the DeePMD-kit package, highlighting its features and technical details. Additionally, this article presents a comprehensive procedure for conducting molecular dynamics as a representative application, benchmarks the accuracy and efficiency of different models, and discusses ongoing developments.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Multi-reward Reinforcement Learning Based Bond-Order Potential to Study Strain-Assisted Phase Transitions in Phosphorene

Here, we introduce a multi-reward reinforcement learning (RL) approach to train a flexible bond-order potential (BOP) for 2D phosphorene based on ab initio training data sets. Our approach is based on a continuous action space Monte Carlo tree search algorithm that is general and scalable and presents an efficient multiobjective optimization scheme for high-dimensional materials design problems. As a proof-of-concept, we deploy this scheme to parametrize multiple structural and dynamical properties of 2D phosphorene polymorphs. Our RL-trained BOP model adequately captures the structure, energetics, transformation barriers, equation of state, elastic constants, and phonon dispersions of various 2D P polymorphs. We use this model to probe the impact of temperature and strain rate on the phase transition from black (α-P) to blue phosphorene (β-P) through molecular dynamics simulations. A decrease in critical strain for this phase transition with increase in temperature is observed, and the underlying atomistic mechanisms are discussed.

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

A Graph Dynamical neural network approach for decoding dynamical states in ferroelectrics.

Ferroelectric materials such as BaTiO 3 show tremendous potential for emerging advances in memory devices, particular neuromorphic type devices. High density of memory can be obtained by stabilising polar domain walls at the nanoscale, regions of discontinuity between the well-defined polarization order parameter, but little is known about what controls their structure and dynamics in real nanoscale materials. Indeed, chiral polar domain walls have been observed in heterogeneous ferroelectrics, such as oxygen-deficient BaTiO 3 , but very little is known about how such polar-domains walls interact with defects. Indeed, a critical understanding of how dynamics of domain-walls depend on point-defects is crucial to create engineered ferroelectric memory devices. For this work, we perform large-scale simulations of nansocale domain-wall dynamics in pristine and defective BaTiO 3 using reactive force-field developed by us earlier (Phys. Chem. Chem. Phys., 2019, 21, 18240–18249), and capture their dynamical dependence on point defects using a graph dynamical neural-network approach, which we adapted to interrogate solids with well-defined order-parameters, and implemented using Pytorch based libraries. Our machine learning (ML) approach goes beyond the traditional post-processing methods to capture both spatial and temporal heterogeneities of large-scale molecular dynamics simulations of complex defective ferroelectric oxide materials. We crucially find that isolated oxygen vacancies introduce very localized spatial regions (~1–2 unit-cell in length) that show slow dipole relaxation due to formation of defect-dipoles, and that these defect-dipoles in turn slow the intrinsic dynamics of domain walls. Further, the roughness of domain walls, also influenced by vacancies, introduce dynamic heterogeneity along the domain-wall. As such we find a novel mechanism by which quenched disorder due to defects introduce dynamic heterogeneity thereby influencing response to external fields (particularly time varying fields) in a ferroelectric. Our study also emphasizes the need for creating digital twins of dynamical quantities to achieve autonomous in operando control of nanoscale switching.

42 ENGINEERING↗
DOE OSTI2023

Dynamics of collision-induced energy transfer

Collisions are of fundamental importance in many relevant physico-chemical phenomena in the gas phase and at the gas-surface interface. Collision-induced energy transfer not only offers a sensitive probe of interatomic interactions, but also impacts an array of processes ranging from transport to reactivity. This perspective highlights some recent advances on collision-induced energy transfer dynamics. In the gas phase, cold collisions revealed remarkable quantum effects and stereodynamics of bimolecular inelastic scattering, which motivated development of new state-to-state quantum scattering methods. For surface processes, experiments of hydrogen atom scattering from various surfaces combined with first-principles molecular dynamics simulations shed valuable light on the relative importance of adiabatic energy transfer to surface phonons and nonadiabatic energy transfer to surface electrons. We outline further challenges in these fields and expect more fruitful interplay between experiment and theory.

Cold collision↗
DOE OSTI2021

Automated discovery of a robust interatomic potential for aluminum

Machine learning, trained on quantum mechanics (QM) calculations, is a powerful tool for modeling potential energy surfaces. A critical factor is the quality and diversity of the training dataset. Here we present a highly automated approach to dataset construction and demonstrate the method by building a potential for elemental aluminum (ANI-Al). In our active learning scheme, the ML potential under development is used to drive non-equilibrium molecular dynamics simulations with time-varying applied temperatures. Whenever a configuration is reached for which the ML uncertainty is large, new QM data is collected. The ML model is periodically retrained on all available QM data. The final ANI-Al potential makes very accurate predictions of radial distribution function in melt, liquid-solid coexistence curve, and crystal properties such as defect energies and barriers. We perform a 1.3M atom shock simulation and show that ANI-Al force predictions shine in their agreement with new reference DFT calculations.

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

Extended Rice–Thomson analysis and atomistic simulations revealing grain boundary effects on fracture in refractory high-entropy alloys

Significance This work serves to extend the fundamental ductile vs. brittle fracture theory, specifically the Rice–Thomson criterion, by introducing a grain boundary ahead of an initiating crack which propagates at an oblique angle to impinge the boundary. Atomistic fracture simulations on two refractory complex concentrated alloys, the brittle NbMoTaW and the ductile Nb 45 Ta 25 Ti 15 Hf 15 , demonstrate qualitative correspondence with the extended Rice–Thomson criterion and experimental observations. Abstract Understanding how grain boundaries mediate fracture remains a critical challenge in designing ductile, high-performance refractory alloys. Here, we extend the Rice–Thomson criterion to account for the angle between cracks and the impinging grain boundaries (GBs), capturing the competition between intergranular fracture and dislocation-mediated plasticity. Using machine learning interatomic potentials, we performed molecular statics simulations to probe fracture mechanisms in nanocrystalline NbMoTaW and Nb 45 Ta 25 Ti 15 Hf 15 , each with two different grain sizes, revealing trends consistent with experimental observations and the extended Rice model. Comparison with averaged R-curves for bulk samples demonstrates that GBs enhance ductility in Nb 45 Ta 25 Ti 15 Hf 15 in both grain sizes investigated. In contrast, GBs only locally improve fracture resistance in NbMoTaW when cracks are temporarily pinned at GBs inclined at high angles from the crack, but generally promote brittle intergranular fracture. These contrasting behaviors are attributed to differences in GB cohesion, reflecting clear alloying trends that align with ab-initio calculations and trends observed experimentally. Our results bridge classical fracture theory, atomistic simulations, and experimental observations, providing a comprehensive understanding of the fracture mechanisms in nanocrystalline refractory complex concentrated alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Molecular Tuning of Ether Cosolvent Chemistry for High-Voltage Sodium-Ion Batteries

Ethers as electrolyte cosolvents in sodium-ion batteries (SIBs) provide favorable Na + solvation and interfacial properties, but their low oxidative stability limits their use in high-voltage SIBs. Herein, we address this limitation via molecular tuning of ether cosolvents for high-voltage (4.2 V) hard carbon || NaNi 0.33 Fe 0.33 Mn 0.33 O 2 full cells. Tetrahydropyran (THP) is functionalized with a nitrile group to form tetrahydropyran-4-carbonitrile (THPCN). To delineate the effect of nitrile functionalization and benchmark ether against a conventional carbonate, THP, THPCN, and diethyl carbonate (DEC) are evaluated as cosolvents with ethylene carbonate. Nitrile functionalization lowers the HOMO energy of the ether, extends the electrolyte stability window, and alters Na⁺ solvation. Spectroscopic techniques and molecular dynamics simulations reveal that THPCN exhibits predominantly aggregate-dominated solvation (95.1 %) with weakened Na + -solvent interactions, producing the most anion-rich environment relative to DEC and THP cosolvents. THPCN-modified solvation promotes the formation of highly conductive, fluorine-enriched interphases that suppress parasitic reactions. Pouch full cells with THPCN sustained ~ 600 cycles at 4.2 V, outperforming THP and DEC. Operando gas analysis reveals that THPCN reduces CO 2 generation by 45% and H 2 generation by 30% relative to THP. Furthermore, the findings demonstrate nitrile functionalization as a molecular design strategy to stabilize ethers and enable high-voltage SIBs.

25 ENERGY STORAGE↗
DOE OSTI2021

Effect of Helium Flux on Near-surface Helium Accumulation in Plasma-Exposed Tungsten

We report results of object kinetic Monte Carlo (OKMC) simulations aimed at understanding the effect of helium flux on the near-surface helium accumulation in plasma-facing tungsten, which is initially defect-free and has a W(100) surface orientation. These OKMC simulations are performed at 933 K for fluxes ranging from 10 22 to 4 × 10 25 He/m 2 s, with 100 eV helium atoms impinging on a W(100) surface up to a maximum fluence of 4×10 19 He/m 2 . In the near-surface region, helium clusters interact elastically with the free surface. The interaction is attractive and results in the drift of mobile helium clusters towards the surface as well as increased trap mutation rates. The associated kinetics and energetics of the above-mentioned processes obtained from molecular dynamics simulations are also considered. The OKMC simulations indicate that as the flux decreases, the retention of implanted helium decreases, and its depth distribution shifts to deeper below the surface in initially pristine tungsten. Furthermore, the fraction of retained helium diffusing into the bulk increases as well, so much so that for 10 22 He/m 2 s, almost all of the retained helium diffused into the bulk with minimal/negligible near-surface helium accumulation. At a given flux, with increasing fluence, the fraction of retained helium initially decreases and then starts to increase after reaching a minimum. The occurrence of the retention minimum shifts to higher fluences as the flux decreases. Although the near-surface helium accumulation spreads deeper into the material with decreasing flux and increasing fluence, the spread appears to saturate at depths between 80 and 100 nm. Finally, we present a detailed analysis of the influence of helium flux on the size and depth distribution of total helium and helium bubbles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Theory and Simulation of Nanoscale Phenomena

The Center for Integrated Nanotechnologies (CINT) is a Department of Energy/Office of Science Nanoscale Science Research Center (NSRC), operating as a national user facility devoted to establishing the scientific principles that govern nanoscale integration. Nanoscale integration is defined as assembling diverse nanoscale materials across length scales to design and achieve new properties and functionality. The CINT Theory and Simulation of Nanoscale Phenomena thrust is the component of CINT dedicated to developing and applying theory to enable nanoscale integration. Our focus is on understanding and simulating the unique behavior of integrated materials and systems with nanoscale structure. This mission is achieved through collaborations with CINT Users, between thrust scientists, and with CINT scientists from other thrusts. Our research is focused on three science directions that together form the basis for integration at the nanoscale, namely (i) Hierarchical structure and dynamics in soft matter, (ii) Excitation and Transport in Nanostructured Systems, and (iii) Emergent phenomena at surfaces and interfaces. A broad spectrum of techniques is developed and applied including continuum fluid theory, atomistic and coarse-grained molecular dynamics simulations, static and dynamic electronic structure calculations, multiscale modeling, low-energy effective Hamiltonian methods, and perturbative and exact quantum many-body approaches. These tools are applied to physical systems of interest to CINT Users, the other CINT thrusts, and the general scientific community with the goals of understanding and controlling the interactions between nanoscale building blocks to assemble specific integrated structures, controlling energy transfer and other interactions over multiple length scales, and designing and exploiting the interactions within assembled structures to achieve new materials functionality.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2022

Bridging confined phase behavior of CH 4 -CO 2 binary systems across scales

Phase behavior of confined fluids may deviate significantly from that of the bulk fluid due to the fluid-wall interactions being a significant portion of all intermolecular interactions under confinement. Despite recent advancements in understanding confined phase behavior of pure fluids, confined phase behavior of mixtures remains an understudied topic. In this work, we examine the confined phase behavior of a CH 4 -CO 2 binary system by combining Monte Carlo (MC) simulations, a cubic equation of state (EoS), and the lattice Boltzmann method (LBM). First, the effects of confinement on density and phase distribution in nano-size pores are established using Gibbs Ensemble MC calculations, which produce precise results of liquid and vapor confined pressures and account for the modification of the phase change location. By comparing the phase envelopes of bulk and confined mixtures at a fixed temperature, here it is observed that the phase envelopes shrink with reductions in pore size. Based on this observation, we extend a modified Peng-Robinson EoS, which was originally developed for pure fluids under confinement, to mixtures via van-der-Waals-type mixing rules and by accounting for shifts in the critical properties of confined CH 4 -CO 2 . The resulting phase envelopes are in good agreement with the MC data. In addition, a local density model is used in combination with the confined EoS to calculate adsorption isotherms of CH 4 -CO 2 mixtures and to characterize the behavior of confined matter in nanopores. Finally, we incorporate this EoS in a multicomponent multiphase LBM that uses a pseudopotential model to represent intermolecular forces. This workflow utilizes multiscale simulation techniques to bridge the behavior of multicomponent systems across scales and to shed light on the confined phase behavior of CH 4 -CO 2 binary systems.

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

Molecular Design of Functional Polymers for Silica Scale Inhibition

Silica polymerization, which involves the condensation reaction of silicic acid, is a fundamental process with wide-ranging implications in biological systems, material synthesis, and scale formation. The formation of a silica-based scale poses significant technological challenges to energy-efficient operations in various industrial processes, including heat exchangers and water treatment membranes. Despite the common strategy of applying functional polymers for inhibiting silica polymerization, the underlying mechanisms of inhibition remain elusive. In this study, we synthesized a series of nitrogen-containing polymers as silica inhibitors and elucidated the role of their molecular structures in stabilizing silicic acids. Polymers with both charged amine and uncharged amide groups in their backbones exhibit superior inhibition performance, retaining up to 430 ppm of reactive silica intact for 8 h under neutral pH conditions. In contrast, monomers of these amine/amide-containing polymers as well as polymers containing only amine or amide functionalities present insignificant inhibition. Molecular dynamics simulations reveal strong binding between the deprotonated silicic acid and a polymer when the amine groups in the polymer are protonated. Notably, an extended chain conformation of the polymer is crucial to prevent proximity between the interacting monomeric silica species, thereby facilitating effective silica inhibition. Furthermore, the hydrophobic nature of alkyl segments in polymer chains disrupts the hydration shell around the polymer, resulting in enhanced binding with ionized silicic acid precursors compared to monomers. Finally, our findings provide novel mechanistic insights into the stabilization of silicic acids with functional polymers, highlighting the molecular design principles of effective inhibitors for silica polymerization.

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

Asymmetric Composition of Ionic Aggregates and the Origin of High Correlated Transference Number in Water-in-Salt Electrolytes

"Water-in-salt" electrolytes open up exciting new avenues for expanding the electrochemical window of aqueous electrolytes. We investigated the solvation structure and dynamics of highly concentrated lithium bis(trifluoromethane)sulfonimide aqueous electrolyte using experimentally corroborated molecular dynamics simulations. The simulations revealed that the heterogeneous structure of the electrolyte comprises percolating networks of ion and water domains/aggregates. Interestingly, the ionic regions are composed of more TFSI - ions than Li + ions. The Li + -ion transport mechanism was further explored. Li + ions can hop along the coordinated TFSI- ions in the ionic aggregates. The calculated correlated transference number of the 20 m electrolyte is similar to 0.32, which is reasonably high for the high concentration due to a weak negative correlation between the motion of cations and anions within the heterogeneous microscopic domains. These molecular dynamics results connect the heterogeneous structure of the electrolyte to the correlated dynamics of the Li + ion and provide a new understanding of the Li + -ion transport mechanism in this novel electrolyte.

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

Time-Evolving Chirality Loss in Molecular Photodissociation Monitored by X-ray Circular Dichroism Spectroscopy

In this work, the ultrafast photoinduced chirality loss of 2-iodobutane is studied theoretically by time- and frequency-resolved X-ray circular dichroism (TRXCD) spectroscopy. Following an optical excitation, the iodine atom dissociates from the chiral center, which we capture by quantum non-adiabatic molecular dynamics simulations. At variable time delays after the pump, the resonant X-ray pulse selectively probes the iodine and carbon atom involved in the chiral dissociation through a selected core-to-valence transition. The TRXCD signal at the iodine L 1 edge accurately captures the timing of C–I photodissociation and thereby chirality loss, c.a 70 fs. The strong electric dipole–electric quadrupole (ED–EQ) response makes this signal particularly sensitive to vibronic coherence at the high X-ray regime. At the carbon K-edges, the signals monitor the molecular chirality of the 2-butyl radical photoproduct and the spin state of the iodine atom. The ED–EQ response is masked under the strong electric dipole-magnetic dipole response, making this signal intuitive for the electronic population. The evolution of the core electronic states and its chiral sensitivity is discussed. Overall, the element-specific TRXCD signal provides a detailed picture of molecular dynamics and offers a unique sensitive window into the time-dependent chirality of molecules.

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

Uncertainty-driven dynamics for active learning of interatomic potentials

Machine learning (ML) models, if trained to data sets of high-fidelity quantum simulations, produce accurate and efficient interatomic potentials. Active learning (AL) is a powerful tool to iteratively generate diverse data sets. In this approach, the ML model provides an uncertainty estimate along with its prediction for each new atomic configuration. If the uncertainty estimate passes a certain threshold, then the configuration is included in the data set. Here we develop a strategy to more rapidly discover configurations that meaningfully augment the training data set. The approach, uncertainty-driven dynamics for active learning (UDD-AL), modifies the potential energy surface used in molecular dynamics simulations to favor regions of configuration space for which there is large model uncertainty. The performance of UDD-AL is demonstrated for two AL tasks: sampling the conformational space of glycine and sampling the promotion of proton transfer in acetylacetone. The method is shown to efficiently explore the chemically relevant configuration space, which may be inaccessible using regular dynamical sampling at target temperature conditions.

97 MATHEMATICS AND COMPUTING↗