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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Conformational restriction shapes the inhibition of a multidrug efflux adaptor protein

Abstract Membrane efflux pumps play a major role in bacterial multidrug resistance. The tripartite multidrug efflux pump system from Escherichia coli , AcrAB-TolC, is a target for inhibition to lessen resistance development and restore antibiotic efficacy, with homologs in other ESKAPE pathogens. Here, we rationalize a mechanism of inhibition against the periplasmic adaptor protein, AcrA, using a combination of hydrogen/deuterium exchange mass spectrometry, cellular efflux assays, and molecular dynamics simulations. We define the structural dynamics of AcrA and find that an inhibitor can inflict long-range stabilisation across all four of its domains, whereas an interacting efflux substrate has minimal effect. Our results support a model where an inhibitor forms a molecular wedge within a cleft between the lipoyl and αβ barrel domains of AcrA, diminishing its conformational transmission of drug-evoked signals from AcrB to TolC. This work provides molecular insights into multidrug adaptor protein function which could be valuable for developing antimicrobial therapeutics.

59 BASIC BIOLOGICAL SCIENCES↗
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 OSTI2023

Simulations of Glass Transition and Mechanical Behavior of Off-Stoichiometric Crosslinked Polymers

Here, this work explores the influence of blend composition, network architecture, and hydrogen bonding on the material properties of crosslinked epoxy networks, focusing on the glass transition temperature (T g ) and Young’s modulus (Y). We used coarse-grained molecular dynamics simulations to simulate varying compositions of stiff and flexible components in epoxy monomer blends with varying excess of curative. We find that, without hydrogen bonding, networks of any composition show a monotonically increasing T g with decreasing excess curative, consistent with theory. In contrast, we find that when hydrogen bonding is introduced, the binary blend networks show significant enhancement in T g for lightly crosslinked systems. This result contributes to an explanation of the anomalous T g behavior observed experimentally in these systems. We further find that Y is generally enhanced by hydrogen bonds, especially below T g , demonstrating that hydrogen bonding has a significant influence on mechanical properties and can allow access to other desirable dynamic behavior, especially self-healing.

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

Using a Coarse-Grained Modeling Framework to Identify Oligomeric Motifs with Tunable Secondary Structure

Coarse-grained modeling can be used to explore general theories that are independent of specific chemical detail. In this paper, we present cg_openmm, a Python-based simulation framework for modeling coarse-grained hetero-oligomers and screening them for structural and thermodynamic characteristics of cooperative secondary structures. cg_openmm facilitates the building of coarse-grained topology and random starting configurations, setup of GPU-accelerated replica exchange molecular dynamics simulations with the OpenMM software package, and features a suite of postprocessing thermodynamic and structural analysis tools. In particular, native contact analysis, heat capacity calculations, and free energy of folding calculations are used to identify and characterize cooperative folding transitions and stable secondary structures. In this work, we demonstrate the capabilities of cg_openmm on a simple 1–1 Lennard-Jones coarse-grained model, in which each residue contains 1 backbone and 1 side-chain bead. By scanning both nonbonded and bonded force-field parameter spaces at the coarse-grained level, we identify and characterize sets of parameters which result in the formation of stable helices through cooperative folding transitions. Furthermore, we show that the geometries and stabilities of these helices can be tuned by manipulating the force-field parameters.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Exploring the formation of gold/silver nanoalloys with gas-phase synthesis and machine-learning assisted simulations

While nanoalloys are of paramount scientific and practical interest, the main processes leading to their formation are still poorly understood. Key structural features in the alloy systems, including the crystal phase, chemical ordering, and morphology, are challenging to control at the nanoscale, making it difficult to extend their use to industrial applications. Here, in this contribution, we focus on the gold/silver system that has two of the most prevalent noble metals and combine experiments with simulations to uncover the formation mechanisms at the atomic level. Nanoparticles were produced using a state-of-the-art inert-gas aggregation source and analyzed using transmission electron microscopy and energy-dispersive X-ray spectroscopy. Machine-learning-assisted molecular dynamics simulations were employed to model the crystallization process from liquid droplets to nanocrystals. Our study finds a preponderance of nanoparticles with five-fold symmetric morphology, including icosahedra and decahedra which is consistent with previous results on mono-metallic nanoparticles. However, we observed that gold atoms, rather than silver atoms, segregate at the surface of the obtained nanoparticles for all the considered alloy compositions. These segregation tendencies are in contrast to previous studies and have consequences on the crystallization dynamics and the subsequent crystal ordering. We finally showed that the underpinning of this surprising segregation dynamics is due to charge transfer and electrostatic interactions rather than surface energy considerations.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Tight-Binding-Based Method for Simulations of Complex Materials with Accompanying Database

