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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2023

Energetics of high temperature degradation of fentanyl into primary and secondary products

Fentanyl is a synthetic opioid used for managing chronic pain. Due to its higher potency (50–100×) than morphine, fentanyl is also an abused drug. A sensor that could detect illicit fentanyl by identifying its thermally degraded fragments would be helpful to law enforcement. While experimental studies have probed the thermal degradation of fentanyl, little theoretical work has been done to understand the mechanism. Here, we studied the thermal degradation pathways of fentanyl using extensive ab initio molecular dynamics simulations combined with enhanced sampling via multiple-walker metadynamics. We calculated the free energy profile for each bond suggested earlier as a potential degradation point to map the thermodynamic driving forces. In conclusion, we also estimated the forward attempt rate of each bond degradation reaction to gain information about degradation kinetics.

Poudel, Bharat↗
DOE OSTI2023

Locking water molecules via ternary O–H⋯O intramolecular hydrogen bonds in perhydroxylated closo -dodecaborate

A multitude of applications related to perhydroxylated closo-dodecaborate B 12 (OH) 12 2− in the condensed phase are inseparable from the fundamental mechanisms underlying the high water orientation selectivity based on the base B 12 (OH) 12 2− . Herein, we directly compare the structural evolution of water clusters, ranging from monomer to hexamer, oriented by functional groups in the bases B 12 H 12 2− , B 12 H 11 OH 2− and B 12 (OH) 12 2− using multiple theoretical methods. A significant revelation is made regarding B 12 (OH) 12 2− : each additional water molecule is locked into the intramolecular hydrogen bond B–O–H ternary ring in an embedded form. This new pattern of water cluster growth suggests that B–(H–O)⋯H–O interactions prevail over the competition from water–hydrogen bonds (O⋯H–O), distinguishing it from the behavior observed in B 12 H 12 2− and B 12 H 11 OH 2− bases, in which competition arises from a mixed competing model involving dihydrogen bonds (B–H⋯H–O), conventional hydrogen bonds (B–(H–O)⋯H–O) and water hydrogen bonds (O⋯H–O). Through aqueous solvation and ab initio molecular dynamics analysis, we further demonstrate the largest water clusters in the first hydrated shell with exceptional thermodynamic stability around B 12 (OH) 12 2− . These findings provide a solid scientific foundation for the design of boron cluster chemistry incorporating hydroxyl-group-modified borate salts with potential implications for various applications.

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

Carbon Kagome nanotubes—quasi-one-dimensional nanostructures with flat bands

In recent years, a number of bulk materials and heterostructures have been explored due their connections with exotic materials phenomena emanating from flat band physics and strong electronic correlation. The possibility of realizing such fascinating material properties in simple realistic nanostructures is particularly exciting, especially as the investigation of exotic states of electronic matter in wire-like geometries is relatively unexplored in the literature. Motivated by these considerations, we introduce in this work carbon Kagome nanotubes (CKNTs)—a new allotrope of carbon formed by rolling up Kagome graphene, and investigate this material using specialized first principles calculations. We identify two principal varieties of CKNTs—armchair and zigzag, and find both varieties to be stable at room temperature, based on ab initio molecular dynamics simulations. CKNTs are metallic and feature dispersionless states (i.e., flat bands) near the Fermi level throughout their Brillouin zone, along with an associated singular peak in the electronic density of states. We calculate the mechanical and electronic response of CKNTs to torsional and axial strains, and show that CKNTs appear to be more mechanically compliant than conventional carbon nanotubes (CNTs). Additionally, we find that the electronic properties of CKNTs undergo significant electronic transitions—with emergent partial flat bands and tilted Dirac points—when twisted. We develop a relatively simple tight-binding model that can explain many of these electronic features. We also discuss possible routes for the synthesis of CKNTs. Overall, CKNTs appear to be unique and striking examples of realistic elemental quasi-one-dimensional materials that may display fascinating material properties due to strong electronic correlation. Distorted CKNTs may provide an interesting nanomaterial platform where flat band physics and chirality induced anomalous transport effects may be studied together.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Topology driven and soft phonon mode enabled Na-ion diffusion in quaternary chalcogenides, Na 3 ZnGaX 4 (X = S, and Se)

