First-Principles and Thermodynamics Study of Compositionally Tuned Complex Metal Oxides: Cation Release from the (001) Surface of Mn-Rich Lithium Nickel Manganese Cobalt Oxide

Synthesis of Single Crystal Li2NpO4 and Li4NpO5 from Aqueous Lithium Hydroxide Solutions under Mild Hydrothermal Conditions

The ternary oxides, Li2NpO4 and Li4NpO5, were synthesized under mild hydrothermal conditions using concentrated LiOH solutions containing NpO2(NO3)2. The reactions resulted in the formation of single crystals of both compounds, enabling the determination of their single crystal structures for the first time. Exploration of the synthetic phase space demonstrates that the resulting neptunate phases are dependent on the concentration of LiOH, transitioning from Li2NpO4, containing a typical octahedral neptunyl geometry with two shorter Np≡O bonds, at lower LiOH concentrations to Li4NpO5 with two long and four short Np-O bonds under saturated solution conditions. Reactions exploring the same synthetic conditions are also reported for uranyl(VI) for comparison. Raman spectra of the compounds were collected and analyzed to evaluate the Np-O bonding in these compounds.

Guiding Principles for the Rational Design of Hybrid Materials: Use of DFT Methodology for Evaluating Non-Covalent Interactions in a Uranyl Tetrahalide Model System

Together with the synthesis and experimental characterization of 14 hybrid materials containing [UO2 X4 ]2- (X=Cl- and Br- ) and organic cations, we report on novel methods for determining correlation trends in their formation enthalpy (ΔHf ) and observed vibrational signatures. ΔHf values were analyzed through isothermal acid calorimetry and a Density Functional Theory+Thermodynamics (DFT+T) approach with results showing good agreement between theory and experiment. Three factors (packing efficiency, cation protonation enthalpy, and hydrogen bonding energy [E H , norm total ${{E}_{H,{\rm { norm}}}^{{\rm { total}}}}$ ]) were assessed as descriptors for trends in ΔHf . Results demonstrated a strong correlation betweenE H , norm total ${E_{{\rm{H}},{\rm{norm}}}^{{\rm{total}}} }$ and ΔHf , highlighting the importance of hydrogen bonding networks in determining the relative stability of solid-state hybrid materials. Lastly, we investigate how hydrogen bonding networks affect the vibrational characteristics of uranyl solid-state materials using experimental Raman and IR spectroscopy and theoretical bond orders and find that hydrogen bonding can red-shift U≡O stretching modes. Overall, the tightly integrated experimental and theoretical studies presented here bridge the trends in macroscopic thermodynamic energies and spectroscopic features with molecular-level details of the geometry and electronic structure. This modeling framework forms a basis for exploring 3D hydrogen bonding as a tunable design feature in the pursuit of supramolecular materials by rational design.

Chemical Transformations of 2D Kaolinic Clay Mineral Surfaces from Sulfuric Acid Exposure

A combined experimental and computational approach is used to investigate the chemical transformations of kaolinite and metakaolin surfaces when exposed to sulfuric acid. These clay minerals are hydrated ternary metal oxides and are shown to be susceptible to degradation by loss of Al as the water-soluble salt Al₂(SO₄)₃, due to interactions between H₂SO₄ and aluminum cations. This degradation process results in a silica-rich interfacial layer on the surfaces of the aluminosilicates, most prominently observed in metakaolin exposed to pH environments of less than 4. Our observations are supported by XPS, ATR-FTIR, and XRD experiments. Concurrently, DFT methodologies are used to probe the interactions between the clay mineral surfaces and H₂SO₄ as well as other sulfur-containing adsorbates. An analysis performed using a DFT + thermodynamics model shows that the surface transformation processes that lead to the loss of Al and SO₄ from metakaolin are favorable at pH below 4; however, such transformations are not favorable for kaolinite, a result that agrees with our experimental efforts. The data obtained from both experimental techniques and computational studies support that the dehydrated surface of metakaolin interacts more strongly with sulfuric acid and provide atomistic insight into the acid-induced transformations of these mineral surfaces.

Periodic Density Functional Theory Calculations of Uranyl Tetrachloride Compounds Engaged in Uranyl-Cation and Uranyl-Hydrogen Interactions: Electronic Structure, Vibrational, and Thermodynamic Analyses

