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 1,387 records · Page 77
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Discovery of Small-Molecule Inhibitors of SARS-CoV-2 Proteins Using a Computational and Experimental Pipeline

A rapid response is necessary to contain emergent biological outbreaks before they can become pandemics. The novel coronavirus (SARS-CoV-2) that causes COVID-19 was first reported in December of 2019 in Wuhan, China and reached most corners of the globe in less than two months. In just over a year since the initial infections, COVID-19 infected almost 100 million people worldwide. Although similar to SARS-CoV and MERS-CoV, SARS-CoV-2 has resisted treatments that are effective against other coronaviruses. Crystal structures of two SARS-CoV-2 proteins, spike protein and main protease, have been reported and can serve as targets for studies in neutralizing this threat. We have employed molecular docking, molecular dynamics simulations, and machine learning to identify from a library of 26 million molecules possible candidate compounds that may attenuate or neutralize the effects of this virus. The viability of selected candidate compounds against SARS-CoV-2 was determined experimentally by biolayer interferometry and FRET-based activity protein assays along with virus-based assays. In the pseudovirus assay, imatinib and lapatinib had IC 50 values below 10 μM, while candesartan cilexetil had an IC 50 value of approximately 67 µM against M pro in a FRET-based activity assay. Comparatively, candesartan cilexetil had the highest selectivity index of all compounds tested as its half-maximal cytotoxicity concentration 50 (CC 50 ) value was the only one greater than the limit of the assay (>100 μM).

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Facet-dependent structure and dissociation of water at pristine IrO 2 /water interfaces

Understanding the microscopic structure of water at metal oxide interfaces is crucial for advancing electrocatalysis. IrO 2 , specifically, has shown exceptional activity for electrochemical water oxidation, but we currently lack a fundamental understanding of how the surface structure of IrO 2 impacts water reactivity. In this work, we developed a machine learning potential trained to first-principles accuracy for modeling IrO 2 /water interfaces across different facets: (110), (100), (101), and (001). Using extensive machine learning molecular dynamics simulations, we investigated the spontaneous dissociation of water molecules at these interfaces. Our results reveal a distinct dissociation probability trend: (110) > (100) ≈ (101) > (001), which we attribute primarily to the reaction thermodynamics of surface water dissociation. A strong correlation is observed between the surface Ir–O bond distances and the dissociation probabilities, highlighting the role of surface geometry in modulating reactivity. As a consequence, the interfacial solvation structures and hydrogen bonding environments are dynamically tuned by the varying water dissociation capabilities across facets. This work elucidates how water dissociation energetics depend on surface orientation and interfacial structure, offering atomistic insights into manipulating reaction chemistry at electrocatalytic interfaces.

organic↗
DOE OSTI2021

Traction-separation response of bilayer graphene interfaces: The role of an intercalated single atomic layer of water molecules and hydroxyl groups

Molecular dynamics simulations of the normal and shear traction-separation response of dry and wet bilayer graphene interfaces were performed. The influence of a single, atomic layer of water molecules sandwiched between graphene layers and hydroxyl functionalization on one or both layers were investigated. Hydroxyl groups provided a shielding effect during normal separation while favoring stick-slip friction under sliding conditions. The single, atomic layer of water molecules confined between hydrophilic and asymmetric hydrophobic/hydrophilic surfaces was found to ease sliding friction. The mechanisms arose from the balance between the attractive and repulsive graphene/water interactions through bridges of single water molecules between the two surfaces and the structure of the water monolayer under confinement for which an out-of-plane water/water hydrogen bond network could not be developed. The single, atomic layer of water molecules reduced the sliding friction the most when both surfaces were hydrophilic and lowered the sliding friction of the asymmetric hydrophobic/hydrophilic interfaces to that when both surfaces were hydrophobic. An optimal hydroxyl surface coverage of less than 10% may exist for efficient water lubrication, with diminishing effect with increasing number of hydroxyl groups. Finally, the findings offer a design strategy for assembling bilayer graphene interfaces with tailored friction and sliding contact properties.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Dynamics of poly[n]catenane melts

