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 343 records · Page 19

Li + Diffusion in Amorphous and Crystalline Al 2 O 3 for Battery Electrode Coatings

Al 2 O 3 is often applied protectively to lithium-ion battery anode and cathode materials to inhibit surface degradation, suppress dendrite formation, and relieve mechanical stresses. Given the very high intrinsic band gap and diffusion barrier of the material, the mechanism that allows Li diffusion through these coatings is not well understood, and widely varying laboratory results indicate that there may be dependencies on morphology and stoichiometry. Using nudged elastic band calculations and ab initio molecular dynamics, we perform a systematic investigation across Al 2 O 3 structures, both crystalline and amorphous, and at various concentrations of Li + to uncover the optimal parameters for maximally diffusive coatings. We find a correlation between the low proximity of Li + to Al 3+ and the low Li + migration barrier. Although barriers are the lowest in the highly diffusive one-dimensional channels of crystalline θ-Al 2 O 3 , the system is structurally delicate and subject to detrimental distortion as the Li + content is increased. The α-Al 2 O 3 lattice is, conversely, highly stable against distortion at all Li + concentrations but disadvantageous for Li + migration. In amorphous systems, unscreened Li + –Li + Coulomb repulsion and pre-emptive occupation of “trapping sites” combine to lower the energy barriers as a function of increasing concentration. One of our most important findings is that Al-deficient materials can sharply increase Li + movement, and we predict that an amorphous material with a combination of high Li + concentration and Al deficiency would enable highly Li + -conductive protective coatings for electrodes.

25 ENERGY STORAGE↗

Tailoring Interfaces in Solid-State Batteries Using Interfacial Thermochemistry and Band Alignment

Solid-state lithium ion batteries have enhanced thermal stability compared to batteries using liquid electrolytes. Although their practical implementation is limited by undesired side reactions between the electrode and electrolyte, different modification strategies have been proposed to stabilize these reactive interfaces. These approaches have been primarily based upon bulk materials properties, however, which may not necessarily be representative of the chemistry at the solid/solid interface. In this work, first principles calculations and X-ray scattering experiments are used to elucidate molecular-level reactivity and develop design principles for tailored interfaces based on surface and interfacial properties. Because of its well-known instability, the interface between the Li metal anode and the lithium lanthanum titanate (LLTO) solid electrolyte is applied as a model system. Density functional theory calculations are used to describe bulk, surface, and interfacial thermochemistry of the Li/LLTO system, and interfacial reconstruction is probed using ab initio molecular dynamics. These simulations of the Li/LLTO interface demonstrate facile decomposition concomitant with reduction of Ti 4+ ions. Based on further insights from computational analysis of the surface band edge positions, La 2 O 3 is proposed as an interlayer coating and is shown to provide an energetic barrier for interfacial charge transfer and decomposition reactions from X-ray reflectivity measurements and theoretical calculations. The findings in this work suggest that design strategies based on surface and interfacial properties can be used to kinetically stabilize electrode/electrolyte interfaces in solid-state batteries.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Solid-State Calcium-Ion Diffusion in Ca 1.5 Ba 0.5 Si 5 O 3 N 6

Rechargeable batteries based on multivalent working ions are promising candidates for next-generation high-energy-density batteries. Development of these technologies, however, is largely limited by the low diffusion rate of multivalent ions in solid-state materials, thereby necessitating a better understanding of the design principles that control multivalent-ion mobility. We report Ca 1.5 Ba 0.5 Si 5 O 3 N 6 as a potential calcium solid-state conductor and investigate its Ca migration mechanism by means of ab initio computations and neutron diffraction. This compound contains partially occupied Ca sites in close proximity to each other, providing a unique mechanism for Ca migration. Nuclear density maps obtained with the maximum entropy method from neutron powder diffraction data provide strong evidence for low-energy percolating one-dimensional pathways for Ca-ion migration. Ab initio molecular dynamics simulations further support a low Ca-ion migration barrier of ~400 meV when Ca vacancies are present and reveal a unique "vacancy-adjacent"concerted ion migration mechanism. This work provides a new understanding of solid-state Ca-ion diffusion and insights into the future design of novel cation configurations that utilize the interactions between mobile ions to enable fast multivalent-ion conduction in solid-state materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Functionalization of Electrodes with Tunable [EMIM] x [Cl] x +1 – Ionic Liquid Clusters for Electrochemical Separations

