1510 J. Phys. Chem. A, Vol. 104, No. 7, 2000
Good et al.
Figure 5. Relative energetics for methyl formate reactions.
CH3OCOH f CH3OCO + H
reaction 1a resulting in the formation of CH3OCO radical. A
reasonable estimate for the uncertainty associated with the
computation of each activation energy is (0.5 kcal mol-1. With
this uncertainty, the ratio k1a/k1 could range from 65% to 97%.
Bartels et al.23 examined the analogous hydrocarbon system,
i.e.,
(2a)
(2b)
CH3OCOH f CH2OCOH +H
The 298 K heat of formation9 of the hydrogen atom is 52.1
kcal mol-1. Thus, the BDE for reactions 2a and 2b is estimated
to be 99.6 and 100.6 kcal mol-1, respectively. Previous
investigators have noted a correlation between a molecule’s
C-H bond dissociation energy and the reaction rate of CH bond
dissociation reactions.27 In addition to being thermodynamically
favored, our analysis suggests that removal of the carbonyl
hydrogen is also kinetically favored over removal of one of the
methyl hydrogen atoms. This suggestion is supported by our
molecular orbital calculations of the activation energies for
pathways 1 and 2. At the QCISD(T)/6-311++G(2d,2p)//UMP2/
6-311++G(2d,2p) level of theory, removal of the carbonyl
hydrogen is calculated to proceed over a reaction barrier of 0.2
kcal mol-1 while the activation energy for removal of a methyl
CH3CHO + Cl f CH3CO + HCl
CH3CHO + Cl f CH2CHO + HCl
(4a)
(4b)
Their investigation determined k4b/k4a to be 0.0753 and thus k4a/
k4 to be over 92%. Thus, for the acetylaldehyde system, removal
of the carbonyl hydrogen is shown to be favored over abstraction
of a methyl hydrogen. The analogous conclusions are suggested
for the methyl formate system.
IV. Conclusions
hydrogen is a substantially higher 1.9 kcal mol-1
.
FTIR analysis of the products from the reaction of chlorine
atoms with methyl formate show that both formic acid anhydride
and formic acid are products of the reaction. Chlorine radical
initiated hydrogen abstraction of methyl formate is found to
2. Comparison of Branching Ratio and Rate Constant from
Theory and Experiment. For reactions 1a and 1b, the branching
ratio is expressed as the ratio of the rate constants, i.e., k1a/k1.
From transition state theory, each rate constant is given by the
following expression:
occur with a rate of 1.4 ( 0.5 × 10-12 cm3 molecule-1 s-1
.
This value is in agreement with previous determinations by
Wallington et al.16 (1.4 ( 0.1 × 10-12 cm3 molecule-1 s-1
)
Qq
LkbT
e-E /k T
(3)
and Notorio et al.17 (1.8 ( 0.2 × 10-12 cm3 molecule-1 s-1).
Ab initio calculations at the QCISD(T)/6-311++G(2d,2p)//
MP2/6-311++G(2d,2p) level of theory predict the rate to be
2.8 × 10-12 cm3 molecule-1 s-1, which is in reasonable
agreement with experimental determinations. Additionally, our
calculations suggest that the reaction between methyl formate
and chlorine occurs predominately (90%) at the carbonyl
hydrogen resulting in the formation of CH3OCO radical.
a
b
k )
h
QClQmethylformate
where Ea is the activation energy, T is the temperature, Q
represents the total partition function incorporating translational,
rotational, vibrational, and electronic terms (QT ) QeQvQrQt),
L is a statistical factor representing the number of equivalent
extractable hydrogen atoms, and kb is Boltzmann’s constant.
Tunneling corrections, which may slightly lower the activation
energy, have not been included in this work. Removal of a
methyl hydrogen is estimated to have a rate constant of 2.8 ×
10-13 cm3 molecule-1 s-1, while removal of the carbonyl
hydrogen is predicted to have a rate constant of 2.5 × 10-12
cm3 molecule-1 s-1 at 298 K. The combined rate constant is
thus 2.8 × 10-12 cm3 molecule-1 s-1. The experimentally
determined rate16,17 has previously been found to range between
1.4 × 10-12 and 1.8 × 10-12 cm3 molecule-1 s-1. Thus, the
298 K rate constant for the reaction of methyl formate with
chlorine atom as determined by our ab initio methodology is a
reasonable estimate. It is also reasonable to expect the branching
ratio of reactions 1a and 1b to be more reliable, as errors
common to both calculations may cancel.
References and Notes
(1) Rouhi, A. M. Chem Eng. News 1985, 44, 37.
(2) Japar, S. M.; Wallington, T. J.; Richert, J. F. O.; Ball, J. C. Int. J.
Chem. Kinet. 1990, 22, 1257.
(3) Jenkin, M. E.; Hayman, G. D.; Wallington, T. J.; Hurley, M. D.;
Ball, J. C.; Nielsen, O. J.; Ellerman, T. J. Phys. Chem. 1993, 97, 11712.
(4) Wallington, T. J.; Hurley, M. D.; Ball, J. C.; Jenkin, M. E. Chem.
Phys. Lett. 1993, 211, 41.
(5) Langer, S.; Ljungstrom, E.; Ellerman, T.; Nielsen, O. J.; Sehested,
J. Chem. Phys. Lett. 1995, 240, 53.
(6) Sehested, J.; Mogelberg, T.; Wallington, T. J.; Kaiser, E. W.;
Nielsen, O. J. J. Phys. Chem. 1996, 100, 17218.
(7) Sehested, J.; Sehested, K.; Platz, J.; Egsgaard, H.; Nielsen, O. J.
Int. J. Chem. Kinet. 1997, 29, 627.
(8) Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Gill, P. M. W.;
Johnson, B. G.; Robb, M. A.; Cheeseman, J. R.; Keith, T.; Peterson, G. A.;
Montgomery, J. A.; Raghavachari, K.; Al-Laham, M. A.; Zakrzewski, V.
G.; Ortiz, J. V.; Foresman, J. B.; Cioslowski, J.; Stefanov, B. B.;
Using the above rate data for reactions 1a and 1b, it is
predicted that 90% of the total reaction proceeds through