Using JCP format - Washington State University

4 downloads 0 Views 500KB Size Report
Nov 8, 2000 - calculations, use of the MP3 method to simulate the effect of ... should be significant since the MP3 method formally scales as N6, while the ...
JOURNAL OF CHEMICAL PHYSICS

VOLUME 113, NUMBER 18

8 NOVEMBER 2000

Approximating the basis set dependence of coupled cluster calculations: Evaluation of perturbation theory approximations for stable molecules Thom H. Dunning, Jr.a) and Kirk A. Petersonb)

Environmental Molecular Sciences Laboratory, Pacific Northwest National Laboratory,c兲 Richland, Washington 99352

共Received 4 February 2000; accepted 17 August 2000兲 The coupled cluster CCSD共T兲 method provides a theoretically sound, accurate description of the electronic structure of a wide range of molecules. To obtain accurate results, however, very large basis sets must be used. Since the computational cost of CCSD共T兲 calculations formally increases with the seventh power of the number of basis functions (N 7 ), the CCSD共T兲 method can only be applied to a restricted range of molecules. In this work we show that the basis set dependence of the CCSD(T) method is well described by perturbation theory. Starting with CCSD共T兲/aug-cc-pVTZ calculations, use of the MP3 method to simulate the effect of increasing the basis set to aug-cc-pV5Z leads to average absolute errors, relative to the full CCSD共T兲/aug-cc-pV5Z calculations, of less than ⫾0.4 kcal/mol (D e ), ⫾0.0002 Å (r e ), ⫾2 cm⫺1 ( ␻ e ), 0.1 kcal/mol (I P e ), and 0.2 kcal/mol (EA e ) for the test set of diatomic molecules considered here. Although the corresponding MP2 approximation does not provide this high level of accuracy, it also should be useful for many molecular studies. When properly implemented, the savings in computer time should be significant since the MP3 method formally scales as N 6 , while the MP2 method scales as only N 5 . © 2000 American Institute of Physics. 关S0021-9606共00兲31942-0兴

I. INTRODUCTION

large basis sets must be used and/or a series of calculations must be extrapolated to the complete basis set limit. This is unfortunate as CCSD共T兲 calculations formally scale as N 7 , where N is the number of basis functions. Thus using a basis set with double the number of functions in the set or doubling the number of atoms in the molecule with the same basis set increases the cost of the calculation by two orders of magnitude. Although work is underway to reduce the dependence of CCSD共T兲 on N through a series of controlled approximations,8 even in the best of circumstances it is unlikely that the dependence will be reduced below N4–5 for many of the molecules of interest in chemistry. This steep dependence of the CCSD共T兲 method on the number of basis functions greatly restricts the range of applicability of this otherwise promising theoretical approach. Before continuing, one additional fact should be acknowledged. To achieve the accuracies stated above, it is necessary to include the effects of correlating the core electrons 共which requires even larger basis sets!兲 as well as scalar and spin-orbit relativistic effects. For first row atoms, these corrections are small relative to the errors introduced by truncation of the basis set. For atoms further down in the periodic chart, these effects, especially relativistic effects, become much more important. Even in these cases, however, truncation of the basis set is still a major problem. It is the errors introduced by basis set truncation that are of interest here. In examining the rate of convergence of a number of molecular properties with basis set, e.g., D e , r e , and ␻ e , we noted that:

During the past decade it has been established that coupled cluster methods provide a theoretically sound as well as rapidly convergent description of the electronic structure of atoms and molecules. The coupled cluster method that includes all single and double excitations plus a perturbative estimate for triple excitations—CCSD共T兲1—provides unrivaled accuracy for a wide range of molecules 共see the review by Lee and Scuseria,2 as well as the more recent papers of Dunning and co-workers,3 Martin and co-workers,4 Helgaker and co-workers,5 and Feller and co-workers6兲. For molecules well described by a Hartree–Fock wave function, the CCSD共T兲 method predicts bond dissociation energies, ionization potentials, and electron affinities to an accuracy of approximately ⫾0.5 kcal/mol, bond lengths accurate to ⫾0.0005 Å, and vibrational frequencies accurate to ⫾5 cm⫺1. Further, Feller7 has recently shown that, for N2 , CO, and HF, the differences between the bond dissociation energies calculated by full CI and the CCSD共T兲 method are ⫺0.4 ⫾0.25, ⫺0.4⫾0.3, and ⫺0.05⫾0.15 kcal/mol. Thus for D e , the results obtained with the CCSD共T兲 method are converged to within the error bars noted above. To obtain high accuracy in CCSD共T兲 calculations, very a兲

Current address: Office of Science, U.S. Department of Energy, 19901 Germantown Road, Germantown, MD 20874. b兲 Also at: Department of Chemistry, Washington State University, Richland, WA 99352. c兲 The Pacific Northwest National Laboratory is operated by Battelle Memorial Institute for the U.S. Department of Energy under Contract No. DEAC06-76RLO 1830. 0021-9606/2000/113(18)/7799/10/$17.00

7799

© 2000 American Institute of Physics

Downloaded 22 Feb 2001 to 192.101.100.146. Redistribution subject to AIP copyright, see http://ojps.aip.org/jcpo/jcpcpyrts.html

7800

J. Chem. Phys., Vol. 113, No. 18, 8 November 2000

T. H. Dunning, Jr. and K. A. Peterson

TABLE I. List of molecules included in the current study.

In the next section we provide the computational details. In the following two sections we first assess the accuracy of the various combinations of 共perturbation theory method, basis set兲 for the test set of molecules, and then determine the accuracy of the resulting ‘‘recommended’’ approach for all of the molecules in the test set. Finally, in the last section, we summarize the most important conclusions that can be drawn from the current work.

Neutrals AH AB

BH, CH, NH, OH, HF C2 , N2 , O2 , F2 , BF, CN, CO, NO

Cations AH ⫹ AB ⫹

BH⫹ , CH⫹ , NH⫹ , OH⫹ , HF⫹ ⫹ ⫹ O⫹ 2 , CO , NO

Anions AH ⫺ AB ⫺

II. COMPUTATIONAL DETAILS CH⫺ , NH⫺ , OH⫺ ⫺ O⫺ 2 , CN

⌬Q L-S i ⫽Q L 关 CCSD共T兲 ]⫺Q S i 关 CCSD共T兲 ] ⬇Q L 共 MPn 兲 ⫺Q S i 共 MPn 兲 ,

共1兲

