Search NASA⌕ Search

SEARCH · Search NASA

Results for “MD 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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Computational studies reveal Fluorine based quinolines to be potent inhibitors for proteins involved in SARS-CoV-2 assembly

World is witnessing one of the worst pandemics of this century caused by SARS-CoV-2 virus which has affected millions of individuals. Despite rapid efforts to develop vaccines and drugs for COVID-19, the disease is still not under control. Chloroquine (CQ) and Hydroxychloroquine (HCQ) are two very promising inhibitors which have shown positive effect in combating the disease in preliminary experimental studies, but their use was reduced due to severe side-effects. Here, we performed a theoretical investigation of the same by studying the binding of the molecules with SARS-COV-2 Spike protein, the complex formed by Spike and ACE2 human receptor and a human serine protease TMPRSS2 which aids in cleavage of the Spike protein to initiate the viral activation in the body. Both the molecules had shown very good docking energies in the range of -6kcal/mol. Subsequently, we did a high throughput screening for other potential quinoline candidates which could be used as inhibitors. From the large pool of ligand candidates, we shortlisted the top three ligands (binding energy -8kcal/mol). We tested the stability of the docked complexes by running Molecular Dynamics (MD) simulations where we observed the stability of the quinoline analogues with the Spike-ACE2 and TMPRSS2 nevertheless the quinolines were not stable with the Spike protein alone. Thus, although the inhibitors bond very well with the protein molecules their intrinsic binding affinity depends on the protein dynamics. Moreover, the quinolines were stable when bound to electronegative pockets of Spike-ACE2 or TMPRSS2 but not with Viral Spike protein. We also observed that a Fluoride based compound: 3-[3-(Trifluoromethyl)phenyl]quinoline helps the inhibitor to bind with both Spike-ACE2 and TMPRSS2 with equal probability. The molecular details presented in this study would be very useful for developing quinoline based drugs for COVID-19 treatment.

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

How alkyl branching shapes structure in imidazolium and pyrrolidinium NTf 2 ionic liquids

High-energy X-ray scattering experiments and molecular dynamics (MD) simulations were carried out on ionic liquids (ILs) consisting of 1-alkyl-3-methylimidazolium and 1-alkyl-1-methylpyrrolidinium cations. These cations were paired with bis(trifluoromethylsulfonyl)amide anions and identical alkyl tails were used for both cationic species. The goal of this work is to investigate how the nanoscale structure of the ionic liquid changes with the length and with the degree of branching of the alkyl tail, for ILs having a common anion. We investigate spatial correlations in the intermolecular region, focusing on the intrinsic charge-charge interactions that characterize all ionic liquids, as well as the nanoscale domain segregation that is present in IL species with significant nonpolar components.

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

Irradiation-enhanced diffusion and diffusion-limited creep in U 3 Si 2

U 3 Si 2 is an advanced fuel candidate due to its relatively high fissile density and attractive thermal properties. Compared to standard UO 2 fuel, there are significant data gaps for the thermophysical and thermomechanical properties of U 3 Si 2 . Point defect concentrations and mobilities under irradiation govern a number of important fuel performance properties, such as creep and fission gas release. In this work, we utilized density functional theory (DFT) data to inform a cluster dynamics framework to predict point defect concentrations in U 3 Si 2 under irradiation. Molecular dynamics (MD) simulations were used to examine the contribution of atomic mixing during ballistic cascades to diffusion, as well as the diffusivity of U and Si at grain boundaries. These atomic scale models for diffusivity were then used to inform a creep model based on bulk (Nabarro-Herring) and grain boundary (Coble) diffusional creep, and climb-limited dislocation creep. The model compares well against available experimental data and has been implemented in the BISON fuel performance code. Finally, a demonstration case using simple power profiles has been carried out, showing that negligible creep occurs due to the low temperatures experienced by U 3 Si 2 in-reactor, a consequence of its high thermal conductivity.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Diffusion of small anti-Schottky clusters in UO 2

