Search NASA⌕ Search

SEARCH · Search NASA

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

Sticking Coefficients of Fusion Reactor Impurities from Molecular Dynamics Simulations for the Design of Cryopumps

A cryopump can be utilized as an impurity removal component of a direct internal recirculation (DIR) system for the fusion fuel cycle. The DIR facilitates a low fuel inventory by continuously pumping unburnt fuel while removing impurities from the fusion exhaust stream. A cryopump can target multiple impurity species by maintaining a temperature lower than the gas triple-point temperature that promotes desublimation. The desublimation/condensation of gases in cryopumps can be characterized by the sticking coefficient, which is defined as the probability for a gas particle to stick to a (cryo-)surface upon collision. The sticking coefficient is one of the important design/operation parameters for cryopumps, and it depends on a variety of surface and gas properties. Here, in this study, molecular dynamics simulations were utilized to estimate the sticking coefficients of typical fusion gas impurity species N 2 , CO 2 , and CH 4 over a Cu surface for a range of gas temperatures and surface coverages. The molecular dynamics study showed that the sticking coefficients for gases decrease with an increase in gas temperature. The presence of a single full monolayer of condensate on the metallic surface showed an adverse effect on the sticking of gases; however the sticking improved with two full monolayers of condensate on the surface. The sticking of gases over the mixed condensate on a surface was more favorable than the condensate of the same species for N 2 and CH 4 , with an exception for CO 2 , which showed a decrease in sticking over the mixed condensate.

cryopump↗
DOE OSTI2024

KCl-UCl 3 molten salts investigated by Ab Initio Molecular Dynamics (AIMD) simulations

Ab Initio Molecular Dynamics (AIMD) simulations are performed on molten KCl-UCl3 salt mixtures to determine energies, heat capacities, and densities. The density-dependent energy correction (DFT-dDsC), Grimme et al.'s DFT-D3, and Langreth & Lundqvist (vdW-cx) models are used for dispersion forces and combined with the Perdew-Burke-Ernzerhof (PBE) exchange-correlation potential with a Hubbard U parameter for the 5f electrons of uranium. After validating predictions for the end-member systems to literature data, KCl-UCl 3 mixtures are studied at select temperatures. Densities and energies both deviate from ideal solution behavior, with the maximum deviation occurring around 36% UCl 3 for mixing energies and slightly lower (29% UCl3) for densities. Compared to the NaCl-UCl 3 system, which was previously investigated using the same simulation methodologies, the KCl-UCl 3 density and mixing energy deviations from ideal solution behavior are larger by almost a factor of two. No deviation from ideal solution behavior for heat capacity was observed. The AIMD predictions for mixing energies and densities agree qualitatively with experimental data, though the spread in data obtained from the various dispersion force models utilized, measurements, and empirical estimates makes strong conclusions difficult. The dependence of thermodynamic and thermophysical properties on composition is correlated with the local chemistry of the solution phase, in particular, the tendency of UCl 3 to form network structures.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Machine Learning‐Augmented Molecular Dynamics Simulations (MD) Reveal Insights Into the Disconnect Between Affinity and Activation of ZTP Riboswitch Ligands

Abstract The challenge of targeting RNA with small molecules necessitates a better understanding of RNA–ligand interaction mechanisms. However, the dynamic nature of nucleic acids, their ligand‐induced stabilization, and how conformational changes influence gene expression pose significant difficulties for experimental investigation. This work employs a combination of computational and experimental methods to address these challenges. By integrating structure‐informed design, crystallography, and machine learning‐augmented all‐atom molecular dynamics simulations (MD), we synthesized, biophysically and biochemically characterized, and studied the dissociation of a library of small molecule activators of the 5‐aminoimidazole–4–carboxamide ribonucleotide triphosphate (ZTP) riboswitch, a ligand‐binding RNA motif that regulates bacterial gene expression. We uncovered key interaction mechanisms, revealing valuable insights into the role of ligand binding kinetics on riboswitch activation. Further, we established that ligand on‐rates determine activation potency as opposed to binding affinity and elucidated RNA structural differences, which provide mechanistic insights into the interplay of RNA structure on riboswitch activation.

Chemistry↗
DOE OSTI2025

