Assessing the performance of the Tao-Mo semilocal density functional in the projector-augmented-wave method

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.

Performance of Density Functionals for Excited-State Properties of Isolated Chromophores and Exciplexes: Emission Spectra, Solvatochromic Shifts, and Charge-Transfer Character

This study assesses the performance of various meta-generalized gradient approximation (meta-GGA), global hybrid, and range-separated hybrid (RSH) density functionals in capturing the excited-state properties of organic chromophores and their excited-state complexes (exciplexes). Motivated by their uses in solar energy harvesting and photoredox CO2 reduction, we use oligo-(p-phenylenes) and their excited-state complexes with triethylamine as model systems. We focus on the fluorescence properties of these systems, specifically emission energies. We also consider solvatochromic shifts and wave function characteristics. The latter is described by using reduced quantities such as natural transition orbitals (NTOs) and exciton descriptors. The functionals are benchmarked against the experimental fluorescence spectra and the equation-of-motion coupled-cluster method with single and double excitations. Both in isolated chromophores and in exciplexes, meta-GGA functionals drastically underestimate the emission energies and exhibit significant exciton delocalization and anticorrelation between electron and hole motion. The performance of global hybrid functionals is strongly dependent on the percentage of exact exchange. Our study identifies RSH GGAs as the best-performing functionals, with ωPBE demonstrating the best agreement with experimental results. RSH meta-GGAs often overestimate emission energies in exciplexes and yield larger hole NTOs. Their performance can be improved by optimally tuning the range-separation parameter.

Density functional theory methods applied to homogeneous and heterogeneous catalysis: a short review and a practical user guide

The application of density functional theory (DFT) methods in catalysis has been growing fast in the last few decades thanks to both the availability of more powerful high computing resources and the development of new efficient approximations and approaches. DFT calculations allow for the understanding of crucial catalytic aspects that are difficult or even impossible to access by experiments, thus contributing to faster development of more efficient and selective catalysts. Depending on the catalytic system and properties under investigation, different approaches should be used. Moreover, the reliability of the obtained results deeply depends on the approximations involved in both the selected method and model. This review addresses chemists, physicists and materials scientists whose interest deals with the application of DFT-based computational tools in both homogeneous catalysis and heterogeneous catalysis. First, a brief introduction to DFT is presented. Then, the main approaches based on atomic centered basis sets and plane waves are discussed, underlining the main differences, advantages and limitations. Eventually, guidance towards the selection of the catalytic model is given, with a final focus on the evaluation of the energy barriers, which represents a crucial step in all catalytic processes. Overall, the review represents a rational and practical guide for both beginners and more experienced users involved in the wide field of catalysis.

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.

First-principle investigation of structural, electronic, and phase stabilities in chalcopyrite semiconductors: insights from Meta-GGA functionals

We undertake a comprehensive first-principles investigation into the factors influencing the optoelectronic efficiencies of PIQIIIR2VIchalcopyrite semiconductors. The structural attributes, electronic properties, and phase stabilities are explored using various meta-GGA exchange-correlation (XC) functionals within the density functional framework. In particular, we assess the relative performance of these XC functionals in obtaining estimates of various relevant parameters. The structural parameteruin chalcopyrite semiconductors is a noteworthy aspect, as it is intrinsically tied to the extent of orbital hybridization between distinct atoms and thereby strongly influences the electronic properties. In general, the application of widely used GGA-PBE XC functional to these chalcopyrites results in unreliable predictions of band gaps and 'u' parameter due to delocalization errors that in turn arise due to the inclusion ofdandfcore electrons. While hybrid functionals offer remarkable accuracy through state-of-the-art methods, their main drawback lies in their computational expense and resource demands. Our findings strongly suggest that in comparison to GGA-PBE, the meta-GGA XC functionals perform quite well and provide results that closely align with experimental values. In particular, ther2SCAN andrMGGAC XC functionals are preferable and superior for investigating chalcopyrites and other solid-state systems. This preference is rooted in their excellent performance and substantially reduced computational costs compared to hybrid functionals.

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.

Density functional applications of jellium with a local gap model correlation energy functional

We develop a realistic density functional approximation for the local gap, which is based on a semilocal indicator that shows good screening properties. The local band model has remarkable density scaling behaviors and works properly for the helium isoelectronic series for the atoms of the Periodic Table, as well as for the non-relativistic noble atom series (up to 2022 e-). Due to these desirable properties, we implement the local gap model in the jellium-with-gap correlation energy, developing the local-density-approximation-with-gap correlation functional (named LDAg) that correctly gives correlation energies of atoms comparable with the LDA ones but shows an improvement for ionization potential of atoms and molecules. Thus, LDAg seems to be an interesting and useful tool in density functional theory.

Efficient and improved prediction of the band offsets at semiconductorheterojunctions from meta-GGA density functionals: a benchmark study