The NRL-Tight-Binding (NRL-TB) method has been applied successfully to many monatomic, binary, and a few ternary systems. The method consists of a fit to Density Functional Theory (DFT) results and has transferability beyond the fitted DFT database including performing tight-binding molecular dynamics simulations. We are working on a new website with a detailed guide of how to use the method in performing the fit to DFT results (Fitting Code) and obtaining non-orthogonal TB parameters expressed in a polynomial form as environment-dependent Slater-Koster parameters (on-site terms), and neighbor distance-dependent Hamiltonian and Overlap terms. These parameters are then used in a Static Code to verify or make predictions for large systems often beyond the computational capabilities of standard DFT codes. Our starting point was a website developed at NRL some years ago which NRL no longer supports. We have now been augmenting that website with new examples of expanding the original capabilities of the method such as incorporating in the codes f-electron states. Our website (esd.cos.gmu.edu/tb) contains two separate branches: One named NRL-TB with the full documentation of the method including source codes and examples. The other branch, named ESD-dB, is a library of electronic structure calculations for all the monatomic elements (a subset of the elements also found in a recent book by the author of this proposal. See reference 22) and an approximately 700 hundred compounds in specific crystal structures. We discuss below the progress made so far in improving and expanding the capabilities of our TB method and the plans for further extensions of our approach.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Simulation insights into wetting properties of hydrogen-brine-clay for hydrogen geo-storage

Hydrogen geo-storage is attracting substantial interdisciplinary interest as a cost-effective and sustainable option for medium- and long-term storage. Hydrogen can be stored underground in diverse formations, including aquifers, salt caverns, and depleted oil and gas reservoirs. The wetting dynamics of the hydrogen-brine-rock system are critical for assessing both structural and residual storage capacities, and ensuring containment safety. Through molecular dynamics simulations, we explore how varying concentrations of cushion gases (CO 2 or CH 4 ) influence the wetting properties of hydrogen-brine-clay systems under geological conditions (15 MPa and 333 K). We employed models of talc and the hydroxylated basal face of kaolinite (kaoOH) as clay substrates. Our findings reveal that the effect of cushion gases on hydrogen-brine-clay wettability is strongly dependent on the clay-brine interactions. Notably, CO 2 and CH 4 reduce the water wettability of talc in hydrogen-brine-talc systems, while exerting no influence on the wettability of hydrogen-brine-kaoOH systems. Detailed analysis of free energy of cavity formation near clay surfaces, clay-brine interfacial tensions, and the Willard-Chandler surface for gas-brine interfaces elucidate the molecular mechanisms underlying wettability changes. Our simulations identify empirical correlations between wetting properties and the average free energy required to perturb a flat interface when clay-brine interactions are less dominant. Here, our thorough thermodynamic analysis of rock-fluid and fluid-fluid interactions, aligning with key experimental observations, underscores the utility of simulated interfacial properties in refining contact angle measurements and predicting experimentally relevant properties. These insights significantly enhance the assessment of gas geo-storage potential. Prospectively, the approaches and findings obtained from this study could form a basis for more advanced multiscale simulations that consider a range of geological and operational variables, potentially guiding the development and improvement of geo-storage systems in general, with a particular focus on hydrogen storage.

25 ENERGY STORAGE↗
DOE OSTI2020

Efficient Treatment of Large Active Spaces through Multi-GPU Parallel Implementation of Direct Configuration Interaction

In this study, we have extended our graphical processing unit (GPU)-accelerated direct configuration interaction program to multiple devices, reducing iteration times for configuration spaces of 165 million determinants to only 3 s using NVIDIA P100 GPUs. Similar improvements in the one- and two-particle reduced density matrix formation allow for fast analytical energy gradients and electronic properties. Our parallel algorithm enables the calculation of arbitrarily large configuration spaces (limited only by available system memory), with iteration times of 13 min for an active space of 18 electrons in 18 orbitals (2.4 billion determinants) using six consumer grade NVIDIA 1080Ti GPUs. These advances enable routine molecular dynamics simulations, geometry optimizations, and absorption spectrum calculations for molecules with large configuration spaces, a task that has heretofore required massive computational effort. In this work, we demonstrate the utility of our program by generating the absorption spectrum for diphenyl acetylene at the floating occupation molecular orbital complete active space configuration interaction level of theory. Lastly, several active spaces were investigated to assess the dependence of spectral features on orbital space dimension.

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

The influence of Holliday junction sequence and dynamics on DNA crystal self-assembly

