Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Molecular dynamics 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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Water dissociation at the water–rutile TiO 2 (110) interface from ab initio-based deep neural network simulations

The interaction of water with TiO 2 surfaces is of crucial importance in various scientific fields and applications, from photocatalysis for hydrogen production and the photooxidation of organic pollutants to self-cleaning surfaces and bio-medical devices. In particular, the equilibrium fraction of water dissociation at the TiO 2 –water interface has a critical role in the surface chemistry of TiO 2 , but is difficult to determine both experimentally and computationally. Among TiO 2 surfaces, rutile TiO 2 (110) is of special interest as the most abundant surface of TiO 2 ’s stable rutile phase. While surface-science studies have provided detailed information on the interaction of rutile TiO 2 (110) with gas-phase water, much less is known about the TiO 2 (110)–water interface, which is more relevant to many applications. In this work, we characterize the structure of the aqueous TiO 2 (110) interface using nanosecond timescale molecular dynamics simulations with ab initio-based deep neural network potentials that accurately describe water/TiO 2 (110) interactions over a wide range of water coverages. Simulations on TiO 2 (110) slab models of increasing thickness provide insight into the dynamic equilibrium between molecular and dissociated adsorbed water at the interface and allow us to obtain a reliable estimate of the equilibrium fraction of water dissociation. We find a dissociation fraction of 22 ± 6% with an associated average hydroxyl lifetime of 7.6 ± 1.8 ns. These quantities are both much larger than corresponding estimates for the aqueous anatase TiO 2 (101) interface, consistent with the higher water photooxidation activity that is observed for rutile relative to anatase.

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

Exploring the Structural, Dynamic, and Functional Properties of Metal–Organic Frameworks through Molecular Modeling

This review spotlights the role of atomic-level modeling in research on metal-organic frameworks (MOFs), especially the key methodologies of density functional theory (DFT), Monte Carlo (MC) simulations, and molecular dynamics (MD) simulations. The discussion focuses on how periodic and cluster-based DFT calculations can provide novel insights into MOF properties, with a focus on predicting structural transformations, understanding thermodynamic properties and catalysis, and providing information or properties that are fed into classical simulations such as force field parameters or partial charges. Classical simulation methods, highlighting force field selection, databases of MOFs for high-throughput screening, and the synergistic nature of MC and MD simulations, are described. By predicting equilibrium thermodynamic and dynamic properties, these methods offer a wide perspective on MOF behavior and mechanisms. Additionally, the incorporation of machine learning (ML) techniques into quantum and classical simulations is discussed. These methods can enhance accuracy, expedite simulation setup, reduce computational costs, as well as predict key parameters, optimize geometries, and estimate MOF stability. By charting the growth and promise of computational research in the MOF field, the aim is to provide insights and recommendations to facilitate the incorporation of computational modeling more broadly into MOF research.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Wettability of Two-Dimensional Carbon Allotropes from Molecular Simulations

Force-field Monte Carlo and Molecular Dynamics simulations are used to compare wetting behaviors of model carbon sheets mimicking neat graphene, its saturated derivative, graphane, and related planar allotropes penta-graphene, γ-graphyne, and ψ-graphene in contact with aqueous droplets or an aqueous film confined between parallel carbon sheets. Atomistic and area-integrated surface/water potentials are found to be essentially equivalent in capturing moderate differences between the wetting free energies of tested substrates. Despite notable differences in mechanical and electric properties of distinct allotropes, the predicted allotrope/water contact angles span a narrow window of weakly hydrophilic values. Contact angles in the range of 80 ± 10° indicate modest hydration repulsion incapable of competing with van der Waals attraction between carbon particles. Poor dispersibility in neat water is hence a common feature of studied materials.

Biochemistry & Molecular Biology↗
DOE OSTI2020

Simulations of the Biodegradation of Citrate Based Polymers for Artificial Scaffolds Using Accelerated Reactive Molecular Dynamics

Here, we explore the reactivity and mechanical properties of poly(1,6-hexanediol-co-Citric acid) via ReaxFF molecular dynamics simulations. We implement an accelerated scheme within the ReaxFF framework to study the hydrolysis reaction of the polymer which is provided with sufficient amount of energy known as the restrain energy after a suitable pre-transition state configuration is obtained to overcome the activation energy barrier and the desired product is obtained. The validity of the ReaxFF force field is established by comparing the ReaxFF energy barriers of ester and ether hydrolysis with benchmark DFT values in literature. We perform chemical and mechanical degradation of polymer chain bundles at 300 K. We find that ester hydrolyzes faster than ether due to lower activation energy barrier of the reaction. The selectivity of the bond-boost scheme has been demonstrated by lowering the boost parameters of the accelerated simulation which almost stops the ether hydrolysis. Mechanical degradation of pre-hydrolyzed and intermittent hydrolyzed polymer bundles is performed along the longitudinal direction at two different strain rates. We find that the tensile modulus of the polymers are increases with increase in strain rates which shows that polymers show a strain dependent behavior. The tensile modulus of the Polyester-ether is higher than Polyester yet reaches yield-stress faster than Polyester. This makes Polyester more ductile than Polyester-ether.

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

