Electron Correlation Methodology - American Chemical Society

et al. concluded that the protonation in aniline occurs on the Ν atom (6). 9. 0. ©. Η Η la lb ... between the Ν and carbon protonated forms being...
2 downloads 0 Views 877KB Size
Chapter 11

The Protonation Site of Aniline Revisited: A 'Torture Test' for Electron Correlation Methods Downloaded by UNIV OF CALIFORNIA SAN DIEGO on March 19, 2016 | http://pubs.acs.org Publication Date: March 13, 2007 | doi: 10.1021/bk-2007-0958.ch011

1

1

2

A. Daniel Boese, Jan M. L. Martin , Frank De Proft , and Paul Geerlings 2

1

2

Department of Organic Chemistry, Weizmann Institute of Sceince, IL-76100 Rehovot, Israel (email: [email protected]) Eenheid Algemene Chemie, Faculteit Wetenschappen, Vrije Universiteit Brussel, Pleinlaan 2,1050 Brussels, Belgium (email: [email protected])

The site of protonation of aniline is revisited using many high level ab initio wave function methodologies and several D F T exchange-correlation functional. The results indicate that the prediction of the protonation site and the relative energy of the Ν-, ο-, and p-protonated species depend critically on the method and basis set used. At the highest level of theory (W1c), aniline is predicted to be a nitrogen base with a proton affinity of 211.0 kcal/mol, in very good agreement with experiment, the p- and o-protonated forms lying respectively 1.1 and 6.2 kcal/mol higher in energy. Despite yielding errors of less than 5 kcal/mol in all cases, all density functional methods considered wrongly predict protonation at the para site.

© 2007 American Chemical Society

Wilson and Peterson; Electron Correlation Methodology ACS Symposium Series; American Chemical Society: Washington, DC, 2007.

183

184

Introduction

Downloaded by UNIV OF CALIFORNIA SAN DIEGO on March 19, 2016 | http://pubs.acs.org Publication Date: March 13, 2007 | doi: 10.1021/bk-2007-0958.ch011

It is a well established fact that the preferred site of protonation of aniline in solution is the nitrogen atom, due to the higher stabilization upon solvatation of the N-protonated species with respect to the ring protonated ion (/). In the gas phase however, the situation is less clear, and the preferred site of protonation has been the subject of many experimental and theoretical studies. Based on proton transfer equilibria at 600 Κ and upon correlation of the proton affinities with σ* constants for p-substituents,

Lau and Kebarle concluded that the

preferred site of protonation is on the aromatic ring (2). This was confirmed by mass spectroscopic studies involving the collision-activated dissociation of partially deuterated aniline ions, suggesting preferential ring protonation under chemical ionization conditions (5-5). From the correlation of the proton affinities of a series of substituted anilines with Ν is ionization energies, Pollack et al. concluded that the protonation in aniline occurs on the Ν atom (6). 9

0

©

Η Η

la

lb

le

Karpas et al. used the ion mobility/mass spectroscopy method to show that under atmospheric pressure, protonation of the aniline yields two isomeric ions, a N-protonated and a ring protonated species. Assuming that the N-protonated species is the least mobile of the two, they concluded that protonation occurs on the nitrogen (7). Smith et al. showed that the kinetically favored site of protonation is Ν (8). The structure of the protonated species was studied by Nold and Wesdemiotis using neutralization-reionization mass spectrometry. Fast atom bombardment ionization was shown to yield predominantly the N-protonated isomer, whereas chemical ionization with a variety of reagents yielded the ring protonated species (9). Minimal basis set Hartree-Fock and semi-empirical calculations have shown that aniline is a nitrogen base, the energy difference between the Ν and carbon protonated forms being 1-3 kcal/mol (6,10). Sjoberg, Murray, Brinck and Politzer used the average local ionization energy and showed that the para carbon is the most reactive toward electrophiles (11). The absolute minimum in the molecular electrostatic potential of the molecule however is found near the nitrogen atom as was concluded by Ritchie (12). Based on M P 2 computations on HF optimized structures, Hillebrand et al. concluded that aniline is protonated on the nitrogen atom (73). Roy, De Proft and Geerlings used D F T based reactivity indices, such as the relative

