High Electron Affinities of SO4 and HSO4
J. Phys. Chem. A, Vol. 104, No. 3, 2000 507
TABLE 4: Relative Energies (kcal/mol) and Number of
Imaginary Frequencies of the Various Symmetries of SO
at B3LYP/TZVP+, B3LYP/6-31+G*, and MP2/6-31+G*
TABLE 5: Energies (kcal/mol) of Various Symmetries of
-
-
4
SO
4
Relative to C2W from CCSD(T)/6-311+G*//B3LYP/
tzvp+ Calculations
B3LYP/tzvp+
B3LYP/6-31+G* MP2/6-31+G*
symmetry MP2
MP4(sdq)
CCSD
CCSD(T)
symmetry energy imaginary energy imaginary energy imaginary
D
2d 2A
2
-3.94
-3.52
0.0
8.39
3.75
0.0
9.68
4.56
0.0
5.88
2.85
0.0
2
C
C
C
3V
B
A
1
D
2d 2A
2
1.98
1.90
0.0
2
2
0
1
3.63
1.90
0.0
2
2
0
1
0.0
0.54
3.08
12.67
0
0
0
1
2
2
2V
2
C
C
C
3V
B
A
1
2A
12.59
6.20
2.80
2.59
2
s
”
2V
2
2A
1.65
1.64
s
”
TABLE 6: Relative Energies (kcal/mol) of the Various
Symmetries of SO at B3LYP/6-31+G* and MP2/6-31+G*
ADE and VDE values of SO4- were in good agreement with
the experimental values. However, the difference between the
calculated ADE and VDE was much larger than that which was
experimentally determined, at the highest level of theory carried
out, suggesting that the measured ADE may only represent an
upper limit of the real adiabatic electron affinity. Theory predicts
4
symmetry B3LYP/tzvp+ B3LYP/6-31+G* MP2/6-31+G*
D2d 1A1
28.20
33.78
0.0
27.11
32.71
0.0
0.0
12.47
17.53
1
C
C
3V
A
A
1
1
2V
2
TABLE 7: Relative Energies (kcal/mol) of the Various
Symmetries of SO
Calculations
-
that both HSO4 and HSO4 have Cs symmetry. The calculated
4
from CCSD(T)/6-311+G*//B3LYP/tzvp+
-
VDE and ADE for HSO4 are in excellent agreement with the
experimental values.
symmetry
MP2
MP4(sdq)
CCSD
CCSD(T)
D
2d 1A
1
-10.12
5.18
43.92
35.19
0.0
36.87
34.07
0.0
15.68
17.89
0.0
Acknowledgment. This work is supported by The U.S.
Department of Energy, Office of Basic Energy Sciences,
Chemical Science Division. Acknowledgment is also made to
the donors of the Petroleum Research Fund, administered by
the American Chemical Society, for partial support of this
research. The work was performed at the W. R. Wiley
Environmental Molecular Sciences Laboratory, a national
scientific user facility sponsored by DOE’s Office of Biological
and Environmental Research and located at Pacific Northwest
National Laboratory, which is operated for DOE by Battelle.
L.S.W. is an Alfred P. Sloan Foundation Research Fellow.
1
C
C
3V
A
A
1
1
2V
1
0.0
three levels of theory were in close agreement with each other,
despite the differing numbers of negative frequencies, we used
only our B3LYP/tzvp+ geometries. The relative CCSD(T)/6-
11+G* energies, including the MP2, MP4(sdq), and CCSD
predictions, also obtained from these single-point calculations,
are shown in Table 5.
Only the MP2/6-311+G* calculation predicts that the D2d
symmetry is lowest in energy. All of the more extensive
treatments of electron correlation predict the lowest energy for
the C2V geometry, which is consistent with the DFT results. The
difference in energy between the various states becomes less
when the triple excitations are included, but CCSD(T) still gives
larger relative energy differences than those given by B3LYP.
From these results, we tentatively concluded that the lowest-
3
Appendix
4
In a recent paper, McKee presented a theoretical study of a
-
variety of sulfur-containing complexes, including SO4 and SO4 .
B3LYP/6-31+G* and MP2/6-31+G* optimizations were re-
-
ported for SO4 and SO4 , with single-point energies calculated
-
at higher levels of theory. Frequencies were also obtained at
the B3LYP/6-31+G* level. The zero-point energies from these
frequency calculations were given, but no other details were
presented. McKee reported one geometry for SO4 (C2V) and three
energy geometry of SO4 has C2V symmetry. We still cannot
reach a firm conclusion, because the energy differences among
the various isomers of SO4 are not great enough (Table 5). In
-
principle, they might even be accessible at our experimental
conditions (room temperature). Frequency calculations at a
higher-order level of correlation (i.e., MP4) would be needed
to verify whether the D2d and C3V geometries are stable points.
Such studies are beyond the scope of the current investigation.
Extending this approach to SO4, we find that B3LYP again
predicts that the lowest energy isomer is that of C2V symmetry,
whereas MP2 predicts that the D2d isomer is the lowest in
-
for SO4 (D2d, C3V, and C2V). We found that the relative energies
-
of the three isomers of SO4 were different at the B3LYP/6-
3
1+G* level from those found at the MP2/6-31+G* level. To
-
obtain the ADE and VDE for SO4 , for comparison to the
experiment, we had to determine which of these three states to
choose as the anion ground state. Our optimizations and
frequency calculations at the B3LYP/tzvp+ level indicated that
the D2d and C3V geometries were not stable points. We thus
decided to duplicate the B3LYP/6-31+G* and MP2/6-31+G*
optimizations and frequency calculations to clarify the situation.
The relative energies and number of imaginary frequencies for
the various states are presented in Table 4. We were also able
to optimize a geometry of Cs symmetry. At the B3LYP/tzvp+
and B3LYP/6-31+G* levels; only the C2V geometry is a
minimum. The C2V molecule is also the lowest in energy. The
other isomers have one or two large imaginary frequencies.
However, at the MP2/6-31+G* level, only the Cs isomer has a
negative frequency. In this case, the D2d geometry is the lowest
in energy. The B3LYP/tzvp+ structures of the four symmetries
-
energy. However, contrary to the results for SO4 , both DFT
and MP2 predict that all the symmetries are stable points. In
addition, during attempts to optimize a Cs geometry, the
molecule reverted to C2V symmetry at both levels of theory.
The relative energies of the various isomers are given in Table
6, and the B3LYP/tzvp+ structures are shown in Figure 3b.
To again clarify the correct symmetry, we obtained CCSD-
(T)/6-311+G* single-point energies for each species. As before,
we used the B3LYP/tzvp+ geometries for the single point
calculations. The calculated energies, relative to those of the
C2V isomer, are shown in Table 7.
Once again, the MP2/6-311+G* single-point energy is the
lowest for the D2d geometry, whereas all the other results
indicate that the C2V geometry is the energy minimum. The
relative energies from the more extensive treatment of electron
correlation are not close to those obtained with DFT; CCSD-
(T) gives predictions of the relative energies that are ∼10 kcal/
-
of SO4 are shown in Figure 3a.
The only isomer that is a minimum at both levels of theory
is that of C2V symmetry. To verify which symmetry was the
lowest energy state, we obtained CCSD(T)/6-311+G* single-
point energies for each species. Although the geometries at all