Occupation matrix control of d- and f-electron localisations using DFT + U

The energy landscape of magnetic materials

Magnetic materials can display many solutions to the electronic-structure problem, corresponding to different local or global minima of the energy functional. In Hartree-Fock or density-functional theory different single-determinant solutions lead to different magnetizations, ionic oxidation states, hybridizations, and inter-site magnetic couplings. The vast majority of these states can be fingerprinted through their projection on the atomic orbitals of the magnetic ions. We have devised an approach that provides an effective control over these occupation matrices, allowing us to systematically explore the landscape of the potential energy surface. We showcase the emergence of a complex zoology of self-consistent states; even more so when semi-local density-functional theory is augmented - and typically made more accurate - by Hubbard corrections. Such extensive explorations allow to robustly identify the ground state of magnetic systems, and to assess the accuracy (or not) of current functionals and approximations.

Hybridization State Transition under Working Conditions: Activity Origin of Single-Atom Catalysts

Single-atom catalysts (SACs) have been widely investigated and have emerged as a transformative approach in electrocatalysis. Despite their clear structure, the origin of their exceptional activity remains elusive. Herein, we elucidate a common phenomenon of the hybridization state transition of metal centers, which is responsible for the activity origin across various SACs for different reactions. Focusing on N-doped carbon-supported Ni SAC (NiN4 SAC) for CO2 reduction reaction (CO2RR), our comprehensive computations successfully clarify the hybridization state transition under working conditions and its relation with the activity. This transition, triggered by the reaction intermediates and applied potential, converts the Ni center from the inert dsp2 hybridization state to the active d2sp3 hybridization state. Importantly, the calculated activity and selectivity of the CO2RR over the d2sp3-hybridized Ni center are consistent with the experimental results, offering strong support for the proposed hypothesis. This work suggests a universal principle of electronic structure evolution in SACs that could revolutionize catalyst design, which also introduces a new paradigm for manipulating electronic states to enhance catalytic performance, with implications for various reactions and catalyst platforms.

Role of Dy 4felectrons on magnetic coupling and reorientation in DyFeO3

Spin reorientation transition is an ubiquitous phenomenon observed in magnetic rare earth orthferrites RFeO3, which has garnered significant attention in recent years due to its potential applications in spintronics or magnetoelectric devices. Although a plenty of experimental works suggest that the magnetic interaction between R3+and Fe3+spins is at the heart of the spin reorientation, but a direct and conclusive theoretical support has been lacking thus far, primarily due to the challenging nature of handling R 4felectrons. In this paper, we explored DyFeO3as an example by means of comprehensive first principles calculations, and compared two different approaches, where the Dy 4felectrons were treated separately as core or valence states, aiming to elucidate the role of Dy 4felectrons, particularly in the context of the spin reorientation transition. The comparison provides a solid piece of evidence for the experimental argument that the Dy3+-Fe3+magnetic interactions play a vital role in triggering spin reorientation of Fe3+moments at low temperatures. The findings revealed here not only extend our understanding on the underlying mechanism for spin reorientation transition in RFeO3, but also highlight the importance of explicit description of R 4felectrons in rationally reproducing their structural, electronic and magnetic properties.

Unexpected Two-Dimensional Polarons Induced by Oxygen Vacancies in Layered Structure MoO3-x

Unveiling the effects of oxygen vacancies on the structural stability of layered α-MoO3 is critical for optimizing its physical and chemical properties. Herein, we present experimental evidence regarding the phase stability of α-MoO3 with ∼2% oxygen vacancy concentrations. Interestingly, we report a previously ignored oxygen-deficient orthorhombic MoO3-x phase in space group Cmcm. Further density functional theory calculations reveal a detailed phase transition mechanism from α-MoO3 to MoO3-x. More importantly, we demonstrate that two-dimensional (2D) large polarons must exist to stabilize the MoO3-x crystal structure. 2D large polarons are suspected to exist in numerous quasi-2D systems but have never been found in layered α-MoO3 or other molybdenum oxides. Our work contributes to a basic understanding of the polaronic behavior in MoO3-x and may broaden the application realm of molybdenum oxides.

Time-Domain View of Polaron Dynamics in Metal Oxide Photocatalysts

The polaron is a fundamental physical phenomenon in transition metal oxides (TMOs), and it has been studied extensively for decades. However, the implication of a polaron on photochemistry is still ambiguous. As such, understanding the fundamental properties and controlling the dynamics of polarons at the atomistic level is desired. In this Perspective, we seek to highlight the recent advances in studying small polarons in TMOs, with a particular focus on nonadiabatic molecular dynamics at the ab initio level, and discuss the implications for photocatalysis from the aspects of the structure, intrinsic physical properties, formation, migration, and recombination of small polarons. Finally, various methods were proposed to advance our understanding of manipulating the small-polaron dynamics, and strategies to design high-performance TMO-based photoelectrodes were discussed.