Accurate theoretical prediction of the band offsets at interfaces of semiconductor heterostructures can of-ten be quite challenging. Although density functional theory has been reasonably successful to carry outsuch calculations and efficient and accurate semilocal functionals are desirable to reduce the computational cost. In general, the semilocal functionals based on the generalized gradient approximation (GGA) significantly underestimate the bulk band gaps. This, in turn, results in inaccurate estimates of the band offsets at the heterointerfaces. In this paper, we investigate the performance of several advanced meta-GGA functionals in the computational prediction of band offsets at semiconductor heterojunctions. In particular, we investigate the performance of r 2 SCAN (revised strongly-constrained and appropriately-normed functional), rMGGAC (revised semilocal functional based on cuspless hydrogen model and Pauli kinetic energy density functional), mTASK (modified Aschebrock and Kümmel meta-GGA functional), and LMBJ (local modified Becke-Johnson) exchange-correlation functionals. Our results strongly suggest that these meta-GGA functionals for supercell calculations perform quite well, especially, when compared to computationally more demanding GW calculations. We also present band offsets calculated using ionization potentials and electron affinities, as well as band alignment via the branch point energies. Overall, our study shows that the aforementioned meta-GGA functionals can be used within the DFT framework to estimate the band offsets in semiconductor heterostructures with predictive accuracy.

A Local Hybrid Exchange Functional Approximation from First Principles

Local hybrid functionals are a more flexible class of density functional approximations allowing for a position-dependent admixture of exact exchange. This additional flexibility, however, comes with a more involved mathematical form and a more complicated design. A common denominator for previously constructed local hybrid funtionals is usage of thermochemical benchmark data to construct these functionals. Herein, we design a local hybrid functional without relying on benchmark data. Instead, we construct it in a more ab initio manner, following the principles of modern meta-generalized gradient approximations and considering theoretical constrains. To achieve this, we make use of the density matrix expansion and a local mixing function based on an approximate correlation length. The accuracy of the developed density functional approximation is assessed for thermochemistry, excitation energies, polarizabilities, magnetizabilities, NMR spin-spincoupling constants, NMR shieldings and shifts, as well as EPR g-tensors and hyperfine coupling constants. Here, the new exchange functional shows a robust performance and is especially well suited for atomization energies, barrier heights, excitation energies, NMR coupling constants, and EPR properties, whereas it looses some ground for the NMR shifts.Therefore, the designed functional is a major step forwards for functionals that have been designed from first principles.

Solid-state performance of a meta-GGA screened hybrid density functional constructed from Pauli kinetic enhancement factor dependent semilocal exchange hole

The semilocal form of the exchange hole is highly useful in developing non-local range-separated hybrid density functionals for finite and extended systems. The way to construct the conventional exact exchange hole model is based on either the Taylor series expansion or the reverse engineering technique from the corresponding exchange energy functional. Although the latter technique is quite popular in context of generalized gradient approximation (GGA) functionals, the same for the meta-GGA functionals is not so much explored. Thus, in this study, we propose a reverse-engineered semilocal exchange hole of a meta-GGA functional, that depends only on the meta-GGA ingredient α (also known as the Pauli kinetic energy enhancement factor). The model is used subsequently to design the short-range-separated meta-GGA hybrid density functional. We show that the present method can be successfully applied for several challenging problems in the context of solids, especially for which the GGA based hybrid fails drastically. This assessment proves that the present functional is quite useful for materials sciences. Finally, we also use this method for several molecular test cases, where the results are also as comparative as its base semilocal functional.

Accurate band gaps from exchange potentials designed from a cuspless hydrogen density-based exchange hole model

The explicit forms of exchange-correlation (XC) potentials, which are not functional derivatives of any XC energy functional, are reasonably efficient in predicting the band gap of materials. The most successful example in this genre is the MBJ [F. Tran et al., Phys. Rev. Lett., 2009 102, 226401] exchange potential, which is based on the asymptotically correct Becke-Roussel (BR) exchange potential. We employ the cuspless hydrogen density and corresponding exchange hole to construct a BR like potential. The modified BR potential is again utilized in the framework of MBJ for band gap calculations. Also, we employ a Laplacian free model of the exchange hole in the framework of MBJ. These methods are analyzed by calculating band gaps of various test sets containing narrow, intermediate, and wide bandgap materials. Besides, we apply these potentials to eighteen ternary semiconductors with a chalcopyrite structure, exciting materials for photovoltaic applications. By comparing them with MBJ, we find that the band gaps obtained using the new potentials are not uniformly larger values than the MBJ potential in all cases. But, in many instances where MBJ overestimates the gap, the new potentials' band gaps are comparatively smaller and closer to the experimental ones. We also show that these potentials can correctly predict the band structure of three-dimensional topological insulators.

Improved electronic structure prediction of chalcopyrite semiconductors from a semilocal density functional based on Pauli kinetic energy enhancement factor