where Q refers to the property, ‘‘L’’ to the largest basis set used, ‘‘S i ’’ to the series of smaller basis sets used, and MPn to nth order Møller–Plesset perturbation theory 共‘‘L’’ can also refer to the complete basis set limit; however, the analysis is then complicated by the differences in the extrapolation errors for the various methods兲. That is, the dependence of the property ‘‘Q’’ on the basis set was found to be approximately the same for perturbation theory methods as for the far more accurate CCSD共T兲 method. This suggests that it might be possible to obtain accurate CCSD共T兲 energies by performing CCSD共T兲 calculations in a small basis set and then correcting those energies by adding the difference between the energies obtained in MPn calculations with the same small basis set and the larger, more accurate basis set, i.e., E L 关 CCSD共T兲 ]⬇E S i 关 CCSD共T兲 ]⫹ 关 E L 共 MPn 兲 ⫺E S i 共 MPn 兲兴 . 共2兲 This approximation is also a key assumption in the Gn methods (n⫽2,3) of Curtiss, Raghavarchari, Pople, and co-workers.9 Since MP2 formally scales as N 5 , MP3 as N 6 , and MP4 as N 7 , use of the approximation in Eq. 共2兲 could result in significant computational savings if either the MP2 or MP3 methods result in acceptable errors. In this paper we examine the use of Eq. 共2兲 to approximate CCSD共T兲 calculations. Specifically, we use the energies obtained from Eq. 共2兲 to calculate the minimum energies (E e ), spectroscopic properties (D e ,r e , ␻ e , ␻ e x e ), ionization potentials (I P e ), and electron affinities (EA e ) of a number of diatomic molecules. This set of molecules was selected to be representative of a wide range of molecule types: ionic and covalent, open-shell and closed-shell, and those well described or poorly described by Hartree–Fock. The molecules included in the set are listed in Table I. For the present exploratory calculations, we take the large basis set 共L兲 to be the aug-cc-pV5Z basis set of Peterson and Dunning10 and determine how well Eq. 共2兲 reproduces the results of CCSD共T兲/aug-cc-pV5Z calculations when the small basis set ranges from aug-cc-pVDZ to aug-cc-pVQZ11 ( 兵 S i 其 ) and MPn ranges from MP2 to MP4.

Calculations on the molecules listed in Table I used the augmented correlation consistent polarized valence (aug-cc-pVnZ) basis sets of Kendall et al.11 and Peterson and Dunning.10 All CCSD共T兲 and closed-shell MPn calculations were carried out with the MOLPRO program package.12 The MPn calculations on the open-shell atoms and molecules were carried out with ACES II.13 The energies of the open-shell atoms and molecules were computed with the spin-restricted coupled cluster methods of Knowles et al.14 and the restricted many-body perturbation theory of Lauerdale et al.15 Spherical harmonics were used for the angular parts of the Gaussian basis functions. The calculations were performed on an SGI Origin 2000 at PNNL and on CRAY J90’s in DOE’s National Energy Research Scientific Computing Center at Lawrence Berkeley National Laboratory. The energies of the diatomic molecules were computed at seven internuclear separations around the equilibrium internuclear distance. Potential energy curves were obtained by fitting polynomials of sixth order through the computed energies. A standard Dunham analysis16 was employed to derive the minimum energy E e , bond dissociation energy D e , equilibrium distance r e , harmonic frequency ␻ e , and anharmonic correction ␻ e x e from the polynomial coefficients. The equilibrium ionization potentials, I P e , and electron affinities, EA e , are simply the differences in E e ’s for the ions and corresponding neutral molecules. CCSD共T兲 calculations were first carried out with the aug-cc-pV5Z set. This defined the reference results. Then, the CCSD共T兲/aug-cc-pV5Z results were approximated using Eq. 共2兲. The difference between these results and the reference results defined the approximation errors: ⌬E e , ⌬D e , ⌬r e , ⌬ ␻ e , ⌬ ␻ e x e , ⌬I P e , and ⌬EA e . We also carried out a corresponding series of calculations with the standard correlation consistent basis sets.17,10 For molecules that are well described by the standard sets, the conclusions are essentially the same as those drawn here for the augmented sets. Therefore, these results are not reported. III. EVALUATION OF PERTURBATION THEORY APPROXIMATIONS: CCSD METHOD

Although only the CCSD共T兲 method provides high accuracy for molecular calculations, it is important to also check the accuracy of the perturbation theory approximations to the CCSD method. The errors resulting from use of the various MPn approximations to describe the basis set dependence of the CCSD method are summarized in Table II for the total energies (⌬E e ) and spectroscopic constants

Downloaded 22 Feb 2001 to 192.101.100.146. Redistribution subject to AIP copyright, see http://ojps.aip.org/jcpo/jcpcpyrts.html

J. Chem. Phys., Vol. 113, No. 18, 8 November 2000

Basis set dependence

7801

TABLE II. Errors in the approximate CCSD energies (⌬E e ) and spectroscopic constants (⌬D e ,⌬r e ,⌬ ␻ e ) for the molecules included in the current study. Errors are relative to results obtained from CCSD calculations with the aug-cc-pV5Z set. MP2 Molecule set

a v dz

MP3 a v tz

a v qz

a v dz

a v tz

MP4 a v qz

a v dz

a v tz

a v qz

1.902

0.485

0.124

10.774

2.532

0.685

⌬E e (mE e ) Average Absolute Errors A AH AB Maximum Absolute Errors A AH AB

4.593 5.749 11.713

2.704 3.299 6.468

1.113 1.357 2.751

1.738 1.638 2.853

0.487 0.481 1.087

0.096 0.091 0.201

5.125 6.220 15.253

3.142 3.783 7.120

1.615 1.860 3.505

1.984 1.888 3.453

0.724 0.727 1.773

0.205 0.184 0.476

2.990

0.891

0.259

14.017

3.368

0.886

⌬D e (kcal/mol) Average Absolute Errors AH AB Maximum Absolute Errors AH AB

0.681 1.747

0.306 0.720

0.126 0.350

0.067 0.484

0.019 0.103

0.008 0.046

4.166

0.925

0.273

1.347 3.139

0.438 1.983

0.189 0.995

0.104 0.849

0.040 0.244

0.013 0.160

6.968

1.746

0.480

⌬r e (mÅ) Average Absolute Errors AH AB Maximum Absolute Errors AH AB

1.41 2.66

0.25 0.91

0.08 0.23

0.33 0.58

0.11 0.21