Inspired by advances in the chemical synthesis of interlocking polymer architectures, extensive molecular dynamics simulations have been conducted to study the dynamical properties of poly[n]catenanes-polymers composed entirely of interlocking rings-in the melt state. Both the degree of polymerization (number of links) and the number of beads per ring are systematically varied, and the results are compared to linear and ring polymers. Herein, a simple Rouse-like model is presented, and its analytical solution suggests a decomposition of the dynamics into "ring-like" and "linear-like" regimes at short and long times, respectively. In agreement with this picture, multiple sub-diffusive regimes are observed in the monomer mean-squared-displacements even though interchain entanglement is not prevalent in the system. However, the Rouse-type model does not account for the topological effects of the mechanical bonds, which significantly alter the dynamics at intermediate length scales both within the rings and at the chain segment scales. The stress relaxation in the system is extremely rapid and may be conveniently separated into ring-like and linear-like contributions, again in agreement with the Rouse picture. However, the viscosity has a non-monotonic dependence on the ring size for long chains, which disagrees strongly with theoretical predictions. This unexpected observation cannot be explained in terms of chain disentanglement and is inconsistent with other measures of polymer relaxation. Possible mechanisms for this behavior are proposed and implications for materials design are discussed.

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

Demonstration of an x-ray Raman spectroscopy setup to study warm dense carbon at the high energy density instrument of European XFEL

We present a proof-of-principle study demonstrating x-ray Raman Spectroscopy (XRS) from carbon samples at ambient conditions in conjunction with other common diagnostics to study warm dense matter, performed at the high energy density scientific instrument of the European x-ray Free Electron Laser (European XFEL). We obtain sufficient spectral resolution to identify the local structure and chemical bonding of diamond and graphite samples, using highly annealed pyrolytic graphite spectrometers. Due to the high crystal reflectivity and XFEL brightness, we obtain signal strengths that will enable accurate XRS measurements in upcoming pump–probe experiments with a high repetition-rate, where the samples will be pumped with high-power lasers. Molecular dynamics simulations based on density functional theory together with XRS simulations demonstrate the potential of this technique and show predictions for high-energy-density conditions. Our setup allows simultaneous implementation of several different diagnostic methods to reduce ambiguities in the analysis of the experimental results, which, for warm dense matter, often relies on simplifying model assumptions. The promising capabilities demonstrated here provide unprecedented insights into chemical and structural dynamics in warm dense matter states of light elements, including conditions similar to the interiors of planets, low-mass stars, and other celestial bodies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1997

Interactions of anesthetics with the water-hexane interface. A molecular dynamics study

The free energy profiles characterizing the transfer of nine solutes across the liquid-vapor interfaces of water and hexane and across the water-hexane interface were calculated from molecular dynamics simulations. Among the solutes were n-butane and three of its halogenated derivatives, as well as three halogenated cyclobutanes. The two remaining molecules, dichlorodifluoromethane and 1,2-dichloroperfluoroethane, belong to series of halo-substituted methanes and ethanes, described in previous studies (J. Chem. Phys. 1996, 104, 3760; Chem. Phys. 1996, 204, 337). Each series of molecules contains structurally similar compounds that differ greatly in anesthetic potency. The accuracy of the simulations was tested by comparing the calculated and the experimental free energies of solvation of all nine compounds in water and in hexane. In addition. the calculated and the measured surface excess concentrations of n-butane at the water liquid-vapor interface were compared. In all cases, good agreement with experimental results was found. At the water-hexane interface, the free energy profiles for polar molecules exhibited significant interfacial minima, whereas the profiles for nonpolar molecules did not. The existence of these minima was interpreted in terms of a balance between the free energy contribution arising from solute-solvent interactions and the work to form a cavity that accommodates the solute. These two contributions change monotonically, but oppositely, across the interface. The interfacial solubilities of the solutes, obtained from the free energy profiles, correlate very well with their anesthetic potencies. This is the case even when the Meyer-Overton hypothesis, which predicts a correlation between anesthetic potency and solubility in oil, fails.

