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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

CO 2 in Ionene–Ionic Liquid Composite Membranes

Abstract Ionene – ionic liquid (IL) composites are promising materials for CO 2 separation, yet a molecular‐level understanding of their structure and its impact on CO 2 speciation, solubility, rotation, and diffusivity remains unclear. Herein, using multimodal nuclear magnetic resonance (NMR), time‐of‐flight secondary ion mass spectrometry (ToF‐SIMS), atomic force microscopy (AFM), and molecular dynamics (MD) simulations, we reveal that the composites contain IL‐rich domains extending across hundreds of nanometres within the ionene matrix, and these bicontinuous domains span the entire membrane depth. CO 2 also absorbs into the ionene matrix, with the distribution between two CO 2 species varying with temperature and time. The rotational correlation times of these two species are on the timescale of 0.1 and 1 ns, respectively. As IL content increases, the ionic domains expand, resulting in higher CO 2 solubility due to enhanced molecular dynamics and increased free volume in both ionene backbones and IL‐rich regions. Although CO 2 diffusion in the membranes is an order of magnitude slower than in bulk IL, the activation energy for CO 2 diffusion remains comparable. Ionene‐IL composites represent a promising platform for designing CO 2 separation membranes, offering enhanced CO 2 diffusion and selectivity through IL‐rich domains, and increased CO 2 solubility and mechanical integrity from the ionene matrix.

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

Dynamic Monkey Bar Mechanism of Superionic Li‐ion Transport in LiTaCl 6

Abstract The LiTaCl 6 solid electrolyte has the lowest activation energy of ionic conduction at ambient conditions (0.165 eV), with a record high ionic conductivity for a ternary compound (11 mS cm −1 ). However, the mechanism has been unclear. We train machine‐learning force fields (MLFF) on ab initio molecular dynamics (AIMD) data on‐the‐fly and perform MLFF MD simulations of AIMD quality up to the nanosecond scale at the experimental temperatures, which allows us to predict accurate activation energy for Li‐ion diffusion (at 0.164 eV). Detailed analyses of trajectories and vibrational density of states show that the large‐amplitude vibrations of Cl − ions in TaCl 6 − enable the fast Li‐ion transport by allowing dynamic breaking and reforming of Li−Cl bonds across the space in between the TaCl 6 − octahedra. We term this process the dynamic‐monkey‐bar mechanism of superionic Li + transport which could aid the development of new solid electrolytes for all‐solid‐state lithium batteries.

Lei, Ming↗
DOE OSTI2024

Fluorophilic Sigma(σ)‐Lock Self‐Healable Copolymers

