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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Vapor-phase self-assembly for generating thermally stable single-atom catalysts

Preparation of thermally stable metal single-atom catalysts (SACs) is a challenge in heterogeneous catalysis, especially on conventional supports that provide a weak metal-support interaction. In this work, we report that a modified support MgAl 2 O 4 can stabilize Pt single atoms by a mechanism of vapor-phase self-assembly in a high-temperature treatment (800°C, air). The experimental results on the formation mechanism and the structure are validated by DFT and ab initio molecular dynamics simulations. We infer that stable triangular K 3 O 3 structures help stabilize Pt single atoms at high temperatures in oxidizing conditions, exhibiting excellent reactivity for methane oxidation. The obtained Pt/K/MgAl 2 O 4 SAC presents excellent stability in methane oxidation after steam treatment at elevated temperatures, whereas the Pt/MgAl 2 O 4 nanocatalyst suffers from rapid deactivation due to Pt nanoparticle growth. Now this work paves the way for preparing thermally stable and highly active SACs using conventional high-surface-area supports, despite the weak metal-support interaction.

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

Computational and experimental analysis of structural and Thermophysical properties of LiX-KX (X = chloride, iodide, or bromide) molten salts

Molten salts containing lithium (Li) and potassium (K) halides, including chlorides (Cl), bromides (Br), and iodides (I), are pivotal in advanced technological applications including nuclear reactor technologies, thermal batteries, and industrial pyrochemical processes. Despite their significance, the availability of the temperature-dependent structural and thermophysical properties of these molten salts remain scarce. Here, this study addresses this knowledge gap by predicting and measuring the key properties such as coordination number, radial distribution function, density, heat capacity, and volumetric thermal expansion of molten LiI-KI, LiBr-KBr and LiCl-KCl eutectic mixtures. Utilizing ab initio molecular dynamics (AIMD) simulations, we predict these properties of LiX-KX (X=Br,I, Cl) across various temperatures, with chlorides as a benchmark. Experimental measurements using Archimedes methods and Differential Scanning Calorimetry validate the densities and heat capacities of these molten salts. Our comprehensive analysis of the structural and thermophysical properties provides critical insights into the behavior of these molten salts at various temperatures, enhancing the understanding necessary for their application in advanced technologies.

36 - MATERIALS SCIENCE↗
DOE OSTI2021

Combined density functional theory/kinetic Monte Carlo investigation of surface morphology during cycling of Li-Cu electrodes

In this work, we employed a combined density functional theory and kinetic Monte Carlo approach in order to investigate the nanoscale behavior of lithium stripping and plating on lithium-copper nanoalloys. Density functional theory calculations were used to develop an atom-by-atom understanding of the nucleation behavior of lithium on copper surfaces with different surface morphologies. We found that nanometric pits resulted in significantly more favorable lithium deposition when compared to pristine or lightly-defective Cu surfaces, which is in agreement with previous experimental studies. We then used density functional theory and constrained ab initio molecular dynamics to develop interaction and reaction parameters for a larger-scale kinetic Monte Carlo model in order to provide insight into the stripping and plating behavior of nanostructured Li-Cu nanoalloy anodes. Substantial losses in capacity and strong indications of morphological change were observed for a wide range of Cu:Li ratios over 100 cycles. We observed a complex relationship between the nanoparticle composition and the properties of the formed solid-electrolyte interphase.

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

Ionic Interdiffusion at Cathode|Solid-Electrolyte Interface: A Machine Learning–Assisted Multiscale Investigation and Mitigation Strategies

Future lithium batteries are expected to use solid electrolytes to achieve higher energy density and fast charge capabilities. However, most solid electrolytes are thermodynamically unstable against layered oxide cathodes. In this study, the stability of LiCoO2 (LCO) cathode with Li10GeP2S12 (LGPS) solid electrolyte is investigated using ab initio molecular dynamics (AIMD) and machine learning molecular dynamics (MLMD). The propensity of ionic interdiffusion, formation of a passivating interphase layer, and corresponding decay in cell performance is addressed using a continuum model. Large-scale MLMD simulations confirm that the LCO|LGPS interface permits interdiffusion of cobalt (Co) and other ionic species, leading to the formation and growth of a resistive interphase and to dramatic capacity fade even in the first cycle. We examine the literature evidence that incorporating a thin layer of LiNb0.5Ta0.5O3 (LNTO) between LCO and LGPS prevents the interdiffusion of ions. Atomistic simulations suggest that substituting lithium (Li) in LNTO with Co is thermodynamically unfavorable, thereby inhibiting ionic interdiffusion. The stable Nb5+/Ta5+ states form a rigid metal-oxide framework, which consequently also prevents the substitution of niobium (Nb) or tantalum (Ta). However, continuum-level analysis suggests that the higher mechanical stiffness of LNTO can lead to interfacial delamination between the LCO and LNTO. This phenomenon reduces the effectiveness of the protective layer. This paper, therefore, highlights the need to develop novel interlayers that balance low ionic interdiffusion with low mechanical stiffness.