NASA Center ARC↗
DOE OSTI2024

Comparing ab initio and quantum-kinetic approaches to electron transport in warm dense matter

Accurate knowledge of the electronic transport properties of warm dense matter is one of the main concerns of research in high-energy-density physics. Three modern approaches with vastly different levels of fidelity are reviewed and compared: the Kubo–Greenwood (KG) approach based on density-functional-theory molecular dynamics simulations (QMD), quantum kinetic theory based on average-atom models, and time-dependent density functional theory. Throughout, emphasis is placed on the connection between static properties of the electrons (e.g., density of states) and transport properties. Overall, it is found that whenever the conduction electrons can be modeled as being nearly free, fair to excellent agreement is found between QMD and kinetic theory approaches. Such a circumstance is required for modeling warm dense matter as a plasma of ions and free electrons, which is assumed in most kinetic theory approaches. The sensitivity of transport properties to the electronic structure is further highlighted by comparing different exchange–correlation approximations in QMD and KG calculations. It is found that the inclusion of exact exchange via thermal hybrid functionals can make a pronounced impact on electrical and thermal conduction in warm dense matter. We also investigate dynamic screening physics via kinetic theory and time-dependent density functional theory calculations of the mean free path of an electron in a hot dense plasma. In sum, we identify three axes along which to make progress in predicting electron transport in warm dense matter.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Dataset for manuscript "Consequences of the failure of equipartition for the p-V behavior of liquid water and the hydration free energy components of a small protein"

Previously, we showed that in the molecular dynamics simulation of a rigid model of water it is necessary to use an integration time-step dt that is less than or equal to 0.5 fs to ensure equipartition between translational and rotational modes. We extended that study in the NVT ensemble to NpT conditions and to an aqueous protein. We study neat liquid water with the rigid, SPC/E model and the protein BBA (PDB ID: 1FME) solvated in the rigid, TIP3P model. We examined integration time-steps ranging from 0.5 fs to 4.0 fs for various thermostat plus barostat combinations. We find that a small time-step, dt, is necessary to ensure consistent prediction of the simulation volume. Hydrogen mass repartitioning alleviates the problem somewhat, but is ineffective for the typical time-step used with this approach. The compressibility, a measure of volume fluctuations, is seen to be sensitive to dt. Using the mean volume estimated from the NpT simulation, we examined the electrostatic and van der Waals contribution to the hydration free energy of the protein in the NVT ensemble. These contributions are also sensitive to dt. In going from a time-step of 2 fs to a time-step of 0.5 fs, the change in the net electrostatic plus van der Waals contribution to the hydration of BBA is already in excess of the folding free energy reported for this protein. The data-set contains the simulation metadata and log files that support the claims noted above.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

High-pressure melt dynamics in shock-compressed titanium

In this work we study the high-pressure melting behavior of titanium using laser-driven shock compression with in situ femtosecond x-ray diffraction and molecular-dynamics simulations based on a machine-learned interatomic potential. The MD simulations predict the solid-liquid coexistence on the Hugoniot in the ∼111−124GPa range. Experimentally, we observe the first evidence of liquid at 86 GPa. We also observe pronounced microstructural changes with pressure, with strong grain refinement associated with the emergence of liquid, within the solid-liquid coexistence (∼110−126GPa). Above 126 GPa, we observe the persistence of residual levels of highly textured crystalline Ti to ∼180GPa, well above the expected melt completion pressure. We discuss the accuracy that current laser-shock experimental platforms have at determining the melt onset and completion pressures.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Anharmonicity and Phase Diagram of Magnesium Oxide in the Megabar Regime