Exploring Li-Ion Transport Properties of Li 3 TiCl 6 : A Machine Learning Molecular Dynamics Study

We performed large-scale molecular dynamics simulations based on a machine-learning force field (MLFF) to investigate the Li-ion transport mechanism in cation-disordered Li 3 TiCl 6 cathode at six different temperatures, ranging from 25°C to 100°C. In this work, deep neural network method and data generated by ab − initio molecular dynamics (AIMD) simulations were deployed to build a high-fidelity MLFF. Radial distribution functions, Li-ion mean square displacements (MSD), diffusion coefficients, ionic conductivity, activation energy, and crystallographic direction-dependent migration barriers were calculated and compared with corresponding AIMD and experimental data to benchmark the accuracy of the MLFF. From MSD analysis, we captured both the self and distinct parts of Li-ion dynamics. The latter reveals that the Li-ions are involved in anti-correlation motion that was rarely reported for solid-state materials. Similarly, the self and distinct parts of Li-ion dynamics were used to determine Haven’s ratio to describe the Li-ion transport mechanism in Li 3 TiCl 6 . Obtained trajectory from molecular dynamics infers that the Li-ion transportation is mainly through interstitial hopping which was confirmed by intra- and inter-layer Li-ion displacement with respect to simulation time. Ionic conductivity (1.06 mS/cm) and activation energy (0.29eV) calculated by our simulation are highly comparable with that of experimental values. Overall, the combination of machine-learning methods and AIMD simulations explains the intricate electrochemical properties of the Li 3 TiCl 6 cathode with remarkably reduced computational time. Thus, our work strongly suggests that the deep neural network-based MLFF could be a promising method for large-scale complex materials.

