# Basis convergence of range-separated density-functional theory

###### Abstract

Range-separated density-functional theory is an alternative approach to Kohn-Sham density-functional theory. The strategy of range-separated density-functional theory consists in separating the Coulomb electron-electron interaction into long-range and short-range components, and treating the long-range part by an explicit many-body wave-function method and the short-range part by a density-functional approximation. Among the advantages of using many-body methods for the long-range part of the electron-electron interaction is that they are much less sensitive to the one-electron atomic basis compared to the case of the standard Coulomb interaction. Here, we provide a detailed study of the basis convergence of range-separated density-functional theory. We study the convergence of the partial-wave expansion of the long-range wave function near the electron-electron coalescence. We show that the rate of convergence is exponential with respect to the maximal angular momentum for the long-range wave function, whereas it is polynomial for the case of the Coulomb interaction. We also study the convergence of the long-range second-order Møller-Plesset correlation energy of four systems (He, Ne, N, and HO) with the cardinal number of the Dunning basis sets cc-p(C)VZ, and find that the error in the correlation energy is best fitted by an exponential in . This leads us to propose a three-point complete-basis-set extrapolation scheme for range-separated density-functional theory based on an exponential formula.

^{†}

^{†}preprint:

## I Introduction

Range-separated density-functional theory (DFT) (see, e.g., Ref. TouColSav-PRA-04, ) is an attractive approach for improving the accuracy of Kohn-Sham DFT HohKoh-PR-64 (); KohSha-PR-65 () applied with usual local or semi-local density-functional approximations. This approach is particularly relevant for the treatment of electronic systems with strong (static) or weak (van der Waals) correlation effects. The strategy of range-separated DFT consists in separating the Coulomb electron-electron interaction into long-range and short-range components, and treating the long-range part by an explicit many-body wave-function method and the short-range part by a density-functional approximation. In particular, for describing systems with van der Waals dispersion interactions, it is appropriate to use methods based on many-body perturbation theory for the long-range part such as second-order perturbation theory AngGerSavTou-PRA-05 (); GerAng-CPL-05b (); GerAng-JCP-07 (); Ang-PRA-08 (); FroJen-PRA-08 (); GolLeiManMitWerSto-PCCP-08 (); JanScu-PCCP-09 (); FroCimJen-PRA-10 (); ChaStoWerLei-MP-10 (); ChaJacAdaStoLei-JCP-10 (); FroJen-JCP-11 (); KulSau-CP-12 (); CorStoJenFro-PRA-13 (), coupled-cluster theory GolWerSto-PCCP-05 (); GolWerStoLeiGorSav-CP-06 (); GolStoThiSch-PRA-07 (); GolWerSto-CP-08 (); GolErnMoeSto-JCP-09 (), or random-phase approximations TouGerJanSavAng-PRL-09 (); JanHenScu-JCP-09 (); JanHenScu-JCP-09b (); JanScu-JCP-09 (); ZhuTouSavAng-JCP-10 (); TouZhuAngSav-PRA-10 (); PaiJanHenScuGruKre-JCP-10 (); TouZhuSavJanAng-JCP-11 (); AngLiuTouJan-JCTC-11 (); IreHenScu-JCP-11 (); GouDob-PRB-11 (); CheMusAngRei-CPL-12 (); MusSzaAng-JCTC-14 ().

Among the advantages of using such many-body methods for the long-range part only of the electron-electron interaction is that they are much less sensitive to the one-electron atomic basis compared to the case of the standard Coulomb interaction. This has been repeatedly observed in calculations using Dunning correlation-consistent basis sets Dun-JCP-89 () for second-order perturbation theory AngGerSavTou-PRA-05 (); GerAng-JCP-07 (); JanScu-PCCP-09 (); KulSau-CP-12 (); CorStoJenFro-PRA-13 (), coupled-cluster theory GolWerSto-PCCP-05 () and random-phase approximations TouGerJanSavAng-PRL-09 (); JanHenScu-JCP-09 (); JanScu-JCP-09 (); TouZhuAngSav-PRA-10 (); IreHenScu-JCP-11 (). The physical reason for this reduced sensitivity to the basis is easy to understand. In the standard Coulomb-interaction case, the many-body wave-function method must describe the short-range part of the correlation hole around the electron-electron coalescence which requires a lot of one-electron basis functions with high angular momentum. In the range-separation case, the many-body method is relieved from describing the short-range part of the correlation hole, which is instead built in the density-functional approximation. The basis set is thus only used to describe a wave function with simply long-range electron-electron correlations (and the one-electron density) which does not require basis functions with very high angular momentum.

In the case of the Coulomb interaction, the rate of convergence of the many-body methods with respect to the size of the basis has been well studied. It has been theoretically shown that, for the ground-state of the helium atom, the partial-wave expansion of the energy calculated by second-order perturbation theory or by full configuration interaction (FCI) converges as where in the maximal angular momentum of the expansion Sch-PR-62 (); Sch-INC-63 (); CarSilMet-JCP-79 (); Hil-JCP-85 (); God-SJMA-09 (). Furthermore, this result has been extended to arbitrary atoms in second-order perturbation theory KutMor-JCP-92 (); GriLud-JPB-02 (). This has motivated the proposal of a scheme for extrapolating the correlation energy to the complete-basis-set (CBS) limit based on a power-law dependence of the correlation energy on the cardinal number of the Dunning hierarchical basis sets HelKloKocNog-JCP-97 (); HALKIER98 (). This extrapolation scheme is widely used, together with other more empirical extrapolation schemes Fel-JCP-92 (); Fel-JCP-93 (); PetKenDun-JCP-93b (); PetDun-JPC-95 (); Mar-CPL-96 (); Tru-CPL-98 (); Fel-JCP-13 (). In the case of range-separated DFT the rate of convergence of the many-body methods with respect to the size of the basis has never been carefully studied, even though the reduced sensitivity to the basis is one of the most appealing feature of this approach.