0.02 0.06

3.77

0.72

0.18

1.92 5.78

0.35 1.88

0.10 0.54

0.48 2.64

0.17 0.91

0.04 0.19

7.95

1.58

0.40

⫺1

⌬ ␻ e (cm Average Absolute Errors AH AB Maximum Absolute Errors AH AB

)

15.47 16.69

4.20 5.23

0.93 1.32

4.79 3.81

2.04 1.52

0.41 0.35

26.15

5.35

1.32

20.65 23.79

5.83 7.18

1.85 1.86

7.54 12.26

3.18 3.99

0.75 1.15

46.72

8.91

2.57

(⌬D e ,⌬r e ,⌬ ␻ e ). Both the average absolute errors and the maximum absolute errors are listed. In the tables we use the shorthand notation: a v nz⫽aug-cc-pVnZ 共n⫽d/D, t/T, q/Q兲. Double excitations and products of double excitations are the dominant contributions to the CCSD wave function. The MP2 and MP3 methods involve only double excitations and, thus may be expected to provide reliable approximations to the CCSD method. The presence of triple excitations in the MP4 wave function can be expected to lead to a loss of accuracy. To check this expectation, the MP4 method was used to approximate the CCSD energies of the AB molecules as well as those of the first row atoms. The results, which are also reported in Table II, confirm this expectation. Errors in the MP4 approximation are usually larger than the errors resulting from use of the MP2 approximation and are always larger than from the MP3 approximation. This difference is most pronounced for the MP3 and MP4 maximum absolute errors. Below we only discuss the results obtained with the MP2 and MP3 approximations. Consider the errors in the total energies—a very demanding criterion for assessing the accuracy of the perturba-

tion theory approximations. First, we note that the error decreases significantly as the basis set is expanded from augcc-pVDZ to aug-cc-pVTZ to aug-cc-pVQZ. This is to be expected—at n⫽5, the error is exactly zero; see Eq. 共2兲. Second, it is quite clear that the MP3 approximation is far more accurate than the MP2 approximation. Differences in the errors of factors of 5⫻ or more are not uncommon. Finally, note that the errors in the atomic energies are approximately the same as those in the AH energies and about half those in the AB energies. Since the signs of the errors are the same, much of the error in the absolute energies will cancel when differences are taken; compare, e.g., the errors in E e versus those in D e . We see similar trends for the approximated spectroscopic constants as for the total energies. In general, we find the errors to be significantly larger for the AB molecules than for the AH molecules. For example, for the AH molecules the average absolute error in r e is 0.11 mÅ for the MP3/augcc-pVTZ approximation, whereas it is twice that value 共0.21 mÅ兲 for the AB molecules. We again find that the error decreases as n (a v nz) increases. For D e , the error in the MP2 approximation for the AB molecules decreases by a

Downloaded 22 Feb 2001 to 192.101.100.146. Redistribution subject to AIP copyright, see http://ojps.aip.org/jcpo/jcpcpyrts.html

7802

J. Chem. Phys., Vol. 113, No. 18, 8 November 2000

T. H. Dunning, Jr. and K. A. Peterson

TABLE III. Errors in the approximate CCSD共T兲 total energies, ⌬E e , for the atoms and molecules included in the current study. Errors are relative to results obtained from CCSD共T兲 calculations with the aug-cc-pV5Z set. Errors in total energies are given in millihartrees. MP2

MP3

MP4

a v dz

a v tz

a v qz

a v dz

a v tz

a v qz

a v dz

a v tz

a v qz

Average Absolute Errors A 2.813 AH 3.227 AB 4.163

2.225 2.683 4.642

0.976 1.188 2.235

0.674 1.146 4.762

0.059 0.135 0.764

0.050 0.095 0.337

0.122 0.303 3.266

0.058 0.073 0.702

0.020 0.017 0.169

Maximum Absolute Errors A 3.880 AH 4.420 AB 8.113

2.508 2.947 4.927

1.366 1.551 2.807

1.335 3.046 6.764

0.084 0.348 1.277

0.067 0.125 0.534

0.189 0.538 6.876

0.083 0.131 1.562

0.033 0.025 0.413

Molecule set

factor of slightly more than 2 from the aug-cc-pVDZ to the aug-cc-pVTZ set and by another factor of 2 from the aug-ccpVTZ to the aug-cc-pVQZ set. A similar trend is observed for the MP3 approximation, although the overall decrease is a factor of 10, rather than a factor of 5. For both r e and ␻ e , the errors tend to decrease more rapidly with increasing n than for D e . The MP3 approximation is also found to be substantially more accurate than the MP2 approximation for the spectroscopic constants. For example, with the aug-cc-pVTZ set, the average absolute error in D e for the MP3 approximation is 16 (AH) and 7 (AB) times smaller than the corresponding errors for the MP2 approximation. The reductions are nearly as large for the maximum absolute errors; the factors are 11 and 8, respectively. The error ratios are smaller for r e and ␻ e , although still significant. The average absolute errors are decreased by factors ranging from 2.1 to 4.3, while the maximum absolute errors are reduced by factors ranging from 1.8 to 2.1. IV. EVALUATION OF PERTURBATION THEORY APPROXIMATIONS: CCSD„T… METHOD

The errors resulting from use of the various MPn approximations to describe the basis set dependence of the CCSD共T兲 method are summarized in Tables III 共total energies: ⌬E e 兲, IV 共spectroscopic constants: ⌬D e , ⌬r e , ⌬ ␻ e , ⌬ ␻ e x e 兲, and V 共ionization potentials and electron affinities: ⌬I P e , ⌬EA e 兲. The average absolute errors in D e , r e , and ␻ e are plotted in Figs. 1 and 2 for the AH and AB sets; the average absolute errors in the I P e ’s and EA e ’s are plotted in Figs. 3 and 4, respectively. Also included in the figures are the average absolute errors associated with the corresponding CCSD共T兲/aug-cc-pVnZ calculations 共heavy solid lines兲 as well as the intrinsic errors relative to experiment discussed in the Introduction 共dashed lines兲. Consider first the errors in the total energies reported in Table III. As can be seen, the errors for the MP2 method are fairly large. Even for the aug-cc-pVQZ basis set, the average absolute error for the AB molecule set exceeds 2 millihartrees (mE h ) and the maximum error is nearly 3 mE h . Again, it should be noted that the average absolute error for the AH molecule set is approximately the same as that for the A atom set 共1.2 vs 1.0 mE h 兲 and the average absolute error for the

