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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Atomistic simulations of nanoindentation on nanoglasses: Effects of grain size and gradient microstructure on the mechanical properties

The use of nanoglass (NG) microstructures has been shown as an effective strategy to improve the ductility of rather brittle metallic glass alloys. To evaluate the effects of grain size and gradient design on the mechanical properties of NGs, we perform molecular dynamics simulations of nanoindentation on Cu 64 Zr 36 NGs. We consider samples with uniform 3 and 7 nm grain sizes as well as gradient microstructures with gradient grain sizes varying from 3 to 7 nm. Here, the results show that the deformation mechanism in NGs with small grain sizes is dominated by the activation and evolution of multiple shear transformation zones. Increasing grain sizes enhances the elastic modulus and hardness at the cost of reduced abrasion resistance. Though the average grain size at the indentation surface plays a crucial role in the deformation behavior of gradient NG models, their strain localization and plastic deformation states are affected by the grain sizes far from the indenter. Desired combinations of mechanical properties can be realized by different NG microstructure designs.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Taming the Pseudoelastic Response of Nitinol Using Ion Implantation

Implantation of Ni 50.5 Ti 49.5 wire with 30 MeV Ni 6+ ions at doses (< 0.1 DPA) typically smaller than employed in the literature is shown to systematically alter the pseudoelastic response, with extrema in Berkovich nanoindentation load (+50%), hysteresis (–60%), and recoverable displacement (–19%) occurring at ~ 3.6 μm below the implantation surface. These extraordinary values are attributed to ~ 10 to 20 nm amorphous clusters that constrain the stress-induced B2-B19' phase transformation. This is substantiated by phase field simulations of crystalline-amorphous composites and molecular dynamics simulations of crystalline-vacancy cluster composites showing the spatial refinement of martensite caused by nm-scale defects. Finally, the results suggest that ion implantation may potentially expand the processing and performance space for NiTi, by creating amorphous defects at smaller length scales than dislocation substructures produced by conventional deformation processing.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Stress-Induced Structural Transformations in Au Nanocrystals

Nanocrystals can exist in multiply twinned structures like icosahedron or single crystalline structures like cuboctahedron. Transformations between these structures can proceed through diffusion or displacive motion. Experimental studies on nanocrystal structural transformations have focused on high-temperature diffusion-mediated processes. Limited experimental evidence of displacive motion exists. We report structural transformation of 6 nm Au nanocrystals under nonhydrostatic pressure of 7.7 GPa in a diamond anvil cell that is driven by displacive motion. X-ray diffraction and transmission electron microscopy were used to detect the structural transformation from multiply twinned to single crystalline. Single crystalline nanocrystals were recovered after unloading, then quickly reverted to the multiply twinned state after dispersion in toluene. Furthermore, the dynamics of recovery was captured using TEM which showed surface recrystallization and rapid twin boundary motion. Molecular dynamics simulations showed that twin boundaries are unstable due to defects nucleated from the interior of the nanocrystal.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

Surfactant-Specific AI-Driven Molecular Design: Integrating Generative Models, Predictive Modeling, and Reinforcement Learning for Tailored Surfactant Synthesis

Molecular design is a critical aspect of various scientific and industrial fields, where the properties of molecules hold significant importance. In this study, a 3-fold methodology design is presented that leverages the power of generative artificial intelligence (AI), predictive modeling, and reinforcement learning to create tailored molecules with desired properties. This model synergistically combines deep learning techniques with Self-Referencing Embedded Strings (SELFIES) molecular representation to build a generative model that generates valid molecules and a graphical neural network model that accurately forecasts molecular properties. The Variational Autoencoder (VAE) coupled with reinforcement learning helps refine molecule generation based on targeted attributes. Data from an experimental study involving surfactants were used to test the framework. A validation of the structural integrity of the molecules generated was conducted, and Tanimoto similarities were used to quantify the similarity and diversity between the original and generated molecular structures. Also, saliency maps for the generated surfactants were produced to identify the features explaining the property values. Lastly, molecular dynamics simulations were used to validate the stability of the generated molecules. The results showed that the proposed framework can effectively produce valid molecules within the set property threshold value.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Amino Acid Residues Controlling Domain Interaction and Interdomain Electron Transfer in Cellobiose Dehydrogenase