In this work, we provide a detailed study of the basis convergence of range-separated DFT. First, we review the theory of range-separated DFT methods (Section II) and we study the convergence of the partial-wave expansion of the long-range wave function near the electron-electron coalescence. We show that the rate of convergence is exponential with respect to the maximal angular momentum (Section III). Second, we study the convergence of the long-range second-order Møller-Plesset (MP2) correlation energy of four systems (He, Ne, N, and HO) with the cardinal number of the Dunning basis sets, and find that the error in the correlation energy is best fitted by an exponential in . This leads us to propose a three-point CBS extrapolation scheme for range-separated DFT based on an exponential formula (Section IV).

Hartree atomic units are used throughout this work.

## Ii Range-separated density-functional theory

In range-separated DFT, the exact ground-state energy of an electronic system is expressed as a minimization over multideterminantal wave functions (see, e.g., Ref. TouColSav-PRA-04, )

(1) |

where is the kinetic-energy operator, is the nuclear–electron interaction operator, is the short-range Hartree–exchange–correlation density functional (evaluated at the density of ), and is the long-range electron-electron interaction operator written in terms of the pair-density operator . In this work, we define the long-range interaction with the error function

(2) |

where is the distance between two electrons and (in bohr) controls the range of the separation, with acting as a smooth cutoff radius. For , the long-range interaction vanishes and range-separated DFT reduces to standard Kohn-Sham DFT. In the opposite limit , the long-range interaction becomes the Coulomb interaction and range-separated DFT reduces to standard wave-function theory. In practical applications, one often uses bohr GerAng-CPL-05 (); FroTouJen-JCP-07 ().

The minimizing wave function in Eq. (1) satisfies the Schrödinger-like equation

where is the short-range Hartree–exchange–correlation potential operator (obtained by taking the functional derivative of ), and is the eigenvalue associated with .

In practice, many-body perturbation theory can be used to solve Eq. (LABEL:HPsiEPsi). An appropriate reference for perturbation theory is the range-separated hybrid (RSH) approximation AngGerSavTou-PRA-05 () which is obtained by limiting the search in Eq. (1) to single-determinant wave functions

(4) |

The corresponding minimizing wave function will be denoted by . The exact ground-state energy is then expressed as

(5) |

where is the long-range correlation energy which is to be approximated by perturbation theory. For example, in the long-range variant of MP2 perturbation theory, the long-range correlation energy is AngGerSavTou-PRA-05 ()

(6) |

where is the first-order correction to the wave function (with intermediate normalization). In the basis of RSH spin orbitals , takes a standard MP2 form

where are the long-range two-electron integrals and are the RSH orbital energies. The long-range correlation energy can also be approximated beyond second-order perturbation theory by coupled-cluster GolWerSto-PCCP-05 () or random-phase TouGerJanSavAng-PRL-09 (); JanHenScu-JCP-09 (); PaiJanHenScuGruKre-JCP-10 (); TouZhuSavJanAng-JCP-11 (); AngLiuTouJan-JCTC-11 () approximations. Beyond perturbation theory approaches, Eq. (LABEL:HPsiEPsi) can be (approximately) solved using configuration interaction LeiStoWerSav-CPL-97 (); PolSavLeiSto-JCP-02 (); TouColSav-PRA-04 () or multiconfigurational self-consistent field PedJen-JJJ-XX (); FroTouJen-JCP-07 (); FroReaWahWahJen-JCP-09 () methods. Alternatively, it has also been proposed to use density-matrix functional approximations for the long-range part of the calculation Per-PRA-10 (); RohTouPer-PRA-10 ().

Since the RSH scheme of Eq. (4) simply corresponds to a single-determinant hybrid DFT calculation with long-range Hartree-Fock (HF) exchange, it is clear that the energy has an exponential basis convergence, just as standard HF theory HALKIER99 (). We will thus focus our study on the basis convergence of the long-range wave function and the long-range MP2 correlation energy .

## Iii Partial-wave expansion of the wave function near electron-electron coalescence

In this section, we study the convergence of the partial-wave expansion of the wave function at small interelectronic distances, i.e. near the electron-electron coalescence, which for the case of the Coulomb interaction determines the convergence of the correlation energy. We first briefly review the well-known case of the Coulomb interaction and then consider the case of the long-range interaction.

### iii.1 Coulomb interaction

For systems with Coulomb electron-electron interaction , the electron-electron cusp condition Kat-CPAM-57 () imposes the wave function to be linear with respect to when KUTZELNIGG85 ()

(8) |

Here and in the rest of this section, we consider only the dependence of the wave function on and we restrict ourselves to the most common case of the two electrons being in a natural-parity singlet state KutMor-JCP-92 () for which . The function

(9) |

thus gives the behavior of the wave function at small interelectronic distances. Writing where is the relative angle between the position vectors and of the two electrons, the function can be written as a partial-wave expansion

(10) |

where are the Legendre polynomials and the coefficients are

(11) |

with and . The coefficients decrease slowly with when and are similar. In particular, for , we have as FraMusLupTou-JJJ-XX-note1 (). Therefore, the approximation of by a truncated partial-wave expansion, ,

(12) |

also converges slowly with near . This is illustrated in Figure 1 (left) which shows 1+ as a function of for bohr for increasing values of the maximal angular momentum . Comparing with the converged value corresponding to [Eq. (10)], it is clear that the convergence near the singularity at is indeed painstakingly slow.

This slow convergence of the wave function near the electron-electron coalescence leads to the slow power-law convergence of the partial-wave increments to the correlation energy Sch-PR-62 (); Sch-INC-63 (); CarSilMet-JCP-79 (); KutMor-JCP-92 (); GriLud-JPB-02 () or, equivalently, to the power-law convergence of the truncation error in the correlation energy Hil-JCP-85 (); God-SJMA-09 ().

### iii.2 Long-range interaction

For systems with the long-range electron-electron interaction , the behavior of the wave function for small interelectronic distances was determined by Gori-Giorgi and Savin GorSav-PRA-06 ()

(13) |

where the function is given by

(14) |

We thus need to study the function

(15) |

For a fixed value of , and for , it yields

(16) |

which exhibits no linear term in , i.e. no electron-electron cusp. On the other hand, for and , we obtain

(17) |

i.e. the Coulomb electron-electron cusp is recovered. The function thus makes the transition between the cuspless long-range wave function and the Coulomb wave function.

As for the Coulomb case, we write as a partial-wave expansion