The aggregation of irradiation-induced defect clusters in UO 2 leads to the formation of dislocation loops and cavities which contribute swelling and ultimately produce deleterious effects on the fuel. Despite their fundamental role in fuel evolution, the kinetics of these defect clusters are not currently well understood. In this work, we investigate the diffusive behavior of interstitial clusters via Molecular Dynamics (MD) simulations. In this work, our investigation considers a range of defect cluster sizes, each composed of N UO 2 units or anti-Schottky defects. We report a complex cluster size - mobility relation; increasing cluster size corresponds to an increase in mobility up until a critical size (N=4) where diffusivity reaches a maximum, after which the trend reverts and further increasing size corresponds to a decrease in mobility. The rapid migration observed correlates well with the very low barriers reported in a few historic experimental studies. Further analysis of the cluster shape and orientation reveals a strong structural preference for near-planar configurations oriented normal to the $\langle$100$\rangle$ direction. Rotational energy barriers are found to be similar in both magnitude and trend to the observed diffusion barriers. Lastly, connections are drawn to the behavior of large defect clusters including the formation of (1/3) $\langle$111$\rangle$ Frank dislocation loops.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Effect of free surfaces on localized plastic deformation in single-crystal nickel containing helium bubbles and radiation-induced self-interstitial atom clusters

This paper presents molecular dynamics (MD) simulations of plastic deformation in single-crystal Ni models containing internal distributions of helium (He) bubbles and self-interstitial atom (SIA) clusters corresponding to experimental characterization of ex-service Inconel X-750 components. Here, we simulate single crystals with and without free surfaces at room temperature and 588K, the latter temperature corresponding to the Inconel X-750 service conditions. Both models are deformed under multiple slip to a maximum deviatoric strain of ~10 %. A dislocation slip band forms in the model without free surfaces, but not in the model with free surfaces. However, the latter exhibits nano-twin nucleation from surfaces. Helium (He) bubbles that overlap with a slip band or a nano-twinned region deform more severely than other bubbles in the model. Nevertheless, the bubbles do not link up and do not form configurations that resemble incipient cracks.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Molecular-dynamics study of diffusional creep in uranium mononitride

Uranium mononitride (UN) is a promising advanced nuclear fuel due to its high thermal conductivity and high fissile density. Yet, many aspects of its mechanical behavior and microstructural features are currently unknown. In this paper, molecular dynamics (MD) simulations are used to study UN's diffusional creep. Nanometer-sized polycrystals are used to simulate diffusional creep and to calculate an effective GB width. It is found that Nabarro-Herring creep is not dominant in the temperature range of 1700–2000 K and that the dominant diffusional creep mechanism is Coble creep with an activation energy of 2.28 ± 0.09 eV. A method is proposed to calculate the diffusional GB width and its temperature dependence in polycrystals. The effective GB width of UN is calculated as 2.69 ± 0.08 nm. This value fits very well with the prefactor of the phenomenological Coble creep formula. It is demonstrated that the most comprehensive thermal creep model for UN can be represented as the combination of our Coble creep model and the dislocation creep model proposed by Hayes et al.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Void growth via dislocation pileup impingement on grain boundary

In this work, the effect of edge dislocations, impinging on grain boundary in fcc Fe near a preexisting void on the volume of this void, was studied using MD simulation. It was found that virtually unlimited growth of the void via this mechanism is possible under certain conditions. To the best of our knowledge this is the first atomistic simulation study, demonstrating feasibility of this void growth mechanism.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Engineering the fracture resistance of 2H-transition metal dichalcogenides using vacancies: An in-silico investigation based on HRTEM images

Vacancy engineering of 2H-transition metal dichalcogenides (2H-TMDs) has recently attracted great attention due to its potential to fine-tune the phonon and opto-electric properties of these materials. From a mechanical perspective, this symmetry-breaking process typically reduces the overall crack resistance of the material and adversely affects its reliability. However, vacancies can trigger the formation of heterogeneous phases that synergistically improve fracture properties. In this study, using MoSe2 as an example, we characterize the types and density of vacancies that can emerge under electron irradiation and quantify their effect on fracture. Molecular dynamic (MD) simulations, employing a re-parameterized Tersoff potential capable of accurately capturing bond dissociation and structural phase changes, reveal that isolated transition metal monovacancies or chalcogenide divacancies tend to arrest the crack tip and hence enhance the monolayer toughness. In contrast, isolated chalcogenide monovacancies do not significantly affect toughness. The investigation further reveals that selenium vacancy lines, formed by high electron dose rates, alter the crack propagating direction and lead to multiple crack kinking. Using atomic displacements and virial stresses together with a continuum mapping, displacement, strain, and stress fields are computed to extract mechanistic information, e.g., conditions for crack kinking and size effects in fracture events. The study also reveals the potential of specific defect patterns, “vacancy engineering,” to improve the toughness of 2H-TMDs materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Effects of rejuvenation modes on the microstructures and mechanical properties of metallic glasses

