Accurate first-principles structures and energies of diversely bonded systems from an efficient density functional

Reworking the Tao-Mo Exchange-Correlation Functional. III. Improved Deorbitalization Strategy and Faithful Deorbitalization

We present a deorbitalization of the recent simplified, regularized Tao-Mo exchange functional ( J. Chem. Phys. 2023, 159, 214102) that is faithful to the parent functional. That is a major gain relative to our earlier deorbitalization which did poorly on molecular heats of formation ( J. Chem. Phys. 2023, 159, 214103). The improvement arises from augmentation of the Mejía-Rodríguez and Trickey deorbitalization strategy ( Phys. Rev. A 2017, 96, 052512) to use a smoothed replacement for the reduced density Laplacian (conventionally denoted q) obtained from that Laplacian itself. The augmentation also rationalizes the improvement obtained from the cutoff of q < 0 that was poorly understood at the time of the previous paper. The new scheme yields deorbitalized chemical region indicators that are much closer to those from the parent, orbital-dependent functional than were obtainable from the previous deorbitalization. It also replicates the good 3d elemental magnetization of the parent functional reasonably well.

Mn2C MXene functionalized by oxygen is a semiconducting antiferromagnet and an efficient visible light absorber

Manganese-based MXenes are promising two-dimensional materials due to the broad palette of their magnetic phases and the possibility of experimental preparation because the corresponding MAX phase was already prepared. Here, we systematically investigated geometrical conformers and spin solutions of oxygen-terminated Mn2C MXene and performed subsequent many-body calculations to obtain reliable electronic and optical properties. Allowing energy-lowering using the correct spin ordering via supercell magnetic motifs is essential for the Mn2CO2 system. The stable ground-state Mn2CO2 conformation is antiferromagnetic (AFM) with zigzag lines of up and down spins on Mn atoms. The AFM nature is consistent with the parent MAX phase and even the clean depleted Mn2C sheet. Other magnetic states and geometrical conformations are energetically very close, providing state-switching possibilities in the material. Subsequent many-body GW and Bethe-Salpeter equation (BSE) calculations provide indirect semiconductor characteristics of AFM Mn2CO2 with a fundamental gap of 2.1 eV (and a direct gap of 2.4 eV), the first bright optical transition at 1.3 eV and extremely strongly bound (1.1 eV) first bright exciton. Mn2CO2 absorbs efficiently the whole visible light range and near ultraviolet range (between 10 and 20%).

Universal scales of electronegativity and ionicity from electron scattering factors

Electronegativity and ionicity are important but difficult concepts. In this work, we present the universal scale of electronegativity by utilizing electron scattering factors. A mapping of our proposed "scattering electronegativity" with Pauling electronegativity and other electronegativity scales is given. Two new ionicity scales based on electron scattering factors are also proposed and compared with Philips's ionicity and Pauling's ionicity. Although electron scattering factors are based on neutral, free, atomic features, scattering electronegativity and ionicity are able to provide rich charge or bonding information of materials in molecules or solid states. Our newly proposed electronegativity scale satisfies the metalloid band criterion (or Si rule). A new Mg rule and a new concept of critical ionicity gap are proposed, successfully distinguishing between the four-fold and six-fold coordinates in the ANB8-N crystals.

Unveiling unconventional CH4-Xe compounds and their thermodynamic properties at extreme conditions

Inert gases (e.g., He and Xe) can exhibit chemical activity at high pressure, reacting with other substances to form compounds of unexpected chemical stoichiometry. This work combines first-principles calculations and crystal structure predictions to propose four unexpected stable compounds of CH4Xe3, (CH4)2Xe, (CH4)3Xe, and (CH4)3Xe2 at pressure ranges from 2 to 100 GPa. All structures are composed of isolated Xe atoms and CH4 molecules except for (CH4)3Xe2, which comprises a polymerization product, C3H8, and hydrogen molecules. Ab initio molecular dynamics simulations indicate that pressure plays a very important role in the different temperature driving state transitions of CH4-Xe compounds. At lower pressures, the compounds follow the state transition of solid-plastic-fluid phases with increasing temperature, while at higher pressures, the stronger Xe-C interaction induces the emergence of a superionic state for CH4Xe3 and (CH4)3Xe2 as temperature increases. These results not only expand the family of CH4-Xe compounds, they also contribute to models of the structures and evolution of planetary interiors.

Hyperfine interactions for small systems including transition-metal elements using self-interaction corrected density-functional theory

The interactions between the electronic magnetic moment and the nuclear spin moment, i.e., magnetic hyperfine (HF) interactions, play an important role in understanding electronic properties of magnetic systems and in realizing platforms for quantum information science applications. We investigate the HF interactions for atomic systems and small molecules, including Ti or Mn, by using Fermi-Löwdin orbital (FLO) based self-interaction corrected (SIC) density-functional theory. We calculate the Fermi contact (FC) and spin-dipole terms for the systems within the local density approximation (LDA) in the FLO-SIC method and compare them with the corresponding values without SIC within the LDA and generalized-gradient approximation (GGA), as well as experimental data. For the moderately heavy atomic systems (atomic number Z ≤ 25), we find that the mean absolute error of the FLO-SIC FC term is about 27 MHz (percentage error is 6.4%), while that of the LDA and GGA results is almost double that. Therefore, in this case, the FLO-SIC results are in better agreement with the experimental data. For the non-transition-metal molecules, the FLO-SIC FC term has the mean absolute error of 68 MHz, which is comparable to both the LDA and GGA results without SIC. For the seven transition-metal-based molecules, the FLO-SIC mean absolute error is 59 MHz, whereas the corresponding LDA and GGA errors are 101 and 82 MHz, respectively. Therefore, for the transition-metal-based molecules, the FLO-SIC FC term agrees better with experiment than the LDA and GGA results. We observe that the FC term from the FLO-SIC calculation is not necessarily larger than that from the LDA or GGA for all the considered systems due to the core spin polarization, in contrast to the expectation that SIC would increase the spin density near atomic nuclei, leading to larger FC terms.

Nonlinear Effects of Hydrophobic Confinement on the Electronic Structure and Dielectric Response of Water

Fundamental studies of the dielectrics of confined water are critical to understand the ion transport across biological and synthetic nanochannels. The relevance of these fundamental studies, however, surmounts the difficulty of probing water's dielectric constant as a function of a fine variation in confinement. In this work, we explore the computational efficiency of machine learning potentials to derive the confinement effects on the dielectric constant, polarization, and dipole moment of water. Our simulations predict an enhancement of the axial dielectric constant of water under extreme confinement, arising from either the formation of ferroelectric structures of ordered water or larger dipole fluctuations facilitated by the disruption of water's H-bond network. Our study highlights the impact of hydrophobic nanoconfinement on the dielectric constant and on the ionic and electronic structure of water molecules, pointing to the importance of geometric flexibility and electronic polarizability to properly model confinement effects on water.

Exploring unconventional ferromagnetism in hole-doped LaCrAsO: insights into charge-transfer and magnetic interactions

