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,477 records · Page 82
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Elevating density functional theory to chemical accuracy for water simulations through a density-corrected many-body formalism

Abstract Density functional theory (DFT) has been extensively used to model the properties of water. Albeit maintaining a good balance between accuracy and efficiency, no density functional has so far achieved the degree of accuracy necessary to correctly predict the properties of water across the entire phase diagram. Here, we present density-corrected SCAN (DC-SCAN) calculations for water which, minimizing density-driven errors, elevate the accuracy of the SCAN functional to that of “gold standard” coupled-cluster theory. Building upon the accuracy of DC-SCAN within a many-body formalism, we introduce a data-driven many-body potential energy function, MB-SCAN(DC), that quantitatively reproduces coupled cluster reference values for interaction, binding, and individual many-body energies of water clusters. Importantly, molecular dynamics simulations carried out with MB-SCAN(DC) also reproduce the properties of liquid water, which thus demonstrates that MB-SCAN(DC) is effectively the first DFT-based model that correctly describes water from the gas to the liquid phase.

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

Amplifying Nanoparticle Reinforcement through Low Volume Topologically Controlled Chemical Coupling

Here, we present a streamlined method to covalently bond hydroxylated carbon nanotubes (CNOH) within a polyphenol matrix, all achieved through a direct, solvent-free process. Employing an extremely small concentration of CNOH (0.01% w/w) along with topologically contrasting linkers led to a maximum of 5-fold increase in modulus and a 25% enhancement in tensile strength compared to the unaltered matrix, an order of magnitude greater reinforcement (w/w) compared to state-of-the-art melt-processed nanocomposites. Through dynamic mechanical analysis, low field solid-state nuclear magnetic resonance spectroscopy, and molecular dynamics simulations, we uncovered the profound influence of linker’s conformational degrees of freedom on the segmental dynamics and therefore the material’s properties.

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

Stereoselective recognition of morphine enantiomers by μ -opioid receptor

Stereospecific recognition of chiral molecules plays a crucial role in biological systems. The μ-opioid receptor (MOR) exhibits binding affinity towards (-)-morphine, a well-established gold standard in pain management, while it shows minimal binding affinity for the (+)-morphine enantiomer, resulting in a lack of analgesic activity. Understanding how MOR stereoselectively recognizes morphine enantiomers has remained a puzzle in neuroscience and pharmacology for over half-a-century due to the lack of direct observation techniques. To unravel this mystery, we constructed the binding and unbinding processes of morphine enantiomers with MOR via molecular dynamics simulations to investigate the thermodynamics and kinetics governing MOR's stereoselective recognition of morphine enantiomers. Our findings reveal that the binding of (-)-morphine stabilizes MOR in its activated state, exhibiting a deep energy well and a prolonged residence time. In contrast, (+)-morphine fails to sustain the activation state of MOR. Furthermore, the results suggest that specific residues, namely D114 2.50 and D147 3.32 , are deprotonated in the active state of MOR bound to (-)-morphine. This work highlights that the selectivity in molecular recognition goes beyond binding affinities, extending into the realm of residence time.

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

Using limited neural networks to assess relative mechanistic influence on shock heating in granular solids

The rapid compaction of granular media results in localized heating that can induce chemical reactions, phase transformations, and melting. However, there are numerous mechanisms in play that can be dependent on a variety of microstructural features. Machine learning techniques such as neural networks offer a ubiquitous method to develop models for physical processes. Limiting what kind of microstructural information is used as an input and assessing normalized changes in network error, the relative importance of different mechanisms can be inferred. Here we utilize binned, initial density information as network inputs to predict local shock heating in a granular high explosive trained from large-scale molecular dynamics simulations. In this study, the spatial extent of the density field used in the network is altered to assess the importance and relevant length scales of the physical mechanisms in play, where different microstructural features result in different predictive capabilities.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Simulations of classical three-body thermalization in one dimension

One-dimensional systems, such as nanowires or electrons moving along strong magnetic field lines, have peculiar thermalization physics. The binary collision of pointlike particles, typically the dominant process for reaching thermal equilibrium in higher-dimensional systems, cannot thermalize a 1D system. We study how dilute classical 1D gases thermalize through three-body collisions. We consider a system of identical classical point particles with pairwise repulsive inverse power-law potential V ij ∝ 1/|x i –x j | n or the pairwise Lennard-Jones potential. Using Monte Carlo methods, we compute a collision kernel and use it in the Boltzmann equation to evolve a perturbed thermal state with temperature T toward equilibrium. We explain the shape of the kernel and its dependence on the system parameters. Additionally, we implement molecular dynamics simulations of a many-body gas and show agreement with the Boltzmann evolution in the low-density limit. For the inverse power-law potential, the rate of thermalization is proportional to ρ 2 ⁢T$\frac{1}{2}$ – $\frac{1}{n}$, where ρ is the number density. Furthermore, the corresponding proportionality constant decreases with increasing n.

