Theoretical Analysis of the Methane Elimination from Oxonium Cations

DOI: 10.1021/jp952453+. Publication Date (Web): February 8, 1996. Copyright © 1996 American Chemical Society. Cite this:J. Phys. Chem. 100, 6, 2089-2...
0 downloads 0 Views 318KB Size
J. Phys. Chem. 1996, 100, 2089-2093

2089

Theoretical Analysis of the Methane Elimination from Oxonium Cations [R3O]+, R ) H, CH3 Minh Tho Nguyen† and Guy Bouchoux*,‡ Department of Chemistry, UniVersity of LeuVen, Celestijnenlaan 200F, 3001 LeuVen, Belgium, and De´ partement de Chimie, Laboratoire des Me´ canismes Re´ actionnels CNRS/URA 1307, Ecole Polytechnique, 91128 Palaiseau Cedex, France ReceiVed: August 21, 1995; In Final Form: NoVember 2, 1995X

Molecular orbital calculations up to the MP4/6-311++G(2df,2p)//MP2/6-31G(d,p) + ZPE level have been undertaken on the isomerization/dissociation processes of oxonium cations [CH3OH2]+, [(CH3)2OH]+, and [(CH3)3O]+. In addition to the classical oxonium cation structure, two ion-neutral complexes were identified: [H2CHO(R)‚‚‚CH3]+ and [H2CdO(R)‚‚‚CH4]+ (R ) H, CH3). The former is only a weakly stabilized structure while the latter (stabilized by 10-15 kJ/mol) is a reaction intermediate during the direct hydride ion abstraction [CH3]+ + CH3OR f CH4 + [CH2OR]+. The barrier height associated with the 1,2-elimination reaction [H2CHO(R)sCH3]+ f [H2CdO(R)]+ + CH4 is calculated to be 267 and 295 kJ/ mol for R ) H and CH3, respectively. This energy barrier is mainly due to the heterolytic CsO bond dissociation energy, and as a result, the corresponding transition structure looks like a [CH3]+‚‚‚CH3OR complex. Calculation has also been used to predict a heat of formation of the trimethyloxonium cation: ∆Hf°300K[O(CH3)3]+ ) 536 ( 15 kJ/mol.

Introduction A growing amount of accurate thermochemical data concerning gas phase ions has been collected in recent years.1 In spite of this rapid progress, the evaluation of activation parameters of elementary reactions involving these chemical species remains a great challenge for the spectroscopists. As far as mass spectrometry is concerned, reports on energy barriers of gas phase reactions of species generated by relevant techniques are rather scarce. Moreover, when available, the accuracy of reported results may be a matter of discussion (see for example refs 2 and 3). In this context it is of primordial interest to carry out a systematic comparison between the known experimental estimates and theoretical results obtained by, for example, highlevel quantum mechanical calculations, in order to establish the limit of confidence of both approaches in determining activation parameters. In the present contribution, we consider the energy barriers for methane elimination from oxonium cations (eqs 1 and 2), for which a theory-experiment comparison could be undertaken.

[(CH3)2OH]+ f CH4 + [CH2OH]+

(1)

[(CH3)3O]+ f CH4 + [CH2OCH3]+

(2)

Deuterium- and 13C-labeling experiments showed that both reactions 1 and 2 are 1,2-eliminations of CH4 involving one methyl group and one hydrogen of a vicinal methyl group; the proton bonded to oxygen in [(CH3)2OH]+ is not lost with the methane molecule. In addition, no hydrogen scrambling seems to occur prior to the decomposition and the carbon atoms play identical roles in each oxonium cation.4,5,8,9 Gas phase ionmolecule reactions of [CH3]+ with methanol and dimethyl ether have been recently carried out in an ion cyclotron resonance spectrometer.9 Formation of ions [CH2OH]+ or [CH2OCH3]+ †

University of Leuven. Ecole Polytechnique. X Abstract published in AdVance ACS Abstracts, January 1, 1996. ‡

0022-3654/96/20100-2089$12.00/0

was observed, suggesting an initial condensation process leading to the oxonium cations [(CH3)2OH]+ and [(CH3)3O]+, respectively. From labeling experiments, however, the following findings were reported: +

[ CH3] + CH3OH 13

