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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Liquid Structure and Transport Properties of the Deep Eutectic Solvent Ethaline

A range of techniques including physical property measurements, neutron scattering experiments, ab initio molecular dynamics, and classical molecular dynamics simulations are used to probe the structural, thermodynamic, and transport properties of a deep eutectic solvent comprised of a 1:2 molar ratio of choline chloride and ethylene glycol. This mixture, known as Ethaline, has many desirable properties for use in a range of applications, and therefore, understanding its liquid structure and transport properties is of interest. Simulation results are able to capture experimental densities, diffusivities, viscosities, and structure factors extremely well. Here, the solvation environment is dynamic and dominated by different hydrogen bonding interactions. Dynamic heterogeneities resulting from hydrogen bonding interactions are quantified. Rotational dynamics of molecular dipole moments of choline and ethylene glycol are computed and found to exhibit a fast and slow mode.

25 ENERGY STORAGE↗
DOE OSTI2021

Specific Ion Solvation and Pairing Effects in Glycerol Carbonate

Identifying the driving forces behind the solvation of inorganic salts by nonaqueous solvents is an important step in the development of green solvents. Here we focus on one promising solvent: glycerol carbonate (GC). Using ab initio molecular dynamics simulations, we build upon our previous work by detailing glycerol carbonate’s interactions with a series of anions, a lithium ion, and the LiF ion pair. Through these investigations, we highlight the changes in solvation behavior as the anion size increases, the competition of binding shown by lithium for the oxygens of GC, and the behavior of the LiF ion pair in a GC solution. These results in this work indicate the importance of the cation’s identity in ion-pairing structure and dynamics and lend insight into the key factors behind the specific ion effects seen in GC.

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

Exploring the Unusual Reactivity of the Hydrated Electron with CO 2

Many questions remain about the reactions of the hydrated electron despite decades of study. Of particular note is the fact that they do not appear to follow the Marcus theory of electron transfer reactions, a feature that has yet to be explained. To investigate these issues, we use ab initio molecular dynamics (AIMD) simulations to investigate a one of the better studied reactions, the hydrated electron reduction of CO 2 . The rate constant for the hydrated electron-CO 2 reaction complex to react to form CO 2 - is, for the first time, estimated from AIMD simulations. Results at 298 and 373 K show the rate constant is insensitive to temperature, consistent with the low measured activation energy for the reaction, and the implications of this behavior are examined. The sampling provided by the simulations yields insight into the reaction mechanism. The reaction is found to involve both solvent reorganization and changes in the carbon dioxide structure. The latter lead to significant vibrational excitation of the bending and symmetric stretch vibrations in the CO 2 - product, indicating the reaction is vibrationally nonadiabatic. The former is estimated from calculation of an approximate collective solvent coordinate and the free energy in this coordinate is determined. Furthermore, these results indicate that AIMD simulations can reasonably estimate hydrated electron reaction activation energies and provide new insight into the mechanism that can help illuminate the features of this unusual chemistry.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2025

Modeling the Behavior of Complex Aqueous Electrolytes Using Machine Learning Interatomic Potentials: The Case of Sodium Sulfate

Understanding the structure and thermodynamics of solvated ions is essential for advancing applications in electrochemistry, water treatment, and energy storage. While ab initio molecular dynamics methods are highly accurate, they are limited by short accessible time and length scales whereas classical force fields struggle with accuracy. Herein, we explore the structure and thermodynamics of complex monovalent-divalent ion pairs using Na 2 SO 4 (aq) as a case study by applying a machine learning interatomic potential (MLIP) trained on density functional theory (DFT) data. Our MLIP-based approach reproduces key bulk properties such as density and radial distribution functions of water. We provide the hydration structure of the sodium and sulfate ions in the 0.1–2 M concentration range and the one-dimensional and two-dimensional potentials of mean force for the sodium–sulfate ion pairing at the low concentration limit (0.1 M), which are inaccessible to DFT. At low concentrations, the sulfate ion is strongly solvated, leading to the stabilization of solvent-separated ion pairs over contact ion pairs. Minimum energy pathway analysis revealed that coordinating two sodium ions with a sulfate ion is a multistep process whereby the sodium ions coordinate to the sulfate ion sequentially. Finally, we demonstrate that MLIPs allow the study of solvated ions beyond simple monovalent pairs with DFT-level accuracy in their low concentration limit (0.1 M) via statistically converged properties from ns-long simulations.