This study systematically investigates the effects of three rejuvenation modes, namely biaxial loading, thermal pressure loading, and cryogenic thermal cycling, on the microstructures and mechanical properties of Cu 64 Zr 36 metallic glasses (MGs) using molecular dynamics (MD) simulations. The simulation results demonstrate that biaxial tension application during loading results in a rejuvenated glassy state with high strength. Specifically, we find that the yield stress, Young's modulus, and shear modulus of the rejuvenated glass obtained by biaxial tensile loading decreased by 2.8%, 2.1%, and 2.7%, respectively, compared to an as-cast MG reference. We also find that MGs achieve the greatest degree of rejuvenation, at the expense of strength, with the use of thermal pressure loading. Interestingly, the application of cryogenic thermal cycling results in a rejuvenated state with a lower strength drop compared to thermal pressure loading. Furthermore, analysis of the potential energy and medium-range-order (MRO) size indicates that the microstructures of MGs achieves a high-energy ordered glass state after thermal pressure loading, but a high-energy disordered glass state after biaxial tensile loading and cryogenic thermal cycling. Furthermore, these findings provide meaningful insights into the effect of rejuvenation modes on the microstructures and mechanical properties of MGs, as well as useful guidelines for designing MGs with high strength and superior plasticity for prospective applications in industry.

36 MATERIALS SCIENCE↗
DOE OSTI2025

An atomistic study connecting underlying dislocation behavior with superior mechanical properties of NiCoCr medium entropy alloy

NiCoCr-based medium-entropy alloy (MEA) with a simple face-centered cubic crystal phase exhibits excellent mechanical properties, often attributed to the synergy of multiple deformation mechanisms. However, the atomistic origin of their outstanding mechanical response, including microstructural evolution and dislocation behavior under varying strain-rates and orientation, remains unclear. In this work, we employ large-scale molecular dynamics (MD) simulations to investigate the changes in deformation mechanisms along three distinct orientations ([110], [111], [100]) under varying strain rates (1 ×10 8 /sec, 1 ×10 10 /sec, 1 ×10 12 /sec) in the NiCoCr MEA. The presence of the stair-rod and the Shockley partial dislocations under uniaxial tensile strain are found to play a key role in the formation of deformation twinning and ε-martensite, which positively correlates with strain-rate dependent dislocation analysis. These findings further establish the role of the dislocations in controlling the superior mechanical response and excellent fracture toughness of the NiCoCr MEA. Systematic transmission-electron microscopy tests performed on the [111]-oriented crystals, deformed at different strain levels, at room temperature provide clear evidence of both the extended stacking-fault and the stair rods, confirming the predicted microstructural features. Finally, this study offers key insights into the complex nucleation mechanisms of deformation twinning and ε-martensite, such as twinning – and transformation–induced plasticity (TWIP-TRIP), providing valuable guidelines for studying similar material classes.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Critical mineral substitutions in IN617: A combined computational and experimental approach to performance evaluation and feasibility

Addressing the escalating demand for critical minerals (CMs) driven by global climate change initiatives, this study explores compositional modifications to Inconel 617 (IN617) by substituting cobalt (Co) with manganese (Mn) across various atomic percentages. Here, we conducted a computational feasibility study employing Molecular Dynamics (MD) simulations to provide strategic guidance for experimental validation. The simulations analyzed tensile strength and corrosion resistance for five modified compositions (M1 to M5) to identify optimal properties. Tensile tests on cubic simulation cells were performed to generate stress-strain curves, revealing the impact of Co replacement with Mn on tensile strength—a metric correlated with hardness. Oxygen penetration simulations were conducted to evaluate corrosion resistance, indicating that reduced oxygen penetration depth corresponds to enhanced resistance. Promising compositions underwent phase diagram calculations for assessing phase stability. The optimal composition (M1), characterized by high tensile strength and minimal oxygen penetration, was chosen for experimental validation using induction melting and friction stir consolidation techniques. The materials were further characterized using SEM-EDS, XRD, and Vickers hardness testing. Our findings suggest that Mn substitution in IN617 can yield mechanical performance at par with high-Co alloys in energy-critical applications.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Weakly solvating ester electrolyte for high voltage sodium-ion batteries