The programmable synthesis of rationally engineered crystal architectures for the precise arrangement of molecular species is a foundational goal in nanotechnology, and DNA has become one of the most prominent molecules for the construction of these materials. In particular, branched DNA junctions have been used as the central building block for the assembly of 3D lattices. Here, crystallography is used to probe the effect of all 36 immobile Holliday junction sequences on self-assembling DNA crystals. Contrary to the established paradigm in the field, most junctions yield crystals, with some enhancing the resolution or resulting in unique crystal symmetries. Unexpectedly, even the sequence adjacent to the junction has a significant effect on the crystal assemblies. Six of the immobile junction sequences are completely resistant to crystallization and thus deemed “fatal,” and molecular dynamics simulations reveal that these junctions invariably lack two discrete ion binding sites that are pivotal for crystal formation. The structures and dynamics detailed here could be used to inform future designs of both crystals and DNA nanostructures more broadly, and have potential implications for the molecular engineering of applied nanoelectronics, nanophotonics, and catalysis within the crystalline context.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Molecular transport enhancement in pure metallic carbon nanotube porins

Nanofluidic channels impose extreme confinement on water and ions, giving rise to unusual transport phenomena strongly dependent on the interactions at the channel–wall interface. Yet how the electronic properties of the nanofluidic channels influence transport efficiency remains largely unexplored. Here we measure transport through the inner pores of sub-1 nm metallic and semiconducting carbon nanotube porins. We find that water and proton transport are enhanced in metallic nanotubes over semiconducting nanotubes, whereas ion transport is largely insensitive to the nanotube bandgap value. Molecular simulations using polarizable force fields highlight the contributions of the anisotropic polarizability tensor of the carbon nanotubes to the ion–nanotube interactions and the water friction coefficient. We also describe the origin of the proton transport enhancement in metallic nanotubes using deep neural network molecular dynamics simulations. Finally, these results emphasize the complex role of the electronic properties of nanofluidic channels in modulating transport under extreme nanoscale confinement.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

Critical fluctuations in liquid–liquid extraction organic phases controlled by extractant and diluent molecular structure

Extractant aggregation in liquid–liquid extraction organic phases impacts extraction energetics and is related to the deleterious efficiency-limiting liquid–liquid phase transition known as third phase formation. Using small angle X-ray scattering, we find that structural heterogeneities across a wide range of compositions in binary mixtures of malonamide extractants and alkane diluents are well described by Ornstein–Zernike scattering. This suggests that structure in these simplified organic phases originates from the critical point associated with the liquid–liquid phase transition. To confirm this, we measure the temperature dependence of the organic phase structure, finding critical exponents consistent with the 3D Ising model. Molecular dynamics simulations were also consistent with this mechanism for extractant aggregation. Due to the absence of water or any other polar solutes required to form reverse-micellar-like nanostructures, these fluctuations are inherent to the binary extractant/diluent mixture. Further, we also show how the molecular structure of the extractant and diluent modulate these critical concentration fluctuations by shifting the critical temperature: critical fluctuations are suppressed by increasing extractant alkyl tail lengths or decreasing diluent alkyl chain lengths. This is consistent with how extractant and diluent molecular structure are known to impact metal and acid loading capacity in many-component LLE organic phases, suggesting phase behavior of practical systems may be effectively studied in simplified organic phases. Overall, the explicit connection between molecular structure, aggregation and phase behavior demonstrated here will enable the design of more efficient separations processes.

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

Polar liquids at charged interfaces: A dipolar shell theory

The structure of polar liquids and electrolytic solutions, such as water and aqueous electrolytes, at interfaces underlies numerous phenomena in physics, chemistry, biology, and engineering. In this work, we develop a continuum theory that captures the essential features of dielectric screening by polar liquids at charged interfaces, including decaying spatial oscillations in charge and mass, starting from the molecular properties of the solvent. The theory predicts an anisotropic dielectric tensor of interfacial polar liquids previously studied in molecular dynamics simulations. We explore the effect of the interfacial polar liquid properties on the capacitance of the electrode/electrolyte interface and on hydration forces between two plane-parallel polarized surfaces. In the linear response approximation, we obtain simple formulas for the characteristic decay lengths of molecular and ionic profiles at the interface.

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

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

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

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Chain and Segmental Dynamics in Polymer–Nanoparticle Composites with High Nanoparticle Loading

Recent experimental and simulation studies have shown that polymer–nanoparticle (NP) composites (PNCs) with ultra-high NP loading (>50%) exhibit remarkable mechanical properties and dramatic increases in polymer glass-transition temperature, viscosity, and thermal stability compared to the bulk polymer. These deviations in macroscopic properties suggest a slowdown in both segmental and chain-scale polymer dynamics due to confinement. In this work, we examine the polymer conformations and dynamics in these PNCs using molecular dynamics simulations of both unentangled and entangled coarse-grained polymers in random-close-packed NP packings with varying polymer fill fractions. We find that the changes in the polymer dynamics depend on the number of NPs in contact with a polymer segment. Using the number of polymer–NP contacts and different polymer chain conformations as criteria for categorization, we further examine the polymer dynamics at multiple length scales to show the high level of dynamic heterogeneity in PNCs with ultra-high NP loading.

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