anions↗
DOE OSTI2025

Temperature-Dependent Speciation of Ni(II) in Molten Chloride, Bromide, and Iodide Salts

Understanding the fundamental speciation and local structure of metal cations─such as nickel, which is present as a corrosion product─in molten salt media is essential for assessing their impact on the thermal and physical properties of salts used for advanced nuclear energy applications. In this study, we employed an integrated approach of combining X-ray absorption spectroscopy (XAS), UV–vis spectroscopy, and ab initio molecular dynamics (AIMD) to investigate the halide-dependent speciation of divalent nickel ions, Ni(II), in LiCl–KCl, LiBr–KBr, and LiI–KI eutectic salt mixtures from room temperature to 600 °C. Our findings show that both halide type and temperature significantly influence the local coordination of nickel in these molten salt systems. A notable aspect of Ni(II) speciation in LiI–KI salts is the formation of a polyiodide species. In conclusion, our findings provide critical insights into metal ion speciation in molten salt systems, enhancing our understanding of the factors that govern the physicochemical properties of these complex fluids.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Phase Transition and Liquid-like Superionic Conduction in Ag 2 S

Superionic conductivity in certain polymorphs of Ag 2 S has inspired numerous concepts for materials applications, but the relationship between the structure and the mobility of silver ions remains poorly explored. In this paper, we report ab initio molecular dynamics simulations for low- (acanthite) and high-temperature (argentite) Ag 2 S polymorphs that reveal the dynamical processes, giving rise to the superionic behavior in the latter. Similarities between their sulfur sublattices enable simulations of silver ion diffusivities and pathways on essentially an equal footing. For the higher temperature polymorph, calculated temperature-dependent mean square displacements and activation energies by the nudged elastic band method show good correspondence with expectations from the experiment. In the superionic state, silver atoms diffuse in a liquid-like behavior with no preferred diffusion pathways, within the relatively stable body-centered cubic sulfur framework. Conversely, conduction in acanthite appears to depend more on the mobilities of electronic charge carriers.

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

Surface Hydration and Hydroxyl Configurations of Gibbsite and Boehmite Nanoplates

Discontinuation of the bulk structure at the interface between metal oxide particles and water leads to altered bonding characteristics and unique facet-dependent molecular environments. Surface hydration and consequent hydroxylation adds further complexity to the interface, details that for metal (oxy)hydroxides are especially difficult to isolate from the background signal of bulk structural hydroxyls. Here we probe for the first time the hydroxylation and hydration structures on basal surfaces of gibbsite (?-Al(OH)3) and boehmite (?-AlOOH) nanoplates at ambient conditions, using the interface-sensitive technique vibrational sum frequency generation spectroscopy (VSFG). Hydroxyl stretching modes at the interfaces with adsorbed water layers compared directly to Raman and infrared bulk modes show that while gibbsite surface frequencies were sharp and nearly identical to those in the bulk, boehmite surface hydroxyls displayed a very different broad spectrum of states. Ab initio molecular dynamics simulations of both basal surfaces with and without hydration waters reveal that gibbsite surface hydroxyls interact only weakly with overlying hydration waters remaining essentially unperturbed, whereasthose on boehmite can interact more strongly facilitated by higher configurational degrees of freedom at the interface. The findings clearly unveil substantial differences in the hydrated interfacial dynamics of these two otherwise similar materials, with implications for their interfacial chemistry and wettability.

Sassi, Michel JPC↗
DOE OSTI2020

Consequences of 131 I Transmutation in Gas Phase Radioiodine Molecules and Adsorbed on Graphite Surface

