Search NASA⌕ Search

SEARCH · Search NASA

Results for “Molecular Dynamic simulations”

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

Shadow Molecular Dynamics and Atomic Cluster Expansions for Flexible Charge Models

Here, a shadow molecular dynamics scheme for flexible charge models is presented where the shadow Born–Oppenheimer potential is derived from a coarse-grained approximation of range-separated density functional theory. The interatomic potential, including the atomic electronegativities and the charge-independent short-range part of the potential and force terms, is modeled by the linear atomic cluster expansion (ACE), which provides a computationally efficient alternative to many machine learning methods. The shadow molecular dynamics scheme is based on extended Lagrangian (XL) Born–Oppenheimer molecular dynamics (BOMD). XL-BOMD provides stable dynamics while avoiding the costly computational overhead associated with solving an all-to-all system of equations, which normally is required to determine the relaxed electronic ground state prior to each force evaluation. To demonstrate the proposed shadow molecular dynamics scheme for flexible charge models using atomic cluster expansion, we emulate the dynamics generated from self-consistent charge density functional tight-binding (SCC-DFTB) theory using a second-order charge equilibration (QEq) model. The charge-independent potentials and electronegativities of the QEq model are trained for a supercell of uranium oxide (UO 2 ) and a molecular system of liquid water. The combined ACE+XL-QEq molecular dynamics simulations are stable over a wide range of temperatures both for the oxide and for the molecular systems and provide a precise sampling of the Born–Oppenheimer potential energy surfaces. Accurate ground Coulomb energies are produced by the ACE-based electronegativity model during an NVE simulation of UO 2 , predicted to be within 1 meV of those from SCC-DFTB on average during comparable simulations.

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

Ultrafast Events in Electrocyclic Ring-Opening Reactions

Electrocyclic reactions are characterized by the concerted formation and cleavage of multiple σ and π bonds in a molecular system and have been extensively studied since they were introduced by Robert Burns Woodward and Roald Hoffmann in 1965. Recent advances and the integration of time-resolved experiments and nonadiabatic quantum molecular dynamics simulations have transformed the traditional understanding of electrocyclic reactions beyond the Woodward–Hoffmann rules. In this review, we focus on recent studies of 1,3-cyclohexadiene and two of its derivatives, α-phellandrene and α-terpinene, to shed light on the underlying mechanisms of electrocyclic photochemical reactions. We highlight recent progress in ultrafast electron diffraction techniques and the simulation approach of ab initio multiple spawning. Together, these approaches can elucidate molecular structure dynamics from femtosecond to picosecond timescales as well as nuclear and electronic responses at conical intersections.

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

A model for trapping and re-solution regarding intra-granular bubbles in UO 2 , linked to atomic-scale simulations

In the literature, a clear definition of the irradiation re-solution frequency of gas from bubbles in the UO 2 fuel is absent. Moreover, for intra-granular bubbles, a detailed calculation of the cumulated displaced gas quantities in function of the distance from the radius of the bubble after a re-solution event has never been published. The assessment of these two elements is very useful if we want to increase the adherence of fission gas release codes to our present knowledge of the behavior of fission gases. Hence, we suggest to link the definition of the re-solution frequency to atomic-scale simulations. Furthermore, we present the cumulated displaced gas quantities obtained from Molecular Dynamics calculations, from which we have derived a re-solution profile that can be exploited to better consider the irradiation re-solution phenomenon inside Fission Gas Release codes. On top of that, we have built a new trapping/re-solution model for intra-granular bubbles linked to Molecular Dynamics simulations that can be easily incorporated into Fission Gas Release codes. In conclusion, we also check that the model is properly built through the comparison of the new model against a reference.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
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

Revisiting Theoretical Tools and Approaches for the Valorization of Recalcitrant Lignocellulosic Biomass to Value-Added Chemicals

