Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “initio molecular dynamics”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Tuning collective anion motion enables superionic conductivity in solid-state halide electrolytes

Halides of the family Li 3 MX 6 (M = Y, In, Sc and so on, X = halogen) are emerging solid electrolyte materials for all-solid-state Li-ion batteries. They show greater chemical stability and wider electrochemical stability windows than existing sulfide solid electrolytes, but have lower room-temperature ionic conductivities. Here, in this work, we report the discovery that the superionic transition in Li 3 YCl 6 is triggered by the collective motion of anions, as evidenced by synchrotron X-ray and neutron scattering characterizations and ab initio molecular dynamics simulations. Based on this finding, we used a rational design strategy to lower the transition temperature and thus improve the room-temperature ionic conductivity of this family of compounds. We accordingly synthesized Li 3 YCl x Br 6-x and Li 3 GdCl 3 Br 3 and achieved very high room-temperature conductivities of 6.1 and 11 mS cm -1 for Li 3 YCl 4.5 Br 1.5 and Li 3 GdCl 3 Br 3 , respectively. These findings open new routes to the design of room-temperature superionic conductors for high-performance solid batteries.

25 ENERGY STORAGE↗
DOE OSTI2021

A two-dimensional type I superionic conductor

Superionic conductors (SICs) possess liquid-like ionic diffusivity in the solid state, finding wide applicability from electrolytes in energy storage to materials for thermoelectric energy conversion. Type I SICs (e.g., AgI, Ag 2 Se, etc.) are defined by a first-order transition to the superionic state and have so far been found exclusively in three-dimensional crystal structures. Here, we reveal a two-dimensional type I SIC, α-KAg 3 Se 2 by scattering techniques and complementary simulations. Quasi-elastic neutron scattering and ab initio molecular dynamics simulations confirm that the superionic Ag + ions are confined to sub-nanometre sheets, with the simulated local structure validated by experimental X-ray powder pair-distribution-function analysis. Finally, we demonstrate that the phase transition temperature can be controlled by chemical substitution of the alkali metal ions that comprise the immobile charge-balancing layers. Our work thus extends the known classes of SICs and will facilitate the design of new materials with tailored ionic conductivities and phase transitions.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Tuning aminopolycarboxylate chelators for efficient complexation of trivalent actinides

The complexation of trivalent lanthanides and minor actinides (Am 3+ , Cm 3+ , and Cf 3+ ) by the acyclic aminopolycarboxylate chelators 6,6'-((ethane-1,2-diylbis–((carboxymethyl)azanediyl))bis–(methylene))dipicolinic acid (H 4 octapa) and 6,6'-((((4-(1-(2-(2-(2-hydroxyethoxy)ethoxy)ethyl)-1H-1,2,3-triazol-4-yl)pyridine-2,6-diyl)bis–(methylene))bis–((carboxymethyl)azanediyl))bis–(methylene)) dipicolinic acid (H 4 pypa-peg) were studied using potentiometry, spectroscopy, competitive complexation liquid–liquid extraction, and ab initio molecular dynamics simulations. Two studied reagents are strong multidentate chelators, well-suited for applications seeking radiometal coordination for in-vivo delivery and f-element isolation. The previously reported H 4 octapa forms a compact coordination packet, while H 4 pypa-peg is less sterically constrained due to the presence of central pyridine ring. The solubility of H 4 octapa is limited in a non-complexing high ionic strength perchlorate media. However, the introduction of a polyethylene glycol group in H 4 pypa-peg increased the solubility without influencing its ability to complex the lanthanides and minor actinides in solution.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Search for a Grotthuss mechanism through the observation of proton transfer