(18) |

and calculate with Mathematica Math9-PROG-12 () the coefficients for each

(19) |

with , , and using the following explicit expression for

(20) |

Since the partial-wave expansion of the first term in in Eq. (16) terminates at , we expect a fast convergence with of , for small enough, and thus also a fast convergence of the truncated partial-wave expansion

(21) |

Plots of this truncated partial-wave expansion for in Figure 1 (right) confirm this expectation. The Coulomb singularity at has disappeared and the approximation converges indeed very fast with , being converged to better than a.u. already at .

We study in detail the dependence of the coefficients of the partial-wave expansion on . We compare two possible convergence behaviors, a power-law form

(22) |

and an exponential-law form

(23) |

where , , and are (-dependent) parameters.

To determine which form best represents , in Figure 2 we plot for as a function of (left) and as a function of (right), for several values of , as well as for the Coulomb case () FraMusLupTou-JJJ-XX-note4 (). A straight line on the plot of vs. indicates a power-law dependence, whereas a straight line on the plot of vs. indicates an exponential-law dependence.

For the Coulomb case (black curve nearly superimposed with the curve for ), we observe that the plot of vs. is linear, whereas the plot of vs. is curved upward. This is expected for a power law form. Moreover, we find as expected from Section III.1. When going from large to small values of , we observe that the plot of vs. becomes more and more curved downward, and the plot of vs. becomes more and more linear. We thus go from a power-law dependence to an exponential-law dependence. Already for , the exponential law is a better description than the power law.

When decreases, the absolute value of the slope of the plot of vs. increases, i.e. the convergence becomes increasingly fast. More precisely, we have found for .

The exponential convergence of the partial-wave expansion of the long-range wave function near the electron-electron coalescence implies a similar exponential convergence for the partial-wave expansion of the corresponding energy. The present study is thus consistent with the approximate exponential convergence of the partial-wave expansion of the energy of the helium atom in the presence of a long-range electron-electron interaction reported in Refs. SIRBU02, ; SIRBU03, . However, no quantitative comparison can be made between the latter work and the present work since the form of the long-range interaction is different.

## Iv Convergence in one-electron atomic basis sets

In this section, we study the convergence of the long-range wave function and correlation energy with respect to the size of the one-particle atomic basis. This problem is closely related to the convergence of the partial-wave expansion studied in the previous section. Indeed, for a two-electron atom in a singlet S state, it is possible to use the spherical-harmonic addition theorem to obtain the partial-wave expansion in terms of the relative angle between two electrons by products of the spherical harmonic part of the one-particle atomic basis functions

(24) |

where with spherical angles , and ,. The partial-wave expansion can thus be obtained from a one-particle atomic basis, provided that the basis saturates the radial degree of freedom for each angular momentum . In practice, of course, for the basis sets that we use, this latter condition is not satisfied. Nevertheless, one can expect the convergence with the maximal angular momentum of the basis to be similar to the convergence of the partial-wave expansion.

For this study, we have analyzed the behavior of He, Ne, N, and HO at the same experimental geometries used in Ref. HALKIER98, ( Å, Å and ). We performed all the calculations with the program MOLPRO 2012 Molproshort-PROG-12 () using Dunning correlation-consistent cc-p(C)VZ basis sets for which we studied the convergence with respect to the cardinal number , corresponding to a maximal angular momentum of for He and for atoms from Li to Ne. We emphasize that the series of Dunning basis sets does not correspond to a partial-wave expansion but to a principal expansion KloBakJorOlsHel-JPB-99 (); HelJorOls-BOOK-02 () with maximal quantum number for He and for Li to Ne. The short-range exchange-correlation PBE density functional of Ref. GolWerStoLeiGorSav-CP-06, (which corresponds to a slight modification of the one of Ref. TouColSav-JCP-05, ) was used in all range-separated calculations.

### iv.1 Convergence of the wave function

We start by analyzing the convergence of the FCI ground-state wave function of the He atom with respect to the cardinal number of the cc-pVZ basis sets. We perform a FCI calculation with the long-range Hamiltonian in Eq. (LABEL:HPsiEPsi) using a fixed RSH density, calculated from the orbitals obtained in Eq. (4), in the short-range Hartree–exchange-correlation potential. To facilitate the extraction of the wave function from the program, we use the Löwdin-Shull diagonal representation of the spatial part of the FCI wave function in terms of the spatial natural orbitals (NO) LowShu-PR-56 (); SzaOst-BOOK-96 ()

(25) |

where the coefficients are related to the NO occupation numbers by the relation . As the signs of are undetermined we have chosen a positive leading coefficient , and we assumed that all the other coefficients are negative for GoeUmr-INC-00 (). Even though it has been shown that, for the case of the Coulomb interaction, there are in fact positive coefficients in the expansion in addition to the leading one, for a weakly correlated system such as the He atom, these positive coefficients appear only in larger basis sets than the ones that we consider here and have negligible magnitude SheMenGriBae-JCP-13 (); GieLee-JCP-13 (); GieLee-JCP-13b ().

In Figure 3 we show the convergence of the FCI wave function with the cardinal number for which corresponds to the Coulomb interaction (left) and for (right). The first electron is fixed at the Cartesian coordinates (measured from the nucleus) and the position of the second electron is varied on a circle at the same distance of the nucleus, . For the Coulomb interaction, we compare with the essentially exact curve obtained with a highly accurate 418-term Hylleraas-type wave function FreHuxMor-PRA-84 (); BakFreHilMor-PRA-90 (); UmrGon-PRA-94 (). The curve of as a function of reveals the angular correlation between the electrons FraMusLupTou-JJJ-XX-note2 (). Clearly, correlation is much weaker for the long-range interaction. Note that a single-determinant wave function , where is a spherically symmetric 1s orbital, does not depend on , and the HF and RSH single-determinant wave functions indeed just give horizontal lines in Figure 3.