Biorefinery processes for converting lignocellulosic biomass to fuels and chemicals proceed via an integrated series of steps. Biomass is first pretreated and deconstructed using chemical catalysts and/or enzymes to liberate sugar monomers and lignin fragments. Deconstruction is followed by a conversion step in which engineered host organisms assimilate the released sugar monomers and lignin fragments, and produce value-added fuels and chemicals. Over the past couple of decades, a significant amount of work has been done to develop innovative biomass deconstruction and conversion processes that efficiently solubilize biomass, separate lignin from the biomass, maximize yields of bioavailable sugars and lignin fragments and convert the majority of these carbon sources into fuels, commodity chemicals, and materials. Herein, we advocate that advanced in silico approaches provide a theoretical framework for developing efficient processes for lignocellulosic biomass valorization and maximizing yields of sugars and lignin fragments during deconstruction and fuel and chemical titers during conversion. This manuscript surveys the latest developments in lignocellulosic biomass valorization with special attention given to highlighting computational approaches used in process optimization for lignocellulose pretreatment; enzyme engineering for enhanced saccharification and delignification; and prediction of the genome modification necessary for desired pathway fine-tuning to upgrade products from biomass deconstruction into value-added products. Physics-based modeling approaches such as density functional theory calculations and molecular dynamics simulations have been most impactful in studies aimed at exploring the molecular level details of solvent-biomass interactions, reaction mechanisms occurring in biomass-solvent systems, and the catalytic mechanisms and engineering of enzymes involved in biomass degradation. More recently, with ever increasing amounts of data from, for example, advanced mutli-omics experiments, machine learning approaches have begun to make important contributions in synthetic biology and optimization of metabolic pathways for production of biofuels and chemicals.

09 BIOMASS FUELS↗
DOE OSTI2024

Library Screening, In Vivo Confirmation, and Structural and Bioinformatic Analysis of Pentapeptide Sequences as Substrates for Protein Farnesyltransferase

Protein farnesylation is a post-translational modification where a 15-carbon farnesyl isoprenoid is appended to the C-terminal end of a protein by farnesyltransferase (FTase). This process often causes proteins to associate with the membrane and participate in signal transduction pathways. The most common substrates of FTase are proteins that have C-terminal tetrapeptide CaaX box sequences where the cysteine is the site of modification. However, recent work has shown that five amino acid sequences can also be recognized, including the pentapeptides CMIIM and CSLMQ. In this work, peptide libraries were initially used to systematically vary the residues in those two parental sequences using an assay based on Matrix Assisted Laser Desorption Ionization–Mass Spectrometry (MALDI-MS). In addition, 192 pentapeptide sequences from the human proteome were screened using that assay to discover additional extended CaaaX-box motifs. Selected hits from that screening effort were rescreened using an in vivo yeast reporter protein assay. The X-ray crystal structure of CMIIM bound to FTase was also solved, showing that the C-terminal tripeptide of that sequence interacted with the enzyme in a similar manner as the C-terminal tripeptide of CVVM, suggesting that the tripeptide comprises a common structural element for substrate recognition in both tetrapeptide and pentapeptide sequences. Molecular dynamics simulation of CMIIM bound to FTase further shed light on the molecular interactions involved, showing that a putative catalytically competent Zn(II)-thiolate species was able to form. Bioinformatic predictions of tetrapeptide (CaaX-box) reactivity correlated well with the reactivity of pentapeptides obtained from in vivo analysis, reinforcing the importance of the C-terminal tripeptide motif. This analysis provides a structural framework for understanding the reactivity of extended CaaaX-box motifs and a method that may be useful for predicting the reactivity of additional FTase substrates bearing CaaaX-box sequences.

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

Chemical and Morphological Structure of Transgenic Switchgrass Organosolv Lignin Extracted by Ethanol, Tetrahydrofuran, and $\mathcal{γ}$-Valerolactone Pretreatments

We report the recalcitrance of lignocellulosic biomass is a challenge in biological-based biorefinery systems due to the complex physicochemical structure of plant cell walls. Pretreatment and genetic modification are two approaches in biomass conversion that have succeeded in modifying the structure of lignocellulose to enable better enzymatic deconstruction. However, the structural differences among pretreatment-solubilized lignin isolated from switchgrass genotypes have not been extensively investigated. Here, three organosolv pretreatment systems-ethanol (EtOH), tetrahydrofuran (THF), and γ-valerolactone (GVL)-were used on wild-type (WT) and two transgenic switchgrasses. All organosolv pretreatments caused a significant reduction in the molecular mass of lignins; particularly, up to ~90% decrease was observed in EtOH-pretreated lignin compared to untreated lignin. The WT EtOH lignin also presented the smallest particle size among all WT lignins. THF pretreated transgenic lignins showed a higher molecular mass, β-O-4 linkages, and aliphatic hydroxyl content compared to EtOH and GVL pretreated lignin. The number of hydrogen bonds between lignin and the organic solvents calculated from the molecular dynamics simulations followed the same trend as the experimentally determined reduction in lignin molecular mass. The results revealed the structural changes of solubilized lignin isolated from wild-type and transgenic switchgrass after different organosolv pretreatments.

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