Itinerant ferromagnetism due to the canonical double exchange (CDE) mechanism always occurs at low doping concentrations. Here we demonstrate the occurrence of robust itinerant ferromagnetism that can persist high doping concentrations. Using experimentally synthesized LaCrAsO as an illustrative example, we study the effects of hole doping via first-principles calculations and observe that the parent G-type antiferromagnetism vanishes quickly at a low doping concentration (∼0.20) and the system becomes a ferromagnetic metal due to the CDE mechanism. As the doping concentration continues to increase, the As 4p orbitals are gradually pushed up to the Fermi level and doped with holes. These ligand holes participate in the exchange interactions and drive the system toward ferromagnetism. Therefore, itinerant ferromagnetism doesn't terminate at an intermediate doping concentration as the CDE mechanism usually predicts. Furthermore, our results reveal that both the nearest and the next-nearest ferromagnetic exchange coupling strengths keep growing with doping concentration monotonically, showing that the emergent ferromagnetism mediated by As 4p orbitals is "stronger" than that of the CDE picture. Our work unlocks a new mechanism of itinerant ferromagnetism and potentially paves the way towards novel magneto-transport properties.

Spin Currents via the Gauge Principle for Meta-Generalized Gradient Exchange-Correlation Functionals

The prominence of density functional theory in the field of electronic structure computation stems from its ability to usefully balance accuracy and computational effort. At the base of this ability is a functional of the electron density: the exchange-correlation energy. This functional satisfies known exact conditions that guide the derivation of approximations. The strongly constrained and appropriately normed (SCAN) approximation stands out as a successful, modern, example. In this Letter, we demonstrate how the SU(2) gauge invariance of the exchange-correlation functional in spin current density functional theory allows us to add an explicit dependence on spin currents in the SCAN functional (here called JSCAN)-and similar meta-generalized-gradient functional approximations-solely invoking first principles. In passing, a spin-current dependent generalization of the electron localization function (here called JELF) is also derived. The extended forms are implemented in a developer's version of the crystal23 program. Applications on molecules and materials confirm the practical relevance of the extensions.

Emerging Conformational-Analysis Protocols from the RTCONF55-16K Reaction Thermochemistry Conformational Benchmark Set

RTCONF55-16K is a new, reactive conformational data set based on cost-efficient methods to assess different conformational analysis protocols. Our reference calculations underpinned the accuracy of the CENSO (Grimme et al. J. Phys. Chem. A, 2021, 125, 4039) procedure and resulted in alternative recipes with different cost-accuracy compromises. Our general-purpose and economical protocols (CENSO-light and zero, respectively) were found to be 10-30 times faster than the original algorithm, adding only 0.4-0.7 kcal/mol absolute error to the relative free energy estimates.

Confinement Effects on Proton Transfer in TiO2 Nanopores from Machine Learning Potential Molecular Dynamics Simulations

Improved understanding of proton transfer in nanopores is critical for a wide range of emerging applications, yet experimentally probing mechanisms and energetics of this process remains a significant challenge. To help reveal details of this process, we developed and applied a machine learning potential derived from first-principles calculations to examine water reactivity and proton transfer in TiO2 slit-pores. We find that confinement of water within pores smaller than 0.5 nm imposes strong and complex effects on water reactivity and proton transfer. Although the proton transfer mechanism is similar to that at a TiO2 interface with bulk water, confinement reduces the activation energy of this process, leading to more frequent proton transfer events. This enhanced proton transfer stems from the contraction of oxygen-oxygen distances dictated by the interplay between confinement and hydrophilic interactions. Our simulations also highlight the importance of the surface topology, where faster proton transport is found in the direction where a unique arrangement of surface oxygens enables the formation of an ordered water chain. In a broader context, our study demonstrates that proton transfer in hydrophilic nanopores can be enhanced by controlling pore size, surface chemistry, and topology.

Disorder induced band gap lowering in kesterite type Cu2ZnSnSe4and Ag2ZnSnSe4: a first-principles and special quasirandom structures investigation

Quaternary chalcogenides, i.e. Cu2ZnSnS4, crystallising in the kesterite crystal structure have already been demonstrated as potential building blocks of thin film solar cells, containing only abundant elements and exhibiting power conversion efficiencies of about 14.9% so far. However, due to the potential presence of several structurally similar polymorphs, the unequivocal identification of their ground state crystal structures required the application of more elaborate neutron diffraction experiments. One particular complication arose from the later identified Cu-Zn disorder, present in virtually all thin film samples. Subsequently, it has been shown experimentally that this unavoidable Cu-Zn disorder leads to a band gap lowering in the respective samples. Additional theoretical investigations, mostly based on Monte-Carlo methods, tried to understand the atomistic origin of this disorder induced band gap lowering. Here, we present theoretical results from first-principles calculations based on density functional theory for the disorder induced band gap lowering in kesterite Cu2ZnSnSe4and Ag2ZnSnSe4, where the Cu-Zn and Ag-Zn disorder is modelled via a supercell approach and special quasirandom structures. Results of subsequent analyses of structural, electronic, and optical properties are discussed with respect to available experimental results, and will provide additional insight and knowledge towards the atomistic origin of the observed disorder induced band gap lowering in kesterite type materials.

Noise reduction of stochastic density functional theory for metals

Density Functional Theory (DFT) has become a cornerstone in the modeling of metals. However, accurately simulating metals, particularly under extreme conditions, presents two significant challenges. First, simulating complex metallic systems at low electron temperatures is difficult due to their highly delocalized density matrix. Second, modeling metallic warm-dense materials at very high electron temperatures is challenging because it requires the computation of a large number of partially occupied orbitals. This study demonstrates that both challenges can be effectively addressed using the latest advances in linear-scaling stochastic DFT methodologies. Despite the inherent introduction of noise into all computed properties by stochastic DFT, this research evaluates the efficacy of various noise reduction techniques under different thermal conditions. Our observations indicate that the effectiveness of noise reduction strategies varies significantly with the electron temperature. Furthermore, we provide evidence that the computational cost of stochastic DFT methods scales linearly with system size for metal systems, regardless of the electron temperature regime.

Advances in Computational Methods for Modeling Photocatalytic Reactions: A Review of Recent Developments

Photocatalysis is a fascinating process in which a photocatalyst plays a pivotal role in driving a chemical reaction when exposed to light. Its capacity to harness light energy triggers a cascade of reactions that lead to the formation of intermediate compounds, culminating in the desired final product(s). The essence of this process is the interaction between the photocatalyst's excited state and its specific interactions with reactants, resulting in the creation of intermediates. The process's appeal is further enhanced by its cyclic nature-the photocatalyst is rejuvenated after each cycle, ensuring ongoing and sustainable catalytic action. Nevertheless, comprehending the photocatalytic process through the modeling of photoactive materials and molecular devices demands advanced computational techniques founded on effective quantum chemistry methods, multiscale modeling, and machine learning. This review analyzes contemporary theoretical methods, spanning a range of lengths and accuracy scales, and assesses the strengths and limitations of these methods. It also explores the future challenges in modeling complex nano-photocatalysts, underscoring the necessity of integrating various methods hierarchically to optimize resource distribution across different scales. Additionally, the discussion includes the role of excited state chemistry, a crucial element in understanding photocatalysis.

Band Gap Engineering of MoxW1-xS2 Alloy Monolayers with Wafer-Scale Uniformity

