Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Molecular 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,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Cryo-EM, Protein Engineering, and Simulation Enable the Development of Peptide Therapeutics against Acute Myeloid Leukemia

Cryogenic electron microscopy (cryo-EM) has emerged as a viable structural tool for molecular therapeutics development against human diseases. However, it remains a challenge to determine structures of proteins that are flexible and smaller than 30 kDa. The 11 kDa KIX domain of CREB-binding protein (CBP), a potential therapeutic target for acute myeloid leukemia and other cancers, is a protein which has defied structure-based inhibitor design. Here, we develop an experimental approach to overcome the size limitation by engineering a protein double-shell to sandwich the KIX domain between apoferritin as the inner shell and maltose-binding protein as the outer shell. To assist homogeneous orientations of the target, disulfide bonds are introduced at the target–apoferritin interface, resulting in a cryo-EM structure at 2.6 Å resolution. We used molecular dynamics simulations to design peptides that block the interaction of the KIX domain of CBP with the intrinsically disordered pKID domain of CREB. The double-shell design allows for fluorescence polarization assays confirming the binding between the KIX domain in the double-shell and these interacting peptides. Further cryo-EM analysis reveals a helix–helix interaction between a single KIX helix and the best peptide, providing a possible strategy for developments of next-generation inhibitors.

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

Large-Scale Atomistic Simulations: Investigating Free Expansion

The experiment investigates free expansion of a supercritical fluid into a two-phase liquid-vapor coexistence region. A huge molecular dynamics simulation (6 billion Lennard-Jones atoms) was run on 5760 GPUs (33% of LLNL Sierra) using LAMMPS/Kokkos software. This improved visualization workflow and started preliminary simulations of aluminum using SNAP machine learning potential.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Understanding Metal–Organic Framework Nucleation from a Solution with Evolving Graphs

A mechanistic understanding of metal–organic framework (MOF) synthesis and scale-up remains underexplored due to the complex nature of the interactions of their building blocks. In this work, we investigate the collective assembly of building units at the early stages of MOF nucleation, using MIL-101(Cr) as a prototypical example. Using large-scale molecular dynamics simulations, we observe that the choice of solvent (water and N,N-dimethylformamide), the introduction of ions (Na+ and F–) and the relative populations of MIL-101(Cr) half-secondary building unit (half-SBU) isomers have a strong influence on the cluster formation process. Additionally, the shape, size, nucleation and growth rates, crystallinity, and short and long-range order largely vary depending on the synthesis conditions. We evaluate these properties as they naturally emerge when interpreting the self-assembly of MOF nuclei as the time evolution of an undirected graph. Solution-induced conformational complexity and ionic concentration have a dramatic effect on the morphology of clusters emerging during assembly. While pure solvents lead to the rapid formation of a small number of large clusters, the presence of ions in aqueous solutions results in smaller clusters and slower nucleation. This diversity is captured by the key features of the graph representation. Principle component analysis on graph properties reveals that only a small number of molecular descriptors is needed to deconvolute MOF self-assembly. Furthermore, descriptors such as the average coordination number between half-SBUs and fractal dimension are of particular interest as they can be can be followed experimentally by techniques like by time-resolved spectroscopy. Ultimately, graph theory emerges as an approach that can be used to understand complex processes revealing molecular descriptors accessible by both simulation and experiment.

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

Extemporaneous Mechanochemistry: Shock-Wave-Induced Ultrafast Chemical Reactions Due to Intramolecular Strain Energy

We report regions of energy localization referred to as hotspots are known to govern shock initiation and the run-to-detonation in energetic materials. Mounting computational evidence points to accelerated chemistry in hotspots from large intramolecular strains induced via the interactions between the shock wave and microstructure. However, definite evidence mapping intramolecular strain to accelerated or altered chemical reactions has so far been elusive. From a large-scale reactive molecular dynamics simulation of the energetic material 1,3,5-triamino-2,4,6-trinitrobenzene, we map decomposition kinetics to molecular temperature and intramolecular strain energy prior to reaction. Both temperature and intramolecular strain are shown to accelerate chemical kinetics. A detailed analysis of the atomistic trajectory shows that intramolecular strain can induce a mechanochemical alteration of decomposition mechanisms. The results in this paper could inform continuum-level chemistry models to account for a wide range of mechanochemical effects.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Species Separation and Hydrogen Streaming upon Shock Release from Polystyrene under Inertial Confinement Fusion Conditions