In the case of nuclear power plant accidents, several radionuclides are being dispersed into the environment. In particular, the release of radioiodine 131 I represents a major health concern because of its absorption by the thyroid gland. To reduce the release of 131 I into the environment, capture materials, such as activated carbon, provide a mean to efficiently retain radioiodine-bound molecular species. Here, we use ab initio molecular dynamic calculations to investigate the effects of 131 I radioactive decay on the molecular stability of I 2 , acid iodine (HI, HIO, HIO 2 , and HIO 3 ), and methyl-iodine species in both gas phase and adsorbed on graphite. In gas phase, we found that I 2 is the least resistant to 131 I transmutation, followed by the HIO x family of molecules, then HI and CH 3 I which are the most resistant. However, the graphite surface is found to greatly favorize the fragmentation of adsorbed molecular species. While the identification of gas phase fragments can offer new chemical route to consider in chemical kinetics models, the formation reactive fragment products at the graphite surface can favorize their reaction with surface carbon atoms which overtime could reduce the number of adsorption sites of capture materials available, and hence, impacting their capture capability.

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

Intercalation of Lithium into Graphite: Insights from First-Principles Simulations

Understanding ion intercalation at electrode–electrolyte interfaces is key to the development of energy storage and water desalination. Here, we investigate Li + kinetics at a prototypical interface between graphite anodes and an organic electrolyte, and we elucidate key factors that determine ion transport, using first-principles methodology coupling ab initio molecular dynamics simulations with a solvation model. We show that surface chemical composition significantly influences the kinetics of ion intercalation from the liquid into graphite. We find that this is partly related to the ion desolvation process, which varies notably for different graphite surface chemical terminations. In addition, interfacial polarization is found to play an important role in determining energy barriers for ion transfer. We also discuss the impact of electrode potentials, which is often neglected in conventional first-principles calculations despite being a key factor in device configurations. Our study provides insights into the coupling of electronic and ionic effects of interfacial chemistry on ion transport at complex electrode–electrolyte interfaces.

25 ENERGY STORAGE↗
DOE OSTI2020

Electron Trapping and Ion Leaching at the Li-Modified Quartz–Water Interface

Electrons and ions at metal oxide-water interfaces have a critical role in many phenomena and applications, thus making their properties of considerable interest. We here investigate the behavior of a lithium impurity at the water interface with Li-doped SiO 2 using hybrid density functional based ab initio molecular dynamics simulations. We find that the excess electron donated by the lithium dopant localizes on a surface silicon atom and is partially solvated at the aqueous interface. While the position of the excess electron does not change during our ~60 ps simulation, the lithium cation diffuses from its initial position in the SiO 2 subsurface toward the interface and eventually leaches out of the oxide surface and becomes solvated by four water molecules, forming an aqueous Li+- electron complex at the interface. The degree of interaction between the localized electron and the lithium ion controls not only the energy level of the excess electron but also the extent of structural distortion and hydrogen-bonding at the interface. The results of our study can be relevant for understanding the initial stages of glass corrosion, where alkali ions are leached from the surface region by interactions with water.

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

Mechanistic Study of the Li–Air Battery with a Co 3 O 4 Cathode and Dimethyl Sulfoxide Electrolyte

The lithium–air battery, a powerful competitor to replace the traditional lithium-ion battery, has attracted increasing attention due to its extremely high theoretical energy density. However, its development is limited by the cathode and electrolyte properties, which should include high stability, conductivity, and electrocatalytic properties in oxygen-rich environments. Here, we employ a systematic first-principles study of Li–O 2 discharge and charge reactions on the Co 3 O 4 -based cathode with the assistance of dimethyl sulfoxide (DMSO) electrolyte. The structure, stability, and electronic properties of different surface reconstructions of the Co 3 O 4 (100) facet are investigated. In addition, the mechanisms and thermodynamic overpotentials of multi-step reactions between Li + /e – and O 2 are provided, where lithium suboxide products (Li 2 O 2 or Li 3 O 2 ) are formed on the different Co 3 O 4 (100) terminations. The solvation shell of Li + components in explicit DMSO solvent is investigated through ab initio molecular dynamics simulations. In general, we find that the Co 3 O 4 (100)-O (oxidized) surface is the most stable one under standard conditions, and the stable Li + solvation structure is found in a tetrahedral Li(DMSO) 4 + shell in the DMSO-based electrolyte. Moreover, in the system of the Co 3 O 4 (100)-O cathode and DMSO electrolyte, the solution model pathway is energetically favorable for the Li–O 2 discharge reaction. It provides a low constant overpotential of 0.17 V during a long-term discharging process, thus causing the final toroid Li 2 O 2 formation on the cathode. During the charging process, an overpotential of 0.36 V is required to rapidly decompose Li 2 O 2 .