Ncube, Musawenkosi K.↗
DOE OSTI2024

Additively manufactured Ni-20Cr to V functionally graded material: Computational predictions and experimental verification of phase formations

Throughout this work, a database for the Cr-Ni-V system was constructed by modeling the binary Cr-V and ternary Cr-Ni-V systems using the CALPHAD approach aided by density functional theory (DFT)-based first-principles calculations and ab initio molecular dynamics (AIMD) simulations. To validate this new database, a functionally graded material (FGM) using Ni-20Cr and elemental V was fabricated using directed energy deposition additive manufacturing (DED AM) and experimentally characterized. The deposited Ni-20Cr was pure fcc phase, while increasing the amount of V across the gradient resulted in the formation of sigma phase, followed by the bcc phase. The experimentally measured phase data was compared with computational predictions made using a Cr-Ni-V thermodynamic database from the literature as well as the database developed in the present work. The newly developed database was shown to better predict the experimentally observed phases due to its accurate modeling of binary systems within the database and the ternary liquid phase, which is critical for accurate Scheil calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Nitrogen: A promising doping strategy for high-performance ovonic threshold switching selectors

The Ovonic Threshold Switching (OTS) selector serves as an essential component in the development of three-dimensional high-density memory integration technology. Nevertheless, the state-of-the-art high-performance OTS materials usually contain toxic elements such as arsenic (As), posing significant risks to both environmental and human health. Nitrogen (N), which belongs to the same group as arsenic (As), has emerged as a highly promising alternative for As doping. However, the underlying mechanisms that govern N-based OTS materials have not yet been extensively investigated. In this study, we delve into the effects of N doping on the structural, bonding, and electronic properties of amorphous GeSe (a-GeNSe) by ab initio molecular dynamics simulations to bridge the knowledge gap. Our findings indicate that upon N doping in a-GeSe, the formation of robust Ge-N bonds, along with N-centered tetrahedral and triangular structures, resulting in the sluggish atomic movement that enhances the thermal stability and endurance of a-GeNSe. The OTS characteristics are significantly influenced by the material’s electronic band structure, and thus the relatively slow performance drift can be attributed to the stabilization of mid-gap states, a result of N doping which effectively slows down the aging process of chalcogenide glass. Moreover, the increased mobility gap in a-GeNSe raises the threshold voltage (V th ), making it more compatible with commercially available phase-change memory materials. Furthermore, our findings reveal the extensive impact of the N element on a typical OTS material and offer valuable perspectives for alternative doping strategies that could potentially supplant As practices.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Adsorption of hydroxamic acid ligands for improved extraction of rare earth elements from monazite ores

Efficient separation of rare earth element (REE) ores via froth flotation requires the development of novel ligands with enhanced adsorption capacity and selectivity. To realize these advances, understanding the mechanisms underlying interactions between the ligand and mineral surfaces is essential. This study systematically evaluates the adsorption behavior of alkyl- and aromatic alkyl-substituted hydroxamic acid ligands on monazite surfaces using complementary spectroscopic techniques, including UV–visible (UV–vis) spectroscopy, Raman spectroscopy, infrared spectroscopy, and vibrational sum frequency generation (SFG) spectroscopy, together with the ab initio molecular dynamics (AIMD) simulations. Among the studied ligands, octanohydroxamic acid (OHA) and 4-ethoxy-N,2-dihydroxybenzamide (EDHBA) exhibit high adsorption capacity under basic pH (8–10) by forming multilayers on the surface. OHA has a higher equilibrium adsorption capacity compared to EDHBA, but it forms a less stable multilayer susceptible to disruption in the presence of interfering ions. AIMD results show that OHA adopts a single stable chelating geometry, while EDHBA exhibits multiple binding modes involving distinct interactions with La surface atoms and phosphate-bound oxygens, resulting in more complex adsorption kinetics. The variations in surface binding and intermolecular interactions observed between alkyl and aromatic molecules influence the differences in adsorption kinetics, equilibrium adsorption capacities on the mineral surface, and their flotation performance. This work provides valuable insight into the adsorption mechanism of ligands at mineral interfaces, which is crucial for guiding the design of new ligands with enhanced separation performance.