Shock release from inertial confinement fusion (ICF) shells poses a great challenge to single-fluid hydrodynamic equations, especially for describing materials composed of different ion species. This has been evidenced by a recent experiment [Haberberger et al., Phys. Rev. Lett. 123, 235001 (2019)], in which low-density plasmas ( 10 19 to 10 20 cm – 3 ) are measured to move far ahead of what radiation-hydrodynamic simulations predict. To understand such experimental observations, in this work we have performed large-scale nonequilibrium molecular-dynamics simulations of shock release in polystyrene (CH) at experimental conditions. These simulations revealed that upon shock releasing from the back surface of a CH foil, hydrogen can stream out of the bulk of the foil due to its mass being lighter than carbon. This released hydrogen, exhibiting a much broader velocity distribution than carbon, forms low-density plasmas moving in nearly constant velocities ahead of the in-flight shell, which is in quantitative agreement with the experimental measurements. Such kinetic effect of species separation is currently missing in single-fluid radiation-hydrodynamics codes for ICF simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Spatial arrangement of dynamic surface species from solid-state NMR and machine learning-accelerated MD simulations

Here, the surface arrangement of motional organic functionalities is explored by experimental dipolar coupling measurements and the prediction of motionally-averaged coupling constant from molecular dynamics simulations. The use of machine learning potentials was key to reaching the timescale required. The distance between dynamic surface species are important in cooperative heterogeneous catalysis.

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

Isolating Chemical Reaction Mechanism as a Variable with Reactive Coarse-Grained Molecular Dynamics: Step-Growth versus Chain-Growth Polymerization

Here we present a general approach to isolate chemical reaction mechanism as an independently controllable variable across chemically distinct systems. Modern approaches to reduce the computational expense of molecular dynamics simulations often group multiple atoms into a single “coarse-grained” interaction site, which leads to a loss of chemical resolution. In this work we convert this shortcoming into a feature and use identical coarse-grained models to represent molecules that share nonreactive characteristics but react by different mechanisms. As a proof of concept, we use this approach to simulate and investigate distinct, yet similar, trifunctional isocyanurate resin formulations that polymerize by either chain- or step-growth. Because the underlying molecular mechanics of these models are identical, all emergent differences are a function of the reaction mechanism only. We find that the microscopic morphologies resemble related all-atom simulations and that simulated mechanical testing reasonably agrees with experiment.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Cooperativity in the Aldol Condensation Using Bifunctional Mesoporous Silica–Poly(styrene) MCM-41 Organic/Inorganic Hybrid Catalysts

This work explores the efficacy of silica/organic hybrid catalysts, where the organic component is built from linear aminopolymers appended to the silica support within the support mesopores. Specifically, the role of molecular weight and polymer chain composition in amine-bearing atom transfer radical polymerization-synthesized poly(styrene-co-2-(4-vinylbenzyl)isoindoline-1,3-dione) copolymers is probed in the aldol condensation of 4-nitrobenzaldehyde and acetone. Controlled polymerization produces protected amine-containing poly(styrene) chains of controlled molecular weight and dispersity, and a grafting-to thiol–ene coupling approach followed by a phthalimide deprotection step are used to covalently tether and activate the polymer hybrid catalysts prior to the catalytic reactions. Site-normalized batch kinetics are used to assess the role of polymer molecular weight and chain composition in the cooperative catalysis. Lower-molecular-weight copolymers are demonstrated to be more active than catalysts built from only molecular organic components or from higher-molecular-weight chains. Molecular dynamics simulations are used to probe the role of polymer flexibility and morphology, whereby it is determined that higher-molecular-weight hybrid structures result in congested pores that inhibit active site cooperativity and the diffusivity of reagents, thus resulting in lower rates during the reaction.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Probing C–I bond fission in the UV photochemistry of 2-iodothiophene with core-to-valence transient absorption spectroscopy