Modulating the band gap of two-dimensional (2D) transition metal dichalcogenide (TMDC) semiconductors is critical for their application in a wider spectral range. Alloying has been demonstrated as an effective method for regulating the band gap of 2D TMDC semiconductors. The fabrication of large-area 2D TMDC alloy films with centimeter-scale uniformity is fundamental to the application of integrated devices. Herein, we report a liquid-phase precursor one-step chemical vapor deposition (CVD) method for fabricating a MoxW1-xS2 alloy monolayer with a large size and an adjustable band gap. Good crystalline quality and high uniformity on a wafer scale enable the continuous adjustment of its band gap in the range of 1.8-2.0 eV. Density functional theory calculations provided a deep understanding of the Raman-active vibration modes of the MoxW1-xS2 alloy monolayer and the change in the conductivity of the alloy with photon energy. The synthesis of large-area MoxW1-xS2 alloy monolayers is a critical step toward the application of 2D layered semiconductors in practical optoelectronic devices.

Local Structures of Ex-Solved Nanoparticles Identified by Machine-Learned Potentials

In this study, we identify the local structures of ex-solved nanoparticles using machine-learned potentials (MLPs). We develop a method for training machine-learned potentials by sampling local structures of heterointerface configurations as a training set with its efficacy tested on the Ni/MgO system, illustrating that the error in interface energy is only 0.004 eV/Å2. Using the developed scheme, we train an MLP for the Ni/La0.5Ca0.5TiO3 ex-solution system and identify the local structures for both exo- and endo-type particles. The established model aligns well with the experimental observations, accurately predicting a nucleation size of 0.45 nm. Lastly, the density functional theory calculations on the established atomistic model verify that the kinetic barrier for the dry reforming of methane are substantially reduced by 0.49 eV on the ex-solved catalysts compared to that on the impregnated catalysts. Our findings offer insights into the local structures, growth mechanisms, and underlying origin of the catalytic properties of ex-solved nanoparticles.

Challenges for density functional theory in simulating metal-metal singlet bonding: A case study of dimerized VO2

VO2 is renowned for its electric transition from an insulating monoclinic (M1) phase, characterized by V-V dimerized structures, to a metallic rutile (R) phase above 340 K. This transition is accompanied by a magnetic change: the M1 phase exhibits a non-magnetic spin-singlet state, while the R phase exhibits a state with local magnetic moments. Simultaneous simulation of the structural, electric, and magnetic properties of this compound is of fundamental importance, but the M1 phase alone has posed a significant challenge to the density functional theory (DFT). In this study, we show none of the commonly used DFT functionals, including those combined with on-site Hubbard U to treat 3d electrons better, can accurately predict the V-V dimer length. The spin-restricted method tends to overestimate the strength of the V-V bonds, resulting in a small V-V bond length. Conversely, the spin-symmetry-breaking method exhibits the opposite trends. Each of these two bond-calculation methods underscores one of the two contentious mechanisms, i.e., Peierls lattice distortion or Mott localization due to electron-electron repulsion, involved in the metal-insulator transition in VO2. To elucidate the challenges encountered in DFT, we also employ an effective Hamiltonian that integrates one-dimensional magnetic sites, thereby revealing the inherent difficulties linked with the DFT computations.

The rise and fall of stretched bond errors: Extending the analysis of Perdew-Zunger self-interaction corrections of reaction barrier heights beyond the LSDA

Incorporating self-interaction corrections (SIC) significantly improves chemical reaction barrier height predictions made using density functional theory methods. We present a detailed orbital-by-orbital analysis of these corrections for three semi-local density functional approximations (DFAs) situated on the three lowest rungs of Jacob's ladder of approximations. The analysis is based on Fermi-Löwdin Orbital Self-Interaction Correction (FLOSIC) calculations performed at several steps along the reaction pathway from the reactants (R) to the transition state (TS) to the products (P) for four representative reactions selected from the BH76 benchmark set. For all three functionals, the major contribution to self-interaction corrections of the barrier heights can be traced to stretched bond orbitals that develop near the TS configuration. The magnitude of the ratio of the self-exchange-correlation energy to the self-Hartree energy (XC/H) for a given orbital is introduced as an indicator of one-electron self-interaction error. XC/H = 1.0 implies that an orbital's self-exchange-correlation energy exactly cancels its self-Hartree energy and that the orbital, therefore, makes no contribution to the SIC in the FLOSIC scheme. For the practical DFAs studied here, XC/H spans a range of values. The largest values are obtained for stretched or strongly lobed orbitals. We show that significant differences in XC/H for corresponding orbitals in the R, TS, and P configurations can be used to identify the major contributors to the SIC of barrier heights and reaction energies. Based on such comparisons, we suggest that barrier height predictions made using the strongly constrained and appropriately normed meta-generalized gradient approximation may have attained the best accuracy possible for a semi-local functional using the Perdew-Zunger SIC approach.

Structure and dynamics of liquid water from ab initio simulations: adding Minnesota density functionals to Jacob's ladder

The accurate representation of the structural and dynamical properties of water is essential for simulating the unique behavior of this ubiquitous solvent. Here we assess the current status of describing liquid water using ab initio molecular dynamics, with a special focus on the performance of all the later generation Minnesota functionals. Findings are contextualized within the current knowledge on DFT for describing bulk water under ambient conditions and compared to experimental data. We find that, contrary to the prevalent idea that local and semilocal functionals overstructure water and underestimate dynamical properties, M06-L, revM06-L, and M11-L understructure water, while MN12-L and MN15-L overdistance water molecules due to weak cohesive effects. This can be attributed to a weakening of the hydrogen bond network, which leads to dynamical fingerprints that are over fast. While most of the hybrid Minnesota functionals (M06, M08-HX, M08-SO, M11, MN12-SX, and MN15) also yield understructured water, their dynamical properties generally improve over their semilocal counterparts. It emerges that exact exchange is a crucial component for accurately describing hydrogen bonds, which ultimately leads to corrections in both the dynamical and structural properties. However, an excessive amount of exact exchange strengthens hydrogen bonds and causes overstructuring and slow dynamics (M06-HF). As a compromise, M06-2X is the best performing Minnesota functional for water, and its D3 corrected variant shows very good structural agreement. From previous studies considering nuclear quantum effects (NQEs), the hybrid revPBE0-D3, and the rung-5 RPA (RPA@PBE) have been identified as the only two approximations that closely agree with experiments. Our results suggest that the M06-2X(-D3) functionals have the potential to further improve the reproduction of experimental properties when incorporating NQEs through path integral approaches. This work provides further proof that accurate modeling of water interactions requires the inclusion of both exact exchange and balanced (non-local) correlation, highlighting the need for higher rungs on Jacob's ladder to achieve predictive simulations of complex biological systems in aqueous environments.

Reductive Dechlorination of Chlorinated Ethenes at the Sulfidated Zero-Valent Iron Surface: A Mechanistic DFT Study