The function of cellobiose dehydrogenase (CDH) in biosensors, biofuel cells, and as a physiological redox partner of lytic polysaccharide monooxygenase (LPMO) is based on its role as an electron donor. Before donating electrons to LPMO or electrodes, an interdomain electron transfer from the catalytic FAD-containing dehydrogenase domain to the electron shuttling cytochrome domain of CDH is required. This study investigates the role of two crucial amino acids located at the dehydrogenase domain on domain interaction and interdomain electron transfer by structure-based engineering. The electron transfer kinetics of wild-type Myriococcum thermophilum CDH and its variants M309A, R698S, and M309A/R698S were analyzed by stopped-flow spectrophotometry and structural effects were studied by small-angle X-ray scattering. The data show that R698 is essential to pull the cytochrome domain close to the dehydrogenase domain and orient the heme propionate group towards the FAD, while M309 is an integral part of the electron transfer pathway – its mutation reducing the interdomain electron transfer 10-fold. Structural models and molecular dynamics simulations pinpoint the action of these two residues on the domain interaction and interdomain electron transfer.

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

The SARS-CoV-2 Spike Variant D614G Favors an Open Conformational State

The COVID-19 pandemic underwent a rapid transition with the emergence of a dominant viral variant (from the “D-form” to the “G-form”) that carried an amino acid substitution D614G in its “Spike” protein. The G-form is more infectious in vitro and associated with increased viral loads. To gain insight into the molecular-level underpinnings of these characteristics, we employed microsecond all-atom simulations. Here we show that changes in the protein energetics favor a higher population of infection-capable (open) states through release of hydrogen bonds of an asymmetry present in the D-form but not the G-form. Thus, the increased infectivity of the G-form is likely due to a higher rate of profitable binding encounters with the host receptor. It is also predicted to be more neutralization sensitive due to enhanced exposure of the receptor binding domain, a key target region for neutralizing antibodies. These results are significant for vaccine design. The Molecular Dynamics Simulations datasets generated in this study, namely that of the soluble form SARS-CoV-2 Spike protein are made available here. Details of the trajectory and supporting files are given in the README.txt file.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

eReaxFF force field development for BaZr 0.8 Y 0.2 O 3-δ solid oxide electrolysis cells applications

The use of solid-oxide materials in electrocatalysis applications, especially in hydrogen-evolution reactions, is promising. However, further improvements are warranted to overcome the fundamental bottlenecks to enhancing the performance of solid-oxide electrolysis cells (SOECs), which is directly linked to the more-refined fundamental understanding of complex physical and chemical phenomena and mass exchanges that take place at the surfaces and in the bulk of electrocatalysis materials. Here, we developed an eReaxFF force field for barium zirconate doped with 20 mol% of yttrium, BaZr 0.8 Y 0.2 O 3-δ (BZY20) to enable a systematic, large-length-scale, and longer-timescale atomistic simulation of solid-oxide electrocatalysis for hydrogen generation. All parameters for the eReaxFF were optimized to reproduce quantum-mechanical (QM) calculations on relevant condensed phase and cluster systems describing oxygen vacancies, vacancy migrations, electron localization, water adsorption, water splitting, and hydrogen generation on the surfaces of the BZY20 solid oxide. Using the developed force field, we performed both zero-voltage (excess electrons absent) and non-zero-voltage (excess electrons present) molecular dynamics simulations to observe water adsorption, water splitting, proton migration, oxygen-vacancy migrations, and eventual hydrogen-production reactions. Based on investigations offered in the present study, we conclude that the eReaxFF force field-based approach can enable computationally efficient simulations for electron conductivity, electron leakage, and other non-zero-voltage effects on the solid oxide materials using the explicit-electron concept. Moreover, we demonstrate how the eReaxFF force field-based atomistic-simulation approach can enhance our understanding of processes in SOEC applications and potentially other renewable-energy applications.

