78
J. Phys. Chem. A, Vol. 103, No. 1, 1999
McClean
9
1 3
3
d 4s a D3 state, any barriers due to the closed s-subshell could
oxygen atoms, are observed and calculated. The binding energy
of NiSO2 estimated in this work is approximately equal to the
be screened. On the basis of the orbital-occupancy argument
presented above, we conclude that the ground state is converted
to the a D3 state, which is reactive. k∞ is close to the gas kinetic
2
binding energy of C2V η O,O Ni-O2 from an ab initio study: Eb
3
-1 22
) 48 ( 7 kcal mol . This observation suggests that the
rate constant, and ko decreases with increasing temperature,
indications of no energy barrier to product formation.
bonding in C2V NiO2 and C2V NiSO2 is similar. Covalent and
ionic components are present in NiO2. The relatively large
limiting low-pressure and high-pressure rate constants for Ni
Physical quenching is the likely depletion channel of the states
6
-
1
+ SO2 and the larger electron affinity of SO2 (than O2) suggest
above 880 cm . An abstraction channel is not accessible, and
-
10
3
-1
an ionic component in NiSO2.
rate constants on the order of 10
cm s for a termolecular
The reaction of Mn(3d 4s a6S5/2) with SO2 was also
5
2
channel at only 10 Torr total pressure is not expected. E-V
energy transfer may take place via the vibrational modes of SO2
investigated in this work to determine if a relatively unreactive
-
14
3
TM atom would react with SO2. An upper limit of ∼10 cm
since the energy differences between the Ni states are close to
-
1
-1
the vibrational modes of SO2.7
molecule
s
was obtained over the temperature range of
2
96-622 K. The electron-transfer mechanism has been offered
To estimate the binding energy Eb of the NiSO2 adduct
as an explanation for the reactions of TM atoms with SO2. The
ionization energy of Mn is less than that of Ni. Thus, the orbital
3
relative to Ni(a D3), simplified RRKM calculations using the
6
9
formalism of Troe were performed at 296 K. In brief, the rate
occupancy of the TM atom is the predominant driving force
constant for the unimolecular dissociation in the low-pressure
limit, kuni, of NiSO2 was first calculated. ko was then calculated
from the equilibrium expression Keq ) kuni/ko, where Keq is the
equilibrium constant for NiSO2 h Ni + SO2. The collision
efficiency, âc, of Ar was assumed to be 0.20. âc can be
determined by comparing the experimental and calculated
(at least when comparing Mn and Ni).
Summary and Conclusions
Results presented here indicate that Ni is very reactive toward
3
3
3
SO2. The three lowest states of Ni, a F4, a D3, and a D2, react
with SO2 via a termolecular mechanism. The 1:1 adduct
formation is efficient. The kinetic results were combined with
DFT and RRKM calculations to provide an estimate of 47 kcal
(strong collision) third-order rate constants. However, the
uncertainties and adjustable binding energy in the calculations
precluded such an approach. The molecular structure and
vibrational frequencies of NiSO2 needed for the RRKM calcula-
-
1
2
mol for the binding energy of NiSO2 (η O,O isomer). A
comparison between NiO and NiSO suggest covalent and ionic
8
tions were calculated from density functional methods. The
2
2
3
3
3
components to the bonding in NiSO2. The a F3, a F2, and a D1
states of Ni deplete at approximately the gas collision rate in
the presence of SO2.
binding energy was varied until agreement was obtained
between the calculated and experimental third-order rate con-
stant. The density functional calculations are not intended to
be exhaustive, but instead serve mainly to provide estimates of
the molecular parameters.
Acknowledgment. This research was supported by the Naval
Academy Research Council and a Cottrell College Science
Award of Research Corporation.
The molecular parameters of NiSO2 were calculated by the
15,16
SPARTAN
suite of programs. Calculations were performed
with the LSDA/pBP86 model and DN** numerical basis sets.
Several triplet isomers were assumed. The “side-on bonded,”
η O,O isomer (C2V geometry) gave the lowest energy. The DFT
results are to be taken with caution; numerical errors are likely
since calculations involve numerical integration steps. Addition-
ally, the relatively large number of low-lying states of Ni and
other TM atoms poses a challenge for the theoretical treatment
of transition metals. Pertinent output from the computations are
listed in the Appendix.
Appendix
2
The DFT suite of SPARTAN programs15,16 was used to
calculate the molecular parameters of NiSO2. For the triplet η O,O
2
isomer of NiSO2 (C2V geometry), we found bond lengths of r(S-
O) ) 1.5847 Å and r(O-Ni) ) 1.9539 Å. Bond angles are
∠ (O-S-O) ) 100.1412°, ∠ (S-O-Ni) ) 91.4730°, and ∠ (O-
Ni-O) ) 76.9128°. Calculated frequencies, in units of cm ,
are 127.12, 301.61, 349.38, 504.62, 836.33, and 841.61.
For the simplified RRKM calculations at 296 K, Eo was taken
as an adjustable parameter. The rate constant for the unimo-
lecular dissociation in the low-pressure limit, kuni, is given by 9
-
1
Agreement between the measured and calculated ko at 296
-
1
K was obtained for Eb ) 47 ( 3 kcal mol ; the uncertainty in
Eb is based only on the overall uncertainty in ko at 296 K. Since
there is no activation barrier, significantly smaller Eb values
would give association rate constants that are too small because
ko is a sensitive function of the binding energy.13 The presence
of several low-lying states in atomic Ni and the possibility of
several low-lying bound states in NiSO2 could lead to additional
pathways, thus affecting the calculated results. The equilibrium
constant, and hence ko, is also affected by the relatively low
level of the DFT calculations. Thus, only an estimate of Eb is
reported. This estimate is reasonable, however. A comparison
of the rates of TM and main-group metals with small oxygen-
containing molecules indicates that a bond energy of 40-50
F(E )RT
Qvib
o
kuni ) âCZLJ
exp(-E /RT)F F F FrotintnFcorr
o anh E rot
where âC is the collisional efficiency of Ar, ZLJ is the Leonard-
Jones rate constant, F(Eo) is the vibrational density of states of
NiSO2 at the threshold energy Eo for dissociation, Qvib is the
vibrational partition function of NiSO2, Fanh is a correction for
vibrational anharmonicity, FE is a correction for the variation
of the density of states, and Frot is the molecular rotational
correction factor. The corrections for internal rotational modes
and the coupling of the F factors, Frot int and Fcorr, were taken
as one.
-
1
kcal mol is in accord with third-order rate constants of the
magnitude measured in this work.1
7-21
-
10
cm3
The binding energies of Na-SO2 and K-SO2 have been
A Lennard-Jones collision frequency of 6.38 × 10
-
1
-1
found to lie within the range of the corresponding binding
molecule was calculated from reasonable values of the
s
energies of the metal-O2 superoxides.1
9,20
C2V η O,O-type
2
Lennard-Jones parameters. âc was taken as 0.20. s ) 6 and m
) 3, where s is the number of vibrational modes in NiSO2 and
structures, where the alkali-metal atoms are side-bonded to the