Sulfidated nano- and microscale zero-valent iron (S-(n)ZVI) has shown enhanced selectivity and reactive lifetime in the degradation of chlorinated ethenes (CEs) compared to pristine (n)ZVI. However, varying effects of sulfidation on the dechlorination rates of structurally similar CEs have been reported, with the underlying mechanisms remaining poorly understood. In this study, we investigated the β-dichloroelimination reactions of tetrachloroethene (PCE), trichloroethene (TCE), cis-1,2-dichloroethene (cis-DCE), and trans-1,2-dichloroethene (trans-DCE) at the S and Fe sites of several S-(n)ZVI surface models by using density functional theory. Dechlorination reactions were both kinetically and thermodynamically more favorable at Fe sites compared to S sites, indicating that maintaining the accessibility of reactive Fe sites is crucial for achieving high S-(n)ZVI reactivity with contaminants. At Fe sites adjacent to S atoms, the reactivity for CE dechlorination followed the order trans-DCE ≈ TCE > cis-DCE > PCE. PCE degradation was hindered at these sites due to the steric effects of S atoms. At the S sites, the energy barriers correlated with the CEs' energy of the lowest unoccupied molecular orbital in the order PCE < TCE < DCE isomers. Our findings reveal that the experimentally observed selectivity of S-(n)ZVI materials for individual CEs can be explained by an interplay of the varying reactivities of Fe and S sites in CE dechlorination reactions.

The Case for a Defect Genome Initiative

The Materials Genome Initiative (MGI) has streamlined the materials discovery effort by leveraging generic traits of materials, with focus largely on perfect solids. Defects such as impurities and perturbations, however, drive many attractive functional properties of materials. The rich tapestry of charge, spin, and bonding states hosted by defects are not accessible to elements and perfect crystals, and defects can thus be viewed as another class of "elements" that lie beyond the periodic table. Accordingly, a Defect Genome Initiative (DGI) to accelerate functional defect discovery for energy, quantum information, and other applications is proposed. First, major advances made under the MGI are highlighted, followed by a delineation of pathways for accelerating the discovery and design of functional defects under the DGI. Near-term goals for the DGI are suggested. The construction of open defect platforms and design of data-driven functional defects, along with approaches for fabrication and characterization of defects, are discussed. The associated challenges and opportunities are considered and recent advances towards controlled introduction of functional defects at the atomic scale are reviewed. It is hoped this perspective will spur a community-wide interest in undertaking a DGI effort in recognition of the importance of defects in enabling unique functionalities in materials.

Highly Potent and Low-Volume Concentration Additives for Durable Aqueous Zinc Batteries: Machine Learning-Enabled Performance Rationalization

The essential virtues of aqueous zinc battery chemistry stem from the energy-dense zinc metal anode and mild aqueous electrolytes. Yet, their incompatibility - as exposed by zinc's corrosion and associated dendrite problem - poses a challenge to achieving improved cycle life under practically relevant parameters. While electrolyte additives are a scalable strategy, additives that can function at low volume concentrations remain elusive. Here, through screening alkanol and alkanediol chemistries, 1,2-butanediol and pentanediol are unveiled as highly potent additives, which operate at a practical 1 volume% concentration owing to their ability to furnish dynamic solid-electrolyte interphase through pronounced interfacial filming. This unique mechanistic action renders effective corrosion and dendrite mitigation, resulting in up to five to twenty-fold zinc cyclability enhancement with a high Coulombic efficiency (up to 99.9%) and improved full-cell performance under demanding conditions, including at elevated temperatures. A machine learning-based analysis is presented to rationalize the additive performance relative to critical physicochemical descriptors, which can pave the way for a rational approach to efficient additive discoveries.

Origin of Multiferroism in VOX2 (X = Cl, Br, I) Monolayers

Based on the proposed microscopic model, we investigate the multiferroic characteristics of VOX2 (X = Cl, Br, I) monolayers using a Green's function method. The dependence of the microscopic parameters of the ferroelectric system (pseudo-spin arrangement and flipping rate) on the magnitude and sign of the exchange magnetic interaction along the b-axis and the value of the Dzyaloshinskii-Moria vector have been investigated and qualitatively explained. The possibility of observing a spin-reorientation transition with a change in the character of spin ordering from antiferromagnetic to ferromagnetic is investigated. It is found that the antisymmetric magnetoelectric interaction may be responsible for the spin-reorientation transition without a change in the ordering of magnetic moments. Changing the sign of the exchange magnetic interaction along the b-axis leads to ferromagnetic ordering without observing a spin-reorientation transition. The dependence of isotropic and antisymmetric magnetic interactions on the microscopic parameters of the ferroelectric system is qualitatively explained. A mechanism for the occurrence of the spin-reorientation transition is presented based on the proposed microscopic model. The obtained results qualitatively coincide with Density Functional Theory calculations.

Nonadiabatic Derivative Couplings through Multiple Franck-Condon Modes Dictate the Energy Gap Law for Near and Short-Wave Infrared Dye Molecules

Near infrared (NIR, 700-1000 nm) and short-wave infrared (SWIR, 1000-2000 nm) dye molecules exhibit significant nonradiative decay rates from the first singlet excited state to the ground state. While these trends can be empirically explained by a simple energy gap law, detailed mechanisms of nearly universal behavior have remained unsettled for many cases. Theoretical and experimental results for two representative NIR/SWIR dye molecules reported here clarify the key mechanism for the observed energy gap law behavior. It is shown that the first derivative nonadiabatic coupling terms serve as major coupling pathways for nonadiabatic decay processes from the first excited singlet state to the ground state for these NIR and SWIR dye molecules and that vibrational modes other than the highest frequency modes also make significant contributions to the rate. This assessment is corroborated by further theoretical comparison with possible alternative mechanisms of intersystem crossing to triplet states and also by comparison with experimental data for deuterated molecules.

Comparing first-principles density functionals plus corrections for the lattice dynamics of YBa2Cu3O6

The enigmatic mechanism underlying unconventional high-temperature superconductivity, especially the role of lattice dynamics, has remained a subject of debate. Theoretical insights have long been hindered due to the lack of an accurate first-principles description of the lattice dynamics of cuprates. Recently, using the r2SCAN meta-generalized gradient approximation (meta-GGA) functional, we have been able to achieve accurate phonon spectra of an insulating cuprate YBa2Cu3O6 and discover significant magnetoelastic coupling in experimentally interesting Cu-O bond stretching optical modes [Ning et al., Phys. Rev. B 107, 045126 (2023)]. We extend this work by comparing Perdew-Burke-Ernzerhof and r2SCAN performances with corrections from the on-site Hubbard U and the D4 van der Waals (vdW) methods, aiming at further understanding on both the materials science side and the density functional side. We demonstrate the importance of vdW and self-interaction corrections for accurate first-principles YBa2Cu3O6 lattice dynamics. Since r2SCAN by itself partially accounts for these effects, the good performance of r2SCAN is now more fully explained. In addition, the performances of the Tao-Mo series of meta-GGAs, which are constructed in a different way from the strongly constrained and appropriately normed (SCAN) meta-GGA and its revised version r2SCAN, are also compared and discussed.

A first-principles machine-learning force field for heterogeneous ice nucleation on microcline feldspar

The formation of ice in the atmosphere affects precipitation and cloud properties, and plays a key role in the climate of our planet. Although ice can form directly from liquid water under deeply supercooled conditions, the presence of foreign particles can aid ice formation at much warmer temperatures. Over the past decade, experiments have highlighted the remarkable efficiency of feldspar minerals as ice nuclei compared to other particles present in the atmosphere. However, the exact mechanism of ice formation on feldspar surfaces has yet to be fully understood. Here, we develop a first-principles machine-learning model for the potential energy surface aimed at studying ice nucleation at microcline feldspar surfaces. The model is able to reproduce with high-fidelity the energies and forces derived from density-functional theory (DFT) based on the SCAN exchange and correlation functional. Our training set includes configurations of bulk supercooled water, hexagonal and cubic ice, microcline, and fully-hydroxylated feldspar surfaces exposed to a vacuum, liquid water, and ice. We apply the machine-learning force field to study different fully-hydroxylated terminations of the (100), (010), and (001) surfaces of microcline exposed to a vacuum. Our calculations suggest that terminations that do not minimize the number of broken bonds are preferred in a vacuum. We also study the structure of supercooled liquid water in contact with microcline surfaces, and find that water density correlations extend up to around 10 Å from the surfaces. Finally, we show that the force field maintains a high accuracy during the simulation of ice formation at microcline surfaces, even for large systems of around 30 000 atoms. Future work will be directed towards the calculation of nucleation free-energy barriers and rates using the force field developed herein, and understanding the role of different microcline surfaces in ice nucleation.