25 ENERGY STORAGE↗
DOE OSTI2023

Cation-Induced Disruption of the Local Structure of Water in Faujasite

Here, the changes in the structure of water confined in faujasite were resolved by a series of ab initio molecular dynamic simulations of explicit water molecules in faujasite with varying Si/Al ratio and counterions present. Radial distribution functions of water oxygen atom pairs, as well as counterions to Al pairs, were calculated to quantify the structure of water in faujasite and the distribution of counterions in the water-filled faujasite pores. Radial distribution functions reveal that confinement without counterions results in more water molecules hydrogen bound to each other. In hydrogen-exchanged faujasite, the increase in aluminum content led to the disruption of the confined water structure as charge-compensating protons dissociated from Al tetrahedra and entered the liquid water in the pores. In addition, a potential of mean force of proton transfer shows that protons are thermodynamically favored to be in the liquid water rather than bound on Al tetrahedra. Exchanged Na + and La 3+ counterions, which remain closer to Al tetrahedra in liquid water than H + , also disordered the confined water structure. Two La 3+ ions disrupted water less than did six Na + ions, suggesting that, for the alkali and lanthanide ions, the number of counterions in FAU has a larger effect on the structure of water than the charge of the ion itself.

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

Proton Dissociation and Delocalization under Stepwise Hydration of Zeolite HZSM-5

The protonation behavior of zeolite Brønsted acid sites (BAS) in the presence of water is important for the performance of these widely used catalysts. Despite extensive study, the number of water molecules necessary for deprotonation is not well understood, in large part because experiments have been unable to access this information. In this work, we report experimental evidence for full deprotonation of the BAS in the presence of two or more water molecules, with a deprotonation energy of 1.6 kcal/mol. Linear IR absorption and 2D IR spectra were measured over a wide range of controlled hydration levels from 0.5 to 8.0 equivalents of H 2 O/Al at a constant temperature. Distinct spectral signatures of the protonated BAS and excess proton are identified, and their hydration dependence is analyzed quantitatively. Using the experiment as a benchmark, ab initio molecular dynamics simulations are reported that reproduce the experimental trends in the protonation state and IR spectra. The proton charge position and delocalization are quantified in clusters of 1–8 H 2 O molecules using the recently developed rCEC method. This analysis provides insight into the proton structure in confined water clusters, showing that the excess charge remains relatively localized between two oxygen atoms across the hydration range.

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

Confinement-Controlled Rearrangements in Dioxolane Upgrading on H-ZSM-5 Revealed by Periodic DFT

Zeolitic Brønsted acid sites catalyze carbocation rearrangements central to upgrading biomass-derived oxygenates. Here we elucidate the mechanism of dioxolane conversion to methyl ethyl ketone and isobutanal on H-ZSM-5 using periodic density functional theory on the MFI model, complemented by ab initio molecular dynamics to probe confinement effects. Dioxolane adsorption at the Brønsted site is followed by protonation-assisted ring opening to form an oxocarbenium intermediate stabilized by the deprotonated framework. From this common intermediate, selectivity is governed by two competing rearrangements, namely, the 1,2-hydride shift with a free-energy barrier of 18.05 kcal mol –1 at 498 K leading toward MEK, and the 1,2-methyl shift with a higher barrier of 25.40 kcal mol –1 leading toward isobutanal. The hydride-shift channel is kinetically preferred over the methyl-shift channel, lowering the isobutanal/MEK ratio below the 3:1 limit expected for equal branching. Adsorption thermodynamics further indicate stronger stabilization of MEK than isobutanal within ZSM-5 channels, suggesting that confinement-controlled binding can bias product distributions in addition to intrinsic rearrangement barriers. These results highlight how Brønsted acidity and pore confinement jointly shape the rearrangement landscape in MFI zeolites.

adsorption↗
DOE OSTI2021

Resolving Configurational Disorder for Impurities in a Low-Entropy Phase