Wilson and Peterson; Electron Correlation Methodology ACS Symposium Series; American Chemical Society: Washington, DC, 2007.

Downloaded by UNIV OF CALIFORNIA SAN DIEGO on March 19, 2016 | http://pubs.acs.org Publication Date: March 13, 2007 | doi: 10.1021/bk-2007-0958.ch011

185 nucleophilicity, to study the site of protonation of aniline and a series of m- and p-substituted anilines (14). They also concluded that the nitrogen atom is the site of protonation. These systems were also revisited by Roy et al. in a study addressing the non-negativity of Fukui functions. It was confirmed by Fuentealba, Perez and Contreras that the electrophilic Fukui function indeed predicted the correct protonation site (15). In contrast, Russo et al. concluded that the orbital Fukui function indices, associated with the hardest orbitals, point to the para carbon as the preferred site of protonation (16). Moreover, they performed proton affinity calculations at the B 3 L Y P , B P , MP4 and G2(MP2) levels of theory. A l l methods indeed predict that the nitrogen and para carbon of the aromatic ring are the thermodynamic most favorable protonation sites. However, the D F T methods and M P 4 point to the /?-carbon as the preferred protonation site, whereas G2(MP2) yields N-protonation: the difference with the p-protonated form is only 0.7 kcal/mol. At the G3(MP2) level, N-protonation is favored by 1.2 kcal/mol over ring protonation; this energy difference increases to 10.1 kcal/mol when applying the I P C M continuum method with an G value of 78.3 (water) (17). Water clusters of protonated substituted anilines generated by an electrospray ion source were investigated using a Fourier transform ion cyclotron resonance mass spectrometer by Lee, Cox, Goddard III and Beauchamp (18). From these studies, it was concluded that the preferred site of protonation in the nitrogen atom. From all this theoretical and experimental work, it can be concluded that the energy difference between the N-protonated form la and the p-protonated species lb is very small. (The o-protonated structure 1c lies somewhat higher in energy.) Moreover, the relative stability of the two species changes when using different calculation levels and can even change when the basis set is altered for a given methodology. In this contribution, the definitive relative stability of the different protonated aniline species is theoretically established and the performance of the different methods is critically investigated.

Methods A l l calculations were carried out on a Compaq GS160 minisupercomputer at Brussels Free University and on a Linux farm consisting of dual-CPU 2GHz and 2.4Ghz Pentium Xeon workstations at the Weizmann Institute of Science. The coupled cluster calculations involved in W l c theory were carried out using M O L P R O 2002.3. (19) (For the open-shell calculations on the constituent atoms, the definition of the open-shell C C S D and CCSD(T) energies in Ref. (20) was employed.) The density functional calculations were carried out using a modified version of Gaussian 98 rev. A l l (2/), as were the Gn theory (22-24) and CBS-n (25) calculations.

Wilson and Peterson; Electron Correlation Methodology ACS Symposium Series; American Chemical Society: Washington, DC, 2007.

Downloaded by UNIV OF CALIFORNIA SAN DIEGO on March 19, 2016 | http://pubs.acs.org Publication Date: March 13, 2007 | doi: 10.1021/bk-2007-0958.ch011