Doubly Screened Coulomb Correction Approach for Strongly Correlated Systems

Strongly correlated systems containing d/f electrons present a challenge to conventional density functional theory such as the local density approximation or generalized gradient approximation. We developed a doubly screened Coulomb correction (DSCC) approach to perform on-site Coulomb interaction correction for strongly correlated materials. The on-site Coulomb interaction between localized d/f electrons is self-consistently determined from a model dielectric function that includes both the static dielectric and Thomas-Fermi screening. We applied DSCC to simulate the electronic and magnetic properties of typical 3d, 4f, and 5f strongly correlated systems. The accuracy of DSCC is comparable to that of hybrid functionals but an order of magnitude faster. In addition, DSCC can reflect the difference in the Coulomb interaction between metallic and insulating situations, similar to the popular but computationally expensive constrained random phase approximation approach. This feature suggests that DSCC is also a promising method for simulating Coulomb interaction parameters.

Synergistic effects of mixing and strain in high entropy spinel oxides for oxygen evolution reaction

Developing stable and efficient electrocatalysts is vital for boosting oxygen evolution reaction (OER) rates in sustainable hydrogen production. High-entropy oxides (HEOs) consist of five or more metal cations, providing opportunities to tune their catalytic properties toward high OER efficiency. This work combines theoretical and experimental studies to scrutinize the OER activity and stability for spinel-type HEOs. Density functional theory confirms that randomly mixed metal sites show thermodynamic stability, with intermediate adsorption energies displaying wider distributions due to mixing-induced equatorial strain in active metal-oxygen bonds. The rapid sol-flame method is employed to synthesize HEO, comprising five 3d-transition metal cations, which exhibits superior OER activity and durability under alkaline conditions, outperforming lower-entropy oxides, even with partial surface oxidations. The study highlights that the enhanced activity of HEO is primarily attributed to the mixing of multiple elements, leading to strain effects near the active site, as well as surface composition and coverage.

Impact of Surface Oxidation on the Morphology of Uranium Dioxide Nanoparticles

The undeniable importance of nanoparticles has led to vast efforts, in many fields of science, to understand their chemical and physical properties. In this paper, the morphology dependence of f-element nanoparticles is correlated to the oxygen environment and the type and coverage of capping ligands. This dependence was evaluated by first-principles calculations of the surface energies of different crystallographic planes (001, 110, and 111) as a function of the relative oxygen chemical potential and under the influence of different ligands. Uranium dioxide nanoparticles were the focus of this study due to their high sensitivity to oxidation compared to thorium dioxide nanoparticles, a homoleptic material but insensitive to oxidation. To fully explain the experimental observations of uranium dioxide nanocrystals, theoretical modeling shows that the consideration of surfaces with different oxidation conditions is necessary. It is shown that, for materials with low oxidation potential, such as uranium dioxide, the oxygen environment and capping ligand concentration are competing factors in determining the nanoparticle morphology.

Systematic Characterization of Electronic Metal-Support Interactions in Ceria-Supported Pt Particles

Electronic metal-support interactions affect the chemical and catalytic properties of metal particles supported on reducible metal oxides, but their characterization is challenging due to the complexity of the electronic structure of these systems. These interactions often involve different states with varying numbers and positions of strongly correlated d or f electrons and the corresponding polarons. In this work, we present an approach to characterize electronic metal-support interactions by means of computationally efficient density functional calculations within the projector augmented wave method. We describe Ce3+ cations with potentials that include a Ce4f electron in the frozen core, overcoming prevalent convergence and 4f electron localization issues. We systematically explore the stability and chemical properties of different electronic states for a Pt8/CeO2(111) model system, revealing the predominant effect of electronic metal-support interactions on Pt atoms located directly at the metal-oxide interface. Adsorption energies and the reactivity of these interface Pt atoms vary significantly upon donation of electrons to the oxide support, pointing to a strategy to selectively activate interfacial sites of metal particles supported on reducible metal oxides.

Monolayer TiSi2P4as a high-performance anode for Na-ion batteries

Exploring anode materials with overall excellent performance remains a great challenge for rechargeable Na-ion battery technologies. Herein, we have identified that monolayer TiSi2P4is just such a prospective anode candidate via first-principles calculations. It is showed to be dynamically, thermally, mechanically, and energetically stable, which provides feasibility for experimental realization. The Na diffusion on the its surface is proved to be ultrafast, with a migration energy barrier as low as 73 meV. Electronic structure confirms that the pristine system undergoes a transition from the semiconductor to metal during the whole sodiation process, which is a significant advantage to the electrode conductivity. More excitingly, monolayer TiSi2P4can accommodate up to double-sided five-layer adatoms, resulting in an ultrahigh theoretical capacity of 1176 mA h g-1and a low average open-circuit voltage of 0.195 V. Moreover, the maximally sodiated electrode monolayer yields rather small in-plane lattice expansion of only 1.40%, which ensures reversible deformation and excellent cycling stability as further corroborated by structural relaxation andab initiomolecular dynamics simulation. Overall, all of these results point to the potential that monolayer TiSi2P4can serve as a promising anode candidate for application in high-performance low-cost Na-ion batteries.

