Applicability of the Strongly Constrained and Appropriately Normed Density Functional to Transition-Metal Magnetism

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.

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.

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.

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.

Testing the r2SCAN Density Functional for the Thermodynamic Stability of Solids with and without a van der Waals Correction

A central aim of materials discovery is an accurate and numerically reliable description of thermodynamic properties, such as the enthalpies of formation and decomposition. The r2SCAN revision of the strongly constrained and appropriately normed (SCAN) meta-generalized gradient approximation (meta-GGA) balances numerical stability with high general accuracy. To assess the r2SCAN description of solid-state thermodynamics, we evaluate the formation and decomposition enthalpies, equilibrium volumes, and fundamental band gaps of more than 1000 solids using r2SCAN, SCAN, and PBE, as well as two dispersion-corrected variants, SCAN+rVV10 and r2SCAN+rVV10. We show that r2SCAN achieves accuracy comparable to SCAN and often improves upon SCAN's already excellent accuracy. Although SCAN+rVV10 is often observed to worsen the formation enthalpies of SCAN and makes no substantial correction to SCAN's cell volume predictions, r2SCAN+rVV10 predicts marginally less accurate formation enthalpies than r2SCAN, and slightly more accurate cell volumes than r2SCAN. The average absolute errors in predicted formation enthalpies are found to decrease by a factor of 1.5 to 2.5 from the GGA level to the meta-GGA level. Smaller decreases in error are observed for decomposition enthalpies. For formation enthalpies r2SCAN improves over SCAN for intermetallic systems. For a few classes of systems-transition metals, intermetallics, weakly bound solids, and enthalpies of decomposition into compounds-GGAs are comparable to meta-GGAs. In total, r2SCAN and r2SCAN+rVV10 can be recommended as stable, general-purpose meta-GGAs for materials discovery.

Re-examining the giant magnetization density in α''-Fe16N2 with the SCAN+U method

We present an in-depth discussion of the magnetic ground state of α''-Fe₁₆N₂ within the framework of the density functional theory (DFT). The exchange-correlation effects are treated using a variety of schemes, including the local-spin-density approximation, the generalized-gradient approximation, and the Strongly-Constrained-and-Appropriately-Normed (SCAN) scheme. We also delineate effects of adding an on-site interaction parameter U on the Fe sites. Among all the schemes considered, only SCAN+U is found to capture the surprisingly large magnetization density in α''-Fe₁₆N₂ that has been observed experimentally. Our study shows how the combination of SCAN and self-interaction corrections applied on different Fe sites through the parameter U can reproduce both the correct equilibrium volume and the giant magnetization density of α''-Fe₁₆N₂.

Performance Enhancement of APW+lo Calculations by Simplest Separation of Concerns

Full-potential linearized augmented plane wave (LAPW) and APW plus local orbital (APW+lo) codes differ widely in both their user interfaces and in capabilities for calculations and analysis beyond their common central task of all-electron solution of the Kohn–Sham equations. However, that common central task opens a possible route to performance enhancement, namely to offload the basic LAPW/APW+lo algorithms to a library optimized purely for that purpose. To explore that opportunity, we have interfaced the Exciting-Plus (“EP”) LAPW/APW+lo DFT code with the highly optimized SIRIUS multi-functional DFT package. This simplest realization of the separation of concerns approach yields substantial performance over the base EP code via additional task parallelism without significant change in the EP source code or user interface. We provide benchmarks of the interfaced code against the original EP using small bulk systems, and demonstrate performance on a spin-crossover molecule and magnetic molecule that are of size and complexity at the margins of the capability of the EP code itself.

Scalable Molecular GW Calculations: Valence and Core Spectra

We present a scalable implementation of the GW approximation using Gaussian atomic orbitals to study the valence and core ionization spectroscopies of molecules. The implementation of the standard spectral decomposition approach to the screened-Coulomb interaction, as well as a contour-deformation method, is described. We have implemented both of these approaches using the robust variational fitting approximation to the four-center electron repulsion integrals. We have utilized the MINRES solver with the contour-deformation approach to reduce the computational scaling by 1 order of magnitude. A complex heuristic in the quasiparticle equation solver further allows a speed-up of the computation of core and semicore ionization energies. Benchmark tests using the GW100 and CORE65 data sets and the carbon 1s binding energy of the well-studied ethyl trifluoroacetate, or ESCA molecule, were performed to validate the accuracy of our implementation. We also demonstrate and discuss the parallel performance and computational scaling of our implementation using a range of water clusters of increasing size.