Defect passivation in methylammonium/bromine free inverted perovskite solar cells using charge-modulated molecular bonding

Molecular passivation is a prominent approach for improving the performance and operation stability of halide perovskite solar cells (HPSCs). Herein, we reveal discernible effects of diammonium molecules with either an aryl or alkyl core onto Methylammonium-free perovskites. Piperazine dihydriodide (PZDI), characterized by an alkyl core-electron cloud-rich-NH terminal, proves effective in mitigating surface and bulk defects and modifying surface chemistry or interfacial energy band, ultimately leading to improved carrier extraction. Benefiting from superior PZDI passivation, the device achieves an impressive efficiency of 23.17% (area ~1 cm2) (low open circuit voltage deficit ~0.327 V) along with superior operational stability. We achieve a certified efficiency of ~21.47% (area ~1.024 cm2) for inverted HPSC. PZDI strengthens adhesion to the perovskite via -NH2I and Mulliken charge distribution. Device analysis corroborates that stronger bonding interaction attenuates the defect densities and suppresses ion migration. This work underscores the crucial role of bifunctional molecules with stronger surface adsorption in defect mitigation, setting the stage for the design of charge-regulated molecular passivation to enhance the performance and stability of HPSC.

AFLOW-CCE for the thermodynamics of ionic materials

Accurate thermodynamic stability predictions enable data-driven computational materials design. Standard density functional theory (DFT) approximations have limited accuracy with average errors of a few hundred meV/atom for ionic materials, such as oxides and nitrides. Thus, insightful correction schemes as given by the coordination corrected enthalpies (CCE) method, based on an intuitive parametrization of DFT errors with respect to coordination numbers and cation oxidation states, present a simple, yet accurate solution to enable materials stability assessments. Here, we illustrate the computational capabilities of our AFLOW-CCE software by utilizing our previous results for oxides and introducing new results for nitrides. The implementation reduces the deviations between theory and experiment to the order of the room temperature thermal energy scale, i.e., ∼25 meV/atom. The automated corrections for both materials classes are freely available within the AFLOW ecosystem via the AFLOW-CCE module, requiring only structural inputs.

Prediction of transition metal carbonitride monolayers MN4C6 (M = Cr, Mn, Fe, and Co) made up of a benzene ring and a planar MN4 moiety

Based on first-principles calculations, we predict a class of graphene-like magnetic materials, transition metal carbonitrides MN4C6 (M = Cr, Mn, Fe, and Co), which are made up of a benzene ring and an MN4 moiety, two common planar units in the compounds. The structural stability is demonstrated by the phonon and molecular dynamics calculations, and the formation mechanism of the planar geometry of MN4C6 is ascribed to the synergistic effect of sp2 hybridization, M-N coordination bond, and π-d conjugation. The MN4C6 materials consist of only one layer of atoms and the transition metal atom is located in the planar crystal field, which is markedly different from most two-dimensional materials. The calculations indicate that MnN4C6, FeN4C6, and CoN4C6 are ferromagnetic while CrN4C6 has an antiferromagnetic ground state. The Curie temperatures are estimated by solving the anisotropic Heisenberg model with the Monte Carlo method.

Self-interaction corrected SCAN functional for molecules and solids in the numeric atom-center orbital framework

Semilocal density-functional approximations (DFAs), including the state-of-the-art SCAN functional, are plagued by the self-interaction error (SIE). While this error is explicitly defined only for one-electron systems, it has inspired the self-interaction correction method proposed by Perdew and Zunger (PZ-SIC), which has shown promise in mitigating the many-electron SIE. However, the PZ-SIC method is known for its significant numerical instability. In this study, we introduce a novel constraint that facilitates self-consistent localization of the SIC orbitals in the spirit of Edmiston-Ruedenberg orbitals [Rev. Mod. Phys. 35, 457 (1963)]. Our practical implementation within the all-electron numeric atom-centered orbitals code FHI-aims guarantees efficient and stable convergence of the self-consistent PZ-SIC equations for both molecules and solids. We further demonstrate that our PZ-SIC approach effectively mitigates the SIE in the meta-generalized gradient approximation SCAN functional, significantly improving the accuracy for ionization potentials, charge-transfer energies, and bandgaps for a diverse selection of molecules and solids. However, our PZ-SIC method does have its limitations. It cannot improve the already accurate SCAN results for properties such as cohesive energies, lattice constants, and bulk modulus in our test sets. This highlights the need for new-generation DFAs with more comprehensive applicability.

Stable room-temperature ferromagnetism and gate-tunable quantum anomalous Hall effect of two-dimensional 5d transition-metal trihalide OsX3 (X = Cl, Br, I) monolayers

5d transition-metal compounds are usually not expected to exhibit distinct magnetic ordering owing to their substantial binding energy associated with 5d electrons. In this study, we demonstrate that two-dimensional (2D) 5d transition-metal Os trihalide OsX3 monolayers can exhibit room-temperature ferromagnetism and quantum anomalous Hall effect (QAHE) by utilizing density functional theory and Monte Carlo simulation. Our calculation results of coexisting Raman and infrared activities of lattice vibration reveal the structural stability of 2D OsX3 (X = Cl, Br, I) and structural instability of 2D OsX3 (X = F). Furthermore, all 2D OsX3 trihalides (X = Cl, Br, I) are half-metals, and their ferromagnetism remains stable under ambient temperature, where 2D OsCl3 and OsBr3 have an in-plane easy axis while 2D OsI3 has an out-of-plane easy axis. Notably, when spin-orbit coupling is included, the gate-tunable QAHE could emerge in ferromagnetic 2D OsI3, while 2D OsCl3 and OsBr3 are topologically trivial. Additionally, the magnon bands of 2D OsX3 (X = Cl, Br, I) possess two spin-wave branches with dispersion similar to that of the Dirac cone in the electronic structure of graphene, which are attributed to the unique ferromagnetic honeycomb sublattice of osmium atoms.

Assessing r2SCAN meta-GGA functional for structural parameters, cohesive energy, mechanical modulus, and thermophysical properties of 3d, 4d, and 5d transition metals

