Search NASA⌕ Search

SEARCH · Search NASA

Results for “Molecular Dynamic simulations”

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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Molecular dynamics of interface rupture

Several situations have been studied in which a fluid-vapor or fluid-fluid interface ruptures, using molecular dynamics simulations of 3000 to 20,000 Lennard-Jones molecules in three dimensions. The cases studied are the Rayleigh instability of a liquid thread, the burst of a liquid drop immersed in a second liquid undergoing shear, and the rupture of a liquid sheet in an extensional flow. The late stages of the rupture process involve the gradual withdrawal of molecules from a thinning neck, or the appearance and growth of holes in a sheet. In all cases, it is found that despite the small size of the systems studied, tens of angstroms, the dynamics is in at least qualitative accord with the behavior expected from continuum calculations, and in some cases the agreement is to within tens of percent. Remarkably, this agreement occurs even though the Eulerian velocity and stress fields are essentially unmeasurable - dominated by thermal noise. The limitations and prospects for such molecular simulation techniques are assessed.

Koplik, Joel↗
DOE OSTI2025

Solvated Electrons Have Multiple Personalities in Molten Salts

P-Crosscut: Solvated Electrons Have Multiple Personalities in Molten Salts [EFRC – MSEE] Alejandro Ramos-Ballesteros;1 Hung H. Nguyen;2 Kazuhiro Iwamatsu;3 Santanu Roy;4 Vyacheslav Bryantsev;4 Michael E. Woods;1 Ruchi Gakhar;1 Phillip Halstenberg;4 Bobby Layne;5 Jay A. LaVerne;6 Claudio J. Margulis:2* and James F. Wishart5* 1Idaho National Laboratory; 2The University of Iowa; 3Hunter College; 4Oak Ridge National Laboratory; 5Brookhaven National Laboratory; 6University of Notre Dame Abstract: The solvated (eS–) is a powerful reducing agent and one of the primary products of molten salt radiolysis. In these extreme high-temperature environments, the eS– can initiate cascades of redox processes that significantly alter the physical and chemical properties of a molten salt, posing challenges for molten salt reactor (MSR) performance and longevity. Consequently, mastering the fundamental behavior of the eS– could enable the design of specific molten salt mixtures with tailored Lewis acidities for controlling the speciation and chemical reactivity of the eS–, thereby mitigating its overall impact on MSR technologies. Here, we use time-resolved electron pulse radiolysis techniques for determining chemical kinetics and transient absorption spectra, combined with ab initio molecular dynamics simulations to explore the influence of multivalent metal cations on the fundamental speciation and distribution of eS– coordination environments.

38 - RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCL↗
DOE OSTI2023

Molecular pathway and structural mechanism of human oncochannel TRPV6 inhibition by the phytocannabinoid tetrahydrocannabivarin

The calcium-selective oncochannel TRPV6 is an important driver of cell proliferation in human cancers. Despite increasing interest of pharmacological research in developing synthetic inhibitors of TRPV6, natural compounds acting at this channel have been largely neglected. On the other hand, pharmacokinetics of natural small-molecule antagonists optimized by nature throughout evolution endows these compounds with a medicinal potential to serve as potent and safe next-generation anti-cancer drugs. Here we report the structure of human TRPV6 in complex with tetrahydrocannabivarin (THCV), a natural cannabinoid inhibitor extracted from Cannabis sativa . We use cryo-electron microscopy combined with electrophysiology, calcium imaging, mutagenesis, and molecular dynamics simulations to identify THCV binding sites in the portals that connect the membrane environment surrounding the protein to the central cavity of the channel pore and to characterize the allosteric mechanism of TRPV6 inhibition. We also propose the molecular pathway taken by THCV to reach its binding site. Our study provides a foundation for the development of new TRPV6-targeting drugs.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Comment on “Brownian motion of droplets induced by thermal noise”

We simulate phase separated fluids using the Cahn-Hillard fluctuating hydrodynamic (CH-FHD) model and measure the statistical properties of capillary waves generated by thermal fluctuations. Our measurements are in good agreement with stochastic lubrication theory and molecular dynamics simulations but differ significantly from recent CH-FHD results by Zhang et al. [Phys. Rev. E 109, 024208 (2024)2470-004510.1103/PhysRevE.109.024208]. Specifically, we find that capillary wave statistics at thermodynamic equilibrium are independent of transport properties, namely viscosity and species diffusion.

Bell, J B↗
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 OSTI2020

Mapping the Interfacial Chemistry and Structure of Partially Fluorinated Bottlebrush Polymers and Their Linear Analogues

