Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecular dynamic (MD) 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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Orientation and Conformation of Proteins at the Air–Water Interface Determined from Integrative Molecular Dynamics Simulations and Sum Frequency Generation Spectroscopy

Understanding the assembly of proteins at the air-water interface (AWI) informs the formation of protein films, emulsion properties, and protein aggregation. Determination of protein conformation and orientation at an interface is difficult to resolve with a single experimental or simulation technique alone. To date, the interfacial structure of even one of the most widely studied proteins, lysozyme, at the AWI remains unresolved. Here in this study, molecular dynamics (MD) simulations are used to determine if the protein adopts a side-on, head-on, or axial orientation at the AWI with two different forcefields, GROMOS-53a6 + SPC/E and a99SB-disp + TIP4P-D. Vibrational sum frequency generation (SFG) spectroscopy experiments and spectral SFG calculations validate consistency between the structure determined from MD and experiments. Overall, we show with strong agreement that lysozyme adopts an axial conformation at pH 7. Further, we provide molecular-level insight as to how pH influences the binding domains of lysozyme resulting in side-on adsorption near the isoelectric point of the lysozyme.

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

All-atom molecular dynamics simulations of weak polyionic brushes: influence of charge density on the properties of polyelectrolyte chains, brush-supported counterions, and water molecules

All atom molecular dynamics (MD) simulations of planar Na+-counterion-neutralized polyacrylic acid (PAA) brushes are performed for varying degrees of ionization (and thereby varying charge density) and varying grafting density. Variation in the PE charge density (or degree of ionization) and grafting density leads to massive changes of the properties of the PE molecules (quantified by the changes in the height and the mobility of the PE brushes) as well as the local arrangement and distribution of the brush-supported counterions and water molecules within the brushes. The effect on the counterions is manifested by the corresponding variation of the counterion mobility, counterion concentration, extent of counterion binding to the charged site of the PE brushes, water-in-salt-like structure formation, and counterion-water–oxygen radial distribution function within the PE brushes. On the other hand, the effect on water molecules is manifested by the corresponding variation of water–oxygen–water–oxygen RDF, local water density, water–water and water–PE functional group hydrogen bond networks, static dielectric constant of water molecules, orientational tetrahedral order parameter, and water mobility. Enforcing such varying degree of ionization of weak polyelectrolytes is possible by changing the pH of the surrounding medium. Thus, our results provide insights into the changes in microstructure (at the atomistic level) of weak polyionic brushes at varying pH. Finally, we anticipate that this knowledge will prove to be vital for the efficient design of several nano-scale systems employing PE brushes such as nanomechanical gates, current rectifiers, etc.

42 ENGINEERING↗
DOE OSTI2021

Investigations of water/oxide interfaces by molecular dynamics simulations

Water/oxide interfaces are ubiquitous on earth and show significant influence on many chemical processes. For example, understanding water and solute adsorption as well as catalytic water splitting can help build better fuel cells and solar cells to overcome our looming energy crisis; the interaction between biomolecules and water/oxide interfaces is one hypothesis to explain the origin of life. However, knowledge in this area is still limited due to the difficulty of studying water/solid interfaces. As a result, research using increasingly sophisticated experimental techniques and computational simulations has been carried out in recent years. Although it is difficult for experimental techniques to provide detailed microscopic structural information, molecular dynamics (MD) simulations have satisfactory performance. In this report, we discuss classical and ab initio MD simulations of water/oxide interfaces. Generally, we are interested in the following questions: How do solid surfaces perturb interfacial water structure? How do interfacial water molecules and adsorbed solutes affect solid surfaces and how do interfacial environments affect solvent and solute behavior? Finally, we discuss progress in the application of neural network potential based MD simulations, which offer a promising future because this approach has already enabled ab initio level accuracy for very large systems and long trajectories.

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

Machine Learning Predictions of Simulated Self-Diffusion Coefficients for Bulk and Confined Pure Liquids

