C O M M U N I C A T I O N S
predicted by SVWN/6-311++G** may be more accurate than the
CCSD(T) and the other DFT distances, which may overestimated.
The overly long (2.296 Å) separation predicted by CASSCF(11,11)/
6-311G** evidently is due to the lack of dynamic electron
correlation. These problems - self-interaction error in DFT, large
spin contamination in CCSD(T), and lack of dynamic correlation
in CASSCF - might, in principle, be overcome at multireference
levels, such as MRSDCI and MRCI-SD/AQCC, but optimization
and frequency computations are not feasible at present.
In summary, we have characterized H2OOO+ to be a complex
involving a 3c-1e bond, instead of a simple ion-molecule
complex. The H2OOO+ cation is stable under visible irradiation; it
only was destroyed by UV photolysis. However, the cation
absorption diminished quickly on annealing, suggesting high
reactivity toward, for example, O2 or H2O.
Figure 2. (A) Geometries calculated at various levels (bond length in Å,
+
bond angle in deg). (B) Interaction between the πy* of O2 and in-plane
water lone pair. The values in parentheses are S2 values.
Table 2. Comparison of Calculated and Observed IR Frequencies
Acknowledgment. The NSFC (20125033), the NKBRSF, and
the University of Georgia supported this work. We thank Hans
Lischka for attempting multireference calculations, and H. F.
Schaefer, III and the referees for helpful discussions.
(in cm-1) for H2OOO+ (Values in Parentheses Are IR Intensities)
OH asy str.
OH sym. str.
O−O str.
H O−O str.
2 2
B3LYP
BHHLYP
BLYP
3723(238)
3861(320)
3587(206)
3615(211)
3624(271)
3836
3621(535)
3746(884)
3490(340)
3516(323)
3521(349)
3738
1762(593)
1860(1423)
1645(240)
1694(219)
1773(186)
1820
679(133)
687(109)
689(167)
711(169)
745(190)
647
Supporting Information Available: Spectra and calculated vibra-
tional frequencies and intensities (PDF). This material is available free
BP86
SVWN
CCSD(T)
observed
3430.7
3341.7
1601.1a
731.0
References
a The D2OOO value. The frequency for H2OOO is expected to be 2-3
(1) Byers Brown, W. Chem. Phys. Lett. 1995, 235, 94.
(2) Cacace, F.; de Petris, G.; Pepi, F.; Troiani, A. Angew. Chem., Int. Ed.
2000, 39, 367.
(3) Wentworth, P., Jr.; Jones, L. H.; Wentworth, A. D.; Zhu, X.; Larsen, N.
A.; Wilson, I. A.; Xu, X.; Goddard, W. A.; Janda, K. D.; Eschenmoser,
A.; Lerner, R. A. Science 2001, 293, 1806.
(4) Engdahl, A.; Nelander, B. Science 2002, 295, 482.
(5) Fehsenfeld, F. C.; Mosesman, M.; Ferguson, E. E. J. Chem. Phys. 1971,
55, 2115.
(6) Howard, C. J.; Bierbaum, V. M.; Rundle, H. W.; Kaufman, F. J. Chem.
Phys. 1972, 57, 3491.
(7) Smith, G. P.; Lee, L. C. J. Chem. Phys. 1978, 69, 5393.
(8) Beyer, R. A.; Vanderhoff, J. A. J. Chem. Phys. 1976, 65, 2313.
(9) McCrumb, J. L. Planet. Space Sci. 1982, 30, 559.
(10) Angel, L.; Stace, A. J. J. Phys. Chem. A 1999, 103, 2999.
(11) Cheng, B. M.; Lee, Y. P.; Ogilvie, J. F. Chem. Phys. Lett. 1988, 151,
109.
(12) Milligan, D. E.; Jacox, M. E. J. Chem. Phys. 1963, 38, 2627.
(13) Smith, D. W.; Andrews, L. J. Chem. Phys. 1974, 60, 81.
(14) (a) Langford, V. S.; McKinley, A. J.; Quickenden, T. I. J. Am. Chem.
Soc. 2000, 122, 12859. (b) Cooper, P. D.; Kjaergaard, H. G.; Langford,
V. S.; McKinley, A. J.; Quickenden, T. I.; Schofield, D. P. J. Am. Chem.
Soc. 2003, 125, 6048.
cm-1 higher.
nature of the bonding.19 The optimized geometries at various levels
of theory are given in Figure 2A, and the computed harmonic
frequencies are compared with the experimental values in Table 2.
At all levels employed, the radical has a doublet ground state (2A′′)
structure with Cs symmetry which correlates with H2O + O2+. As
+
shown in Figure 2B, the empty πy* orbital of O2 (the singly
occupied πy* in the ground state of O2+) interacts with an in-plane
water lone pair orbital favorably. The single electron occupying
this orbital glues O2+ and H2O together. The Wiberg bond indices
of O1O2 and O1O3 are 0.26 and 0.21, and the O2 and H2O fragment
charges are +0.59 and +0.41e, respectively. Therefore, this radical
is a complex involving a 3c-1e bond, rather than a simple ion-
molecule complex, O2+‚H2O. The calculated isotopic frequency
ratios (Table 2 of Supporting Information) are in good agreement
with the experimental values, which add strong support to the
experimental assignment.
(15) Andrews, L.; Ault, B. S.; Grzybowski, J. M.; Allen, R. O. J. Chem. Phys.
1975, 62, 2461.
(16) Chertihin, G. V.; Andrrews, L. J. Chem. Phys. 1998, 108, 6404.
(17) (a)Thompson, W. E.; Jacox, M. E. J. Chem. Phys. 1989, 91, 3826. (b)
Chertihin, G. V.; Saffel, W.; Yustein, J. T.; Andrews, L.; Neurock, M.
Ricca, A.; Bauschlicher, C. W. J. Phys. Chem. 1996, 100, 5261.
(18) Zhou, M. F.; Andrews, L. J. Am. Chem. Soc. 1998, 120, 11499.
(19) Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M.
A.; Cheeseman, J. R.; Zakrzewski, V. G.; Montgomery, J. A., Jr.;
Stratmann, R. E.; Burant, J. C.; Dapprich, S.; Millam, J. M.; Daniels, A.
D.; Kudin, K. 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. A.; Ayala, P. Y.; Cui, Q.; Morokuma, K.;
Malick, D. K.; Rabuck, A. D.; Raghavachari, K.; Foresman, J. B.;
Cioslowski, J.; Ortiz, J. V.; Baboul, A. G.; Stefanov, B. B.; Liu, G.;
Liashenko, A.; Piskorz, P.; Komaromi, I.; Gomperts, R.; Martin, R. L.;
Fox, D. J.; Keith, T.; Al-Laham, M. A.; Peng, C. Y.; Nanayakkara, A.;
Gonzalez, C.; Challacombe, M.; Gill, P. M. W.; Johnson, B.; Chen, W.;
Wong, M. W.; Andres, J. L.; Gonzalez, C.; Head-Gordon, M.; Replogle,
E. S.; Pople, J. A. Gaussian 98, revision A.7; Gaussian, Inc.: Pittsburgh,
PA, 1998.
The O1-O2 stretching vibrational frequency is an important
characteristic of the interaction between the two fragments. Mainly
due to the neglect of anharmonicity, computed vibrational frequen-
cies are generally higher than the experimental values. While the
OH symmetric stretching, OH antisymmetric stretching, and O2-
O3 stretching frequencies follow this generalization (Table 2), the
O1-O2 stretching frequency is underestimated. Thus, the computed
hybrid B3LYP (679 cm-1) and BHHLYP (687 cm-1) values are
smaller than the experimental value of 731 cm-1. The gradient-
corrected BLYP and BP86 values of 689 and 711 cm-1, respec-
tively, are closer. The CCSD(T) “panacea” did not improve the
agreement. We attribute this to the large spin contamination: the
S2 ) 1.68 value for the doublet is quite large. Surprisingly, the
simpler SVWN/6-311++G** method, which generally is less
accurate than hybrid and gradient-corrected DFT, performs well
and gives a 745 cm-1 frequency. As has been discussed in several
papers,20,21 DFT methods suffer from self-interaction errors in
describing radical cations with odd electron bonds (e.g., hemibonds).
However, the simpler LSD method is known to perform better for
hemi-bonded species.21 Therefore, the O1O2 bond length (1.965 Å)
(20) (a) Braida, B.; Lauvergnat, D.; Hiberty, P. C. J. Chem. Phys. 2001, 115,
90. (b) Gruning, M.; Gritsenko, O. V.; van Gisbergen, S. J. A.; Baersends,
E. J. J. Phys. Chem. A 2001, 105, 9211. (c) Ghanty, T. K.; Ghosh, S. K.
J. Phys. Chem. A 2002, 106, 11815. (d) Sodupe, M.; Bertran, J.; Bodriguez-
Santiago, L.; Baerends, E. J. J. Phys. Chem. A 1999, 103, 166. (e) Xie,
Y. M.; Schaefer, H. F.; Fu, X. Y.; Liu, R. Z. J. Chem. Phys. 1999, 111,
2532. (f) Jaramillo, J.; Scuseria, G. E.; Ernzerhof, M. J. Chem. Soc. 2003,
118, 1068.
(21) Braida, B.; Hiberty, P. C.; Savin, A. J. Phys. Chem. A 1998, 102, 7872.
JA037125Z
9
J. AM. CHEM. SOC. VOL. 125, NO. 38, 2003 11513