08 HYDROGEN↗
DOE OSTI2025

Dynamics and lipid membrane coupling of the RAS-RAF complex revealed via multiscale simulations

To gain molecular and mechanistic insights into initiation of the RAS-RAF signaling cascade, we developed and used a combination of multiscale simulation and experimental approaches. The influence and impact of the membrane on RAS and RAF proteins is a factor we are just beginning to understand and appreciate in more detail. Molecular simulation is an ideal methodology to further study this complicated relationship between the membrane and associated proteins. Our previous work using Multiscale Machine-learned Modeling Infrastructure investigated different lipid compositions solely around the KRAS4b protein and the interplay between protein behavior and these membrane environments. Multiscale Machine-learned Modeling Infrastructure uses machine learning to couple adjacent simulation scales and has been efficiently scaled across some of the world’s largest high-performance computers. Recently, we have expanded this multiresolution framework to include the all-atom simulation scale and to incorporate the RAF RBDCRD domains. Here, we present the overall analysis results from this new simulation campaign comprising a mixture of RAS and RAF RBDCRD proteins. Approximately 35,000 coarse-grained and 10,000 all-atom molecular dynamics simulations were completed, sampled from a variety of protein/lipid composition configurations that were generated from a micron-scale continuum simulation containing hundreds of copies of the proteins. Our studies suggest that orientations of the RAS-RBDCRD complex on the membrane occupy distinct configurational states, and the spatial patterns of lipid arrangements around these different protein states are unique to each state. The extent and size of lipid “fingerprints” imposed on the membrane by the RAS-RBDCRD protein complex are significantly larger than observed for just the RAS protein on its own. These protein complexes strongly associate, but we do not observe statistically significant preferred protein-protein orientations. These observations indicate that spatial colocalization of RAS-RBDCRD proteins in the same vicinity may be assisted by specific membrane environments, acting to increase the probability of signaling complex formation.

