Technical Comments

CH5+ Stability and Mass Spectrometry

Science  05 Nov 1999:
Vol. 286, Issue 5442, pp. 1051
DOI: 10.1126/science.286.5442.1051a

The report by E. T. White et al. (1) and the accompanying Perspective by D. Marx and M. Parrinello (2) both refer extensively to theoretical papers (3) which “indicate that the five protons in CH5 + are well bound to the central carbon … so that the C–H stretching potential is well defined but the angles between the five C–H bonds can change freely. The five equivalent protons are effectively swarming around the central carbon atom” (1, p. 136). The complexity of the system led Whiteet al. to present the infrared spectrum of CH5 + (near 2939 cm−1) without any specific assignments.

The theoretical models in these papers, however, appear to be in conflict with the Mass Spectral study of Heck, de Koning, and Nibbering (4), which showed that in the absence of intermolecular collisions, CH4D+ and CD4H+ are not rearranging but are stable. That these ions are relatively stable complexes of [CH3···HD]+ and [CD3···HD]+ was inferred from quenching reactions with NH3.

Assuming that the lifetime in the mass spectrometry experiments is of the order of milliseconds implies a significant activation energy, possibly exceeding 10 kcal/mol for the intramolecular scrambling process. Thus, CH5 + appears to be a particularly stable van der Waals complex of a methyl cation and a hydrogen molecule. It may well be described as a (3c-2e) system, but it is not scrambling. It may be a classic example of an ion-molecule interaction, but the existence of static isomeric CS structures (1, 2) with energies close to that of the (3c-2e) entity must be considered questionable.


CH5+ Stability and Mass Spectrometry

Response: We regret not to have sufficiently discussed in our report (1) earlier experimental results on the structure of CH5+, especially of the proton/deuteron transfer experiments (2, 3, 4), the cluster stability measurements by Hiraoka and his colleagues (5), and the predissociation spectroscopy of clusters by Boo and Lee (6). In particular, the paper by Heck, de Koning, and Nibbering (4) gives definitive experimental results indicating that the two (3c-2e) C–H bonds are isolated from the other three C–H bonds and do not scramble, a conclusion incompatible with the most recent ab initio theory (7, 8).

The high resolution of their mass spectrometer allowed Heck et al. (4) to extract CH4D+ and CD4H+that were uncontaminated with other ions with the same nominal mass. They then examined results of the proton (deuteron) transfer reactionEmbedded Image(1-1)With the use of a low-pressure CH5 + plasma source, they observed a 1:1 ratio of NH4 +/NH3D+ for both CH4D+ and CD4H+, and concluded that the two (3c-2e) C–H bonds were isolated from the rest and that all chemical reactions occurred via those two bonds. When the pressure of the CH5 + source was increased, the NH4 +/NH3D+ ratio approached to statistical values of 4/1 for CH4D+ and 1/4 for CD4H+.

The CH4D+ (CD4H+) was produced (4) from a 1:1 mixture of CH4and CD4 through the ion-neutral reactionsEmbedded Image Embedded ImageIn order to accept the conclusion made by Heck et al., it is necessary to assume that (i), CH4D+ (CD4H+) is produced by deuteron (proton) transfer and that the transferred deuteron (proton) ends up in one of the (3c-2e) bonds and not in the other three C–H bonds, without causing scrambling, and (ii) subsequent CH4D+ (CD4H+) collisions with CH4 or CD4 scramble isotopes. The collisions are weak, without exchange of proton and deuterons, because Heck et al. did not observe CH3D2 + or CH2D3 +. We find it difficult to accept these two results together in view of the fact that reaction (2) is exothermic by approximately 5 to 8 kcal/mol and that the collisions are “weak interactions between the protonated methane and methane” with energies of perhaps much less than 1 kcal/mol (4).

In view of our lack of knowledge of their experimental method and the limited space for text and references in our report (1), we chose not to go into the discussion given above. The same applies to discussions on works by Hiraoka and his colleagues (5), and Boo and Lee (6).

We stated that “our results are not inconsistent with the theoretical predictions of Schreiner et al. [7] and Müller et al. [8].” If the conclusion of Heck et al. is correct, the spectrum will be composed of five C–H stretch vibration bands with vibration-rotation interaction between them. Such a spectrum will also be complicated. The definitive spectroscopic conclusion will have to wait until we understand the reported CH5 +spectrum.