With density functional molecular dynamics simulations, we computed the phase diagram of MgO from 50 to 2000 GPa up to 20000 K. Via thermodynamic integration (TDI), we derive the Gibbs free energies of the B1, B2, and liquid phases and determine their phase boundaries. With TDI and a pseudo-quasi-harmonic approach, we show that anharmonic effects are important and stabilize the B1 phase in particular. As a result, the B1-B2 transition boundary in the pressure-temperature plane exhibits a steep slope. We predict the B1-B2-liquid triple point to occur at approximately T= 10000 K and P = 370 GPa, which is higher in pressure than was inferred with quasi-harmonic methods alone. We predict the principal shock Hugoniot curve to enter the B2 phase stability domain but only over a very small range of parameters. Here, this may render it difficult to observe this phase with shock experiments because of kinetic effects.

3-dimensional systems↗
DOE OSTI2022

Cost-effective carbon fiber precursor selections of polyacrylonitrile-derived blend polymers: carbonization chemistry and structural characterizations

Blending polyacrylonitrile (PAN) with plastic wastes and bio-based polymers provides a convenient and inexpensive method to realize cost-effective carbon fiber (CF) precursors. In this work, PAN-based blend precursors are investigated using ReaxFF reactive molecular dynamics simulations with respect to the formation of all-carbon rings, the evolutions of oxygen-containing and nitrogen-containing species, and the migration of carbon atoms to form turbostratic graphene. From these simulations, we identify that PAN/cellulose (CL) blend manifests the highest carbon yield and the most substantial all-carbon ring formation. This ReaxFF-based finding is confirmed by Raman and TEM experiments indicating high crystallinity for PAN/CL-derived blend CFs. We trace the pathway of gasification and carbonization of PAN/CL to elaborate the mechanism of the formation of all-carbon ring networks. We discover that the acetals of CL can catalyze the cyclization of the blend precursor, allowing for the search for CL derivatives or the other kinds of bio-based polymers with similar functionalities as alternative blends. In addition, we examine the structural characteristics using the carbon–carbon (C–C) radial distribution functions, C–C bond length distributions, and sp2 C atom ratios for the four representative precursors, i.e., PAN, oxidized PAN, PAN/nylon 6,6, and PAN/CL. Our simulation results show the most extensive all-carbon ring cluster and graphitic structure growths for PAN/CL. Here, we propose PAN/CL as a cost-effective alternative CF precursor, since (a) CL is naturally abundant and eco-friendly for production, (b) the blend precursor PAN/CL does not require oxidation treatment, (c) PAN/CL has a high carbon yield with substantial all-carbon ring formation, and (d) PAN/CL based CFs potentially provide a mechanical property enhancement.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Elucidating the Role of Hydrogen at c‐In 2 O 3 /a‐In 2 O 3− x Interface

Ab initio molecular dynamics simulations and hybrid functional electronic structure calculations are employed to determine the formation, the structural and electronic properties, and the dynamics of covalent (In—OH) and ionic (In—H—In) hydrogen defects at crystalline-In 2 O 3 /amorphous-In 2 O 3−x interface. This comprehensive computational study considers (i) various interstitial and substitutional hydrogen site locations within the crystalline, amorphous, and interfacial regions; (ii) several oxygen-to-hydrogen ratios; and (iii) possible defect charged states. The results reveal hydrogen's inability to fully passivate the undercoordinated under-shared in atoms in amorphous highly substoichiometric oxide, giving rise to the formation of deep electron traps even in net-charge neutral structures. These trap defects are found to be sensitive to photoexcitation, in contrast to In—OH with H's electronic states located below the valence band and to In—H—In, where H is found to maintain its charge state upon illumination. Nevertheless, H plays a critical role in photoinduced conductivity and its relaxation by promoting In—O coordination transformations at the interfacial region, including deterioration of the crystalline layer. Finally, the results help identify mobile H species and metastable H defect complexes (such as In—H—H—In, In—OH–H—In, and In—OH–O—In) that are responsible for long relaxation times of the conductivity decay.

H defects↗
DOE OSTI2025

The Martini 3 Lipidome: Expanded and Refined Parameters Improve Lipid Phase Behavior