AB molecule set is approximately twice that of the A atom set 共2.2 vs 1.0 mE h 兲. Thus the errors, although large, will again cancel when energy differences are taken. There is a dramatic improvement in both the average and maximum absolute errors listed in Table III when the MP3 method is used. For the AB molecule set, the average absolute error for the MP3 method with the aug-cc-pVTZ set is three times smaller than for the MP2 method with the augcc-pVQZ set. Use of the aug-cc-pVQZ set reduces the error by another factor of 2. This trend continues with the MP4 method, although the decrease is in general less dramatic. One point that should be noted is that for both the MP3 and MP4 methods, the errors in the atomic calculations are markedly smaller than half the errors in the AB calculations. Thus there is less tendency for the errors to cancel when taking energy differences. One further point to note in Table III is that the maximum absolute errors are rarely more than twice the average absolute errors. The largest difference is for the calculations on the AH molecule set with the aug-cc-pVDZ basis set, where the ratio is a factor of 2.7. This is encouraging. A method that yields maximum absolute errors not much larger than the average absolute errors is obviously superior to a method that yields maximum absolute errors that are much larger than the average absolute errors even if the average absolute errors are smaller. Now, let us compare the average absolute errors from CCSD共T兲/aug-cc-pVnZ calculations for D e , r e , and ␻ e with those obtained by using the various MPn approximations to describe the basis set dependence of the CCSD共T兲 method. That is, we will compare the results obtained using only the first term on the right hand side of Eq. 共2兲 with the results obtained using the complete expression on the right hand side. For the AH molecules 共see Fig. 1兲, use of the MPn approximation reduces the error by approximately an order of magnitude in the worst case and by nearly two orders of magnitude in the best case. For D e , use of the MP2 approximation yields less accurate results than use of the MP3 approximation, which, in turn, yields less accurate results than the MP4 approximation. This seems intuitively reasonable— the MP3 method includes the effects of interactions between double excitations as is the case in the CCSD part of the CCSD共T兲 method 共see the last section兲, while the MP4

Downloaded 22 Feb 2001 to 192.101.100.146. Redistribution subject to AIP copyright, see http://ojps.aip.org/jcpo/jcpcpyrts.html

J. Chem. Phys., Vol. 113, No. 18, 8 November 2000

Basis set dependence

7803

TABLE IV. Errors in the approximate CCSD共T兲 spectroscopic constants (⌬D e ,⌬r e ,⌬ ␻ e ,⌬ ␻ e x e ) for the molecules included in the current study. Errors are relative to results obtained from CCSD共T兲 calculations with the aug-cc-pV5Z set. MP2 a v dz

Molecule set

MP3 a v tz

a v qz

a v dz

MP4

a v tz

a v qz

a v dz

a v tz

a v qz

⌬D e (kcal/mol) Average Absolute Errors AH AB Maximum Absolute Errors AH AB

0.301 0.862

0.237 0.163

0.110 0.129

0.530 2.854

0.089 0.486

0.029 0.152

0.117 1.873

0.031 0.425

0.010 0.119

0.369 1.481

0.369 0.317

0.171 0.279

1.073 3.675

0.166 0.902

0.052 0.280

0.272 4.097

0.052 1.019

0.014 0.288

⌬r e (mÅ) Average Absolute Errors AH AB Maximum Absolute Errors AH AB

0.80 1.36

0.14 0.71

0.04 0.17

0.43 1.30

0.02 0.23

0.02 0.09

0.32 2.62

0.08 0.52

0.02 0.12

1.29 2.23

0.19 1.19

0.05 0.35

0.83 2.77

0.03 0.49

0.03 0.20

0.48 4.97

0.12 0.85

0.03 0.22

⌬ ␻ e (cm⫺1 ) Average Absolute Errors AH AB Maximum Absolute Errors AH AB

7.86 8.34

2.65 3.40

0.57 0.73

3.34 8.65

0.47 1.96

0.22 0.73

5.28 18.81

1.29 3.67

0.26 0.81

12.93 13.66

3.78 4.40

1.08 1.11

9.19 15.47

1.04 5.46

0.43 1.82

8.32 36.71

2.22 7.24

0.51 2.11

⌬ ␻ e x e (cm⫺1 ) Average Absolute Errors AH AB Maximum Absolute Errors AH AB

0.61 0.25

0.17 0.03

0.11 0.02

0.46 0.08

0.13 0.02

0.03 0.01

0.38 0.45

0.11 0.08

0.03 0.02

0.97 1.09

0.35 0.07

0.16 0.04

0.96 0.10

0.24 0.05

0.09 0.01

0.85 1.39

0.25 0.18

0.06 0.06

TABLE V. Errors in the approximate CCSD共T兲 ionization potentials (⌬I P e ) and electron affinities (⌬EA e ) for the molecules included in the current study. Errors are relative to results obtained from CCSD共T兲 calculations with the aug-cc-pV5Z set. MP2 Molecule set

a v dz

a v tz

MP3 a v qz

a v dz

a v tz

MP4 a v qz

a v dz

a v tz

a v qz

⌬I P e 共in kcal/mol兲 Average Absolute Errors 0.440 AH AB 0.146 Maximum Absolute Errors 1.611 AH AB 0.353

0.270 0.211

0.130 0.130

0.609 0.042

0.113 0.034

0.040 0.018

0.089 0.345

0.021 0.089

0.007 0.048

0.809 0.348

0.250 0.167

1.271 0.084

0.185 0.057

0.051 0.024

0.235 0.438

0.029 0.136

0.013 0.087

⌬EA e 共in kcal/mol兲 Average Absolute Errors 0.492 AH AB 0.453 Maximum Absolute Errors 0.553 AH AB 0.587

0.326 0.351

0.150 0.193

1.014 0.107

0.267 0.088

0.112 0.030

0.401 0.193

0.096 0.053

0.033 0.032

0.452 0.420

0.208 0.236

1.584 0.154

0.360 0.093

0.145 0.031

0.704 0.227

0.163 0.079

0.064 0.046

Downloaded 22 Feb 2001 to 192.101.100.146. Redistribution subject to AIP copyright, see http://ojps.aip.org/jcpo/jcpcpyrts.html

7804

J. Chem. Phys., Vol. 113, No. 18, 8 November 2000

T. H. Dunning, Jr. and K. A. Peterson