The fact that Figure 3 resembles Figure 1 confirms that the convergence with respect to is similar to the convergence of the partial-wave expansion with respect to , and thus corroborates the relevance of the study of Section III for practical calculations. As for the partial-wave expansion, the convergence with of the Coulomb wave function near the electron-electron cusp is exceedingly slow. The long-range wave function does not have an electron-electron cusp and converges much faster with , the differences between the curves obtained with the cc-pV5Z and cc-pV6Z basis being smaller than 0.03 mhartree. Note, however, that the convergence of the long-range wave function seems a bit less systematic than the convergence of the Coulomb wave function, with the difference between the cc-pVQZ and cc-pV5Z basis being about 3 times larger than the difference between the cc-pVTZ and cc-pVQZ basis. This may hint to the fact the Dunning basis sets have been optimized for the Coulomb interaction and are not optimal for the long-range interaction. Finally, we note that we have found the same convergence behavior with the short-range exchange-correlation LDA density functional of Ref. PazMorGorBac-PRB-06, .

### iv.2 Convergence of the correlation energy

Coulomb interaction | ||||||||

He | Ne | N | HO | |||||

Basis | ||||||||

cc-pVDZ | -25.828 | 11.549 | -185.523 | 134.577 | -306.297 | 114.903 | -201.621 | 98.479 |

cc-pVTZ | -33.138 | 4.239 | -264.323 | 55.777 | -373.683 | 47.517 | -261.462 | 38.638 |

cc-pVQZ | -35.478 | 1.899 | -293.573 | 26.527 | -398.749 | 22.451 | -282.798 | 17.302 |

cc-pV5Z | -36.407 | 0.970 | -306.166 | 13.934 | -409.115 | 12.085 | -291.507 | 8.593 |

cc-pV6Z | -36.807 | 0.570 | -311.790 | 8.310 | -413.823 | 7.377 | -295.205 | 4.895 |

CBS limit | -37.377 | -320(1) | -421(2) | -300(1) | ||||

Long-range interaction | ||||||||

He | Ne | N | HO | |||||

Basis | ||||||||

cc-pVDZ | -0.131 | 0.227 | -0.692 | 1.963 | -20.178 | 3.316 | -6.462 | 3.532 |

cc-pVTZ | -0.262 | 0.096 | -1.776 | 0.879 | -22.663 | 0.830 | -8.956 | 1.038 |

cc-pVQZ | -0.322 | 0.036 | -2.327 | 0.328 | -23.263 | 0.231 | -9.626 | 0.367 |

cc-pV5Z | -0.346 | 0.012 | -2.557 | 0.098 | -23.430 | 0.064 | -9.901 | 0.092 |

cc-pV6Z | -0.358 | -2.655 | -23.494 | -9.993 | ||||

Taken from Ref. PatCenJezJezSza-JPCA-07, where it was obtained by a Gaussian-type geminal MP2 calculation. | ||||||||

Taken from Ref. HALKIER98, where it was estimated from R12-MP2 calculations. |

We also study the basis convergence of the long-range MP2 correlation energy, given in Eq. (LABEL:MP2orb), calculated with RSH orbitals for He, Ne, N and HO.

In Table 1 we show the valence MP2 correlation energies and their errors as a function of the cardinal number of the cc-pVZ basis sets for X . We compare the long-range MP2 correlation energies at and the standard Coulomb MP2 correlation energies corresponding to . For the case of the Coulomb interaction, the error is calculated as where is the MP2 correlation energy in the estimated CBS limit taken from Refs. PatCenJezJezSza-JPCA-07, ; HALKIER98, . For the range-separated case we do not have an independent estimate of the CBS limit of the long-range MP2 correlation energy for a given value of . Observing that the difference between the long-range MP2 correlation energies for and is below 0.1 mhartree for , we choose the cc-pV6Z result as a good estimate of the CBS limit. Of course, the accuracy of this CBS estimate will deteriorate for larger values of , but in practice this is a good estimate for the range of values of in which we are interested, i.e. FraMusLupTou-JJJ-XX-note3 (). The error on the long-range correlation energy is thus calculated as .

Coulomb interaction | ||||||||
---|---|---|---|---|---|---|---|---|

Power law | Exponential law | |||||||

He | 2 | 6 | 2.749 | 81.84 | 99.82 | 0.749 | 44.04 | 98.55 |

3 | 6 | 2.902 | 104.00 | 99.97 | 0.669 | 29.51 | 99.19 | |

Ne | 2 | 6 | 2.543 | 843.36 | 99.59 | 0.696 | 479.76 | 99.00 |

3 | 6 | 2.754 | 1169.90 | 99.93 | 0.636 | 355.28 | 99.37 | |

N | 2 | 6 | 2.513 | 697.74 | 99.71 | 0.686 | 397.67 | 98.76 |

3 | 6 | 2.693 | 923.76 | 99.98 | 0.621 | 286.90 | 99.16 | |

HO | 2 | 6 | 2.742 | 717.04 | 99.52 | 0.751 | 391.27 | 99.09 |

3 | 6 | 2.988 | 1051.28 | 99.92 | 0.690 | 288.61 | 99.39 | |

Long-range interaction | ||||||||

Power law | Exponential law | |||||||

He | 2 | 5 | 3.128 | 2.35 | 96.71 | 0.974 | 1.68 | 99.77 |

3 | 5 | 3.997 | 8.16 | 99.11 | 1.028 | 2.13 | 99.95 | |

Ne | 2 | 5 | 3.189 | 22.06 | 95.04 | 0.998 | 15.97 | 99.18 |

3 | 5 | 4.257 | 101.51 | 98.27 | 1.098 | 24.57 | 99.65 | |

N | 2 | 5 | 4.257 | 73.26 | 98.63 | 1.313 | 44.49 | 99.97 |

3 | 5 | 4.997 | 211.13 | 99.44 | 1.283 | 39.08 | 99.99 | |

HO | 2 | 5 | 3.861 | 60.27 | 97.30 | 1.198 | 39.24 | 99.72 |

3 | 5 | 4.686 | 196.20 | 97.63 | 1.210 | 41.50 | 99.33 |

The first observation to be made is that the long-range MP2 correlation energies only represent about 1 to 5 % of the Coulomb MP2 correlation energies. Although the long-range correlation energy may appear small, it is nevertheless essential for the description of dispersion interactions for instance. The errors on the long-range MP2 correlation energies are also about two orders of magnitude smaller than the errors on the Coulomb MP2 correlation energies.