Understanding Viscoelasticity of an Entangled Silicone Copolymer via Coarse-Grained Molecular Dynamics Simulations

Entangled dynamics is important for understanding rheological properties of long-chain polymers. For entangled homopolymers, the classic tube-reptation model and its refinements have been successfully applied to quantify properties like diffusion coefficient and zero-rate viscosity. However, the application of such models to copolymers has been limited despite scientific and industrial importance. Here, we study the entangled melt dynamics of poly(dimethyl-co-diphenyl)siloxane random copolymer for a range of mean-composition-ratio ϕ of the diphenyl component via long-term molecular dynamics simulation with a recently developed coarse-grained model. We found that the segmental relaxation is heterogeneous at the monomeric level because of compositional fluctuations. However, at the chain-entanglement level and higher length scales, the viscoelastic response is homogeneous with compositional dependence only through the overall diphenyl fraction ϕ. The relaxation modulus of the entangled copolymer melt conforms to the Likhtman–McLeish model, and the viscosity predicted using our current coarse-grained parameters is in good quantitative agreement with experimental data.

Copolymers↗
DOE OSTI2025

Deep potential molecular dynamics simulations of ion-enhanced etching of silicon by atomic chlorine

The continued development of plasma-assisted processing techniques requires a fundamental understanding of plasma-surface interactions. Molecular dynamics (MD) simulations have been employed to complement experimental studies and better understand the properties of such systems. Recently, machine learning (ML) methods have enabled the development of ab initio-based interatomic potentials, which can be generalized to complex combinations of multiple atom types. In this work, we use ML potentials developed using the Deep Potential Molecular Dynamics (DeepMD) framework to provide a model of ion-enhanced etching of Si by Cl atoms. We demonstrate the importance of proper selection of the training data set to the accuracy of the DeepMD model and compare our results to MD results using empirical potentials, as well as to experimental measurements. Exposure of undoped Si at 300 K to thermal Cl atoms yields a steady-state Cl coverage of 1.25 monolayers, which is slightly lower than the value obtained in previous experimental studies. Predictions of Si etch yields by simultaneous Cl atom and Ar + ion impacts as a function of ion energy, neutral to ion flux ratio, and angle of incidence of the ions are in reasonably good agreement with classical MD results and experimental measurements. Finally, etch yields and SiCl x mixed layer thicknesses during simultaneous bombardment of the Si(100) surface by Cl atoms and Cl + ions are in good agreement with experimental data. In conclusion, the present work is a necessary condition for the extension of the DeepMD procedure to more complex systems of interest in plasma-surface interactions.

Artificial neural networks↗
DOE OSTI2024

Microscale mechanical property variations of Al-substituted LLZO: insights from compression testing and molecular dynamics simulations

Ceramic solid electrolytes based on LLZO (Li 7 La 3 Zr 2 O 12 ) are promising candidates for all-solid-state batteries due to their high ionic conductivity and good apparent stability vs. lithium metal, however they are prone to mechanical failure. Lithium metal intrusions, alongside cell stack pressure, transition polycrystalline solid electrolyte grains into a compressed state that promotes crack propagation and fracture. Here this work examines the mechanical response of Al-substituted LLZO to compressive forces by measuring ultimate strength under pillar compression with a flat punch tip. Failure modes characterized by in situ scanning electron microscopy show diverse splitting patterns arising from internal porosity, grain boundaries, and slip planes. Large correlated variations in compressive strength (0.93–2.63 GPa) and Young's modulus (72.1–150.97 GPa) are observed across microscale regions of the solid electrolyte. Molecular dynamics simulations of LLZO with different porosities describe the variation of compressive strength and Young's modulus, and enable a microscale porosity model to be fit accounting for Young's modulus reduction across the solid electrolyte. Overall, the results indicate the importance of microscale mechanical testing of ceramic solid electrolytes to identify preferential sites for mechanical degradation and Li intrusion, and ensure the robust design of all-solid-state lithium metal batteries.

25 ENERGY STORAGE↗
DOE OSTI2025

Si–Cl 2 –Ar + Atomic Layer Etching Window: A Fundamental Study Using Molecular Dynamics Simulations and a Reduced Order Model

