Yanjun Haoab,
Yaoming Xieb and
Henry F. Schaefer III*b
aCollege of Physical Science and Technology, Sichuan University, Chengdu, Sichuan 610064, P. R. China
bCenter for Computational Quantum Chemistry, University of Georgia, Athens, GA 30602, USA. E-mail: qc@uga.edu
First published on 17th September 2014
The SiOOH potential energy surface has become central to the understanding of recent experiments (Science 2013, 342, 463) by Chakraborty associated with nebular meteorite formation. The entrance complex, transition states, and exit complex for the title reaction SiO + OH→ SiO2 + H have been studied using the CCSD(T) method with correlation consistent basis sets as large as cc-pV(Q+d)Z. Reported here are characteristics of the reactants, products, six transition states, and four intermediate complexes for this reaction. These show four previously undiscovered stationary point geometries. The entrance complex OH⋯OSi is predicted to lie 28.6 kJ mol−1 below the separated reactants. The classical barriers cis-TS1 and trans-TS1 are predicted to lie 21.8 kJ mol−1 and 6.8 kJ mol−1, respectively, below the reactants. The exit complex HSiO2 is bound by 115.3 kJ mol−1 relative to the separated products. After zero-point vibrational energy corrections, the reaction energy is predicted to be −1.4 kJ mol−1. Vibrational frequencies of the stationary points are reported and compared with the limited available experimental results. The SiOOH potential surface is found to be very different from that for COOH, contrary to the analogy drawn by Chakraborty. Notwithstanding, the assumption of Chakraborty appears justified, because all the stationary points for the SiO + OH reaction have lower relative energies than known for the analogous carbon system.
In 1996, Chagger, et al.9 studied silica formation from hexamethyldisiloxane oxidation by means of a CH4–N2/air opposed diffusion flame technique, and the reaction mechanism involved the reaction noted above, SiO + OH → SiO2 + H. These authors also listed the available reaction enthalpy and rate constant. In 2009, the kinetics of the reaction SiO + OH → SiO2 + H were investigated by Martin, Blitz and Plane10 using the pulsed laser photolysis of a Si atom precursor in the presence of O2, followed by time-resolved non-resonant laser-induced fluorescence. They measured the rate constant to be (5.7 ± 2.0) × 10−12 cm3 per molecule per second at 293 K and 4–20 Torr.
Theoretically, in 1995, Zachariah and Tsang11 predicted features of the potential energy surface for the reaction SiO + OH → SiO2 + H with the MP4/6-311(2df,p)//MP2/6-31(d) method. Recently, Martin, Blitz and Plane10 obtained a potential energy profile of this reaction using CBS-Q single point energies at geometries optimized with the B3LYP/6-311+G(2d,p) method. As Martin et al. discussed, the low pressure limiting rate coefficient at room temperature (SiO + OH → HOSiO → SiO2 + H) obtained by Zachariah and Tsang11 is 15 times smaller than their experimental determination and Rice–Ramsperger–Kassel–Markus (RRKM) fit; while it is 10 times larger than their estimate for the cis-HOSiO → SiO2 + H pathway. It is well known that deduced rate constants depend strongly on the quality of the potential energy surface adopted.12 Although the research of ref. 10 and 11 is very valuable, higher level theoretical studies are now possible and appropriate. In the present study, we will adopt very high level theoretical methods up to CCSD(T)/cc-pV(Q+d)Z to predict a more complete potential energy surface, which will give new possible pathways and more reliable activation energies for the SiO + OH → SiO2 + H reaction. We hope that the newer results will provide a solid base for the future kinetic studies.
In our research, coupled cluster computations were carried out to explore the stationary points on the SiO + OH → SiO2 + H potential energy surface, including the reactants, products, transition states, and intermediate complexes. All geometries are fully optimized using single and double excitation coupled-cluster theory with perturbative triple excitations CCSD(T).16,17 The Dunning correlation-consistent cc-pVnZ (n = D, T, Q) basis sets were adopted for the hydrogen and oxygen atoms,18 while the cc-pV(n+d)Z (n = D, T, Q) basis sets19 were utilized for the silicon atom. With the coupled-cluster methods, the 1s-like molecular orbital for oxygen is frozen, while the 1s2s2p-like molecular orbitals for silicon are frozen. Harmonic vibrational frequencies predicted at the same level were used to evaluate the zero point vibrational energies (ZPVE) and for the characterization of the stationary points. All open-shell Hartree–Fock computations were unrestricted (UHF), and the coupled cluster methods adopted were also unrestricted, UCCSD(T). All coupled cluster computations were carried out using the CFOUR program.20
Fig. 1 shows that there are two possible pathways for the SiO + OH → SiO2 + H reaction, and for each pathway there is more than one transition state. However, the two pathways share same entrance complex and exit complex. The entrance complex OH⋯OSi is of linear geometry with C∞v symmetry, which was not reported in previous theoretical studies.10,11 This complex is predicted to lie 28.6 kJ mol−1 below the separated reactants (SiO + OH) at the CCSD(T)/cc-pV(Q+d)Z level of theory.
For the valence isoelectronic HO + CO → H + CO2 reaction, Ma, Li, and Guo in 2012 reported two entrance complexes.24 Beside the OHCO complex on the cis pathway, there is another linear complex OH⋯CO on the trans pathway. It may be noted that the electronegativity of silicon is less than that of carbon, and a generic Si⋯H interaction may be weaker than the comparable C–H interaction. The greater importance of polarization functions (d orbitals) for silicon may also be pertinent here. In fact, we have attempted to optimize the OHSiO structure analogous to OHCO, but it has at least two imaginary vibrational frequencies, and it eventually collapses to the OHOSi complex seen on the left of Fig. 1. The predicted OH⋯OSi distance for our single entrance complex is 1.90 Å, which lies in the typical range of hydrogen bonding distances. It is significantly shorter than the analogous OH⋯OC distance (2.30 Å) predicted via the CCSD(T)-F12b/aug-cc-pVTZ method by Xie, Li, Xie, and Guo.25 The exit complex HSiO2 (Fig. 1) for the SiO + OH→ SiO2 + H reaction is predicted to be a planar structure with C2v symmetry, lying 10.5 − (−104.8) = 115.3 kJ mol−1 below the products SiO2 + H. This exit complex was reported by Martin, Blitz, and Plane to lie below separated SiO2 + H by 96.1 kJ mol−1,10 but no such a complex was reported in Zachariah and Tsang's paper.11 The Si–O and H–Si bond lengths for C2v symmetrical HSiO2 are 1.565 and 1.467 Å, respectively, which are not inconsistent with the B3LYP/6-311+G(2d,p) optimized results (1.523 Å for Si–O and 1.481 Å for H–Si) of Plane and coworkers.10
From the entrance complex to the exit complex, there are two pathways. The first pathway (red in Fig. 1) proceeds via a cis-HOSiO intermediate, and the other pathway (blue in Fig. 1) involves a trans-HOSiO intermediate. The cis and trans isomers have similar bond distances, bond angles, and relative energies. The cis isomer lies below the reactants by 265.0 kJ mol−1, with the trans isomer is lower than reactants by 262.4 kJ mol−1.
Between the entrance complex and the cis or trans intermediates, there are transition states, i.e., cis-TS1 and trans-TS1, which were not reported before.10,11 Each of these has one imaginary vibrational frequency, namely, 81i cm−1 ((Q+d)Z) for cis-TS1 and 123i cm−1 for trans-TS1. Structures cis-TS1 and trans-TS1 are predicted to lie 21.8 kJ mol−1 and 6.8 kJ mol−1, respectively, below the reactants at the CCSD(T)/cc-pV(Q+d)Z level of theory. These two transition states are very different in geometry, although the nonbonded OSi⋯O distances in the two structures are comparable (∼3.0 Å). Specifically, trans-TS1 has a chain structure, while the cis transition state has a cyclic structure, with an O⋯H hydrogen bond facilitated by the cis geometry. Accordingly, the cis isomer has a lower energy than the trans-TS1, by 15.0 kJ mol−1. Compared with the unique entrance complex OHOSi, the barriers of these two transition states are [(−6.8) − (−28.6)] = 21.8 (trans) and [(−21.8) − (−28.6)] = 6.8 (cis) kJ mol−1, respectively. After zero-point vibrational energy corrections, the barriers (with respect to the entrance complex) for trans-TS1 and cis-TS1 decrease to [(−5.0) − (−22.3)] = 17.3 and [(−16.1) − (−22.3)] = 6.2 kJ mol−1 (Fig. 2).
![]() | ||
Fig. 2 Stationary points on the SiO + OH potential energy surface after zero-point vibrational energy (ZPVE) corrections. |
From the cis and trans intermediates to the unique exit complex (HSiO2), there are corresponding cis and trans transition states, i.e., cis-TS2 and trans-TS2. Each of these is predicted to display one imaginary vibrational frequency, i.e., 1271i cm−1 for cis-TS2 and 1179i cm−1 for trans-TS2 with the cc-pV(Q+d)Z basis sets. The corresponding normal modes reveal Si–O bond formation and O–H bond breaking as the reaction proceeds toward products SiO2 + H. Since the cis-TS2 and trans-TS2 transition states have rather different geometries, they are predicted to have divergent relative energies, 33 kJ mol−1 for cis-TS2 and −60 kJ mol−1 for trans-TS2. The latter prediction is consistent with the previous theoretical value of 63 kJ mol−1 by Martin, Blitz and Plane.10 The former prediction is comparable to the theoretical value of 22 kJ mol−1,10 (note that 12.13 in Fig. 10 of ref. 10 may be a typographical error for 22.13 if their ESI data are correct).
The isomerization between the two intermediates (trans-HOSiO to cis-HOSiO) involves two transition states (trans-TS3 and gauche-TS3, both green in Fig. 1), which are neither on cis nor on trans reaction pathway. The barrier of trans-TS3 is much higher than that of gauche-TS3. This is because the trans-TS3 (Cs symmetry) structure is confined to the plane, and the isomerization is by way of the Si–O–H angle change in the plane. In contrast, gauche-TS3 (C1 symmetry) is through the lower energy internal-rotation around the Si–O single bond. The barrier for gauche-TS3 (13.5 kJ mol−1 from trans-HOSiO side) is in good agreement with the barrier of 11.8 kJ mol−1 from Martin, Blitz and Plane,10 and with the 10 kJ mol−1 reported by Zachariah and Tsang.11 The trans-TS3 transition state was not mentioned in previous theoretical studies.10,11
The theoretical relative energies (ΔE) with respect to the reactants (OH + SiO), harmonic vibrational frequencies (cm−1), and the zero-point vibrational energies (kJ mol−1) for the stationary points on the SiO + OH→ SiO2 + H potential energy surface are reported in Table 1. It is noted that the predicted results for SiO2, SiO, and OH are in good agreement with the experimental vibrational frequencies.21,26
ΔE (kJ mol−1) | ZPVE (kJ mol−1) | Vibrational frequencies (cm−1) | |||||||
---|---|---|---|---|---|---|---|---|---|
ω1 | ω2 | ω3 | ω4 | ω5 | ω6 | ω7 | |||
a Fundamental frequencies.26b Harmonic vibrational frequencies.21 | |||||||||
OHOSi | −28.6 | 36.12 | 3641 | 1244 | 459 | 459 | 163 | 37 | 37 |
cis-TS1 | −21.8 | 35.51 | 3706 | 1226 | 477 | 361 | 165 | 81i | |
cis-HOSiO | −265.0 | 44.15 | 3831 | 1201 | 893 | 772 | 389 | 294 | |
cis-TS2 | 32.8 | 17.91 | 1395 | 957 | 285 | 282 | 75 | 1271i | |
HSiO2 | −104.8 | 32.78 | 2323 | 1028 | 895 | 567 | 364 | 304 | |
trans-HOSiO | −262.4 | 44.28 | 3875 | 1218 | 860 | 789 | 353 | 307 | |
trans-TS1 | −6.8 | 31.68 | 3742 | 1235 | 171 | 101 | 47 | 123i | |
trans-TS2 | −59.5 | 28.28 | 1883 | 1302 | 893 | 350 | 299 | 1179i | |
gauche-TS3 | −248.9 | 41.31 | 3897 | 1207 | 868 | 621 | 314 | 384i | |
trans-TS3 | −6.2 | 31.68 | 3748 | 1235 | 185 | 87 | 42 | 96i | |
SiO | 7.40 | 1237 | |||||||
OH | 22.43 | 3750 | |||||||
SiO + OH | 0.0 | 29.83 | |||||||
SiO2 | 17.95 | 1436 | 989 | 288 | 288 | ||||
SiO2 + H | 10.5 | 17.95 | |||||||
![]() |
|||||||||
Experiment | |||||||||
SiO2a | 17.78 | 1416 | 1010 | 273 | 273 | ||||
SiOb | 7.43 | 1242 | |||||||
OHb | 22.36 | 3738 |
A sketch of the CCSD(T) potential energy surface after zero-point vibrational energy (ZPVE) corrections is presented in Fig. 2. Relative to the reactants (SiO + OH), the energies for the intermediate complexes and transition states, except for cis-TS2 and trans-TS2, slightly increase, compared with those values without ZPVE corrections.
The SiO + OH → SiO2 + H reaction is predicted to be endothermic by 10.5 kJ mol−1 at the CCSD(T)/cc-pV(Q+d)Z level of theory without ZVPE corrections. After the ZPVE corrections, the SiO + OH endothermicity with the cc-pV(D+d)Z, cc-pV(T+d)Z, and cc-pV(Q+d)Z basis sets is predicted to be 39.4, 4.8, and −1.4 kJ mol−1, respectively. The best predicted reaction energy with ZPVE correction (−1.4 kJ mol−1 with the cc-pV(Q+d)Z basis set) is smaller in magnitude than the previous (exothermic) theoretical values of −7 kJ mol−1 derived from their theoretical heats of formation at 298.15 K by Zachariah and Tsang11 and −6.6 kJ mol−1 (CBS-Q) by Martin, Blitz and Plane.10 Our theoretical reaction energy (−1.4 kJ mol−1) is somewhat different from the experimental value (−24.9 kJ mol−1 at 0 K, and −26.0 kJ mol−1 at 298.15 K) derived from the NIST heats of formation.27 However, there is another experimental reaction enthalpy −0.4 kJ mol−1 at 298.15 K reported by Chagger et al.,9 who adopted different heat of formation of −279.9 ± 17 kJ mol−1 for SiO2 reported by that of Hildenbrand et al. in 1994.28 It is notable that there is a significant uncertainty in the heat of formation of SiO2, so this comparison with experimental data may not be entirely helpful. Actually Hildenbrand and Lau also report a heat of formation for a HOSiO species,28 which is far from all theoretical values.
For this reaction, there are two possible pathways. One pathway may be called the cis pathway (red in Fig. 1), i.e., SiO + OH (reactants) → entrance complex → cis-TS1 → cis-HOSiO (intermediate) → cis-TS2 → exit complex → SiO2 + H (products). In this pathway the highest barrier is 32.8 kJ mol−1 (cis-TS2). Another pathway may be called trans pathway (blue in Fig. 1), i.e., SiO + OH (reactants) → entrance complex → trans-TS1 → trans-HOSiO (intermediate) → trans-TS2 → exit complex → SiO2 + H (products). In this pathway the highest barrier is −6.8 kcal mol−1 (the trans-TS1 structure). Clearly, the trans pathway is more favorable than the cis pathway. However, with the aid of the transition states between cis and trans intermediates (green in Fig. 1), we can have a slightly more complicated pathway with even lower barrier. That is, SiO + OH (reactants) → entrance complex → cis-TS1 → cis-HOSiO (intermediate) → gauche-TS3 → trans-HOSiO (intermediate) → trans-TS2 → exit complex → SiO2 + H (products). In this cis–trans combination pathway, the highest lying transition state lies −21.8 kJ mol−1 (cis-TS1) without ZPVE or −16.1 kJ mol−1 with ZPVE.
Compared with the previous theoretical PES of ref. 10 and 11, we have found some new stationary points, i.e., the entrance complex OHOSi, trans-TS1, cis-TS1, and trans-TS3. We suggest that both cis and trans pathway proceed via the exit complex HSiO2. The transition state cis-TS2 does not directly connect to the products because there is no barrier (lower than 115.3 kJ mol−1) from the products to the exit complex. Based on our more complete and reliable PES, we suggest a possible pathway through cis-TS1 (−21.8 kJ mol−1) and trans-TS2 (−59.5 kJ mol−1).
In Fig. 3, we have redrawn (to match our Fig. 2 as closely as possible) the reliable energetics of the OH + CO system predicted by Li, Wang, Jiang, Ma, Dawes, Xie, Bowman, and Guo.29 One sees very quickly that SiOOH and COOH potential surfaces differ in many respects. First, the endothermic SiO + OH → SiO2 + H reaction is barrierless. Other than the products and reactants, the highest lying stationary point, cis-TS1, lies 16.1 kJ mol−1 below the reactants. For the exothermic CO + OH → CO2 + H, there is a barrier (cis-TS2) lying 10.1 kJ mol−1 above the reactants. Secondly, the trans and cis HOSiO minima lie much lower (−248.0 and −250.7 kJ mol−1) than the analogous equilibrium structures of trans and cis HOCO (−123.8 and −116.4 kJ mol−1). One should also note that the HCO2 minimum is separated by a barrier from products CO2 + H, while the HSiO2 equilibrium goes directly uphill to products SiO2 + H. Another way of saying this is that SiO2 + H reacts without a barrier, while the CO2 + H reaction faces a significant (58.2 kJ mol−1) barrier, consistent with the remarkable stability of CO2.
![]() | ||
Fig. 3 The CO + OH reaction pathways analogous to our SiO + OH system. Redrawn from Fig. 1 of Bowman, Guo, and coworkers.29 |
Despite the major differences between the COOH and SiOOH energy landscapes, the assumption of Chakraborty, Yanchulova, and Thiemens6 seems to be justified. Indeed the transition states and intermediates of the SiO + OH system are all more accessible energetically than those for the CO + OH system.
This journal is © The Royal Society of Chemistry 2014 |