FIG. 1. Average absolute errors for the AH molecule set for MPn approximations to CCSD共T兲/aug-cc-pV5Z calculations: 共a兲 ⌬D e ; 共b兲 ⌬r e ; 共c兲 ⌬ ␻ e . Also plotted are the average absolute errors for CCSD共T兲/aug-cc-pVnZ calculations 共heavy solid lines兲 as well as the intrinsic errors for the CCSD共T兲 method 共dashed lines兲; see the text.

method includes the effects of triple excitations, which are also included in the CCSD共T兲 method. For r e and ␻ e , on the other hand, the results obtained with the MP3 method are nearly as accurate as, and often more accurate than, those obtained with the MP4 method. For the AB molecule set, use of the MPn approximation also reduces the error by approximately an order of magnitude. For D e , the MP2 approximation actually leads to the most accurate results with the MP3 and MP4 approximations

yielding very similar results 共the MP4 method provides slightly more accurate results than the MP3 method兲. For r e and ␻ e , the MP3 approximation yields the most accurate results with the MP2 and MP4 methods often providing similar results. It is surprising, and perhaps even counterintuitive, that for both the AH and AB molecule sets, the MP3 method yields results for r e and ␻ e as good as and often better than the MP4 method. However, in a study of the convergence of

FIG. 2. Average absolute error for the AB molecule set for MPn approximations to CCSD共T兲/aug-cc-pV5Z calculations: 共a兲 ⌬D e ; 共b兲 ⌬r e ; 共c兲 ⌬ ␻ e . Also plotted are the average absolute errors for CCSD共T兲/aug-cc-pVnZ calculations 共heavy solid lines兲 as well as the intrinsic errors for the CCSD共T兲 method 共dashed lines兲; see text.

Downloaded 22 Feb 2001 to 192.101.100.146. Redistribution subject to AIP copyright, see http://ojps.aip.org/jcpo/jcpcpyrts.html

J. Chem. Phys., Vol. 113, No. 18, 8 November 2000

Basis set dependence

7805

FIG. 3. Average absolute errors, ⌬I P e , for the 共a兲 AH and 共b兲 AB molecule sets for MPn approximations to CCSD共T兲/aug-cc-pV5Z calculations 共in kcal/mol兲. Also plotted are the average absolute errors for CCSD共T兲/aug-cc-pVnZ calculations 共heavy solid lines兲 as well as the intrinsic errors for the CCSD共T兲 method 共dashed lines兲; see text.

the perturbation expansion for molecular spectroscopic constants, Dunning and Peterson18 also noted a marked difference in the behavior of D e and 兵 r e , ␻ e 其 . Equation 共2兲 can also be recast as E L 关 CCSD共T兲 ]⬇E S i 关 CCSD兴 ⫹ 关 E L 共 MPn 兲 ⫺E S i 共 MPn 兲兴 ⫹ 兵 E S i 关 CCSD共T兲兴 ⫺E S i 共 CCSD兲 其 ;

共3兲

that is, we can regard the perturbation theory approximation to the basis set dependence of the energy as a correction to the CCSD energy, which is then corrected for the effect of triple excitations. If changes in the CCSD energy with basis set are substantially larger than the corresponding changes in the triples correction, one would expect the MP3 approximation to be more accurate than the MP4 approximation as shown in the last section. In Eq. 共3兲, proper choice of S i is

required to ensure that both the MPn approximation and the triples correction are converged to the desired accuracy. For both the AH and AB molecule sets, there is a nearly exponential decrease in the average absolute error with increasing basis set size. The MP2 method tends to behave less monotonically than the MP3 and MP4 methods, although nonmonotonic behavior can also be observed for the MP3 method 共see, e.g., ⌬r e for the AH set兲. Many of the same conclusions can be drawn about the calculated ionization potentials and electron affinities 共see Figs. 3 and 4兲. Here, however, with the exception of the EA e ’s for the AH molecule set calculated with the aug-ccpVQZ basis set, the MP2 approximation is markedly less accurate than the MP3 and MP4 approximations. In addition, the MP2 approximation often does not show significant im-

FIG. 4. Average absolute errors, ⌬EA e , for the 共a兲 AH and 共b兲 AB molecule sets for MPn approximations to CCSD共T兲/aug-cc-pV5Z calculations 共in kcal/mol兲. Also plotted are the average absolute errors for CCSD共T兲/aug-cc-pVnZ calculations 共heavy solid lines兲 as well as the intrinsic errors for the CCSD共T兲 method 共dashed lines兲; see text.

Downloaded 22 Feb 2001 to 192.101.100.146. Redistribution subject to AIP copyright, see http://ojps.aip.org/jcpo/jcpcpyrts.html

7806

J. Chem. Phys., Vol. 113, No. 18, 8 November 2000

provement with increasing basis set size 共increasing n兲. Since the errors in the CCSD共T兲/aug-cc-pVnZ calculations decrease steadily with increasing basis set size, the improvement resulting from use of the MP2 method to approximate the basis set dependence also decreases steadily with basis set. For the approximation in Eq. 共2兲 to be most useful, the errors introduced by the approximation should be comparable to or even less than the intrinsic errors in the CCSD共T兲 method itself. The intrinsic errors 共D e , I P e , EA e : 0.5 kcal/ mol; r e : 0.0005 Å; ␻ e : 5 cm⫺1兲 are plotted as dashed lines in Figs. 1–4. For the AH molecule set, essentially all of the errors except those obtained with the aug-cc-pVDZ set fall below the dashed lines. For the aug-cc-pVDZ set, the approximation errors lie slightly above the intrinsic errors for both r e , ␻ e , and EA e . The same trend holds for the AB molecule set, although the errors resulting from use of the aug-cc-pVDZ set are now well above the intrinsic errors for D e , r e , and ␻ e . The reason for the poor performance of the aug-cc-pVDZ set is a result of the errors associated with the CCSD共T兲/aug-cc-pVDZ calculations. The average absolute errors for CCSD共T兲 calculations with the aug-cc-pVDZ basis set are three to four times larger than those from calculations with the aug-cc-pVTZ set. The use of perturbation theory to approximate the basis set dependence of the CCSD共T兲 method is simply not able to compensate for this increased error.

V. RESULTS OBTAINED WITH „MP3Õaug-cc-pVTZ… APPROXIMATION