The UV photochemistry of small heteroaromatic molecules serves as a testbed for understanding fundamental photo-induced chemical transformations in moderately complex compounds, including isomerization, ring-opening, and molecular dissociation. Here, a combined experimental-theoretical study of 268 nm UV light-induced dynamics in 2-iodothiophene (C 4 H 3 IS) is performed. The dynamics are experimentally monitored with a femtosecond extreme ultraviolet (XUV) probe that measures iodine N-edge 4d core-to-valence transitions. Experiments are complemented by density functional theory calculations of both the pump-pulse induced valence excitations and the XUV probe-induced core-to-valence transitions. Possible intramolecular relaxation dynamics are investigated by ab initio molecular dynamics simulations. Furthermore, gradual absorption changes up to ~0.5 to 1 ps after excitation are observed for both the parent molecular species and emerging iodine fragments, with the latter appearing with a characteristic rise time of 160 ± 30 fs. Comparison of spectral intensities and energies with the calculations identifies an iodine dissociation pathway initiated by a predominant π → π* excitation. In contrast, initial excitation to a nearby n ⟂ → σ* state appears unlikely based on a significantly smaller oscillator strength and the absence of any corresponding XUV absorption signatures. Excitation to the π → π* state is followed by contraction of the C–I bond, enabling a nonadiabatic transition to a dissociative π→σ C-I * state. For the subsequent fragmentation, a relatively narrow bond-length region along the C–I stretch coordinate between 230 and 280 pm is identified, where the transition between the parent molecule and the thienyl radical + iodine atom products becomes prominent in the XUV spectrum due to rapid localization of two singly occupied molecular orbitals on the two fragments.

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

Efficient end-to-end simulation of time-dependent coherent X-ray scattering experiments

Physical optics simulations for beamlines and experiments allow users to test experiment feasibility and optimize beamline settings ahead of beam time in order to optimize valuable beam time at synchrotron light sources like NSLS-II. Further, such simulations also help to develop and test experimental data processing methods and software in advance. The Synchrotron Radiation Workshop ( SRW ) software package supports such complex simulations. We demonstrate how recent developments in SRW significantly improve the efficiency of physical optics simulations, such as end-to-end simulations of time-dependent X-ray photon correlation spectroscopy experiments with partially coherent undulator radiation (UR). The molecular dynamics simulation code LAMMPS was chosen to model the sample: a solution of silica nanoparticles in water at room temperature. Real-space distributions of nanoparticles produced by LAMMPS were imported into SRW and used to simulate scattering patterns of partially coherent hard X-ray UR from such a sample at the detector. The partially coherent UR illuminating the sample can be represented by a set of orthogonal coherent modes obtained by simulation of emission and propagation of this radiation through the coherent hard X-ray (CHX) scattering beamline followed by a coherent-mode decomposition. GPU acceleration is added for several key functions of SRW used in propagation from sample to detector, further improving the speed of the calculations. The accuracy of this simulation is benchmarked by comparison with experimental data.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Probing Carbon Mineralization Mechanisms in Pore and Bulk Fluids by Harnessing Architected Calcium Silicates

The ability to synthesize materials with well-controlled pore structures gives us unprecedented control over probing fluid interactions with reactive interfaces and advancing calibrated insights into coupled chemo-morphological interactions. One of the primary challenges in developing crystalline silicate materials lies in achieving ordered pore structures. Existing approaches of producing amorphous mesoporous metal silicates via sol–gel methods and heat-treatment of these materials to produce crystalline phases cause the pore structures in the amorphous phases to collapse. To overcome this challenge, carbon coating of amorphous mesoporous calcium silicate particles is carried out to retain the pore structure, while the material is heated to produce crystalline calcium silicate with calcium sulfate inclusions. The pore diameter in these materials is about 3.9 nm, with a surface area and a pore volume of 28.75 m 2 /g and 0.092 cm 2 /g, respectively. The mechanisms of carbon mineralization are investigated by reacting architected calcium silicate with 1 M Na 2 CO 3 and monitoring the evolution in the structural phases using operando wide-angle X-ray scattering (WAXS) measurements. Formation of stable calcium carbonate polymorph or calcite and metastable calcium carbonate polymorph or vaterite in pore and bulk fluids, respectively, resulting from the reaction between Na 2 CO 3 and CaSiO 3 , are noted. The mechanisms associated with carbon mineralization are delineated using ReaxFF molecular dynamics (MD) simulations. The surface dissolution reaction is initiated by 2H + ions that replace a Ca 2+ ion in the Ca–silicate matrix. Ca 2+ ions in the solution initially react with water to form calcium hydroxide and eventually form calcium (bi)carbonate. A slow and gradual increase in the formation of sodium silicate in the solution resulting from the reactions of silicic acid or the silicon dioxide reaction with sodium hydroxide is noted. When carbon mineralization occurs in environments bearing interfacial fluids, calcite is the dominant calcium carbonate polymorph, as determined using experiments with pore fluids and molecular-scale simulations. In conclusion, these studies provide fundamental insights into the mechanisms underlying the carbon mineralization of calcium silicate informed by experiments and molecular-scale simulations.