1-dimensional systems↗
DOE OSTI2024

Phase Diagram of Nuclear Pastas in Neutron Star Crusts

Two neural networks were trained to predict, respectively, the Euler characteristic and the curvature of nuclear pastas in neutron star crust conditions generated by molecular dynamics simulations of neutron star matter with 0.1 < x < 0.5, 0.040 fm−3 < ρ < 0.085 fm−3 (0.68 × 1014 g/cm3 < ρ < 1.43 × 1014 g/cm3), and 0.2 MeV < T < 4.0 MeV, where x is proton content, the density is ρ, and the temperature is T. The predictions of the two networks were combined to determine the nuclear pasta phase that is thermodynamically stable at a given x, ρ, and T, and a three-dimensional phase diagram that extrapolated slightly the regions of existing molecular dynamics data was computed. The jungle gym and anti-jungle gym structures are prevalent at high temperature and low density, while the anti-jungle gym and anti-gnocchi structures dominate at high temperature and high density. A diversity of structures exist at low temperatures and intermediate density and proton content. The trained models used in this work are open access and available at a public repository to promote comparison to pastas obtained with other models.

Muñoz, Jorge A. (ORCID:0000000294870079)↗
DOE OSTI2021

Crystal Structure and Preparation of Li7La3Zr2O12 (LLZO) Solid-State Electrolyte and Doping Impacts on the Conductivity: An Overview

As an essential part of solid-state lithium-ion batteries, solid electrolytes are receiving increasing interest. Among all solid electrolytes, garnet-type Li7La3Zr2O12 (LLZO) has proven to be one of the most promising electrolytes because of its high ionic conductivity at room temperature, low activation energy, good chemical and electrochemical stability, and wide potential window. Since the first report of LLZO, extensive research has been done in both experimental investigations and theoretical simulations aiming to improve its performance and make LLZO a feasible solid electrolyte. These include developing different methods for the synthesis of LLZO, using different crucibles and different sintering temperatures to stabilize the crystal structure, and adopting different methods of cation doping to achieve more stable LLZO with a higher ionic conductivity and lower activation energy. It also includes intensive efforts made to reveal the mechanism of Li ion movement and understand its determination of the ionic conductivity of the material through molecular dynamic simulations. Nonetheless, more insightful study is expected in order to obtain LLZO with a higher ionic conductivity at room temperature and further improve chemical and electrochemical stability, while optimal multiple doping is thought to be a feasible and promising route. This review summarizes recent progress in the investigations of crystal structure and preparation of LLZO, and the impacts of doping on the lithium ionic conductivity of LLZO.

Raju, Md Mozammal (ORCID:0000000319040523)↗
DOE OSTI2025

Equilibrium Separation of Siloxanes in Metal–Organic Frameworks

We present an in silico assessment of metal–organic frameworks (MOFs) for the equilibrium separation of linear and cyclic siloxanes. Using a combination of configurational bias/continuous fractional component Monte Carlo (CB/CFCMC) simulations and Ideal Adsorbed Solution Theory (IAST), we investigated the adsorption of both equimolar and nonequimolar mixtures of linear and cyclic siloxanes in a selection of synthesizable MOFs with medium to large pore volumes. We showed that configurational entropy effects drive the preferential adsorption of linear siloxanes over cyclic siloxanes. Based on synthesizability metrics we identified ZIF-70 as a promising adsorbent for the separation of linear and cyclic siloxanes. Self-diffusivities of linear and cyclic siloxanes in ZIF-70 calculated from molecular dynamics simulations show that equilibrium can be reached on reasonable time scales. We explored vacuum-temperature swing adsorption (VTSA) as a potential method for the recovery of adsorbed linear siloxanes from ZIF-70, achieving a 75% reduction in adsorbed phase concentration. Additionally, we demonstrated that supercritical CO 2 offers an alternative desorption strategy by displacing adsorbed linear siloxanes in ZIF-70 at pressures above 200 bar driven by entropy effects.

Adsorption↗
DOE OSTI2022

Dissociating the phononic, magnetic and electronic contributions to thermal conductivity: a computational study in alpha-iron