Solid-state uranyl hybrid structures are often formed through unique intermolecular interactions occurring between a molecular uranyl anion and a charge-balancing cation. In this work, solid-state structures of the uranyl tetrachloride anion engaged in uranyl-cation and uranyl-hydrogen interactions were studied using density functional theory (DFT). As most first-principles methods used for systems of this type focus primarily on the molecular structure, we present an extensive benchmarking study to understand the methods needed to accurately model the geometric properties of these systems. From there, the electronic and vibrational structures of the compounds were investigated through projected density of states and phonon analysis and compared to the experiment. Lastly, we present a DFT + thermodynamics approach to calculate the formation enthalpies (ΔHf) of these systems to directly relate to experimental values. Through this methodology, we were able to accurately capture trends observed in experimental results and saw good quantitative agreement in predicted ΔHf compared to the value calculated through referencing each structure to its standard state. Overall, results from this work will be used for future combined experimental and computational studies on both uranyl and neptunyl hybrid structures to delineate how varying intermolecular interaction strengths relates to the overall values of ΔHf.

Metal Release Mechanism and Electrochemical Properties of Lix(Ni1/3Mn1/3Co1/3)O2

Complex metal oxides (CMOs) are used broadly in applications including electroreactive forms found in lithium-ion battery technology. Computational chemistry can provide unique information about how the properties of CMO cathode materials change in response to changes in stoichiometry, for example, changes of the lithium (Li) content during the charge–discharge cycle of the battery. However, this is difficult to measure experimentally due to the small cross-sectional area of the cations. Outside of operational conditions, the Li content can influence the transformations of the CMO when exposed to the environment. For example, metal release from CMOs in aqueous settings has been identified as a cross-cutting mechanism important to CMO degradation. Computational studies investigating metal release from CMOs show that the thermodynamics depend on the oxidation states of lattice cations, which is expected to vary with the lithium content. In this work, computational studies track changes in metal release trends as a function of Li content in Lix(Ni1/3Mn1/3Co1/3)O2 (NMC). The resulting dataset is used to construct a random forest tree (RFT) machine learning (ML) model. A modeling challenge in delithiation studies is the large configurational space to sample. Through investigating multiple configurations at each lithium fraction, we find structural features associated with favorable energies to chemically guide the identification of relevant structures and adequately predict voltage values.

Density Functional Theory (DFT) as a Nondestructive Probe in the Field of Art Conservation: Small-Molecule Adsorption on Aragonite Surfaces

Humans have incorporated minerals in objects of cultural heritage importance for millennia. The surfaces of these objects, which often long outlast the humans that create them, are undeniably exposed to a diverse mixture of chemicals throughout their lifetimes. As of yet, the art conservation community lacks a nondestructive, accurate, and inexpensive flexible computational screening method to evaluate the potential impact of chemicals with art, as a complement to experimental studies. In this work, we propose periodic density functional theory (DFT) studies as a way to address this challenge, specifically for the aragonite phase of calcium carbonate, a mineral that has been used in pigments, marble statues, and limestone architecture since ancient times. Computational models allow art conservation scientists to better understand the atomistic impact of small-molecule adsorbates on common mineral surfaces across a wide variety of environmental conditions. To gain insight into the surface adsorption reactivity of aragonite, we use DFT to investigate the atomistic interactions present in small-molecule-surface interfaces. Our adsorbate set includes common solvents, atmospheric pollutants, and emerging contaminants. Chemicals that significantly disrupt the surface structure such as carboxylic acids and sulfur-containing molecules are highlighted. We also focus on comparing adsorption energies and changes in surface bonds, which allows for the identification of key features in the electronic structure presented in a projected-density-of-state analysis. The trends outlined here will guide future experiments and allow art conservators to gain a better understanding of how a wide range of molecules interact with an aragonite surface under variable conditions and in different environments.

Eliminating Delocalization Error to Improve Heterogeneous Catalysis Predictions with Molecular DFT + U