Silicon (Si) atomic layer etching (ALE) by alternating exposure to chlorine gas (Cl 2 ) and argon ions (Ar + ) is studied by using molecular dynamics (MD) simulations and a reduced order model (ROM). Here, the purpose of this study is to elucidate the properties of the ALE window, a range of ion energies where the amount of Si etched over a series of cycles is nonzero and nearly independent of ion energy. Experimental studies of the Si–Cl 2 –Ar + ALE system report contradictory results related to the ALE window’s ion energy range. Both MD simulations and the ROM show that there is an ALE window present from approximately 15 to 20 eV for normal incidence argon ions. The Si–Cl 2 –Ar + system, therefore, exhibits a narrow ALE window. The amount of Si etched per cycle is less than one atomic layer because of the higher etch yield of Cl atoms relative to atomic Si and silicon chlorides. A modified version of the ROM with an artificially increased Si physical sputtering threshold energy expands the ALE window, illustrating the importance of the difference in chemical and physical sputtering threshold energies in the ALE window energy range. The ROM is also used to examine the dependence of the EPC on the Ar + ion fluence.

energy↗
DOE OSTI2024

Molecular dynamics simulations of uranyl and plutonyl cations in a task-specific ionic liquid

Ionic liquids (ILs) are a unique class of solvents with potential applications in advanced separation technologies relevant to the nuclear industry. ILs are salts with low melting points and a wide range of tunable physical properties, such as viscosity, hydrophobiciy, conductivity, and liquidus range. ILs have negligible vapor pressure, are often non-flammable, and can have high thermal stability and a wide electrochemical window, making them attractive for use in separations processes relevant to the nuclear industry. Metal salts generally have a low solubility in ILs; however, by incorporating new functional groups into the IL cation or anion that promote complexation with the metal, the solubility can be greatly increased. One such task-specific ionic liquid (TSIL) is 1-carboxy-N, N, N-trimethylglycine bis(trifluoromethylsulfonyl)imide ([Hbet][Tf 2 N]). Water, which is detrimental for electrochemical separations, is a common impurity in ILs and can coordinate with actinyl cations, particularly in ILs containing only weakly coordinating components. Understanding the behavior of actinides in TSIL/water mixtures on a molecular level is vital for designing improved separations processes. Classical molecular dynamics simulations of uranyl(VI) and plutonyl(VI) in 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ([EMIM][Tf 2 N]) with deprotonated Hbet (betaine) and water have been performed to understand the coordination and dynamics of the actinyl cations. We find that betaine is a much stronger ligand than water and prefers to coordinate the metal in a bidentate manner. Potential of mean force simulations yield a relative free energy for betaine coordination of approximately -120 to -90 kJ/mol in mixtures with water. As the amount of betaine coordinated to the actinide increases, the diffusion coefficient of the actinyl cation decreases. Moreover, the betaine ligand is able to bridge between two metal centers, resulting in dimeric complexes with actinide–actinide distances of ~5 Å. Potential of mean force simulations show that these structures are stable, with relative free energies of up to -40 kJ/mol. The crystal structure for [(UO 2 ) 2 (bet) 6 (H 2 O) 2 ][Tf 2 N] 4 shows that the betaine bridges between two uranium atoms to form dimeric complexes similar to those found in our simulations.

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

Coarse-grained molecular dynamics simulation of solvent-dependent cellulose nanofiber interactions

Associations between cellulose are important both in biofuel production and in the use of cellulose for biomaterials. Cellulose nanofibers (CNFs) are sustainable, strong, light-weight alternatives to traditional materials in manufacturing, but are challenging to obtain due to irreversible aggregation in solution during preparative fibrillation. Therefore, it is imperative to understand the underlying factors driving aggregation with a view to designing solvents that can effectively compete with interfiber interactions, hence reducing aggregation. Molecular dynamics (MD) simulation at atomic detail can provide useful information on local interactions. However, the length and timescales accessible are too short to fully capture association processes. Here, we provide a method for accessing the longer length and timescales required using coarse-grained (CG) MD simulations with a MARTINI force field to calculate the interaction behavior of CNFs in three selected solvents: NaOH-urea-water, acetone, and neat water. The CG results are consistent with our prior all-atom MD and with previous experimental results. While acetone is found not to be an effective solvent, urea and ionic moieties in NaOH-urea-water not only solvate the fibrils but also improve the confinement of water molecules around them as shown by the solvent residence times and mean-square displacements. Overall, the presence of urea and ions reduces the likelihood of aggregation in multi-CNF systems relative to neat water irrespective of whether the hydrophobic or hydrophilic CNF surfaces are interacting. In conclusion, the CG method shows clear promise for selecting potential high-performance solvents for experimental prioritization in bioenergy and biomaterials research in a relatively fast manner as well as for understanding the aggregation and rheological behavior of CNF-solvent systems.