Potential-Dependent Active Moiety of Fe-N-C Catalysts for the Oxygen Reduction Reaction

The real active moiety of Fe-N-C single-atom catalysts (SACs) during the oxygen reduction reaction (ORR) depends on the applied potential. Here, we examine the ORR activity of various SAC active moieties (Fe-N₄, Fe-(OH)N₄, Fe-(O₂)N₄, and Fe-(OH₂)N₄) over a wide potential window ranging from -0.8 to 1.0 V (vs. SHE) using constant potential density functional theory calculations. We show that the ORR activity of the Fe-N₄ moiety is hindered by the slow *OH protonation, while the Fe-(OH₂)N₄ (0.4 V ≤ U ≤ 1.0 V), *O₂-assisted Fe-N₄ (-0.6 V ≤ U ≤ 0.2 V), and Fe-(OH)N₄ (U = -0.8 V) moieties dominate the ORR activity of the Fe-N-C catalysts at different potential windows. These oxygenated species modified the single-atom Fe sites and can promote *OH protonation by regulating the electron occupancy of the Fe 3d_z² (spin-up) and Fe 3d_xz (spin-down) orbitals. Overall, our findings provide guidance for understanding the active moieties of SACs.

Tuning the carrier injection barrier of hybrid metal-organic interfaces on rare earth-gold surface compounds

Magnetic hybrid metal-organic interfaces possess a great potential in areas such as organic spintronics and quantum information processing. However, tuning their carrier injection barriers on-demand is fundamental for the implementation in technological devices. We have prepared hybrid metal-organic interfaces by the adsorption of copper phthalocyanine CuPc on REAu₂ surfaces (RE = Gd, Ho and Yb) and studied their growth, electrostatics and electronic structure. CuPc exhibits a long-range commensurability and a vacuum level pinning of the molecular energy levels. We observe a significant effect of the RE valence of the substrate on the carrier injection barrier of the hybrid metal-organic interface. CuPc adsorbed on trivalent RE-based surfaces (HoAu₂ and GdAu₂) exhibits molecular level energies that may allow injection carriers significantly closer to an ambipolar injection behavior than in the divalent case (YbAu₂).

Room-temperature valence transition in a strain-tuned perovskite oxide

Cobalt oxides have long been understood to display intriguing phenomena known as spin-state crossovers, where the cobalt ion spin changes vs. temperature, pressure, etc. A very different situation was recently uncovered in praseodymium-containing cobalt oxides, where a first-order coupled spin-state/structural/metal-insulator transition occurs, driven by a remarkable praseodymium valence transition. Such valence transitions, particularly when triggering spin-state and metal-insulator transitions, offer highly appealing functionality, but have thus far been confined to cryogenic temperatures in bulk materials (e.g., 90 K in Pr_1-xCa_xCoO₃). Here, we show that in thin films of the complex perovskite (Pr_1-yY_y)_1-xCa_xCoO_3-δ, heteroepitaxial strain tuning enables stabilization of valence-driven spin-state/structural/metal-insulator transitions to at least 291 K, i.e., around room temperature. The technological implications of this result are accompanied by fundamental prospects, as complete strain control of the electronic ground state is demonstrated, from ferromagnetic metal under tension to nonmagnetic insulator under compression, thereby exposing a potential novel quantum critical point.

Role of the Magnetic Anisotropy in Atomic-Spin Sensing of 1D Molecular Chains

One-dimensional metal-organic chains often possess a complex magnetic structure susceptible to modification by alteration of their chemical composition. The possibility to tune their magnetic properties provides an interesting playground to explore quasi-particle interactions in low-dimensional systems. Despite the great effort invested so far, a detailed understanding of the interactions governing the electronic and magnetic properties of the low-dimensional systems is still incomplete. One of the reasons is the limited ability to characterize their magnetic properties at the atomic scale. Here, we provide a comprehensive study of the magnetic properties of metal-organic one-dimensional (1D) coordination polymers consisting of 2,5-diamino-1,4-benzoquinonediimine ligands coordinated with Co or Cr atoms synthesized under ultrahigh-vacuum conditions on a Au(111) surface. A combination of integral X-ray spectroscopy with local-probe inelastic electron tunneling spectroscopy corroborated by multiplet analysis, density functional theory, and inelastic electron tunneling simulations enables us to obtain essential information about their magnetic structures, including the spin magnitude and orientation at the magnetic atoms, as well as the magnetic anisotropy.