Hematite (α-Fe 2 O 3 ) exerts a strong control over the transport of minor but critical metals in the environment and is used in multiple industrial applications; the photocatalysis community has explored the properties of hematite nanoparticles over a wide range of transition metal dopants. Nonetheless, simplistic assumptions are used to rationalize the local coordination environment of impurities in hematite. Here, we use ab initio molecular dynamics (AIMD)-guided structural analysis to model the extended X-ray absorption fine structure (EXAFS) of Cu 2+ - and Zn 2+ -doped hematite nanoparticles. Specific defect–impurity associations were identified, and the local coordination environments of Cu and Zn both displayed considerable configurational disorder that, in aggregate, approached Jahn–Teller-like distortion for Cu but, in contrast, maintained hematite-like symmetry for Zn. This study highlights the role of defects in accommodating impurities in a nominally low-entropy phase and the limits to traditional shell-by-shell fitting of EXAFS for dopants/impurities in unprecedented bonding environments.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Hydrogen Bonds and H 3 O + Formation at the Water Interface with Formic Acid Covered Anatase TiO 2

Carboxylic acid-modified TiO 2 surfaces in aqueous environment are of widespread interest, yet atomic-scale understanding of their structure is limited. In this work, we investigate formic acid (FA) on anatase TiO 2 (101) (A-101) in contact with water using density functional theory (DFT) and ab-initio molecular dynamics (AIMD). Isolated FA molecules adsorbed in a deprotonated bridging bidentate (BD) form on A-101 are found to remain stable at the interface with water, with the acid proton transferred to a surface oxygen to form a surface bridging hydroxyl (O br H). With increasing FA coverage, adsorbed monolayers of only BD and successively of alternating monodentate (MD) and BD species give rise to a higher concentration of surface O br H’s. Simulations of these adsorbed monolayers in water environment show that some protons are released from the surface O br H’s to water resulting in a negatively charged surface with nearby solvated H 3 O + ions. These results provide insight into the complex acid-base equilibrium between an oxide surface, adsorbates and water and can also help obtain a better understanding of the wetting properties of chemically modified TiO 2 surfaces.

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

Spectral Signatures of Ultrafast Excited-State Intramolecular Proton Transfer from Computational Multi-edge Transient X-ray Absorption Spectroscopy

Excited state intramolecular proton transfer (ESIPT) is a fundamental chemical process with several applications. Ultrafast ESIPT involves coupled electronic and atomic motions and has been primarily studied using femtosecond optical spectroscopy. X-ray spectroscopy is particularly useful because it is element-specic and enables direct, individual probes of the proton donating and accepting atoms. Here, we report a computational study to resolve the ESIPT in 10-Hydroxybenzo[h]quinoline (HBQ), an intramolecularly hydrogen bonded compound, using time-dependent density functional theory combined with ab initio molecular dynamics (AIMD) and time-resolved X-ray absorption spectroscopy (XAS) computations to track the ultrafast excited-state dynamics. Overall, our results reveal clear X-ray spectral signatures of coupled electronic and atomic motions during and following ESIPT at the oxygen and nitrogen K-edge paving the way for future experiments at X-ray free electron lasers.

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

Birth of the Hydrated Electron via Charge-Transfer-to-Solvent Excitation of Aqueous Iodide

A primary means to generate hydrated electrons in laboratory experiments is excitation to the charge-transfer-to-solvent (CTTS) state of a solute such as I – (aq), but this initial step in the genesis of e – (aq) has never been simulated directly using ab initio molecular dynamics. Here, we report the first such simulations, combining ground- and excited-state simulations of I – (aq) with a detailed analysis of fluctuations in the Coulomb potential experienced by the nascent solvated electron. What emerges is a two-step picture of the evolution of e – (aq) starting from the CTTS state: I – (aq) + hν → I –* (aq) → I • (aq) + e – (aq). Notably, the equilibrated ground state of e – (aq) evolves from I –* (aq) without any nonadiabatic transitions, simply as a result of solvent reorganization. The methodology used here should be applicable to other photochemical electron transfer processes in solution, an important class of problems directly relevant to photocatalysis and energy transfer.

Carter-Fenk, Kevin↗