3
688 J. Am. Chem. Soc., Vol. 118, No. 15, 1996
Samuni et al.
secondary ozonides calculated by ab initio methods13,17 were
found to agree satisfactorily with experimental data.
Table 1. Parameters of the Atom-Atom Pairwise Lennard-Jones
a
6
-12 Potentials Used in the Molecular Dynamics Simulation
atom pair
σ (Å)
ꢀ (cm-1)
In this paper we report that the primary and secondary
ozonides of ethylene can be formed at much lower temperatures
than previously observed, and that at a given temperature the
reaction’s progress depends on the matrix surrounding the
reacting pair. Thus, no reaction was found to take place in an
argon matrix at temperatures up to 35 K, while in a CO2 matrix
reaction products were identified at temperatures as low as 25
K. Molecular dynamics (MD) simulations of the trapping sites
in an argon matrix are reported, and it is shown that their
structure is unsuitable for ozonide formation. The apparently
more open structure of amorphous CO2 appears to allow the
reaction to proceed at temperatures as low as 25 K.
In all experiments both the primary and secondary ozonides
of ethylene were observed by their infrared spectrum, whenever
reaction took place. The assignment of the spectra was aided
by an ab initio calculation of the vibrational frequencies and
infrared intensities of these two compounds. The implication
of the simultaneous appearance of the ozonides (and the failure
to observe directly the Criegee intermediate) on the applicability
of the Criegee mechanism to this system is briefly discussed.
Ar-Ar
C-Ar
C-Ar
H-Ar
H-Ar
O-Ar
H-O
3.40
83
3.357
3.375
3.207
3.105
3.175
2.88
3.15
4.03
3.69
3.42
42 (for ethylene)
54 (for POZ)
18 (for ethylene)
22 (for POZ)
59
16
39
154
74
30
81
C-O
Xe-Xe
C-Xe
H-Xe
O-Xe
3.49
a
Sources: Reference 19. Hallam, H. E. Vibrational Spectroscopy
of Trapped Species; Wiley: London, 1973; p 45. The potential
parameters for unlike atoms, were derived from the literature atom-
19
atom parameters using the usual combination rules.
approach the surface under the action of the interaction potentials.
Periodic boundary conditions were imposed on the crystal surfaces
perpendicular to the growing surface, thus simulating an infinite crystal.
After the molecules were deposited, more argon atoms were deposited,
at least as many as used for the template. The total number of argon
atoms used in a typical simulation was 600. Many repetitions of the
simulations were performed, in which the initial locations of the
deposited atoms and molecules were varied at random, but always
chosen from a 300 K Boltzmann distribution. Different trapping sites
could be generated in different runs, and the most frequently encoun-
tered ones were considered the most stable ones, as discussed in more
detail elsewhere.18
Experimental and Computational Details
Ozone was prepared from a purified oxygen (98.5%, main impurity
N
2
) by electric discharge. Ethylene (99.5%) from Aldrich was used
as received; each gas was mixed separately with a carrier gas (argon
or CO ) to the desired concentration using standard manometric
2
techniques. An all-glass vacuum line was used for ozone, and a
stainless steel one for ethylene. The gases were codeposited (at a rate
of 2-10 mmol/h) on a KBr window attached to the cold tip of a closed
cycle helium cryostat (Air Products Model CS202). The infrared
spectrum of the resulting matrix was recorded by a Fourier transform
Ab initio calculations were carried out using the Gaussian 92 program
21
package at the HF or HF-MP2 levels. Several basis sets were used:
preliminary calculations were performed with a 6-31G, and more
22
sophisticated ones with a 6-311G* basis set. Complete structural
-
1
spectrometer (Nicolet Model 520, 0.5-cm resolution) at any desired
temperature higher than 12 K. The temperature was controlled by a
Lake Shore Cryogenics Temperature Controller (Model 330 Autotuning)
to within (0.5 K. Three characteristic temperatures were used in the
optimization was performed for each species, and vibrational frequen-
cies were calculated from the Hessian matrix using the harmonic
approximation at the optimized geometries. No imaginary frequencies
were found in the calculation, showing that the optimized structure
was indeed a minimum. The resulting frequencies were uniformly
experiments:
Tdep, the deposition temperature; Treac, the reaction
temperature; and Tspec, the temperature at which the spectra were
recorded. The reaction temperature was defined as the highest
temperature at which a reaction could be observed in a sample by
monitoring its infrared spectrum. The time intervals, treac, at which
the matrix was exposed to Treac varied between a few minutes and
23
scaled by 0.9427 for the smaller basis set and by 0.954 for the larger.
Results
Figure 1 shows a portion of spectra obtained in an argon
-
1
several hours. In a typical experiment, ozone in CO
2
(1:265) and
matrix between 920 and 980 cm . Ozone has no absorption
bands in this range, and ethylene was found to exhibit two
ethylene in CO (1:310) were codeposited at 15 K. The matrix was
2
warmed to 26.5 K in 128 min; during the heating period, spectra were
-1
absorption features at 946.9 and 959.7 cm , which are
recorded in the temperature intervals 16-17.4, 19-20.4, 22-23.4, and
apparently due to two different trapping sites of monomeric
2
5-26 K. The temperature was held at 26.5 K for 18 h; most of the
24
ethylene in an argon matrix. These two bands were observed
spectral changes occurred in the first 2 h, and after 10 h further changes
were minimal. The matrix was then warmed to 27 K in 10 min, and
held at that temperature for 12 h, without showing any further changes.
Warming gradually to higher temperatures at a very slow rate did not
for highly diluted ethylene in argon matrices (argon/ethylene
ratio ∼ 6000). A new band, which appeared only when both
ozone and ethylene were codeposited, is clearly seen at 951.1
-
1
cause further reaction until the CO
2
began to evaporate.
cm ; this band was not due to an oligomer of either ethylene
or ozone, as verified by depositing the pure compounds at
Molecular dynamics simulations were performed as previously
18
described, except that the molecules were treated as rigid bodies using
the RATTLE algorithm.1 Lennard-Jones pairwise atom-atom interac-
tion potentials were used, with the potential parameters listed in Table
(21) Frisch, M. J.; Trucks, G. W.; Gordon-Head, M.; Gill, P. M. W.;
9
Wong, M. W.; Foresman, J. B.; Johnson, B. G.; Schlegel, H. B.; Robb, M.
A.; Replogle, E. S.; Gomperts, R.; Andres, J. L.; Raghavachari, K.; Binkley,
J. S.; Gonzalez, C.; Martin, R. L.; Fox, D. J.; Defrees, D. J.; Baker, J.;
Stewart, J. J. P.; Pople, J. A. Gaussian 92, Revision F2; Gaussian, Inc.:
Pittsburg, PA, 1992.
1
. Full technical details of the simulations can be found in a separate
2
0
paper; briefly, the geometry of the molecules was taken from
microwave experimental data. The simulation was begun by construct-
ing a small crystal lattice, consisting of a few hundred atoms, that serves
as a template for further growth. Argon atoms were allowed to
approach the surface of the template, which is kept at a low constant
temperature. The molecules to be deposited were also allowed to
(
22) Hehre, W. J.; Ditchfield, R.; Pople, J. A. J. Chem. Phys. 1972, 56,
2
257.
(
23) Pople, J. A.; Scott, A. P.; Wong, M. W.; Radom, L. Isr. J. Chem.
1
993, 33, 345.
(24) Surprisingly little information is available on the spectra of ethylene
in an argon matrix. To our knowledge, the most dilute spectra (argon/
ethylene ratio ) 999) were reported by: Rytter, E.; Gruen, D. M.
Spectrochim. Acta 1979, 35A, 199. They found an incompletely resolved
(
(
(
17) McKee, M. L.; Rohlfig, C. L. J. Am. Chem. Soc. 1989, 111, 2497.
18) Fraenkel, R.; Haas, Y. Chem. Phys. 1994, 186, 185.
19) Allen, M. P.; Tildesley, D. J. Computer Simulations of Liquids;
-1
doublet at 948 and 959.5 cm . As Figure 1 shows, we find a well-resolved
-
1
Oxford University: New York, 1987.
20) Fraenkel, R.; Haas, Y. To be submitted for publication.
pair at 946.9 and 959.7 cm (these values vary somewhat with deposition
conditions).
(