Functionalization of electrodes with clusters of hydrophobic molecules may improve the energy efficiency and selectivity of electrochemical separations by modulating the desolvation process occurring at the interface. Ionic liquids (IL), which have a broad range of compositions and properties, are potential candidates for controlling the transport, desolvation, and adsorption of target ions at electrochemical interfaces. We report a joint experimental and theoretical investigation of the structure, stability, and selective adsorption properties of the IL clusters 1-ethyl-3-methylimidazolium chloride [EMIM] x [Cl] x+1 - (x = 1 – 10) and demonstrate their ability to adsorb and separate ions from solution. The structure and stability of the IL clusters are determined experimentally using high-mass-resolution electrospray ionization mass spectrometry, collision-induced dissociation, and negative ion photoelectron spectroscopy. Global optimization theory and ab initio molecular dynamics simulations provide molecular-level insight into the bonding and structural fluxionality of these species. Ion soft landing is used to selectively functionalize the surface of highly oriented pyrolytic graphite (HOPG) working electrodes with [EMIM] 1 [Cl] 2 - , [EMIM] 3 [Cl] 4 - , and [EMIM] 5 [Cl] 6 - clusters. Kelvin probe microscopy provides insight into the relative stability of the clusters on HOPG and their effect on the work function of IL-functionalized electrodes. Cyclic voltammetry measurements reveal irreversible adsorption of Fe(CN) 6 4-/3- anions during redox cycling, while electrochemical impedance spectroscopy indicates a substantial decrease in the electron transfer resistance of the IL-functionalized electrodes due to adsorption of Fe(CN) 6 4-/3- . Overall, our findings demonstrate that IL clusters with different size and stoichiometry may be used to increase the efficiency of electrochemical separations, opening new horizons in selective electrode functionalization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Medium-Range Ordering in the Ionic Glass Electrolytes LiPON and LiSiPON

In this work, we provide an in-depth structural characterization of the amorphous ionic glasses LiPON and LiSiPON with high Li content. Based on ab initio molecular dynamics simulations, the structure of these materials is an inverted structure with either isolated polyanion tetrahedra or polyanion dimers in a Li+ matrix. Based on neutron scattering data, this type of inverted structure leads to a significant amount of medium-range ordering in the structure, as demonstrated by two sharp diffraction peaks and a periodic structural oscillation in the density function G(r). While this medium-range ordering is commonly observed in liquids and metallic glasses, it has not previously been observed in oxides. On a local scale, adding N and Si increases the number of anion bridges and polyanion dimer structures, leading to higher ionic conductivity. In the medium-range ordering, the addition of Si leads to more disorder in the polyanion substructure but a significant increase in the ordering of the O substructure. Finally, we demonstrate that this inverted structure with medium-range ordering results in a glassy material that is both mechanically stiff and ductile on the nanoscale.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Semiconducting Zn x Mo 3 S 13 -GO Chalcocarbogel: A High-Capacity and Stable Sulfur-Equivalent Conversion-Based Electrode for Lithium-Ion Batteries

Lithium–sulfur batteries with a sulfur electrode offer a theoretical capacity of ∼1672 mAh g –1 , but rapid capacity loss mainly constrains their practical application. This work introduces a semiconducting and amorphous Zn x Mo 3 S 13 -GO (x = 0.5) chalcocarbogel sulfur-equivalent electrode with superior capacity and stability for lithium-ion batteries (LIBs). The Zn x Mo 3 S 13 -GO is synthesized in solution under ambient conditions, and its local structure contains S–S, M-Q (M = Mo, Zn; Q = S, O), C–S, and Mo–Mo bonding motifs with Mo coordination environment closely related to Mo 3 S 13 anions, as determined by X-ray photoelectron spectroscopy, synchrotron X-ray scattering, X-ray absorption spectroscopy, and ab initio molecular dynamics simulations. The Li/Zn x Mo 3 S 13 -GO cell offers an initial discharge capacity of 1019 mAh g –1 at a rate of C/3. After the activation cycles, the Li/Zn x Mo 3 S 13 -GO cell demonstrates good cycling stability, retaining a discharge capacity of 519.4 mAh g –1 after 250 cycles with ∼99.98% Coulombic efficiency and excellent rate capabilities. Moreover, it provides an initial discharge capacity of ∼574 mAh g –1 and maintains a retention capacity of 279 mAh g –1 at 1C after 625 cycles. The Lewis acidic Zn 2+ ion enhances the Lewis basic polysulfide anchoring ability and reduces the dissolution of polysulfides produced during the redox process through Zn–S covalent interaction, while the semiconducting and amorphous structure of the chalcocarbogel increases the electrical and ionic conductivity. Furthermore, this work highlights chalcocarbogels’ potential for developing high-capacity and stable electrodes for LIBs.