Lipid membranes are central to cellular life. Complementing experiments, computational modeling has been essential in unraveling complex lipid-biomolecule interactions, crucial in both academia and industry. The Martini model, a coarse-grained force field for efficient molecular dynamics simulations, is widely used to study membrane phenomena but has faced limitations, particularly in capturing realistic lipid phase behavior. Here, we present refined Martini 3 lipid models with a mapping scheme that distinguishes lipid tails that differ by just two carbon atoms, enhancing the structural resolution and thermodynamic accuracy of model membrane systems including ternary mixtures. The expanded Martini lipid library includes thousands of models, enabling simulations of complex and biologically relevant systems. These advancements establish Martini as a robust platform for lipid-based simulations across diverse fields.

Lipids↗
DOE OSTI2024

Shock Hugoniot calculations using on-the-fly machine learned force fields with ab initio accuracy

We present a framework for computing the shock Hugoniot using on-the-fly machine learned force field (MLFF) molecular dynamics simulations. In particular, we employ an MLFF model based on the kernel method and Bayesian linear regression to compute the free energy, atomic forces, and pressure, in conjunction with a linear regression model between the internal and free energies to compute the internal energy, with all training data generated from Kohn–Sham density functional theory (DFT). We verify the accuracy of the formalism by comparing the Hugoniot for carbon with recent Kohn–Sham DFT results in the literature. In so doing, we demonstrate that Kohn–Sham calculations for the Hugoniot can be accelerated by up to two orders of magnitude, while retaining ab initio accuracy. We apply this framework to calculate the Hugoniots of 14 materials in the FPEOS database, comprising 9 single elements and 5 compounds, between temperatures of 10 kK and 2 MK. We find good agreement with first principles results in the literature while providing tighter error bars. In addition, we confirm that the inter-element interaction in compounds decreases with temperature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Towards sharp error analysis of extended Lagrangian molecular dynamics

The extended Lagrangian molecular dynamics (XLMD) method provides a useful framework for reducing the computational cost of a class of molecular dynamics simulations with constrained latent variables. The XLMD method relaxes the constraints by introducing a fictitious mass ε for the latent variables and solving a set of singularly perturbed ordinary differential equations. While favorable numerical performance of XLMD has been demonstrated in several different contexts in the past decade, mathematical analysis of the method remains scarce. Here, we propose the first error analysis of the XLMD method in the context of a classical polarizable force field model. While the dynamics with respect to the atomic degrees of freedom are general and nonlinear, the key mathematical simplification of the polarizable force field model is that the constraints on the latent variables are given by a linear system of equations. We prove that when the initial value of the latent variables is compatible in a sense that we define, XLMD converges as the fictitious mass ε is made small with $\mathscr{O}$(ε) error for the atomic degrees of freedom and with $\mathscr{O}$($\sqrt{ε}$) error for the latent variables, when the dimension of the latent variable d' is 1. Furthermore, when the initial value of the latent variables is improved to be optimally compatible in a certain sense, we prove that the convergence rate can be improved to $\mathscr{O}$(ε) for the latent variables as well. Numerical results verify that both estimates are sharp not only for d'=1, but also for arbitrary d'. In the setting of general d', we do obtain convergence, but with the non-sharp rate of $\mathscr{O}$($\sqrt{ε}$) for both the atomic and latent variables.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Highly tail-asymmetric lipids interdigitate and cause bidirectional ordering