+

[ CH3] + (CH3)2O 13

f 13CH4 + [CH2OH]+

43%

f CH4 + [13CH2OH]+

57%

f 13CH4 + [CH2OCH3]+

50%

f CH4 + [13CH2OCH3]+

50%

The reaction of [13CH3]+ with methanol leads to a 57/43 ratio of 13CH4/CH4 elimination whereas a statistical 1/1 ratio is expected from the oxonium cation. Similarly, the reaction of [13CH3]+ with dimethyl ether leads to a 50/50 ratio of 13CH4/ CH4 elimination whereas a statistical 1/2 ratio is expected. These results are readily explainable by the existence of two parallel reaction channels: the first one involves the expected oxonium cation while a second route is responsible for the exclusive 13CH elimination. For the latter process (which participates 4 to the composition to 14% and 25% for CH3OH and (CH3)2O, respectively), Audier et al.9 propose the intermediacy of an ionneutral complex, but it is clear that a direct hydride ion abstraction may also account for the observations. Recently, we have reported a theoretical estimate for the classical energy barrier of the CH4 loss from [(CH3)2OH]+ (reaction 1) using ab initio molecular orbital calculations.6 At the MP4/6-311++G(2df,2p)//MP2/6-31G(d,p) + ZPE level, the barrier height of reaction 1 has been calculated to be 267 ( 15 kJ/mol. In a subsequent spectrometric study, Holmes7 deduced a value of 243 ( 10 kJ/mol from the appearance energy measurement of [CH2OH]+ ions coming from metastable [(CH3)2OH]+ ions. The latter have been generated by dissociative ionization of 1-methoxy-2-hydroxypropene (eq 3). © 1996 American Chemical Society

2090 J. Phys. Chem., Vol. 100, No. 6, 1996

