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

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.

Semilocal Meta-GGA Exchange-Correlation Approximation from Adiabatic Connection Formalism: Extent and Limitations

The incorporation of a strong-interaction regime within the approximate semilocal exchange-correlation functionals still remains a very challenging task for density functional theory. One of the promising attempts in this direction is the recently proposed adiabatic connection semilocal correlation (ACSC) approach [Constantin, L. A.; Phys. Rev. B 2019, 99, 085117] allowing one to construct the correlation energy functionals by interpolation of the high and low-density limits for the given semilocal approximation. The current study extends the ACSC method to the meta-generalized gradient approximations (meta-GGA) level of theory, providing some new insights in this context. As an example, we construct the correlation energy functional on the basis of the high- and low-density limits of the Tao-Perdew-Staroverov-Scuseria (TPSS) functional. Arose in this way, the TPSS-ACSC functional is one-electron self-interaction free and accurate for the strictly correlated and quasi-two-dimensional regimes. Based on simple examples, we show the advantages and disadvantages of ACSC semilocal functionals and provide some new guidelines for future developments in this context.

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.

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.

Benchmark test of a dispersion corrected revised Tao-Mo semilocal functional for thermochemistry, kinetics, and noncovalent interactions of molecules and solids

In the density functional theory, dispersion corrected semilocal approximations are often used to benchmark weekly interacting finite and extended systems. Here, the focus is on providing a broad overview of the performance of D3 dispersion corrected revised Tao-Mo (revTM) semilocal functionals [A. Patra et al., J. Chem. Phys. 153, 084 117 (2020)] for thermochemistry and kinetics of molecules, molecular crystals, ice polymorphs, metal-organic systems, atom/molecular adsorption on solids, water interacting with nano-materials, binding energies of layered materials, and properties of weekly and strongly bonded solids. We show that the most suitable "optimized power" function for the revTM functional needs a modification to make it suitable for properties related to the diverse nature of finite and extended systems. The present work is an extension of the previously proposed revTM+D3 method with the motivation to design and benchmark the dispersion corrected cost-effective method based on this semilocal approximation. We show that the revised revTM+D3 functional provides various general purpose molecular and solid properties with the closest to experimental findings than its predecessor. The present assessment and benchmarking can be practically useful for performing cost-effective method based simulations of various molecular and solid-state properties.

Bandgap of two-dimensional materials: Thorough assessment of modern exchange-correlation functionals

The density-functional theory (DFT) approximations that are the most accurate for the calculation of bandgap of bulk materials are hybrid functionals, such as HSE06, the modified Becke-Johnson (MBJ) potential, and the GLLB-SC potential. More recently, generalized gradient approximations (GGAs), such as HLE16, or meta-GGAs, such as (m)TASK, have also proven to be quite accurate for the bandgap. Here, the focus is on two-dimensional (2D) materials and the goal is to provide a broad overview of the performance of DFT functionals by considering a large test set of 298 2D systems. The present work is an extension of our recent studies [T. Rauch, M. A. L. Marques, and S. Botti, Phys. Rev. B 101, 245163 (2020); Patra et al., J. Phys. Chem. C 125, 11206 (2021)]. Due to the lack of experimental results for the bandgap of 2D systems, G₀W₀ results were taken as reference. It is shown that the GLLB-SC potential and mTASK functional provide the bandgaps that are the closest to G₀W₀. Following closely, the local MBJ potential has a pretty good accuracy that is similar to the accuracy of the more expensive hybrid functional HSE06.

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.

Density Functional Theory for Molecule-Metal Surface Reactions: When Does the Generalized Gradient Approximation Get It Right, and What to Do If It Does Not

While density functional theory (DFT) is perhaps the most used electronic structure theory in chemistry, many of its practical aspects remain poorly understood. For instance, DFT at the generalized gradient approximation (GGA) tends to fail miserably at describing gas-phase reaction barriers, while it performs surprisingly well for many molecule-metal surface reactions. GGA-DFT also fails for many systems in the latter category, and up to now it has not been clear when one may expect it to work. We show that GGA-DFT tends to work if the difference between the work function of the metal and the molecule's electron affinity is greater than ∼7 eV and to fail if this difference is smaller, with sticking of O₂ on Al(111) being a spectacular example. Using dynamics calculations we show that, for this system, the DFT problem may be solved as done for gas-phase reactions, i.e., by resorting to hybrid functionals, but using screening at long-range to obtain a correct description of the metal. Our results suggest the GGA error in the O₂ + Al(111) barrier height to be functional driven. Our results also suggest the possibility to compute potential energy surfaces for the difficult-to-treat systems with computationally cheap nonself-consistent calculations in which a hybrid functional is applied to a GGA density.

Electron-electron correlations and structural, spectral and polarization properties of tetragonal BaTiO3

To analyze the role of electron-electron correlation effects in structural (local-geometry), spectral and polarization properties of tetragonal BaTiO3we apply DFT +Uapproach. We demonstrate that the system properties drastically change when the value of the local Coulomb repulsionUcrosses the critical valueUc≈ 7 eV. In particular, the correlation effects cause a change of the ratio of the in-plane and inter-plane Ti-O bond lengths, which results in a flip of the order of the Tid-bands and change of the polarizability of the system. Since the consensus value ofUin BaTiO3is unknown, we discuss how the obtained results may be revealed in experimental data, especially in the optical response and ultrafast charge dynamics, where effectiveUis dynamically tuned.

Unveiling the Physics Behind Hybrid Functionals

We show that accurate exchange-correlation hybrid functionals give very physically optimized effective-correlation potentials, capable of correctly describing the quantum oscillations of atoms and molecules. Based on this analysis and on understanding the error cancellation between semilocal exchange and correlation functionals, we propose a very simple, semilocal correlation potential model compatible with the exact exchange of density functional theory, which performs remarkably well for charge densities and orbital energies.