25 ENERGY STORAGE↗

Thermodynamic and Kinetic Mechanisms Governing the Synthesis of Nickel-Poor Cathodes

A deeper understanding of the thermodynamics and kinetics governing the lithiation and layering mechanisms of NMC cathode materials (LiNixMnyCozO 2 , where x + y + z = 1) offers valuable insights for enhancing synthesis methods and improving cathode performance. By employing atomistic and mesoscale approaches informed by in situ powder X-ray diffraction (PXRD) experiments, critical parameters for comprehending lithiation and layering processes and reaction rates were identified. The mesoscale approach captured the evolution of the phases and crystallite size observed in the in situ PXRD, revealing the differences in reaction rates with the use of different lithium salts and starting precursors. Ab initio molecular dynamics (AIMD) underscored the importance of vacancies and structural defects in promoting ion mobility and facilitating the nucleation of a layered domain. This nucleation disrupts the symmetry of disordered phases, ultimately creating a strained phase that serves as a buffer between layered and disordered regions. The lithiation and layering processes reflect a dynamic balance between the thermodynamic drive for a low-energy layered structure and the kinetic of diffusion, which is influenced by temperature and lithium vacancy concentration. Overall, reaction mechanisms are driven by the inherent defects of the intermediate phase that differ for NMC cathode materials. The lithium salts impact the rates of lithiation and layering, with a much slower process for Li 2 CO 3 .

25 ENERGY STORAGE↗

Understanding Superionic Conductivity in Lithium and Sodium Salts of Weakly Coordinating Closo -Hexahalocarbaborate Anions

Solid-state ion conductors based on closo-polyborate anions combine high ionic conductivity with a rich array of tunable properties. Cation mobility in these systems is intimately related to the strength of the interaction with the neighboring anionic network and the energy for reorganizing the coordination polyhedra. In this paper, we explore such factors in solid electrolytes with two anions of the weakest coordinating ability, [HCB 11 H 5 Cl 6 ] – and [HCB 11 H 5 Br 6 ] – , and a total of 11 polymorphs are identified for their lithium and sodium salts. Our approach combines ab initio molecular dynamics, synchrotron X-ray powder diffraction, differential scanning calorimetry, and AC impedance measurements to investigate their structures, phase-transition behavior, anion orientational mobilities, and ionic conductivities. We find that M(HCB 11 H 5 X 6 ) (M = Li, Na, X = Cl, Br) compounds exhibit order–disorder polymorphic transitions between 203 and 305 °C and display Li and Na superionic conductivity in the disordered state. Through detailed analysis, we illustrate how cation disordering in these compounds originates from a competitive interplay among the lattice symmetry, the anion reorientational mobility, the geometric and electronic asymmetry of the anion, and the polarizability of the halogen atoms. These factors are compared to other closo-polyborate-based ion conductors to suggest guidelines for optimizing the cation–anion interaction for fast ion mobility. This study expands the known solid-state poly(carba)borate-based materials capable of liquid-like ionic conductivities, unravels the mechanisms responsible for fast ion transport, and provides insights into the development of practical superionic solid electrolytes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Influence of Cation Size on the Local Atomic Structure and Electronic Properties of Ta Perovskite Oxynitrides

Partial anion substitution in transition metal oxides provides rich opportunities to control and tune physical and chemical properties, for example, combining the merits of oxides and nitrides. In addition, the possibility of resulting anion sublattice order provides a means to target polar and chiral structures based on a wide array of accessible structural archetypes by design. Here, in this work, we investigate the local structures of a family of perovskite tantalum oxynitrides—ATaO 2 N (A = Ba, Sr, and Ca)—using a combination of experimental and theoretical approaches including neutron total scattering, density functional theory (DFT), and ab initio molecular dynamics (AIMD) simulations. We present the first experimental study of chemical short-range order (CSRO) in CaTaO 2 N, confirming local cis N ordering of the anion sub-lattice. Our systematic exploration of a local structure across the A cation size series (from the larger Ba to the smaller Ca) reveals a perovskite motif increasingly distorted with respect to long-range average structures. DFT and AIMD simulations support the observed trends. Overall, structures with cis ordering of the nitrogen anions in each TaO 4 N 2 octahedron are favored over those with trans ordering. With diminishing A cation size, local cis ordering and Ta off-centering play decreasing roles in overall lattice stability, overshadowed by the stabilizing effects of octahedral tilting. The influence of these factors on local dipole formation and frustrated dipole ordering are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Solvation Structure of 237 Np 4+ in a Noncomplexing Environment