The transport of protons is critical in a variety of bio- and electro-chemical processes and technologies. The Grotthuss mechanism is considered to be the most efficient proton transport mechanism, generally implying a transfer of protons between ‘chains’ of host molecules via elementary reactions within the hydrogen bonds. Although Grotthuss proposed this concept more than 200 years ago, only indirect experimental evidence of the mechanism has been observed. Here we report the first experimental observation of proton transfer between the molecules in pure and 85% aqueous phosphoric acid. Employing dielectric spectroscopy, quasielastic neutron, and light scattering, and ab initio molecular dynamic simulations we determined that protons move by surprisingly short jumps of only ~0.5–0.7 Å, much smaller than the typical ion jump length in ionic liquids. Our analysis confirms the existence of correlations in these proton jumps. However, these correlations actually reduce the conductivity, in contrast to a desirable enhancement, as is usually assumed by a Grotthuss mechanism. Furthermore, our analysis suggests that the expected Grotthuss-like enhancement of conductivity cannot be realized in bulk liquids where ionic correlations always decrease conductivity.

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

A hybrid calorimetry-simulation model of mixing enthalpy for molten salt

Calorimetric determination of enthalpies of mixing (ΔH mix ) in multicomponent molten salts is often interpreted using empirical models that lack physically meaningful parameters. However, for improving pyrochemical separation of spent nuclear fuel, where lanthanides are major fission products and critical elements, a deeper thermodynamic understanding of the link between excess thermodynamic properties and solvation structure is critically needed. In this work, we implement a hybrid and physics-informed framework, MIVM+Calorimetry+AIMD, which integrates experimentally measured ΔH mix (via high temperature drop calorimetry) with solvation structures from ab initio molecular dynamics (AIMD). This approach is demonstrated using LaCl 3 mixed with eutectic LiCl-KCl (58 mol% – 42 mol%) at 873 K and 1133 K. MIVM-derived parameters enable extrapolation of excess Gibbs energy and La 3+ activity across compositions. In contrast, direct ΔH mix predictions from AIMD and polarizable ion model simulations deviate significantly. By incorporating experimentally benchmarked solvation structures into an interpretable thermodynamic model, the MIVM+Calorimetry+AIMD formalism achieves higher accuracy and generalizable method for studying molten salts, offering a robust path for understanding and optimizing molten salt chemistry relevant to nuclear fuel cycles and separation science.