From the results summarized in the preceding sections it would appear that a reasonable compromise between accuracy and computational cost in approximating the basis set dependence of the CCSD共T兲 method is to use the MP3 method with an aug-cc-pVTZ basis set; that is, ˜E a v 5z 关 CCSD共T 兲 ]⫽E a v tz 关 CCSD共T 兲 ] ⫹ 关 E a v 5z 共 MP3 兲 ⫺E a v tz 共 MP3 兲兴 . 共4兲 This approximation leads to average absolute errors that are less than the intrinsic errors associated with the CCSD共T兲 method itself. It replaces an N 7 CCSD共T兲 calculation with the aug-cc-pV5Z basis set 共127 functions兲 with an N 7 CCSD共T兲 calculation with the aug-cc-pVTZ set 共46 functions兲 plus two N 6 calculations, one with the aug-cc-pV5Z set and the other with the aug-cc-pVTZ set. CCSD共T兲 calculations with the aug-cc-pVTZ set are as much as three orders of magnitude cheaper than such calculations with the aug-ccpV5Z set and, so, more than compensate for the cost of the two additional N 6 calculations. The errors resulting from the use of Eq. 共4兲, relative to those obtained from CCSD共T兲/aug-cc-pV5Z calculations, are listed in Tables VI (⌬D e ,⌬r e ,⌬ ␻ e ,⌬ ␻ e x e ) and VII (⌬I P e ,⌬EA e ) for all of the molecules considered here. Also included in Tables VI and VII are the errors for the base CCSD共T兲/aug-cc-pVTZ calculations. Comparing the columns labeled ‘‘CCSD共T兲’’ and ‘‘CCSD共T兲⫹MP3’’ allows

T. H. Dunning, Jr. and K. A. Peterson

the impact of the MP3 approximation in Eq. 共4兲 to be quickly quantified. The improvements are dramatic. For D e , the errors are reduced by a factor of 5 (C2 ) to nearly a factor of 40 共NH兲; for I P e , the improvements span nearly the same scale, from a factor of 5 共NH兲 to 30 共BH兲. For EA e , the improvements are not as large, but are still significant, varying from a factor of 3 共CH兲 to more than 15 共CN兲. For both r e and ␻ e , substantial improvements are found, varying from a low of 4 关 ␻ e (CN) 兴 to a high of more than 300 关 r e (OH) 兴 . All in all, there can be little doubt that the use of Eq. 共3兲 dramatically decreases the errors in the uncorrected CCSD共T兲/aug-ccpVTZ calculations. Examining Table VI more carefully, we see that the largest errors resulting from use of the MPn approximation are associated with molecules that are poorly described by a Hartree–Fock wave function. For D e , the largest error is for C2 共⫺0.9 kcal/mol兲, a molecule that is known to have substantial multireference character in the wave function. C2 is followed by O2 , which has an error of 0.65 kcal/mol, another difficult molecule to describe with a single reference wave function. For both r e and ␻ e , the largest errors are for CN 共⫺0.49 mÅ, 5.45 cm⫺1兲, which, like C2 , is best described by a multireference wave function. In both cases CN is followed by O2 共⫺0.45 mÅ, 3.36 cm⫺1兲. The coupled cluster and perturbation theory methods used here are both based on Hartree–Fock wave functions. However, past work has clearly established that the range of applicability of coupled cluster methods far exceeds that of perturbation theory methods. Thus it is not surprising that the MPn approximation to the basis set dependence of coupled cluster theory is less accurate for molecules that are poorly described by a Hartree–Fock wave function. Although the above recommends use of the MP3 method to correct for basis set truncation in CCSD共T兲 calculations, it should be noted that the MP2 method, when combined with the aug-cc-pVTZ set, also improves the accuracy of the base CCSD共T兲/aug-cc-pVTZ calculations. The average accuracies of the MP2 and MP3 approximations can be readily compared in Tables IV and V. From these tables, we see that the MP3 approximation is on the average significantly more accurate than the MP2 approximation. However, with the single exception of ⌬r e for the AB molecule set, the average absolute errors obtained using the MP2 approximation also fall below the intrinsic errors associated with the CCSD共T兲 method 共there are, however, larger variations in the errors themselves兲. The MP2 approximation, while clearly not as accurate as the MP3 approximation, may be adequate for many computational studies. The MP2 method scales as N 5 and approaches have already been developed that reduce this scaling to N 3 , although with some loss of accuracy for smaller basis sets19 共see also Ref. 20兲.

VI. CONCLUSIONS

The CCSD共T兲 method—a coupled cluster method that includes all single and double excitations plus a perturbative estimate of triple excitations—provides an accurate description of the electronic structure of a broad range of molecules. To achieve high accuracy, however, large basis sets must be

Downloaded 22 Feb 2001 to 192.101.100.146. Redistribution subject to AIP copyright, see http://ojps.aip.org/jcpo/jcpcpyrts.html

J. Chem. Phys., Vol. 113, No. 18, 8 November 2000

Basis set dependence

7807

TABLE VI. Errors in the spectroscopic constants (⌬D e ,⌬r e ,⌬ ␻ e ,⌬ ␻ e x e ) obtained from CCSD共T兲 calculations with the aug-cc-pVTZ set compared to the errors obtained from the same calculations with an MP3 correction. Errors are relative to results obtained from CCSD共T兲 calculations with the aug-cc-pV5Z set. ⌬D e (kcal/mol)

⌬ ␻ e (cm⫺1 )

⌬r e (mÅ)

⌬ ␻ e x e (cm⫺1 )

Molecule

CCSD共T兲

CCSD共T兲⫹MP3

CCSD共T兲

CCSD共T兲⫹MP3

CCSD共T兲

CCSD共T兲⫹MP3

CCSD共T兲

CCSD共T兲⫹MP3

BH CH NH OH HF C2 N2 O2 F2 BF CN CO NO

⫺0.96 ⫺1.31 ⫺1.85 ⫺1.93 ⫺1.93 ⫺4.57 ⫺7.60 ⫺4.02 ⫺1.64 3.32 ⫺5.74 ⫺5.51 ⫺5.90

⫺0.06 ⫺0.06 ⫺0.05 ⫺0.11 ⫺0.17 ⫺0.90 ⫺0.45 ⫺0.65 ⫺0.28 ⫺0.26 ⫺0.40 ⫺0.47 ⫺0.47

2.70 2.24 2.55 3.12 3.66 5.96 4.44 6.02 7.03 7.23 4.82 5.08 5.08

⫺0.03 ⫺0.02 ⫺0.02 0.01 ⫺0.02 0.16 ⫺0.25 ⫺0.45 0.03 0.21 ⫺0.49 ⫺0.07 ⫺0.17