Angular Jahn-Teller Effect and Photoluminescence of the Tetrahedral Coordinated Mn2+ Activators in Solids-A First-Principles Study

First-principles calculations based on density functional theory have been performed to investigate the electronic structure, excited-state Jahn-Teller distortion, and photoluminescence of the multielectron d⁵ system of the strongly covalent tetrahedral coordinated Mn²⁺ activator in solids. The electronic structure of the ⁴T₁ and ⁴A₁/⁴E excited states is analyzed, and Slater's transition-state method and occupation matrix control methodology are applied to deal with the spin contamination in the lower-spin excited states, which is due to the mixing of the ground state of the same spin projection number. In a series of covalent tetrahedral coordinations, the ⁶A₁ → ⁴T₁ and ⁴A₁/⁴E excitations and the ⁴T₁ → ⁶A₁ emission energies are obtained and compared to the reported experimental results. The nephelauxetic effect follows O^2- < S^2- ≈ Se^2- < N^3-, and the larger nephelauxetic effect and crystal field strength lead to the red-shifted emission of nitride phosphors. The Jahn-Teller distortion of the ⁴T₁ states is dominated by the e-type angular distortion of the [MnL₄] moiety (L being the ligand), which accounts for the small Stokes shift of tetrahedral coordinated Mn²⁺. The results show that the ground- and excited-state electronic and geometric structures and the luminescent property of tetrahedral coordinated Mn²⁺ can be reliably predicted. The method can be further explored to interpret and discriminate the luminescent properties of materials containing a variety of different Mn²⁺ sites and complexes and even other transition metals.

DFT + U Study of Uranium Dioxide and Plutonium Dioxide with Occupation Matrix Control

DFT + U with occupation matrix control (OMC) is applied to study computationally bulk UO₂ and PuO₂, the latter for the first time. Using the PBESol functional in conjunction with OMC locates AFM and NM ground states for UO₂ and PuO₂, respectively, in agreement with experimental findings. By simulating the lattice parameter, magnetic moment, band gap, and densities of states, U = 4.0 eV is recommended for AFM UO₂, yielding data close to experiments for all considered properties. U = 4.5 and 4.0 eV are recommended for NM and AFM PuO₂, respectively, though much larger U values (c. 10 eV) are required to yield the most recently reported PuO₂ band gap. For both oxides, several excited states have similar properties to the ground state, reinforcing the need to employ OMC wherever possible.

Role of Polarons in Single-Atom Catalysts: Case Study of Me1 [Au1, Pt1, and Rh1] on TiO2(110)

The local environment of metal-oxide supported single-atom catalysts plays a decisive role in the surface reactivity and related catalytic properties. The study of such systems is complicated by the presence of point defects on the surface, which are often associated with the localization of excess charge in the form of polarons. This can affect the stability, the electronic configuration, and the local geometry of the adsorbed adatoms. In this work, through the use of density functional theory and surface-sensitive experiments, we study the adsorption of Rh1, Pt1, and Au1 metals on the reduced TiO2(110) surface, a prototypical polaronic material. A systematic analysis of the adsorption configurations and oxidation states of the adsorbed metals reveals different types of couplings between adsorbates and polarons. As confirmed by scanning tunneling microscopy measurements, the favored Pt1 and Au1 adsorption at oxygen vacancy sites is associated with a strong electronic charge transfer from polaronic states to adatom orbitals, which results in a reduction of the adsorbed metal. In contrast, the Rh1 adatoms interact weakly with the excess charge, which leaves the polarons largely unaffected. Our results show that an accurate understanding of the properties of single-atom catalysts on oxide surfaces requires a careful account of the interplay between adatoms, vacancy sites, and polarons.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11244-022-01651-0.

The Role of Local DFT+U Minima in the First-Principles Modeling of the Metal-Insulator Transition in Vanadium Dioxide

The DFT+U method is frequently employed to improve the first-principles description of strongly correlated materials. However, it is prone to deliver metastable electronic minima. While these local minima of the DFT+U method are often considered to be computational artifacts, their physical meaning and relationship to true excited states remains unclear. In this work, the possibility of theoretically modeling transformations in the solid state that require thermal or optical excitations of electrons is explored, taking into account the metastable states of the computationally undemanding DFT+U formalism. For this purpose, we choose to examine the example of the VO₂ metal-insulator transition. Metastable states that are located on different electronic potential energy surfaces are found to correspond to experimentally observed VO₂ phases. The identified metastable electronic states can be used to model the collapse of the VO₂ band gap at elevated temperatures and upon photoexcitation as well as other monoclinic-monoclinic phase transformations. The results suggest that local DFT+U minima can indeed carry physical meaning, while they remain under-reported in theoretical literature on transition metal oxides like VO₂.

Oxygen vacancy distributions and electron localization in a CeO2(100) nanocube