Approximate semilocal density functional theory (DFT) is known to underestimate surface formation energies yet paradoxically overbind adsorbates on catalytic transition-metal oxide surfaces due to delocalization error. The low-cost DFT + U approach only improves surface formation energies for early transition-metal oxides or adsorption energies for late transition-metal oxides. In this work, we demonstrate that this inefficacy arises due to the conventional usage of metal-centered atomic orbitals as projectors within DFT + U. We analyze electron density rearrangement during surface formation and O atom adsorption on rutile transition-metal oxides to highlight that a standard DFT + U correction fails to tune properties when the corresponding density rearrangement is highly delocalized across both metal and oxygen sites. To improve both surface properties simultaneously while retaining the simplicity of a single-site DFT + U correction, we systematically construct multi-atom-centered molecular-orbital-like projectors for DFT + U. We demonstrate this molecular DFT + U approach for tuning adsorption energies and surface formation energies of minimal two-dimensional models of representative early (i.e., TiO₂) and late (i.e., PtO₂) transition-metal oxides. Molecular DFT + U simultaneously corrects adsorption energies and surface formation energies of multilayer models of rutile TiO₂(110) and PtO₂(110) to resolve the paradoxical description of surface stability and surface reactivity of semilocal DFT.

Adsorption of small organic acids and polyphenols on hematite surfaces: Density Functional Theory + thermodynamics analysis

HYPOTHESIS

The interactions of organic molecules with mineral surfaces are influenced by several factors such as adsorbate speciation, surface atomic and electronic structure, and environmental conditions. When coupled with thermodynamic techniques, energetics from atomistic modeling can provide a molecular-level picture of which factors determine reactivity. This is paramount for evaluating the chemical processes which control the fate of these species in the environment.

EXPERIMENTS

Inner-sphere adsorption of oxalate and pyrocatechol on (001), (110), and (012) α-Fe₂O₃ surfaces was modeled using Density Functional Theory (DFT). Unique bidentate binding modes were sampled along each facet to study how different adsorbate and surface factors govern site preference. Adsorption energetics were then calculated using a DFT + thermodynamics approach which combines DFT energies with tabulated data and Nernst-based corrective terms to incorporate different experimental parameters.

FINDINGS

Instead of a universal trend, each facet displays a unique factor that dominates site preference based on either strain (001), functional groups (110), or topography (012). Adsorption energies predict favorable inner-sphere adsorption for both molecules but opposite energetic trends with varying pH. Additionally, vibrational analysis was conducted for each system and compared to experimental IR data. The work presented here provides an effective, computational methodology to study numerous adsorption processes occurring at the surface-aqueous interface.

Density functional theory and thermodynamics analysis of MAl12 Keggin substitution reactions: Insights into ion incorporation and experimental confirmation

Polyaluminum cations, such as the MAl₁₂ Keggin, undergo atomic substitutions at the heteroatom site (M), where nanoclusters with M = Al³⁺, Ga³⁺, and Ge⁴⁺ have been experimentally studied. The identity of the heteroatom M has been shown to influence the structural and electronic properties of the nanocluster and the kinetics of ligand exchange reactions. To date, only three ε-analogs have been identified, and there is a need for a predictive model to guide experiment to the discovery of new MAl₁₂ species. Here, we present a density functional theory (DFT) and thermodynamics approach to predicting favorable heteroatom substitution reactions, alongside structural analyses on hypothetical ε-MAl₁₂ nanocluster models. We delineate trends in energetics and geometry based on heteroatom cation properties, finding that Al³⁺-O bond lengths are related to heteroatom cation size, charge, and speciation. Our analyses also enable us to identify potentially isolable new ε-MAl₁₂ species, such as FeAl₁₂ ⁷⁺. Based upon these results, we evaluated the Al³⁺/Zn²⁺/Cr³⁺ system and determined that substitution of Cr³⁺ is unfavorable in the heteroatom site but is preferred for Zn²⁺, in agreement with the experimental structures. Complimentary experimental studies resulted in the isolation of Cr³⁺-substituted δ-Keggin species where Cr³⁺ substitution occurs only in the octahedral positions. The isolated structures Na[AlO₄Al_9.6Cr_2.4(OH)₂₄(H₂O)₁₂](2,6-NDS)₄(H₂O)₂₂ (δ-Cr_nAl_13-n-1) and Na[AlO₄Al_9.5Cr_2.5(OH)₂₄(H₂O)₁₂](2,7-NDS)₄(H₂O)_18.5 (δ-Cr_nAl_13-n-2) are the first pieces of evidence of mixed Al³⁺/Cr³⁺ Keggin-type nanoclusters that prefer substitution at the octahedral sites. The δ-Cr_nAl_13-n-2 structure also exhibits a unique placement of the bound Na⁺ cation, which may indicate that Cr³⁺ substitution can alter the surface reactivity of Keggin-type species.