Diffusion properties of bulk fluids have been predicted using empirical expressions and machine learning (ML) models, suggesting that predictions of diffusion also should be possible for fluids in confined environments. The ability to quickly and accurately predict diffusion in porous materials would enable new discoveries and spur development in relevant technologies such as separations, catalysis, batteries, and subsurface applications. Here in this work, we apply artificial neural network (ANN) models to predict the simulated self-diffusion coefficients of real liquids in both bulk and pore environments. The training data sets were generated from molecular dynamics (MD) simulations of Lennard-Jones particles representing a diverse set of 14 molecules ranging from ammonia to dodecane over a range of liquid pressures and temperatures. Planar, cylindrical, and hexagonal pore models consisted of walls composed of carbon atoms. Our simple model for these liquids was primarily used to generate ANN training data, but the simulated self-diffusion coefficients of bulk liquids show excellent agreement with experimental diffusion coefficients. ANN models based on simple descriptors accurately reproduced the MD diffusion data for both bulk and confined liquids, including the trend of increased mobility in large pores relative to the corresponding bulk liquid.

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

De Novo Design of Proteins That Bind Naphthalenediimides, Powerful Photooxidants with Tunable Photophysical Properties

De novo protein design provides a framework to test our understanding of protein function and build proteins with cofactors and functions not found in nature. Here, we report the design of proteins designed to bind powerful photooxidants and the evaluation of the use of these proteins to generate diffusible small-molecule reactive species. Because excited-state dynamics are influenced by the dynamics and hydration of a photooxidant’s environment, it was important to not only design a binding site but also to evaluate its dynamic properties. Thus, we used computational design in conjunction with molecular dynamics (MD) simulations to design a protein, designated NBP (NDI Binding Protein), that held a naphthalenediimide (NDI), a powerful photooxidant, in a programmable molecular environment. Solution NMR confirmed the structure of the complex. We evaluated two NDI cofactors in this de novo protein using ultrafast pump–probe spectroscopy to evaluate light-triggered intra- and intermolecular electron transfer function. Moreover, we demonstrated the utility of this platform to activate multiple molecular probes for protein labeling.

carbonyls↗
DOE OSTI2020

Application of First-Principles-Based Artificial Neural Network Potentials to Multiscale-Shock Dynamics Simulations on Solid Materials

The use of artificial neural network (ANN) potentials trained with first-principles calculations has emerged as a promising approach for molecular dynamics (MD) simulations encompassing large space and time scales while retaining first-principles accuracy. To date, however, the application of ANN-MD has been limited to near-equilibrium processes. Here we combine first-principles-trained ANN-MD with multiscale shock theory (MSST) to successfully describe far-from-equilibrium shock phenomena. Overall, our ANN-MSST-MD approach describes shock-wave propagation in solids with first-principles accuracy but a 5000 times shorter computing time. Accordingly, ANN-MD-MSST was able to resolve fine, long-time elastic deformation at low shock speed, which was impossible with first-principles MD because of the high computational cost. This work thus lays a foundation of ANN-MD simulation to study a wide range of far-from-equilibrium processes.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Applying molecular dynamics simulations to the study of micro cold spray of yttria

Micro cold spray (MCS) is a novel manufacturing process with proven success depositing dense, thick films of technical ceramics without the need for elevated processing temperatures or chemical precursors. The process transforms a loose sub-micrometer powder into an adherent film with mechanical properties approaching that of the bulk material through a phenomenon known as room temperature impact consolidation. While this phenomenon has been demonstrated for a wide variety of ceramic materials and film substrates, yttria stands out as a material that is both well-characterized in experiments and of significant interest for industrial applications due to its excellent refractory properties. This dissertation seeks to investigate the fundamental mechanisms that enable room temperature impact consolidation by using molecular dynamics (MD) simulations of yttria particle impacts and characterizing the deformation behavior of the system.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Computational Design of Miniproteins as SARS-CoV-2 Therapeutic Inhibitors

