Exploring the Mechanisms of LiNiO 2 Cathode Degradation by the Electrolyte Interfacial Deprotonation Reaction
Not Available
25 ENERGY STORAGE↗
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.
Characterization of Hole States at the Zn-Doped Hematite/Water Interface from Ab Initio Simulations
Not Available
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
A Perspective on Multiscale Modeling of Explicit Solvation-Enabled Simulations of Catalysis at Liquid–Solid Interfaces
Catalysis at liquid-solid interfaces is profoundly influenced by the interfacial solvent structure, which affects catalytic activity, selectivity, and reaction pathways. This perspective discusses state-of-the-art multiscale modeling methods that integrate quantum mechanics and molecular mechanics approaches to apply explicit solvent molecules to capture these interfacial phenomena. Specifically, the construction of multiscale models, the importance of capturing the interfacial solvent structure, and the computational strategies used to achieve this are explored, and the challenges in balancing chemical accuracy with computational expense are highlighted. Additionally, this perspective addresses the limitations of current methods. Opportunities for integrating machine learning are proposed. Here, by advancing the efficiency and user friendliness of multiscale modeling, it is argued that deeper insights into heterogeneous catalysis in liquid phases can be provided, which will ultimately contribute to the development of more efficient catalytic processes.
Ab initio molecular dynamics↗
Reduced Energy Cost of Ammonia Synthesis Via RF Plasma Pulsing
Not Available
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
X-ray and Optical Circular Dichroism as Local and Global Ultrafast Chiral Probes of [12]Helicene Racemization
Optical circular dichroism (OCD) is a well-established technique for probing molecular chirality typically in the UV-visible window. Although molecular chirality is usually related to a chiral center, it can alternatively describe a global structure. The electronic transitions are then delocalized and OCD is said to probe the global chirality. Recent advances controlling the polarization and spatial structure of X-ray beams have led to novel spectroscopic techniques such as X-ray circular dichroism (XCD), which can exploit the localized nature of the X-ray electronic transitions. XCD is thus sensitive to local structures and the relevant chirality probed with it can be referred as local. Here, we use the racemization mechanism of the [12]Helicene as a model to assess the capabilities of OCD and XCD as probes of global and local chiralities. In conclusion, our simulations suggest that XCD is a reliable spectroscopic technique for probing the local chirality of molecules, suggesting the further development of this experimental technique.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Transient Covalency in Molten Uranium(III) Chloride
Not Available
38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
Transferable Deep Learning Potential Reveals Intermediate-Range Ordering Effects in LiF–NaF–ZrF 4 Molten Salt
Not Available
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NWChem: Past, present, and future
Specialized computational chemistry packages have permanently reshaped the landscape of chemical sciences by providing tools to support and guide the experimental effort and for prediction of chemical and materials properties. In this regard, a special role has been played by electronic structure packages where complex chemical and materials processes can be modeled using first-principle-driven methodologies. Over the last few decades, the rapid development of computing technologies and a tremendous increase in computational power has offered a unique chance to study complex chemical transformations using sophisticated and predictive many-body techniques to describe correlated behavior of electrons in molecular and condensed phase systems at different levels of theory. In enabling these simulations, a critical role has been played by novel parallel algorithms capable of taking advantage of computational resources to address polynomial scaling of electronic structure methods. NWChem was among the first electronic structure codes that focused on delivering scalable performance for electronic structure simulations. Herein, we briefly review the NWChem suite of computational codes including its history, design principles, parallel tools, current capabilities, outreach and outlook.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Localized high concentration electrolytes decomposition under electron-rich environments
Localized high concentration electrolytes have been proposed as an effective route to construct stable solid–electrolyte interphase (SEI) layers near Li-metal anodes. However, there is still a limited understanding of decomposition mechanisms of electrolyte components during SEI formation. In this work, we investigate reactivities of lithium bis(fluorosulfonyl)imide (LiFSI, salt), 1,2-dimethoxyethane (DME, solvent) and tris(2,2,2-trifluoroethyl)orthoformate (TFEO, diluent) species in DME+TFEO mixed solvents and 1M LiFSI/DME/TFEO solutions. By supplying excess of electrons into the simulation cell, LiFSI is initially reduced via a 4-electron charge transferring reaction yielding F - and N(SO 2 ) 2 3- . The local solvation environment has little effect on the subsequent TFEO reaction, which typically requires 6 |e| to decompose into F - , HCOO - , CH2CF and - OCH 2 CF 3 . Besides, TFEO dehydrogenation reaction mechanism under the attack of anions is also identified. Unlike salt and diluent, DME shows good stability with any excess of electrons. Here, the energetics of most relevant reactions are characterized. Most reactions are thermodynamically favorable with low activation barriers.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
The behavior of methane–water mixtures under elevated pressures from simulations using many-body potentials
Non-polarizable empirical potentials have been proven to be incapable of capturing the mixing of methane–water mixtures at elevated pressures. Although density functional theory-based ab initio simulations may circumvent this discrepancy, they are limited in terms of the relevant time and length scales associated with mixing phenomena. Here, we show that the many-body MB-nrg potential, designed to reproduce methane–water interactions with coupled cluster accuracy, successfully captures this phenomenon up to 3 GPa and 500 K with varying methane concentrations. Two-phase simulations and long time scales that are required to fully capture the mixing, affordable due to the speed and accuracy of the MBX software, are assessed. Constructing the methane–water equation of state across the phase diagram shows that the stable mixtures are denser than the sum of their parts at a given pressure and temperature. As a result, we find that many-body polarization plays a central role, enhancing the induced dipole moments of methane by 0.20 D during mixing under pressure. Overall, the mixed system adopts a denser state, which involves a significant enthalpic driving force as elucidated by a systematic many-body energy decomposition analysis.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
GPU acceleration of local and semilocal density functional calculations in the SPARC electronic structure code
We present a Graphics Processing Unit (GPU)-accelerated version of the real-space SPARC electronic structure code for performing Kohn–Sham density functional theory calculations within the local density and generalized gradient approximations. In particular, we develop a modular math-kernel based implementation for NVIDIA architectures wherein the computationally expensive operations are carried out on the GPUs, with the remainder of the workload retained on the central processing units (CPUs). Here, using representative bulk and slab examples, we show that relative to CPU-only execution, GPUs enable speedups of up to 6× and 60× in node and core hours, respectively, bringing time to solution down to less than 30 s for a metallic system with over 14 000 electrons and enabling significant reductions in computational resources required for a given wall time.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Extreme sensitivity of higher-order interatomic force constants and thermal conductivity to the energy surface roughness of exchange-correlation functionals
In this Letter, we report that the fourth-order interatomic force constants (4th-IFCs) are significantly sensitive to the energy surface roughness of exchange-correlation (XC) functionals in density functional theory calculations. This sensitivity, which is insignificant for the second- (2nd-) and third-order (3rd-) IFCs, varies for different functionals in different materials and can cause misprediction of thermal conductivity by several times of magnitude. As a result, when calculating the 4th-IFCs using the finite difference method, the atomic displacement needs to be taken large enough to overcome the energy surface roughness, in order to accurately predict phonon lifetime and thermal conductivity. We demonstrate this phenomenon on a benchmark material (Si), a high-thermal conductivity material (BAs), and a low thermal conductivity material (NaCl). For Si, we find that the LDA, PBE, and PBEsol XC functionals are all smooth to the 2nd- and 3rd-IFCs but all rough to the 4th-IFCs. This roughness can lead to a prediction of nearly one order of magnitude lower thermal conductivity. For BAs, all three functionals are smooth to the 2nd- and 3rd-IFCs, and only the PBEsol XC functional is rough for the 4th-IFCs, which leads to a 40% underestimation of thermal conductivity. For NaCl, all functionals are smooth to the 2nd- and 3rd-IFCs but rough to the 4th-IFCs, leading to a 70% underprediction of thermal conductivity at room temperature. In conclusion, with these observations, we provide general guidance on the calculation of 4th-IFCs for an accurate thermal conductivity prediction.
42 ENGINEERING↗
Palladium at high pressure and high temperature: A combined experimental and theoretical study
Palladium is one of the most important technological materials, yet its phase diagram remains poorly understood. At ambient conditions, its solid phase is face-centered cubic (fcc). However, another solid phase of Pd, body-centered cubic (bcc), was very recently predicted in two independent theoretical studies to occur at high pressures and temperatures. In this work, we report an experimental study on the room-temperature equation of state (EOS) of Pd to a pressure of 80 GPa, as well as a theoretical study on the phase diagram of Pd including both fcc-Pd and bcc-Pd. Our theoretical approach consists in ab initio quantum molecular dynamics (QMD) simulations based on the Z methodology which combines both direct Z method for the simulation of melting curves and inverse Z method for the calculation of solid–solid phase transition boundaries. We obtain the melting curves of both fcc-Pd and bcc-Pd and an equation for the fcc–bcc solid–solid phase transition boundary as well as the thermal EOS of Pd which is in agreement with experimental data and QMD simulations. We uncover the presence of another solid phase of Pd on its phase diagram, namely, random hexagonal close-packed (rhcp), and estimate the location of the rhcp-bcc solid–solid phase transition boundary and the rhcp–fcc–bcc triple point. We also discuss the topological similarity of the phase diagrams of palladium and silver, the neighbor of Pd in the periodic table. We argue that Pd is a reliable standard for shock-compression studies and present the analytic model of its principal Hugoniot in a wide pressure range.
36 MATERIALS SCIENCE↗
Tunable noninteracting free-energy density functionals for high-energy-density physics applications
In this work, we introduce the concept of a tunable noninteracting free-energy density functional and present two examples realized: (i) via a simple one-parameter convex combination of two existing functionals and (ii) via the construction of a generalized gradient approximation (GGA) enhancement factor that contains one free parameter and is designed to satisfy a set of incorporated constraints. Functional (i), constructed as a combination of the local Thomas–Fermi and a pseudopotential-adapted GGA for the noninteracting free-energy, has already demonstrated its practical usability for establishing the high temperature end of the equation of state of deuterium [Phys. Rev. B 104, 144104 (2021)] and CHON resin [Phys. Rev. E 106, 045207 (2022)] for inertial confinement fusion applications. Hugoniot calculations for liquid deuterium are given as another example of how the application of computationally efficient orbital-free density functional theory (OF-DFT) can be utilized with the employment of the developed functionals. Once the functionals have been tuned such that the OF-DFT Hugoniot calculation matches the Kohn–Sham solution at some low-temperature point, agreement with the reference Kohn–Sham results for the rest of the high temperature Hugoniot path is very good with relative errors for compression and pressure on the order of 2% or less.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
How the choice of exchange–correlation functional affects DFT-based simulations of the hydrated electron
Hydrated electrons are anionic species that are formed when an excess electron is introduced into liquid water. Building an understanding of how hydrated electrons behave in solution has been a long-standing effort of simulation methods, of which density functional theory (DFT) has come to the fore in recent years. The ability of DFT to model the reactive chemistry of hydrated electrons is an attractive advantage over semi-classical methodologies; however, relatively few density functional approximations (DFAs) have been used for the hydrated electron simulations presented in the literature. Here, we simulate hydrated electron systems using a series of exchange–correlation (XC) functionals spanning Jacob’s ladder. We calculate a variety of experimental and other observables of the hydrated electron and compare the XC functional dependence for each quantity. We find that the formation of a stable localized hydrated electron is not necessarily limited to hybrid XC functionals and that some hybrid functionals produce delocalized hydrated electrons or electrons that react with the surrounding water at an unphysically fast rate. Here, we further characterize how different DFAs impact the solvent structure and predicted spectroscopy of the hydrated electron, considering several methods for calculating the hydrated electron’s absorption spectrum for the best comparison between structures generated using different density functionals. None of the dozen or so DFAs that we investigated are able to correctly predict the hydrated electron’s spectroscopy, vertical detachment energy, or molar solvation volume.
Ab-initio molecular dynamics↗
GPU acceleration of hybrid functional calculations in the SPARC electronic structure code
We present a Graphics Processing Unit (GPU)-accelerated version of the real-space SPARC electronic structure code for performing hybrid functional calculations in generalized Kohn–Sham density functional theory. In particular, we develop a batch variant of the recently formulated Kronecker product-based linear solver for the simultaneous solution of multiple linear systems. We then develop a modular, math kernel based implementation for hybrid functionals on NVIDIA architectures, where computationally intensive operations are offloaded to the GPUs, while the remaining workload is handled by the central processing units (CPUs). Considering bulk and slab examples, we demonstrate that GPUs enable up to 8× speedup in node-hours and 80× in core-hours compared to CPU-only execution, reducing the time to solution on V100 GPUs to around 300 s for a metallic system with over 6000 electrons, and significantly reducing the computational resources required for a given wall time.
Kohn-Sham density functional theory↗
Temperature-induced structural evolution in liquid Ag-Ga alloys
Temperature-dependent atomic structural evolutions of liquid Ag 60 Ga 40 and Ag 70 Ga 30 alloys have been studied by in situ high-energy x-ray-diffraction (HEXRD) experiments combined with ab initio molecular-dynamics simulations. The experimental data show a reversible structural crossover at about 1050 ~ 1100 (± 50) K in both liquid Ag 60 Ga 40 and Ag 70 Ga 30 alloys. Obvious changes of the electrical resistivity, absolute thermoelectric power, and atomic diffusivity around the similar temperature range for both Ag-Ga liquids strongly support the HEXRD results. The origin of the liquid-to-liquid crossover in both Ag-Ga liquids was suggested to link with the rearrangements of Ag and Ga atoms, i.e., Ag and Ga atoms prefer to associate with themselves in the higher temperature range above 1100 K, consistent with the accelerated increase of the strong covalently bonded Ga–Ga dimers in both Ag-Ga liquids. In addition, more studies from the energy aspect are still desirable to understand the rearrangements of Ag and Ga atoms in the higher temperature range in both Ag-Ga liquids.
36 MATERIALS SCIENCE↗
Topological electride phase of sodium at high pressures and temperatures
Sodium, a textbook example of a nearly-free-electron metal, exhibits unforeseen pressure-induced behavior including the stabilization of numerous polymorphs—some possessing extremely complex unit cells—as well as a metal-to-insulator transition to the iconic hP4 phase. However, until recently, most of the experimental and theoretical studies on solid sodium have been restricted to the low-temperature regime. In this work, using ab initio evolutionary structure searches coupled with quasiharmonic calculations, we discover seven new phases of sodium that are more stable than the known hP4 phase at pressure-temperature conditions that were recently attained in ramp-compression experiments. From these, our calculations suggest that P6 3 /m Na is the ground state between ≈250 GPa at 710 K and ≈350 GPa at 1900 K. Electronic structure calculations show that this phase is a topological semimetal with a Dirac nodal surface that is protected by a nonsymmorphic symmetry S 2z and an electride owing to its non-nuclear charge localized within one-dimensional honeycomb channels and zero-dimensional cages. Our results highlight the complexity of dense sodium's free-energy landscape and intricate electronic structure that emerges at finite temperatures and conditions where ionic cores overlap.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