aggregation↗
NASA NTRS2018

Molecular Dynamics Simulations of Liquid and Polymer Electrolytes for Energy Storage Devices

Advancing beyond current lithium-ion technology is necessary in order to enable energy storage devices for electric airplanes. Electrolyte stability is a key limiting factor, yet the design of improved electrolytes remains a formidable challenge. Molecular dynamics (MD) simulations are a powerful tool for studying electrolytes, since they can be used to evaluate structural, thermodynamic, and transport properties, and can provide molecular-level detail often inaccessible to experimental techniques. Our computational materials groups at the NASA Ames Research Center has developed models and methods to accurately simulate both liquid and polymer electrolytes.We report the results from atomistic MD simulations of several electrolyte materials, with lithium salts dissolved in ionic liquids, dimethoxyethane (DME), and polyethylene oxide (PEO). For improved accuracy, we employ polarizable models, where each atom is given an environment-dependent atomic dipole. The simulations accurately predict bulk transport properties, including viscosity, diffusion, and ionic conductivity, in quantitative agreement with available experimental data. Moreover, the simulations provide important insights into the solvation structure of the lithium ions.We also report the results from coarse-grained MD simulations of polyanion electrolytes. In order to more efficiently capture the longer length- and time-scales of these systems, we employ a generic bead-spring model. These simulations provide important insight into how the polymer chain architecture and ionic interaction strengths affect the ionic aggregation behavior and cation dynamics. Despite the simplicity of the model, the simulations yield qualitative agreement with experimental data for similar systems.

Abbott, Lauren J.↗
NASA NTRSDate not supplied

Cardiovascular Responses to Simulated Spaceflight: Molecular Signatures and Surrogate Outputs to Measure CVD Risk

During extended space missions beyond low Earth orbit, astronauts will encounter prolonged periods of weightlessness and low dose space radiation. Previous studies have shown that exposure to small doses of high LET radiation (< 50 cGy) can lead to both short-term and long-term alterations in heart function, structure and underlying molecular mechanisms. In this study, we aim to identify the molecular signature associated with the cardiovascular response to simulated galactic cosmic radiation (5-ion GCR) alone or in combination with simulated weightlessness at time intervals relevant to mission length and recovery. Additionally, we aim to determine whether sex impacts cardiovascular responses to these spaceflight factors. Our overarching goal is to enhance our understanding of the cardiovascular risks associated with extended space missions and the clinical endpoints they suggest. We hypothesize that exposure to simulated space radiation leads to enduring alterations in the transcriptome, redox signaling and cytokine environment of cardiovascular tissue, some which have known links with reduced cardiovascular performance, aging, and increased risk of cardiovascular disease (CVD). Furthermore, we posit that simulated space radiation exposure in combination with simulated microgravity exacerbates cardiovascular deficits compared to single factor exposure. Female and male C57BL/6J mice, aged 23-24 weeks, were exposed to a single dose of 5, 15, or 50 cGy of 5-ion GCR, or sham-treated (0 cGy). Euthanasia was performed at 14 days and ~4 months post-irradiation. Hearts, aorta and blood plasma were collected shortly thereafter. RNA-sequencing of left ventricles at ~4 months post-GCR exposure revealed sex differences in the heart transcriptome with a few genes showing radiation-dependent changes in expression levels. Notably, some of the differentially expressed genes in 15 and 50 cGy GCR groups are known to play roles in the development of CVD. Analysis of protein levels of a subset of inflammatory cytokines in the heart indicated sex differences but no differences between sham and 50 cGy groups. Results also showed correlations among differentially expressed genes and a subset of inflammatory cytokines, with some correlations altered by GCR exposure. These findings suggest that GCR exposure can modify protein and gene networks linked to inflammation and CVD progression. In the aorta, telomere lengths were comparable across treatment groups sexes. Mitochondrial copy number is a biomarker for mitochondrial function with decreased copy numbers associated with cardiometabolic disease traits. Mitochondrial copy numbers of aorta also showed no sex nor dose differences. In a second study, mice underwent one week of simulated microgravity by hindlimb unloading (HU) and then exposed to a single dose of 15 cGy of 5-ion GCR. HU was conducted for an additional two weeks following GCR exposure. Single factor exposure groups (HU or GCR only) also were included in the study. Euthanasia was then performed and the same tissues were collected. Protein levels of select inflammatory cytokines in the heart showed sex-dependent differences in expression. In the aorta, telomere lengths and mitochondrial copy number also showed sex differences. In summary, our results indicate differences between sexes in biomarkers related to cardiovascular health. Exposure to 5-ion GCR or HU, alone or in combination, did not result in changes in most of the cardiovascular biomarkers that were examined. However, in the heart, simulated space radiation at doses of 15 and 50 cGy led to long-term alterations in the expression levels of a small group of genes known to be associated with the progression of CVD. The long-term transcriptomic changes resulting from exposure to simulated space radiation should be carefully investigated to mitigate adverse cardiovascular events during and after deep space missions. Our results also highlight the importance of sex-specific strategies in monitoring and maintaining cardiovascular health during and after deep space missions.