A rational therapeutic strategy is urgently needed for combating SARS-CoV-2 infection. Viral infection initiates when the SARS-CoV-2 receptor-binding domain (RBD) binds to the ACE2 receptor, and thus, inhibiting RBD is a promising therapeutic for blocking viral entry. In this study, the structure of lead antiviral candidate binder (LCB1), which has three alpha-helices (H1, H2, and H3), is used as a template to design and simulate several miniprotein RBD inhibitors. LCB1 undergoes two modifications: structural modification by truncation of the H3 to reduce its size, followed by single and double amino acid substitutions to enhance its binding with RBD. We use molecular dynamics (MD) simulations supported by ab initio density functional theory (DFT) calculations. Complete binding profiles of all miniproteins with RBD have been determined. The MD investigations reveal that the H3 truncation results in a small inhibitor with a -1.5 kcal/mol tighter binding to RBD than original LCB1, while the best miniprotein with higher binding affinity involves D17R or E11V + D17R mutation. DFT calculations provide atomic-scale details on the role of hydrogen bonding and partial charge distribution in stabilizing the minibinder:RBD complex. This study provides insights into general principles for designing potential therapeutics for SARS-CoV-2.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Xe Recovery from Nuclear Power Plants Off-Gas Streams: Molecular Simulations of Gas Permeation through DD3R Zeolite Membrane

Recent experimental work has shown zeolite membrane-based separation as a promising potential technology for Kr/Xe gas mixtures due to its much lower energy requirements in comparison to cryogenic distillation, the conventional separation method for such mixtures. Such a separation is also economically rewarding because Xe is in high demand, as a valuable product for many applications/processes. In this work, we have used Molecular Dynamics (MD) simulations to study the effects of different conditions, i.e., temperature, pressure, and gas feed composition, on Kr/Xe separation performance via DD3R zeolite membranes. We provide a comprehensive study of the permeation of the different gas species, density profiles, and diffusion coefficients. Molecular simulations show that if the feed is changed from pure Kr/Xe to an equimolar mixture, the Kr/Xe separation factor increases, which agrees with experiments. In addition, when Ar is introduced as a sweep gas, the adsorption of both Kr and Xe increases, while the permeation of pure Kr increases. A similar behavior is observed with equimolar mixtures of Kr/Xe with Ar as the sweep gas. High-separation Kr/Xe selectivity is observed at 50 atm and 425 K but with low total permeation rates. Changing pressure and temperature are found to have profound effects on optimizing the separation selectivity and the permeation throughput.

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

Thermal and concentration effects on 1 H NMR relaxation of Gd 3+ -aqua using MD simulations and measurements

We report gadolinium-based contrast agents are key in clinical MRI for enhancing the longitudinal NMR relativity (r 1 ) of hydrogen nuclei ( 1 H) in water and improving the contrast among different tissues. The importance of MRI in clinical practice cannot be gainsaid, yet the interpretation of MRI relies on models with severe assumptions, reflecting a poor understanding of the molecular-scale relaxation processes. In a step towards building a clearer understanding of the relaxation processes, here we investigate thermal and concentration effects on r 1 of the Gd 3+ -aqua complex using both semi-classical molecular dynamics (MD) simulations and measurements. We follow the MD simulation approach recently introduced by [Singer et al., Phys. Chem. Chem. Phys., 2021, 23, 20974], in which no NMR relaxation model or free-parameter is assumed to predict r 1 , thereby bringing new insights into the physics of r 1 on a molecular scale. We expand the autocorrelation function G(t) in terms of molecular modes and determine the thermal activation energies of the two largest modes, both of which are consistent with the range of literature values for rotational diffusion. We also determine the activation energies for translational diffusion and low-field electron-spin relaxation, both of which are consistent with the literature. Furthermore, we validate the MD simulations at human body temperature and concentrations of the paramagnetic ion used in clinical MRI, and we quantify the uncertainties in both simulations and measurements.

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

