Semilocal exchange hole with an application to rangeseparation density functional
Abstract
Exchangecorrelation hole is a central concept in density functional theory. It not only provides justification for an exchangecorrelation energy functional, but also serves as a local ingredient in nonlocal rangeseparation density functional. However, due to the nonlocal nature, modelig the conventional exact exchange hole presents a great challenge to density functional theory. In this work, we propose a semilocal exchange hole underlying the TaoPerdewStaroverovScuseria (TPSS) metaGGA functional. The present model is distinct from previous models at small separation between an electron and the hole around the electron. It is also different in the way it interpolates between the rapidly varying isoorbital density and the slowly varying density, which is determined by the wave vector analysis based on the exactly solvable infinite barrier model for jellium surface. Our numerical tests show that the exchange hole generated from this model mimics the conventional exact exchange hole quite well for atoms. Finally, as a simple application, we apply the hole model to construct a TPSSbased rangeseparation functional. Our tests show that this TPSSbased rangeseparation functional can substantially improve TPSS band gaps and barrier heights, without losing much accuracy of molecular atomization energies.
I Introduction
KohnSham density functional theory (DFT) ks65 (); Parr89 (); RMDreizler90 () is a mainstream electronic structure theory, due to the useful accuracy and high computational efficiency. Formally it is exact, but in practice the exchangecorrelation energy component, which accounts for all manybody effects, has to be approximated as a functional of the electron density. Development of reliable exchangecorrelation energy functionals for a wide class of problems has been the central task of DFT. Many density functionals have been proposed PW86 (); B88 (); LYP88 (); BR89 (); B3PW91 (); B3LYP (); PBE96 (); VSXC98 (); HCTH (); PBE0 (); HSE03 (); AE05 (); MO6L (); TPSS03 (); revTPSS (); PBEsol (); SCAN15 (); Kaup14 (); TaoMo16 (), and some of them have achieved remarkable accuracy in condensed matter physics or quantum chemistry or both.
According to the local ingredients, density functionals can be classified into two broad categories: semilocal and nonlocal. Semilocal functionals make use of the local electron density, density derivatives, and/or the orbital kinetic energy density as inputs, such as the LSDA (local spindensity approximation) VWN80 (); PW92 (), GGA (generalizedgradient approximation) PBE96 (); PW91 (); ZWu06 () and metaGGA TPSS03 (); MO6L (); SCAN15 (); TaoMo16 (). Due to the simplicity in theoretical development, relatively easy numerical implementation, and cheap computational cost, semilocal functionals have been widelyused in electrunic structure calculations CJCramer09 (); Quest12 (); Yangreview (); Becke14 (). Indeed, semilocal DFT can give a quick and often accurate prediction of many properties such as atomization energies VNS03 (); PHao13 (); LGoerigk10 (); LGoerigk11 (), bond lengths CAdamo10 (); YMo16 (), lattice constants Csonka09 (); PBlaha09 (); FTran16 (), cohesive energies VNS04 (), etc.
Semilocal DFT has achieved high level of sophistication and practical success for many problems in chemistry, physics, and materials science, but it encounters difficulty in the prediction of reaction barrier heights, band gaps, charge transfer, and excitation energies. A proper description of these properties requires electronic nonlocality information PSTS08 (), which is absent in semilocal functionals. Nonlocality effect can be accounted for via mixing some amount of exact exchange into a semilocal DFT. This leads to the development of hybrid PBE0 (); B3PW91 (); VNS03 (); JJaramillo03 () and rangeseparation functionals HSE03 (); Truhlar11 (). The former involves the exact exchange energy or energy density, while the latter involves the exact as well as approximate semilocal exchange holes.
There are three general ways to develop the exchange hole: One is from paradigm systems such as the slowly varying density PW86 (); PBY96 (); Taobook10 () and oneelectron density BR89 (), another is from the reverse engineering approach EP98 (); tpsshole (); Lucian13 (), and the third is from the density matrix expansion TaoMo16 (). However, semilocal exchange holes developed with the reverse engineering approach may not be in the gauge of the conventional exchange hole. In the construction of the semilocal exchange hole, one must impose the hole to recover the underlying semilocal exchange energy density, which is usually not in the same gauge of the conventional exchange energy density, due to the integration by parts performed in the development of semilocal DFT. Examples include the PBE GGA EP98 () and TPSS metaGGAs tpsshole (); Lucian13 () exchange holes. Many rangeseparation functionals have been proposed LKronik11 (); RBaer10 (); JCTC16 (); Kresse11 (), and some of them have obtained great popularity in electronic structure calculations.
Semilocal exchange hole in the gauge of the conventional exchange is of special interest in DFT. In this work, we aim to develop an exchange hole, which reproduces the TPSS exchange functional. To ensure that our model hole is in the conventional gauge, we not only impose the exact constraints in the conventional gauge (e.g., correct shortrange behaviour without integration by parts) on the hole model, but also alter the TPSS exchange energy density by adding the TaoPerdewStaroverovScuseria gauge function tsspg08 () with a modification for improving the density tail behaviour of the original gauge function. The change in the local exchange energy density does not alter the integrated TPSS exchange energy, but it largely improves the agreement of the model hole with the conventional exact exchange hole. Furthermore, the hole model can generate the exact systemaveraged exchange hole by replacing the TPSS exchange energy density with the gaugecorrected TPSS exchange energy density but with the TPSS part replaced by the conventional exact exchange energy density. Finally, as a simple application, we apply our semilocal exchange hole to construct a rangeseparation functional. This TPSSbased rangeseparation functional yields band gaps of semiconductors and barrier heights in much better agreement with experimental values, without losing much accuracy of molecular atomization energies.
Ii Conventional exact exchange hole
For simplicity, we consider a spinunpolarized density (). For such a density, the exchange (x) energy can be written as
(1)  
where is the total electron density, is the conventional exchange energy per electron, or loosely speaking exchange energy density, and is the exchange hole at around an electron at . It is conventionally defined by
(2) 
Here is the KohnSham singleparticle density matrix defined by
(3) 
with being the occupied KohnSham orbitals. According to the expression (1), we can regard the exchange energy as the electrostatic interaction between reference electrons and their associated exchange hole. Therefore, an exchange energy functional cannot be fully justified unless the underlying exchange hole has been found.
The exchange hole for a spinunpolarized density can be generalized to any spin polarization with the spinscaling relation OP79 ()
(4) 
Therefore, in the development of the exchange hole, we only need to consider a spincompensated density. Performing the spherical average of the exchange hole over the direction of separation vector , the exchange part of Eq. (1) may be rewritten as
(5) 
where is the spherical average of the exchange hole defined by
(6) 
This suggests that the exchange energy does not depend on the detail of the associated hole. Rearranging Eq. (5) leads to a simple expression
(7) 
where is the system average of the exchange hole defined by
(8) 
Here is the number of electrons of a system.
Although the conventional exchange hole of Eq. (2) satisfies the sum rule,
(9) 
the most important property of the exchange hole, the exchange hole under a general coordinate transformation does not. Nevertheless, its system average always satisfies the sum rule
(10) 
This is a constraint that has been imposed in the development of semilocal exchange hole. While the exchange energy is uniquely defined, the exchange energy density and the exchange hole are not. For example, these local quantities can be altered by a general coordinate transformation, or arbitrarily adding the Laplacian of the electron density, without changing the exchange energy tsspg08 (); KBurke98 ().
Iii Constraints on the exchange hole
The conventional exchange hole is related to the pair distribution function by
(11) 
In general, a semilocal exchange hole can be written as
(12) 
where is the shape function that needs to be constructed, with being the dimensionless reduced density gradient, being the Fermi wave vector, , and . Here is the von Weizäscker kinetic energy density, and is the KohnSham orbital kinetic energy density given by
(13) 
Here are the occupied KohnSham orbitals with spin summed over.
iii.1 Constraints on the shape function
We will seek for a shape function that satisfies the following constraints:

Ontop value
(14) 
In the uniformgas limit,
(15) 
Normalization
(16) 
Negativity
(17) 
Energy constraint
(18) 
Small behaviour
(19) where will be discussed below. This is a constraint for the conventional exchange hole. That means no integration by parts is performed.

In the largegradient limit, the TPSS enhancement factor approaches PBE enhancement factor. In this limit, the TPSS shape function should also approach the PBE shape function, i.e.,
(20)
Among these constraints, (vi) is for the conventional exchange hole, while (vii) is a constraint used in the development of TPSS functional. These two constraints will be discussed in detail below. In previous works tpsshole (); Lucian13 (), constraint (vi) was used with integration by parts and thus is not for the conventional exchange hole, and (vii) was not considered.
iii.2 Small behaviour and largegradient limit
Expanding the sphericallyaveraged exchange hole up to second order in yields
(21)  
Since the Laplacian of the density tends to negative infinity at a nucleus, the negativity of the exchange hole in the small expansion of Eq. (21) is not ensured. Therefore, we must eliminate it. In previous works, the Laplacian of the density was eliminated by integration by parts tpsshole (). In order to model the conventional exchange hole, instead of performing integration by parts, here we eliminate it with the secondorder gradient expansion of the kinetic energy density,
(22) 
which is exact for slowly varying densities. This technique has been used in the development of TPSS TPSS03 () and other semilocal functionals PKZB (); revTPSS () as well as in the construction of electron localization indicator TLZR15 ().
Substituting Eqs. (21) into Eq. (12) and eliminating the Laplacian via (22) yields the small expansion of the shape function
(23)  
leading to
(24) 
For one or twoelectron densities, reduces to
(25) 
while for the uniform gas, . Note that .
In the largegradient limit, the TPSS shape function should recover the PBE shape function. This requires that must be merged smoothly with the PBE small behaviour,
(26) 
We can achieve this with
(27)  
where is the complementary error function defined by
(28) 
is a switching parameter that defines the point at which the small behaviour smoothly changes from TPSS to PBE. This choice of ensures that the small behaviour of our shape function is essentially determined by Eq. (23), while it merges into the PBE shape function in the largegradient limit.
Iv Shape function for the TPSS exchange hole
iv.1 TPSS shape function
The shape function for the TPSS exchange hole is assumed to take the following form
(29)  
where , , are determined by the conditions for the uniform electron gas, while the functions , and are determined by constraints (iii), (v) and (vi). They can be analytically expressed in terms of as
(30)  
(31)  
(32)  
where and . Following the procedure of Constantin, Perdew, and Tao tpsshole () for the construction of the original TPSS shape function, here we determine dependence of by fitting to the twoelectron hydrogenlike density, while dependence is determined by the wave vector analysis of the surface energy.
iv.2 Gradient dependence of
In soorbital regions where (e.g., core and density tail regions), we assume that the function takes the form of
(33) 
In the largegradient regime, the function should recover the PBE function Henderson08 ()
(34) 
For any density between the two regimes, we take the interpolation formula,
(35)  
Finally, we insert Eq. (35) into Eqs. (30)(32) and perform the fitting procedure by minimizing the following quantity
(36) 
where is the system average of the exchange hole defined by Eq. (8). It can be expressed in terms of the shape function as . For numerical convenience, we replace the integral with discretsized summation. All the parameters for and are listed in Table 1.
Figure 1 shows the systemaveraged exchange hole for the twoelectron exponential density evaluated with different hole models, compared to the exact one. We can observe from Fig. 1 that the present TPSS hole is slightly closer to the conventional exact hole than the original TPSS hole, but it is much closer than the PBE GGA and LSDA holes.
iv.3 Infinite barrier model and wave vector analysis for surface energy
For isoorbital regions, the dependence of is determined by fitting the model hole to the conventional exact exchange hole for the twoelectron exponential density. In the uniformgas limit, our exchange hole should correctly reduce to the LDA. This requires that vanishes in this limit. To fulfill these considerations, we assume that
(37) 
where is an integer. In order to determine , we follow the procedure of Ref. tpsshole () to study the wavevector analysis (WVA) of the surface energy. But instead of using the jellium surface model with linearly increasing barrier, here we employ the exactly solvable infinite barrier model (IBM). Since the singleparticle density matrix and hence the electron density of IBM is analytically known, this allows us to obtain insight into the dependence of from this model more easily.
Let us consider a uniform gas of noninteracting electrons subject to the infinite potential barrier perpendicular to the axis ( for ). The oneparticle density matrix is given by LMiglio81 (); IDMoore76 ()
(38)  
where is a step function, with for and for . is the average bulk valence electron density, , , , and
(39) 
with being the firstorder spherical Bessel function. The electron density can be obtained from the singleparticle density matrix by taking in Eq. (38). This yields
(40) 
The WVA for the surface exchange energy density is given by tpsshole ()
(41) 
where
(42) 
The exchange hole of IBM can be obtained from the oneparticle density matrix of Eq. (38). With some algebra, we can express the WVA surface exchange energy as Langreth82 ()
(43) 
where , and is given by
(44)  
and
(45) 
Here , with being defined by Eq. (3.18) of Ref. Langreth82 ().
Figure LABEL:Gamma shows of Eq. (44). The area under the curve is proportional to the surface exchange energy. From the electron density and density matrix of IBM given by Eqs. (38) and (40), the exact surface exchange energy can be calculated with the WVA of Eq. (43). Langreth and Perdew Langreth82 () reported that the value of is a.u., where is Seitz radius. This value is slightly smaller than the value obtained earlier by Harris and Jones Harris74 () and Ma and Sahni CQMa79 () ( a.u.). Our present work gives a.u., which might be the most accurate one.
iv.4 dependence of
The dependence of [Eq. (37)] can be determined by fitting the TPSS hole to the wave vector analysis. We start with the specific expressions for the local ingredients of the hole model in IBM.
From the electron density of Eq. (40), the reduced density gradient can be explicitly expressed as
(46) 
The kinetic energy density can be obtained from the singleparticle density matrix of Eq. (38). This yields
(47)  
Finally, the von Weizäscker kinetic energy density can be expressed as
(48) 
Next, we calculate from the TPSS hole. Starting with of Eq.(42), we obtain
(49)  
(Note that is implicit on the electron density.) Substituting Eq. (49) into Eq. (41), we obtain
(50) 
where
(51)  
and . Rearrangement of Eq. (50) leads to the final expression
(52) 
where
(53)  
(54) 
Figure 3 shows the comparison of with different choices of to the exact one. From Fig. 3, we see that the best choice is . Figure 4 shows that, compared to the LDA, PBE, and original TPSS holes, the WVA of the present model is in best agreement with the exact one. To further understand these two models, we plot the TPSS shape function of the present and the original models in IBM at , as shown by Figs. (5) and (6), respectively. From Figs. (5) and (6), we observe that while the present model hole is always negative, the original TPSS hole can be positive at some values of and .
To check our wave vector analysis for the surface exchange energy, we have computed from
(55) 
The results are shown in Table 2. From Table 2, we can see that the surface energy from the WVA of the PBE exchange hole shows some inconsistency with the surface energy evaluated directly from the PBE exchange functional with Eq. (55), but the surface energy from the WVA of the TPSS hole (both original and the present version) agrees very well with the surface energy calculated from TPSS exchange functional with Eq. (55). Furthermore, TPSS surface energy is more closer to the exact value than those of the LDA and PBE. The LDA significantly overestimates the surface exchange energy. All these observations are consistent with those evaluated from jellium surface model. It is interesting to note that even the original TPSS shape function in certain range is positive, the surface energy from the original TPSS hole is the same as that from the present model. This result is simply due to the cancellation of the original hole model between positive values and toonegative values at certain values, as seen from the comparison of Fig. (6) to (5). The IBM surface energy presents a great challenge to semilocal DFT. It is more difficult to get it right than the surface energy of jellium model with finite linear potential, because the electron density at the surface of IBM is highly inhomogeneous, due to immediate cutoff at surface, and is too far from the slowly varying regime, in which semilocal DFT can be exact (e.g., TPSS functional).
Eq .(55)  WVA integration  
LDA  6.318   
PBE  2.576  2.67 
TPSS  2.945  2.95 (original hole) 
2.95 (present hole) 
Figure 7 shows the comparison of the LSDA, PBE, and the original and present TPSS exchange hole models to the conventional exact exchange hole of the Ne atom, in which is in general different from 0 (slowly varying density) and 1 (isoorbital density). The PBE and LSDA curves are plotted with the hole models of Ref. Henderson08 (). From Fig. 7 we can see that the present TPSS hole model is generally closer to the exact one than the original TPSS hole model, while the original TPSS hole model is very close to PBE, but both of them are closer to the exact one than the LSDA.
V TPSS hole in the gauge of the conventional exact exchange
The shape function explicitly depends on the enhancement factor via the energy constraint of Eq. (18). The latter may be altered by arbitrarily adding any amount of the Laplacian of the density without changing the exchange energy. This ambiguity of the exchange energy density KBurke98 () leads to the ambiguity of the semilocal exchange hole. Our primary goal of this work is to develop a semilocal exchange hole in the gauge of the conventional exact exchange. This is partly motivated by the fact that in the development of rangeseparationals, the exact exchange part is usually provided in the conventional gauge.
The exact exchange energy density in the conventional gauge can be conveniently evaluated with the Della SalaGörling (DSG) DSG01 () identity resolution
(56) 
where is spin block of DSG matrix tsspg08 (). However, many semilocal exchange energy densities or enhancement factors of Eq. (18) are not in the gauge of the conventional exact exchange, due to the constraints such as the LiebOxford bound and the slowlyvarying gradient expansion imposed on the enhancement factor. For example, for the twoelectron exponential density, the conventionally defined exact enhancement factor is less than 1 near the nucleus, while the TPSS enhancement factor is by design. To construct the TPSS exchange hole in the conventional gauge, we can replace the original energy density constraint [Eq. (18)], as used in the construction of the original TPSS exchange hole tpsshole (), with the TPSS exchange energy density or enhancement factor in the conventional gauge. In this gauge, the TPSS exchange energy density can be written as tsspg08 ()
(57) 
where . Based on the uniform and nonuniform coordinate scaling properties of the exact exchange energy density, Tao, Staroverov, Scuseria, and Perdew (TSSP) tsspg08 () proposed a gauge function
(59)  
(60) 
Here and are determined by a fit to the conventional exact exchange energy density for the H atom, and is an integer which is chosen to be 4. is the exact exchange energy density in the conventional gauge. This gauge function is integrated to zero, i.e., , as required. It satisfies the correct uniform coordinate scaling relation, , and nonuniform scaling relation .
However, in the density tail () limit of an atom, the exact exchange energy density in the conventional gauge decays as , but the TSSP gauge function decays as . As a result, the exchange energy density in this gauge becomes positive. In order to fix this deficiency, we impose a constraint on the density tail,
(61) 
This can be achieved by requiring that in the limit, decays as with . Here we choose and take the same form of the TSSP gauge function, but with abd given by
(62)  
(63) 
Here and are determined by fitting the systemaveraged hole generated from the present TPSS shape function with the exact enhancement factor in the conventional gauge to the exact systemaveraged hole of the twoelectron exponential density. The fitting procedure is the same as that in the determination of the function. The parameter remains the same as that in the original version of Eq. (60).
Figure 8 shows the comparison of the present TPSS systemaveraged exchange hole and the systemaveraged exchange hole calculated from the present TPSS hole model with the conventional exact exchange energy density with and without gauge correction of Eqs. (62) and (63 to the exact one. From Fig. 8 we can observe that the exact systemaveraged exchange hole generated from the present TPSS hole model simply by replacing the TPSS enhancement factor with the conventional exact exchange enhancement factor without gauge correction significantly deviates from the exact systemaveraged hole. However, the agreement can be significantly improved with our gauge correction of Eqs. (62) and (63.
Figure 9 shows the comparison of the TPSS exchange energy density evaluated with the TPSS functional without and with gauge correction to the exact conventional exchange energy density for the twoelectron exponential density. From Fig. 9, we can observe that gauge correction is important, leading to a better agreement of the gaugecorrected TPSS exchange energy density with the exact one.
Vi Appication to rangeseparation density functional
As a simple application, we apply the present TPSS hole model to construct a rangeseparation functional. In general, there are two ways to construct rangeseparation functionals, depending on the need. For example, we may employ a semilocal DFT for the longrange part, while the exact exchange is used for the shortrange part, as pioneered by Heyd, Scuseria and Ernzerhof HSE03 (). This kind of rangeseparation functional is developed for solids, particularly useful for metallic solids, because usual hybrids requires much larger momentum cutoff for metallic systems with electrons nonlocalized. We may employ a semilocal DFT for the shortrange part, while the exact exchange is used for the longrange part, as developed by Henderson et al. Henderson08 () on the basis of the PBE hole. This kind of rangeseparation functionals are usually developed for molecular calculations, because the improved longrange part of the exchange hole greatly benefits the description of molecular properties, without significant increase of computational cost. Many rangeseparation functionals have been proposed hse06 (); Henderson07 (); AVKrukau08 (); Henderson09 (); GalliPRB14 (); Truhlar11 (). In the following, we will explore the TPSS holebased range separation functional with TPSS exchange functional being the longrange (LR) part and the exact exchange being the shortrange (SR) part, aiming to improve the TPSS description of twosmall band gaps and reaction barrier heights.
The idea of the construction of our TPSSbased rangeseparation functional is rooted in the constrcution of usual oneparameter hybrid functionals, which, in general, can be written as
(64) 
where is the mixing parameter that controls the amount of exact exchange mixed into a semilocal (sl) functional.
LSDA  PBE  TPSS  HSE  PW  Exp  
C  4.17  4.2  4.24  5.43  5.48  5.48 
CdSe  0.31  0.63  0.85  1.48  1.82  1.90 
GaAs  0.04  0.36  0.6  1.11  1.44  1.52 
GaN  2.15  2.22  2.18  3.48  3.5  3.50 
GaP  1.56  1.74  1.83  2.39  2.53  2.35 
Ge    0.13  0.32  0.8  0.99  0.74 
InAs  0    0.08  0.57  0.85  0.41 
InN  0  0  0  0.72  0.75  0.69 
InSb    0    0.47  0.73  0.23 
Si  0.53  0.62  0.71  1.2  1.31  1.17 
ZnS  2.02  2.3  2.53  3.44  3.78  3.66 
ME  0.14  
MAE  0.89  0.86  0.76  0.15  0.17 
LSDA  PBE  TPSS  TPSSh  HSE  PW  Expt  
SiH  347.4  313.2  333.7  333.6  314.5  333.6  322.4 
SiO  223.9  195.7  186.7  182.0  182.1  175.4  192.1 
S  135.1  114.8  108.7  105.9  106.3  101.9  101.7 
CH  802.1  721.2  707.5  704.4  705.9  699.9  704.8 
CHO  754.9  665.1  636.0  628.0  635.3  616.4  633.4 
CH  1304  1168  1156  1154  1152  1152  1149 
ME  77.4  12.4  4.1  0.75  
MAE  77.4  15.5  5.9  6.1  4.8  8.8 
Following the prescription of Heyd, Scuseria and Ernzerhof (HSE) HSE03 (), we write the TPSSbased rangeseparation functional as
(65) 
where is the HartreeFock (HF) exchange serving as part of the shortrange contribution, while is the TPSS exchange that provides the rest of the shortrange contribution. is TPSS correlation. The longrange contribution is provided fully by the TPSS exchange . They are given, respectively, by
(66) 
(67) 
(68) 
where is a rangeseparation parameter, and is the error function defined by Eq. (28). From Eqs. (65)(68), we can see that the amount of exact exchange mixing is controlled by two parameters and . Determination of them is discussed below. To test this functional, we have implemented it into the development version of Gaussian 09 gaussian09 ().
PBE  TPSS  TPSSh  HSE  PW  Ref  
OH + CH CH + HO  1.50  1.96  4.86  6.54(f)  
8.95  9.90  11.79  13.9  14.3  19.6(r)  
H + OH O + H  3.69  7.06  1.75  10.5(f)  
4.73  6.90  5.93  9.89  12.9(r)  
H + HS H + HS  1.03  3.55(f)  
9.40  12.72  13.4  12.4  14.4  17.3(r)  
ME  
MAE  9.37  8.34  6.76  4.66  4.63 
In the TPSSbased hybrid functional TPSSh, was fitted to molecular properties. In other words, in TPSSH, the optimized of is 0.1. So, in the TPSSbased rangeseparation functional of Eq. (65), the best value of will be 0, if is chosen. Since in the rangeseparation functional, some amount of the exact exchange (here the longrange part) in TPSSh is replaced by semilocal TPSS functional, to compensate for this, we need a value of larger than 0.1. Then we can find the best rangeseparation parameter by fitting to some electronic properties. To avoid possible overfitting, here we choose , a value that was recommended by Perdew, Ernzerhof, and Burke Perdew96 () and adopted with PBE0 functional PBE0 (). The parameter is determined by a fit to the band gap of diamond (C). This yields . Then we apply this rangeseparation functional to calculate the band gaps of 10 semiconductors. The results are listed in Table 3. From Table 3, we see that the band gaps of this rangeseparation functional is remarkably accurate, with a mean absolute deviation from experiments of only eV, about the same accuracy of HSE functional. We can also see from Table 3 that this TPSSbased rangeseparation functional should provide more accurate description for large bandgap materials, and therefore it provides an alternative choice for bandgap and other solidstate calculations.
Next, we apply our rangeseparation functional to calculate atomization energies of AE6 molecules. The results are listed in Tables 4. From Tables 4, we can see that, our rangeseparation functional only worsens the TPSS atomization energies of TPSS functional for this special set by about 3 kcal/mol. This error is still smaller than many other DFT methods such as LSDA, PBE, and PBEsol.
Reaction barrier heights are decisive quantity in the study of chemical kinetics. However, semilocal functionals tend to underestimate this quantity. As an interesting application, we apply our rangeseparation functional to calculate representative reaction barrier heights of BH6, which consists of three forward (f) and three reverse (r) barrier heights. The results are listed in Table 5. For comparison, we also calculated these barrier heights using PBE, TPSS, TPSSh, and HSE. From Table 5, we obserce that our rangeseparation functional provides the most accurate description of BH6, compared to other functionals considered.
Vii Conclusion
In conclusion, we have developed a semilocal exchange hole underlying TPSS exchange functional. The hole is exact in the uniformgas limit and accurate for isoorital densities. Then is is interpolated between these two densities through WVA on the surface of infinite barrier model. Our numerical test on the Ne atom shows that the hole mimics the conventional exact exchange hole quite accurately. In particular, with our new gauge function correction, the hole model can generate the systemaveraged hole from the conventional exact exchange energy density in good agreement with the present TPSS hole and the exact hole for the twoelectron exponential density.
Based on the TPSS hole, we have constructed a rangeseparation functional. This functional yields very accurate band gaps and reaction barrier heights, without losing much accuracy on atomization energies. This accuracy makes this functional an attractive alternative in solidstate calculations and molecular kinetics.
References
 (1) W. Kohn and L.J. Sham, Phys. Rev. 140, A1133 (1965).
 (2) R.G. Parr and W. Yang, Density Functional Theory of Atoms and Molecules (Oxford University Press, New York, 1989).
 (3) R.M. Dreizler and E.K.U. Gross, Density Functional Theory (Springer, Berlin, 1990).
 (4) J.P. Perdew and Y. Wang, Phys. Rev. B 33, 8800 (1986).
 (5) A.D. Becke, Phys. Rev. A 38, 3098 (1988).
 (6) C. Lee, W. Yang, and R.G. Parr, Phys. Rev. B 37, 785 (1988).
 (7) A.D. Becke and M. R. Roussel, Phys. Rev. A 39, 3761 (1989).
 (8) A.D. Becke, J. Chem. Phys. 104, 1040 (1996).
 (9) P.J. Stephens, F.J. Devlin, C.F. Chabalowski, and M.J. Frisch, J. Phys. Chem. 98, 11623 (1994).
 (10) J.P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. 77, 3865 (1996).
 (11) T.V. Voorhis and G.E. Scuseria, J. Chem. Phys. 109, 400 (1998).
 (12) F.A. Hamprecht, A.J. Cohen, D.J. Tozer, and N.C. Handy, J. Chem. Phys. 109, 6264 (1998).
 (13) M. Ernzerhof and G.E. Scuseria, J. Chem. Phys. 110, 5029 (1999).
 (14) J. Heyd, G.E. Scuseria, and M. Ernzerhof, J. Chem. Phys. 118, 8207 (2003).
 (15) R. Armiento and A.E. Mattsson, Phys. Rev. B 72, 085108 (2005).
 (16) Y. Zhao and D.G. Truhlar, J. Chem. Phys. 125, 194101 (2006).
 (17) J. Tao, J.P. Perdew, V.N. Staroverov, and G.E. Scuseria, Phys. Rev. Lett. 91, 146401 (2003).
 (18) J.P. Perdew, A. Ruzsinszky, G.I. Csonka, L.A. Constantin, and J. Sun, Phys. Rev. Lett. 103, 026403 (2009).
 (19) J.P. Perdew, A. Ruzsinszky, G.I. Csonka, O.A. Vydrov, G.E. Scuseria, L.A. Constantin, X. Zhou, and K. Burke, Phys. Rev. Lett. 100, 136406 (2008).
 (20) J. Sun, A. Ruzsinszky, and J.P. Perdew, Phys. Rev. Lett. 115, 036402 (2015).
 (21) A.V. Arbuznikov and M. Kaupp, J. Chem. Phys. 141, 204101 (2014).
 (22) J. Tao and Y. Mo, Phys. Rev. Lett. 117, 073001 (2016).
 (23) S.H. Vosko, L. Wilk, and M. Nusair, Can. J. Phys. 58, 1200 (1980).
 (24) J.P. Perdew and Y. Wang, Phys. Rev. B 45, 13244 (1992).
 (25) J.P. Perdew, J.A. Chevary, S.H. Vosko, K.A. Jackson, M.R. Pederson, D.J. Singh, and C. Fiolhais, Phys. Rev. B 46, 6671 (1992).
 (26) Z. Wu and R.E. Cohen, Phys. Rev. B 73, 235116 (2006).
 (27) C.J. Cramer and D.G. Truhlar, Phys.Chem.Chem.Phys. 11, 10757 (2009).
 (28) R. Peverati and D.G. Truhlar, Phys. Phil.Trans. R. Soc. A 372, 20120476 (2011).
 (29) A.J. Cohen, P. MoriSanchez, and W. Yang, Chem.Rev. 112, 289 (2012).
 (30) A.D. Becke, J. Chem. Phys. 140, 18A301 (2014).
 (31) V.N. Staroverov, G.E. Scuseria, J. Tao, and J.P. Perdew, J. Chem. Phys. 119, 12129 (2003); 121, 11507(E) (2004).
 (32) P. Hao, J. Sun, B. Xiao, A. Ruzsinszky, G.I. Csonka, J. Tao, S. Glindmeyer, and J.P. Perdew, J. Chem. Theory Comput. 9, 355 (2013).
 (33) L. Goerigk and S. Grimme, J. Chem. Theory Comput. 6, 107 (2010).
 (34) L. Goerigk and S. Grimme, J. Chem. Theory Comput. 7, 291 (2011).
 (35) D. Jacquemin and C. Adamo, J. Chem. Theory Comput. 7, 369 (2011).
 (36) Y. Mo, G. Tian, R. Car, V.N. Staroverov, G.E. Scuseria, and J. Tao, submitted.
 (37) G.I. Csonka, J.P. Perdew, A. Ruzsinszky, P.H.T. Philipsen, S. Lebégue, J. Paier, O.A. Vydrov, and J.G. Ángyán Phys. Rev. B 79, 155107 (2009).
 (38) P. Haas, F. Tran, and P. Blaha, Phys. Rev. B 79, 085104 (2009).
 (39) F. Tran, J. Stelzl, and P. Blaha, J. Chem. Phys. 144, 204120 (2016).
 (40) V.N. Staroverov, G.E. Scuseria, J. Tao, and J.P. Perdew, Phys. Rev. B 69, 075102 (2004).
 (41) J.P. Perdew, V.N. Staroverov, J. Tao, and G.E. Scuseria, Phys. Rev. A 78, 052513 (2008).
 (42) J. Jaramillo, G.E. Scuseria, and M. Ernzerhof, J. Chem. Phys. 118, 1068 (2003).
 (43) R. Peverati and D.G. Truhlar, Phys. Chem. Lett. 2, 2810 (2011).
 (44) Phys. Rev. B 54, 16533 (1996).
 (45) J. Tao, Density Functional Theory of Atoms, Molecules, and Solids (VDM Verlag, Germany, 2010).
 (46) M. Ernzerhof and J.P. Perdew J. Chem. Phys. 109, 3313 (1998).
 (47) L.A. Constantin, J.P. Perdew, and J. Tao, Phys. Rev. B 73, 205104 (2006).
 (48) L.A. Constantin, E. Fabiano, and F. Della Sala Phys. Rev. B 88, 125112 (2013).
 (49) Phys. Rev. B 84, 075144 (2011).
 (50) Ann. Rev. Phys. Chem. 61, 85 (2010).
 (51) M. Modrzejewski, M. Hapka, G. Chalasinski, and M.M. Szczesniak, J. Chem. Theory Comput., 12, 3662 (2016).
 (52) L. Schimka, J. Harl, and G. Kresse, J. Chem. Phys. 134, 024116 (2011).
 (53) J. Tao, J.P. Perdew, V.N. Staroverov, and G.E. Scuseria, Phys. Rev. A 77, 012509 (2008).
 (54) G.L. Oliver and J.P. Perdew, Phys. Rev. A 20, 397 (1979).
 (55) K. Burke, F.G. Cruz, and K.C. Lam, J. Chem. Phys. 109, 8161 (1998).
 (56) J.P. Perdew, S. Kurth, A. Zupan, and P. Blaha, Phys. Rev. Lett. 82, 2544 (1999).
 (57) J. Tao, S. Liu, Z. Fan, and A.M. Rappe, Phys. Rev. B 92, 060401(R) (2015).
 (58) T.M. Henderson, B.G. Janesko, and G.E. Scuseria, J. Chem. Phys. 128, 194105 (2008).
 (59) L. Miglio, M.P. Tosi, and N.H. March, Sur. Sci. 111, 119 (1981).
 (60) I.D. Moore and N.H. March, Ann. Phys. 97, 136 (1976).
 (61) D.C. Langreth and J.P. Perdew, Phys. Rev. B 26, 2810 (1982).
 (62) J. Harris and R.O. Jones, J. Phys. F 4, 1170 (1974).
 (63) C.Q. Ma and V. Sahni, Phys. Rev. B 20, 2291 (1979).
 (64) F. Della Sala and A. Görling, J. Chem. Phys. 115, 5718 (2001).
 (65) A.V. Krukau, O.A. Vydrov, A.F. Izmaylov, and G.E. Scuseria, J. Chem. Phys. 125, 224106 (2006).
 (66) T.M. Henderson, A.F. Izmaylov, G.E. Scuseria, and A. Savin, J. Chem. Phys. 127, 221103 (2007).
 (67) A.V. Krukau, G.E. Scuseria, J.P. Perdew, and A. Savin, J. Chem. Phys. 129, 124103 (2008).
 (68) T.M. Henderson, B.G. Janesko, G.E. Scuseria, and A. Savin, Int. J. Quantum Chem. 109, 2023 (2009).
 (69) J.H. Skone, M. Govoni, and G. Galli, Phys. Rev. B 89, 195112 (2014).
 (70) Y. Zhao, J. Pu, B.J. Lynch, and D.G. Truhlar, Phys. Chem. Chem. Phys. 6, 673 (2004).
 (71) J. Paier, B.G. Janesko, T.M. Henderson, G.E. Scuseria, A. Güneis, and G. Kresse, J. Chem. Phys. 132, 094103 (2010).
 (72) B.J. Lynch and D.G. Truhlar, J. Phys. Chem. A 107, 8996 (2003).
 (73) Gaussian Development Version, Revision B.01, M.J. Frisch et al., Gaussian, Inc., Wallingford CT (2009).
 (74) J.P. Perdew, M. Ernzerhof, and K. Burke, J. Chem. Phys. 105, 9982 (1996).