Carpenter, Timothy S. [Lawrence Livermore National↗
DOE OSTI2022

Imaging Wavepackets in Real-Space & Disentangling Ultrafast Solvation Dynamics with High-Energy Ultrafast X-ray Scattering

The microscopic information of solute and solute-solvent motions can be measured using ultrafast diffuse x-ray scattering, and compared to molecular dynamics simulations, as demonstrated in studies of solvated molecular systems done in X-ray free-electron lasers. However, for the typical photon energies used in these experiments, the photon momentum transfer range was limited to lower Q values where the signals arising from different parts of the studied system overlap significantly. Recently, this limitation was greatly alleviated by extending the photon energy at the Linac Coherent Light-Source (LCLS), significantly increasing the accessible momentum transfer range. This improvement is transformative for ultrafast diffuse x-ray scattering measurements, enabling for the first time to disentangle the solute scattering difference signal, that persists to higher Q values, from the bulk solvent and solute-solvent cross-terms difference signals. In addition, it enhances the prediction ability of modeling and simulation methods such as the hybrid QM/MM approach. The extended Q range also opens the way to resolve in real-space details regarding coherent wavepackets motions beyond their center-of-mass positions. Here, we present the first results on high-energy (18keV) diffuse x-ray scattering in solution, demonstrating high fidelity time-resolved scattering and analysis of the photoexcited model photocatalysts PtPOP ([Pt2(POP)4]4-), and IrDimen ([Ir2(dimen)4]2+) in several solvents. These complexes provide ideal systems for demonstrating the ability of high-Q ultrafast scattering and QM/MM simulations to decompose ultrafast solvation dynamics into specific changes in the solute-solvent pair distribution function with a particular focus on how electronically excited states change the interaction between the solvent and photo-catalytically active metal sites. We introduce an approach to further utilize the high-energy capability and demonstrate a single-shot ultra-wide-angle X-ray Scattering modality using two perpendicular detectors, spanning a scattering angle range of more than 100 degrees, allowing to extend the accessible momentum transfer range up to Q~14 Å-1. We develop a model-free method to invert the x-ray scattering signals and enable the recovery of multiple pair-density motions that happen simultaneously, allowing us to directly measure in real-space nuclear wavepacket motions. . [1] Natan, Adi. "Real-Space Inversion and Super-Resolution of Ultrafast X-ray Scattering using Natural Scattering Kernels." arXiv preprint arXiv:2107.05576 (2021)

Natan, Adi↗
DOE OSTI2024

Evaluating Cryo–TEM Reconstruction Accuracy of Self–Assembled Polymer Nanostructures

Cryogenic transmission electron microscopy (cryo–TEM) combined with single particle analysis (SPA) is an emerging imaging approach for soft materials. However, the accuracy of SPA–reconstructed nanostructures, particularly those formed by synthetic polymers, remains uncertain due to potential packing heterogeneity of the nanostructures. In this study, the combination of molecular dynamics (MD) simulations and image simulations is utilized to validate the accuracy of cryo–TEM 3D reconstructions of self–assembled polypeptoid fibril nanostructures. Using CryoSPARC software, image simulations, 2D classifications, ab initio reconstructions, and homogenous refinements are performed. By comparing the results with atomic models, the recovery of molecular details is assessed, heterogeneous structures are identified, and the influence of extraction location on the reconstructions is evaluated. In conclusion, these findings confirm the fidelity of single particle analysis in accurately resolving complex structural characteristics and heterogeneous structures, exhibiting its potential as a valuable tool for detailed structural analysis of synthetic polymers and soft materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Composition-dependent fracture energy in metallic glasses

The interplay between metallic glasses (MGs) mechanical properties, fracture energy (G), and glass-forming ability (GFA) and their dependence on alloy composition remain poorly understood. Here, we perform molecular dynamics simulations to investigate the intrinsic composition-dependence of G in Cu x Zr 100-x MGs (x = 20, 30, 40, 50, 64). The results indicate that the value of G increases with Cu content. In addition, it is revealed that MGs with higher G values display higher Poisson’s ratio (ν) and GFA, suggesting a close correlation between fracture toughness, mechanical properties, and GFA. This correlation between G, ν and GFA can be understood based on the fragility (m) of supercooled liquids, which is directly related to the structural heterogeneity in MGs. Larger m values are related to dynamic slowdown and supercooled liquid stabilization, which enhance GFA and the formation of pronounced structural heterogeneity, comprised of loosely packed regions that favor β-relaxation and the activation of shear transformation zones. Those concurrently promote the expansion of the plastic zone at the crack tip, enhancing the observed value of G. Furthermore, these simulation results shed light on the intrinsic relationship between fracture toughness, mechanical properties, and alloy composition in MGs.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Effect of Grain Boundary Misorientation on Spall Strength in Ta via Shock-Free Simulations with Relatively Few Atoms

A suite of 37 molecular dynamics simulations is conducted at two system sizes to systematically characterize the role of grain boundary (GB) misorientation on spall strength in pure BCC tantalum (Ta). The systems studied consist of bicrystals with a single [110] symmetric tilt grain boundary. Two loading conditions are compared: (i) homogeneous extension under uniaxial strain simulated in this study and (ii) piston/flyer impact of sample, which induces heterogeneous deformation via shockwave propagation along the length of the sample. The piston/flyer impact is taken from the literature and run on the same set of GB misorientation angles using LAMMPS. The major finding here is that both methods result in similar spall strength predictions, but the homogeneous extension method generally requires two to three orders of magnitude fewer atoms and similar reductions in computational costs. Spall strength results systematically overpredict using this method, by about 10% for the dataset three orders of magnitude smaller than piston/flyer simulations, and 5% for the dataset two orders of magnitude smaller. Lastly, the effect of system size and pre-compression magnitude on spall strength is systematically characterized.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Polyphenols Weaken Pea Protein Gel by Formation of Large Aggregates with Diminished Noncovalent Interactions

Polyphenols are well-known native cross-linkers and gel strengthening agents for many animal proteins. However, their role in modifying plant protein gels remains unclear. In this study, multiple techniques were applied to unravel the influence of green tea polyphenols (GTP) on pea protein gels and the underlying mechanisms. We found that the elasticity and viscosity of pea protein gels decreased with increased GTP. The protein backbone became less rigid when GTP was present based on shortened T 1ρ H in relaxation solid-state NMR measurements. Electron microscopy and small-angle X-ray scattering showed that gels weakened by GTP possessed disrupted networks with the presence of large protein aggregates. Solvent extraction and molecular dynamic simulation revealed a reduction in hydrophobic interactions and hydrogen bonds among proteins in gels containing GTP. The current findings may be applicable to other plant proteins for greater control of gel structures in the presence of polyphenols, expanding their utilization in food and biomedical applications.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Phase evolution in two-phase alloys during severe plastic deformation

Herein, phase evolution in FCC metals with strongly interacting alloy components during severe plastic deformation is investigated using molecular dynamics simulations. Specifically, we study the alloy microstructure in steady state, nucleation and growth of precipitates in supersaturated alloys, and the decomposition of precipitates in undersaturated alloys. The results are compared to a modified effective temperature model, providing a physical understanding for the atomic processes underlying the model and a perspective on its strengths and weaknesses. Key observations in this work are nucleation and growth of precipitates during SPD at a temperature of 100 K; Gibbs-Thomson-like behavior relating solubility to precipitate size under steady-state shearing; a direct relationship between the effective temperature and the shear modulus; and the importance of cluster agglomeration during precipitate growth. The study also reveals that the mechanisms of forced chemical mixing depends on precipitate size, adding complications for effective temperature models describing inhomogeneous systems. The simulations are shown to provide good semiquantitative agreement with experimental findings reported in the literature.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A molecular dynamics survey study of impact-induced amorphization in yttria nanoparticles

The impact of a single yttria (Y 2 O 3 ) nanoparticle onto a yttria substrate was studied as a function of particle velocity (300–1200 m/s) and diameter (5–12 nm) for three particle orientations and two substrate orientations using molecular dynamics simulations. Impacts resulted in minimal dislocation motion but extensive disordering of the crystal lattice without approaching the melting temperature. The resulting amorphous microstructure was both precisely coincident with plastic strain and retained until the end of the simulation, suggesting a deformation mechanism of amorphization-facilitated viscous flow. Further analysis revealed that the degree of amorphization depended primarily on the impact velocity and secondarily on the orientations of the particle and substrate, with significant second-order effects from variable interactions.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Probing ferroelectricity by x-ray absorption spectroscopy in molecular crystals

We carry out x-ray absorption spectroscopy experiment at the oxygen $K$ edge in croconic acid (C 5 H 2 O 5 ) crystal as a prototype of ferroelectric organic molecular solid, whose electric polarization is generated by proton transfer. The experimental spectrum is well reproduced by the electron-hole excitation theory simulations from configuration generated by ab initio molecular dynamics simulation. When inversion symmetry is broken in the ferroelectric state, the hydrogen bonding environment on the two bonded molecules become inequivalent. Such a difference is sensitively probed by the bound excitation in the pre-edge, which is strongly localized on the excited molecules. Our analysis shows that a satellite peak in the pre-edge will emerge at higher excitation energy, which serves as a clear signature of ferroelectricity in the material.

36 MATERIALS SCIENCE↗
DOE OSTI2021

First principles molecular dynamics maps out complete mineral surface acidity landscape

In “Interfacial structures and acidity constants (pKa) of goethite from first principles molecular dynamics simulations,” authors Y. Zhang, X. Lui, J. Cheng, and X. Lu apply First Principles molecular dynamics (FPMD, also called Density Functional Theory MD, DFT/MD, or ab initio MD, AIMD), to evaluable the complete set of acidity constants (pK a ) of the hydroxyl groups on the most prominent goethite crystal facets. The pK a of these OH and OH$^+_2$ groups are compared with available data from the multisite complexation (MUSIC) model traditionally used to estimate pK a on mineral surfaces. The authors have presented eloquent rationale for the importance and implications of understanding goethite acidity constants in room temperature geochemistry settings. In this paper, I focus on the computational aspects, the strengths of FPMD, and its possibilities.

58 GEOSCIENCES↗