Modelling dislocation-graphene interactions in a BCC Fe matrix by molecular dynamics simulations and gradient plasticity theory

We report graphene nanosheets (GNS) can enhance the strength and ductility of metal-based composites as they can obstruct the propagation of dislocations. The present article employs Molecular Dynamics (MD) simulations to investigate dislocation-GNS interaction mechanisms and possible influencing factors, including the number of GNS layers, the thickness of the metallic amorphous layer and the C - C bond strength. The results indicated that the shear strength of the metal/GNS interface and the bending stiffness of GNS determined the ability of GNS to block dislocation transmission. A physically based phenomenological parameter that can capture such dislocation-GNS interactions is the mechanical interface energy that has been put forth within gradient plasticity. By fitting the theoretical expressions to the simulation data, it was possible to obtain estimates for the mechanical interface energy for the GNS. It was found that increasing the GNS layers and adding an amorphous layer resulted in a strengthening in the stress–strain response and increased the value of this interfacial parameter. This indicates that the mechanical interfacial energy can be a unified measure for capturing and tuning the strength of various interfaces such as grain boundaries, GNS, amorphous-crystalline interface and bimetal interfaces.

36 MATERIALS SCIENCE↗
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 OSTI2023

Sensitivity of Dislocation-GB interactions to simulation setups in atomistic models

Dislocation-grain boundary (GB) interactions play an integral role in strengthening of crystalline materials, and are dependent on external loading conditions and atomic arrangements at the grain boundary. While molecular dynamics (MD) simulations can provide critical insights into relationships between these parameters and the exact interaction, the outcomes are sensitive to the computational setup. Here, in this work, we explore the effect of computational setup (system size, GB orientation and boundary conditions) on stress-controlled dislocation-grain boundary (DGB) reactions using MD simulations for a well-studied copper symmetric $\langle$1 1 2$\rangle$ tilt boundary. The study demonstrates that the DGB reactions (pinning, absorption, transmission, etc.) are sensitive to the system size and the orientation of the GB to the applied stress. Additionally, the current work shows that there is a critical system size for the dislocation-grain boundary reaction stresses to converge, one which was found to be much larger than the setups in comparable MD studies. This was attributed to the specific stress states and the influence of surfaces, fixed or free, as these computational setups are modified. Finally, a stress-controlled setup with minimal model artifacts is examined, and the effect of coupled Schmid and non-Schmid stress components on the dislocation-GB reactions are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Si–Cl 2 –Ar + Atomic Layer Etching Window: A Fundamental Study Using Molecular Dynamics Simulations and a Reduced Order Model

Silicon (Si) atomic layer etching (ALE) by alternating exposure to chlorine gas (Cl 2 ) and argon ions (Ar + ) is studied by using molecular dynamics (MD) simulations and a reduced order model (ROM). Here, the purpose of this study is to elucidate the properties of the ALE window, a range of ion energies where the amount of Si etched over a series of cycles is nonzero and nearly independent of ion energy. Experimental studies of the Si–Cl 2 –Ar + ALE system report contradictory results related to the ALE window’s ion energy range. Both MD simulations and the ROM show that there is an ALE window present from approximately 15 to 20 eV for normal incidence argon ions. The Si–Cl 2 –Ar + system, therefore, exhibits a narrow ALE window. The amount of Si etched per cycle is less than one atomic layer because of the higher etch yield of Cl atoms relative to atomic Si and silicon chlorides. A modified version of the ROM with an artificially increased Si physical sputtering threshold energy expands the ALE window, illustrating the importance of the difference in chemical and physical sputtering threshold energies in the ALE window energy range. The ROM is also used to examine the dependence of the EPC on the Ar + ion fluence.

energy↗
DOE OSTI2023