Inspired by Ref. HALKIER99, , we compare two possible forms of convergence for the correlation energy: a power-law form

(26) |

and an exponential-law form

(27) |

where is the CBS limit of the long-range correlation energy and , , , and are parameters depending on , as in Section III. In practice, we actually make linear fits of the logarithm of the error for the two forms:

(28) |

(29) |

In Table 2, we show the results of the fits for the Coulomb interaction and the long-range interaction at using either the complete range of or excluding the value for . We use the squared Pearson correlation coefficient as a measure of the quality of the fit. For the Coulomb interaction and for all the systems studied, the best fit is achieved with the power law with , which is roughly what is expected HALKIER98 (). We note however that the fits to the exponential law are also very good with when the value is excluded. This explains why extrapolations of the total energy based on an exponential formula have also been used for the case of the Coulomb interaction Fel-JCP-92 (); PetDun-JPC-95 (). For the case of the long-range interaction, the difference between the fits to the power law and to the exponential law is much bigger. The best fit is by far obtained for the exponential law with for all systems and ranges of considered. This exponential convergence of the long-range correlation energy with respect to is in accordance with the exponential convergence of the partial-wave expansion of the long-range wave-function observed in Section III.

We have also performed fits for several other values of between and and always obtained an exponential convergence of the long-range valence MP2 correlation energy with respect to . However, contrary to what was observed for the partial-wave expansion, we found that when decreases also decreases a bit for the four systems considered. In other words, when the interaction becomes more long range, the convergence of the long-range correlation energy becomes slower. This surprising result may be due to the fact that the cc-pV6Z result may not be as good an estimate of the CBS limit when increases. When decreases, the prefactor decreases and goes to zero for , as expected. Moreover, we have also checked that we obtain very similar results for the long-range all-electron MP2 correlation energy (including core excitations) with cc-pCVZ basis sets.

We note that Prendergast et al. PreNolFilFahGre-JCP-01 () have argued that the removal of the electron-electron cusp in a small region around the coalescence point does not significantly improve the convergence of the energy in the millihartree level of accuracy. At first sight, their conclusion might appear to be in contradiction with our observation of the exponential convergence of the long-range correlation energy with . There are however important differences between the two studies: (1) their form of long-range interaction is different from ours, (2) they consider interelectronic distance “cutoffs” of bohr whereas we consider larger “cutoffs” bohr, (3) they do not investigate exponential-law versus power-law convergence.

Finally, in the Appendix we provide a complement analysis of the basis convergence of the correlation energy of the He atom for truncated configuration-interaction (CI) calculations in natural orbitals. The analysis shows that, contrary to the case of the Coulomb interaction, the convergence of the long-range correlation energy is no longer limited by the truncation rank the CI wave function but by the basis convergence of the natural orbitals themselves. This result is consistent with an exponential basis convergence of the long-range correlation energy.

### iv.3 Extrapolation scheme

For the long-range interaction case, since both the RSH energy and the long-range correlation energy have an exponential convergence with respect to the cardinal number , we propose to extrapolate the total energy to the CBS limit by using a three-point extrapolation scheme based on an exponential formula. Suppose that we have calculated three total energies , , for three consecutive cardinal numbers , , . If we write

(30) |

(31) |

(32) |

and eliminate the unknown parameters and , we obtain the following estimate of the CBS-limit total energy

(33) |

In Table 3, we report the errors on the RSH+lrMP2 total energy, , obtained with the three-point extrapolation formula using either , , () or , , (), and we compare with the errors obtained with each cc-pVZ basis set from to (). Here again the errors are calculated with respect to the cc-pV6Z total energy, for several values of the range-separation parameter , , , and only valence excitations are included in the MP2 calculations. For all the systems studied the errors are less than 1.5 mhartree. For Ne, N, and HO, these errors are significantly smaller (by a factor of about 3 to 15) than the errors obtained with the largest basis used for the extrapolation, and are overall comparable with the errors . Thus, the three-point extrapolation formula with , , provides a useful CBS extrapolation scheme for range-separated DFT. Except for He, the errors are negative (i.e. , the extrapolation overshoots the CBS limit) and larger than the errors . Thus, the three-point extrapolation scheme with , , does not seem to be useful.

These conclusions extend to calculations including core excitations with cc-pCVZ basis sets, which are presented in Table 4. All the errors are smaller than for the valence-only calculations. The errors are now less than 0.9 mhartree, and are smaller or comparable to the errors . The errors are always negative and are overall larger than the errors .

He | 0.1 | 8.508 | 0.772 | 0.261 | 0.089 | 0.224 | 0.003 |
---|---|---|---|---|---|---|---|

0.5 | 8.488 | 0.781 | 0.245 | 0.078 | 0.205 | 0.002 | |

1.0 | 8.258 | 0.924 | 0.259 | 0.078 | 0.192 | 0.011 | |

Ne | 0.1 | 72.999 | 20.215 | 5.842 | 0.716 | 0.464 | -2.127 |

0.5 | 74.523 | 20.337 | 5.763 | 0.751 | 0.401 | -1.876 | |

1.0 | 79.311 | 20.962 | 5.726 | 0.803 | 0.342 | -1.548 | |

N | 0.1 | 47.061 | 13.026 | 4.136 | 0.853 | 0.993 | -1.069 |

0.5 | 51.581 | 13.406 | 4.090 | 0.810 | 1.083 | -0.972 | |

1.0 | 61.053 | 15.108 | 4.513 | 0.868 | 1.337 | -1.043 | |

HO | 0.1 | 54.861 | 15.229 | 5.005 | 0.857 | 1.451 | -1.975 |

0.5 | 55.850 | 14.736 | 4.499 | 0.726 | 1.105 | -1.475 | |

1.0 | 61.013 | 15.212 | 4.423 | 0.724 | 1.099 | -1.206 |

Ne | 0.1 | 70.932 | 18.941 | 4.929 | 0.522 | -0.240 | -1.501 |
---|---|---|---|---|---|---|---|

0.5 | 72.501 | 18.990 | 4.831 | 0.537 | -0.263 | -1.333 | |

1.0 | 77.497 | 19.517 | 4.775 | 0.554 | -0.250 | -1.140 | |