Zhou, Muchu [ORNL] (ORCID:0000000182650215)↗
DOE OSTI2024

Structural and electronic characteristics of amorphous Ge 8 Sb 2 Te 11

GeTe-rich phase-change materials can be utilized in rewriteable optical memory due to the large contrast in reflectivity between amorphous and crystalline phases. Here we explored the structure and electronic properties of amorphous Ge 8 Sb 2 Te 11 using ab initio molecular dynamics simulations. The results indicate that amorphous Ge 8 Sb 2 Te 11 is dominantly composed of Ge-, Sb- and Te-centered octahedrons with distortions, while 30.4 % of Ge-centered clusters are in the form of tetrahedrons which are randomly distributed. The 5-fold rings possess a large proportion, and the Gesingle bondTe and Sbsingle bondTe bonds present larger formation energies than other bonds, leading to the ABAB bonding sequence (A: Ge and Sb, B: Te). The lone pair electrons locating at the opposite direction of bonds possess a large fraction of 14.8 %, which may enhance the distortions of local clusters. Importantly, these structural properties lead to the high stability of amorphous Ge 8 Sb 2 Te 11 and thus long data retention in the information storage.

36 MATERIALS SCIENCE↗
DOE OSTI2024

First-principles investigation of the thermophysical properties of NaCl, PuCl 3 , and NaCl-PuCl 3 Molten salts

Molten salts have a variety of applications that span the nuclear and solar industries, and which involve thermal storage and heat transfer. There is a present knowledge gap in the thermophysical properties of molten salts, which limits the readiness level of molten salt applications. This is especially pertinent for molten salt reactors where the fissile material is dissolved within the molten salt. Ab initio molecular dynamics is a common method employed to investigate the structural and thermophysical properties at elevated temperatures. NaCl-PuCl 3 is a candidate fuel salt in molten salt reactors. The scope of this study is to investigate the seven unique compositions and a range of temperatures of NaCl-PuCl 3 and calculate the density, heat capacity, compressibility, enthalpy of mixing, and coefficient of thermal expansion. Within this work, the van der Waals (vdW) interactions are primarily handled with the vdW-DF2 functional but spot-checked with both the dDsC and DFT-D3 methods. The calculated densities are in good agreement with the NaCl and the eutectic compositions, and no experimental values exist for PuCl 3 at these temperatures. The densities are fit to a first-order Redlich-Kister expansion as a function of both composition and temperature. The heat capacity in the literature is scattered but the calculated values agree well with one of the experimental results. Here, the heat capacity increases as a linear function with respect to concentration at a rate of about 7.5 J/mol-K per 10 mol% PuCl 3 .

36 MATERIALS SCIENCE↗
DOE OSTI2022

Enhanced ionic conductivity and lack of paddle-wheel effect in pseudohalogen-substituted Li argyrodites

Superionic conductors are key to the development of safe and high-energy-density, all-solid-state batteries. Using a combined theoretical and experimental approach, we explore the feasibility of increasing the ionic conductivity through pseudohalogen substitution in the Li argyrodite structure. Under the guidance of calculated thermodynamic stability, BH 4 -substituted Li argyrodite, Li 5.91 PS 4.91 (BH 4 ) 1.09 , was successfully synthesized via a mechanochemical method. As-synthesized BH 4- substituted Li argyrodite displays an ionic conductivity of 4.8 mS/cm at 25°C. Ab initio molecular dynamics simulation trajectory analysis was used to investigate how BH 4 facilitates Li-ion diffusion and indicates only a weak correlation with the B–H bond motion. Here, we find that the enhanced conductivity mainly originates from the weak interaction between Li and BH 4 and find no evidence of a paddle-wheel effect from the polyanion. This work provides insight into how cluster ions enhance Li diffusion and systematically describes how to explore superionic conductors with pseudohalogen substitution.

25 ENERGY STORAGE↗
DOE OSTI2020

Fluorinated co-solvent promises Li-S batteries under lean-electrolyte conditions