186 Aside from the well-known L D A (local density approximation (26, 27)\ B L Y P (Becke-Lee-Yang-Parr (28, 29)) and B 3 L Y P (29, 30) functional, we considered the more recent B97-1 functional (which is a reparametrization (31) of Becke's 1997 hybrid functional), the B97-2 functional (32) (a variation of B97-1 which includes a kinetic energy density term), and the HCTH-407 (33) functional of Boese and Handy (arguably the best G G A functional in existence at the time of writing). The rationale behind the W l c (Weizmann-1 cheap) approach is extensively discussed elsewhere. (34,35) For the sake of selfcontainedness of the paper, we briefly summarize the steps involved for the specific system discussed here: • the geometry is optimized at the B 3 L Y P level with a Dunning correlation consistent (36) cc-pVTZ basis set; • the zero-point energy was obtained from B3LYP/cc-pVTZ harmonic frequencies scaled by 0.985, and rigid rotor-harmonic oscillator (RRHO) thermal corrections were derived from the B3LYP/cc-pVTZ rotational constants and harmonic frequencies; • CCSD(T)/cc-pVnZ+aug(N) (n=D,T) calculations are carried out using a conventional (disk-based) algorithm, where aug(N) stands for the use of the diffuse function augmented aug-cc-pVnZ basis set on the nitrogen atom; • an integral-direct (37) CCSD/cc-pVQZ+aug(N) is carried out; • the SCF energy is extrapolated to the infinite basis limit using £«, = E + (E - E - )l((L/L-\f-\\ where α = 5 and L = {3,4} for {VTZ,VQZ} respectively; L

L

L

X

• the C C S D valence correlation energy is extrapolated using the same expression, but using α = 3.22; • the (T) contribution is extrapolated in the same manner but using L = {2,3} for { V D Z , V T Z } ; • the inner-shell correlation and scalar relativistic contributions are obtained using the Martin-Sundermann-Fast-Truhlar (MSFT) (38) bondadditivity model The 'ultrafine* pruned (99,590) integration grid was used in all D F T calculations.

Wilson and Peterson; Electron Correlation Methodology ACS Symposium Series; American Chemical Society: Washington, DC, 2007.

187

Table I: Computed total atomization energy at 0 Κ of aniline and proton affinities at 298 Κ for protonation at different sites. Proton affinities for N protonation of aniline. All values are given in kcal/mol. TAEo PA 98 PA PA298 N-prot. para vs. N-prot." para ortho para vs. ortho Expt. (39) 1467.7±0.7 210.9 Wlc 1468.7 211.0 209.9 204.9 1.1 5.1 Gl 1462.7 210.8 209.2 204.3 4.9 1.6 G2 1464.0 211.4 210.0 205.1 1.4 4.9 G2MP2 1463.0 211.5 209.9 205.0 1.6 4.9 G3 1467.0 211.4 209.8 5.0 204.8 1.6 G3B3 1467.3 211.4 209.9 204.8 1.5 5.1 CBS-Q 1468.8 210.5 207.9 204.2 3.6 2.7 CBS-QB3 1466.2 211.4 203.2 208.3 3.1 5.1 SVWN3 1738.2 203.1 204.3 4.7 199.6 -1.2 BLYP 1459.7 209.9 212.5 4.1 208.5 -2.7 B3LYP 1461.2 211.0 212.8 208.3 -1.8 4.5 B97-1 1468.7 212.5 214.0 4.4 209.5 -1.5 B97-2 1472.2 213.7 215.0 210.6 -1.4 4.5 HCTH/407 1465.0 212.6 214.9 210.8 4.1 -2.3 HF/VDZ+ 215.1 219.8 214.1 -4.8 5.8 HF/VTZ+ 215.6 220.7 214.7 5.9 -5.1 HF/VQZ+ 215.4 220.5 214.6 -5.2 5.9 H F limit 215.3 220.5 214.5 -5.2 5.9 CCSD/VDZ+ 213.9 214.3 209.1 -0.5 5.3 CCSD/VTZ+ 213.2 213.5 208.0 -0.4 5.5 CCSD/VQZ+ 212.6 212.9 207.4 -0.2 5.5 C C S D limit 212.3 212.5 207.0 5.5 -0.1 MP2/VTZ+ 210.0 203.8 198.8 6.2 5.0 CCSD/VTZ+ 213.2 213.5 208.0 -0.4 5.5 CCSD(T)/VTZ+ 212.1 211.4 206.3 0.7 5.1 CCSD[T]/VTZ+ 211.9 211.7 206.6 0.3 5.1

Downloaded by UNIV OF CALIFORNIA SAN DIEGO on March 19, 2016 | http://pubs.acs.org Publication Date: March 13, 2007 | doi: 10.1021/bk-2007-0958.ch011

2

a

2 9 8

includes correction for internal rotation of the - NH kcal/mol)

+ 3

group (RT/2, +0.3

Wilson and Peterson; Electron Correlation Methodology ACS Symposium Series; American Chemical Society: Washington, DC, 2007.

188

Downloaded by UNIV OF CALIFORNIA SAN DIEGO on March 19, 2016 | http://pubs.acs.org Publication Date: March 13, 2007 | doi: 10.1021/bk-2007-0958.ch011

Results and Discussion A l l relevant results are gathered in Table I. The relative energies apparently exhibit a rather weak basis set dependence: the equilibrium between the N - and p-protonated forms is driven by electron correlation, while the ortho-para equilibrium is apparently quite well reproduced even at the HartreeFock level. We may safely argue that the W l c results are converged with respect to the basis set. The W l c total atomization energy at 0 Κ of aniline, 1468.7 kcal/mol, is in satisfying agreement with the value obtained from heats of formation in the NIST WebBook (59), 1467.7±0.7 kcal/mol. (Most of the uncertainty derives from the heat of vaporization of graphite.) The various contributions to this result are (in kcal/mol): SCF limit 1144.4, valence C C S D correlation energy limit 359.0, connected triple excitations 31.7, inner shell correlation 7.6, scalar relativistic effects -1.2, atomic spin-orbit coupling -0.5 kcal/mol. Extrapolations account for 0.6, 12.1, and 2.5 kcal/mol, respectively, out of the three first contributions. As the M S F T (38) additivity model predicts only very weak core correlation and scalar relativistic contributions to the proton affinity, we have not attempted their explicit (and very expensive) calculation. The B 3 L Y P B97-n, and HCTH/407 compare quite well to experiment and W l c : it should be recalled that all four fiinctionals were wholly or partially fitted to atomization energies. The B97-1 result in fact is identical to W l c , while the 'pure D F T ' HCTH/407 functional does surprisingly well. The relative energies of the three protonated species are well reproduced by all methods from the Gn family. This can largely be explained by (a) the fact that all these methods involved CCSD(T) or QCISD(T) steps (and apparently triple excitations are quite important here); (b) the relatively rapid basis set convergence noted above, which means that it is not really an issue that the CCSD(T) and QCISD(T) steps are carried out in relatively small basis sets. CBS-QB3 likewise reproduces the relative energetics quite well. As has been pointed out in the past (e.g. concerning the linear-cyclic equilibrium in C a n d C carbon clusters (40)), Hartree-Fock underestimates the resonance stabilization of aromatic relative to non-aromatic systems (in the case at hand, between the N - and /^-protonated isomers) and M P 2 overcorrects. The structures are found to be nearly isoenergetic at the C C S D level; inclusion of connected triple excitations favors the N-protonated ion. The direction of the effect of connected quadruples is somewhat unclear, and a CCSD(TQ) or CCSDT(Q) calculation impossible on systems this size, but the contribution will anyhow be much smaller in absolute magnitude than that of connected triple excitations, particularly for systems like these which are dominated by a single reference determinant. We may therefore infer that at the full CI limit, the N protonated species will be slightly more stable than its /?-protonated counterpart. 6

1 0

Wilson and Peterson; Electron Correlation Methodology ACS Symposium Series; American Chemical Society: Washington, DC, 2007.

189

Downloaded by UNIV OF CALIFORNIA SAN DIEGO on March 19, 2016 | http://pubs.acs.org Publication Date: March 13, 2007 | doi: 10.1021/bk-2007-0958.ch011

On the other hand, the PAs at the para and ortho positions are manifestly more sensitive to the electron correlation treatment than the P A at the nitrogen. From the point of view of the 7] diagnostic (a measure for the importance of nondynamical correlation), aniline (7i=0.0113) and N-anilinium (71=0.0096) are very similar (and basically purely single-reference), while the /?- and oprotonated species exhibit very mild multireference character (0.0149 and 0.0157, respectively). Since protonation at these sites thus involves a noticeable change in 71, the P A is expected to be more sensitive to the correlation method. Let us now turn to the density functional methods. A l l of them correctly predict the para-ortho ordering, but considering that this is the case even for a Hartree-Fock treatment this is a somewhat hollow victory. Without exception, all functionals wrongly predict p-protonation. Arguably, this small energy difference falls within the error margin of any type of calibration for (semi-) 'empirical' DFT functionals. In an attempt at finding the 'root cause' of this failure of DFT, we first estimate the conjugation stabilization of the cations by means of isodesmic reaction energies at the CCSD(T)/cc-pVTZ+aug(N) and B3LYP/ccpVTZ+aug(N) levels. (The same qualitative trends are exhibited by other D F T methods.) The isodesmic reactions +

N-prot +7 C H — C H N H + 3 C H + 3 C H 4

3

3

/?-prot + 7 C H — + C H N H 4

3

2

+ 2

4

2

+2 C H +4 C H 2

4

2

(1)

6

6

(2)

have stabilization energies of 65.9 and 67.0 kcal/mol at the CCSD(T) level of theory, respectively. This result is somewhat counterintuitive, since one might expect the conjugation energy of the (aromatic) N-protonated aniline to be the greater of the two. Still, B 3 L Y P overestimates reaction (1) by 3.2 and (2) by 3.8 kcal/mol, hence favoring the /7-protonated aniline by an additional 0.6 kcal/mol. Nevertheless, this does not fully account for the 2.5 kcal/mol discrepancy between B 3 L Y P and CCSD(T) for the energy difference between N - and /?protonated aniline. A second potential source of error may be overstabilization of the 'tautomeric' form of aniline, 2,5-cyclohexadien-l-imine. At the CCSD(T)/ccpVTZ+aug(N) level, the energy difference is 23.9 kcal/mol. This, incidentally, implies a very high proton affinity for the tautomeric form of 235.3 kcal/mol. However, while B 3 L Y P overestimates the tautomerization energy by 3.5 kcal/mol compared to CCSD(T), it does so even more (4.9 kcal/mol) for the imine proton affinity, which is enough to tip the scales to the wrong protonation site. Considering the proton affinities of some small model compounds at the B3LYP/cc-pVTZ+aug(N) and CCSD(T)/cc-pVTZ+aug(N) levels, we find an overestimate by B 3 L Y P of 1.7 kcal/mol for C H N H , and 3.1 kcal/mol for vinyl imine. The 'Schiff base tautomer' of aniline has two conjugated branches similar to vinyl imine, and exhibits an even larger overestimate. (Similar findings are seen with other exchange-correlation functionals.) 2

Wilson and Peterson; Electron Correlation Methodology ACS Symposium Series; American Chemical Society: Washington, DC, 2007.

190 Overall, we see that imperfect cancellation between two (apparently systematic) errors leads to D F T consistently predicting the wrong protonation site.

Downloaded by UNIV OF CALIFORNIA SAN DIEGO on March 19, 2016 | http://pubs.acs.org Publication Date: March 13, 2007 | doi: 10.1021/bk-2007-0958.ch011

Conclusions In this contribution, the protonation site of aniline has been revisited. From our results, it can be inferred that at the full CI limit, the N-protonated species will be slightly more stable than the /?-protonated one. For this problem, W l c theory was found to be converged in basis set and pointed to an N-protonation of the base, the p- and o-protonated forms lying respectively 1.1 and 6.2 kcal/mol higher in energy. The resonance stabilization of N-protonated relative to the pprotonated species is underestimated by Hartree-Fock theory, whereas M P 2 overestimates it. At the C C S D level both forms are nearly isoenergetic, the inclusion of connected triples favoring protonation on the nitrogen. Furthermore, considering the T diagnostic, aniline and the N-protonated conjugate acid are classified as having essentially single-reference character, while the pprotonated and o-protonated forms exhibit very mild multireference character. Finally, all the D F T methods considered severely overestimate the proton affinity of 2,5-cyclohexadien-l-imine at nitrogen, hence erroneously predicting protonation of aniline at the para site rather than on the nitrogen. x

Acknowledgments J M L M is a member of the Helen and Martin Kimmel Center for Molecular Design. A D B acknowledges a postdoctoral fellowship from the Feinberg Graduate School. Research at the Weizmann Institute was partially supported by the Minerva Foundation, Munich, Germany.

Supporting Information B 3 L Y P / c c - p V T Z geometries and vibrational frequencies, as well as total energies at all levels of theory, are available on the World Wide Web at http://theochem.weizmann.ac.il/web/papers/aniline.html

Literature Cited 1. 2.

McMurry, J. Organic Chemistry, Fifth Edition, Brooks/Cole, 2000. Lau, Y . K . ; Kebarle, P. J. Am. Chem. Soc. 1976, 98, 7452.

Wilson and Peterson; Electron Correlation Methodology ACS Symposium Series; American Chemical Society: Washington, DC, 2007.

191 3. 4. 5.

Downloaded by UNIV OF CALIFORNIA SAN DIEGO on March 19, 2016 | http://pubs.acs.org Publication Date: March 13, 2007 | doi: 10.1021/bk-2007-0958.ch011

6. 7. 8. 9. 10. 11. 12. 13. 14. 15. 16. 17. 18. 19.

20. 21.

Wood, Κ. V . ; Burinsky, D. J.; Cameron, D. ; Cooks, R. G . J. Org. Chem. 1983, 48, 5236. Pachuta, S. J.; Isern-Flecha, I.; Cooks, R. G. Org. Mass. Spectrom. 1986, 27, 1. Maquestiau, Α.; Van Haverbeke, Y . ; Mispreuve, H . ; Flammang, R.; Harris, J. Α.; Howe, I.; Benon, J. H . Org. Mass. Spectrom. 1980, 15, 144. Pollack, S. K.; Devlin III, J. L.; Summerhays, K . D.; Taft, R. W.; Hehre, W. J. J Am. Chem. Soc. 1977, 99, 4583. Karpas, Z.; Berant, Z.; Stimac, R. M. Struct. Chem. 1990, 1, 201. Smith, R. L . ; Chyall, L . J.; Beaseley, B . J.; Kenttamaa, Η. I. J. Am. Chem. Soc. bf 1995, 117, 201. Nold, M. J.; Wesdemiotis, C. J. Mass. Spectrom. 1996, 31, 1169. Dewar, M. J. S.; Dieter, K . M. J. Am. Chem. Soc. 1986, 108, 8075. Sjoberg, P.; Murray, J. S.; Brinck, T.; Politzer, P. Can. J Chem. 1990, 68, 1440. Ritchie, J. P. J. Mol. Struct. (Theochem) 1992, 225, 297. Hillebrand, C.; Klessinger, M.; Eckert-Maksié, M . ; Maksié, Ζ. Β. J. Phys. Chem. 1996, 100, 9698. Roy, R. K.; De Proft, F.; Geerlings, P. J. Phys. Chem. A 1998, 102, 7035. Fuentealba, P.; Perez, P.; Contreras, R. J. Chem. Phys. 2000, 113, 2544. Russo, N.; Toscano, M.; Grand, Α.; Mineva, T. J. Phys. Chem. A 2000, 104, 4017. Bagno, Α.; Terrier, F. J. Phys. Chem. A 2001, 105, 6537. Lee, S.-W.; Cox, H . ; Goddard III, W. Α.; Beauchamp, J. L . J. Am. Chem. Soc. 2000, 122, 9201. M O L P R O is a package of ab initio programs written by H.-J. Werner and P. J. Knowles, with contributions from J. Almlöf, R. D. Amos, A . Berning, D. L . Cooper, M. J. O. Deegan, A . J. Dobbyn, F. Eckert, S. T. Elbert, C. Hampel, R. Lindh, A . W. Lloyd, W. Meyer, A . Nicklass, K . Peterson, R. Pitzer, A . J. Stone, P. R. Taylor, M. E. Mura, P. Pulay, M. Schütz, H . Stoll, and T. Thorsteinsson. Watts, J. D.; Gauss, J.; Bartlett, R. J.; J. Chem. Phys. 1993, 98, 8718. Frisch, M . J.; Trucks, G . W.; Schlegel, H . B . ; Scuseria, G . E.; Robb, Μ. Α.; Cheeseman, J. R.; Zakrzewski, V . G.; Montgomery, Jr., J. Α.; Stratmann, R. E.; Burant, J. C.; Dapprich, S.; Millam, J. M.; Daniels, A . D.; Kudin, Κ. N.; Strain, M. C.; Farkas, O.; Tomasi, J.; Barone, V.; Cossi, M.; Cammi, R.; Mennucci, B.; Pomelli, C.; Adamo, C.; Clifford, S.; Ochterski, J.; Petersson, G. Α.; Ayala, P. Y . ; Cui, Q.; Morokuma, K . ; Salvador, P.; Dannenberg, J. J.; Malick, D. K . ; Rabuck, A . D.; Raghavachari, K . ; Foresman, J. B.; Cioslowski, J.; Ortiz, J. V.; Baboul, A . G.; Stefanov, Β. B . ; Liu, G . ; Liashenko, Α.; Piskorz, P.; Komaromi, I.; Gomperts, R.; Martin, R. L.; Fox, D. J.; Keith, T.; Al-Laham, Μ. Α.; Peng, C. Y . ; Nanayakkara, Α.; Challacombe, M.; Gill, P. M. W.; Johnson, B . G.; Chen, W.; Wong, M. W.;

Wilson and Peterson; Electron Correlation Methodology ACS Symposium Series; American Chemical Society: Washington, DC, 2007.

192

22. 23. 24.

Downloaded by UNIV OF CALIFORNIA SAN DIEGO on March 19, 2016 | http://pubs.acs.org Publication Date: March 13, 2007 | doi: 10.1021/bk-2007-0958.ch011

25. 26. 27. 28. 29. 30. 31. 32. 33. 34. 35.

36.

37. 38. 39.

40.

Andres, J. L.; Gonzalez, C.; Head-Gordon, M.; Replogle, F. S.; Pople, J. Α.; Gaussian 98, Revision A.]] (Gaussian, Inc., Pittsburgh, PA, 1998). Curtiss, L . Α.; Raghavachari, K . ; Trucks, G . W.; Pople, J. Α.; J. Chem. Phys. 1991, 94, 7221. Curtiss, L . Α.; Raghavachari, K . ; Redfern, P. C.; Rassolov, V . ; Pople, J. Α.; J. Chem. Phys. 1998, 109, 7764. Curtiss, L . Α.; Redfern, P. C.; Raghavachari, K . ; Pople, J. Α.; J. Chem. Phys. 2001, 114, 108. Montgomery, Jr., J. Α.; Frisch, M. J.; Ochterski, J. W.; Petersson, G . Α.; J. Chem. Phys. 1999, 110, 2822 and references therein. Slater, J. C.; Phys. Rev. 1951, 81, 385. Vosko, S.H.; Wilk, L.; Nusair, M.; Can J. Phys. 1980, 58, 1200. A . D. Becke, Phys. Rev. A 1988, 38, 3098. Lee, C.; Yang, W.; Parr, R. G.; Phys. Rev. Β 1988, 37, 785. Becke, A . D.; J. Chem. Phys. 1993, 98, 5648. Hamprecht, F. Α.; Cohen, A . J.; Tozer, D. J.; Handy, N. C.; J. Chem. Phys. 1998, 109, 6264. Wilson, P. J.; Bradley, T. J.; Tozer, D . J.; J. Chem. Phys. 2001, 115, 9233. Boese, A . D.; Handy, N. C.; J. Chem. Phys. 2001, 114, 5497. Martin, J. M. L.; de Oliveira, G . J. Chem. Phys. 1999, 111, 1843. Martin, J. M . L . ; Parthiban, S. in: Quantum-mechanical prediction of thermochemical data (ed. J. Cioslowski), Understanding Chemical Reactivity, vol. 22 (Kluwer Academic Publishers, Dordrecht, 2001), pp. 3 1 65. Dunning Jr., T . H . ; Peterson, K . A . ; Woon, D.E.; "Correlation consistent basis sets for molecular calculations", in Encyclopedia of Computational Chemistry, ed. P. von Ragué Schleyer, (Wiley & Sons, Chichester, 1998), pp. 88-115. Schütz, M.; Lindh, R.; Werner, H.-J.; Mol. Phys. 1999, 96, 719 Martin, J. M. L . ; Sundermann, Α.; Fast, P. L . ; Truhlar, D. G.; J. Chem. Phys. 2000, 113, 1348 Afeefy, H.Y.; Liebman, J.F.; Stein, S.F.; "Neutral Thermochemical Data", and Hunter, E.P.; Lias, S. G.;, "Proton Affinity Data", in NIST Chemistry WebBook, NIST Standard Reference Database Number 69, Mallard, W.G.; Linstrom, P.G.; Eds.; February 2000, National Institute of Standards and Technology, Gaithersburg, M D , 20899; http://webbook.nist.gov. Martin, J. M. L.; Taylor, P. R. J. Phys. Chem. 1996, 100, 6047

Wilson and Peterson; Electron Correlation Methodology ACS Symposium Series; American Chemical Society: Washington, DC, 2007.