Here, the solvation structure of an Np 4+ ion in an aqueous, noncomplexing and nonoxidizing environment of trifluoromethanesulfonic (triflic) acid was investigated with X-ray absorption spectroscopy (XAS) combined with ab initio molecular dynamics (AIMD) and time-dependent density functional theory (TDDFT) calculations. Np L III -edge X-ray absorption near-edge structure (XANES) and extended X-ray absorption fine structure (EXAFS) data were collected for Np 4+ in 1, 3, and 7 M triflic acid using a laboratory-scale spectrometer and separately at a synchrotron facility, producing data sets in excellent agreement. TDDFT calculations revealed a weak pre-edge feature not previously reported for Np L III -edge XANES. AIMD modeling results showed differences in the hydration shell of the Np 4+ ion at different concentrations of triflic acid; these results are supported by the experiment. EXAFS fit models to the experiment resulted in similar coordination of Np 4+ in noncomplexing aqueous media as reported in the literature for 1 M perchloric acid but, together with calculations, revealed more than one distance between Np and O atoms in 7 M triflic acid. These results imply monodentate coordination with sulfonate groups in 7 M triflic acid and suggest the possibility of proto-neptunyl species in relatively low-concentration Np 4+ acid solutions.

Boglaienko, Daria V. [Pacific Northwest National L↗

Modeling Molecules under Pressure with Gaussian Potentials

The computational modeling of molecules under high pressure is a growing research area that augments experimental high-pressure chemistry. Here, a new electronic structure method for modeling atoms and molecules under pressure, Gaussians On Surface Tesserae Simulate HYdrostatic Pressure (GOSTSHYP) approach, is introduced. In this method, a set of Gaussian potentials is distributed evenly on the van der Waals surface of the investigated chemical system, leading to a compression of the electron density and the atomic scaffold. Since no parameters other than pressure need to be specified, GOSTSHYP allows straightforward geometry optimizations and ab initio molecular dynamics simulations of chemical systems under pressure for nonexpert users. In this work, calculated energies, bond lengths, and dipole moments under pressure fall within the range of established computational methods for high-pressure chemistry. A Diels-Alder reaction and the cyclotrimerization of acetylene showcase the ability of GOSTSHYP to model pressure-induced chemical reactions. The connection to mechanochemistry is pointed out.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Assessing the Limitations of Self-Interaction-Corrected Functionals for Describing the Hydrated Electron

Simulating the hydrated electron using density functional theory is challenging due to the prevalence of self-interaction error in standard functionals. Hybrid functionals like PBEh(40) can reasonably describe the chemistry of an excess electron in water and partially mitigate self-interaction error by incorporating exact Hartree–Fock exchange, but they are computationally expensive making them impractical for large-scale and long-time ab initio molecular dynamics simulations. Explicit self-interaction correction schemes that are applied on an orbital-by-orbital basis offer a potential alternative when the correction is limited to the singly occupied molecular orbital obtained with a generalized gradient approximation functional. Here, we examine whether the Perdew–Zunger self-interaction correction scheme applied to the revPBE functional can provide a computationally efficient and physically sensible alternative to PBEh(40) for the hydrated electron. We find that functionals incorporating a self-interaction correction scheme should be viewed with caution when applied to the hydrated electron and its reactivity. Furthermore, we show that it is critical to consider extensive sampling and diverse chemical environments when validating their performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The CP-PAW Code Package for First-Principles Calculations from a User’s Perspective

CP-PAW is a combined electronic structure and ab initio molecular dynamics code to perform mixed quantum and classical simulations of atomistic condensed phase systems, such as solids, liquids, and molecular systems. As the name suggests, the CP-PAW code unifies the all-electron projector augmented-wave (PAW) method with the Car–Parrinello (CP) approach to determine not only the electronic and nuclear ground states of condensed matter but also to study their properties and dynamics. In addition to briefly outlining the underlying theory, the focus will be on the unique aspects of CP-PAW and how to correctly employ them as a user. How to install CP-PAW using the new build system will also be briefly mentioned.

Blöchl, Peter E [Institute for Theoretical Physic↗

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↗

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↗

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↗

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↗

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↗