Modeling the partitioning of amphiphilic molecules and co-solvents in biomembranes

We report amphiphilic co-solvents can have a significant impact on the structure, organization and physical properties of lipid bilayers. Describing the mutual impact of partitioning and induced structure changes is therefore a crucial consideration for a range of topics such as anesthesia and other pharmacokinetic effects, as well as microbial solvent tolerance in the production of biofuels and other fermentation products, where molecules such as ethanol, butanol or acetic acid might be generated. Small-angle neutron scattering (SANS) is a key method for studying lipid and polymer bilayer structures, with many models for extracting bilayer structure (thickness, area per lipid etc.) from scattering data in use today. However, the molecular details of co-solvent partitioning are conflated with induced changes to bilayer structure, making interpretation and modeling of the scattering curves a challenge with the existing set of models. To address this, a model of a bilayer structure is presented which invokes a two-term partition constant accounting for the localization of the co-solvent within the bilayer. This model was validated using a series of SANS measurements of lipid vesicles in the presence of the co-solvent tetrahydrofuran (THF), showing several strategies of how to deploy the two-parameter partition constant model to describe scattering data and extract both structure and partitioning information from the data. Molecular dynamics simulations are then used to evaluate assumptions of the model, provide additional molecular scale details and illustrate its complementary nature to the data fitting procedure. This approach results in estimates of the partition coefficient for THF in 1,2-dimyristoyl-sn-glycero-3-phosphocholine at 35°C, along with an estimate of the fraction of THF residing in the hydrophobic core of the membrane. The authors envision that this model will be applicable to a wide range of other bilayer/amphiphile interactions and provide the associated code needed to implement this model as a fitting algorithm for scattering data in the SasView suite.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

The importance of long-timescale simulations for driven systems: An example of He bubble growth at a W GB

Accelerated Molecular Dynamics (AMD) is used to study complex systems under realistic conditions by extending the timescales accessible by Molecular Dynamics. However, some studies rely instead on driving atomic systems harder with higher temperature, faster growth, etc. Here, we study He bubble growth at a W grain boundary as an illustration of harnessing AMD methods to avoid consequences of over-driving the system. The growth mechanisms observed for a He bubble grown under realistic conditions are compared to bubbles-grown orders of magnitude faster, at rates typical of conventional molecular dynamics simulations. We find that progressive growth mechanisms and bubble structures depend on the rate at which the bubble is grown providing further evidence that care must be taken when simulating the dynamics of driven systems such as this one.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Post-Marcus Theory and Simulation of Interfacial Charge Transfer Dynamics in Organic Semiconducting Materials (Final Report)

The research program developed, validated and applies predictive computational protocols for calculating charge transfer (CT) rates in complex molecular systems, including molecular dyads and triads in liquid solution and solid-state organic semiconducting (OSC) materials. We have established a transformative computational scheme that goes beyond widely used simplifications, to achieve realistic descriptions of CT processes. The approach properly addresses the contribution of molecular environment at ambient conditions to CT processes. Our approach achieves unique insight on investigated CT processes in relevant experimental efforts. The collaborative team included three principal investigators (PIs), with complimentary expertise in classical molecular dynamics simulations and data science (Cheung), state-of-the-art electronic structure calculations (Dunietz), and cutting-edge theory and simulation techniques for modeling energy, charge and coherence transfer dynamics in molecular systems (Geva).

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

Quantum Ornstein-Zernike theory for two-temperature two-component plasmas