The compounds Na 3 ZnGaX 4 (X = S, Se) are potential solid electrolyte materials in sodium-based batteries, which have certain advantages over oxide materials and have shown significant ionic conductivity at ambient temperature. In this paper, we bring out atomic-level features of the diffusion process in these new materials using the microscopic techniques of inelastic neutron scattering (INS), quasielastic neutron scattering (QENS), and ab initio molecular dynamics (AIMD) simulations. The insights obtained from these techniques are unique and not available from other macroscopic experiments. Neutron scattering experiments have been performed at temperatures from 100 to 700 K. The simulations have been carried out up to 900 K. We have calculated the phonon spectra and the space–time correlation functions and found good agreement with the results of the neutron scattering experiments. The simulations enable detailed analysis of the atomic-site dependent dynamical information. We observe low-energy phonon modes of ~6 meV involving the vibrations of certain Na atoms in the lattice. This reveals that the Na at 32g Wyckoff sites (Na2) has sufficiently shallow potential among the two available crystallographic sites. This shallow potential facilitates diffusion. Furthermore, the specific structural topology of the network of interconnected zig-zag chains of the Na2 atomic sites provides the low-barrier energy pathways for diffusion. A small fraction of vacancy defects appears essential for diffusion. We further observe that the Na2 atoms undergo jump-like diffusion to the vacant next or the 2nd next neighbour sites at ~4 Å. While the QENS experiments reveal the jump-like diffusion and its time scale, detailed analysis of the AIMD simulations shows that the jumps appear mostly along zig-zag chains of the Na2 sites in the tetragonal ab-plane, as well as between the chains along the c-axis. In conclusion, the net diffusion is essentially 3-dimensional, with little anisotropy despite the anisotropy of the tetragonal crystal structure.

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

Structure, bonding and electronic characteristics of amorphous Se

Ovonic threshold switching (OTS) selector can effectively improve the storage density and suppress the leakage current of advanced phase-change memory. As a prototypical OTS material, amorphous GeSe is widely investigated. But the attentions paid to amorphous Se (i.e., the functional constituent in amorphous GeSe) is very limited up to now. Here, in this work, we have explored the structure, bonding and electronic characteristics of amorphous Se using ab initio molecular dynamics simulations. The results reveal that the Se atoms in amorphous Se tend to form the 2-coordinated configurations, and they connect with each other to form the long chains. The fraction of vibrational density of state located in high frequency range is relatively large, and the formation energy of Se-Se bond is as large as 4.44 eV, hinting that the Se-Se bonds in chains possess a high stability. In addition, the mid-gap state related to the OTS behavior is also found in the amorphous Se despite of the small proportion. Our findings enrich the knowledge of the amorphous Se, which contributes to the applications of Se-based OTS selectors.

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

Modeling separation of lanthanides via heterogeneous ligand binding

Individual lanthanide elements have physical/electronic/magnetic properties that make each useful for specific applications. Several of the lanthanides cations (Ln 3+ ) naturally occur together in the same ores. They are notoriously difficult to separate from each other due to their chemical similarity. Predicting the Ln 3+ differential binding energies (ΔΔE) or free energies (ΔΔG) at different binding sites, which are key figures of merit for separation applications, will help design of materials with lanthanide selectivity. We apply ab initio molecular dynamics (AIMD) simulations and density functional theory (DFT) to calculate ΔΔG for Ln 3+ coordinated to ligands in water and embedded in metal–organic frameworks (MOFs), and ΔΔE for Ln 3+ bonded to functionalized silica surfaces, thus circumventing the need for the computational costly absolute binding (free) energies ΔG and ΔE. Perturbative AIMD simulations of water-inundated simulation cells are applied to examine the selectivity of ligands towards adjacent Ln 3+ in the periodic table. Static DFT calculations with a full Ln 3+ first coordination shell, while less rigorous, show that all ligands examined with net negative charges are more selective towards the heavier lanthanides than a charge-neutral coordination shell made up of water molecules. Amine groups are predicted to be poor ligands for lanthanide-binding. Finally, we also address cooperative ion binding, i.e., using different ligands in concert to enhance lanthanide selectivity.

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