CH3-CH(OH)-CH2-OCH3 f [(CH3)2OH]+ + OCH3 (3) Although the experimental value is smaller than our theoretical estimate, the agreement between both theoretical and experimental values can be regarded as satisfactory as long as their respective lower and upper limits are overlapping each other. Concerning reaction 2, Wang et al.8 estimated an activation energy of 357 ( 19 kJ/mol on the basis of the observed [CH2OCH3]+ ion appearance curve from collision-induced dissociation of [(CH3)3O]+ under single-collision conditions in a flowing afterglow apparatus. From a crude evaluation based on the trend of [CH3]+ dissociation energy following methyl substitutions in the oxonium cations, we also suggested an upper limit of about 328 ( 15 kJ/mol for the energy barrier of reaction 2. This estimate is smaller than the experimental value of Wang et al.,8 but again, the theory-experiment comparison for both estimates is acceptable if we consider that upper and lower limits coincide. However, the theory-experiment difference for both oxonium ions is apparently not systematic. In this context, it is desirable to carry out ab initio calculations on the [(CH3)3O]+ system in order to obtain a more reliable theoretical estimate. It is the main purpose of the present work. As a necessary calibration, the simpler system [CH3OH2]+ has also been considered. On the other hand, we have searched for possible weakly-bonded complexes and their role in the methane elimination reaction. Method of Calculation All ab initio molecular orbital calculations were performed using a local version of the GAUSSIAN 90 set of programs.10 Stationary points have initially been located by geometry optimization using single-reference Hartree-Fock (HF) wave functions with the split valence 3-21G basis set and characterized by harmonic vibrational analysis at this level. Geometrical parameters of relevant structures have been reoptimized at the second-order Moller-Plesset perturbation theory (MP2) using a dp-polarized 6-31G(d,p) basis set. Improved relative energies have subsequently been obtained from single-point electronic energy calculations which utilize MP2-optimized geometries, larger 6-311++G(d,p) and 6-311++G(2df,2p) basis sets, and methods incorporating valence-electron correlation up to full fourth-order perturbation theory (MP4SDTQ). Calculated total and relative energies are summarized in Table 1. To facilitate comparison, we have used the same notation for each of the three oxonium systems considered, namely 1 for the oxonium cation, COM1 for the complex of [CH3]+, TS1 for the transition structure converting COM1 to 1, TS for the transition structure for H2 or CH4 elimination, and COM2 for the complex of CH4. Results and Discussion [CH3OH2]+ System. We first consider the methyloxonium cation, [CH3OH2]+. Since structures and dissociation pathways of this cation have been the subject of two previous studies,11,12 we only mention here the essentials of our results summarized in Figure 1. The two low-energy fragmentation modes of [CH3OH2]+ are losses of water and molecular hydrogen. It appears from Figure 1 that extension of the one-electron basis function, in particular the polarization space, induces rather small changes in relative energies. Comparison with available experimental values is satisfactory. Note that the proton affinity of methanol is calculated to be 763 and 752 kJ/mol at the MP4/6-311++G(d,p)//MP2/6-31G(d,p) and MP4/6-311++G(2df,2p)//MP2/6-

Nguyen and Bouchoux 31G(d,p) + ZPE levels, respectively. The corresponding experimental value is 760 kJ/mol at 300 K,14 i.e. 753 kJ/mol at 0 K, as deduced from thermal correction.12 Our results confirm two important points. The first is that the initial hydrogenbonded complex COM1 ([H2CH‚‚‚H2O]+) does exist as a local minimum but the energy barrier for converting this complex to the methyloxonium cation, which corresponds to a small reorientation of both fragments, is tiny (about 1 kJ/mol). The result suggests that such complexes are only transient species during ion-molecule reactions between [CH3]+ and H2O. The second observation is that both the fragments ([CH3]+ + H2O) and the transition structure for H2 loss, TS, seem to possess comparable potential energy. Within the expected uncertainty of ab initio methods, it is rather hard to calculate accurately the small energy difference between them. It may be recalled however that a value of 250 kJ/mol has been experimentally determined by Huntress and Bowers15 for the barrier [CH3OH2]+ f H2 + [CH2OH]+, in satisfactory agreement with our MP4/6-311++G(2df,2p)//MP2/6-31G(d,p) + ZPE and Uggerud’s12 estimate of 265 kJ/mol. In the same experimental study15 it has also been demonstrated that the loss of H2O needs more energy than the H2 elimination and occurs at its thermochemical threshold. [(CH3)2OH]+ System. The dissociation of the dimethyloxonium cation [(CH3)2OH]+ has been examined in detail in a recent paper.6 Figure 2 summarizes and completes the potential energy profile related to this ion. Additional results obtained in the present work must be mentioned: (i) A hydrogen-bonded complex between [CH3]+ and CH3OH, COM1, persists. The complexation energy amounts to 85 kJ/mol. Nevertheless, this complex disappears upon slight reorientation of both fragments to form the covalent C-O bond. This weak complex lies 12 kJ/mol (MP4/6-311++G(2df,2p)/ /MP2/6-31G(d,p)+ZPE) below the transition structure for CH4 loss, a value much smaller than that of about 60 kJ/mol in the methyloxonium ion (Figure 1). (ii) Dimethyloxonium cation is also the protonated form of dimethyl ether. The corresponding proton affinity is calculated to be 800 and 791 kJ/mol using MP4/6-311++G(d,p)//MP2/ 6-31G(d,p) and MP4/6-311++G(2df,2p)//MP2/6-31G(d,p) methods with zero point corrections, respectively. The latter value compares quite well with the experimental estimate for the 0 K proton affinity of dimethyl ether, PA((CH3)2O) ) 788 kJ/mol, which may be deduced from the experimental 300 K value of 793 kJ/mol14 and a correction of 5 kJ/mol for thermal energies.16 (iii) Starting from [CH3]+ and CH3OH, a hydride transfer from the methyl group of methanol to the methyl cation is also a possible reaction channel yielding methane. Due to the large difference in energy between ([CH3]+ + CH3OH) and ([CH2OH]+ + CH4), the hydride transfer occurs without an energy barrier, resulting in a lone complex (COM2, Figure 2). In fact, this complex has been identified several years ago by Swanton et al. during a theoretical investigation of the decomposition pathways of ethyloxonium.16 The hydride transfer corresponds thus to a situation having a single potential energy well. It is likely that this mode of CH4 formation, exclusively from the incoming [CH3]+ reactant, contributes to the observation of a nonstatistical distribution of the carbon atoms in the eliminated methane. [(CH3)3O]+ System. Selected optimized geometries of species relevant to the [(CH3)3O]+ system (1, COM1, TS1, TS, and COM2) are recorded in Figure 3. Parameters of the four fragments are omitted but are available upon request from the authors. A schematic potential energy profile is illustrated in Figure 4.

Methane Elimination from Oxonium Cations [R3O]+

J. Phys. Chem., Vol. 100, No. 6, 1996 2091

TABLE 1: Calculated Total (hartrees), Zero Point Vibrational (kJ/mol), and Relative Energies (kJ/mol) of the Oxonium Systems Considered

method experimentale MP2(F)/6-31G(d,p)b MP4/6-311++G(d,p) MP4/6-311++G(2df,2p) ZPEc

method

total energya of oxonium cation 1 [CH3OH2]+ 0 -115.697 36 -115.779 03 -115.840 29 total energya of oxonium cation 1 [(CH3)2OH]+

e

experimental MP2(F)/6-31G(d,p)b MP4/6-311++G(d,p) MP4/6-311++G(2df,2p) ZPEc

method experimentale MP2(F)/6-31G(d,p)b MP4/6-311++G(d,p) ZPEc

-154.884 62 -154.992 01 -155.076 40 total energya of oxonium cation 1 [(CH3)3O]+ -194.071 73 -194.204 39

1

COM1

TS1

TS

0 0 0 158

215 195 201 136

216 196 202 136

250 282 275 265 140

1

COM1

TS1

TS

0 0 0 0 232

266 247 254 212

267 248 255 212

243 280 264 267 214

1

COM1

TS1

TS

0 0 0 306

294 279 285

295 280 285

357 308 295 287

relative energyd [CH3]+ + H2O 283 296 264 268 129 relative energyd COM2 [CH3]+ + CH3OH 67 60 59 217

338 352 328 332 207

relative energyd COM2 [CH3]+ + (CH3)2O 51 43 288

386 382 364 280

[CH2OH]+ + H2 135 132 124 124 126

[CH2OH]+ + CH4 75 80 74 75 214

[CH2OCH3]+ + CH4 60 54 286

a Based on MP2(F)/6-31G(d,p) optimized geometries. Core orbitals are frozen in MP4SDTQ calculations. b Using full sets of molecular orbitals. Zero point vibrational energies from HF/3-21G calculations and scaled by 0.9. d All relative energies are corrected for ZPEs. e Experimental ∆H°, see text for references and discussion. c

Figure 1. Schematic potential energy profile for the two main reactions of [CH3OH2]+. The energy entries are MP4/6-311++G(d,p) (upper) and MP4/6-311++G(2df,2p) (lower) values, both based on MP2(F)/ 6-31G(d,p) geometries. In parentheses, are experimental values derived from available heats of formation at 300 K (ref 1 except ∆Hf°300([CH3OH2]+ ) 568 kJ/mol deduced from the PA(CH3OH) value of ref 14).

The trimethyloxonium ion 1 possesses a C3V symmetry in its most stable conformation and a strongly pyramidal oxygen atom. The O-C distance is consistently shortened upon methylation of the [H3O]+ ion, namely in going from 1.51 Å in [CH3OH2]+ to 1.49 Å in [(CH3)2OH]+ and to 1.48 Å in [(CH3)3O]+. The COC angle is also closed up: 116.2° in [(CH3)2OH]+ and 113.1° in [(CH3)3O]+. Calculations (at the MP4/6-311++G(d,p)//MP2/ 6-31G(d,p) + ZPE level) indicate that trimethyloxonium ion 1 is 364 and 54 kJ/mol below ([CH3]+ + O(CH3)2) and (CH4 + [CH2OCH3]+), respectively. The calculated difference in energy between the latter products is thus equal to 310 kJ/mol, a value

Figure 2. Schematic potential energy profile for the two main reactions of [(CH3)2OH]+. The energy entries are MP4/6-311++G(d,p) (upper) and MP4/6-311++G(2df,2p) (lower) values, both based on MP2(F)/ 6-31G(d,p) geometries. In parentheses are experimental values derived from available heats of formation at 300 K (ref 1 except ∆Hf°300([(CH3)2OH]+) ) 553 kJ/mol deduced from the PA((CH3)2O) value of ref 14).

corresponding with the experimental estimate of about 309 kJ/ mol (values at 0 K1,13). Together with available experimental heats of formation of the fragments,1,13 the above relative energies allow the heat of formation of trimethyl oxonium cation to be evaluated: ∆Hf°0K[O(CH3)3]+ ) 568 ( 2 kJ/mol and ∆Hf°300K[O(CH3)3]+ ) 537 ( 8 kJ/mol. Another way to estimate the heat of formation of the [O(CH3)3]+ cation is to consider isodesmic reactions 4 and 5:

2092 J. Phys. Chem., Vol. 100, No. 6, 1996

Nguyen and Bouchoux calculated at the same level, i.e. -76.564 006 hartrees, one obtains ∆E4 ) -7 kJ/mol and ∆E5 ) 2 kJ/mol, and consequently, if we consider that for isodesmic processes the difference in heat capacities between reactants and products is negligible, ∆Hf°300K[O(CH3)3]+ ) 535 ( 5 kJ/mol. Finally, we note that the two preceding estimates are in correct agreement with a 300 K value of 530 ( 8 kJ/mol, which may be calculated by a group equivalent analysis17 using the known heat of formation1,13 and proton affinity14 of H2O, CH3OH, and (CH3)2O. Consequently, if we put trust in MP4/6-311++G(d,p)//MP2/6-31G(d,p) calculations and assign an uncertainty of (15 kJ/mol, our suggested ∆Hf°300K[O(CH3)3]+ value is 536 ( 15 kJ/mol. Concerning the dissociation processes of oxonium cations 1, it is interesting to visualize the effect of methyl group substitution on the heterolytic bond energies associated with the C-O bond cleavage. The following dissociation energies (∆E at 0 K) have been obtained from MP4/6-311++G(d,p)//MP2/6-31G(d,p) + ZPE calculations:

Figure 3. Selected MP2(F)/6-31G(d,p) optimized geometries for structures relevant to the [(CH3)3O]+ system. Bond lengths are given in angstroms and bond angles in degrees.

Figure 4. Schematic potential energy profile for the two main reactions of [(CH3)3O]+. Relative energies are obtained at the MP4/6-311++G(d,p)//MP2(F)6-31G(d,p) level.

[(CH3)3O]+ + [H3O]+ f [CH3OH2]+ + [(CH3)2OH]+

∆E4 (4)

[(CH3)3O]+ + [CH3OH2]+ f 2[(CH3)2OH]+ ∆E5 (5) Using the MP4/6-311++G(d,p)//MP2/6-31G(d,p) total energies reported in Table 1 together with the total energy of [H3O]+

[CH3OH2]+ f [CH3]+ + H2O

∆E ) 264 kJ/mol

[(CH3)2OH]+ f [CH3]+ + CH3OH

∆E ) 328 kJ/mol