Computational tools to study thermodynamic properties of magnetic materials have, until recently, been limited to phenomenological modeling or to small domain sizes limiting our mechanistic understanding of thermal transport in ferromagnets. Herein, we study the interplay of phonon and magnetic spin contributions to the thermal conductivity in a-iron utilizing non-equilibrium molecular dynamics simulations. It was observed that the magnetic spin contribution to the total thermal conductivity exceeds lattice transport for temperatures up to two-thirds of the Curie temperature after which only strongly coupled magnon-phonon modes become active heat carriers. Characterizations of the phonon and magnon spectra give a detailed insight into the coupling between these heat carriers, and the temperature sensitivity of these coupled systems. Comparisons to both experiments and ab initio data support our inferred electronic thermal conductivity, supporting the coupled molecular dynamics/spin dynamics framework as a viable method to extend the predictive capability for magnetic material properties.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Investigation of Low-Energy Lattice Dynamics and Their Role in Superionic Na Diffusion and Ultralow Thermal Conductivity of Na 3 PSe 4 as a Solid-State Electrolyte

The atomic dynamics of Na 3 PSe 4 were investigated using a combination of neutron scattering experiments and ab initio and machine-learned molecular dynamics simulations to probe the interplay of fast ionic diffusion with atomic vibrations (phonons) of the host lattice. Our results reveal the existence of low-energy vibrational modes, simultaneously involving motions of Na + ions and framework polyanion subunits, and show that these modes become strongly overdamped in the superionic regime as they couple with the Na + hopping process. In particular, the Na+ migration energy landscape is strongly impacted by low-energy phonons derived from a soft acoustic branch of the host lattice, which modulates the diameter of the Na + diffusion channel at the bottleneck. We find that an additional factor for the enhanced Na + conductivity in Na 3 PSe 4 is the presence of Na-vacancies, which also affect the low-frequency dynamics and thermal vibration amplitudes, pointing to an interplay between Na + vacancies and host dynamics, jointly enhancing ionic diffusivity. Finally, we investigate the origin of ultralow thermal conductivities in Na 3 PSe 4 and Na 3 PS 4 using Green–Kubo simulations and find that low-energy acoustic phonon modes of the overall crystal framework provide a dominant contribution to the thermal conductivity.

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

Surface Hydroxyls of Imogolite Nanotubes Drive Distinct Structures and Mobility Differentiation of Nanoconfined Water

Abstract Nanoconfined fluids, particularly water, govern subsurface geochemistry, yet the molecular-level mechanisms by which mineral surfaces dictate confined water structure and mobility remain poorly resolved. Here, we combine solution- and solid-state proton (1H) NMR spectroscopy, NMR relaxometry, modulated-gradient spin–echo (MGSE) NMR diffusometry, infrared spectroscopy, and molecular dynamics simulations to demonstrate that surface hydroxyls drive the structural and dynamic differentiation of water in imogolite nanotubes. In saturated suspensions, we observe the coexistence of distinct water 1H environments, each exhibiting significantly reduced mobility, which we attribute to strong interactions with the imogolite surfaces. As more mobile water is removed, long-range water diffusivity in the fibrous solid samples slows to 1.6 × 10–10 m2·s–1 while local fluctuations increase to 3.4 × 10–8 m2·s–1, indicating a significant increase in molecular restriction through surface interactions. Together, our experiments reveal three coexisting populations with varied structures and mobilities, all distinct from liquid bulk water. We resolve a persistent solid-like interfacial layer strongly bound to surface hydroxyls, characterized by an unusually short T2 (<1 ms), and a very close effective 1H–1H distance of ∼1.55 Å between water and the inner-surface silanol group. An inner-core population occupying the inner cavity exhibits an intermediate dynamic regime with restricted axial diffusion, while outer-surface water associated with aluminum hydroxyls retains relatively higher mobility. These results experimentally substantiate a hierarchy of water populations in imogolite and show how surface chemistry dictates the structures and dynamics of confined water.