Exploring Cr and molten salt interfacial interactions for molten salt applications

Molten salts play an important role in various energy-related applications such as high-temperature heat transfer fluids and reaction media. However, the extreme molten salt environment causes the degradation of materials, raising safety and sustainability challenges. A fundamental understanding of material–molten salt interfacial evolution is needed. This work studies the transformation of metallic Cr in molten 50/50 mol% KCl–MgCl 2 via multi-modal in situ synchrotron X-ray nano-tomography, diffraction and spectroscopy combined with density functional theory (DFT) and ab initio molecular dynamics (AIMD) simulations. Notably, in addition to the dissolution of Cr in the molten salt to form porous structures, a δ-A15 Cr phase was found to gradually form as a result of the metal–salt interaction. This phase change of Cr is associated with a change in the coordination environment of Cr at the interface. DFT and AIMD simulations provide a basis for understanding the enhanced stability of δ-A15 Cr vs. bcc Cr, by revealing their competitive phase thermodynamics at elevated temperatures and probing the interfacial behavior of the molten salt at relevant facets. This study provides critical insights into the morphological and chemical evolution of metal–molten salt interfaces. Finally, the combination of multimodal synchrotron analysis and atomic simulation also offers an opportunity to explore a broader range of systems critical to energy applications.

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

Nature of molybdenum carbide surfaces for catalytic hydrogen dissociation using machine-learned potentials: an ensemble-averaged perspective

Molybdenum carbides with an electronic structure similar to noble metals have gained attention as a promising low-cost catalyst for biomass valorization and the hydrogen evolution reaction. However, our fundamental understanding of the catalyst surface and how different phases of these catalysts behave at varying reaction conditions is limited to ground state density functional theory calculations as ab initio molecular dynamics (AIMD) is computationally prohibitive at relevant length and time scales. Here, in this work, we train a multi-atomic cluster expansion (MACE) machine-learned interatomic potentials (MLIP) to study hydrogen dissociation and dynamics over Mo, δ-MoC, α-Mo 2 C, and β-Mo 2 C surfaces at varying temperatures and hydrogen partial pressures. Our simulations identify unique and different molecular and atomic hydrogen adsorption sites on different surfaces that do not depend on the temperature. At low hydrogen pressures, the surface coverage is monolayer, which transitions to two-layer adsorption at higher pressures. We find that atomic hydrogen diffusion and recombinations are preferred over molybdenum atom hollow sites, while the diffusion over carbon-terminated facets was negligible, signifying particularly strong C–H interactions. In contrast, molecular hydrogen adsorption occurs mostly atop Mo or the bridging sites. At a comparable hydrogen loading, β-Mo 2 C (001) is the most active surface for hydrogen dissociation reaction. This work provides insights into the dynamic nature of the hydrogen dissociation chemistry and the diversity of hydrogen adsorption sites on molybdenum carbides.

08 HYDROGEN↗
DOE OSTI2024

Revealing the structure and electronic characteristics of Te-rich threshold switching materials for high-density integration