[(CH3)3O]+ f [CH3]+ + CH3OCH3

∆E ) 365 kJ/mol

Thus, an increase of 100 kJ/mol in dissociation energy follows two successive methyl substitutions. This remarkable effect is due to an increase of the number of (polarizable) methyl groups on the ionic reactant. In contrast, but for a similar reason, while the C-O bond cleavage becomes more difficult to achieve, the methane expulsion is getting less endothermic in going from 74 kJ/mol in reaction 1 to 54 kJ/mol in reaction 2. For these two reactions the stabilizing effect of the polarizable methyl group occurs on both the reactant and ionic product. A stronger ion-induced dipole interaction seems to be operative on the latter species, probably due to a more concentrated positive charge. Similar to the two previous cases, a hydrogen-bonded [H2CH‚‚‚O(CH3)2]+ complex (COM1) is found to lie 85 kJ/mol below the [CH3]+ + O(CH3)2 fragments, but it rearranges almost barrier-free, giving oxonium ion 1 (Figure 4). In many aspects the transition structure TS for 1,2-elimination of CH4 from 1 closely resembles that from [(CH3)2OH]+.6 As postulated by Wang et al.,8 TS can be regarded as a loose interacting complex between [CH3]+ and CH3OCH3. In structure TS (Figure 3), one carbon-oxygen bond is virtually broken (C3-O1 ) 2.71 Å) while the hydrogen atom H5 is not yet transferred to C3 (C3-H5 ) 2.40 Å) but prefers to remain attached to C2 (C2-H5 ) 1.10 Å). Again there appears a short distance between O1 and H4 (1.72 Å), indicating a certain hydrogen bond of the type H2C+H‚‚‚O(CH3)2. The latter is expected to bring about a stabilization partly compensating for the energy cost due to the C-O bond breaking. Note that unconstrained hydrogen bonding involving O1 and H4 interaction leads to complex COM1, for which a stabilization energy of 86 kJ/mol is calculated. In the case of transition structure TS a stabilization energy of ca. 69 kJ/mol with respect to ([CH3]+ + O(CH3)2) is calculated. In fact, TS is close in energy to the complex COM1, being only 15 kJ/mol above the latter. When TS is passed, the supersystem goes downhill, rearranging to the complex COM2 before dissociating to CH4 + [CH2OCH3]+. Our best estimate (MP4/6-311++G(d,p)//MP2/6-31G(d,p) + ZPE calculations) suggests an energy barrier of 295 ( 15 kJ/ mol for methane loss from trimethyloxonium cation 1, a value smaller than our earlier crude estimate of 328 ( 15 kJ/mol,

Methane Elimination from Oxonium Cations [R3O]+ which was expected to be only an upper limit.6 The present calculated value turns out to be 30 kJ/mol higher than the barrier height of 264 kJ/mol calculated for dimethyloxonium ion (reaction 1, Figure 2). This is consistent with the difference mentioned above in the heterolytic C-O bond dissociation energies in both cations which amount to 36 kJ/mol. These results lend support for the view that the barrier height for the elimination process from oxonium ions originates essentially from the heterolytic dissociation of the relevant bond, achieved by a specific hydrogen bond interaction which is almost equal to the fully developing hydrogen bond established in complex COM1. On the other hand, it is clearly shown that the experimental value of 357 ( 19 kJ/mol8 is definitely overestimated. Note also that the heterolytic bond dissociation energy of 1 determined by Wang et al.,8 i.e. 386 ( 19 kJ/mol, is also overestimated. The activation entropy accompanying the methane departure remains small but positive, ∆S°q ) 14 J/mol‚K including a negative entropic contribution of about 4 kJ/mol to the free energy of activation at 298 K. Again, a single hydride transfer from ether to methyl cation is found to be possible, yielding the complex COM2, which foreshadows the formation of CH4 (Figure 4). We believe that this transfer constitutes the second reaction route to form CH4 directly and exclusively from the reacting methyl cation in the bimolecular process [CH3]+ + (CH3)2O f CH4 + [CH2OCH3]+. Conclusion In summary, MO calculations reported in the present paper bring several pieces of information about the mechanism and the energetics of 1,2-elimination reactions of the type [H2CHsO(R)sCH3]+ f [H2CdO(R)]+ + CH4: (i) In addition to the classical oxonium cation structure 1, two ion-neutral complexes may be identified: [H2CHsO(R)‚‚‚CH3]+ (COM1) and [H2CdO(R)‚‚‚CH4]+ (COM2). (ii) COM1 is only a weakly stabilized structure; a slight reorientation of the two partners of the complex leads with an essentially negligible activation energy to structure 1. (iii) The responsibility of COM1 in the preferred loss of a methane molecule containing the initial methyl cation in the bimolecular process [CH3]+ + (CH3)2O f CH4 + [CH2OCH3]+ is only conceivable if cation 1 produced by the reaction COM1 f 1 dissociates before internal energy redistribution (nonergodic behavior). (iv) The direct hydride abstraction involving COM2 as reaction intermediate seems to provide a more classical interpretation of the nonstatistical distribution of labeled products. It is noteworthy that a larger participation of the hydride ion abstraction occurs for the system in which COM2 and 1 are closer in energy. This is in keeping with the intuitive idea that a competition between two association reactions is (partly) guided by the exothermicity of both processes. (v) The transition structure TS associated with the methane elimination from 1 looks like a [CH3]+‚‚‚CH3OH (or CH3OCH3)

J. Phys. Chem., Vol. 100, No. 6, 1996 2093 complex. Its binding energy is equal to 64 and 69 kJ/mol, respectively. In fact, the energy barrier for the methane elimination from oxonium cation 1 (calculated values: 267 and 295 kJ/mol, respectively) is due to the C-O bond dissociation energy; this energy demanding process is compensated by a favorable interaction between one hydrogen of the methyl cation and the oxygen atom. The binding energy of TS1 compares well with the stabilization energies of the ion-neutral complex COM1 (78 and 85 kJ/mol, respectively) in direct relation with the basicity of the neutral (PA ) 760 and 793 kJ/mol, respectively14). (vi) The heat of formation of trimethyloxonium cation [O(CH3)3]+ is estimated to be 567 kJ/mol at 0 K and 536 kJ/ mol at 300 K with a probable error of (15 kJ/mol. Acknowledgment. The authors thank the Belgian National Fund for Scientific Research (NFWO), Belgian Government (DWTC), KU Computing Centre, and CNRS for continuing support. References and Notes (1) (a) Lias, S. G.; Liebman, J. F.; Holmes, J. L.; Levin, R. D.; Mallard, W. G. Positive Ion Energetics Database, SRD Database 19A, 1993. (b) Bartmess, J. E. Negative Ion NIST Database, SRD Database 19B, 1994. (2) Nguyen, M. T.; Rademakers, J.; Martin, J. M. L. Chem. Phys. Lett. 1994, 221, 149. (3) Burgers, P. C.; Holmes, J. L.; Terlouw, J. K. J. Am. Chem. Soc. 1984, 106, 2762. (4) Farcasiu, D.; Pancirov, R. G. Int. J. Mass Spectrom. Ion Processes 1986, 74, 207. (5) Jarrold, M. F.; Illis, A. J.; Kirchner, N. J.; Bowers, M. T. Org. Mass Spectrom. 1983, 18, 388. (6) Nguyen, M. T.; Vanquickenborne, L. E.; Bouchoux, G. Int. J. Mass Spectrom. Ion Processes 1993, 124, R11. (7) (a) Holmes, J. L. 41st ASMS Conference on Mass Spectrometry, San-Francisco, 1993. (b) Sirois, M.; George, M.; Holmes, J. L. Org. Mass Spectrom. 1994, 29, 11. (8) Wang, D.; Squires, R. R.; Farcasiu, D. Int. J. Mass Spectrom. Ion Processes 1991, 107, R7. (9) (a) Audier, H. E.; McMahon, T. B. 41st ASMS Conference on Mass Spectrometry, San-Francisco, 1993. (b) Audier, H. E.; Koyanagi, G. K.; McMahon, T. B.; Tholman, D. To be published. (10) Frish, M. J.; Head-Gordon, M.; Trucks, G. W.; Foresman, J. B.; Schlegel, H. B.; Raghavachari, K.; Robb, M.; Binkley, J. S.; Gonzalez, C.; DeFrees, D. F.; Fox, D. J.; Whiteside, R. A.; Seeger, R.; Melius, C. F.; Baker, J.; Martin, R. L.; Kahn, L. R.; Stewart, J. J. P.; Topiol, S.; Pople, J. A. GAUSSIAN 90; Gaussian Inc.: Pittsburgh, PA, 1990. (11) Nobes, R. H.; Radom, L. Org. Mass Spectrom. 1982, 17, 340. (12) Uggerud, E. J. Am. Chem. Soc. 1994, 116, 6873. (13) Lias, S. G.; Bartmess, J. E.; Liebman, J. F.; Holmes, J. L.; Levin, R. D.; Mallard, W. G. J. Phys. Chem. Ref. Data 1988, 17, Suppl. 1. (14) Szulejko, J. E.; McMahon, T. B. J. Am. Chem. Soc. 1993, 115, 7839. (15) Huntress, W. T.; Bowers, M. T. Int. J. Mass Spectrom. Ion Phys. 1973, 12, 1. (16) Swanton, D. J.; Marsden, D. C. J.; Radom, L. Org. Mass Spectrom. 1991, 26, 227. (17) Szulejko, J. E.; Fisher, J. J.; McMahon, T. B.; Wronka, J. Int. J. Mass Spectrom. Ion Processes 1988, 183, 147.

JP952453+