Oxygen vacancy distributions in a 5 nm CeO2 nanocube were determined using the Reverse Monte Carlo method. The oxygen vacancies tend to be located on the surface of the CeO2 nanocube, with far fewer in subsurface and internal regions.

Exploring metastable states in UO2using hybrid functionals and dynamical mean field theory

A detailed exploration of thef-atomic orbital occupancy space for UO₂is performed using a first principles approach based on density functional theory (DFT), employing a full hybrid functional within a systematic basis set. Specifically, the PBE0 functional is combined with an occupancy biasing scheme implemented in a wavelet-based algorithm which is adapted to large supercells. The results are compared with previous DFT +Ucalculations reported in the literature, while dynamical mean field theory is also performed to provide a further base for comparison. This work shows that the computational complexity of the energy landscape of a correlatedf-electron oxide is much richer than has previously been demonstrated. The resulting calculations provide evidence of the existence of multiple previously unexplored metastable electronic states of UO₂, including those with energies which are lower than previously reported ground states.

Oxygen chemistry of halogen-doped CeO2(111)

We study substitutional fluorine, chlorine and bromine impurities at CeO₂(111), and their effects on the oxygen chemistry of the surface, using density functional theory. We find that impurity formation results in a halide ion and one Ce³⁺ ion for all three halogens, although the formation energy depends strongly on the identity of the halogen; however, once formed, all three halogens exhibit a similar propensity to form impurity-impurity pairs. Furthermore, while the effects of halogen impurities on oxygen vacancy formation are marginal, they are more significant for oxygen molecule adsorption, due to electron transfer from the Ce³⁺ ion which results in an adsorbed superoxide molecule. We also consider the displacement of a halide ion on to the surface by half of an oxygen molecule, and find that the energy required to do so depends strongly not only on the identity of the halogen, but also on whether or not a second halogen impurity, with its associated Ce³⁺ ion, is present; if it is, then the process is greatly facilitated. Overall, our results demonstrate the existence of a rich variety of ways in which the oxygen chemistry of CeO₂(111) may be modified by the presence of halogen dopants.

Interplay between London Dispersion, Hubbard U, and Metastable States for Uranium Compounds

High-throughput computational studies of lanthanide and actinide chemistry with density-functional theory are complicated by the need for Hubbard U corrections, which ensure localization of the f-electrons, but can lead to metastable states. This work presents a systematic investigation of the effects of both Hubbard U value and metastable states on the predicted structural and thermodynamic properties of four uranium compounds central to the field of nuclear fuels: UC, UN, UO₂, and UCl₃. We also assess the impact of the exchange-hole dipole moment (XDM) dispersion correction on the computed properties. Overall, the choice of Hubbard U value and inclusion of a dispersion correction cause larger variations in the computed geometric properties than result from metastable states. The weak dependence of structure optimization on metastable states should simplify future high-throughput calculations on actinides. Conversely, addition of the dispersion correction is found to offset the repulsion introduced by the Hubbard U term and provides greatly improved agreement with experiment for both cell volumes and heats of formation. The XDM dispersion correction is largely invariant to the chosen U value, making it a robust dispersion correction for actinide systems.

Computational studies of the electronic structure of copper-doped ZnO quantum dots

Copper-doped ZnO quantum dots (QDs) have attracted substantial interest. The electronic structure and optical and magnetic properties of Cu³⁺(d⁸)-, Cu²⁺(d⁹)-, and Cu⁺(d¹⁰)-doped ZnO QDs with sizes up to 1.5 nm are investigated using the GGA+U approximation, with the +U corrections applied to d (Zn), p(O), and d(Cu) orbitals. Taking +Us parameters, as optimized in previous bulk calculations, we obtain the correct band structure of ZnO QDs. Both the description of electronic structure and thermodynamic charge state transitions of Cu in ZnO QDs agree with the results of bulk calculations due to the strong localization of Cu defect energy levels. Atomic displacements around Cu are induced by strong Jahn-Teller distortion and affect Kohn-Sham energies and thermodynamic transition levels. The average bond length of Cu-O and the defect structure are crucial factors influencing the electronic properties of Cu in ZnO QDs. The analysis of the optical properties of Cu in ZnO QDs is reported. The GGA+U results, compared with the available experimental data, support Dingle's model [Phys. Rev. Lett. 23, 579 (1969)], in which the structured green luminescence observed in bulk and nanocrystals originates from the [(Cu⁺, hole) → Cu²⁺] transition. We also examine the magnetic interaction between the copper pair for two charge states: 0 and +2, and four positions relative to the center of QDs. Ferromagnetic interaction between ions is obtained for every investigated configuration. The magnitude of ferromagnetism increases for positive charge defects due to the strong hybridization of the d(Cu) and p(O) states.

CO2 adsorption on the pristine and reduced CeO2 (111) surface: Geometries and vibrational spectra by first principles simulations

