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

Development of a High-Rate Lithium-Air Battery Using a Gaseous CO 2 Reactant

Li-air batteries are considered a potential alternative to Li-ion batteries for transportation applications due to their high theoretical specific energy. Most works in this area focus on use of O 2 as the reactant. However, newer concepts for using gaseous reactants (such as CO 2 , which has a theoretical specific energy density of 1,876 Wh/kg) provide opportunities for further exploration. The main objective of this project was the development of a novel strategy that enables operation of Li-CO 2 batteries at high-capacity and high-rate, with a long-cycle-life. The team was able to: (1) Synthesize two novel transition metal chalcogenide (TMC) catalysts that work in synergy with ionic liquid-based electrolytes to enhance the efficiency of reactions during discharge and charge processes; (2) Fabricate high-porosity cathode electrodes with 3D printing to increase electrode surface area and gas permeability; (3) Develop a multiscale modeling framework that integrates Density Functional Theory (DFT), Ab-Initio Molecular Dynamics (AIMD), classical Molecular Dynamics (MD), and Finite Element Analysis (FEA) to investigate atomic and cell-level properties of Li-CO 2 batteries; (4) Assemble a stackable Li-CO 2 pouch-cell able to deliver a capacity of >200 mAh. These achievements were realized through an integrated approach based on materials synthesis, testing, characterization, analysis, and computation. This project produced a thorough understanding of key chemical, electronic, and kinetic parameters that govern the operation of Li- CO 2 batteries in realistic conditions. The methodologies employed, and the insight generated, will be valuable beyond advancing the field of Li-CO 2 batteries

25 ENERGY STORAGE↗
DOE OSTI2026

Hydrated Metal and Metal-Nitrate Complexes in Water: Full Lanthanide(III) Series plus Miscellaneous Metal Ions

This is a dataset of hydrated metal complexes and metal–nitrate hydrated complexes intended for public use, reproducibility, and downstream structural analysis. A key feature is coverage across the full lanthanide(III) series (La–Lu), enabling systematic comparisons of coordination motifs and bonding trends across the entire lanthanide sequence. In addition to the lanthanides, the dataset also includes other metal ions such as UO2(VI), Fe(II), and Fe(III). The dataset provides optimized geometries for hydrated and nitrate-containing hydrated complexes, together with representative ab initio molecular dynamics (AIMD) trajectories saved in standard XYZ formats. The accompanying NWChem input decks enable reproduction of the reported calculations and provide a starting point for extending the simulations to related coordination environments. Computationally, DFT calculations employ the B3LYP functional with DFT-D3BJ dispersion corrections and a COSMO continuum solvent model (dielectric constant 78.4) to represent solvation beyond the explicitly treated first hydration shell. AIMD simulations are performed with the NWChem qmd module at 298 K, integrating nuclear motion with the velocity-Verlet algorithm and controlling temperature using a Nosé–Hoover thermostat. Trajectories are approximately 4.8 ps in length and are used primarily to assess short-time stability of candidate coordination motifs, including (for lanthanides) differences between 8- versus 9-water coordination and comparisons between nitrate-bound and nitrate-free hydrated complexes.