Laboratory plasma production almost always preferentially heats either the ions or electrons, leading to a two-temperature state. In this state, density functional theory molecular dynamic simulation is the state of the art for modeling bulk material properties. We construct a statistical mechanics model for the two temperature limit that is theoretically consistent with the molecular dynamics method. We proceed to derive the electron-ion multi-temperature quantum Ornstein-Zernike equations for the first time. This allows the construction of a two-temperature two-component plasma model using the average atom from which we can compute bulk material properties at a fraction of the computation time of the two-temperature density functional theory simulation. The accuracy of the model is benchmarked against ion pair correlation and self-diffusion results from ab initio simulation. Here, we proceed to compute the viscosity and ion thermal conductivity as a function of both ion and electron temperature.

Ab initio molecular dynamics↗
DOE OSTI2021

SARS-CoV2 Protein-Ligand Simulation Dataset: Layer 2 (Lowest Temperature Extracted Protein Coordinate Trajectories)

Protein-only trajectories of the lowest-temperature window (310K) extracted from temperature replica-exchange molecular dynamics simulations of 23 different SARS CoV-2 systems, including S Protein ACE2-receptor binding-domain, MPro, PLPro, NSP3 ADRP (X-domain/macrodomain/phosphatase), NSP15 (endoribonuclease), NSP9, NSP10, NSP16, and N-protein N-terminus. Systems were prepared with charmm-gui and simulated using the GROMACS simulation software suite. Trajectories are provided in compressed dcd format with accompanying coordinate/topology files in pdb and psf formats. This data supplements the data release, DOI: 10.13139/OLCF/1650650 ('SARS-CoV2 Protein-Ligand Simulation Dataset: Layer 1 (Simulation Initial Conditions and Parameters)')

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

200-fold increase in dynamic strength of platinum at 430 GPa without a phase transformation

Here, the high-pressure strength of solid platinum is studied at the National Ignition Facility via Rayleigh-Taylor ripple growth under ramp compression. Laser-driven experiments reached pressures up to 430 GPa, collecting velocimetry and radiography data. Despite no phase transformation, the observed ripple growth aligns with hydrodynamic simulations using an amplified Steinberg-Guinan model, resulting in a 200-fold increase in strength. Molecular dynamics simulations of platinum ramp compression yield similar von Mises stress at comparable high pressure and high strain rate conditions. These findings represent the highest-pressure strength measurements ever achieved for platinum and offer critical insight into its significant strengthening behavior under high pressure and strain rate conditions.