Polymer interfaces are key to a range of applications including membranes for chemical separations, hydrophobic coatings, and passivating layers for antifouling. While important, challenges remain in probing the interfacial monolayer where the molecular ordering and orientation can change depending on the chemical makeup or processing conditions. In this work, we leverage surface specific vibrational sum frequency generation (SFG) and the associated dependence on molecular symmetry to elucidate the ordering and orientations of key functional groups for poly(2,2,2-trifluoroethyl methacrylate) bottlebrush polymers and their linear polymer analogues. These measurements were framed by atomistic molecular dynamic simulations to provide a complementary physical picture of the gas–polymer interface. Simulations and SFG measurements show that methacrylate backbones are buried beneath a layer of trifluoroethyl containing side groups that result in structurally similar interfaces regardless of the polymer molecular weight or architecture. The average orientational angles of the trifluoroethyl containing side groups differ depending on polymer linear and bottlebrush architectures, suggesting that the surface groups can reorient via available rotational degrees of freedom. Results show that the surfaces of the bottlebrush and linear polymer samples do not strongly depend on molecular weight or architecture. As such, one cannot rely on increasing the molecular weight or altering the architecture to tune surface properties. This insight into the polymer interfacial structure is expected to advance the design of new material interfaces with tailored chemical/functional properties.

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

Probing the state of hydrogen in δ − AlOOH at mantle conditions with machine learning potential

Hydrous and nominally anhydrous minerals are a fundamental class of solids of enormous significance to geophysics. They are the water carriers in the deep geological water cycle and impact structural, elastic, plastic, and thermodynamic properties and phase relations in Earth's forming aggregates (rocks). They play a critical role in the geochemical and geophysical processes that shape the planet. Their complexity has prevented predictive calculations of their properties, but progress in materials simulations ushered by machine-learning potentials is transforming this state of affairs. Here, we adopt a hybrid approach that combines deep learning potentials (DPs) with the strongly constrained and appropriately normed meta-generalized gradient approximation functional to simulate a prototypical hydrous system. We illustrate the success of this approach to simulate δ − AlOOH ( δ ) , a phase capable of transporting water down to near the core-mantle boundary of the Earth ( ∼ 2900 km depth and ∼ 135 GPa ) in subducting slabs. A high-throughput sampling of phase space using molecular dynamics simulations with DPs sheds light on the hydrogen-bond behavior and proton diffusion at geophysical conditions. These simulations provide a pathway for a deeper understanding of these crucial components that shape Earth's internal state. Published by the American Physical Society 2024

Luo, Chenxing (ORCID:0000000341166851)↗
DOE OSTI2022

Concentration Dependent Interfacial Chemistry of the NaOH(aq):Gibbsite Interface

Caustic conditions are often employed for dissolution of a wide variety of minerals, where ion sorption, surface diffusion, and interfacial organization impact surface reactivity. In the case of gibbsite, gamma-Al(OH)3, the chemistry at the NaOH(aq) interface is deeply intertwined with industrial processing of aluminum, including metal production and the disposition of Al-containing wastes. To date, little is known about the structure, speciation, and dynamic behavior of gibbsite interfaces (and that of many other minerals) with NaOH(aq)—particularly as a function of ionic strength. Yet concentration-dependent interfacial organization and dynamics are a critical starting point to develop a fundamental understanding of the factors that influence dissolution. This work reports equilibrium molecular dynamics simulations of the gamma-Al(OH)3:NaOH(aq) interface, revealing the sorption behavior and speciation of ions from 0.5–10 M [NaOH]. As inner-sphere complexes, Na+ primarily coordinates to the side of the gibbsite hexagonal cavities, while OH accepts hydrogen-bonding from the surface-OH groups. The mobility of inner-sphere Na+ and OH ions is significantly reduced due to a strong surface affinity in comparison to previous reports of NaCl, CaCl2, or BaCl2 electrolytes. At high [NaOH], contact ion pairing that is observed in the bulk solution is partially disrupted upon sorption to the gibbsite surface by the individual ion–surface interactions. The molecular-scale changes to surface speciation and competition between ion–surface vs. ion–ion interactions influence surface characterization of gibbsite and potential dissolution processes, providing a valuable baseline for starting conditions needed within future reactive molecular simulations.

alkaline solutions, IDREAM, gibbsite, molecular dy↗
DOE OSTI2023

Disentangling the effects of non-adiabatic interactions upon ion self-diffusion within warm dense hydrogen

Warm dense matter is a material state in the region of parameter space connecting condensed matter to classical plasma physics. In this intermediate regime, we investigate the significance of non-adiabatic electron-ion interactions upon ion dynamics. To disentangle non-adiabatic from adiabatic electron-ion interactions, we compare the ion self-diffusion coefficient from the non-adiabatic electron force field computational model with an adiabatic, classical molecular dynamics simulation. A classical pair potential developed through a force-matching algorithm ensures the only difference between the models is due to the electronic inertia. Here, we implement this new method to characterize non-adiabatic effects on the self-diffusion of warm dense hydrogen over a wide range of temperatures and densities. Ultimately we show that the impact of non-adiabatic effects is negligible for equilibrium ion dynamics in warm dense hydrogen.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Isotropic-to-Nematic Transition in Salt-Free Polyelectrolyte Coacervates from Coarse-Grained Simulations