The correct treatment ofdelectrons is of prime importance in order to predict the electronic properties of the prototype chalcopyrite semiconductors. The effect ofdstates is linked with the anion displacement parameteru, which in turn influences the bandgap of these systems. Semilocal exchange-correlation functionals which yield good structural properties of semiconductors and insulators often fail to predict reasonableubecause of the underestimation of the bandgaps arising from the strong interplay betweendelectrons. In the present study, we show that the meta-generalized gradient approximation (meta-GGA) obtained from the cuspless hydrogen density (MGGAC) (2019Phys. Rev.B 100 155140) performs in an improved manner in apprehending the key features of the electronic properties of chalcopyrites, and its bandgaps are comparative to that obtained using state-of-art hybrid methods. Moreover, the present assessment also shows the importance of the Pauli kinetic energy enhancement factor,α= (τ-τ^W)/τ^unifin describing thedelectrons in chalcopyrites. The present study strongly suggests that the MGGAC functional within semilocal approximations can be a better and preferred choice to study the chalcopyrites and other solid-state systems due to its superior performance and significantly low computational cost.

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.

Efficient Band Structure Calculation of Two-Dimensional Materials from Semilocal Density Functionals

The experimental and theoretical realization of two-dimensional (2D) materials is of utmost importance in semiconducting applications. Computational modeling of these systems with satisfactory accuracy and computational efficiency is only feasible with semilocal density functional theory methods. In the search for the most useful method in predicting the band gap of 2D materials, we assess the accuracy of recently developed semilocal exchange-correlation (XC) energy functionals and potentials. Though the explicit forms of exchange-correlation (XC) potentials are very effective against XC energy functionals for the band gap of bulk solids, their performance needs to be investigated for 2D materials. In particular, the LMBJ [J. Chem. Theory Comput.2020, 16, 2654] and GLLB-SC [Phys. Rev. B82, 2010, 115106] potentials are considered for their dominance in bulk band gap calculation. The performance of recently developed meta generalized gradient approximations, like TASK [Phys. Rev. Res.1, 2019, 033082] and MGGAC [Phys. Rev. B. 100, 2019, 155140], is also assessed. We find that the LMBJ potential constructed for 2D materials is not as successful as its parent functional, i.e., MBJ [Phys. Rev. Lett.102, 2009, 226401] in bulk solids. Due to a contribution from the derivative discontinuity, the band gaps obtained with GLLB-SC are in a certain number of cases, albeit not systematically, larger than those obtained with other methods, which leads to better agreement with the quasi-particle band gap obtained from the GW method. The band gaps obtained with TASK and MGGAC can also be quite accurate.

Generalizing Double-Hybrid Density Functionals: Impact of Higher-Order Perturbation Terms

Connections between the Görling-Levy (GL) perturbation theory and the parameters of double-hybrid (DH) density functional are established via adiabatic connection formalism. Moreover, we present a more general DH density functional theory, where the higher-order perturbation terms beyond the second-order GL2 one, such as GL3 and GL4, also contribute. It is shown that a class of DH functionals including previously proposed ones can be formed using the present construction. Based on the proposed formalism, we assess the performance of higher-order DH and long-range corrected DH formed on the Perdew-Burke-Ernzerhof (PBE) semilocal functional and second-order GL2 correlation energy. The underlying construction of DH functionals based on the generalized many-body perturbation approaches is physically appealing in terms of the development of the non-local forms using more accurate and sophisticated semilocal functionals.

Improved transition metal surface energies from a generalized gradient approximation developed for quasi two-dimensional systems

Nonuniform density scaling in the quasi-two-dimensional (quasi-2D) regime is an important and challenging aspect of the density functional theory. Semilocal exchange-correlation energy functionals, developed by solving the dimensional crossover criterion in the quasi-2D regime, have great theoretical and practical importance. However, the only semilocal generalized gradient approximation (GGA) that has been designed to satisfy this criterion is the Q2D-GGA [L. Chiodo et al., Phys. Rev. Lett. 108, 126402 (2012)]. Here, we establish the applicability, broadness, and accuracy of the Q2D-GGA functional by performing an extensive assessment of this functional for transition metal surface energies. The important characteristic of the surface density localization and oscillation due to the rearrangement of the d electrons is also shown for different metal surfaces.

Screened range-separated hybrid by balancing the compact and slowly varying density regimes: Satisfaction of local density linear response

Due to their quantitative accuracy and ability to solve several difficulties, screened range-separated hybrid exchange-correlation functionals are now a standard approach for ab initio simulation of condensed matter systems. However, the screened range-separated hybrid functionals proposed so far are biased either toward compact or slowly varying densities. In this paper, we propose a screened range-separated hybrid functional, named HSEint, which can well describe these density regimes, achieving good accuracy for both molecular and solid-state systems. The semilocal part of the proposed functional is based on the PBEint generalized gradient approximation [E. Fabiano et al., Phys. Rev. B 82, 113104 (2010)], constructed for hybrid interfaces. To improve the functional performance, we employ exact or nearly exact constraints in the construction of range-separated hybrid functional, such as recovering of the local density linear response and semiclassical atom linear response.