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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

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

Hydrated Anions: From Clusters to Bulk Solution with Quasi-Chemical Theory

The interactions of hydrated ions with molecular and macromolecular solution and interface partners are strong on a chemical energy scale. Here we recount the foremost ab initio theory for the evaluation of the hydration free energies of ions, namely, quasi-chemical theory (QCT). We focus on anions, particularly halides but also the hydroxide anion, because they have been outstanding challenges for all theories. For example, this work supports understanding the high selectivity for F – over Cl – in fluoride-selective ion channels despite the identical charge and the size similarity of these ions. QCT is built by the identification of inner-shell clusters, separate treatment of those clusters, and then the integration of those results into the broader-scale solution environment. Recent work has focused on a close comparison with mass-spectrometric measurements of ion-hydration equilibria. We delineate how ab initio molecular dynamics (AIMD) calculations on ion-hydration clusters, elementary statistical thermodynamics, and electronic structure calculations on cluster structures sampled from the AIMD calculations obtain just the free energies extracted from the cluster experiments. That theory–experiment comparison has not been attempted before the work discussed here, but the agreement is excellent with moderate computational effort. This agreement reinforces both theory and experiment and provides a numerically accurate inner-shell contribution to QCT. The inner-shell complexes involving heavier halides display strikingly asymmetric hydration clusters. Asymmetric hydration structures can be problematic for the evaluation of the QCT outer-shell contribution with the polarizable continuum model (PCM). Nevertheless, QCT provides a favorable setting for the exploitation of PCM when the inner-shell material shields the ion from the outer solution environment. For the more asymmetrically hydrated, and thus less effectively shielded, heavier halide ions clustered with waters, the PCM is less satisfactory. We therefore investigate an inverse procedure in which the inner-shell structures are sampled from readily available AIMD calculations on the bulk solutions. This inverse procedure is a remarkable improvement; our final results are in close agreement with a standard tabulation of hydration free energies, and the final composite results are independent of the coordination number on the chemical energy scale of relevance, as they should be. Finally, a comparison of anion hydration structure in clusters and bulk solutions from AIMD simulations emphasize some differences: the asymmetries of bulk solution inner-shell structures are moderated compared with clusters but are still present, and inner hydration shells fill to slightly higher average coordination numbers in bulk solution than in clusters.

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

The Role of Cation Coordination in the Electrical and Optical Properties of Amorphous Transparent Conducting Oxides

Amorphous oxide semiconductor materials have demonstrated numerous advantages without compromise of electrical properties as compared to their crystalline counterparts, yet understanding of the fundamental principles allowing this has remained elusive. To study the origins of enhanced optoelectronic properties, we apply high-throughput, combinatorial sputtering, structural and spectral mapping, and computationally intensive ab initio molecular dynamics simulations with density functional theory to a ternary, post-transition metal oxide system, namely, zinc tin oxide. The deposited thin films exhibit a high figure of merit, achieving carrier densities in the range of 1019 to 1020 cm–3 and carrier mobilities up to 35 cm2/Vs. These results highlight the role of local distortions and cation coordination in determining the microscopic origins of carrier generation and transport. In particular, we identify the strong likelihood of Sn undercoordination in both Zn-poor and Zn-rich phases leading to the high carrier concentrations observed. This not only diverges from the still widespread historical indictment of oxygen vacancies controlling carrier population in crystalline oxides but also provides a comprehensive framework to describe the unique structure–property relationships using specific structural and electronic descriptors in disordered phase materials.

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS,M↗
DOE OSTI2020

Engineering Electrical Conductivity in Stable Zirconium-Based PCN-222 MOFs with Permanent Mesoporosity

Electrically conductive metal–organic frameworks (MOFs) featuring zirconium-based nodes are of great interest for electrochemical and optoelectronic applications due to their exceptional thermal and chemical stability, although the number of such MOFs remains limited. Here in this paper, we demonstrate that electron deficient molecules such as transition metal(IV) bis(dicarbollide) (TM IV (C 2 B 9 H 11 ) 2 ; TM IV CB for short, where TM = Ni, Pd, and Pt) or C 60 , spatially infiltrated in the microporous channels of the zirconium-based porphyrinic MOFs PCN-222 and PCN-222-Zn, generate highly stable and electrically conductive frameworks in which the mesopores remain accessible to other guests. Solid-state density functional theory (DFT) and ab initio molecular dynamics (AIMD) simulations indicate incorporation of Ni IV CB and C 60 in the MOF micropores is energetically and structurally feasible. Interestingly, in contrast to the free species, strong donor–acceptor interactions between the MOF and Ni IV CB restrain it to a single conformation. Calculated electronic structures and charge-hopping conduction probabilities illustrate that efficient charge-transfer (CT) from photoexcited linkers to the guest molecules facilitates charge hopping in the framework, making the MOFs electrically conductive. The donor–acceptor conjugates also enhance exciton dissociation at their heterojunctions, fostering the formation of long-lived electron-trapped states with potential utility for photo- and electrochemical devices.

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

Anionic Oxygen Redox in the High-Lithium Material Li 8 SnO 6

Lithiated transition metal oxide cathode materials that combine oxygen-anion redox [O 2 – /(O 2 ) n– ] with cation redox offer substantial capacities and higher voltage than cathodes without oxygen redox, but oxygen-anion redox is usually accompanied by the formation of peroxo or superoxo species that may cause oxygen release. To clarify the relationship between the oxygen release and peroxide and superoxide dimers without the complication of transition-metal redox, we performed density functional theory calculations to study the layered Li 8 SnO 6 cathode because it has only oxygen-anion redox during delithiation. Features with O–O distances of ~1.5 Å (corresponding to the formation of a peroxide) were observed at 75% lithium concentration, and further delithiation to 62.5% lithium concentration shortened the O–O bond distance to ~1.3 Å (corresponding to the formation of superoxide). The solely anionic oxygen redox stabilizes the delithiated structure of Li 8–x SnO 6 , and it does not cause structural disorder or release oxygen gas. Ab initio molecular dynamics calculations in the supercell at 300 K indicate that oxygen dimerization may occur at 75% lithium concentration. On the (0001) surface with a high lithium concentration (98%), oxygen dimerization is unfavorable. The calculations by HSE06 show that oxygen redox in Li 8 SnO 6 offers an initial voltage (vs Li/Li + ) greater than 4.2 V. The GGA+U calculations indicate that the Li + -polaron in Li 95 Sn 12 O 72 migrates with a migration barrier of only 0.43 eV. Here, this work also quantifies a solely anionic oxygen redox mechanism in Li 8 SnO 6 , and it provides a deeper understanding of anionic oxygen redox in Li-excess oxide cathode materials.

25 ENERGY STORAGE↗
DOE OSTI2021

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

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

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

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

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

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

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

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

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

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

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

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

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↗