⫺11.0 ⫺14.0 ⫺20.2 ⫺24.0 ⫺17.6 ⫺17.7 ⫺19.8 ⫺25.0 ⫺10.6 ⫺11.6 ⫺21.2 ⫺19.3 ⫺20.1

0.05 ⫺0.21 ⫺0.64 ⫺1.04 ⫺0.41 ⫺0.95 2.86 3.36 ⫺0.19 ⫺0.74 5.45 0.79 1.34

⫺0.15 1.16 1.25 2.31 2.47 0.09 0.03 ⫺0.14 0.58 ⫺0.04 0.05 ⫺0.11 0.15

0.01 0.12 ⫺0.20 ⫺0.17 ⫺0.17 0.00 ⫺0.02 ⫺0.04 0.01 0.00 ⫺0.08 0.00 0.04

used. Since the computational cost of the CCSD共T兲 method increases as N 7 , where N is the number of basis functions used in the calculations, CCSD共T兲 calculations becomes prohibitively expensive for large molecules. In this study we examined the use of perturbation theory to predict the effect of increasing the size of the basis set in CCSD共T兲 calculations, using the augmented correlation consistent sets of Dunning and co-workers.10,11 We found that the spectroscopic properties (D e ,r e , ␻ e , ␻ e x e ) for a set of diatomic molecules calculated with the aug-cc-pV5Z basis set, as well as the associated ionization potentials (I P e ) and electron affinities (EA e ), were well reproduced by: ˜E a v 5z 关 CCSD共T 兲 ]⫽E a v tz 关 CCSD共T 兲 ]⫹ 关 E a v 5z 共 MP3 兲 ⫺E a v tz 共 MP3 兲兴 . The average absolute errors for the (AH⫹AB) molecule set resulting from the use of this equation, relative to the full CCSD共T兲/aug-cc-pV5Z calculations, are: 0.33 kcal/mol (D e ), 0.15 mÅ (r e ), 1.4 cm⫺1 ( ␻ e ), 0.06 cm⫺1 ( ␻ e x e ), 0.08 kcal/mol (I P e ), and 0.20 kcal/mol (EA e ). These errors are TABLE VII. Errors in the ionization potential (⌬I P e ) and electron affinity (⌬EA e ) obtained from CCSD共T兲 calculations with the aug-cc-pVTZ set compared to the errors obtained from the same calculations with an MP3 correction. Errors are relative to results obtained from CCSD共T兲 calculations with the aug-cc-pV5Z set. ⌬I P e (kcal/mol) Molecule

CCSD共T兲

CCSD共T兲⫹MP3

BH CH NH OH HF O2 CN CO NO

⫺0.60 ⫺0.84 ⫺0.85 ⫺2.51 ⫺2.35 ⫺0.78

⫺0.02 ⫺0.09 ⫺0.16 ⫺0.11 ⫺0.18 ⫺0.00

⫺1.38 ⫺0.80

⫺0.06 0.04

⌬EA e (kcal/mol) CCSD共T兲

less than the intrinsic errors in the CCSD共T兲 method itself, namely, ⌬D e , ⌬I P e , ⌬EA e : ⫾0.5 kcal/mol; ⌬r e : ⫾0.0005 Å, and ⌬ ␻ e : ⫾5 cm⫺1. Since MP3 calculations scale as N 6 , use of the above approximation should significantly increase the size of molecules that can be treated by the CCSD共T兲 method. The MP3 method also provides an excellent approximation of the basis set dependence of the CCSD method. The MP4 method, which includes triple excitations, does not provide as accurate a representation of the CCSD/aug-cc-pV52 energy as the MP3 method. This is at least partly responsible for the poor performance of the MP4 method in representing the basis set dependence of the CCSD共T兲 method. Although the MP2 approximation when combined with the aug-cc-pVTZ set provides a less accurate description of the energy than does use of the MP3 approximation, it also significantly reduces the errors in the CCSD共T兲/aug-ccpVTZ calculations. The MP2 method scales as N 5 , and variants with reduced scaling requirements currently exist.19 The 共MP2, aug-cc-pVTZ兲 approximation may be adequate for many computational studies. Finally, it should be noted that, by carrying out a sequence of calculations with the augmented correlation consistent basis sets, the second term on the right hand side can be extrapolated to the complete basis set limit. In this case, the above equation leads to an approximation to E CBS 关 CCSD共T)], the CCSD共T兲 energy at the complete basis set limit. ACKNOWLEDGMENTS

CCSD共T兲⫹MP3

⫺0.53 ⫺1.92 ⫺1.71

⫺0.19 ⫺0.25 ⫺0.36

⫺0.99 ⫺1.24

⫺0.09 0.08

This work was supported by the Division of Chemical Sciences in the Office of Basis Energy Sciences of the U.S. Department of Energy at Pacific Northwest National Laboratory, a multiprogram national laboratory operated by Battelle Memorial Institute, under Contract No. DE-AC0676RLO 1830. Computer resources were provided by the Division of Chemical Sciences and by the Office of Scientific Computing, at the National Energy Research Scientific Computing Center 共NERSC兲 at Lawrence Berkeley National

Downloaded 22 Feb 2001 to 192.101.100.146. Redistribution subject to AIP copyright, see http://ojps.aip.org/jcpo/jcpcpyrts.html

7808

J. Chem. Phys., Vol. 113, No. 18, 8 November 2000