Ethyl acetate (EA) was identified as a promising electrolyte solvent for sodium-ion batteries (SIBs), exhibiting low viscosity, cost-effectiveness, and low toxicity. Despite a significant portion of aggregation being linked to the weak solvation of Na + /EA as revealed by molecular dynamics (MD) simulations, pulsed-field gradient nuclear magnetic resonance (pfg-NMR) analysis identified a noteworthy Na+ diffusion coefficient of 3.95×10 -10 m 2 s -1 at 25°C in the presence of 1 m NaPF 6 salt. Employing fluoroethylene carbonate (FEC) as a film-forming additive to create electrode-electrolyte interphase, this electrolyte surprisingly made ~210 mAh Na 0.97 Ca 0.03 [Mn 0.39 Fe 0.31 Ni 0.22 Zn 0.08 ]O 2 (NCMFNZO)/hard carbon (HC) pouch cells achieve a lengthy cycling lifetime of 250 cycles with ~80 % capacity retention, cycled up to 4.0 V at 40°C. X-ray photoelectron spectroscopy (XPS) revealed increasing interphasial organic species over cycling, augmenting charge transfer resistance on both cathode and anode, particularly during fast charging or low temperatures (<10°C), promoting Na plating. Finally, gas chromatography-mass spectrometry combined with density functional theory identified CO 2 as the major gas generated from charged cathode/electrolyte interactions, exhibiting temperature/voltage dependence.

25 ENERGY STORAGE↗
DOE OSTI2021

Experimental and computational investigations of ethane and ethylene kinetics with copper oxide particles for Chemical Looping Combustion

In this work, reaction pathways for the oxidation of methane, ethane, and ethylene with CuO was obtained by ReaxFF Molecular Dynamics (MD) simulations between temperatures of 1000 K and 2000 K. Experiments in a fixed-bed flow reactor were preformed with methane, ethane, and ethylene at temperatures ranging from 500 K to 1000 K with time-dependent species measurements from an Electron-Ionization Molecular Beam Mass Spectrometer (MBMS), and species validation with Gas Chromatography (GC) for detection of complete and intermediate combustion products. The MBMS and GC allow for the detection of oxygenated species and larger species produced from radical reformation. The simulation and experiment agree on the production of such species as CH 3 CHO, CH 2 O, CO, and H 2 O, which allow for the creation of simple C1 and C2 reaction pathways, which can be used in kinetic models of C2 species and larger fuels such as biofuels, which inherently depend on C1 and C2 kinetics and reaction pathway. The simulation and experiment disagree on the formation of C 2 H 2 , CH 3 OH, and CO 2 with large amounts of C 2 H 2 being measured in the ethylene oxidation simulations and CH 3 OH being formed in methane oxidation simulations, while neither species were experimentally found. In the case of CO 2 large amounts of CO 2 are rapidly produced in experiments with C2 fuels at 800 K, while little-to-no CO 2 was observed in simulations. This is believed to be resulting from the extremely short timescale of the simulations, preventing total oxidation of the fuel. Here, the differences in products produced between simulation and experiment allow for the potential to modify the ReaxFF potential functions to more accurately model the experimental products of Cu–H–O–C reaction kinetics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Autonomous alloy composition optimization using molecular dynamics guided by a large language model

Here, we present an autonomous materials discovery framework that couples a large language model (LLM) with molecular dynamics (MD) simulations to optimize Fe–Cr–Mn alloy compositions for tensile strength. Starting from six distinct compositions, the LLM operated as an intelligent agent, iteratively proposing changes based on prior simulation results and constraints. Over 50 iterations per case, the LLM adaptively explored the composition space, identifying high-strength regions, not easily accessible by conventional methods. The highest strength, 18.7 GPa, was achieved with Fe 71 Cr 25 Mn 4 composition, identified from a Fe 75 Cr 20 Mn 5 starting point. The LLM autonomously adjusted its strategy in real time, demonstrating closed-loop decision-making using commodity hardware. This approach showcases the potential of LLMs as scientific co-pilots, capable of accelerating materials discovery and generalizable to other domains like biology and drug design.

Autonomy↗
DOE OSTI2024

Mapping structural and dynamic divergence across the MBOAT family

Membrane-bound O-acyltransferases (MBOATs) are membrane-embedded enzymes that catalyze acyl chain transfer to a diverse group of substrates, including lipids, small molecules, and proteins. MBOATs share a conserved structural core, despite wide-ranging functional specificity across both prokaryotes and eukaryotes. The structural basis of catalytic specificity, regulation and interactions with the surrounding environment remain uncertain. Here, we combine comparative molecular dynamics (MD) simulations with bioinformatics to assess molecular and interactional divergence across the family. In simulations, MBOATs differentially distort the bilayer depending on their substrate type. Additionally, we identify lipid binding sites surrounding reactant gates in the surrounding membrane. Complementary bioinformatic analyses reveal a conserved role for re-entrant loop-2 in MBOAT fold stabilization and a key hydrogen bond bridging DGAT1 dimerization. Finally, we predict differences in MBOAT solvation and water gating properties. These data are pertinent to the design of MBOAT-specific inhibitors that encompass dynamic information within cellular mimetic environments.

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