N | 0.1 | 43.528 | 10.237 | 2.334 | 0.459 | -0.126 | -0.125 |

0.5 | 48.079 | 10.451 | 2.285 | 0.413 | 0.021 | -0.144 | |

1.0 | 57.942 | 12.118 | 2.677 | 0.467 | 0.227 | -0.209 | |

HO | 0.1 | 52.875 | 13.897 | 4.132 | 0.680 | 0.868 | -1.208 |

0.5 | 53.936 | 13.350 | 3.614 | 0.534 | 0.541 | -0.891 | |

1.0 | 59.290 | 13.789 | 3.527 | 0.521 | 0.539 | -0.724 |

We have also tested a more flexible extrapolation scheme where the RSH energy and the long-range MP2 correlation energy are exponentially extrapolated independently but we have not found significant differences. On the contrary, one may want to use a less flexible two-point extrapolation formula using a predetermined value for . The difficulty with such an approach is to choose the value of , which in principle should depend on the system, on the range-separated parameter , and on the long-range wave-function method used. For this reason, we do not consider two-point extrapolation schemes.

## V Conclusions

We have studied in detail the basis convergence of range-separated DFT. We have shown that the partial-wave expansion of the long-range wave function near the electron-electron coalescence converges exponentially with the maximal angular momentum . We have also demonstrated on four systems (He, Ne, N, and HO) that the long-range MP2 correlation energy converges exponentially with the cardinal number of the Dunning basis sets cc-p(C)VZ. This contrasts with the slow convergence of the correlation energy for the standard case of the Coulomb interaction. Due to this exponential convergence, the extrapolation to the CBS limit is less necessary for range-separated DFT than for standard correlated wave function methods. Nevertheless, we have proposed a CBS extrapolation scheme for the total energy in range-separated DFT based on an exponential formula using calculations from three cardinal numbers . For the systems studied, the extrapolation using gives an error on the total energy with respect to the estimated CBS limit which is always smaller than the error obtained with a single calculation at , and which is often comparable or smaller than the error obtained with a calculation at .

We expect the same convergence behavior for range-separated DFT methods in which the long-range part is treated by configuration interaction, coupled-cluster theory, random-phase approximations, or density-matrix functional theory. Finally, it should be pointed out that this rapid convergence is obtained in spite of the fact that the Dunning basis sets have been optimized for the case of the standard Coulomb interaction. The construction of basis sets specially optimized for the case of the long-range interaction may give yet a faster and more systematic convergence.

Coulomb interaction | |||||

rank | VDZ | VTZ | VQZ | V5Z | V6Z |

1s2s | -14.997 | -15.806 | -16.087 | -16.204 | -16.212 |

1s2s2p | -32.434 | -35.256 | -35.664 | -35.794 | -35.808 |

1s2s2p3s | -35.909 | -36.425 | -36.595 | -36.618 | |

1s2s2p3s3p | -37.448 | -38.068 | -38.261 | -38.291 | |

1s2s2p3s3p3d | -39.079 | -39.807 | -39.999 | -40.028 | |

Long-range interaction | |||||

rank | VDZ | VTZ | VQZ | V5Z | V6Z |

1s2s | -0.018 | -0.036 | -0.038 | -0.039 | -0.039 |

1s2s2p | -0.155 | -0.329 | -0.415 | -0.449 | -0.468 |

1s2s2p3s | -0.329 | -0.415 | -0.450 | -0.468 | |

1s2s2p3s3p | -0.329 | -0.415 | -0.450 | -0.469 | |

1s2s2p3s3p3d | -0.329 | -0.416 | -0.451 | -0.470 |

Coulomb interaction | ||||||
---|---|---|---|---|---|---|

Power law | Exponential law | |||||

rank | ||||||

1s2s | 5.081 | 65.641 | 87.44 | 1.618 | 43.148 | 94.39 |

1s2s2p | 5.765 | 241.923 | 95.45 | 1.797 | 131.516 | 98.80 |

1s2s2p3s | 6.616 | 1219.794 | 95.63 | 1.716 | 140.022 | 98.11 |

1s2s2p3s3p | 6.440 | 1168.927 | 96.46 | 1.668 | 140.605 | 98.65 |

1s2s2p3s3p3d | 6.771 | 1874.080 | 97.16 | 1.751 | 200.085 | 99.07 |

Long-range interaction | ||||||

Power law | Exponential law | |||||

rank | ||||||

1s2s | 3.974 | 0.321 | 99.45 | 1.206 | 0.188 | 97.51 |

1s2s2p | 3.058 | 3.108 | 96.28 | 0.953 | 2.249 | 99.64 |

1s2s2p3s | 3.956 | 11.283 | 99.01 | 1.018 | 2.991 | 99.93 |

1s2s2p3s3p | 3.927 | 10.979 | 99.02 | 1.010 | 2.938 | 99.93 |

1s2s2p3s3p3d | 3.892 | 10.701 | 99.00 | 1.001 | 2.898 | 99.93 |

## Acknowledgements

We thank Andreas Savin, Trygve Helgaker, Éric Cancès, and Gabriel Stoltz for insightful discussions. This work was supported by French state funds managed by CALSIMLAB and the ANR within the Investissements d’Avenir program under reference ANR-11-IDEX-0004-02.

## Appendix A Convergence of the correlation energy for truncated CI calculations

In this Appendix, we explore the basis convergence of the correlation energy of the He atom for truncated CI calculations for both the Coulomb and long-range interactions. For a given basis set and interaction, we start by performing a FCI calculation and generating the corresponding natural orbitals. We then use these natural orbitals in truncated CI calculations for increasing orbital active spaces 1s2s, 1s2s2p, 1s2s2p3s, 1s2s2p3s3p, and 1s2s2p3s3p3d. Table 5 shows the Coulomb and long-range correlation energies for the different basis sets and truncation ranks. For the Coulomb interaction, the correlation energy for a fixed rank converges rapidly with the basis size, while the convergence with respect to the rank is much slower. For the long-range interaction, the correlation energy jumps by one order of magnitude when including the 2p natural orbital, which is consistent with the fact that the long-range interaction brings in first angular correlation effects PolColLeiStoWerSav-IJQC-03 (). The long-range correlation energy is essentially converged at rank 1s2s2p, and the overall convergence is now determined by the basis convergence of the natural orbitals.