Laboratory. We thank the referee for stressing the importance of the CCSD extrapolations as well as the CCSD共T兲 extrapolations reported herein. G. D. Purvis and R. J. Bartlett, J. Chem. Phys. 76, 1910 共1982兲; K. Raghavachari, G. W. Trucks, J. A. Pople, and M. Head-Gordon, Chem. Phys. Lett. 157, 479 共1989兲. See also J. F. Stanton, ibid. 281, 130 共1997兲. 2 T. J. Lee and G. Scuseria, in Quantum Mechanical Electronic Structure Calculations with Chemical Accuracy, edited by S. Langhoff 共Kluwer, Dordrecht, 1995兲. 3 K. A. Peterson and T. H. Dunning, Jr., J. Chem. Phys. 102, 2032 共1995兲; K. A. Peterson and T. H. Dunning, Jr., J. Phys. Chem. 99, 3898 共1995兲; D. E. Woon, T. H. Dunning, Jr., and K. A. Peterson, J. Chem. Phys. 104, 5883 共1996兲; S. S. Xantheas, T. H. Dunning, Jr., and A. Mavridis, ibid. 106, 3280 共1997兲; K. A. Peterson and T. H. Dunning, Jr., J. Mol. Struct.: THEOCHEM 400, 93 共1997兲; T. van Mourik and T. H. Dunning, J. Chem. Phys. 107, 2451 共1997兲; K. A. Peterson, A. K. Wilson, D. E. Woon, and T. H. Dunning, Jr., Theor. Chem. Acc. 97, 251 共1997兲; D. E. Woon, K. A. Peterson, and T. H. Dunning, Jr., J. Chem. Phys. 109, 2233 共1998兲; T. van Mourik, A. K. Wilson, and T. H. Dunning, Jr., Mol. Phys. 96, 529 共1999兲; T. van Mourik and T. H. Dunning, Jr., J. Chem. Phys. 111, 9248 共1999兲. 4 J. M. L. Martin and P. R. Taylor, Chem. Phys. Lett. 248, 336 共1996兲; J. M. L. Martin and T. J. Lee, ibid. 258, 129 共1996兲; ibid. 258, 136 共1996兲; J. M. L. Martin, ibid. 259, 669 共1996兲; ibid. 259, 679 共1996兲; J. M. L. Martin, in Computational Thermochemistry. Prediction and Estimation of Molecular Thermodynamics, edited by K. K. Irikura and D. J. Frurip 共American Chemical Society, Washington, D.C., 1996兲, p. 212; J. M. L. Martin and P. R. Taylor, J. Chem. Phys. 106, 8620 共1997兲; J. M. L. Martin, Theor. Chem. Acc. 97, 227 共1997兲; J. M. L. Martin, Chem. Phys. Lett. 273, 98 共1997兲; ibid. 292, 411 共1998兲; J. M. L. Martin, T. J. Lee, and P. R. Taylor, J. Chem. Phys. 108, 676 共1998兲. 5 A. Halkier, O. Christiansen, D. Sundholm, and P. Pyykkoo, Chem. Phys. Lett. 271, 273 共1997兲; A. Halkier, P. Jorgensen, J. Gauss, and T. Helgaker, ibid. 274, 235 共1997兲; T. Helgaker, J. Gauss, P. Jørgensen, and J. Olsen, J. Chem. Phys. 106, 6430 共1997兲; T. Helgaker, W. Klopper, H. Koch, and J. Noga, ibid. 106, 9639 共1997兲; A. Halkier, H. Koch, P. Jorgensen, O. Christiansen, I. M. Beck Nielsen, and T. Helgaker, Theor. Chem. Acc. 97, 150 共1997兲; A. Halkier and P. R. Taylor, Chem. Phys. Lett. 285, 133 1

T. H. Dunning, Jr. and K. A. Peterson 共1998兲; A. Halkier, T. Helgaker, P. Jorgensen, W. Klopper, H. Koch, J. Olsen, and A. K. Wilson, ibid. 286, 243 共1998兲; A. Halkier, S. Coriani, and P. Jorgensen, ibid. 294, 292 共1998兲; A. Halkier, H. Larsen, J. Olsen, and P. Jorgensen, J. Chem. Phys. 110, 734 共1999兲; A. Halkier and S. Coriani, Chem. Phys. Lett. 303, 408 共1999兲; A. Halkier, W. Klopper, T. Helgaker, P. Jorgensen, and P. R. Taylor, J. Chem. Phys. 111, 9157 共1999兲. 6 D. Feller and K. A. Peterson, J. Chem. Phys. 108, 154 共1998兲. 7 D. Feller, J. Chem. Phys. 111, 4373 共1999兲. 8 C. Hampel and H. J. Werner, J. Chem. Phys. 104, 6286 共1996兲; G. E. Scuseria and P. Ayala, ibid. 111, 8330 共1999兲. 9 L. A. Curtiss, K. Raghavachari, G. W. Trucks, and J. A. Pople, J. Chem. Phys. 94, 7221 共1991兲; L. A. Curtiss, J. E. Carpenter, K. Raghavarchari, and J. A. Pople, ibid. 96, 9030 共1992兲; L. A. Curtiss, K. Raghavachari, and J. A. Pople, ibid. 98, 1293 共1993兲; ibid. 103, 4192 共1995兲; L. A. Curtiss, K. Raghavarchi, P. C. Redfern, V. Rassolov, and J. A. Pople, ibid. 109, 7764 共1998兲. 10 K. A. Peterson and T. H. Dunning, Jr., J. Chem. Phys. 99, 1930 共1993兲. 11 R. A. Kendall, T. H. Dunning, Jr., and R. J. Harrison, J. Chem. Phys. 96, 6796 共1992兲. 12 MOLPRO is a package of ab initio programs written by H.-J. Werner and P. J. Knowles with contributions from J. Almlo¨f, R. D. Amos, A. Berning et al. 13 ACES II is a program product of the Quantum Theory Project, University of Florida. Authors: J. F. Stanton, J. Gauss, J. D. Watts et al., Integral packages included are VMOL 共J. Almlo¨f and P. R. Taylor兲; VPROPS 共P. Taylor兲; ABACUS 共T. Helgaker, H. J. Aa. Jensen, P. Jørgensen, J. Olsen, and P. R. Taylor兲. 14 P. J. Knowles, C. Hampel, and H.-J. Werner, J. Chem. Phys. 99, 5219 共1994兲. 15 W. J. Lauderdale, J. F. Stanton, J. Gauss, J. D. Watts, and R. J. Bartlett, Chem. Phys. Lett. 187, 21 共1991兲. 16 J. L. Dunham, Phys. Rev. 41, 713 共1932兲. 17 T. H. Dunning, Jr., J. Chem. Phys. 90, 1007 共1989兲. 18 T. H. Dunning, Jr. and K. A. Peterson, J. Chem. Phys. 108, 4761 共1998兲. 19 M. Feyereisen, G. Fitzgerald, and A. Kormornicki, Chem. Phys. Lett. 208, 359 共1993兲; O. Vahtras, J. Almlo¨f, M. W. Feyereisen, ibid. 213, 514 共1993兲; R. A. Kendall and H. Fru¨chtl, Theor. Chim. Acta 97, 159 共1997兲. 20 G. Hetzer, P. Pulay, and H. J. Werner, Chem. Phys. Lett. 290, 143 共1998兲.

Downloaded 22 Feb 2001 to 192.101.100.146. Redistribution subject to AIP copyright, see http://ojps.aip.org/jcpo/jcpcpyrts.html