cardiovascular↗
NASA NTRS1997

Molecular Dynamics Simulation of a Multi-Walled Carbon Nanotube Based Gear

We used molecular dynamics to investigate the properties of a multi-walled carbon nanotube based gear. Previous work computationally suggested that molecular gears fashioned from (14,0) single-walled carbon nanotubes operate well at 50-100 gigahertz. The gears were formed from nanotubes with teeth added via a benzyne reaction known to occur with C60. A modified, parallelized version of Brenner's potential was used to model interatomic forces within each molecule. A Leonard-Jones 6-12 potential was used for forces between molecules. The gear in this study was based on the smallest multi-walled nanotube supported by some experimental evidence. Each gear was a (52,0) nanotube surrounding a (37,10) nanotube with approximate 20.4 and 16,8 A radii respectively. These sizes were chosen to be consistent with inter-tube spacing observed by and were slightly larger than graphite inter-layer spacings. The benzyne teeth were attached via 2+4 cycloaddition to exterior of the (52,0) tube. 2+4 bonds were used rather than the 2+2 bonds observed by Hoke since 2+4 bonds are preferred by naphthalene and quantum calculations by Jaffe suggest that 2+4 bonds are preferred on carbon nanotubes of sufficient diameter. One gear was 'powered' by forcing the atoms near the end of the outside buckytube to rotate to simulate a motor. A second gear was allowed to rotate by keeping the atoms near the end of its outside buckytube on a cylinder. The ends of both gears were constrained to stay in an approximately constant position relative to each other, simulating a casing, to insure that the gear teeth meshed. The stiff meshing aromatic gear teeth transferred angular momentum from the powered gear to the driven gear. The simulation was performed in a vacuum and with a software thermostat. Preliminary results suggest that the powered gear had trouble turning the driven gear without slip. The larger radius and greater mass of these gears relative to the (14,0) gears previously studied requires a smaller rotation rate and multiple rows of teeth to avoid excessive force on the gear teeth resulting, in slip and failure of the driven gear to turn. We hope that studies such as these will eventually lead to synthesis of components that can be assembled into atomically precise fullerene machines. These machines, in turn, may someday be used in machine-phase fullerene materials with remarkable properties.

Han, Jie↗
DOE OSTI2025

From Electronic Structure to Ion Transport: Photoelectron Spectroscopy and Molecular Dynamics Simulations Reveal the Role of Anions in Lithium Battery Electrolytes