Finally, we compare two possible forms for the convergence of the truncated CI correlation energies with the cardinal number , the power law Eq. (26) and the exponential law Eq. (27). Using as reference the results obtained with the cc-pV6Z basis set, we have calculated for the different truncation ranks the correlation energy errors for the Coulomb interaction, , and for the long-range interaction, , and performed logarithmic fits as in Section IV.2. Table 6 shows the results of the fits. For both the Coulomb and long-range interactions, for the rank 1s2s2p and larger, the best fit is achieved with the exponential law . Thus, in comparison with the Coulomb interaction, the long-range interaction does not significantly change the basis convergence of the correlation energy at a fixed truncation rank. However, for the long-range interaction, this exponential convergence at a fixed truncation rank becomes the dominant limitation to the overall basis convergence.

## References

- (1) J. Toulouse, F. Colonna, and A. Savin, Phys. Rev. A 70, 062505 (2004).
- (2) P. Hohenberg and W. Kohn, Phys. Rev. 136, B 864 (1964).
- (3) W. Kohn and L. J. Sham, Phys. Rev. 140, A1133 (1965).
- (4) J. G. Ángyán, I. C. Gerber, A. Savin, and J. Toulouse, Phys. Rev. A 72, 012510 (2005).
- (5) I. C. Gerber and J. G. Ángyán, Chem. Phys. Lett. 416, 370 (2005).
- (6) I. C. Gerber and J. G. Ángyán, J. Chem. Phys. 126, 044103 (2007).
- (7) J. G. Ángyán, Phys. Rev. A 78, 022510 (2008).
- (8) E. Fromager and H. J. A. Jensen, Phys. Rev. A 78, 022504 (2008).
- (9) E. Goll, T. Leininger, F. R. Manby, A. Mitrushchenkov, H.-J. Werner, and H. Stoll, Phys. Chem. Chem. Phys. 10, 3353 (2008).
- (10) B. G. Janesko and G. E. Scuseria, Phys. Chem. Chem. Phys. 11, 9677 (2009).
- (11) E. Fromager, R. Cimiraglia, and H. J. A. Jensen, Phys. Rev. A 81, 024502 (2010).
- (12) S. Chabbal, H. Stoll, H.-J. Werner, and T. Leininger, Mol. Phys. 108, 3373 (2010).
- (13) S. Chabbal, D. Jacquemin, C. Adamo, H. Stoll, and T. Leininger, J. Chem. Phys. 133, 151104 (2010).
- (14) E. Fromager and H. J. A. Jensen, J. Chem. Phys. 135, 034116 (2011).
- (15) O. Kullie and T. Saue, Chem. Phys. 395, 54 (2012).
- (16) Y. Cornaton, A. Stoyanova, H. J. A. Jensen, and E. Fromager, Phys. Rev. A 88, 022516 (2013).
- (17) E. Goll, H.-J. Werner, and H. Stoll, Phys. Chem. Chem. Phys. 7, 3917 (2005).
- (18) E. Goll, H.-J. Werner, H. Stoll, T. Leininger, P. Gori-Giorgi, and A. Savin, Chem. Phys. 329, 276 (2006).
- (19) E. Goll, H. Stoll, C. Thierfelder, and P. Schwerdtfeger, Phys. Rev. A 76, 032507 (2007).
- (20) E. Goll, H.-J. Werner, and H. Stoll, Chem. Phys. 346, 257 (2008).
- (21) E. Goll, M. Ernst, F. Moegle-Hofacker, and H. Stoll, J. Chem. Phys. 130, 234112 (2009).
- (22) J. Toulouse, I. C. Gerber, G. Jansen, A. Savin, and J. G. Ángyán, Phys. Rev. Lett. 102, 096404 (2009).
- (23) B. G. Janesko, T. M. Henderson, and G. E. Scuseria, J. Chem. Phys. 130, 081105 (2009).
- (24) B. G. Janesko, T. M. Henderson, and G. E. Scuseria, J. Chem. Phys. 131, 034110 (2009).
- (25) B. G. Janesko and G. E. Scuseria, J. Chem. Phys. 131, 154106 (2009).
- (26) W. Zhu, J. Toulouse, A. Savin, and J. G. Ángyán, J. Chem. Phys. 132, 244108 (2010).
- (27) J. Toulouse, W. Zhu, J. G. Ángyán, and A. Savin, Phys. Rev. A 82, 032502 (2010).
- (28) J. Paier, B. G. Janesko, T. M. Henderson, G. E. Scuseria, A. Grüneis, and G. Kresse, J. Chem. Phys. 132, 094103 (2010).
- (29) J. Toulouse, W. Zhu, A. Savin, G. Jansen, and J. G. Ángyán, J. Chem. Phys. 135, 084119 (2011).
- (30) J. G. Ángyán, R.-F. Liu, J. Toulouse, and G. Jansen, J. Chem. Theory Comput. 7, 3116 (2011).
- (31) R. M. Irelan, T. M. Henderson, and G. E. Scuseria, J. Chem. Phys. 135, 094105 (2011).
- (32) T. Gould and J. F. Dobson, Phys. Rev. B 84, 241108 (2011).
- (33) E. Chermak, B. Mussard, J. G. Ángyán, and P. Reinhardt, Chem. Phys. Lett. 550, 162 (2012).
- (34) B. Mussard, P. G. Szalay, and J. G. Ángyán, J. Chem. Theory Comput. 10, 1968 (2014).
- (35) T. H. Dunning, J. Chem. Phys. 90, 1007 (1989).
- (36) C. Schwartz, Phys. Rev. 126, 1015 (1962).
- (37) C. Schwartz, in Methods in Computational Physics, Vol. 2, edited by B. Alder, S. Fernbach, and M. Rotenberg, pages 241–266, Academic Press, New York and London, 1963.
- (38) D. P. Carroll, H. J. Silverstone, and R. M. Metzger, J. Chem. Phys. 71, 4142 (1979).
- (39) R. N. Hill, J. Chem. Phys. 83, 1173 (1985).
- (40) B. D. Goddard, Siam J. Math. Anal. 41, 77 (2009).
- (41) W. Kutzelnigg and J. D. Morgan III, J. Chem. Phys. 96, 4484 (1992).
- (42) G. F. Gribakin and J. Ludlow, J. Phys. B 35, 339 (2002).
- (43) T. Helgaker, W. Klopper, H. Koch, and J. Noga, J. Chem. Phys. 106, 9639 (1997).
- (44) A. Halkier, T. Helgaker, P. Jørgensen, W. Klopper, H. Koch, J. Olsen, and A. K. Wilson, Chem. Phys. Lett. 286, 243 (1998).
- (45) D. Feller, J. Chem. Phys. 96, 6104 (1992).
- (46) D. Feller, J. Chem. Phys. 98, 7059 (1993).
- (47) K. A. Peterson, R. A. Kendall and T. H. Dunning Jr., J. Chem. Phys. 99, 9790 (1993).
- (48) K. A. Peterson and T. H. Dunning Jr., J. Phys. Chem. 99, 3898 (1995).
- (49) J. M. L. Martin, Chem. Phys. Lett. 259, 669 (1996).
- (50) D. G. Truhlar, Chem. Phys. Lett. 294, 45 (1998).
- (51) D. Feller, J. Chem. Phys. 138, 074103 (2013).
- (52) I. Gerber and J. G. Ángyán, Chem. Phys. Lett. 415, 100 (2005).
- (53) E. Fromager, J. Toulouse, and H. J. A. Jensen, J. Chem. Phys. 126, 074111 (2007).
- (54) T. Leininger, H. Stoll, H.-J. Werner, and A. Savin, Chem. Phys. Lett. 275, 151 (1997).
- (55) R. Pollet, A. Savin, T. Leininger, and H. Stoll, J. Chem. Phys. 116, 1250 (2002).
- (56) J. K. Pedersen and H. J. A. Jensen, A second order MCSCF-DFT hybrid algorithm (unpublished).
- (57) E. Fromager, F. Réal, P. Wåhlin, U. Wahlgren, and H. J. A. Jensen, J. Chem. Phys. 131, 054107 (2009).
- (58) K. Pernal, Phys. Rev. A 81, 052511 (2010).
- (59) D. R. Rohr, J. Toulouse, and K. Pernal, Phys. Rev. A 82, 052502 (2010).
- (60) A. Halkier, T. Helgaker, P. Jørgensen, W. Klopper, and J. Olsen, Chem. Phys. Lett. 302, 437 (1999).
- (61) T. Kato, Comm. Pure Appl. Math. 10, 151 (1957).
- (62) W. Kutzelnigg, Theor. Chem. Acta 68, 445 (1985).
- (63) Note that the partial-wave expansion of any odd power of contains an infinite number of terms. In contrast, the partial-wave expansion of any even power of terminates at a finite .
- (64) P. Gori-Giorgi and A. Savin, Phys. Rev. A 73, 032506 (2006).
- (65) Wolfram Research, Inc., Mathematica, Version 9.0, Champaign, IL (2012).
- (66) For , for and for .
- (67) I. Sirbu and H. F. King, J. Chem. Phys. 117, 6411 (2002).
- (68) I. Sirbu and H. F. King, Int. J. Quantum Chem. 92, 433 (2003).
- (69) H.-J. Werner, P. J. Knowles, G. Knizia, F. R. Manby, M. Schütz, and others, Molpro, version 2012.1, a package of ab initio programs, Cardiff, UK, 2012, see http://www.molpro.net.
- (70) W. Klopper, K. L. Bak, P. Jørgensen, J. Olsen, and T. Helgaker, J. Phys. B 32, R103 (1999).
- (71) T. Helgaker, P. Jørgensen, and J. Olsen, Molecular Electronic-Structure Theory, Wiley, Chichester, 2002.
- (72) J. Toulouse, F. Colonna, and A. Savin, J. Chem. Phys. 122, 014110 (2005).
- (73) D. E. Freund, B. D. Huxtable, and J. D. Morgan III, Phys. Rev. A 29, 980 (1984).
- (74) J. D. Baker, D. E. Freund, R. N. Hill, and J. D. Morgan III, Phys. Rev. A 41, 1247 (1990).
- (75) C. J. Umrigar and X. Gonze, Phys. Rev. A 50, 3827 (1994).
- (76) P.-O. Löwdin and H. Shull, Phys. Rev. 101, 1730 (1956).
- (77) A. Szabo and N. S. Ostlund, Modern Quantum Chemistry: Introduction to Advanced Electronic Structure Theory, Dover, New York, 1996.
- (78) S. Goedecker and C. J. Umrigar, in Many-Electron Densities and Reduced Density Matrices, edited by J. Cioslowski, pages 165–181, Kluwer Academic, Dordrecht/New York, 2000.
- (79) X. W. Sheng, Ł. M. Mentel, O. V. Gritsenko, and E. J. Baerends, J. Chem. Phys. 138, 164105 (2013).
- (80) K. J. H. Giesbertz and R. van Leeuwen, J. Chem. Phys. 139, 104109 (2013).
- (81) K. J. H. Giesbertz and R. van Leeuwen, J. Chem. Phys. 139, 104110 (2013).
- (82) For a discussion of angular and radial correlation and the shape of the correlation hole in the He atom with the long-range interaction, see Ref. PolColLeiStoWerSav-IJQC-03, .
- (83) S. Paziani, S. Moroni, P. Gori-Giorgi, and G. B. Bachelet, Phys. Rev. B 73, 155111 (2006).
- (84) K. Patkowski, W. Cencek, M. Jeziorska, B. Jeziorski, and K. Szalewicz, J. Phys. Chem. A 111, 7611 (2007).
- (85) For the difference between the long-range MP2 correlation energies for and is still as small as about 0.2 mhartree.
- (86) D. Prendergast, M. Nolan, C. Filippi, S. Fahy, and J. C. Greer, J. Chem. Phys. 115, 1626 (2001).
- (87) R. Pollet, F. Colonna, T. Leininger, H. Stoll, H.-J. Werner, and A. Savin, Int. J. Quantum. Chem. 91, 84 (2003).