Calcium↗
DOE OSTI2022

Method to determine the electron–ion temperature relaxation rate from test particle distributions

Here a method to calculate the electron–ion energy exchange rate from the forces on and velocities of an ensemble of test particles is described. The essence of the method is that the energy exchange rate is related to the covariance between the distribution of velocities of test particles and the forces acting on them. The energy exchange rate is related to the electron–ion temperature relaxation rate in the limit of a low test particle speed. A proof of concept is conducted with first-principles molecular dynamics simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Structural and thermodynamic characterization of CuZr metallic glass nanoparticles: Insights from atomistic simulations

Metallic crystalline nanoparticles (NPs) have been shown to display intriguing size-dependent properties. However, the properties of metallic glass (MG) NPs remain largely unexplored. Here, using molecular dynamics simulations, we produce Cu 64 Zr 36 MG NPs ranging in size from 1 to 20 nm by cooling molten systems at a relatively slow rate of 10 9 K/s. Results indicate that NPs are coated with a Cu-rich layer with a thickness that increases with NP size. We employ various simulated rates of heating, from 10 11 to 10 13 K/s, to determine the melting points of these NPs. The results show a significant decrease in the solidus temperature for NPs smaller than 10 nm. Furthermore, as NP size decreases, the fraction of Cu full icosahedra within the NPs increases, suggesting that smaller NPs are stronger and harder. These findings provide insights for designing heterogeneous metallic nanoglass materials that leverage the size-dependent properties of MG NPs.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Large scale polymer toughening of two-dimensional materials revealed by in situ TEM fracture tests and multiscale simulations

Two-dimensional (2D) materials offer significant potential for applications in energy-harvesting devices, batteries, sensors, and transistors. However, their intrinsic brittleness makes them prone to mechanical failure, limiting their practical use. In this work, we perform in situ transmission electron microscopy (TEM) fracture tests on monolayer MoSe2 and uncover an extrinsic toughening effect induced by an ultrathin adsorbed polystyrene adlayer. This adlayer substantially enhances the fracture resistance of the 2D flakes. Through a combination of molecular dynamics simulations and finite element analysis, we elucidate the molecular mechanism behind this toughening effect. Further, it arises from the active crack-bridging behavior of entangled polymer chains and the formation of a fracture process zone that stabilizes crack propagation and increases the energy required for crack extension. The proposed toughening mechanism offers a pathway to improving the mechanical reliability of 2D material-based devices by mitigating the risk of sudden failure.

2D materials↗
DOE OSTI2025

Machine learning force field model for kinetic Monte Carlo simulations of itinerant Ising magnets

Here, we present a scalable machine learning (ML) framework for large-scale kinetic Monte Carlo (kMC) simulations of itinerant electron Ising systems. As the effective interactions between Ising spins in such itinerant magnets are mediated by conducting electrons, the calculation of energy change due to a local spin update requires solving an electronic structure problem. Such repeated electronic structure calculations could be overwhelmingly prohibitive for large systems. Assuming the locality principle, a convolutional neural network (CNN) model is developed to directly predict the effective local field and the corresponding energy change associated with a given spin update based on Ising configuration in a finite neighborhood. As the kernel size of the CNN is fixed at a constant, the model can be directly scalable to kMC simulations of large lattices. Our approach is reminiscent of the ML force field models widely used in first-principles molecular dynamics simulations. Applying our ML framework to a square-lattice double-exchange Ising model, we uncover unusual coarsening of ferromagnetic domains at low temperatures. Our work highlights the potential of ML methods for large-scale modeling of similar itinerant systems with discrete dynamical variables.