Practical lithium-sulfur batteries stipulate the use of a lean-electrolyte and a high Coulombic efficiency of the lithium metal anode. In this work, we employ 1,1,2,2-tetrafluoroethyl 2,2,2-trifluoroethyl ether as a co-solvent in the electrolyte of Li-S batteries to meet the demands. The co-solvent fosters the formation of a LiF-rich solid electrolyte interphase on lithium metal anode, as revealed by Ab initio molecular dynamics and Femtosecond Stimulated Raman spectroscopy. The co-solvent results in an average Coulombic efficiency of 99.4% for lithium plating/stripping cycles. Full cells with a capacity ratio of 2 between lithium anode and the S@PAN cathode (3 gs cm –2 ) exhibited a stable cycle life over 100 cycles at an electrolyte/sulfur ratio of 2 µL mg –1 , validating the high Coulombic efficiency of the lithium anode and demonstrating the compatibility of the electrolyte with both electrodes. To enhance the energy density, we prepared a hybrid cathode composed of 45 wt.% VS 2 mixed with ZnS-coated Li 2 S@graphene as the cathode. Based on the mass of both electrodes and the electrolyte, the full cell delivers an energy density of 483 Wh kg –1 , which demonstrates viable Li-S batteries with the lean-electrolyte.

25 ENERGY STORAGE↗
DOE OSTI2021

Computational study of the electrostatic potential and charges of multivalent ionic liquid molecules

The electrostatic interactions in ionic liquids (ILs) can be difficult to quantify, but this information is very valuable for improving intermolecular potentials for performing molecular-level simulations. Although polarizable forcefields can provide improved accuracy, the computational efficiency of fixed charge molecular models is often preferred for larger systems. Here, we propose a multi-scale screening approach for analyzing the charge scaling behavior and other electrostatic properties of ILs, which is more attainable than previous analysis methods that rely on the IL crystalline structures. Based on isolated molecules and small clusters, we study the electrostatic potential properties and charges of multivalent IL molecules using density-functional theory calculations and ab initio molecular dynamics simulations. The charges of the IL molecules are calculated and compared using several different charge estimation techniques. Although there are some differences among the estimated partial charges on the atoms, the total anion and cation charges are very similar. In conclusion, the charge scaling factors are substantial, and they are only slightly affected by the IL cluster size model.

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

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. These data are necessary for designing, constructing, and licensing future commercial molten-salt reactors. Thermophysical properties, such as melting temperature, salt stability, density, and heat capacity were ascertained using NaCl-PuCl3 (36 mol% PuCl3) and a more sodium rich composition containing 25 mol% PuCl3. The NaCl-PuCl3 salt mixture was synthesized for this study and contained 63.4 mol% NaCl, 36.3 mol% PuCl3 and was 99.7% pure. Upon heating, the NaCl–PuCl3 eutectic was stable at temperatures up to 800 °C. The onset of melting occurred at 451 ± 3 °C, and the enthalpy of fusion was determined to be 23.0 ± 1.4 kJ/mol. Heat capacity measurements in the liquid phase ranged from 107.7 to 91.3 J/mol.K, with an average value of 104.6 ± 11.4 J/mol.K between 500 and 650 °C. Three independent trials of the molten NaCl-PuCl3 salt found the density to be ?(T) = 3.8589 – 9.5342·10-4 T(°C). In addition, ab initio molecular dynamic simulations to calculate density and heat capacity values are included.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Temperature as a control knob on spin-orbit coupling

Spin-orbit coupling (SOC) governs many physical phenomena. Through ab initio molecular dynamics simulations, Lu and Sun demonstrated that structural disorder at elevated temperatures substantially reduces the effective SOC contribution that stabilizes band inversion, driving a topological-to-normal-insulator transition in Bi2Se3. Their work identifies temperature as a meaningful control parameter for SOC-mediated topology and other properties.

Liang, Liangbo [ORNL] (ORCID:0000000311990049)↗
DOE OSTI2023

The stability of FeH x and hydrogen transport at Earth’s core mantle boundary

Iron hydride in Earth’s interior can be formed by the reaction between hydrous minerals (water) and iron. Studying iron hydride improves our understanding of hydrogen transportation in Earth’s interior. Our high-pressure experiments found that face-centered cubic (fcc) FeH x (x≤1) is stable up to 165 GPa, and our ab initio molecular dynamics simulations predicted that fcc FeH x transforms to a superionic state under lower mantle conditions. In the superionic state, H-ions in fcc FeH become highly diffusive-like fluids with a high diffusion coefficient of ~3.7 × 10 -4 cm 2 s -1 , which is comparable to that in the liquid Fe-H phase. The densities and melting temperatures of fcc FeH x were systematically calculated. Similar to superionic ice, the extra entropy of diffusive H-ions increases the melting temperature of fcc FeH. Further, the wide stability field of fcc FeH enables hydrogen transport into the outer core to create a potential hydrogen reservoir in Earth’s interior, leaving oxygen-rich patches (ORP) above the core mantle boundary (CMB).