Abstract Although F‐Containing molecules and macromolecules are often used in molecular biology to increase the binding with Lewis acidic groups by introducing favorable C−F dipoles, there is virtually no experimental evidence and limited understanding of the nature of these interactions, especially their role in synthetic polymeric materials. These studies elucidate the molecular origin of inter‐ and intra‐Chain interactions responsible for self‐healing of F‐Containing copolymers composed of pentafluorostyrene and n‐butyl acrylate units (p(PFS/nBA). Guided by dynamic surface oscillating force (SOF) and spectroscopic measurements supported by molecular dynamics (MD) simulations, these studies show that the reformation of σ‐σ orbitals in −C−F of PFS and CH 3 CH 2 − of nBA units enables the recovery of entropic energy via fluorophilic‐σ‐lock van der Waals forces when PFS/nBA molar ratios are ~50/50. The strength of these interactions determined experimentally for self‐healable PFS/nBA compositions is in the order ~0.3 kcal/mol which primarily comes from fluorophilic‐σ‐lock (~70 %) contributions. These interactions are significantly diminished for non‐self‐healable counterparts. Strongly polarized −C−F σ orbitals create lateral dipolar forces enhancing the affinity towards −C−H orbitals, facilitating energetically favorable interactions. Entropic recovery driven by non‐Covalent bonding offers a valuable tool in designing materials with unique functionalities, particularly self‐healable batteries and energy storage devices.

Gaikwad, Samruddhi↗
DOE OSTI2023

Dynamic Monkey Bar Mechanism of Superionic Li–ion Transport in LiTaCl 6

Here, the LiTaCl 6 solid electrolyte has the lowest activation energy of ionic conduction at ambient conditions (0.165 eV), with a record high ionic conductivity for a ternary compound (11 mS cm –1 ). However, the mechanism has been unclear. We train machine-learning force fields (MLFF) on ab initio molecular dynamics (AIMD) data on-the-fly and perform MLFF MD simulations of AIMD quality up to the nanosecond scale at the experimental temperatures, which allows us to predict accurate activation energy for Li-ion diffusion (at 0.164 eV). Detailed analyses of trajectories and vibrational density of states show that the large-amplitude vibrations of Cl – ions in TaCl 6 – enable the fast Li-ion transport by allowing dynamic breaking and reforming of Li–Cl bonds across the space in between the TaCl 6 – octahedra. We term this process the dynamic-monkey-bar mechanism of superionic Li + transport which could aid the development of new solid electrolytes for all-solid-state lithium batteries.

25 ENERGY STORAGE↗
DOE OSTI2024

Fluorophilic Sigma(σ)‐Lock Self‐Healable Copolymers

Although F-Containing molecules and macromolecules are often used in molecular biology to increase the binding with Lewis acidic groups by introducing favorable C−F dipoles, there is virtually no experimental evidence and limited understanding of the nature of these interactions, especially their role in synthetic polymeric materials. These studies elucidate the molecular origin of inter- and intra-Chain interactions responsible for self-healing of F-Containing copolymers composed of pentafluorostyrene and n-butyl acrylate units (p(PFS/nBA). Guided by dynamic surface oscillating force (SOF) and spectroscopic measurements supported by molecular dynamics (MD) simulations, these studies show that the reformation of σ-σ orbitals in −C−F of PFS and CH 3 CH 2 − of nBA units enables the recovery of entropic energy via fluorophilic-σ-lock van der Waals forces when PFS/nBA molar ratios are ~50/50. The strength of these interactions determined experimentally for self-healable PFS/nBA compositions is in the order ~0.3 kcal/mol which primarily comes from fluorophilic-σ-lock (~70 %) contributions. These interactions are significantly diminished for non-self-healable counterparts. Strongly polarized −C−F σ orbitals create lateral dipolar forces enhancing the affinity towards −C−H orbitals, facilitating energetically favorable interactions. Entropic recovery driven by non-Covalent bonding offers a valuable tool in designing materials with unique functionalities, particularly self-healable batteries and energy storage devices.

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

Synthesis of Gold Hydride at High Pressure and High Temperature

Gold is an unreactive metal and its chemical interactions with hydrogen have only recently been explored. Here, in this study, we report the formation of gold hydride above 40 GPa and 2200 K in X-ray free electron laser heated diamond anvil cells using various hydrocarbons as hydrogen sources. Above 40 GPa, a hexagonal phase emerges close to the gold melting point, corresponding to a hydride with stoichiometry Au 2 H x , with x increasing from 0 to near 1 with pressure from 40 to 80 GPa. This is a high-temperature phase which reverts to face centered cubic gold on cooling to 295 K. Accompanying DFT-MD simulations are in excellent agreement with experiment and reveal the structure to consist of an hexagonal close packed gold lattice with atomic hydrogen disordered in the interstices. The hydrogen is superionic and exhibits high diffusivity through the crystalline gold lattice. Our results present the first solid-state binary compound of gold and hydrogen.

Ab initio calculations↗
DOE OSTI2025

Self‐Healable Poly(ionic liquid) Copolymers Driven by Polar and Dipolar Forces

Abstract Commodity aliphatic and aromatic acrylic‐based copolymers self‐heal due to ubiquitouskey‐and‐lock,ring‐and‐lock, andfluorophilic‐σ‐lockvan der Waals (vdW) interactions. However, the role of these interactions in the presence of covalently copolymerized ionic liquid (IL) is not known. This study is driven by the hypothesis that covalently incorporated cation–anion pairs to form poly(ionic liquid) copolymers (PILCs) can perturb inter‐ or intra‐chain vdW interactions reflected in mechanical and electrical responses. To test this hypothesis, we synthesized a series of PILCs comprising of pentafluorostyrene (PFS) and imidazolium‐based IL monomers with variable‐length aliphatic tails (methyl and butyl). Using a combination of 2D 1 H‐ 1 H and 19 F ‐ 19 F NOESY NMR and FTIR measurements supplemented by molecular dynamic (MD) simulations, these studies demonstrate that preferentially alternating/random PILCs topologies facilitate self‐healing. The introduction of cation–anion moieties modifies thefluorophilic‐σ‐lockinteractions and, along with longer aliphatic tails ─(CH 2 ) 3 CH 3 covalently attached to the imidazolium cation, enhances cation‐anion mobility, thus faster recovery from mechanical damage occurs. These findings underline how precise control over dipolar and ionic interactions through copolymer composition enables self‐healing in PILCs. These insights may open pathways for designing sustainable, mechanically resilient materials for applications in energy storage and energy harvesting.

Chemistry↗
DOE OSTI2022

Theoretical Investigation of the Oxygen Reduction Reaction over Platinum Catalysts Supported by Multi-Edged Vertically Aligned Carbon Nanofiber for Electrocatalyst Preparation

Vertically aligned carbon nanofibers (VACNFs) are promising supports for oxygen reduction reaction (ORR) electrocatalysts in fuel cells. Although experimentally these catalytic systems have shown great potential, there is lack of molecular understanding of the catalytic sites and reaction mechanisms. This work investigated the origin of the ORR reactivities of the platinum catalysts on multi-edged VACNFs (Pt/VACNF) using a multiscale modeling approach combining Density Functional Theory (DFT) and classical Molecular Dynamics (MD) simulations. Based on the ReaxFF potential, here all nanoscale Pt particles (Pt 55 , P 100 , and Pt 147 ) are stabilized by the open edges located axially along the VACNF walls. The calculated first-shell coordination numbers, CN, of surface Pt atoms are 6.63, 7.27, and 7.85, respectively, suggesting that the percentage of low coordination sites increases as the particle size decreases. The adsorption energies of OOH, O, and OH on Pt-55 were systematically probed using DFT calculations. These adsorption energies retain a linear correlation against the generalized coordination numbers ($\overline{CN}$). For Pt nanoparticles supported on VACNF, we found that the OOH and OH bind stronger than on Pt (111) by 0.14 and 0.17 eV, respectively, which can hinder the ORR activity with lower limiting potential than Pt (111). Our theoretical prediction is in good agreement with the linear sweeping voltammetry that revealed a left shift of the half-wave potential.

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

Structure- and Interaction-Based Design of Anti-SARS-CoV-2 Aptamers

Aptamer selection against novel infections is a complicated and time-consuming approach. Synergy can be achieved by using computational methods together with experimental procedures. In this study, we aim to develop a reliable methodology for a rational aptamer in silico et vitro design. The new approach combines multiple steps: (1) Molecular design, based on screening in a DNA aptamer library and directed mutagenesis to fit the protein tertiary structure; (2) 3D molecular modeling of the target; (3) Molecular docking of an aptamer with the protein; (4) Molecular dynamics (MD) simulations of the complexes; (5) Quantum-mechanical (QM) evaluation of the interactions between aptamer and target with further analysis; (6) Experimental verification at each cycle for structure and binding affinity by using small-angle X-ray scattering, cytometry, and fluorescence polarization. By using a new iterative design procedure, structure- and interaction-based drug design (SIBDD), a highly specific aptamer to the receptor-binding domain of the SARS-CoV-2 spike protein, was developed and validated. The SIBDD approach enhances speed of the high-affinity aptamers development from scratch, using a target protein structure. The method could be used to improve existing aptamers for stronger binding. This approach brings to an advanced level the development of novel affinity probes, functional nucleic acids. It offers a blueprint for the straightforward design of targeting molecules for new pathogen agents and emerging variants.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Direct Synthesis of Layer-Tunable and Transfer-Free Graphene on Device-Compatible Substrates Using Ion Implantation Toward Versatile Applications

Direct synthesis of layer-tunable and transfer-free graphene on technologically important substrates is highly valued for various electronics and device applications. State of the art in the field is currently a two-step process: a high-quality graphene layer synthesis on metal substrate through chemical vapor deposition (CVD) followed by delicate layer transfer onto device-relevant substrates. Here, we report a novel synthesis approach combining ion implantation for a precise graphene layer control and dual-metal smart Janus substrate for a diffusion-limiting graphene formation to directly synthesize large area, high quality, and layer-tunable graphene films on arbitrary substrates without the post-synthesis layer transfer process. Carbon (C) ion implantation was performed on Cu–Ni film deposited on a variety of device-relevant substrates. A well-controlled number of layers of graphene, primarily monolayer and bilayer, is precisely controlled by the equivalent fluence of the implanted C-atoms (1 monolayer ~4 × 10 15 C-atoms/cm 2 ). Upon thermal annealing to promote Cu-Ni alloying, the pre-implanted C-atoms in the Ni layer are pushed toward the Ni/substrate interface by the top Cu layer due to the poor C-solubility in Cu. As a result, the expelled C-atoms precipitate into a graphene structure at the interface facilitated by the Cu-like alloy catalysis. After removing the alloyed Cu-like surface layer, the layer-tunable graphene on the desired substrate is directly realized. The layer-selectivity, high quality, and uniformity of the graphene films are not only confirmed with detailed characterizations using a suite of surface analysis techniques but more importantly are successfully demonstrated by the excellent properties and performance of several devices directly fabricated from these graphene films. Molecular dynamics (MD) simulations using the reactive force field (ReaxFF) were performed to elucidate the graphene formation mechanisms in this novel synthesis approach. With the wide use of ion implantation technology in the microelectronics industry, this novel graphene synthesis approach with precise layer-tunability and transfer-free processing has the promise to advance efficient graphene-device manufacturing and expedite their versatile applications in many fields.

36 MATERIALS SCIENCE↗
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

Co-Assembly of Carbon Nanotube Porins into Biomimetic Peptoid Membranes

Robust and cost-effective membrane-based separations are essential to solving many global crises, such as the lack of clean water. Even though the current polymer-based membranes are widely used for separations, their performance and precision can be enhanced by using a biomimetic membrane architecture that consists of highly permeable and selective channels embedded in a universal membrane matrix. Researchers have shown that artificial water and ion channels, such as carbon nanotube porins (CNTPs), embedded in lipid membranes can deliver strong separation performance. However, their applications are limited by the relative fragility and low stability of the lipid matrix. In this work, we demonstrate that CNTPs can co-assemble into two-dimension (2D) peptoid membrane nanosheets, opening up a way to produce highly programmable synthetic membranes with superior crystallinity and robustness. A combination of molecular dynamics (MD) simulations, Raman spectroscopy, X-ray diffraction (XRD), and atomic force microscopy (AFM) measurements to verify the co-assembly of CNTP and peptoids are used and show that it does not disrupt peptoid monomer packing within the membrane. In conclusion, these results provide a new option for designing affordable artificial membranes and highly robust nanoporous solids.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Ionic Liquid‐Driven Modulation of DNA Brush Morphology on Nanoparticle Surfaces

Here, the morphology of DNA is strongly influenced by its surrounding environment, including factors such as pH, salt type and valency, and the presence of polymers. Inorganic salts are known to reduce the DNA chain length through mechanisms like electrostatic screening and ion bridging. In contrast, ionic liquids, a new class of organic salts, have previously been found to increase the DNA chain length, indicating a distinct mode of interaction between the ionic liquid and DNA chains. This study utilizes self-assembled DNA-AuNPs as a model system to examine changes in the DNA chain morphology and the nanoscale interaction mechanisms in an ionic liquid environment. The DNA chain lengths are measured in solution using X-ray scattering measurements at varying concentrations of two imidazolium ([$BMIM$] acetate and [$EMIM$] acetate) based ionic liquids. Additionally, Molecular Dynamics (MD) simulations are performed mimicking the experimental system. Our results suggest an interplay of electrostatic and groove-binding interactions governing the DNA chain morphology, which depends on IL concentration and the composition of the DNA chains. It has been found that for DNA chains with majority ssDNA, electrostatic interaction dominate, however with increasing composition of double strands, the DNA chains exhibit compaction due to a non-electrostatic hydrophobic groove-binding mechanism.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2022

Structure/property relationship of semicrystalline polymers during tensile deformation: a molecular dynamics approach

Molecular dynamics (MD) simulation with an ab initio numerical model is conducted to study the structure–property relationship of semicrystalline polymers during uniaxial deformation. The effects of chain length, temperature, and strain rate on mechanical properties are discussed. The influences of microstructural evolution such as bond length, entanglement density, and chain orientation are also studied quantitatively. The temperature will greatly affect the chain conformation in amorphous domains, and the results revealed that the interaction of amorphous and crystalline domains played a crucial role during stretching. Accordingly, the yielding of semicrystalline polymers follows different mechanisms at temperatures above and below T g . A melt-recrystallization scheme is observed during yielding at higher temperatures, while destruction of crystal structures is observed at lower temperatures at the yield point. The correlated effects of different temperatures and strain rates on mechanical properties are examined. This research is part of our efforts to develop a digital twin of a real experiment for efficient optimization of polymer material properties.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Relaxation dynamics of water in the vicinity of cellulose nanocrystals

Water molecules near cellulose nanocrystals (CNCs; produced via the sulfuric acid-catalyzed hydrolysis of wood pulp) are believed to relax slower than those in the bulk liquid, which may result in unique properties of CNC aqueous dispersions. Herein this study analyzed the polarization behavior of water molecules in CNC aqueous dispersions and other reference samples using a dielectric relaxation spectroscopy (DRS) technique in the microwave frequency range (0.2–20 GHz). As the CNC concentration increases, two slow relaxation components become prominent. The comparison with DRS data of aqueous dispersions of nanoporous silica, polyvinyl alcohol (PVA), and hairy CNC (HCNC) with amorphous chains protruding from both ends suggested that these slow relaxation modes of water near CNC surfaces cannot be attributed to direct hydrogen bonding interactions with the hydroxyl (OH) groups exposed and immobilized at the solid surface. Instead, they are similar to the water molecules interacting with OH groups attached to flexible polymer chains. Molecular dynamics (MD) simulations of the polarization behavior of water near the (110) facet of cellulose Iβ crystals confirmed that the interactions of water molecules with the cellulose crystal surface do not cause slower relaxations in the frequency range studied via the DRS. These results indicated that the CNC surface cannot be depicted with the crystallographic facets of cellulose Iβ; instead, it resembles a polymer-brush surface on which the short glucan residues or fragments of the strong acid-catalyzed hydrolysis process are swollen and extended into the aqueous phase.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Deciphering the non-linear impact of Al on chemical durability of silicate glass

The role of Al in aluminosilicate glasses remains somewhat a mystery: at low concentrations, it increases the resistance to hydrolysis of the glass, whereas at high concentrations an opposite effect is observed. To understand the origin of the phenomenon on a fundamental atomistic scale, we performed 577 MD simulations and applied potential mean force (PMF) calculations to estimate the activation barriers for hydrolysis and to statistically correlate them with local structural features of the glass. Models of pure silicate and aluminosilicate glasses are constructed and investigated. PMF simulation results are further validated by the experimental measurements and revealed that Al is very easy to dissociate, but it also increases the glass chemical durability through significantly increasing both the strength of Si and network connectivity of the glass. In contrast, at high Al concentration, preferential dissolution of Al weakens the silicate network, which it supposes to strengthen, and so the glass resistance becomes poor. Through PMF calculations, we evaluated the activation barriers for dissociating bonds around Al as 0.49 eV, which is less than a half of the energy to dissociate bonds around Si in pure silicate (1.22 eV) and around Si in aluminosilicate glass (1.34 eV), all these energy differences being statistically significant. Molecular structural level investigation revealed that Si with Al as a second neighbor in the glass network has a significantly higher activation energy for dissociation than Si in pure silicate glass. Furthermore, the proposed approach opens the way to the development of quantitative predictive models of glass durability.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Transition in helium bubble strengthening of copper from quasi-static to dynamic deformation

Damage from low-temperature irradiation and the subsequent degradation of materials performance pose significant challenges for the storage of radioactive materials and for peripheral components in some nuclear reactor designs. Fully understanding the mechanical behavior of such materials requires test data for strain rates in both the quasi-static (< 10/s) and dynamic (>> 10/s) regimes. While dynamic testing has generally been avoided in the past for neutron irradiated (contamination concerns) and ion irradiated (insufficient volume) materials, surface-sensitive Richtmyer-Meshkov instability (RMI) tests were used in the present work to overcome these limitations. Here, nanopillar compression, nanoindentation, and RMI testing data from a helium implanted surface layer (~10 µm thick) were compiled to explore the effects of helium bubbles on the materials strength of high-purity copper at strain rates of 0.001/s – 10 8 /s. While nano-mechanical testing revealed increases in yield strength and hardness with increasing helium dose from 1000 to 4000 appm He, RMI indicated no significant changes in strength as compared to unimplanted copper. Here, this discrepancy in behavior was rationalized through a combination of recent literature and follow-on molecular dynamics (MD) simulations, leading to the conclusion that the nanoscale helium bubbles acting as dispersed barriers to dislocation motion at quasi-static strain rates collapse under shock loading and cease to be effective barriers at high strain rates.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Evolution of dislocation loops in irradiated α-Uranium: An atomistically-informed cluster dynamics investigation

An atomistically-informed mean field cluster dynamics model has been developed to investigate the nucleation and growth of prismatic loops in irradiated a-Uranium. TEM analysis of neutron irradiated a-Uranium shows the evolution of self-interstitial atom and vacancy loops on (010) and (100) crystallographic planes, respectively, resulting in an anisotropic lattice swelling of its face-centered orthorhombic crystal. To provide model parameters, the crystallography of loops and the binding energy of point defects to these loops were studied using an angular dependent EAM potential and classical molecular dynamics (MD) simulations. Furthermore, using the bond-boost hyperdynamics method, the anisotropic diffusion of interstitials and vacancies in a-Uranium was investigated. Here, the mechanisms of point defect diffusion and the associated migration energies were reported and compared with previous DFT studies. The energetics and kinetic quantities mentioned above were used in the cluster dynamics model to predict the cluster density as a function of dose rate, dose and temperature and the results were compared to the reported neutron irradiation experiments. The model predictions reveal an accumulation of small sized vacancy loops along with a population of large and growing self-interstitial loops, which closely corresponds to the TEM observations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