Goncharov, Vitaliy G. [Washington State Univ., Pul↗
DOE OSTI2025

Pausing ultrafast melting by timed multiple femtosecond-laser pulses

An intense femtosecond-laser excitation of a solid induces highly nonthermal conditions. In materials like silicon, laser-induced bond-softening leads to a highly incoherent ionic motion and eventually nonthermal melting. But is this outcome an inevitable consequence, or can it be controlled? Here, we performed ab initio molecular dynamics simulations of crystalline silicon after timed multiple femtosecond-laser pulse excitations with fluence above the nonthermal melting threshold. Our results demonstrate an excitation mechanism that pauses nonthermal melting and creates a metastable state instead, with an electronic structure similar to the ground state. This mechanism can be generalized to other materials, potentially enabling structural and/or electronic transitions to metastable phases in the high-excitation regime. In addition, our approach could be used to switch off nonthermal contributions in experiments, allowing reliable electron-phonon coupling constants to be obtained more easily.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

CHGNet as a pretrained universal neural network potential for charge-informed atomistic modelling

Large-scale simulations with complex electron interactions remain one of the greatest challenges for atomistic modelling. Although classical force fields often fail to describe the coupling between electronic states and ionic rearrangements, the more accurate ab initio molecular dynamics suffers from computational complexity that prevents long-time and large-scale simulations, which are essential to study technologically relevant phenomena. Here we present the Crystal Hamiltonian Graph Neural Network (CHGNet), a graph neural network-based machine-learning interatomic potential (MLIP) that models the universal potential energy surface. CHGNet is pretrained on the energies, forces, stresses and magnetic moments from the Materials Project Trajectory Dataset, which consists of over 10 years of density functional theory calculations of more than 1.5 million inorganic structures. The explicit inclusion of magnetic moments enables CHGNet to learn and accurately represent the orbital occupancy of electrons, enhancing its capability to describe both atomic and electronic degrees of freedom. We demonstrate several applications of CHGNet in solid-state materials, including charge-informed molecular dynamics in Li x MnO 2 , the finite temperature phase diagram for Li x FePO 4 and Li diffusion in garnet conductors. We highlight the significance of charge information for capturing appropriate chemistry and provide insights into ionic systems with additional electronic degrees of freedom that cannot be observed by previous MLIPs.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Hydration structure and water exchange kinetics at xenotime–water interfaces: implications for rare earth minerals separation

Hydration of surface ions gives rise to structural heterogeneity and variable exchange kinetics of water at complex mineral–water interfaces. In this work, we employ ab initio molecular dynamics (AIMD) simulations and water adsorption calorimetry to examine the aqueous interfaces of xenotime, a phosphate mineral that contains predominantly Y 3+ and heavy rare earth elements. Consistent with natural crystal morphology, xenotime is predicted to have a tetragonal prismatic shape, dominated by the {100} surface. Hydration of this surface induces multilayer interfacial water structures with distinct OH orientations, which agrees with recent crystal truncation rod measurements. The exchange kinetics between two adjacent water layers exhibits a wide range of underlying timescales (5–180 picoseconds), dictated by ion–water electrostatics. Adsorption of a bidentate hydroxamate ligand reveals that {100} xenotime surface can only accommodate monodentate coordination with water exchange kinetics strongly depending on specific ligand orientation, prompting us to reconsider traditional strategies for selective separation of rare-earth minerals.

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

Molecular growth upon ionization of van der Waals clusters containing HCCH and HCN is a pathway to prebiotic molecules

The growth mechanisms of organic molecules in an ionizing environment such as the interstellar medium are not completely understood. Here we examine by means of ab initio molecular dynamics (AIMD) simulations and density functional theory (DFT) computations the possibility of bond formation and molecular growth upon ionization of van der Waals clusters of pure HCN clusters, and mixed clusters of HCN and HCCH, both of which are widespread in the interstellar medium. Ionization of van der Waals clusters can potentially lead to growth in low temperature and low-density environments. Our results show, that upon ionization of the pure HCN clusters, strongly bound stable structures are formed that contain N∙∙∙H bonds, and growth beyond pairwise HCN molecules is seen only in a small percentage of cases. In contrast, mixed clusters, where HCCH is preferentially ionized over HCN, can grow up to 3 or 4 units long with new carbon–carbon and carbon–nitrogen covalent bonds. Moreover, cyclic molecules formed, such as the radical cation of pyridine, which is a prebiotic molecule. The results presented here are significant as they provide a feasible pathway for molecular growth of small organic molecules containing both carbon and nitrogen in cold and relatively denser environments such as in dense molecular clouds but closer to the photo-dissociation regions, and protoplanetary disks. In the mechanism we propose, first, a neutral van der Waals cluster is formed. Once the cluster is formed it can undergo photoionization which leads to chemical reactivity without any reaction barrier.

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

Theoretical evidence of water serving as a promoter for lithium superoxide disproportionation in Li-O 2 batteries

Experimental evidence has demonstrated that the presence of water in non-aqueous electrolytes significantly affects Li-O 2 electrochemistry. Understanding the reaction mechanism for Li 2 O 2 formation in the presence of water impurities is important to understand Li-O 2 battery performance. A recent experiment has found that very small amounts of water (as low as 40 ppm) can significantly affect the product formation in Li-O 2 batteries as opposed to essentially no water (1 ppm). Although experimental as well as theoretical work has proposed mechanisms of Li 2 O 2 formation in the presence of much larger amounts of water, none of the mechanisms provide an explanation for the observations for very small amounts of water. In this work, density functional theory (DFT) was utilized to obtain a mechanistic understanding of the Li-O 2 discharge chemistry in a dimethoxyethane (DME) electrolyte containing an isolated water and no water. The reaction pathways for Li 2 O 2 formation from LiO 2 on a model system were carefully evaluated with different level of theories, i.e. PBE (PW), B3LYP/6-31G(2df,p), B3LYP/6-311++G(2df,p) and G4MP2. The results indicate that the LiO 2 disproportionation reaction to Li 2 O 2 can be promoted by the water in DME electrolyte, which explains why there is a significant difference compared to when no water is present in the experimentally observed discharge product distributions. Ab initio molecular dynamics calculations were also used to investigate the disproportionation of LiO 2 dimer in explicit DME. This work adds to the fundamental understanding of the discharge chemistry of a Li-O 2 battery.

25 ENERGY STORAGE↗
DOE OSTI2020

Thermodynamic and kinetic properties of layered-CaCo2O4 for the Ca-ion batteries: a systematic first-principles study

One of the more promising directions in multivalent energy storage is systems based on Ca ion intercalation due to the potential for high voltage and capacity. A major challenge for enabling such a battery is to find cathode materials capable of fast ionic diffusion and reversible insertion of Ca ions. Here, on the basis of first-principles calculations, we have demonstrated that layered CaCo 2 O 4 exhibits favorable thermodynamic and kinetic properties that should enable topotactic Ca ion intercalation reactions. The P3-type layered Ca x Co 2 O 4 (0 < x < 1) with either of space groups of P 1 or P 2 1 / m are stable at multiple Ca concentrations and show a smooth voltage plateau higher than 3 V up to X = 0.5. The energy barriers of the single Ca ion migration are as low as 0.36 eV and 0.27 eV at the dilute and high vacancy concentration limits, respectively. Therefore, although varying the vacancy environments of the diffusing atom influences the migration barriers, they do not exceed 0.6 eV. Stochastic analysis of Ca hopping events performed by ab initio molecular dynamics (AIMD) simulation has shown that the migration barriers are lower than 0.32 eV. Therefore, the Ca diffusivity at room temperature extrapolated from the AIMD results is comparable to Li diffusivity (>10 -10 cm 2 s -1 ) in conventional Li cathode materials, suggesting the feasibility of layered Ca x Co 2 O 4 as multivalent cathode materials. Finally, the structural factors that enable fast diffusion are discussed.

25 ENERGY STORAGE↗
DOE OSTI2021

Efficient diffusion of superdense lithium via atomic channels for dendrite-free lithium–metal batteries

The non-uniform aggregation of fast-diffused Li on an anode surface would aggravate its tip-effect-induced nucleation/growth, leading to the notorious dendrite growth in Li metal batteries (LMBs). Tuning the Li diffusion on the anode surface has been regarded previously as a mainstream method to induce its uniform deposition, while the diffusion of Li in the anode bulk is usually ignored. In this study, conceptually different from the classic surface modification, we propose a molecular tunnelling strategy to construct atomic channels in graphite bulk, which enables the fast diffusion of superdense Li. Density functional theory calculations and ab initio molecular dynamics simulations prove that the bulk diffusion through atomic channels could become a new and dominating path. Its reversible and efficient diffusion has been further visualized by in situ transmission electron microscopy. As a result, when coupled with high-loading LiFePO 4 cathodes (20 mg cm –2 ), a high areal capacity and 100% capacity retention are achieved over 370 cycles. Through this work a new strategy is developed based on the bulk-diffusion of superdense Li for dendrite-free LMBs, which can be pervasive in other high-performance energy storage systems.

25 ENERGY STORAGE↗
DOE OSTI2021

A solid electrolyte interphase to protect the sulfurized polyacrylonitrile (SPAN) composite for Li–S batteries: computational approach addressing the electrolyte/SPAN interfacial reactivity

This study addresses the reactivity of multiple solvents and lithium bis(fluorosulfonyl)imide (LiFSI) at the interface with sulfurized polyacrylonitrile (SPAN) in multiple stages of lithiation via ab initio molecular dynamics simulations. Here, both ether 1,3-dioxolane (DOL) and dimethyl carbonate (DMC) proved stable on the lithiated SPAN surface regardless of the lithium content, meaning that neither of these species likely contributes to growing a solid electrolyte interphase (SEI) coating to protect the SPAN composite. Conversely, cyclic carbonates, ethylene carbonate (EC) and fluoroethylene carbonate (FEC) are shown to be very active. The EC reduction occurs only on a highly lithiated surface with a 3.0 Li/S molar ratio, the equivalent of the SPAN composite in an over-discharge regime with voltages close to 0.0 V vs. Li/Li + . The FEC reduction starts with a 2.0 Li/S molar ratio and above, suggesting that FEC could act as a useful additive in electrolyte formulations with EC. Both EC and FEC follow multiple reduction mechanisms to produce complex reduction products and LiF in the FEC case. We provide a mechanistic description for each detected decomposition path. The LiFSI salt also proves reactive against the lithiated SPAN surface. The FSI – defluorination is the dominant reduction path. However, the SO 2 NSO 2 F – and SO 2 NSO 2 2– species proved stable against S–N cleavage. This behavior makes the LiFSI salt a potential candidate for SPAN-based Li–S batteries because it produces LiF without releasing SO 2 .

25 ENERGY STORAGE↗
DOE OSTI2021

Materials design of sodium chloride solid electrolytes Na3MCl6 for all-solid-state sodium-ion batteries

All-solid-state sodium-ion batteries have attracted increasing attention owing to the low cost of sodium and the enhanced safety compared to conventional Li-ion batteries. Recently, halides have been considered as promising solid electrolytes (SEs) due to their favorable combination of high ionic conductivity and chemical stability against high-voltage cathode materials. Although a wide variety of lithium chloride SEs, Li 3 MCl 6 , have been developed for high-voltage all-solid-state batteries, only a limited number of sodium chloride SEs have been reported. This study aims to offer a material design insight for the development of sodium chloride SEs through systematic assessment of the phase stability, electrochemical stability, and transport properties of novel Na 3 MCl 6 SEs. Structural calculations indicate that Na 3 MCl 6 exhibits trigonal $P\bar{3}$1c, monoclinic P2 1 /n, and trigonal $R\bar{3}$ phases, and the stable phase of Na 3 MCl 6 is dependent on the type and ionic radius of M. Na 3 MCl 6 typically exhibits a high oxidation potential, demonstrating good electrochemical stability against cathodes. The bond-valence site energy and ab initio molecular dynamics calculations revealed that Na 3 MCl 6 with P2 1 /n and $R\bar{3}$ phases showed low ionic conductivity, while the $P\bar{3}$1c phase slightly improved the ionic conductivity of Na 3 MCl 6 . The formation of Na vacancies by aliovalent substitution considerably increased the ionic conductivity up to four orders of magnitude for pristine Na 3 MCl 6 , exhibiting ~10 –5 S cm –1 for trigonal $P\bar{3}$1c and $R\bar{3}$ phases. The formation of defects could further enhance the ionic conductivity of Na 3 MCl 6 , and the optimization of defect type and ratio can be helpful in developing superionic Na chloride SEs. The material design of Na 3 MCl 6 in this study will provide fundamental guidelines for the development of novel sodium halide SEs for all-solid-state sodium-ion batteries.

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

Hydrogen in disordered titania: connecting local chemistry, structure, and stoichiometry through accelerated exploration

Hydrogen incorporation in native surface oxides of metal alloys often controls the onset of metal hydriding, with implications for materials corrosion and hydrogen storage. A key representative example is titania, which forms as a passivating layer on a variety of titanium alloys for structural and functional applications. These oxides tend to be structurally diverse, featuring polymorphic phases, grain boundaries, and amorphous regions that generate a disparate set of unique local environments for hydrogen. Here, we introduce a workflow that can efficiently and accurately navigate this complexity. First, a machine learning force field, trained on ab initio molecular dynamics simulations, was used to generate amorphous configurations. Density functional theory calculations were then performed on these structures to identify local oxygen environments, which were compared against experimental observations. Second, to classify subtle differences across the disordered configuration space, we employ a graph-based sampling procedure. Finally, local hydrogen binding energies and hopping kinetics are computed using exhaustive density functional theory calculations on representative configurations. Here, we leverage this methodology to show that hydrogen binding energetics are described by local oxygen coordination, which in turn is affected by stoichiometry, and form the basis of hopping kinetics and diffusion. Together these results imply that hydrogen incorporation and transport in TiO x can be tailored through compositional engineering, with implications for improving the performance and durability of titanium-derived alloys in hydrogen environments.

36 MATERIALS SCIENCE↗