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 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↗
DOE OSTI2021

Assessing density functionals for describing methane dissociative chemisorption on Pt(110)-(2×1) surface

In this work, we explore the suitability of several density functionals with the generalized gradient approximation (GGA) and beyond for describing the dissociative chemisorption of methane on the reconstructed Pt(110)-(2×1) surface. The bulk and surface structures of the metal, methane adsorption energy, and dissociation barrier are used to assess the functionals. A van der Waals corrected GGA functional (optPBE-vdW) and a metaGGA functional with van der Waals correction (MS PBEl-rVV10) are selected for ab initio molecular dynamics calculations of the sticking probability. Our results suggest that the use of these two functionals may lead to a better agreement with existing experimental results, thus serving as a good starting point for future development of reliable machine-learned potential energy surfaces for the dissociation of methane on the Pt(110)-(2×1) surface.

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

Disentangling the evolution of electrons and holes in photoexcited ZnO nanoparticles

The evolution of charge carriers in photoexcited room temperature ZnO nanoparticles in solution is investigated using ultrafast ultraviolet photoluminescence spectroscopy, ultrafast Zn K-edge absorption spectroscopy, and ab initio molecular dynamics (MD) simulations. The photoluminescence is excited at 4.66 eV, well above the band edge, and shows that electron cooling in the conduction band and exciton formation occur in <500 fs, in excellent agreement with theoretical predictions. The x-ray absorption measurements, obtained upon excitation close to the band edge at 3.49 eV, are sensitive to the migration and trapping of holes. They reveal that the 2 ps transient largely reproduces the previously reported transient obtained at 100 ps time delay in synchrotron studies. In addition, the x-ray absorption signal is found to rise in ∼1.4 ps, which we attribute to the diffusion of holes through the lattice prior to their trapping at singly charged oxygen vacancies. Indeed, the MD simulations show that impulsive trapping of holes induces an ultrafast expansion of the cage of Zn atoms in <200 fs, followed by an oscillatory response at a frequency of ∼100 cm−1, which corresponds to a phonon mode of the system involving the Zn sub-lattice.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Raman studies of hydrogen trapped in As 4 O 6 ·2H 2 at high pressure and low temperature

In this work, Raman spectroscopic measurements of the arsenolite–hydrogen inclusion compound As 4 O 6 ·2H 2 were performed in diamond anvil cells at high pressure and variable temperature down to 80 K. The experimental results were complemented by ab initio molecular dynamics simulations and phonon calculations. Observation of three hydrogen vibrons in As 4 O 6 ·2H 2 is reported in the entire temperature and pressure range studied (up to 24 GPa). While the experiments performed with protium and deuterium at variable temperatures allowed for the assignment of two vibrons as Q 1 (1) and Q 1 (0) transitions of ortho and para spin isomers of hydrogen trapped in the inclusion compound, the origin of the third vibron could not be unequivocally established. Low-temperature spectra revealed that the lowest-frequency vibron is actually composed of two overlapping bands of A g and T 2g symmetries dominated by H 2 stretching modes as predicted by our previous density functional theory calculations. We observed low-frequency modes of As 4 O 6 ·2H 2 vibrations dominated by H 2 “librations,” which were missed in a previous study. A low-temperature fine structure was observed for the J = 0 → 2 and J = 1 → 3 manifolds of hydrogen trapped in As 4 O 6 ·2H 2 , indicating the lifting of degeneracy due to an anisotropic environment. A non-spherical distribution was captured by molecular dynamics simulations, which revealed that the trajectory of H 2 molecules is skewed along the crystallographic $\langle$111$\rangle$ direction. Last but not least, low-temperature synchrotron powder x-ray diffraction measurements on As 4 O 6 ·2H 2 revealed that the bulk structure of the compound is preserved down to 5 K at 1.6 GPa.

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

Structural disorder, sublattice melting, and thermo-elastic properties of anti-perovskite Li 3 OBr under high pressure and temperature

Here, in situ synchrotron x-ray diffraction experiments were conducted in the system Li-O-Br at pressures up to 6.5 GPa and temperatures up to the melting points. The thermal equation of state determined from the P–V–T measurements indicates that crystalline anti-perovskite Li 3 OBr exhibits similar thermo-elastic properties to common salts. At temperatures that are 50–75 °C below the crystal melting, a pre-melting behavior was inferred based on weakened Bragg intensity and a substantial increase in the diffuse scattering. Concurrent with the onset of pre-melting, an anomalous increase in the lattice volume was observed and is interpreted as anion sublattice disorder. These findings support previous ab initio molecular dynamics simulations based on a simple anion interchange mechanism, indicating that structural disorder below crystal melting temperature can be a main driving force for Li-sublattice melting and superionic transition in the Li 3 OBr anti-perovskite.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Two-temperature warm dense hydrogen as a test of quantum protons driven by orbital-free density functional theory electronic forces

We consider a steady-state (but transient) situation in which a warm dense aggregate is a two-temperature system with equilibrium electrons at temperature Te, ions at Ti, and Te ≠ Ti. Such states are achievable by pump–probe experiments. For warm dense hydrogen in such a two-temperature situation, we investigate nuclear quantum effects (NQEs) on structure and thermodynamic properties, thereby delineating the limitations of ordinary ab initio molecular dynamics. We use path integral molecular dynamics (PIMD) simulations driven by orbital-free density functional theory (OFDFT) calculations with state-of-the-art noninteracting free-energy and exchange-correlation functionals for the explicit temperature dependence. We calibrate the OFDFT calculations against conventional (explicit orbitals) Kohn–Sham DFT. We find that when the ratio of the ionic thermal de Broglie wavelength to the mean interionic distance is larger than about 0.30, the ionic radial distribution function is meaningfully affected by the inclusion of NQEs. Moreover, NQEs induce a substantial increase in both the ionic and electronic pressures. This confirms the importance of NQEs for highly accurate equation-of-state data on highly driven hydrogen. For Te > 20 kK, increasing Te in the warm dense hydrogen has slight effects on the ionic radial distribution function and equation of state in the range of densities considered. In addition, we confirm that compared with thermostatted ring-polymer molecular dynamics, the primitive PIMD algorithm overestimates electronic pressures, a consequence of the overly localized ionic description from the primitive scheme.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