Search NASA⌕ Search

SEARCH · Search NASA

Results for “ab 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 91 records · Page 5

Nanoscale Editing of Multi and Single Layer Tungsten Disulfide via Gas‐Assisted Focused Electron Beam Induced Etching for Device Prototyping

Focused electron beam induced etching (FEBIE) with XeF 2 (xenon difluoride) precursor is conducted on multi-layer exfoliated WS 2 (tungsten disulfide) and monolayer WS 2 grown by chemical vapor deposition (CVD). The films are characterized by atomic force microscopy (AFM) and Raman and photoluminescence (PL) spectroscopy post-etching. The etch rates/efficiencies are reported as a function of electron beam energy, current, dwell time, and XeF 2 pressure. Bulk film Raman spectra are unchanged post-FEBIE, indicating minimal subsurface damage. Monolayer WS 2 shows a decrease in Raman and PL intensity post-FEBIE, with a dose-to-clear of ≈2 nC µm −2 . The study reveals regimes affected by the various mass transport contributions such as refresh time and the ratio of electrons/XeF 2 . Spontaneous etching was discovered during FEBIE of large patterned areas due to the long frame/refresh times. Density functional theory and ab initio molecular dynamics simulations compares desorption of SF x and WF x molecules from pristine WS 2 basal planes and pore edges, revealing the spontaneous etching is consistent with etching of partially etched monolayers during each frame. Single-line etching width of 21 nm, and patterning flakes into 100 nm wide channels are demonstrated. In conclusion, this work demonstrates the possibility of editing WS 2 flakes into electronic devices of arbitrary dimensions for semiconductor applications.

2D materials↗

Breakthrough Conductivity Enhancement in Deep Eutectic Solvents via Grotthuss–Type Proton Transport

There is an increasing demand for the development of ion-conducting electrolytes for energy storage systems. Much attention is directed toward deep eutectic solvents as potential candidates. In the search for highly conductive systems, the possibility of designing deep eutectic solvents with Grotthuss-type proton transport is widely overlooked. Herein, ethaline, a mixture of choline chloride and ethylene glycol is used in a 1:2 molar ratio, to induce a significant conductivity increase with the addition of water and sulfuric acid (H 2 SO 4 ). The achieved breakthrough conductivity is analyzed experimentally and simulated with ab initio molecular dynamics (AIMD). At sufficient water content, an H-bonding network is formed that leads to a significant breakthrough conductivity based on H 2 SO 4 -derived proton transfer following the long-established Grotthuss proton transport mechanism. This result is substantiated by the positive deviation from the ideal KCl line in the Walden plot. Specifically, the data series positioned above the reference line indicates a Grotthuss mechanism in action. The AIMD simulations demonstrate proton transfer between water and ethylene glycol, supported by simulation frames captured at various times.

36 MATERIALS SCIENCE↗

AIMD‐Based Protocols for Modeling Exciplex Fluorescence Spectra and Inter‐System Crossing in Photocatalytic Chromophores

ABSTRACT This study introduces a computational protocol for modeling the emission spectra of exciplexes using excited‐state ab initio molecular dynamics (AIMD) simulations. The protocol is applied to a model exciplex formed by oligo‐p‐phenylenes (OPPs) and triethylamine (TEA), which is of interest in the context of photocatalytic reduction of . AIMD facilitates efficient sampling of the conformational space of OPP3 and OPP4 exciplexes with TEA, offering a dynamic alternative to previously employed static methods. The AIMD‐based protocol successfully reproduces experimental emission spectra for OPP‐TEA exciplexes, agreeing with previous computational and experimental findings. The results show that AIMD simulations provide an efficient means of sampling the conformational space of these exciplexes, requiring less user input and, in some instances, fewer computational resources than multiple excited‐state optimizations initiated from user‐specified initial structures. The study also evaluates the yield of intersystem crossing (ISC) using AIMD and Landau‐Zener probability. The results suggest that ISC is a minor decay channel for OPP3 and OPP4. This work provides new insights into the structural flexibility and emission characteristics of OPP‐TEA photoredox catalyst systems, potentially contributing to improved design strategies for organic chromophores in reduction applications.