First principles simulations of carbon dioxide adsorbed on the ceria (CeO₂) (111) surface are discussed in terms of structural features, stability, charge transfer, and vibrational modes. For this purpose, different density functional theory methods, such as Perdew-Burke-Ernzerhof (PBE) PBE and Hubbard correction, hybrid functionals, and different basis sets have been applied and compared. Both the stoichiometric and the reduced (111) surfaces are considered, where the electronic structure of the latter is obtained by introducing oxygen vacancies on the topmost or the subsurface oxygen layer. Both the potential energy surfaces of the reduced ceria surface and the adsorbate-surface complex are characterized by numerous local minima, of which the relative stability depends strongly on the electronic structure method of choice. Bent CO₂ configurations in close vicinity to the surface oxygen vacancy that partially re-oxidize the reduced ceria surface have been identified as the most probable stable minima. However, the oxygen vacancy concentration on the surface turns out to have a direct impact on the relative stability of possible adsorption configurations. Finally, the vibrational analyses of selected adsorbed species on both the stoichiometric and reduced surfaces show promising agreement with previous theoretical and experimental results.

Eu-Mn Charge Transfer and the Strong Charge-Spin-Electronic Coupling Behavior in EuMnO3

Based on first-principles calculations with the DFT + U method, the couplings of lattice, charge, spin, and electronic behaviors underlying the Eu-Mn charge transfer in a strongly correlated system of EuMnO₃ were investigated. The potential valence transition from Eu³⁺/Mn³⁺ to Eu²⁺/Mn⁴⁺ was observed in a compressed lattice with little distortions, which is achieved under hydrostatic pressure and external strain. The intraplane antiferromagnetism (AFM) of Mn is proved to be instrumental in the emergence of Eu²⁺. Furthermore, we calculated the magnetic exchange interactions within two equilibrium structures of Eu³⁺Mn³⁺O₃ and Eu²⁺Mn⁴⁺O₃. Mn-Mn ferromagnetic exchange in the ab-plane is enhanced strongly in the Eu²⁺Mn⁴⁺O₃ structure, contributing to the existence of mixed states. The versatile electronic structures were obtained within the Eu²⁺Mn⁴⁺O₃ phase by imposing different magnetic configurations on the Eu and Mn sublattice, attributed to the coupling of charge transfer and magnetic orderings. It is found that the intraplane ferromagnetic ordering of Mn leads to a metallic electronic structure with the coexistence of Eu²⁺ and Eu³⁺, while the intraplane AFM Mn spin ordering leads to insulating states only with Eu²⁺. Notably, a half-metallic characteristic emerges at the magnetic ground state of CF ordering (C-type AFM for the Eu sublattice and ferromagnetic for the Mn sublattice), which makes such a supposed phase more intriguing than the conventional experimental phase. Additionally, the mixture of delocalized 4f with 5d states of Eu in the background of Mn 3d and O 2p orbitals implies a pathway of Eu 4f 5d ↔ O 2p ↔ Mn 3d for charge transfer between Eu and Mn. Our calculation shows that the Eu-Mn charge transfer could be expected in compressed EuMnO₃ and the introduction of Eu²⁺ 4f states near the Fermi level plays an important role in manipulating the interlinks of charge and spin together with electronic behaviors.

Formation and stability of small polarons at the lithium-terminated Li4Ti5O12 (LTO) (111) surface

Zero strain insertion, high cycling stability, and a stable charge/discharge plateau are promising properties rendering Lithium Titanium Oxide (LTO) a possible candidate for an anode material in solid state Li ion batteries. However, the use of pristine LTO in batteries is rather limited due to its electronically insulating nature. In contrast, reduced LTO shows an electronic conductivity several orders of magnitude higher. Studying bulk reduced LTO, we could show recently that the formation of polaronic states can play a major role in explaining this improved conductivity. In this work, we extend our study toward the lithium-terminated LTO (111) surface. We investigate the formation of polarons by applying Hubbard-corrected density functional theory. Analyzing their relative stabilities reveals that positions with Li ions close by have the highest stability among the different localization patterns.

Water on Actinide Dioxide Surfaces: A Review of Recent Progress

The fluorite structured actinide dioxides (AnO2), especially UO2, are the most common nuclear fuel materials. A comprehensive understanding of their surface chemistry is critical because of its relevance to the safe handling, usage, and storage of nuclear fuels. Because of the ubiquitous nature of water (H2O), its interaction with AnO2 has attracted significant attention for its significance in studies of nuclear fuels corrosion and the long-term storage of nuclear wastes. The last few years have seen extensive experimental and theoretical studies on the H2O–AnO2 interaction. Herein, we present a brief review of recent advances in this area. We focus on the atomic structures of AnO2 surfaces, the surface energies, surface oxygen vacancies, their influence on the oxidation states of actinide atoms, and the adsorption and reactions of H2O on stoichiometric and reduced AnO2 surfaces. Finally, a summary and outlook of future studies on surface chemistry of AnO2 are given. We intend for this review to encourage broader interests and further studies on AnO2 surfaces.