Phospholipids form structurally and compositionally diverse membranes. A less studied type of compositional diversity involves phospholipid tail variety. Some phospholipids contain two acyl tails which differ in length. These tail-asymmetric lipids are shown to contribute to temperature sensitivity, oxygen adaptability, and membrane fluidity. Membranes of a highly virulent intracellular bacterium, Francisella tularensis, contain highly tail-asymmetric 1-lignoceroyl-2-decanoyl-sn-glycero-3-phosphatidylethanolamine (XJPE) lipids which were previously shown to inhibit inflammatory responses in host cells. XJPE tails have unusually high asymmetry, and how they contribute to membrane properties on a molecular level is unknown. Here, we use small angle X-ray scattering and molecular dynamics simulations to investigate how varying XJPE ratios alters properties of simple membranes. Our results demonstrate that at high concentration they promote liquid-to-gel transition in otherwise liquid membranes, while at low concentration they are tolerated well, minimally altering membrane properties. In liquid membranes, XJPE lipids dynamically adopt two main conformations; with the long tail extended into the opposing leaflet or bent-back residing in its own leaflet. When added to both leaflets XJPE primarily adopts an extended confirmation, while asymmetric addition results in more bent-back orientations. The former increases tail ordering and the latter decreases it. XJPE tails adopt different conformations that induce composition- and leaflet-dependent bidirectional effect on membrane fluidity and this suggests that Francisella tularensis could use tail asymmetry to facilitate vesicle fusion and destabilize host cells. The effect of tail-asymmetric lipids on complex membranes should be further investigated to reveal the regulatory roles of high tail asymmetry.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1998

Folding and translocation of the undecamer of poly-L-leucine across the water-hexane interface. A molecular dynamics study

The undecamer of poly-L-leucine at the water-hexane interface is studied by molecular dynamics simulations. This represents a simple model relevant to folding and insertion of hydrophobic peptides into membranes. The peptide, initially placed in a random coil conformation on the aqueous side of the system, rapidly translocates toward the hexane phase and undergoes interfacial folding into an alpha-helix in the subsequent 36 ns. Folding is nonsequential and highly dynamic. The initially formed helical segment at the N-terminus of the undecamer becomes transiently broken and, subsequently, reforms before the remainder of the peptide folds from the C-terminus. The formation of intramolecular hydrogen bonds during the folding of the peptide is preceded by a dehydration of the participating polar groups, as they become immersed in hexane. Folding proceeds through a short-lived intermediate, a 3(10)-helix, which rapidly interconverts to an alpha-helix. Both helices contribute to the equilibrium ensemble of folded structures. The helical peptide is largely buried in hexane, yet remains adsorbed at the interface. Its preferred orientation is parallel to the interface, although the perpendicular arrangement with the N-terminus immersed in hexane is only slightly less favorable. In contrast, the reversed orientation is highly unfavorable, because it would require dehydration of C-terminus carbonyl groups that do not participate in intramolecular hydrogen bonding. For the same reason, the transfer of the undecamer from the interface to the bulk hexane is also unfavorable. The results suggest that hydrophobic peptides fold in the interfacial region and, simultaneously, translocate into the nonpolar side of the interface. It is further implied that peptide insertion into the membrane is accomplished by rotating from the parallel to the perpendicular orientation, most likely in such a way that the N-terminus penetrates the bilayer.

NASA Discipline Exobiology↗
DOE OSTI2025

The structure of liquid carbon elucidated by in situ X-ray diffraction

Carbon has a central role in biology and organic chemistry, and its solid allotropes provide the basis of much of our modern technology. However, the liquid form of carbon remains nearly uncharted, and the structure of liquid carbon and most of its physical properties are essentially unknown. But liquid carbon is relevant for modelling planetary interiors and the atmospheres of white dwarfs, as an intermediate state for the synthesis of advanced carbon materials, inertial confinement fusion implosions, hypervelocity impact events on carbon materials and our general understanding of structured fluids at extreme conditions. Here we present a precise structure measurement of liquid carbon at pressures of around 1 million atmospheres obtained by in situ X-ray diffraction at an X-ray free-electron laser. Our results show a complex fluid with transient bonding and approximately four nearest neighbours on average, in agreement with quantum molecular dynamics simulations. The obtained data substantiate the understanding of the liquid state of one of the most abundant elements in the universe and can test models of the melting line. The demonstrated experimental abilities open the path to performing similar studies of the structure of liquids composed of light elements at extreme conditions.

74 ATOMIC AND MOLECULAR PHYSICS↗