Fleming, Xander B. [Earth and Environmental Scienc↗
DOE OSTI2022

Argon equation of state data to 1 TPa: Shock compression experiments and simulations

Argon is the most abundant noble gas on Earth and its noble, atomic fluid nature makes it an excellent candidate for comparison of experiment and theory at extreme conditions. For this study, we performed a combined computational and experimental study on shock compressed cryogenic liquid argon. Using Sandia's Z machine, we shock compressed liquid argon to 600 GPa and reshock states up to 950 GPa. Laser shock experiments at the Omega Laser facility extend the principal Hugoniot to 1000 GPa and provided temperature data along the principal Hugoniot. The plate impact experiments and laser shock experiments used well-characterized impedance matching standards and demonstrate consistent results between the two platforms over a common range. Density functional theory based molecular dynamics simulations provided additional data on the Hugoniot to 600 GPa. The combined experimental data and simulation results provide constraints on the development of new equation of state models at extreme conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Modular Approach for the Synthesis of Bottlebrush Diblock Copolymers from Poly(Glycidyl Methacrylate)-block-Poly(Vinyldimethylazlactone) Backbones

A strategy that enables the facile synthesis of bottlebrush block copolymers with flexible backbones was developed. A demonstration of the strategy’s utility was carried out by grafting end-functionalized polymethylmethacrylate (PMMA) and polystyrene (PS) to the dually reactive block copolymer, poly(glycidyl methacrylate)-block-poly(vinyldimethylazlactone) (PGMA-b-PVDMA). Five different bottlebrush diblock copolymers were investigated by size-exclusion chromatography (SEC), 1H NMR, Fourier transform infrared (FT-IR), differential scanning calorimetry (DSC), X-ray scattering methods, atomic force microscopy (AFM), rheology and computational simulations using molecular dynamics (MD), and self-consistent field theory (SCFT). A relationship between the glass transition temperature and the fraction of chain ends was demonstrated by DSC and highlights the potential of this synthetic method to tailor polymer properties. All five samples were found to be in a disordered phase exhibiting multiscale structures revealed by two broad peaks in small-angle X-ray scattering (SAXS) that can be attributed to graft-to-graft and backbone-to-backbone density correlations using MD simulations. The SCFT-based simulations justify the observation of a disordered phase due to its stabilization by the grafts. Additionally, this modular approach can be easily extended to other grafts, including responsive, conducting, and charged polymers with the prerequisite end groups. The versatility and ease of assembling these functional bottlebrushes constitute a powerful “toolbox” method for the rapid and scalable synthesis of novel bottlebrush block copolymers with desired properties.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Substrate binding-induced conformational transitions in the omega-3 fatty acid transporter MFSD2A

Major Facilitator Superfamily Domain containing 2 A (MFSD2A) is a transporter that is highly enriched at the blood-brain and blood-retinal barriers, where it mediates Na + -dependent uptake of ω-3 fatty acids in the form of lysolipids into the brain and eyes, respectively. Despite recent structural insights, it remains unclear how this process is initiated, and driven by Na + . Here, we perform Molecular Dynamics simulations which demonstrate that substrates enter outward facing MFSD2A from the outer leaflet of the membrane via lateral openings between transmembrane helices 5/8 and 2/11. The substrate headgroup enters first and engages in Na + -bridged interactions with a conserved glutamic acid, while the tail is surrounded by hydrophobic residues. This binding mode is consistent with a “trap-and-flip” mechanism and triggers transition to an occluded conformation. Furthermore, using machine learning analysis, we identify key elements that enable these transitions. These results advance our molecular understanding of the MFSD2A transport cycle.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Gold–Thiolate Nanocluster Dynamics and Intercluster Reactions Enabled by a Machine Learned Interatomic Potential

Monolayer protected metal clusters comprise a rich class of molecular systems and are promising candidate materials for a variety of applications. While a growing number of protected nanoclusters have been synthesized and characterized in crystalline forms, their dynamical behavior in solution, including prenucleation cluster formation, is not well understood due to limitations both in characterization and first-principles modeling techniques. Recent advancements in machine-learned interatomic potentials are rapidly enabling the study of complex interactions such as dynamical behavior and reactivity on the nanoscale. Here, we develop an Au-S-C-H atomic cluster expansion (ACE) interatomic potential for efficient and accurate molecular dynamics simulations of thiolate-protected gold nanoclusters (Au n (SCH 3 ) m ). Trained on more than 30,000 density functional theory calculations of gold nanoclusters, the interatomic potential exhibits ab initio level accuracy in energies and forces and replicates nanocluster dynamics including thermal vibration and chiral inversion. Long dynamics simulations (up to 0.1 μs time scale) reveal a mechanism explaining the thermal instability of neutral Au 25 (SR) 18 clusters. Specifically, we observe multiple stages of isomerization of the Au 25 (SR) 18 cluster, including a chiral isomer. Additionally, we simulate coalescence of two Au 25 (SR) 18 clusters and observe series of clusters where the formation mechanisms are critically mediated by ligand exchange in the form of [Au-S] n rings.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