machine learning↗
DOE OSTI2023

Crystal structure of dihydrofolate reductase from the filarial nematode W. bancrofti in complex with NADPH and folate

Lymphatic filariasis is a debilitating illness with an estimated 50 million cases as of 2018. The majority of cases are caused by the parasitic worm W. bancrofti and additional cases by the worms B. malayi and B. timori. Dihydrofolate reductase (DHFR) is an established target in the treatment of cancer, bacterial, and protozoal infections and may be a potential target for drugs targeting parasitic worm infections, including filariasis. Recent studies have shown that known antifolate compounds, including methotrexate, inhibit the activity of W. bancrofti DHFR (WbDHFR). However, the absence of structural information for filarial DHFRs has limited the study of more in-depth structure-function relationships. We report the structure of WbDHFR complexed with NADPH and folate using X-ray diffraction data measured to 2.47 Å resolution. The structure of WbDHFR reveals the usual DHFR fold and is currently only the second nematode DHFR structure in the Protein Data Bank. The equilibrium dissociation constants for NADPH (90 ± 29 nM) and folate (23 ± 4 nM) were determined by equilibrium titrations. The interactions of known antifolates with WbDHFR were analyzed using molecular docking programs and molecular dynamics simulations. Antifolates with a hydrophobic core and extended linker formed favorable interactions with WbDHFR. These combined data should now facilitate the rational design of filarial DHFR inhibitors, which in turn can be used to determine whether DHFR is a viable drug target for filariasis and whether existing antifolates may be repurposed for its treatment.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

An analytic and complete equation of state for condensed phase materials

Analytic equations of state (EOS) are intended to reproduce theoretical and experimental data in a single phase portion of the thermodynamic space. We devise a complete and thermodynamically consistent model with four distinct features: (1) a reference isotherm that remains thermodynamically stable, (2) a flexible specific heat model based on a fourth-order rational polynomial, (3) a Grüneisen parameter that depends on specific volume and temperature, and (4) pressure and internal energy functions that can be inverted analytically in temperature. The model aims to improve the accuracy of existing equations of state while remaining computationally efficient. To demonstrate its features, we include calibrations for single-crystal pentaerythritol tetranitrate (PETN), liquid nitromethane (NM), and hexagonal close-packed beryllium (Be) metal. The parameter optimization uses the specific heat capacity, Grüneisen parameter, and static compression curves obtained from density functional theory for the crystalline solids and molecular dynamics simulations for liquid NM. We also present a velocity autocorrelation function that yields accurate phonon densities of states for the EOS calibration from the molecular dynamics trajectories. Each of the three calibrations is constrained to enforce the ambient state from experimental measurements and validated against experimental Hugoniot data from multiple sources. We also include one-dimensional hydrodynamic simulations of the isentropic compression experiments for beryllium conducted at the Z facility.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Atomistically-informed modeling of point defect clustering and evolution in irradiated ThO 2

A cluster dynamics (CD) model has been developed to investigate the nucleation and growth of point defect clusters, i.e., interstitial prismatic loops and nanoscale and sub-nanoscale voids, in ThO 2 during irradiation by energetic particles. The model considers cluster off-stoichiometry due to the asymmetry of point defect generation on the O and Th sublattices under irradiation, as well as the point defect diffusivities and the defect binding energies to clusters. The energies were established using detailed molecular dynamics simulations considering the statistical variability of cluster configuration. A high-order adaptive time-integration has been used to solve the model. The predicted loop density and their average size is in good agreement with reported experimental observations for proton irradiated ThO 2 at 600°C. In conclusion, the model did not predict void evolution due to the sluggish kinetics of cation vacancies, explaining the absence of voids in proton irradiated ThO 2 (and other oxides) at relatively low temperatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