Ovonic threshold switching (OTS) selectors play an important role in the integration of advanced three-dimensional memory. Selectors based on tellurium (Te)-containing materials exhibit significant promise due to their low threshold voltages and superior consistency. Here we have theoretically studied the structure and electronic properties of a typical OTS material, amorphous GeTe 6 , to explore the switching mechanisms using ab initio molecular dynamics simulations. The results indicate that Ge atoms tend to bond with Te atoms, forming stable chemical bonds. The Te-centered clusters are predominantly in the form of distorted octahedrons, while the Ge-centered clusters are in the form of both octahedrons and tetrahedrons. Notably, the proportion of tetrahedrons within the 4-coordinated Ge-centered clusters reaches an impressive 66.9%. These tetrahedrons are randomly dispersed throughout the simulated cell, leading to a stable amorphous configuration. The mid-gap state observed in the mobility bandgap originates from the atomic chain composed of both over-coordinated Ge and Te atoms. It is the inherent stability of the chemical environment within amorphous GeTe 6 that enables it to maintain its amorphous phase under a repeated threshold voltage, a characteristic that distinguishes it from non-OTS materials, such as amorphous tellurium. Furthermore, our findings provide an in-depth understanding of the structure and electronic characteristics of amorphous GeTe 6 , which can promote the design and application of the Te-rich threshold switching materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Structural incorporation into goethite fractionates rare earth elements

Up to 20% of rare earth elements (REEs) in ion adsorption deposits (IADs) are associated with iron oxide minerals, primarily goethite. Often termed “non-extractable”, goethite-hosted REEs are thought to be structurally incorporated into the mineral lattice. The large mismatch in size and charge density between REEs (ionic radius, r = 0.86–1.03 Å) and Fe 3+ (r = 0.65 Å), however, makes direct substitution energetically unfavorable. To determine REE compatibility with and incorporation into goethite on the atomic level, we used X-ray pair distribution function analysis (PDF) and LIII-edge extended X-ray absorption fine structure (EXAFS) spectroscopy. The compatibility of REEs with goethite and the precursor ferrihydrite (FH) is Lu ≥ Yb ≫ Dy > Nd for both phases. Nd and Dy primarily form secondary amorphous phases, with <30% Nd and Dy incorporated into goethite. In the FH precursor at pH 6.8, Yb and Lu assumed a local REE-OOH like structure with next nearest neighbor Fe. The Yb, Lu, and Nd-FH samples were also matured at ambient conditions for 100 days; despite the presence of only ∼5% goethite, Lu and Yb were 42% and 100% in goethite-like structural environments, respectively, whereas the PDF and EXAFS of Nd showed little evidence of any incorporation. Using ab initio molecular dynamics (AIMD) to model the EXAFS, we determined the presence of protonated Fe vacancies, edge-sharing with structural Lu and Yb, likely helped accommodate these REEs in the goethite structure. Incorporation into Fe oxyhydroxides thus potentially fractionates the REEs during weathering associated with formation of lateritic and ion adsorption deposits.

Mergelsberg, Sebastian T. [Pacific Northwest Natio↗
DOE OSTI2020

New interaction potentials for borate glasses with mixed network formers

Here, we adapt and apply a recently developed optimization scheme used to obtain effective potentials for aluminosilicate glasses to include the network former boron into the interaction parameter set. As input data for the optimization, we used the radial distribution functions of the liquid at high temperature generated by ab initio molecular dynamics simulations, and density, coordination, and elastic modulus of glass at room temperature from experiments. The new interaction potentials are shown to reliably reproduce the structure, coordination, and mechanical properties over a wide range of compositions for binary alkali borates. Furthermore, the transferability of these new interaction parameters allows mixing to reliably reproduce the properties of various boroaluminate and borosilicate glasses.

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

Pressure-induced transformation of Nb 2 O 5 under shock compression from first principles

Ab initio molecular dynamics (AIMD) simulations were carried out to investigate the equation of state of Nb 2 O 5 and its pressure-density relationship under shock conditions. Here, the focus of this study is on the monoclinic B-Nb 2 O 5 (C2/c) polymorph. Enthalpy calculations from AIMD trajectories at 300 K show that the pressure-induced transformation between the thermody-namically most stable crystalline monoclinic parent phase H-Nb 2 O 5 (P2/m) and B-Nb 2 O 5 occurs at ~1.9 GPa. This H→B transition is energetically more favorable than the H→L(Pmm2) pressure-induced transition recently observed at ~5.9-9.0 GPa. The predicted shock properties of Nb 2 O 5 polymorphs are also compared to their Nb and NbO 2 counterparts to assess the impact of niobium oxidation on shock response.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