Righi, Gaia [Lawrence Livermore National Laborator↗
DOE OSTI2024

Mean force emission theory for classical bremsstrahlung in strongly coupled plasmas

This work presents mean force emission theory, which extends the classical theory of bremsstrahlung emission to strongly coupled plasmas. In the high-frequency limit, the theory reduces to solving for the electron trajectory during a binary collision, but where the electron–ion interactions occur through the potential of mean force. In the low-frequency limit, it uses an autocorrelation formalism that captures effects of multiple collisions and strongly correlated motion. The predictions are benchmarked by comparison with first-principles classical molecular dynamics simulations of a fully ionized hydrogen plasma in which all interactions are repulsive. The comparison shows good agreement up to Coulomb coupling strengths of Γ∼30. The theory improves upon traditional models by including strong coupling effects and systematically including the effect of multiple collisions. Furthermore, mean force emission theory provides evidence that the Drude correction factor commonly used in quantum calculations of optical quantities may not be adequate at strong coupling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Ab Initio Direct Dynamics

The reactivity and dynamics of molecular systems can be explored computationally by classical trajectory calculations. The traditional approach involves fitting a functional form of a potential energy surface (PES) to the energies from a large number of electronic structure calculations and then integrating numerous trajectories on this fitted PES to model the molecular dynamics. The ever-decreasing cost of computing and continuing advances in computational chemistry software have made it possible to use electronic structure calculations directly in molecular dynamics simulations without first having to construct a fitted PES. In this “on-the-fly” approach, every time the energy and its derivatives are needed for the integration of the equations of motion, they are obtained directly from quantum chemical calculations. This approach started to become practical in the mid-1990s as a result of increased availability of inexpensive computer resources and improved computational chemistry software. The application of direct dynamics calculations has grown rapidly over the last 25 years and would require a lengthy review article. The present Account is limited to some of our contributions to methods development and various applications. To improve the efficiency of direct dynamics calculations, we developed a Hessian-based predictor-corrector algorithm for integrating classical trajectories. Hessian updating made this even more efficient. Furthermore, this approach was also used to improve algorithms for following the steepest descent reaction paths. For larger molecular systems, we developed an extended Lagrangian approach in which the electronic structure is propagated along with the molecular structure. Strong field chemistry is a rapidly growing area, and to improve the accuracy of molecular dynamics in intense laser fields, we included the time-varying electric field in a novel predictor-corrector trajectory integration algorithm. Since intense laser fields can excite and ionize molecules, we extended our studies to include electron dynamics. Specifically, we developed code for time-dependent configuration interaction electron dynamics to simulate strong field ionization by intense laser pulses. Our initial application of ab initio direct dynamics in 1994 was to CH 2 O → H 2 + CO; the calculated vibrational distributions in the products were in very good agreement with experiment. In the intervening years, we have used direct dynamics to explore energy partitioning in various dissociation reactions, unimolecular dissociations yielding three fragments, reactions with branching after the transition state, nonstatistical dynamics of chemically activated molecules, dynamics of molecular fragmentation by intense infrared laser pulses, selective activation of specific dissociation channels by aligned intense infrared laser fields, angular dependence of strong field ionization, and simulation of sequential double ionization.

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

Three-Dimensional Mass Spectrometric Imaging of Biological Structures Using a Vacuum-Compatible Microfluidic Device

Three-dimensional (3D) molecular imaging of biological structures is important for a wide range of research. In recent decades, secondary ion mass spectrometry (SIMS) has been recognized as a powerful technique for both two-dimensional (2D) and 3D molecular imaging. Sample fixations (e. g., chemical fixation and cryogenic fixation methods) are necessary to adapt biological samples to the vacuum condition in the SIMS chamber, which has been demonstrated to be non-trivial and less controllable, thus limiting the wider application of SIMS on 3D molecular analysis of biological samples. Our group recently developed in situ liquid SIMS that offers great opportunities for the molecular study of various liquids and liquid interfaces. In this work, we demonstrate that a further development of the vacuum-compatible microfluidic device used in in situ liquid SIMS provides a convenient freeze-fixation of biological samples and leads to more controllable and convenient 3D molecular imaging. The special design of this new vacuum-compatible liquid chamber allows an easy determination of sputter rates of ice, which is critical for calibrating the depth scale of frozen biological samples. Sputter yield of a 20 keV Ar 1800 + ion on ice has been determined as 1500 (± 8%) water molecules per Ar 1800 + ion, consistent with our results from molecular dynamics simulations. Moreover, using the information of ice sputter yield, we successfully conduct 3D molecular imaging of frozen homogenized milk and observe network structures of interesting organic and inorganic species. Finally, taken together, our results will significantly benefit various research fields relying on 3D molecular imaging of biological structures.

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

Self-Assembled Oligomers Facilitate Amino Acid-Driven CO 2 Capture at the Air–Aqueous Interface

Direct air capture of CO 2 using amino acid absorbents, such as glycine or sarcosine, is constrained by the relatively slow mass transfer of CO 2 through the air–aqueous interface. Our recent study showed a marked improvement in CO 2 capture by introducing CO 2 -permeable oligo-dimethylsiloxane (ODMS-MIM + ) oligomers with cationic (imidazolium, MIM + ) headgroups. Here, in this work, we have employed all-atom molecular dynamics simulations in combination with subensemble analysis using network theory to provide a detailed molecular picture of the behavior of CO 2 and the glycinate anions (Gly – ) at the ODMS-MIM + decorated air–aqueous interfaces. We show that the cationic head groups of the surfactants enhance the concentration and lifetime of Gly – in the interfacial region, while ODMS tails promote the physisorption of CO 2 in the interfacial region. Together, these two factors increase the effective region of contact and the probability of interactions between CO 2 and Gly – compared to that of the pure air–aqueous interface. The fundamental insights gained in this work establish essential foundations for developing hybrid systems with oligomer-decorated interfaces to maximize the overall CO 2 capture rates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