The recent development of accurate and efficient semilocal density functionals on the third rung of Jacob's ladder of density functional theory, such as the revised regularized strongly constrained and appropriately normed (r2SCAN) density functional, could enable rapid and highly reliable prediction of the elasticity and temperature dependence of thermophysical parameters of refractory elements and their intermetallic compounds using the quasi-harmonic approximation (QHA). Here, we present a comparative evaluation of equilibrium cell volumes, cohesive energy, mechanical moduli, and thermophysical properties (Debye temperature and thermal expansion coefficient) for 22 transition metals using semilocal density functionals, including the local density approximation (LDA), Perdew-Burke-Ernzerhof (PBE) and PBEsol generalized gradient approximations (GGAs), and the r2SCAN meta-GGA. PBEsol and r2SCAN deliver the same level of accuracies for structural, mechanical, and thermophysical properties. PBE and r2SCAN perform better than LDA and PBEsol for calculating cohesive energies of transition metals. Among the tested density functionals, r2SCAN provides an overall well-balanced performance for reliably computing cell volumes, cohesive energies, mechanical properties, and thermophysical properties of various 3d, 4d, and 5d transition metals using QHA. Therefore, we recommend that r2SCAN could be employed as a workhorse method to evaluate thermophysical properties of transition metal compounds and alloys in high throughput workflows.

Balance between Physical Interpretability and Energetic Predictability in Widely Used Dispersion-Corrected Density Functionals

We assess the performance of different dispersion models for several popular density functionals across a diverse set of noncovalent systems, ranging from the benzene dimer to molecular crystals. By analyzing the interaction energies and their individual components, we demonstrate that there exists variability across different systems for empirical dispersion models, which are calibrated for reproducing the interaction energies of specific systems. Thus, parameter fitting may undermine the underlying physics, as dispersion models rely on error compensation among the different components of the interaction energy. Energy decomposition analyses reveal that, the accuracy of revPBE-D3 for some aqueous systems originates from significant compensation between dispersion and charge transfer energies. However, revPBE-D3 is less accurate in describing systems where error compensation is incomplete, such as the benzene dimer. Such cases highlight the propensity for unpredictable behavior in various dispersion-corrected density functionals across a wide range of molecular systems, akin to the behavior of force fields. On the other hand, we find that SCAN-rVV10, a targeted-dispersion approach, affords significant reductions in errors associated with the lattice energies of molecular crystals, while it has limited accuracy in reproducing structural properties. Given the ubiquitous nature of noncovalent interactions and the key role of density functional theory in computational sciences, the future development of dispersion models should prioritize the faithful description of the dispersion energy, a shift that promises greater accuracy in capturing the underlying physics across diverse molecular and extended systems.

My life in science: Lessons for yours?

Because of an acquired obsession to understand as much as possible in a limited but important area of science and because of optimism, luck, and help from others, my life in science turned out to be much better than I or others could have expected or planned. This is the story of how that happened, and also the story of the groundstate density functional theory of electronic structure, told from a personal perspective.

Improved ion adsorption capacities and diffusion dynamics in surface anchored MoS2⊥Mo4/3B2 and MoS2⊥Mo4/3B2O2 heterostructures as anodes for alkaline metal-ion batteries

First-principles calculations were performed to analyze the atomic structures and electrochemical energy storage properties of novel MoS2⊥boridene heterostructures by anchoring MoS2 nanoflakes on Mo4/3B2 and Mo4/3B2O2 monolayers. Both thermodynamic and thermal stabilities of each heterostructure were thoroughly evaluated from the obtained binding energies and through first-principles molecular dynamics simulations at room temperature, confirming the high formability of the heterostructures. The electrochemical properties of MoS2⊥Mo4/3B2 and MoS2⊥Mo4/3B2O2 heterostructures were investigated for their potential use as anodes for alkaline metal ion batteries (Li+, Na+ and K+). It was revealed that Li+ and Na+ can form multiple stable full adsorption layers on both heterostructures, while K+ forms only a single full adsorption layer. The presence of a negative electron cloud (NEC) contributes to the stabilization of a multi-layer adsorption mechanism. For all investigated alkaline metal ions, the predicted ion diffusion dynamics are relatively sluggish for the adsorbates in the first full adsorption layer on MoS2⊥boridene heterostructures due the relatively large migration energies (>0.50 eV), compared to those of second or third full adsorption layers (<0.30 eV). MoS2⊥Mo4/3B2O2 exhibited higher onset and mean open circuit voltages as anodes for alkaline metal-ion batteries than MoS2⊥Mo4/3B2 hybrids because of enhanced interactions between the adsorbate and the Mo4/3B2O2 monolayer with the presence of O-terminations. Tailoring the size and horizontal spacing between two neighboring MoS2 nano-flakes in heterostructures led to high theoretical capacities for LIBs (531 mA h g-1), SIBs (300 mA h g-1) and PIBs (131 mA h g-1) in the current study.

Structural, electronic and magnetic properties of greigite Fe3S4by GGA and GGA+Uversus SCAN meta-GGA density functionals

The performance of exchange-correlation functional of density functional theory represented in generalized gradient approximation (GGA) and in the strongly constrained and appropriately normed (SCAN) meta-GGA scheme to study structural, electronic, and magnetic properties of greigite (Fe3S4) was investigated. The effects of inclusion of strong electron correlations represented by on-site Hubbard correctionU, and nonlocality of the long-range van der Waals (vdW) interactions were also considered. Geometry optimization yielded the inverse spinel structure and lattice parameter of greigite in good agreement with experimental data. Calculated electronic structure revealed a half-metallic nature of the greigite bands for the applied functionals except for GGA, which predicts metallic behavior. Antiferromagnetic coupling of iron ions in tetrahedral and octahedral coordinations makes the overall crystal structure ferrimagnetic. In general the GGA+Uand SCAN show comparable performance in prediction physical properties of greigite. Inclusion of the vdW correction does not change the character of the bands.

Ion-Type Dependence of DNA Electronic Excitation in Water under Proton, α-Particle, and Carbon Ion Irradiation: A First-Principles Simulation Study

Understanding how the electronic excitation of DNA changes in response to different high-energy particles is central to advancing ion beam cancer therapy and other related approaches, such as boron neutron capture therapy. While protons have been the predominant ions of choice in ion beam cancer therapy, heavier ions, particularly carbon ions, have drawn significant attention over the past decade. Carbon ions are expected to transfer larger amounts of energy according to linear response theory. However, molecular-level details of the electronic excitation under heavier ion irradiation remain unknown. In this work, we use real-time time-dependent density functional theory simulations to examine the quantum-mechanical details of DNA electronic excitations in water under proton, α-particle, and carbon ion irradiation. Our results show that the energy transfer does indeed increase for the heavier ions, while the excitation remains highly conformal. However, the increase in the energy transfer rate, measured by electronic stopping power, does not match the prediction by the linear response model, even when accounting for the velocity dependence of the irradiating ion's charge. The simulations also reveal that while the number of holes generated on DNA increases for heavier ions, the increase is only partially responsible for the larger stopping power. Larger numbers of highly energetic holes formed from the heavier ions also contribute significantly to the increased electronic stopping power.

Simulating Crystal Structure, Acidity, Proton Distribution, and IR Spectra of Acid Zeolite HSAPO-34: A High Accuracy Study