Surface Transformations of Lead Oxides and Carbonates Using First-Principles and Thermodynamics Calculations

Lead (Pb)-containing solids find widespread commercial use in batteries, piezoelectrics, and as starting materials for synthesis. Here, we combine density functional theory (DFT) and thermodynamics in a DFT + solvent ion model to compare the surface reactivity of Pb oxides and carbonates, specifically litharge, massicot, and cerussite, in contact with water. The information provided by this model is used to delineate structure-property relationships for surfaces that are able to release Pb as Pb²⁺. We find that Pb²⁺ release is dependent on pH and chemical bonding environment and go on to correlate changes in the surface bonding to key features of the electronic structure through a projected density of states analysis. Collectively, our analyses link the atomistic structure to i) specific electronic states and ii) the thermodynamics of surface transformations, and the results presented here can be used to guide synthetic efforts of Pb²⁺-containing materials in aqueous media or be used to better understand the initial steps in solid decomposition.

Advanced material modulation of nutritional and phytohormone status alleviates damage from soybean sudden death syndrome

Customized Cu3(PO4)2 and CuO nanosheets and commercial CuO nanoparticles were investigated for micronutrient delivery and suppression of soybean sudden death syndrome. An ab initio thermodynamics approach modelled how material morphology and matrix effects control the nutrient release. Infection reduced the biomass and photosynthesis by 70.3 and 60%, respectively; the foliar application of nanoscale Cu reversed this damage. Disease-induced changes in the antioxidant enzyme activity and fatty acid profile were also alleviated by Cu amendment. The transcription of two dozen defence- and health-related genes correlates a nanoscale Cu-enhanced innate disease response to reduced pathogenicity and increased growth. Cu-based nanosheets exhibited a greater disease suppression than that of CuO nanoparticles due to a greater leaf surface affinity and Cu dissolution, as determined computationally and experimentally. The findings highlight the importance and tunability of nanomaterial properties, such as morphology, composition and dissolution. The early seedling foliar application of nanoscale Cu to modulate nutrition and enhance immunity offers a great potential for sustainable agriculture.

Density Functional Theory and Thermodynamics Modeling of Inner-Sphere Oxyanion Adsorption on the Hydroxylated α-Al2O3(001) Surface

The inner-sphere adsorption of AsO₄^3-, PO₄^3-, and SO₄^2- on the hydroxylated α-Al₂O₃(001) surface was modeled with the goal of adapting a density functional theory (DFT) and thermodynamics framework for calculating the adsorption energetics. While DFT is a reliable method for predicting various properties of solids, including crystalline materials comprised of hundreds (or even thousands) of atoms, adding aqueous energetics in heterogeneous systems poses steep challenges for modeling. This is in part due to the fact that environmentally relevant variations in the chemical surroundings cannot be captured atomistically without increasing the system size beyond tractable limits. The DFT + thermodynamics approach to this conundrum is to combine the DFT total energies with tabulated solution-phase data and Nernst-based corrective terms to incorporate experimentally tunable parameters such as concentration. Central to this approach is the design of thermodynamic cycles that partition the overall reaction (here, inner-sphere adsorption proceeding via ligand exchange) into elementary steps that can either be fully calculated or for which tabulated data are available. The ultimate goal is to develop a modeling framework that takes into account subtleties of the substrate (such as adsorption-induced surface relaxation) and energies associated with the aqueous environment such that adsorption at mineral-water interfaces can be reliably predicted, allowing for comparisons in the denticity and protonation state of the adsorbing species. Based on the relative amount of experimental information available for AsO₄^3-, PO₄^3-, and SO₄^2- adsorbates and the well-characterized hydroxylated α-Al₂O₃(001) surface, these systems are chosen to form a basis for assessing the model predictions. We discuss how the DFT + thermodynamics results are in line with the experimental information about the oxyanion sorption behavior. Additionally, a vibrational analysis was conducted for the charge-neutral oxyanion complexes and is compared to the available experimental findings to discern the inner-sphere adsorption phonon modes. The DFT + thermodynamics framework used here is readily extendable to other chemical processes at solid-liquid interfaces, and we discuss future directions for modeling surface processes at mineral-water and environmental interfaces.