Electrolyte anions are pivotal for lithium battery performance, yet their fundamental electronic structural properties are not well understood. In this work, we employ a combination of negative-ion photoelectron spectroscopy (NIPES), ab initio calculations, and molecular dynamics (MD) simulations to investigate the electronic structures of three representative electrolyte anions. This multiscale approach enables us to elucidate how their intrinsic electronic properties govern anion–solvent interactions in gas-phase clusters, as well as lithium-ion (Li + ) solvation structures and ion transport behavior in the condensed phase. NIPES reveals that difluoro(oxalato)borate (DFOB – ), bis(fluorosulfonyl)imide (FSI – ), and bis(oxalato)borate (BOB – ) all exhibit high electron binding energies, with vertical/adiabatic detachment energies increasing from DFOB – (6.09/5.70 eV) to FSI – (6.80/6.10 eV) to BOB – (6.82/6.40 eV), correlating with enhanced oxidation stability. Ab initio calculations reveal that DFOB – /FSI – –solvent complexes bind Li + ∼ 10 kcal/mol stronger than BOB – series, aligning with the strength of a Li + –anion model. DFOB – exhibits pronounced charge localization on both oxygen and fluorine atoms, enabling their involvement in Li + coordination. In contrast, fluorine atoms in FSI – are largely electron-depleted and remain excluded from direct Li + binding. MD simulations further demonstrate that LiDFOB and LiFSI systems exhibit Li + diffusion coefficients three and five times higher than those of LiBOB across four common solvents. Notably, LiFSI salt in acetonitrile (AN) exhibits the fastest Li + diffusion among 12 electrolyte systems, highlighting the synergistic effect of FSI – and AN in promoting ion mobility. In conclusion, these findings provide a molecular-level understanding of the critical roles of anion and its microsolvation in optimizing Li + diffusion dynamics, once again emphasizing the positioning of FSI – and DFOB – as prime candidates for next-generation electrolytes.

25 ENERGY STORAGE↗
DOE OSTI2025

Molecular Dynamics Simulations of Supercritical Carbon Dioxide and Water using TraPPE and SWM4-NDP Force Fields

The increased levels of carbon dioxide (CO 2 ) emissions due to the combustion of fossil fuels and the consequential impact on global climate change have made CO 2 capture, storage, and utilization a significant area of focus for current research. In most electrochemical CO 2 applications, water is used as a proton donor due to its high availability and mobility and use as a polar solvent. Additionally, supercritical CO 2 is a promising avenue for electrochemical applications due to its unique chemical and physical properties. Consequently, understanding the interactions between water and supercritical CO 2 is of great importance for future electrochemical applications. Molecular dynamics (MD) simulation is a powerful tool that enables atomistic-resolution dynamics of molecular systems, which can complement and guide future experimental investigations. This study employed atomistic MD to study the cosolubilities, codiffusivities, and structure of supercritical CO 2 and water systems, with a polarizable water model (SWM4-NDP) and a nonpolarizable CO 2 model (TraPPE). Additionally, ab initio MD simulations were used to better understand how atomistic polarizable/nonpolarizable models compare to explicit modeling of electron densities. The polarizable water model exhibited substantial improvement in water-associated properties. In conclusion, we anticipate the development of a compatible polarizable CO 2 model to yield similar improvement, providing a pathway for realizing novel high-pressure electrochemical systems.

25 ENERGY STORAGE↗
NASA NTRS2023

Molecular Dynamics Simulations of Austenite-Martensite Interfaces in NiTi Shape Memory Alloys

The unique properties of shape memory alloys (SMAs) arise from a reversible martensitic transformation. The nucleation and migration of austenite-martensite interfaces are the key to understanding the SMA properties. Molecular dynamics (MD) simulations can provide important atomic-scale information about these aspects of the transformation, but their time scales prevent the interface formation under near-equilibrium conditions relevant to experiment. We present a new MD methodology which allows for the natural formation of energetically preferred austenite-martensite interfaces under near-equilibrium conditions. Our simulation demonstrates that the interfaces in NiTi are semi-coherent, composed of a series of terrace planes and structural disconnections, and they migrate rapidly through single crystals with only a small thermodynamic driving force. In bi-crystals and polycrystals, the migration of these same interfaces is significantly impeded by grain boundaries and stored elastic energy. This behavior can result in SMA hysteresis via several mechanisms associated with nucleation and non-elastic strain accommodation.

Gabriel Plummer↗
DOE OSTI2026

Iron surface corrosion in supercritical CO2 at atomic scale investigated by molecular dynamics simulations