Accurate density functional made more versatile

We propose a one-electron self-interaction-free correlation energy functional compatible with the order-of-limit problem-free Tao-Mo (TM) semilocal functional (regTM) [J. Tao and Y. Mo, Phys. Rev. Lett. 117, 073001 (2016) and Patra et al., J. Chem. Phys. 153, 184112 (2020)] to be used for general purpose condensed matter physics and quantum chemistry. The assessment of the proposed functional for large classes of condensed matter and chemical systems shows its improvement in most cases compared to the TM functional, e.g., when applied to the relative energy difference of MnO₂ polymorphs. In this respect, the present exchange-correction functional, which incorporates the TM technique of the exchange hole model combined with the slowly varying density correction, can achieve broad applicability, being able to solve difficult solid-state problems.

Electronic-structure methods for materials design

The accuracy and efficiency of electronic-structure methods to understand, predict and design the properties of materials has driven a new paradigm in research. Simulations can greatly accelerate the identification, characterization and optimization of materials, with this acceleration driven by continuous progress in theory, algorithms and hardware, and by adaptation of concepts and tools from computer science. Nevertheless, the capability to identify and characterize materials relies on the predictive accuracy of the underlying physical descriptions, and on the ability to capture the complexity of realistic systems. We provide here an overview of electronic-structure methods, of their application to the prediction of materials properties, and of the different strategies employed towards the broader goals of materials design and discovery.

Describing adsorption of benzene, thiophene, and xenon on coinage metals by using the Zaremba-Kohn theory-based model

Semilocal (SL) density functional approximations (DFAs) are widely applied but have limitations due to their inability to incorporate long-range van der Waals (vdW) interaction. Non-local functionals (vdW-DF, VV10, and rVV10) or empirical methods (DFT+D, DFT+vdW, and DFT+MBD) are used with SL-DFAs to account for such missing interaction. The physisorption of a molecule on the surface of the coinage metals (Cu, Ag, and Au) is a typical example of systems where vdW interaction is significant. However, it is difficult to find a general method that reasonably describes both adsorption energy and geometry of even the simple prototypes of cyclic and heterocyclic aromatic molecules such as benzene (C₆H₆) and thiophene (C₄H₄S), respectively, with reasonable accuracy. In this work, we present an alternative scheme based on Zaremba-Kohn theory, called DFT+vdW-dZK. We show that unlike other popular methods, DFT+vdW-dZK and particularly SCAN+vdW-dZK give an accurate description of the physisorption of a rare-gas atom (xenon) and two small albeit diverse prototype organic molecules on the (111) surfaces of the coinage metals.

Theoretical understanding of the properties of stepped iron surfaces with van der Waals interaction corrections

The stepped surfaces in nanoscale zero-valent iron (nZVI) play an essential role for environmental application. However, there is still currently a deficiency in the atomic understanding of stepped surface properties due to the limitation of the computational methodology. In this study, stepped Fe(210) and (211) surfaces were theoretically investigated using density functional theory (DFT) computations in terms of the flat Fe(110) surface. Our results suggest that the consideration of van der Waals (vdW) interaction correction is beneficial for the DFT study on Fe-based systems. The DF-cx method is found to be the most promising vdW correction method. The DF-cx results reveal that the stepped Fe(210) and Fe(211) surfaces experience significant surface relaxation and abnormal trends in their work function. Their electronic properties and reactivities of the surface atoms are strongly affected by the Fe coordination numbers and the strong adsorption strengths of oxygen on the surfaces are dependent on both the coordination number of the adsorbed atoms and the geometry of the adsorption sites.

Assessing the effect of regularization on the molecular properties predicted by SCAN and self-interaction corrected SCAN meta-GGA