Mobile Small Polarons Qualitatively Explain Conductivity in Lithium Titanium Oxide Battery Electrodes

Lithium titanium oxide Li₄Ti₅O₁₂ is an intriguing anode material promising particularly long-life batteries, due to its remarkable phase stability during (dis)charging of the cell. However, its usage is limited by its low intrinsic electronic conductivity. Introducing oxygen vacancies can be one method for overcoming this drawback, possibly by altering the charge carrier transport mechanism. We use Hubbard corrected density functional theory to show that polaronic states in combination with a possible hopping mechanism can play a crucial role in the experimentally observed increase in electronic conductivity. To gauge polaronic charge mobility, we compute the relative stabilities of different localization patterns and estimate polaron hopping barrier heights.

Carrier localization in perovskite nickelates from oxygen vacancies

Point defects, such as oxygen vacancies, control the physical properties of complex oxides, relevant in active areas of research from superconductivity to resistive memory to catalysis. In most oxide semiconductors, electrons that are associated with oxygen vacancies occupy the conduction band, leading to an increase in the electrical conductivity. Here we demonstrate, in contrast, that in the correlated-electron perovskite rare-earth nickelates, RNiO₃ (R is a rare-earth element such as Sm or Nd), electrons associated with oxygen vacancies strongly localize, leading to a dramatic decrease in the electrical conductivity by several orders of magnitude. This unusual behavior is found to stem from the combination of crystal field splitting and filling-controlled Mott-Hubbard electron-electron correlations in the Ni 3d orbitals. Furthermore, we show the distribution of oxygen vacancies in NdNiO₃ can be controlled via an electric field, leading to analog resistance switching behavior. This study demonstrates the potential of nickelates as testbeds to better understand emergent physics in oxide heterostructures as well as candidate systems in the emerging fields of artificial intelligence.

Optimizing the orbital occupation in the multiple minima problem of magnetic materials from the metaheuristic firefly algorithm

We present the use and implementation of the firefly algorithm to help in scanning the multiple metastable minima of orbital occupations in density functional theory (DFT) plus Hubbard U correction and to identify the ground state occupations in strongly correlated materials. We show the application of this implementation with the Abinit code on KCoF3 and UO2 crystals, which are typical d and f electron systems with numerous occupation minima. We demonstrate the validity and performance of the method by comparing with previous methodologies. The method is general and can be applied to any code using constrained occupation matrices.

The Structure of Oxygen Vacancies in the Near-Surface of Reduced CeO2 (111) Under Strain

Strain has been widely recognized as important for tuning the behavior of defects in metal oxides since properties such as defect configuration, electronic structure, excess charge localization, and local atomic distortions may be affected by surface strain. In CeO₂, the most widely used promoter in three-way catalysts and solid state electrolyte in fuel cells, the behaviors of oxygen vacancies, and associated Ce³⁺ polarons are crucial in applications. Recent STM and AFM investigations as well as DFT-based calculations have indicated that in the near-surface of CeO₂ (111), at low temperatures and vacancy concentrations, subsurface oxygen vacancies are more stable than surface ones, and the Ce³⁺ ions are next-nearest neighbors to both types of vacancies, which can be explained by the better ability of the system to relax the lattice strain induced by vacancy formation as well as by the excess charge localization. The results also revealed that the interaction between first-neighbor vacancies is repulsive. In this work, the relative stability of surface and subsurface oxygen vacancies at the CeO₂ (111) surface under in-plane strain is investigated by means of DFT+U calculations. The tensile strain favors isolated surface vacancies with next nearest neighbor polarons, whereas isolated subsurface vacancies with nearest neighbor polarons are energetically favored under compressive strain. In addition, the formation of both surface and subsurface dimers is favored over having corresponding isolated species under compressive strain, which implies the possibility of controlling the formation of vacancy clusters using strain. In many applications, ceria is employed as a supported thin film or within a heterostructure in which ceria can be strained, and this study shows that strain can be a useful handle to tune properties of such materials.

Origin of band gaps in 3d perovskite oxides

With their broad range of properties, ABO₃ transition metal perovskite oxides have long served as a platform for device applications and as a testing bed for different condensed matter theories. Their insulating character and structural distortions are often ascribed to dynamical electronic correlations within a universal, symmetry-conserving paradigm. This view restricts predictive theory to complex computational schemes, going beyond density functional theory (DFT). Here, we show that, if one allows symmetry-breaking energy-lowering crystal symmetry reductions and electronic instabilities within DFT, one successfully and systematically recovers the trends in the observed band gaps, magnetic moments, type of magnetic and crystallographic ground state, bond disproportionation and ligand hole effects, Mott vs. charge transfer insulator behaviors, and the amplitude of structural deformation modes including Jahn-Teller in low temperature spin-ordered and high temperature disordered paramagnetic phases. We then provide a classification of the four mechanisms of gap formation and establish DFT as a reliable base platform to study the ground state properties in complex oxides.