Recent interest in complex coacervation between oppositely charged polyelectrolytes (PEs) has been fueled by its relevance to biology in the context of membraneless organelle formation within living cells. For PEs with limited flexibility (such as double-stranded DNA), theoretical treatments and recent experiments have reported the emergence of liquid crystalline order (LCO) within the resulting coacervate phases. In this work, we study the underlying physics of this phenomenon using coarse-grained molecular dynamics simulations of symmetric semiflexible-semiflexible and asymmetric semiflexible-flexible coacervates. By comparing coacervates with the corresponding semidilute solutions of neutral polymers, we demonstrate that the presence of Coulomb interactions in coacervates facilitates orientational ordering, in agreement with theoretical predictions. Quantitative comparisons between our simulations and theory indicate that, for asymmetric nematic coacervates, the strong orientational ordering of stiff polyanions induces a weak ordering of the flexible polycations-an effect that was not anticipated by available theoretical studies. Simulations reveal that, for nematic coacervates, the preferred orientation of the PE chains at the liquid-liquid coacervate-supernatant interface is parallel, and the alignment of semiflexible PEs is homogeneous. The results presented here provide new molecular-level insights into the intracoacervate LCO and will help motivate further experimental and theoretical activities in this area.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Probing electrolyte–silica interactions through simulations of the infrared spectroscopy of nanoscale pores

The structural and dynamical properties of nanoconfined solutions can differ dramatically from those of the corresponding bulk systems. Understanding the changes induced by confinement is central to controlling the behavior of synthetic nanostructured materials and predicting the characteristics of biological and geochemical systems. A key outstanding issue is how the molecular-level behavior of nanoconfined electrolyte solutions is reflected in different experimental, particularly spectroscopic, measurements. This is addressed here through molecular dynamics simulations of the OH stretching infrared (IR) spectroscopy of NaCl, NaBr, and NaI solutions in isotopically dilute HOD/D 2 O confined in hydroxylated amorphous silica slit pores of width 1–6 nm and pH ~2. In addition, the water reorientation dynamics and spectral diffusion, accessible by pump–probe anisotropy and two-dimensional IR measurements, are investigated. In this work, the aim is to elucidate the effect of salt identity, confinement, and salt concentration on the vibrational spectra. It is found that the IR spectra of the electrolyte solutions are only modestly blue-shifted upon confinement in amorphous silica slit pores, with both the size of the shift and linewidth increasing with the halide size, but these effects are suppressed as the salt concentration is increased. This indicates the limitations of linear IR spectroscopy as a probe of confined water. However, the OH reorientational and spectral diffusion dynamics are significantly slowed by confinement even at the lowest concentrations. The retardation of the dynamics eases with increasing salt concentration and pore width, but it exhibits a more complex behavior as a function of halide.

Mesoporous material↗
DOE OSTI2023

Solvent effects on extractant conformational energetics in liquid–liquid extraction: a simulation study of molecular solvents and ionic liquids

Extractant design in liquid–liquid extraction (LLE) is a research frontier of metal ion separations that typically focuses on the direct extractant–metal interactions. However, a more detailed understanding of energetic drivers of separations beyond primary metal coordination is often lacking, including the role of solvent in the extractant phase. In this work, we propose a new mechanism for enhancing metal-complexant energetics with nanostructured solvents. Using molecular dynamics simulations with umbrella sampling, we find that the organic solvent can reshape the energetics of the extractant's intramolecular conformational landscape. Here, we calculate free energy profiles of different conformations of a representative bidentate extractant, n-octyl(phenyl)-N,N-diisobutyl carbamoyl methyl phosphinoxide (CMPO), in four different solvents: dodecane, tributyl phosphate (TBP), and dry and wet ionic liquid (IL) 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ([EMIM][Tf 2 N]). By promoting reorganization of the extractant molecule into its binding conformation, our findings reveal how particular solvents can ameliorate this unfavorable step of the metal separation process. In particular, the charge alternating nanodomains formed in ILs substantially reduce the free energy penalty associated with extractant reorganization. Importantly, using alchemical free energy calculations, we find that this stabilization persists even when we explicitly include the extracted cation. These findings provide insight into the energetic drivers of metal ion separations and potentially suggest a new approach to designing effective separations using a molecular-level understanding of solvent effects.

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

Molecular simulation of glycerol-derived triether podands for lithium ion solvation