CH5+ Stability and Mass Spectrometry

Response by Marx and Parinello: Various evidence, particularly from theory, seems to suggest that protonated methane, CH5+, is a fluxional molecule [see (1)]. Contrary to this, mass spectroscopic reactivity experiments (2) on deuteronated methane CH4D+ and protonated per-deuteromethane CD4H+ are interpreted to show scrambling only if it is induced by intermolecular interactions, that is, each isotopomer “exhibits chemically distinguishable hydrogens” (2) without external perturbations. As recently stressed by Kramer, this implies that “in the absence of intermolecular collisions CH4D+ and CD4H+ are not rearranging but are stable” (3). As for some of the assumptions underlying the interpretation of these measurements, we draw attention to the objections put forward by Oka and White (4). In addition, the details of the reaction dynamics of both the formation and detection processes (such as preferred collision geometries and dipole locking) might yield crucial clues toward an understanding of the experimental data (2). Even weakly interacting ligands, such as several additional H2molecules leading to CH5+(H2)n complexes, can freeze the hydrogen scrambling motion in CH5+(5).

In an effort to shed light on this puzzle through theoretical means, we have undertaken further ab initio path integral simulations similar to the previous one performed for CH5 +(6). In particular, we have studied the isotopomer CH4D+, where the initial configuration was the optimized Cs ground-state structure of CH5 + (which may be pictured as consisting of a H2 moiety attached to a CH3 tripod by means of a three-center two-electron bond [see figure 1, left panel (1)], with one proton of the H2 moiety being substituted by a deuteron, that is, symbolically [CH3·HD]+. As for CH5 + (6),Cs -like configurations contribute most to the overall appearance of CH4D+ in its ground state. Most important, the protons, and in particular the deuteron, undergo scrambling (although the run was not long enough to lead to complete scrambling). These rearrangements lead to configurations—contributions to the nuclear density matrix in the sense of statistical sampling—where the deuteron is found in the CH2D tripod, thus leaving two protons in the three-center bond forming a H2 moiety, symbolically [CH2D·H2]+. The onset of scrambling was also observed in a much shorter run for a similarly prepared CD4H+ molecule. In conclusion, our calculations show hydrogen scrambling in CH5 +, CH4D+, and CD4H+. We cannot infer the time-scale of this process from our path integral calculations. We also point out other limitations of these calculations: the approximate nature of the density-functional-based potential energy surface; the neglect of quantum statistics (the Fermi-Dirac and Bose-Einstein exchange statistics for H and D); and the use of the Born-Oppenheimer approximation.

In a fluxional molecule, harmonic analysis is a severely limited tool. Nevertheless, it is instructive to compare—within the Born-Oppenheimer approximation—the harmonic zero-point vibrational energies (ZPE) of various isotopomers in the ground-state Cs structure (7). In the case of CH4D+, substitution of one proton by a deuteron in the CH3 tripod of CH5 + leads to isotopomers that have a lower ZPE in the range of ≈0.25 to 0.45 kcal/mol than those where substitution takes place in the three-center bonded H2 moiety. Exactly the reverse is true for isotope labelling in CD5 +. Here the two isotopomers that possess a mixed HD moiety are preferred by about 0.25 to 0.45 kcal/mol over the other two, where the single proton is located in the CH2H tripod.

These energy differences, although very small, are nonnegligible on the energy scale set by the Born-Oppenheimer energy barrier of only ≈0.6 kcal/mol (6)—the best literature value being ≈0.8 kcal/mol (8)—to the closest transition state ofC2v symmetry [see figure 1, right panel in (1)] leading to hydrogen scrambling. Anharmonicities, tunneling, thermal excitations, and rotational contributions will alter the reported ZPE differences. Nevertheless, it is to be expected that these differences will affect the probability distribution of the various isotopomers by biasing the ZPE-favored substitution sites.

In view of these arguments, it is likely that experiments performed systematically for a family of CH5 + isotopomers will lead to crucial novel insights into this fascinating molecular ion. As we concluded in our perspective (1), “CH5 + will certainly continue to challenge many groups in various fields of expertise for some time to come.”


Navigate This Article