Dinpajooh, Mohammadhasan [Pacific Northwest Nation↗
DOE OSTI2024

Functionalization of Cathode–Electrolyte Interface with Ionic Liquids for High-Performance Quasi-Solid-State Lithium–Sulfur Batteries: A Low-Sulfur Loading Study

We introduce a quasi-solid-state electrolyte lithium-sulfur (Li–S) battery (QSSEB) based on a novel Li-argyrodite solid-state electrolyte (SSE), Super P–Sulfur cathode, and Li-anode. The cathode was prepared using a water-based carboxymethyl cellulose (CMC) solution and styrene butadiene rubber (SBR) as the binder while Li6PS5F0.5Cl0.5 SSE was synthesized using a solvent-based process, via the introduction of LiF into the argyrodite crystal structure, which enhances both the ionic conductivity and interface-stabilizing properties of the SSE. Ionic liquids (IL) were prepared using lithium bis(trifluoromethyl sulfonyl)imide (LiTFSI) as the salt, with pre-mixed pyrrolidinium bis(trifluoromethyl sulfonyl)imide (PYR) as solvent and 1,3-dioxolane (DOL) as diluent, and they were used to wet the SSE–electrode interfaces. The effect of IL dilution, the co-solvent amount, the LiTFSI concentration, the C rate at which the batteries are tested and the effect of the introduction of SSE in the cathode, were systematically studied and optimized to develop a QSSEB with higher capacity retention and cyclability. Interfacial reactions occurring at the cathode–SSE interface during cycling were also investigated using electrochemical impedance spectroscopy, cyclic voltammetry, and X-ray photoelectron spectroscopy supported by ab initio molecular dynamics simulations. This work offers a new insight into the intimate interfacial contacts between the SSE and carbon–sulfur cathodes, which are critical for improving the electrochemical performance of quasi-solid-state lithium–sulfur batteries.

25 ENERGY STORAGE↗
DOE OSTI2021

Free Energies of Hydrated Halide Anions: High Through-Put Computations on Clusters to Treat Rough Energy-Landscapes

With a longer-term goal of addressing the comparative behavior of the aqueous halides F−, Cl−, Br−, and I− on the basis of quasi-chemical theory (QCT), here we study structures and free energies of hydration clusters for those anions. We confirm that energetically optimal (H2O)nX clusters, with X = Cl−, Br−, and I−, exhibit surface hydration structures. Computed free energies, based on optimized surface hydration structures utilizing a harmonic approximation, typically (but not always) disagree with experimental free energies. To remedy the harmonic approximation, we utilize single-point electronic structure calculations on cluster geometries sampled from an AIMD (ab initio molecular dynamics) simulation stream. This rough-landscape procedure is broadly satisfactory and suggests unfavorable ligand crowding as the physical effect addressed. Nevertheless, this procedure can break down when n≳4, with the characteristic discrepancy resulting from a relaxed definition of clustering in the identification of (H2O)nX clusters, including ramified structures natural in physical cluster theories. With ramified structures, the central equation for the present rough-landscape approach can acquire some inconsistency. Extension of these physical cluster theories in the direction of QCT should remedy that issue, and should be the next step in this research direction.

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

Mechanochemical Synthesis of Sustainable Ternary and Quaternary Nanostructured Cu2SnS3, Cu2ZnSnS4, and Cu2ZnSnSe4 Chalcogenides for Thermoelectric Applications

Copper-based chalcogenides have emerged as promising thermoelectric materials due to their high thermoelectric performance, tunable transport properties, earth abundance and low toxicity. We have presented an overview of experimental results and first-principal calculations investigating the thermoelectric properties of various polymorphs of Cu2SnS3 (CTS), Cu2ZnSnS4 (CZTS), and Cu2ZnSnSe4 (CZTSe) synthesized by high-energy reactive mechanical alloying (ball milling). Of particular interest are the disordered polymorphs of these materials, which exhibit phonon-glass–electron-crystal behavior—a decoupling of electron and phonon transport properties. The interplay of cationic disorder and nanostructuring leads to ultra-low thermal conductivities while enhancing electronic transport. These beneficial transport properties are the consequence of a plethora of features, including trap states, anharmonicity, rattling, and conductive surface states, both topologically trivial and non-trivial. Based on experimental results and computational methods, this report aims to elucidate the details of the electronic and lattice transport properties, thereby confirming that the higher thermoelectric (TE) performance of disordered polymorphs is essentially due to their complex crystallographic structures. In addition, we have presented synchrotron X-ray diffraction (SR-XRD) measurements and ab initio molecular dynamics (AIMD) simulations of the root-mean-square displacement (RMSD) in these materials, confirming anharmonicity and bond inhomogeneity for disordered polymorphs.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Violation of the T –1 Relationship in the Lattice Thermal Conductivity of Mg 3 Sb 2 with Locally Asymmetric Vibrations

Most crystalline materials follow the guidelines of T –1 temperature-dependent lattice thermal conductivity (κ L ) at elevated temperatures. Here, we observe a weak temperature dependence of κ L in Mg 3 Sb 2 , T –0.48 from theory and T –0.57 from measurements, based on a comprehensive study combining ab initio molecular dynamics calculations and experimental measurements on single crystal Mg 3 Sb 2 . These results can be understood in terms of the so-called “phonon renormalization” effects due to the strong temperature dependence of the interatomic force constants (IFCs). The increasing temperature leads to the frequency upshifting for those low-frequency phonons dominating heat transport, and more importantly, the phonon-phonon interactions are weakened. In-depth analysis reveals that the phenomenon is closely related to the temperature-induced asymmetric movements of Mg atoms within MgSb 4 tetrahedron. With increasing temperature, these Mg atoms tend to locate at the areas with relatively low force in the force profile, leading to reduced effective 3 rd -order IFCs. The locally asymmetrical atomic movements at elevated temperatures can be further treated as an indicator of temperature-induced variations of IFCs and thus relatively strong phonon renormalization. The present work sheds light on the fundamental origins of anomalous temperature dependence of κ L in thermoelectrics.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Suppressed paramagnetism in amorphous Ta 2 O 5 − x oxides and its link to superconducting-qubit performance

Amorphous-oxide layers in thin-film capacitors are linked to reduced transmon-qubit T 1 coherence times. Ta -based capacitors outperform Nb -based ones, suggesting that amorphous Ta 2 O 5 − x is less lossy than Nb 2 O 5 − x . We investigate the microscopic features of these amorphous oxides using ab initio molecular dynamics and density functional theory, revealing the origins of the superior performance of Ta 2 O 5 − x . We establish that oxygen deficiency is less likely to occur in amorphous Ta 2 O 5 − x than in Nb 2 O 5 − x for 0 ≤ x ≤ 0.25 and that for a given oxygen deficiency x , metal Ta — Ta bond formation is enhanced. Such bonds, which are accommodated by structural flaws in the amorphous network, capture electrons better than in amorphous Nb 2 O 5 − x . These thermochemical differences quench or highly suppress magnetic moments in amorphous Ta 2 O 5 − x and eliminate a potential source of quasiparticles and magnetic flux noise. We also show that hyperfine couplings between Nb nuclei and local magnetic moments in Nb 2 O 5 − x can form “two-level systems” (TLSs) or “two-level fluctuators” with energy splittings of 100–1000 MHz or higher. This reveals a TLS mechanism in amorphous Nb 2 O 5 − x oxide layers that is likely inactive in Ta 2 O 5 − x . Our work provides a fundamental understanding of the materials chemistry and limitations imposed by native oxides of superconducting qubits that can be used to guide materials selection and processing.

Pritchard, P. Graham [Northwestern U.] (ORCID:0000↗
DOE OSTI2022

Synthesis and Thermophysical Property Determination of NaCl-PuCl3 Salts

Currently, a knowledge gap exists in the available data and understanding of thermophysical properties relating to fresh fuel salts, especially those containing plutonium. Determination of this data is necessary for the design and construction of test reactors, as well as the licensing of future commercial molten salt reactors. Thermophysical properties such as melting temperature, salt stability, density, and heat capacity were determined on synthesized eutectic NaCl-PuCl3 and a more sodium rich composition containing 25 mol% PuCl3. These measurements document the baseline properties of the salt as a function of temperature for future experiments on irradiated fuel salt which will provide a holistic perspective on the change of thermophysical properties during reactor operations. It was determined that the NaCl-PuCl3 ingot synthesized for this study contained 63.4 mol% NaCl, 36.3 mol% PuCl3, and was 99.7% pure. Upon heating it was shown that the NaCl- PuCl3 eutectic was stable at temperatures of 800°C. The onset of melting occurred at 541°C and the enthalpy of fusion was determined to be 140.7 ± 8.4 J/g. Specific heat capacity measurements showed a slightly decreasing trend with respect to temperature in the liquid phase ranging from 0.67 to 0.57 J/g·K, with an average value of 0.637 ± 0.03 J/g·K (104 ± 5 J/mol·K) between 500 to 720°C. Three independent trials of the molten NaCl-PuCl3 eutectic salt found the density to be ?(T) = 3.8589 – 9.5342·10-4 T(°C), validated between 500 to 800°C. In addition to salt synthesis and experimentally determining thermodynamic properties, Ab Initio molecular dynamic (AIMD) simulations were used to calculate density and heat capacity values.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Demonstrate Improved Ag Diffusion and Describe the Basis for Pd Penetration Modeling in SiC

In past work, an effective diffusivity coefficient was determined for Ag transport through the silicon carbide layer of a tristructural isotropic fuel particle. The effective diffusivity coefficient accounts for the microstructure of the silicon carbide and includes both bulk diffusion and grain boundary diffusion of Ag. In this report, the model has been improved by accounting for the enhanced concentration of vacancies in the bulk due to irradiation, which substantially influence bulk diffusivity at low temperatures. To improve the BISON model and make it fission rate dependent, effective diffusivity calculations have been performed that incorporate the radiation modified bulk diffusivity. The microstructure and irradiation-dependent effective diffusivity has also been implemented into BISON, and its predictions for Ag release from tristructural isotropic fuel have been successfully compared to AGR-1 post irradiation measurements. Moreover, a new feature has been developed in the Multiphysics Object-Oriented Simulation Environment (MOOSE) to account for different grain boundary types. The Ag diffusivity in 5 (210)/[001] grain boundaries has been computed and was found to be greater than in random high-angle grain boundaries. The presence of the fission product Pd can also have an important effect on the properties of the silicon carbide layer in tristructural isotropic particles. The penetration of Pd into the silicon carbide layer causes a corrosion reaction that can lead to the failure of the silicon carbide layer; however, this corrosion reaction is not well understood. To enable an improved understanding of the mechanism, ab-initio molecular dynamics simulations of Pd interaction with bulk silicon carbide have been performed. The improved understanding of the reaction will form a basis for future improvements to the BISON’s Pd penetration failure model

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Predicting Hugoniot equation of state in erythritol with ab initio and reactive molecular dynamics

Erythritol has been proposed as an inert surrogate for developing theoretical and computational models to study aging in energetic materials. In this work, we present a comparison of mechanical and shock properties of erythritol computed using the ReaxFF reactive force field and from ab initio calculations employing density functional theory (DFT). We screened eight different ReaxFF parameterizations, of which the CHO parameters developed for hydrocarbon oxidation provide the most accurate predictions of mechanical properties and the crystal structure of erythritol. Further validation of the applicability of this ReaxFF parameterization for modeling erythritol is demonstrated by comparing predictions of the elastic constants, crystal structure, vibrational density of states, and Hugoniot curves against DFT calculations. The ReaxFF predictions are in close agreement with the DFT simulations for the elastic constants and shock Hugoniot when the crystal is loaded along its c axis but show as much as 30% disagreement in the elastic constants in the ab plane and 12% difference in shock pressures when shocked along the a or b crystal axes. Last, we compare thermomechanical properties predicted from classical molecular dynamics with those calculated using the quasi-harmonic approximation and show that quantum mechanical effects produce large discrepancies in the computed values of heat capacity and thermal expansion coefficients compared with classical assumptions. Combining classical molecular dynamics predictions of mechanical behavior with phonon-based calculations of thermal behaviors, we show that predicted shock-induced temperatures for pressures up to 6.5 GPa do not exceed the pressure-dependent melting point of erythritol.

Hu, Jing↗
DOE OSTI2023

Probing the self-ionization of liquid water with ab initio deep potential molecular dynamics

The chemical equilibrium between self-ionized and molecular water dictates the acid–base chemistry in aqueous solutions, yet understanding the microscopic mechanisms of water self-ionization remains experimentally and computationally challenging. Herein, Density Functional Theory (DFT)–based deep neural network (DNN) potentials are combined with enhanced sampling techniques and a global acid–base collective variable to perform extensive atomistic simulations of water self-ionization for model systems of increasing size. The explicit inclusion of long-range electrostatic interactions in the DNN potential is found to be crucial to accurately reproduce the DFT free energy profile of solvated water ion pairs in small (64 and 128 H 2 O) cells. The reversible work to separate the hydroxide and hydronium to a distance S is found to converge for simulation cells containing more than 500 H 2 O, and a distance of ∼ 8 Å is the threshold beyond which the work to further separate the two ions becomes approximately zero. The slow convergence of the potential of mean force with system size is related to a restructuring of water and an increase of the local order around the water ions. Calculation of the dissociation equilibrium constant illustrates the key role of long-range electrostatics and entropic effects in the water autoionization process.

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

Predicting thermodynamic and thermophysical properties of molten chloride salts from ab-initio and classical molecular dynamics simulations

Molten salt reactors (MSRs) are among the advanced concepts pursued under the generation IV nuclear energy technology umbrella. However, the basic concept is not new and was first developed as part of the effort to power aircrafts with nuclear energy in the 1950’s. Later in the 1960’s, Oak Ridge National Laboratory (ORNL) built and operated the Molten-Salt Reactor Experiment (MSRE). This reactor used a fluoride salt with uranium as fuel. Fluorides salts are still highly relevant and proposed in several designs. In addition, chloride salts are being considered for MSRs operating in the fast neutron spectrum. This report focuses on chloride salts.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

PySpawn: Software for Nonadiabatic Quantum Molecular Dynamics

The ab initio multiple spawning (AIMS) method enables nonadiabatic quantum molecular dynamics simulations in an arbitrary number of dimensions, with potential energy surfaces provided by electronic structure calculations performed on-the-fly. However, the intricacy of the AIMS algorithm complicates software development, deployment on modern shared computer resources, and post-simulation data analysis. PySpawn is a nonadiabatic molecular dynamics software package that addresses these issues. Here, the program is designed to be easily interfaced with electronic structure software, and an interface to the TeraChem software package is described here. PySpawn introduces a task-based reorganization of the AIMS algorithm, allowing fine-grained restart capability and setting the stage for efficient parallelization in a future release. PySpawn includes a user-friendly and interactive Python analysis module that will enable novice users to painlessly adopt AIMS. As a demonstration of PySpawn’s simulation capability and analysis module, we report complete active space self-consistent field–based AIMS simulations of the 1,2- dithienyl-1,2-dicyanoethene molecule, a promising molecular photoswitch.

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

Thermal disorder and phonon softening in the ferroelectric phase transition of lead titanate

We report a molecular dynamics study of ab initio quality of the ferroelectric phase transition in crystalline PbTi⁢O3. We model anharmonicity accurately in terms of potential energy and polarization surfaces trained on density functional theory data with modern machine learning techniques. Our simulations demonstrate that the transition has a strong order-disorder character, in agreement with diffraction experiments, and provide fresh insight into the approach to equilibrium across the phase transition. We find that the emergence and disappearance of the macroscopic polarization is driven by dipolar switching at the nanometer scale. We also computed the infrared optical absorption spectra in both the ferroelectric and the paraelectric phases, finding good agreement with the experimental Raman frequencies. Often, the almost ideal displacive character of the soft mode detected by Raman scattering in the paraelectric phase has been contrasted with the order-disorder character of the transition suggested by diffraction experiments. We settle this issue by showing that the soft mode coexists with a strong Debye relaxation associated with thermal disordering of the dipoles. The Debye relaxation feature is centered at zero frequency and appears near the transition temperature in both the ferroelectric and the paraelectric phases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Affinity for OH – Produces Four-Coordinated Zn 2+ Impurities in Hydrated Amorphous Calcium Carbonate

Using ab initio based molecular dynamics and electronic structure calculations, we show that Zn impurities in hydrated amorphous calcium carbonate (ACC) have a much lower coordination number than other divalent impurities due to covalent interactions between the 3d Zn shell and the oxygen atoms of the carbonate and water groups. Further, the local structure around Zn in ACC, including the predicted low coordination number, is confirmed by X-ray absorption spectroscopy of synthetic Zn-bearing ACC. The strong Zn–O chemical interaction leads to substantial water dissociation and slightly disrupts the hydrogen bonding network. Implications of Zn 2+ incorporation for ACC stability are discussed.

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

SEI formation mechanisms and Li + dissolution in lithium metal anodes: Impact of the electrolyte composition and the electrolyte-to-anode ratio

The lithium metal battery is one of today's most promising high-energy-density storage devices. Its full-scale implementation depends on solving operational and safety issues intrinsic to the Li metal high reactivity leading to uncontrolled electrolyte decomposition and uneven Li deposition. In this work, we study the spontaneous formation of the solid electrolyte interphase (SEI) upon contact of Li metal with the electrolyte and describe the heterogeneous SEI morphological features. Multiple electrolyte formulations based on lithium bis(fluorosulfonyl)imide (LiFSI), dimethoxyethane (DME), dimethyl carbonate (DMC), 1,2,2-tetrafluoroethyl-2,2,3,3-tetrafluoropropyl ether (TTE) and bis(2,2,2-trifluoroethyl) ether (BTFE) are used. Findings include the description of the SEI evolution from dispersed LiO, LiS, LiN, and LiF clusters to a continuous and compact inorganic phase in which the LiO and LiF content depend on the presence of fluorine diluents. Further, the role of the DME ether solvent helping the growth of a “wet-SEI” is compared to that of the highly unstable carbonate DMC, which decomposes into complex radical oligomers that might contribute to further electrolyte decomposition. The impact of the electrolyte-to-anode ratio on LiFSI decomposition is highlighted. Finally, we suggest the existence of a critical LiFSI concentration and electrolyte-to-anode ratio that could potentially balance the rate of electrolyte depletion and lithium consumption.

25 ENERGY STORAGE↗