Giudetti, Goran [Department of Chemistry Universit↗

Interactive multiscale modeling to bridge atomic properties and electrochemical performance in Li-CO 2 battery design

Li-CO 2 batteries are promising energy storage systems due to their high theoretical energy density and CO 2 fixation capability, relying on reversible Li 2 CO 3 /C formation during discharge/charge cycles. Here, we present a multiscale modeling framework integrating Density Functional Theory (DFT), Ab-Initio Molecular Dynamics (AIMD), classical Molecular Dynamics (MD), and Finite Element Analysis (FEA) to investigate atomic and cell-level properties. The considered Li-CO 2 battery consists of a lithium metal anode, an ionic liquid electrolyte, and a carbon cloth cathode with Sb 0.67 Bi 1.33 Te 3 catalyst. DFT and AIMD determined the electrical conductivities of Sb 0.67 Bi 1.33 Te 3 and Li 2 CO 3 using the Kubo–Greenwood formalism and studied the CO 2 reduction mechanism on the cathode catalyst. MD simulations calculated the CO 2 diffusion coefficient, Li + transference number, ionic conductivity, and Li + solvation structure. The FEA model, parameterized with atomistic simulation data, reproduced the available experimental voltage–capacity profile at 1 mA/cm 2 and revealed spatio-temporal variations in Li 2 CO 3 /C deposition, porosity, and CO 2 concentration dependence on discharge rates in the cathode. Accordingly, Li 2 CO 3 can form large and thin film deposits, leading to dispersed and local porosity changes at 0.1 mA/cm 2 and 1 mA/cm 2 , respectively. The capacity decreases exponentially from 81,570 mAh/g at 0.1 mA/cm 2 to 6200 mAh/g at 1 mA/cm 2 , due to pore clogging from excessive discharge product deposition that limits CO 2 transport to the cathode interior. Therefore, the performance of Li-CO 2 batteries can be improved by enhancing CO 2 transport, regulating Li 2 CO 3 deposition, and optimizing cathode architecture.

Battery performance↗

Strain-modulated intercalated phases of Pb monolayer with dual periodicity in SiC(0001)-graphene interface

Intercalation of metal atoms at the SiC(0001)-graphene (Gr) interface can provide confined 2D metal layers with interesting electronic properties. The intercalated Pb monolayer (ML) has shown the coexistence of the Gr(10 x 10)-moiré and a stripe phase, which still lacks understanding. Using density functional theory calculation and thermal annealing with ab initio molecular dynamics as motivated by experiment, we have studied the formation energy of Gr/Pb/SiC(0001) for different Pb coverages. Near the coverage of a Pb(111)-like ML mimicking the (10 x 10)-moiré, we find a slightly more stable stripe structure, where one half of the structure has compressive strain with Pb occupying the Si-top sites and the other half has tensile strain with Pb off the Si-top sites. This stripe structure along the Gr zigzag direction has a periodicity of 2.3 nm across the [1$\overline{2}$10] direction agreeing with the previous observations using scanning tunneling microscopy. Analysis with electron density difference and density of states show the tensile region has a more metallic character than the compressive region, while both are dominated by charge transfer from Pb ML to SiC(0001). As a result, the small energy difference between the stripe and Pb(111)-like structures means the two phases are almost degenerate and can coexist, which explains the experimental observations.

36 MATERIALS SCIENCE↗

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

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

36 - MATERIALS SCIENCE↗

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

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

Ncube, Musawenkosi K.↗

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

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

36 MATERIALS SCIENCE↗

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

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

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

Temperature as a control knob on spin-orbit coupling

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

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

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

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

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mechanisms of proton transport in aqueous acid solutions

The fundamental understanding of proton transport, specifically the Grotthuss-like mechanism, is critical for many technological applications. The present study presents a comprehensive analysis of proton transport in aqueous solutions of sulfuric and phosphoric acids using quasielastic neutron and light scattering, broadband dielectric spectroscopy, rheology, pulsed-field gradient-NMR measurements, and ab initio molecular dynamic simulations. Results show that in all systems, proton transport occurs through short “jumps” of ~0.5 Å that are faster than structural relaxation. However, proton hopping appears to be coupled to structural relaxation in aqueous solutions of sulfuric acid, while these processes are decoupled in phosphoric acid. Neutron scattering indicates that all protons have the same fast mobility in phosphoric acid systems, while at least one proton per sulfuric molecule remains slower than other protons. Analysis reveals correlated proton jumps, but these correlations suppress conductivity, suggesting that expected Grotthuss-like enhancement of conductivity is unlikely in bulk liquids.

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↗

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↗

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↗

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↗