Solvate ionic liquids (ILs) are promising candidates for several applications due to their stability, high coulombic efficiency, and low volatility. In this work, we investigate the solvation of lithium-bistriflimide by different glycerol-derived triether solvents, using molecular dynamics simulations. Very strong interactions between Li + and the solvent oxygen sites are found, leading to significant conformational changes in the solvent. Here, by comparing the conformation of the neat solvents with their IL mixtures at different concentrations and temperatures, we find that the presence of Li + induces a distinct crown-like structure in the solvent molecules. The Li + cations and the surrounding solvent form a podand complex, which is stable even at elevated temperatures. These glycerol-derived solvents exhibit distinct interactions with Li + cations which may be exploited in electrolytic applications or lithium recovery processes.

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↗
NASA NTRS1992

Terraced spreading of simple liquids on solid surfaces

We have studied the spreading of liquid drops on a solid surface by molecular-dynamics simulations of coexisting three-phase Lennard-Jones systems of liquid, vapor, and solid. We consider both spherically symmetric atoms and diatomic molecules, and a range of interaction strengths. As the attraction between liquid and solid increases we observe a smooth transition in spreading regimes, from partial to complete to terraced wetting. In the terraced case, where distinct monomolecular layers spread with different velocities, the layers are ordered but not solid, with substantial molecular diffusion both within and between layers. The quantitative behavior resembles recent experimental findings, but the detailed dynamics differ. In particular, the layers exhibit an unusual spreading law, where their radii vary in time as R-squared approximately equal to log10t, which disagrees with experiments on polymeric liquids as well as recent calculations.

Yang, Ju-Xing↗
DOE OSTI2020

Phase-field model of deformation twin-grain boundary interactions in hexagonal systems

Here, a phase-field model for describing dynamic interactions between deformation twins and grain boundaries in a hexagonal close packed (HCP) metal is established. It is applied to simulating the coupled evolution mechanisms of grains and twins in magnesium (Mg) by probing the transmission of $\{10\bar{1}2\}$ deformation twins across grain boundaries. We analyze the effect of the strain relaxation near a grain boundary, which is related to the geometric compatibility arising from the misorientation angle between adjoining grains. Phase-field simulations demonstrate that twin transmission across a grain boundary into a neighboring grain leads to grain boundary migration towards the neighboring grain with a reduced grain boundary width. The preferred nucleation site of new twin variant in the neighboring grain is related to the elastic interaction energy distribution. These predicted transmission behaviors agree well with existing experimental observations and molecular dynamics simulations. By analyzing set of systematic phase-field simulation results, we establish a rotation-angle-related twin variant selection rule for analyzing the transmission of $\{10\bar{1}2\}$ deformation twins across grain boundaries in Mg and its alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Structural basis of agonist specificity of α 1A -adrenergic receptor

α 1 -adrenergic receptors (α 1 -ARs) play critical roles in the cardiovascular and nervous systems where they regulate blood pressure, cognition, and metabolism. However, the lack of specific agonists for all α 1 subtypes has limited our understanding of the physiological roles of different α 1 -AR subtypes, and led to the stagnancy in agonist-based drug development for these receptors. Here we report cryo-EM structures of α 1A -AR in complex with heterotrimeric G-proteins and either the endogenous common agonist epinephrine or the α 1A -AR-specific synthetic agonist A61603. These structures provide molecular insights into the mechanisms underlying the discrimination between α 1A -AR and α 1B -AR by A61603. Guided by the structures and corresponding molecular dynamics simulations, we engineer α 1A -AR mutants that are not responsive to A61603, and α 1B -AR mutants that can be potently activated by A61603. Together, these findings advance our understanding of the agonist specificity for α 1 -ARs at the molecular level, opening the possibility of rational design of subtype-specific agonists.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Atoms to fibers: Identifying novel processing methods in the synthesis of pitch-based carbon fibers

Understanding and optimizing the key mechanisms used in the synthesis of pitch-based carbon fibers (CFs) are challenging, because unlike polyacrylonitrile-based CFs, the feedstock for pitch-based CFs is chemically hetero- geneous, resulting in complex fabrication leading to inconsistency in the final properties. In this work, we use molecular dynamics simulations to explore the processing and chemical phase space through a framework of CF models to identify their effects on elastic performance. The results are in excellent agreement with experiments. We find that density, followed by alignment, and functionality of the molecular constituents dictate the CF mechanical properties more strongly than their size and shape. Last, we propose a previously unexplored fabrication route for high-modulus CFs. Unlike graphitization, this results in increased sp 3 fraction, achieved via generating high-density CFs. In addition, the high sp 3 fraction leads to the fabrication of CFs with isometric compressive and tensile moduli, enabling their potential applications for compressive loading.

01 COAL, LIGNITE, AND PEAT↗