It is a challenge to characterize the acid properties of microporous materials in either experiments or theory. This study presents the crystal structure, acid site, acid strength, proton siting, and IR spectra of HSAPO-34 from the SCAN + rVV10 method. The results indicate: the crystal structures of various acid sites of HSAPO-34 deviate from the space group of R3¯; the acid strength inferred from the DPE value likely decreases with the proton binding sites at O(2), O(4), O(1),and O(3), contrary to the stability order in view of the internal energy; the calculated ensemble-averaged DPE is about 1525 kJ/mol at 673.15 K; and the proton siting and the proton distribution are distinctly influenced by the temperature: at low temperatures, the proton is predominantly located at O(3), while it prefers O(2) at high temperatures, and the proton at O(4) assumedly has the least distribution at 273.15-773.15 K. In line with the neutron diffraction experiment, a correction factor of 0.979 is needed to correct for the calculated hydroxyl stretching vibration (ν(O-H)) of HSAPO-34. It seems that the SCAN meta-GGA method, compensating for some drawbacks of the GGA method, could provide satisfying results regarding the acid properties of HSAPO-34.

Reworking the Tao-Mo exchange-correlation functional. I. Reconsideration and simplification

The revised, regularized Tao-Mo (rregTM) exchange-correlation density functional approximation (DFA) [A. Patra, S. Jana, and P. Samal, J. Chem. Phys. 153, 184112 (2020) and Jana et al., J. Chem. Phys. 155, 024103 (2021)] resolves the order-of-limits problem in the original TM formulation while preserving its valuable essential behaviors. Those include performance on standard thermochemistry and solid data sets that is competitive with that of the most widely explored meta-generalized-gradient-approximation DFAs (SCAN and r2SCAN) while also providing superior performance on elemental solid magnetization. Puzzlingly however, rregTM proved to be intractable for de-orbitalization via the approach of Mejía-Rodríguez and Trickey [Phys. Rev. A 96, 052512 (2017)]. We report investigation that leads to diagnosis of how the regularization in rregTM of the z indicator functions (z = the ratio of the von-Weizsäcker and Kohn-Sham kinetic energy densities) leads to non-physical behavior. We propose a simpler regularization that eliminates those oddities and that can be calibrated to reproduce the good error patterns of rregTM. We denote this version as simplified, regularized Tao-Mo, sregTM. We also show that it is unnecessary to use rregTM correlation with sregTM exchange: Perdew-Burke-Ernzerhof correlation is sufficient. The subsequent paper shows how sregTM enables some progress on de-orbitalization.

Reworking the Tao-Mo exchange-correlation functional. II. De-orbitalization

In Paper I [H. Francisco, A. C. Cancio, and S. B. Trickey, J. Chem. Phys. 159, 214102 (2023)], we gave a regularization of the Tao-Mo exchange functional that removes the order-of-limits problem in the original Tao-Mo form and also eliminates the unphysical behavior introduced by an earlier regularization while essentially preserving compliance with the second-order gradient expansion. The resulting simplified, regularized (sregTM) functional delivers performance on standard molecular and solid state test sets equal to that of the earlier revised, regularized Tao-Mo functional. Here, we address de-orbitalization of that new sregTM into a pure density functional. We summarize the failures of the Mejía-Rodríguez and Trickey de-orbitalization strategy [Phys. Rev. A 96, 052512 (2017)] when used with both versions. We discuss how those failures apparently arise in the so-called z' indicator function and in substitutes for the reduced density Laplacian in the parent functionals. Then, we show that the sregTM functional can be de-orbitalized somewhat well with a rather peculiarly parameterized version of the previously used deorbitalizer. We discuss, briefly, a de-orbitalization that works in the sense of reproducing error patterns but that apparently succeeds by cancelation of major qualitative errors associated with the de-orbitalized indicator functions α and z, hence, is not recommended. We suggest that the same issue underlies the earlier finding of comparatively mediocre performance of the de-orbitalized Tao-Perdew-Staroverov-Scuseri functional. Our work demonstrates that the intricacy of such two-indicator functionals magnifies the errors introduced by the Mejía-Rodríguez and Trickey de-orbitalization approach in ways that are extremely difficult to analyze and correct.

An autonomous laboratory for the accelerated synthesis of novel materials

To close the gap between the rates of computational screening and experimental realization of novel materials1,2, we introduce the A-Lab, an autonomous laboratory for the solid-state synthesis of inorganic powders. This platform uses computations, historical data from the literature, machine learning (ML) and active learning to plan and interpret the outcomes of experiments performed using robotics. Over 17 days of continuous operation, the A-Lab realized 41 novel compounds from a set of 58 targets including a variety of oxides and phosphates that were identified using large-scale ab initio phase-stability data from the Materials Project and Google DeepMind. Synthesis recipes were proposed by natural-language models trained on the literature and optimized using an active-learning approach grounded in thermodynamics. Analysis of the failed syntheses provides direct and actionable suggestions to improve current techniques for materials screening and synthesis design. The high success rate demonstrates the effectiveness of artificial-intelligence-driven platforms for autonomous materials discovery and motivates further integration of computations, historical knowledge and robotics.

Comprehensive Study of Lithium Diffusion in Si/C-Layer and Si/C3N4 Composites in a Faceted Crystalline Silicon Anode for Fast-Charging Lithium-Ion Batteries

By using silicon (Si) as an anode of lithium-ion batteries, the capacity can be significantly increased, but relatively large volume expansion limits the application as an efficient anode material. Huge volume expansion of the silicon anode during lithiation, however, leads to cracking and losing its connection with the current collector. This shortcoming can be improved by the deposition of a nanometric carbon- or nitrogen-doped carbon coating on the silicon surface, resulting in Si/C-layer and Si/C3N4 interfaces. In this work, Li+ diffusion in Si/C-layer and Si/C3N4 composite materials along three Si surfaces and various ion pathways were carefully analyzed by using density functional theory and ab initio molecular dynamic (AIMD) simulations. Both Si/C and Si/C3N4 interfaces and three Si surfaces of (100), (110), and (111) were investigated. The formation of nitrogen holes and monatomic carbon binders in the composite increases ion diffusivity and limits volume expansion. Furthermore, the Bader analysis shows that the type and orientation of the surfaces have important effects on ion distribution. The results indicated that the C3N4 composite increases Li+ diffusion in Si (100) from 7.82 × 10-5 to 3.17 × 10-4 cm2/s. The presented results provide a guide for the appropriate design of stable and safe high-energy-density batteries.

Enhanced immobilization of Arsenic(III) and Auto-oxidation to Arsenic(V) by titanium oxide (TiO2), due to Single-Atom vacancies and oxyanion formation

Pollution control of As(III), a naturally occurring carcinogen, has recently gained a global attention, while due to the dominance of neutral H3AsO3 over a wide pH range, As(III) immobilization by most minerals is not efficient as As(V) immobilization. TiO2 shows promise for controlling As(III) pollution, and herein, a comprehensive study about As(III) adsorption by TiO2 and oxyanion formation is conducted by means of DFT + D3 methods. Both anatase and rutile are effective for As(III) adsorption, while As(III) adsorption affinities differ significantly and are -1.48 and -3.79 eV for pristine surfaces, ascend to -3.85 and -5.08 eV for O vacancies, and further to -5.37 and -5.26 eV for Ti vacancies, respectively. The bidentate binuclear complexes dominate for pristine surfaces, and O vacancies prefer OAs insertion into TiO2 lattice, while for Ti vacancies, all As(III) centers are auto-oxidized to As(V). Ti-3d, O-2p or/and As-4p rather than other orbitals contribute significantly to As adsorption, and O and Ti vacancies promote adsorption through stronger orbital hybridization. The superior adsorption for Ti vacancies originates from As(V) formation instead of bonding interactions. The formation of As oxyanions, which may occur spontaneously at pristine surfaces and is greatly promoted by O and Ti vacancies, enhances As(III) adsorption pronouncedly and becomes a viable strategy for As(III) immobilization. H2AsO3- and HAsO32- dominate for pristine surfaces and O vacancies, and for Ti vacancies, H2AsO4- and HAsO42- dominate over anatase whereas AsO43- also makes an important contribution over rutile. Results rationalize experimental observations available, and provide significantly new insights about the migration, bioavailability and fate of As(III) over TiO2 surfaces that facilitate the exploration of scavengers for As and other pollutants.