Selvaraj, Selva Chandrasekaran (ORCID:000000029023↗
DOE OSTI2023

Investigation of the Impact of Lipid Acyl Chain Saturation on Fusion Peptide Interactions with Lipid Bilayers

The interaction of many peptides with lipid bilayer membranes strongly depends on the lipid composition. Here, a study of the impact of unsaturated lipid acyl chains on the interaction of a derivative of the HIV-1 fusion peptide with lipid bilayer vesicles is presented. Lipid bilayer vesicles composed of mixtures of lipids with two saturated acyl chains and lipids and one saturated and one unsaturated acyl chain, but identical head groups, were studied. The dependence of the peptide conformation on the unsaturated lipid content was probed by circular dichroism spectroscopy, while the impact of the peptide on the bilayer structure was determined by small-angle neutron scattering. The impact of the peptide on the lipid bilayer vesicle dynamics was investigated using neutron spin echo spectroscopy. Molecular dynamics simulations were used to characterize the behavior of the systems studied to determine if there were clear differences in their physical properties. The results reveal that the peptide–bilayer interaction is not a simple function of the unsaturated lipid acyl chain content of the bilayer. Instead, the peptide behavior is more consistent with that seen for the bilayer containing only unsaturated lipids, which is supported by lipid-specific interactions revealed by the simulations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Scale and rate in CdS pressure-induced phase transition

Here, we describe recent efforts to improve our predictive modeling of rate-dependent behavior at, or near, a phase transition using molecular dynamics simulations. Cadmium sulfide (CdS) is a well-studied material that undergoes a solid-solid phase transition from wurtzite to rock salt structures between 3 and 9 GPa. Atomistic simulations are used to investigate the dominant transition mechanisms as a function of orientation, size and rate. We found that the final rock salt orientations were determined relative to the initial wurtzite orientation, and that these orientations were different for the two orientations and two pressure regimes studied. The CdS solid-solid phase transition is studied, for both a bulk single crystal and for polymer-encapsulated spherical nanoparticles of various sizes.

42 ENGINEERING↗
DOE OSTI2020

Imogolite in water: Simulating the effects of nanotube curvature on structure and dynamics

Imogolite is a fascinating inorganic nanotube that is found in nature or synthesized in a laboratory. The synthesis process is carried out in liquid media, and leads to the formation of almost monodisperse diameter nanotubes. Here we investigate, employing classical molecular dynamics simulations, the interaction of water and imogolite for nanotubes of several radii. We established that water penetrates the pores of N = 9 and larger nanotubes, and adopts a coaxial arrangement in it. Also, while water molecules can diffuse along the center of the nanotube, the molecules next to the inner imogolite walls have very low mobility. At the outer nanotube wall, an increase of water density is observed, this effect extends up to 1 nm, beyond which water properties are bulk-like. Both phenomena are affected by the imogolite curvature.

58 GEOSCIENCES↗
DOE OSTI2023

Building a DFT+U machine learning interatomic potential for uranium dioxide

Despite uranium dioxide (UO 2 ) being a widely used nuclear fuel, fuel performance models rely extensively on empirical correlations of material behavior, leveraging the historical operating experience of UO 2 . Mechanistic models that consider an atomistic understanding of the processes governing fuel performance (such as fission gas release and creep) will enable a better description of fuel behavior under non-prototypical conditions such as in new reactor concepts or for modified UO 2 fuel compositions. To this end, molecular dynamics simulation is a powerful tool for rapidly predicting physical properties of proposed fuel candidates. However, the reliability of these simulations depends largely on the accuracy of the atomic forces. Traditionally, these forces are computed using either a classical force field (FF) or density functional theory (DFT). While DFT is relatively accurate, the computational cost is burdensome, especially for f-electron elements, such as actinides. By contrast, classical FFs are computationally efficient but are less accurate. For these reasons, we report a new accurate machine learning interatomic potential (MLIP) for UO 2 that provides high-fidelity reproduction of DFT forces at a similar low cost to classical FFs. We employ an active learning approach that autonomously augments the DFT training data set to iteratively refine the MLIP. To further improve the quality of our predictions, we utilize transfer learning to retrain our MLIP to higher-accuracy DFT+U data. We validate our MLIPs by comparing predicted physical properties (e.g., thermal expansion and elastic properties) with those from existing classical FFs and DFT/DFT+U calculations, as well as with experimental data when available.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

PySpawn: Software for Nonadiabatic Quantum Molecular Dynamics

The ab initio multiple spawning (AIMS) method enables nonadiabatic quantum molecular dynamics simulations in an arbitrary number of dimensions, with potential energy surfaces provided by electronic structure calculations performed on-the-fly. However, the intricacy of the AIMS algorithm complicates software development, deployment on modern shared computer resources, and post-simulation data analysis. PySpawn is a nonadiabatic molecular dynamics software package that addresses these issues. Here, the program is designed to be easily interfaced with electronic structure software, and an interface to the TeraChem software package is described here. PySpawn introduces a task-based reorganization of the AIMS algorithm, allowing fine-grained restart capability and setting the stage for efficient parallelization in a future release. PySpawn includes a user-friendly and interactive Python analysis module that will enable novice users to painlessly adopt AIMS. As a demonstration of PySpawn’s simulation capability and analysis module, we report complete active space self-consistent field–based AIMS simulations of the 1,2- dithienyl-1,2-dicyanoethene molecule, a promising molecular photoswitch.

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

State predictive information bottleneck

We report the ability to make sense of the massive amounts of high-dimensional data generated from molecular dynamics simulations is heavily dependent on the knowledge of a low-dimensional manifold (parameterized by a reaction coordinate or RC) that typically distinguishes between relevant metastable states, and which captures the relevant slow dynamics of interest. Methods based on machine learning and artificial intelligence have been proposed over the years to deal with learning such low-dimensional manifolds, but they are often criticized for a disconnect from more traditional and physically interpretable approaches. To deal with such concerns, in this work we propose a deep learning based state predictive information bottleneck approach to learn the RC from high-dimensional molecular simulation trajectories. We demonstrate analytically and numerically how the RC learnt in this approach is connected to the committor in chemical physics and can be used to accurately identify transition states. A crucial hyperparameter in this approach is the time delay or how far into the future the algorithm should make predictions about. Through careful comparisons for benchmark systems, we demonstrate that this hyperparameter choice gives useful control over how coarse-grained we want the metastable state classification of the system to be. We thus believe that this work represents a step forward in systematic application of deep learning based ideas to molecular simulations.

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

Ions Tune Interfacial Water Structure and Modulate Hydrophobic Interactions at Silica Surfaces

The structure and ultrafast dynamics of the electric double layer (EDL) are central to chemical reactivity and physical properties at solid/aqueous interfaces. While the Gouy–Chapman–Stern model is widely used to describe EDLs, it is solely based on the macroscopic electrostatic attraction of electrolytes for the charged surfaces. Structure and dynamics in the Stern layer are, however, more complex because of competing effects due to the localized surface charge distribution, surface–solvent–ion correlations, and the interfacial hydrogen bonding environment. Here, we report combined time-resolved vibrational sum frequency generation (TR-vSFG) spectroscopy with ab initio DFT-based molecular dynamics simulations (AIMD/DFT-MD) to get direct access to the molecular-level understanding of how ions change the structure and dynamics of the EDL. We show that innersphere adsorbed ions tune the hydrophobicity of the silica–aqueous interface by shifting the structural makeup in the Stern layer from dominant water–surface interactions to water–water interactions. This drives an initially inhomogeneous interfacial water coordination landscape observed at the neat interface toward a homogeneous, highly interconnected in-plane 2D hydrogen bonding (2D-HB) network at the ionic interface, reminiscent of the canonical, hydrophobic air–water interface. This ion-induced transformation results in a characteristic decrease of the vibrational lifetime (T 1 ) of excited interfacial O–H stretching modes from T 1 ~ 600 fs to T 1 ~ 250 fs. Hence, we propose that the T 1 determined by TR-vSFG in combination with DFT-MD simulations can be widely used for a quantitative spectroscopic probe of the ion kosmotropic/chaotropic effect at aqueous interfaces as well as of the ion-induced surface hydrophobicity.

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

Molecular-Level Insights into the NMR Relaxivity of Gadobutrol Using Quantum and Classical Molecular Simulations

MRI is an indispensable diagnostic tool in modern medicine; however, understanding the molecular-level processes governing NMR relaxation of water in the presence of MRI contrast agents remains a challenge, hindering the molecular-guided development of more effective contrast agents. By using quantum-based polarizable force fields, the first-of-its-kind molecular dynamics (MD) simulations of Gadobutrol are reported where the 1 H NMR longitudinal relaxivity r 1 of the aqueous phase is determined without any adjustable parameters. The MD simulations of r 1 dispersion (i.e., frequency dependence) show good agreement with measurements at frequencies of interest in clinical MRI. Importantly, the simulations reveal key insights into the molecular level processes leading to r 1 dispersion by decomposing the NMR dipole–dipole autocorrelation function G(t) into a discrete set of molecular modes, analogous to the eigenmodes of a quantum harmonic oscillator. The molecular modes reveal important aspects of the underlying mechanisms governing r 1 , such as its multiexponential nature and the importance of the second eigenmodal decay. By simply analyzing the MD trajectories on a parameter-free approach, the Gadobutrol simulations show that the outer-shell water contributes ∼50% of the total relaxivity r 1 compared to the inner-shell water, in contrast to simulations of (nonchelated) gadolinium-aqua where the outer shell contributes only ∼15% of r 1 . The deviation between simulations and measurements of r 1 below clinical MRI frequencies is used to determine the low-frequency electron-spin relaxation time for Gadobutrol, in good agreement with independent studies.

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

Transient nature of fast relaxation in metallic glass

Metallic glasses exhibit fast mechanical relaxations at temperatures well below the glass transition, one of which shows little variation with temperature known as nearly constant loss (NCL). Despite the important implications of this phenomenon to deformation, the origin of the relaxation is unclear. Through molecular dynamics simulations of a model metallic glass, Cu 64.5 Zr 35.5 , we implement molecular dynamics dynamical mechanical spectroscopy (MD-DMS) with system stress decomposed into atomic-level stresses to identify the group of atoms responsible for NCL. This work demonstrates that NCL relaxation is due to transient groups of atoms that revert to the typical atomic-level viscoelastic behavior over picosecond timescales. They are homogenously distributed throughout the glass and have no outstanding features, rather than having defect-like local structure as previously reported.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Stable nanovesicles formed by intrinsically planar bilayers

Quatsome nanovesicles, formed through the self-assembly of cholesterol (CHOL) and cetyltrimethylammonium bromide (CTAB) in water, have shown long-term stability in terms of size and morphology, while at the same time exhibiting high CHOL-CTAB intermolecular binding energies. We hypothesize that CHOL/CTAB quatsomes are indeed thermodynamically stable nanovesicles, and investigate the mechanism underlying their formation. A systematic study was performed to determine whether CHOL/CTAB quatsomes satisfy the experimental requisites of thermodynamically stable vesicles. Coarse-grain molecular dynamics simulations were used to investigate the molecular organization in the vesicle membrane, and the characteristics of the simulated vesicle were corroborated with experimental data obtained by cryo–electron microscopy, small- and wide-angle X-ray scattering, and multi-angle static light scattering. CHOL/CTAB quatsomes fulfill the requisites of thermodynamically stable nanovesicles, but they do not exhibit the classical membrane curvature induced by a composition asymmetry between the bilayer leaflets, like catanionic nanovesicles. Instead, CHOL/CTAB quatsomes are formed through the association of intrinsically planar bilayers in a faceted vesicle with defects, indicating that distortions in the organization and orientation of molecules can play a major role in the formation of thermodynamically stable nanovesicles.

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

Correlations between Molecular Structure, Solvation Topology, and Transport Properties of Aqueous Organic Flow Battery Electrolyte Solutions

Aqueous organic redox flow batteries (AORFBs) are considered promising technologies for storing energy generated from renewable resources. However, designing organic electrolyte molecules is limited by gaps between fundamental understanding of coupling between solvation structure and dynamics, and macroscopic transport properties like viscosity. Herein, we used molecular dynamics simulations to understand correlations between ionic molecular structures, ion clustering, and transport properties in 2,3-dihydrophenazine (2,3-DHP), a promising AORFB anolyte. We show that experimentally measured viscosity can be reproduced from simulations at relevant concentrations and that the asymmetric structure of 2,3-DHP leads to unique inhomogeneity in the solvation topology. However, order parameters and metrics need to be developed for better correlations over spatiotemporal scales, with careful consideration of the inhomogeneity of organic anolyte molecules. In conclusion, we show that the increased size and asymmetry of the anolyte leads to breakdown of assumptions within methods for determining ion transport mechanisms previously developed for Li-ion batteries.

25 ENERGY STORAGE↗
DOE OSTI2026

Efficient simulation of optical spectra via machine learning and physical decomposition of environmental effects

Simulations of optical spectra can provide key insights to aid experimental interpretation of electronic excitation phenomena. For chromophores in the condensed phase, these spectra, which incorporate the coupling between electronic excitation and molecular and solvent nuclear motions, can be simulated using excitation energies obtained from molecular dynamics simulations of the chromophore and solvent. Here, we present a hybrid scheme that exploits machine learning and physically informed spectral densities to show that as few as 25 ground and excited state energetic gradient calculations can be used to construct models that accurately predict environment-influenced vibronic coupling in optical spectra. We demonstrate our approach for the green fluorescent protein chromophore in water and the cresyl violet chromophore in methanol. We show that our hybrid approach, employing a machine learning model for the high-frequency spectral density and an ab initio parameterized Debye spectral density for the low-frequency, results in a systematic improvement of the optical absorption lineshape, leading to a simple machine learning scheme that can be used for the simulation of spectral densities and optical spectra.

Snider, Andrew [Univ. of California, Merced, CA (U↗
DOE OSTI2023

A High-Speed Rotational Diamond Anvil Cell for In Situ Analysis of Hierarchical Microstructural Evolution of Metallic Alloys during Extreme Shear Deformation

High speed shear deformation is ubiquitous in engineering applications, ranging from material processing methods such as friction stir processing/extrusion and in tribological contacts. However, analyzing the microstructural evolution of materials while they are undergoing high speed shear deformation have been a long-standing challenge. This led to predominant reliance on ex situ microscopy before and after shear deformation. But ex situ microscopy lacks the ability to analyze dynamic and transient hierarchical microstructural evolution mechanisms that could occur during shear deformation of materials. Therefore, to better understand the dynamic mechanisms of mass and energy transfer in materials under shear deformation, we developed a first of its kind high-speed rotational diamond anvil cell (HS-RDAC) for synchrotron-based in situ high-energy x-ray diffraction (XRD). We studied the time resolved lattice strain evolution, XRD peak broadening and changes in spatial variation of shear deformation induced alloying in pure metal and metal alloy sheets and powder mixture using the HS-RDAC. These in situ results were combined with detailed ex situ microstructural characterization before and after the shear deformation using transmission electron microscopy and atom probe tomography, which revealed the different stages of evolution of a shear deformation induced hierarchical nanostructure. Multiscale computational simulations including computational fluid dynamics, crystal plasticity, molecular dynamic simulation and density functional theory uncovered the mechanisms behind morphological changes, evolution of defect structures and changes in driving force for shear deformation induced intermixing. In conclusion, this in situ HS-RDAC capability, in combination with ex situ microstructural characterization and computational simulations, can provide new insights into the hierarchical microstructural evolution pathway during shear deformation.

36 MATERIALS SCIENCE↗