Ga3+ Incorporation into Al13 Keggin Polyoxometalates and the Formation of δ-(GaAl12)7+ and (Ga2.5Al28.5)19+ Polycations

Keggin-type polyaluminum species (ε-Al₁₃, δ-Al₁₃, Al₂₆, Al₃₀, Al₃₂) can form upon partial hydrolysis of Al³⁺-bearing solutions and are important species for water purification and contaminant transport. While the structural features for the major Al³⁺ polyaluminum species have been delineated, much less is known regarding heteroatom substitution and resultant structures other than the previously identified ε-GaAl₁₂⁷⁺ and ε-GeAl₁₂⁸⁺ cations. Single-atom substitution within polyaluminum species can change the surface reactivity within water treatment scenarios; thus, it is important to understand heteroatom incorporation within this system. The present work describes the synthesis and characterization of two novel Ga³⁺-substituted Keggin-type polyaluminum species. Na[GaO₄Al₁₂(OH)₂₄(H₂O)₁₂](2,6-NDS)₄(H₂O)_20.5 (δ-GaAl₁₂) and [Ga₂O₈Al_28.5Ga_0.5(OH)₅₈(H₂O)₂₇(SO₄)₂](SO₄)₄Cl₇(H₂O)_8.5 (Ga_2.5Al_28.5) were crystallized from a thermally aged, partially hydrolyzed Ga³⁺/Al³⁺ solution. Structural refinement from single-crystal X-ray diffraction indicated fully occupied Ga³⁺ within tetrahedral site(s) of both isolated species. Partial substitution was observed for octahedral sites for the larger Ga_2.5Al_28.5 cluster. The chemical compositions of both clusters were confirmed by inductively coupled plasma mass spectrometry (ICP-MS). Density functional theory (DFT) calculations corroborated the structural refinement, with the energetics of Ga³⁺ substitution suggesting preferential substitution within tetrahedral sites for both species. Additional theoretical work suggests that the rotated trimer in δ-GaAl₁₂ is highly reactive, which can serve as the driving force in the formation of the Ga_2.5Al_28.5 cluster.

Energy Storage Materials as Emerging Nano-contaminants

New and emerging nanotechnologies are increasingly using nanomaterials that undergo significant chemical reactions upon exposure to environmental conditions. The rapid advent of lithium ion batteries for energy storage in mobile electronics and electric vehicles is leading to rapid increases in the manufacture of complex transition metal oxides that incorporate elements such as Co and Ni that have the potential for significant adverse biological impact. This Perspective summarizes some of the important technological drivers behind complex oxide materials and highlights some of the chemical transformations that need to be understood in order to assess the overall environmental impact associated with energy storage technologies.

Insights into the Cathode-Electrolyte Interphases of High-Energy-Density Cathodes in Lithium-Ion Batteries

We present a comprehensive study of cycled high-Ni (LiNi_1-xM_xO₂, M = metals), Li-rich (Li_1+xMn_yM_1-x-yO₂), and high-voltage spinel (LiMn_1.5Ni_0.5O₄) electrodes with time-of-flight secondary ion mass spectrometry (TOF-SIMS) and X-ray photoelectron spectroscopy in conjunction with electrochemical techniques to better understand their evolving cathode-electrolyte interphase structure during cycling. TOF-SIMS provides fragment-specific information regarding the surface film content for each of the electrodes. High-Ni cathodes show thick surface films initially containing Li₂CO₃, later developing oxidized organic carbonates throughout cycling. Li-rich electrode surface films develop strong characteristics during their first activation cycles, where released O₂ oxidizes organic carbonates to form polymeric carbons and decomposes LiPF₆. High-voltage spinel electrodes operate outside the standard electrolyte stability window, generating reactive oxidized electrolyte species that further decompose LiPF₆. The distribution and concentration of these different chemical fragments measured by TOF-SIMS are finally summarized by color-coded high-resolution images of cycled high-Ni, Li-rich, and high-voltage spinel electrodes.