Crystal Structures and Phase Stability of the Li2S-P2S5 System from First Principles

The Li2S-P2S5 pseudo-binary system has been a valuable source of promising superionic conductors, with α-Li3PS4, β-Li3PS4, HT-Li7PS6, and Li7P3S11 having excellent room-temperature Li-ion conductivity >0.1 mS/cm. The metastability of these phases at ambient temperature motivates a study to quantify their thermodynamic accessibility. Through calculating the electronic, configurational, and vibrational sources of free energy from first principles, a phase diagram of the crystalline Li2S-P2S5 space is constructed. New ground-state orderings are proposed for α-Li3PS4, HT-Li7PS6, LT-Li7PS6, and Li7P3S11. Well-established phase stability trends from experiments are recovered, such as polymorphic phase transitions in Li7PS6 and Li3PS4, and the instability of Li7P3S11 at high temperature. At ambient temperature, it is predicted that all superionic conductors in this space are indeed metastable but thermodynamically accessible. Vibrational and configurational sources of entropy are shown to be essential toward describing the stability of superionic conductors. New details of the Li sublattices are revealed and are found to be crucial toward accurately predicting configurational entropy. All superionic conductors contain significant configurational entropy, which suggests an inherent correlation between fast Li diffusion and thermodynamic stability arising from the configurational disorder.

Consistent density functional theory-based description of ion hydration through density-corrected many-body representations

Delocalization error constrains the accuracy of density functional theory in describing molecular interactions in ion-water systems. Using Na+ and Cl- in water as model systems, we calculate the effects of delocalization error in the SCAN functional for describing ion-water and water-water interactions in hydrated ions, and demonstrate that density-corrected SCAN (DC-SCAN) predicts n-body and interaction energies with an accuracy approaching coupled cluster theory. The performance of DC-SCAN is size-consistent, maintaining an accurate description of molecular interactions well beyond the first solvation shell. Molecular dynamics simulations at ambient conditions with many-body MB-SCAN(DC) potentials, derived from the many-body expansion, predict the solvation structure of Na+ and Cl- in quantitative agreement with reference data, while simultaneously reproducing the structure of liquid water. Beyond rationalizing the accuracy of density-corrected models of ion hydration, our findings suggest that our unified density-corrected MB formalism holds great promise for efficient DFT-based simulations of condensed-phase systems with chemical accuracy.

Inelastic Neutron Scattering Study of Phonon Density of States of Iodine Oxides and First-Principles Calculations

Iodine oxides I2Oy (y = 4, 5, 6) crystallize into atypical structures that fall between molecular- and framework-base types and exhibit high reactivity in an ambient environment, a property highly desired in the so-called "agent defeat materials". Inelastic neutron scattering experiments were performed to determine the phonon density of states of the newly synthesized I2O5 and I2O6 samples. First-principles calculations were carried out for I2O4, I2O5, and I2O6 to predict their thermodynamic properties and phonon density of states. Comparison of the INS data with the Raman and infrared measurements as well as the first-principles calculations sheds light on their distinctive, anisotropic thermomechanical properties.

Short-range screened density matrix functional for proper descriptions of thermochemistry, thermochemical kinetics, nonbonded interactions, and singlet diradicals

Common one-electron reduced density matrix (1-RDM) functionals that depend on Coulomb and exchange-only integrals tend to underestimate dynamic correlation, preventing reduced density matrix functional theory (RDMFT) from achieving comparable accuracy to density functional theory in main-group thermochemistry and thermochemical kinetics. The recently developed ωP22 functional introduces a semi-local density functional to screen the erroneous short-range portion of 1-RDM functionals without double-counting correlation, potentially providing a better treatment of dynamic correlation around equilibrium geometries. Herein, we systematically evaluate the performance of this functional model, which consists of two parameters, on main-group thermochemistry, thermochemical kinetics, nonbonded interactions, and more. Tests on atomization energies, vibrational frequencies, and reaction barriers reveal that the ωP22 functional model can reliably predict properties at equilibrium and slightly away from equilibrium geometries. In particular, it outperforms commonly used density functionals in the prediction of reaction barriers, nonbonded interactions, and singlet diradicals, thus enhancing the predictive power of RDMFT for routine calculations of thermochemistry and thermochemical kinetics around equilibrium geometries. Further development is needed in the future to refine short- and long-range approximations in the functional model in order to achieve an excellent description of properties both near and far from equilibrium geometries.

Probing a Defect-Site-Specific Electronic Orbital in Graphene with Single-Atom Sensitivity

Visualization of individual electronic states ascribed to specific unoccupied orbitals at the atomic scale can reveal fundamental information about chemical bonding, but it is challenging since bonding often results in only subtle variations in the whole density of states. Here, we utilize atomic-resolution energy-loss near-edge fine structure (ELNES) spectroscopy to map out the electronic states attributed to specific unoccupied p_{z} orbital around a fourfold coordinated silicon point defect in graphene, which is further supported by theoretical calculations. Our results illustrate the power of atomic-resolution ELNES towards the probing of defect-site-specific electronic orbitals in monolayer crystals, providing insights into understanding the effect of chemical bonding on the local properties of defects in solids.

Stable and metastable structures of tin (IV) oxide at high pressure

We have analysed [Formula: see text] with a combination of synchrotron X-ray diffraction and X-ray absorption spectroscopy across a pressure range of [Formula: see text] GPa with thermal annealing by a [Formula: see text] laser allowing access to all of the known high-density polymorphs of [Formula: see text], and here report their crystallographic information. The metastability of the post-rutile [Formula: see text]-[Formula: see text] and [Formula: see text] structures in [Formula: see text] are investigated by experiment and PW-DFT simulations, revealing a complex energetic landscape and suggesting a significant dependence of the observed phases on the pressure-temperature pathway taken in experiment. This article is part of the theme issue 'Exploring the length scales, timescales and chemistry of challenging materials (Part 1)'.

One-dimensional carbon chains encapsulated in hollandite

One-dimensional carbon chains are highly reactive allotropes that are stabilized inside the protective environment of carbon nanotubes. Here we show that carbon chains can be encapsulated in metal oxides containing open structural channels, exemplified by hollandite α-MnO2. The α-MnO2 channels stabilize cumulene chains due to their structural commensurability, whereas the triple bonds in polyyne chains exhibit excessive steric repulsion to the oxide ions bordering the channel. Cumulene exhibits an interaction energy of only 0.065 eV per carbon atom, obtained by first-principles calculations, which is significantly more favorable than for encapsulation in a similarly sized carbon nanotube. Encapsulation of carbon chains is associated with lateral expansion of the α-MnO2 channel and polarization of the manganese and oxygen charge densities adjacent to the chains. Accordingly, the interaction energy is governed by a balance between van der Waals attraction and steric repulsion between the materials.