58 GEOSCIENCES↗
DOE OSTI2026

Adsorption, charge transfer and a coverage-driven transition of alkali metals on rutile TiO 2 (110)

The interaction of alkali metals with metal oxide surfaces is central to tuning surface reactivity in heterogeneous catalysis and photocatalysis. Here we present a comprehensive DFT+U study of the adsorption of alkali metals (Li, Na, K, Rb, Cs) on the (110) surface of rutile TiO 2 . At low coverage (θ = 1/8), all alkali metals bind preferentially to bridging oxygen sites with adsorption energies in the range −4.06 to −3.33 eV, transferring nearly one full electron (0.92–0.99 |e|) to the substrate and inducing Ti 4+ → Ti 3+ reduction. The excess charge localizes preferentially at subsurface Ti sites in the form of small polarons. Diffusion barriers indicate facile motion along bridging-oxygen rows, whereas inter-row hopping is strongly hindered. Coverage effects were examined systematically for potassium: adsorption energy and charge transfer per K atom decrease monotonically with increasing θ. Strikingly, a sharp energy discontinuity occurs between θ = 4/8 and θ = 5/8 (ΔE ≈ 1 eV per atom), which we identify as a coverage-driven structural transition arising from steric packing constraints and enhanced K–K electrostatic repulsion once every (1×1) surface cell is occupied. This structural transition perfectly correlates with a dramatic drop in the work function down to an ultra-low minimum of 0.84 eV at θ=5/8, followed by a metallization- driven recovery at higher coverages. Ab initio molecular dynamics simulations confirm zigzag K arrangements at moderate coverage (θ = 1/3), while at high coverage (θ = 2/3) short-range K–K correlations emerge without long-range order. These results provide atomistic insight into the structure–activity relationships underlying alkali promotion effects on oxide-supported catalysts.

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

Molten-salt-mediated carbon dioxide capture and superequilibrium utilization with ethane oxidative dehydrogenation

Existing CO 2 -mediated oxidative dehydrogenation (CO 2 -ODH) of ethane has yet to demonstrate >60% single-pass CO yield due to the intrinsic equilibrium limitations. We report a unique approach with mixed molten carbonates as a reaction medium for CO 2 -ODH, which strategically partitions the CO 2 -ODH reactions into gas and molten-salt phases and facilitates integrated CO 2 capture from power plant flue gases. An 89% CO yield was achieved at 770°C, doubling the equilibrium limitation of conventional CO 2 -ODH. The high CO yield in turn enhances ethylene formation. Further characterizations confirmed that molten-salt mediated ODH (MM-ODH) proceeds through a gas-phase cracking and molten-salt mediated reverse water-gas-shift reaction pathway. Based on this understanding, thermodynamic analysis and ab initio molecular dynamics simulations were conducted to develop general principles to optimize the molten-salt reaction medium. Process analyses confirm that MM-ODH has the potential to be significantly more efficient for CO 2 capture and utilization than conventional CO 2 -ODH.

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

An ab initio free energy study of the reaction mechanism and rate-limiting steps of CO 2 capture by aqueous glycine

Aqueous amino acids are promising absorbents for direct air capture (DAC) of CO 2 . Herein, we investigate the possibility of kinetic control of CO 2 absorption with aqueous anionic glycine (GLY - ) by employing extensive ab initio molecular dynamics simulations, free energy analysis, and reaction rate theory. We find that first GLY - binds to CO 2 by overcoming a barrier (7.4 kcal/mol) to form a zwitterion intermediate, which then releases a proton by overcoming a similar barrier. Despite the similarity in the barrier, zwitterion formation appears to be the rate-limiting step because it is two orders of magnitude slower (microseconds) than the proton release step. This is predominantly due to stronger nonequilibrium solvent effects for the former that cause many barrier-recrossing events and effectively slow down the reaction rate. Such a detailed fundamental understanding of the amino acid-based CO 2 -absorption mechanism and rates is key to improving the kinetic efficiency of DAC technology.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