Modes of strain accommodation in Cu-Nb multilayered thin film on indentation and cyclic shear

Two-phase layered thin films with a high density of semi-coherent interfaces exhibit excellent mechanical properties and thermal stability. Here, in this study, a magnetron-sputtered Cu-Nb dual-layered thin film (~500 nm for Cu and ~150 nm for Nb) having an amorphous interface between Cu and Nb with a high density of aligned growth twins in Cu is subjected to severe surface deformation. The material is loaded using indentation and cyclic shear under tribological testing. The strain accommodation in the subsurface microstructure after deformation varies based on the local structure and deformation mode. Grain refinement and crack formations in the stressed region of the Nb layer and localized crystallization of the amorphous interface are observed after indentation and scratch testing. Pronounced detwinning of growth twins in the Cu layer under the cyclic shear strain leaves large dislocations sites and loops which are observed both by high-resolution transmission electron microscopy and experiment-guided molecular dynamic (MD) simulations. Our simulations provided insights into understanding the pathway for the detwinning process under cyclic shear loading.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Analysis of the molecular determinants for furin cleavage of the spike protein S1/S2 site in defined strains of the prototype coronavirus murine hepatitis virus (MHV)

We analyzed the spike protein S1/S2 cleavage of selected strains of a prototype coronavirus, mouse hepatitis virus (MHV) by the cellular protease furin, in order to understand the structural requirements underlying the sequence selectivity of the scissile segment. The probability of cleavage of selected MHV strains was first evaluated from furin cleavage scores predicted by the ProP computer software, and then cleavage was measured experimentally with a fluorogenic peptide cleavage assay consisting of S1/S2 peptide mimics and purified furin. We found that in vitro cleavability varied across MHV strains in line with predicted results—but with the notable exception of MHV-A59, which was not cleaved despite a high score predicted for its sequence. Using the known X-Ray structure of furin in complex with a substrate-like inhibitor as an initial structural reference, we carried out molecular dynamics (MD) simulations to learn the modes of binding of the peptides in the furin active site, and the suitability of the complex for initiation of the enzymatic cleavage. We identified the 3D structural requirements of the furin active site configuration that enable bound peptides to undergo cleavage, and the way in which the various strains tested experimentally are fulfilling these requirements. We find that despite some flexibility in the organization of the peptide bound to the active site of the enzyme, the presence of a histidine at P2 of MHV-A59 fails to properly orient the sidechain of His194 of the furin catalytic triad and therefore produces a distortion that renders the peptide/complex structural configuration in the active site incompatible with requirements for cleavage initiation. The Ser/Thr in P1 of MHV-2 and MHV-S has a similar effect of distorting the conformation of the furin active site residues produced by the elimination of the canonical salt-bridge formed by arginine in P1 position. This work informs a study of coronavirus infection and pathogenesis with respect to the function of the viral spike protein, and suggests an important process of viral adaptation and evolution within the spike S1/S2 structural loop.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Chloride Insertion Enhances the Electrochemical Oxidation of Iron Hydroxide Double-Layer Hydroxide into Oxyhydroxide in Alkaline Iron Batteries

Rechargeable alkaline iron batteries that constitute environmentally benign electrolytes and earth-abundant industrial materials are desirable green solutions for large-scale energy storage. As one of the most abundant metal elements in the earth’s crust, iron (Fe) can satisfy nearly all criteria for low-cost and safe battery electrodes. However, challenges in achieving reversible Fe redox impede their extensive implementation in modern energy supply systems. Here, this study revealed that Cl-anion insertion into Fe(OH) 2 layered double hydroxide (LDH) formed a green rust intermediate phase with the formula [Fe 2 2+ Fe 1 3+ (HO – ) 6 ] + [Cl] – , which assisted a high Fe(OH) 2 /FeOOH conversion reaction (64.7%) and improved cycling stability. This new iron redox chemistry was validated by operando X-ray diffraction, electrochemical testing, X-ray absorption spectroscopy (XAS), X-ray photoelectron spectroscopy (XPS) analysis, scanning transmission electron microscopy–energy-dispersive X-ray spectroscopy (STEM-EDS) mapping, and molecular dynamics (MD) simulations. Our study provides new insight into designing LDH materials for high-capacity alkaline iron batteries.

36 MATERIALS SCIENCE↗