Recent regularization of the SCAN meta-GGA functional (rSCAN) has simplified the numerical complexities of the SCAN functional, alleviating SCAN's stringent demand on the numerical integration grids to some extent. The regularization of rSCAN, however, results in the breaking of some constraints such as the uniform electron gas limit, the slowly varying density limit, and coordinate scaling of the iso-orbital indicator. Here, we assess the effects of regularization on the electronic, structural, vibrational, and magnetic properties of molecules by comparing the SCAN and rSCAN predictions. The properties studied include atomic energies, atomization energies, ionization potentials, electron affinities, barrier heights, infrared intensities, dissociation and reaction energies, spin moments of molecular magnets, and isomer ordering of water clusters. Our results show that rSCAN requires less dense numerical grids and gives very similar results to those of SCAN for all properties examined with the exception of atomization energies, which are worsened in rSCAN. We also examine the performance of self-interaction-corrected (SIC) rSCAN with respect to SIC-SCAN using the Perdew-Zunger (PZ) SIC method. The PZSIC method uses orbital densities to compute one-electron self-interaction errors and places an even more stringent demand on numerical grids. Our results show that SIC-rSCAN gives marginally better performance than SIC-SCAN for almost all properties studied in this work with numerical grids that are on average half or less as dense as that needed for SIC-SCAN.

Improved description of hematite surfaces by the SCAN functional

Controversies on the surface termination of α-Fe₂O₃ (0001) focus on its surface stoichiometry dependence on the oxygen chemical potential. Density functional theory (DFT) calculations applying the commonly accepted Perdew-Burke-Ernzerhof (PBE) exchange-correlation functional to a strongly correlated system predict the best matching surface termination, but would produce a delocalization error, resulting in an inappropriate bandgap, and thus are not applicable for comprehensive hematite system studies. Besides, the widely applied PBE+U scheme cannot provide evidence for existence of some of the successfully synthesized stoichiometric α-Fe₂O₃ (0001) surfaces. Hence, a better scheme is needed for hematite DFT studies. This work investigates whether the strongly constrained and appropriately normed (SCAN) approximation reported by Perdew et al. could provide an improved result for the as-mentioned problem, and whether SCAN can be applied to hematite systems. By comparing the results calculated with the PBE, SCAN, PBE+U, and SCAN+U schemes, we find that SCAN and SCAN+U improves the description of the electronic structure of different stoichiometric α-Fe₂O₃ (0001) surfaces with respect to the PBE results, and that they give a consistent prediction of the surface terminations. Besides, the bulk lattice constants and the bulk density of states are also improved with the SCAN functional. This study provides a general characterization of the α-Fe₂O₃ (0001) surfaces and rationalizes how the SCAN approximation improves the results of hematite surface calculations.

Realization of Lieb lattice in covalent-organic frameworks with tunable topology and magnetism

Lieb lattice has been predicted to host various exotic electronic properties due to its unusual Dirac-flat band structure. However, the realization of a Lieb lattice in a real material is still unachievable. Based on tight-binding modeling, we find that the lattice distortion can significantly determine the electronic and topological properties of a Lieb lattice. Importantly, based on first-principles calculations, we predict that the two existing covalent organic frameworks (COFs), i.e., sp2C-COF and sp2N-COF, are actually the first two material realizations of organic-ligand-based Lieb lattice. Interestingly, the sp2C-COF can experience the phase transitions from a paramagnetic state to a ferromagnetic one and then to a Néel antiferromagnetic one, as the carrier doping concentration increases. Our findings not only confirm the first material realization of Lieb lattice in COFs, but also offer a possible way to achieve tunable topology and magnetism in organic lattices.

Localization in the SCAN meta-generalized gradient approximation functional leading to broken symmetry ground states for graphene and benzene

SCAN over localizes orbitals leading to spin symmetry broken ground states in graphene and benzene.

Regularized SCAN functional

We propose modifications to the functional form of the Strongly Constrained and Appropriately Normed (SCAN) density functional to eliminate numerical instabilities. This is necessary to allow reliable, automatic generation of pseudopotentials (including projector augmented-wave potentials). The regularized SCAN is designed to match the original form very closely, and we show that its performance remains comparable.

Prediction of intrinsic two dimensional ferromagnetism realized quantum anomalous Hall effect

The monolayer of FeX₃ (X = Cl, Br, I) possesses a quantum anomalous Hall insulating phase generated by the honeycomb lattice of iron atoms.