It is often stated that first principles studies of transition metal oxides require dynamically correlated methods to correctly produce gap formation, magnetism and structural distortions. Varignon et al. show instead that static correlations are sufficient to capture these features in the ABO₃ oxide series.

Biaxial strain effect on the electronic structure and valleytronic properties of a MoS2/CoO(111) heterostructure

Spin splitting, valley splitting and Berry curvature at the K and K′ valleys of a MoS₂/CoO(111) heterostructure can be tuned continually by biaxial tensile strain.

Metastable electronic states in uranium tetrafluoride

The DFT+U approach, where U is the Hubbard-like on-site Coulomb interaction, has successfully been used to improve the description of transition metal oxides and other highly correlated systems, including actinides. The secret of the DFT+U approach is the breaking of d or f shell orbital degeneracy and adding an additional energetic penalty to non-integer occupation of orbitals. A prototypical test case, UO2, benefits from the +U approach whereby the bare LDA method predicts UO2 to be a ferromagnetic metal, whereas LDA+U correctly predicts UO2 to be insulating. However, the concavity of the energetic penalty in the DFT+U approach can lead to a number of convergent "metastable" electronic configurations residing above the ground state. Uranium tetrafluoride (UF4) represents a more complex analogy to UO2 in that the crystal field has lower symmetry and the unit cell contains two symmetrically distinct U atoms. We explore the metastable states in UF4 using several different methods of selecting initial orbital occupations. Two methods, a "pre-relaxation" method wherein an initial set of orbital eigenvectors is selected via the self-consistency procedure and a crystal rotation method wherein the x, y, z axes are brought into alignment with the crystal field, are explored. We show that in the case of UF4, which has non-collinearity between its crystal axes and the U atoms' crystal field potentials, the orbital occupation matrices are much more complex and should be analyzed using a novel approach. In addition to demonstrating a complex landscape of metastable electronic states, UF4 also shows significant hybridization in U-F bonding, which involves non-trivial contributions from s, p, d, and f orbitals.

Assessment of van der Waals inclusive density functional theory methods for layered electroactive materials

This work provides solid guidance for the selection of accurate and robust vdW-inclusive methods for high-throughput computational screening of layered electroactive materials.

Modification of Charge Trapping at Particle/Particle Interfaces by Electrochemical Hydrogen Doping of Nanocrystalline TiO2

Particle/particle interfaces play a crucial role in the functionality and performance of nanocrystalline materials such as mesoporous metal oxide electrodes. Defects at these interfaces are known to impede charge separation via slow-down of transport and increase of charge recombination, but can be passivated via electrochemical doping (i.e., incorporation of electron/proton pairs), leading to transient but large enhancement of photoelectrode performance. Although this process is technologically very relevant, it is still poorly understood. Here we report on the electrochemical characterization and the theoretical modeling of electron traps in nanocrystalline rutile TiO₂ films. Significant changes in the electrochemical response of porous films consisting of a random network of TiO₂ particles are observed upon the electrochemical accumulation of electron/proton pairs. The reversible shift of a capacitive peak in the voltammetric profile of the electrode is assigned to an energetic modification of trap states at particle/particle interfaces. This hypothesis is supported by first-principles theoretical calculations on a TiO₂ grain boundary, providing a simple model for particle/particle interfaces. In particular, it is shown how protons readily segregate to the grain boundary (being up to 0.6 eV more stable than in the TiO₂ bulk), modifying its structure and electron-trapping properties. The presence of hydrogen at the grain boundary increases the average depth of traps while at the same time reducing their number compared to the undoped situation. This provides an explanation for the transient enhancement of the photoelectrocatalytic activity toward methanol photooxidation which is observed following electrochemical hydrogen doping of rutile TiO₂ nanoparticle electrodes.

TiO2 Band Restructuring by B and P Dopants

An examination of the effect of B- and P-doping and codoping on the electronic structure of anatase TiO2 by performing density functional theory calculations revealed the following: (i) B- or P-doping effects are similar to atomic undercoordination effects on local bond relaxation and core electron entrapment; (ii) the locally entrapped charge adds impurity levels within the band gap that could enhance the utilization of TiO2 to absorb visible light and prolong the carrier lifetime; (iii) the core electron entrapment polarizes nonbonding electrons in the upper edges of the valence and conduction bands, which reduces not only the work function but also the band gap; and (iv) work function reduction enhances the reactivity of the carriers and band gap reduction promotes visible-light absorption. These observations may shed light on effective catalyst design and synthesis.