Understanding the corrosion behavior of steels in supercritical carbon dioxide (S-CO2) is essential for ensuring the safe application of S-CO2 as a heat-transfer fluid in high-temperature energy systems, including advanced nuclear reactors. In this work, molecular dynamics (MD) simulations using ReaxFF potential are performed to explore the atomic-scale corrosion mechanisms of body-centered cubic iron (BCC-Fe) in S-CO2. The results show that CO2 molecules in S-CO2 decompose at the Fe surface, generating free C and O atoms that form Fe-C and Fe-O bonds and subsequently produce oxides and carbides. Concurrently, Fe atoms dissolve from the surface and diffuse into the S-CO2 region, resulting in interdiffusion of Fe, C and O atoms at the interface. The corrosion-layer thickness calculations show that high pressure and temperature induced by S-CO2 have stronger effects than surface orientation on the corrosion process. In addition, surface Fe atoms undergo substantial displacement under S-CO2 exposure, further accelerating corrosion. When a radiation-induced void is introduced near the Fe surface, the corrosion is enhanced. The void-matrix interface expands the reaction surface area and simultaneously induces corrosion reactions inside the bulk, resulting in a deeper penetration of C and O and thicker corrosion layers. All these results indicate that high-temperature, high-pressure and radiation-induced voids can seriously affect the corrosion of Fe in S-CO2, and must be considered to better use S-CO2 in nuclear facilities.

Li, Wenhua↗
DOE OSTI2024

Shock compression of crystalline TeO 2 to the high-pressure fluid regime: Insights from ab initio molecular dynamics simulations

The shock response of fully-dense and porous crystalline tellurium dioxide (TeO 2 ⁠) to the high-pressure and high-temperature fluid regime was investigated within the framework of density functional theory with Mermin’s generalization to finite temperatures. The principal and porous shock Hugoniot curves were predicted from canonical ab initio molecular dynamics (AIMD) simulations, with the phase space sampled along isotherms up to 80 000 K, for densities ranging from ρ = 3 to 17 g/cm 3 . The polymorphs investigated are α-TeO 2 paratellurite (⁠P4 1 2 1 2), TeO 2 cotunnite (⁠Pnma⁠), and TeO 2 post-cotunnite (⁠P2 1 /m⁠). Based on the discontinuity found in the calculated U s – u p slope of TeO 2 post-cotunnite at a shock velocity of U s ≃ 8.35km/s and a particle velocity of u p ≃ 3.64 km/s, the shock melting temperature and pressure are predicted to be ≃ 6500 K and ≃ 170 GPa. Results from the AIMD simulations are in line with the static compression data of TeO 2 paratellurite and cotunnite, and with the recent shock Hugoniot data for single-crystal α- TeO 2 for pressures up to 85 GPa, obtained using the inclined-mirror method and the velocity interferometer system for any reflector combined with powder gun and two-stage light-gas gun.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Unraveling the Nature of Vibrational Dynamics in CsPbI 3 by Inelastic Neutron Scattering and Molecular Dynamics Simulations

Cesium lead iodide, CsPbI 3 , is an optoelectronic material of large interest for various technological applications; however, fundamental questions surrounding the vibrational dynamics of this material, especially regarding its role in structural phase transitions, remain to be elucidated. Here, in a combined variable temperature inelastic neutron scattering (INS) and machine-learning based molecular dynamics (MD) simulation study, we show that the stable phase at room temperature, i.e., the nonperovskite δ-phase, exhibits phonon modes with weak anharmonicity with only a weak temperature dependence from 10 K all the way up to the transition to the cubic perovskite α-phase at approximately 600 K. In contrast, the α-phase features anharmonic and damped vibrational dynamics, mainly associated with overdamped tilting motions of the PbI 6 octahedra. Crucially, these overdamped tilting modes, which relate to the tetragonal and orthorhombic distorted perovskite phases (β- and γ-phase, respectively) formed at lower temperatures, stay overdamped by more than 100 K above the respective phase transition. This suggests a flat energy landscape of octahedral tilting motions in α-CsPbI3 and with structural fluctuations on the picosecond time scale with tilting patterns that locally resemble the structure of the β- and γ-phases. The vibrational dynamics of α-CsPbI 3 are also characterized by pronounced anharmonic motions with large thermal displacements of the Cs + ions, but these modes remain underdamped at 600 K.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