The impact of physiologically relevant temperatures on physical properties of thylakoid membranes: a molecular dynamics study

Thylakoid membranes are energy-converting membranes with a unique lipid composition. Though the membranes are primarily composed of proteins, their photosynthetic function is strongly influenced by the lipid constituents. Here we characterize, with molecular dynamics (MD) simulations, lipid bilayers with compositions representative of plant thylakoid membranes. We determine, in a wide range of temperatures, the physical parameters of the model membranes which are relevant for the photosynthetic function. We found a marked impact of temperature on membrane permeability due to a combination of increased compressibility and curvature of the membrane at elevated temperatures. With increasing temperatures, we observed increasingly smeared transmembrane density profiles of the membrane forming lipid headgroups predicting increased membrane flexibility. The diffusion coefficient of the lipids increased with temperature without apparent specificity for lipid species. Instead of a comprehensive experimental dataset in the relevant temperature range, we quantitatively compared and validated our MD results with MD simulations on a dipalmitoylphosphatidylcholine model system.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

MISPR : an open-source package for high-throughput multiscale molecular simulations

Computational tools provide a unique opportunity to study and design optimal materials by enhancing our ability to comprehend the connections between their atomistic structure and functional properties. However, designing materials with tailored functionalities is complicated due to the necessity to integrate various computational-chemistry software (not necessarily compatible with one another), the heterogeneous nature of the generated data, and the need to explore vast chemical and parameter spaces. The latter is especially important to avoid bias in scattered data points-based models and derive statistical trends only accessible by systematic datasets. Here, we introduce a robust high-throughput multi-scale computational infrastructure coined MISPR (Materials Informatics for Structure–Property Relationships) that seamlessly integrates classical molecular dynamics (MD) simulations with density functional theory (DFT). By enabling high-performance data analytics and coupling between different methods and scales, MISPR addresses critical challenges arising from the needs of automated workflow management and data provenance recording. The major features of MISPR include automated DFT and MD simulations, error handling, derivation of molecular and ensemble properties, and creation of output databases that organize results from individual calculations to enable reproducibility and transparency. In this work, we describe fully automated DFT workflows implemented in MISPR to compute various properties such as nuclear magnetic resonance chemical shift, binding energy, bond dissociation energy, and redox potential with support for multiple methods such as electron transfer and proton-coupled electron transfer reactions. The infrastructure also enables the characterization of large-scale ensemble properties by providing MD workflows that calculate a wide range of structural and dynamical properties in liquid solutions. MISPR employs the methodologies of materials informatics to facilitate understanding and prediction of phenomenological structure–property relationships, which are crucial to designing novel optimal materials for numerous scientific applications and engineering technologies.

36 MATERIALS SCIENCE↗
NASA NTRS2016

Interfacial Structure and Capacitance of Li-doped Ionic Liquid Electrolytes from Molecular Simulations

Ionic liquids have been proposed as candidate electrolytes for high-energy density, rechargeable batteries, supercapacitors, and hybrid energy storage devices. Though Li-salt is often present in these systems, its influence on interfacial properties is largely uncharacterized. We, thereby, present an extensive computational analysis, supported by experimental comparisons, of the properties of a representative set of these electrolytesat an ideal carbon interface as a function of Li-salt doping and voltage. We have performed polarizable molecular (MD) dynamics simulations, using the APPLEP force field, to evaluate electric double layer (EDL) capacitance and distribution of Li+ in the EDL. Differential capacitance exhibits the characteristic camel profile and is insensitive to Li-doping. Li+ localizes in the second molecular layer of the EDL, which is a result of confinement from free energy barriers associated with ion layering. Joint MDelectronic structure computations show the electrochemical window of the electrolytes to be a weak function of Li-doping. Estimates of supercapacitor specific energy are made using the computed window and capacitance. The magnitude and trends in specific energy are in good agreement with experiment.

batteries↗