Ab-Initio Simulation Studies of Cr Solvation in Fluoride Molten Salts

Understanding molten salt chemistry is essential in ongoing research of the molten salt nuclear reactor (MSR). In this context, detailed understanding of the mechanisms underlying selective oxidation of metal species, such as Cr, is required to guide the design of effective corrosion mitigation strategies in molten salts. An important starting point for such mechanistic understanding is knowledge of the solvation structure and its role in controlling metal speciation. In this work, we use ab initio molecular dynamics simulations to study the short-range (on the scale of the nearest-neighbor bond lengths) and medium-range (over length scales of several neighbor spacings) structure in three different fluoride melts with and without Cr addition; namely, 2KF-NaF, 2LiF-BeF 2 , and 3LiF-AlF3. We find that Cr 0 ,Cr 2+ ,Cr 3+ can each be coordinated by different numbers of F-, with the variance in coordination number decreasing as oxidation state increases, and that these coordination geometries are largely independent of solvent. The manner by which Cr changes the medium-range structure, however, is found to be solvent-dependent. While 2KF-NaF melts show short and medium range order that is highly dynamic, 2LiF-BeF 2 and 3LiF-AlF 3 are characterized by molecular associates that are relatively long-lived that organize into oligomer structures on larger length scales. Rather than being solvated by F- ions alone, we find that Cr can incorporate into and be solvated within this oligomer structure. Fluoroacidity, alone, may therefore prove too simple a metric for assessing the corrosivity of molten fluorides. As our work suggests, the ability of Cr to solvate must be understood in the context of the short- and medium-range structure of the solvent.

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

Electro-reduction of organics on metal cathodes: A multiscale-modeling study of benzaldehyde on Au (111)

Here, we present a multiscale modeling study on the electro-reduction of benzaldehyde (BZY) on Au(111) under realistic electrochemical conditions. To model the electrochemical cell, we adopt a capacitor model in which complex solvents are confined between a Au cathode and a carbon anode. Classical molecular dynamics simulations reveal that electrode charge density and the presence of alcohol show strong effects on the density, adsorption geometry and dynamics of benzaldehyde on the Au electrode. Under charging conditions, the surface concentration of benzaldehyde on the Au electrode decreases, while the content of other species increases. Finally, we proposed a scheme that correlates the electric current running through the Au/solvent interface with the applied bias. This study provides a molecular level understanding of how solvent composition, in this case water/alcohol content, controls the activity of electrocatalytic hydrogenation

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

Resolving Heterogeneous Dynamics of Excess Protons in Aqueous Solution with Rate Theory

Rate theories have found great utility across the chemical sciences by providing a physically transparent way to analyze dynamical processes. Here we demonstrate the benefits of using transition state theory and Marcus theory to study the rate of proton transfer in HCl solutions. By using long ab initio molecular dynamics simulations we show that good agreement is obtained between these two different formulations of rate theory and how they can be used to study the pathways and life-time of proton transfer in aqueous solution. Since both rate theory formulations utilize identical sets of molecular data, this provides a self-consistent theoretical picture of the rates of aqueous phase proton transfer. Specifically, we isolate and quantify the rates of proton transfer, ion-pair dissociation, and solvent exchange in concentrated HCl solutions. Our analysis predicts a concentration dependence to both proton transfer and ion-pairing. Moreover, our estimate of the life-time for the Zundel species is 0.8 ps and 1.3 ps for 2M and 8M HCl, respectively. In this paper, we demonstrate that concentration effects, can indeed be quantified through the combination of state-of-the-art simulation and theory and provides a picture of the important correlations between the cation (hydronium) and the counter-ion in acid solutions.

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

Elucidating the Mechanism of Helium Evaporation from Liquid Water

Here, we investigate the evaporation of trace amounts of helium solvated in liquid water using molecular dynamics simulations and theory. Consistent with experimental observations, we find a super-Maxwellian distribution of kinetic energies of evaporated helium. This excess of kinetic energy over typical thermal expectations is explained by an effective continuum theory of evaporation based on a Fokker–Planck equation, parametrized molecularly by a potential of mean force and position-dependent friction. Using this description, we find that helium evaporation is strongly influenced by the friction near the interface, which is anomalously small near the Gibbs dividing surface due to the ability of the liquid–vapor interface to deform around the gas particle. Our reduced description provides a mechanistic interpretation of trace gas evaporation as the motion of an underdamped particle in a potential subject to a viscous environment that varies rapidly across the air–water interface. From it we predict the temperature dependence of the excess kinetic energy of evaporation, which is yet to be measured.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