Collision-induced Energy Transfer and Bond Dissociation in Toluene by H<sub>2</sub>/D<sub>2</sub>
Collision-induced Energy Transfer and Bond Dissociation in Toluene by H2/D2
Bulletin of the Korean Chemical Society. 2013. Dec, 34(12): 3641-3648
Copyright © 2013, Korea Chemical Society
  • Received : August 15, 2013
  • Accepted : September 12, 2013
  • Published : December 20, 2013
Export by style
Cited by
About the Authors
Jongbaik Ree
Yoo Hang Kim
Department of Chemistry, Inha University, Incheon 402-751, Korea
Hyung Kyu Shin
Department of Chemistry, University of Nevada, Reno, Nevada 89557, USA

Energy transfer and bond dissociation of C-H methyl and C-H ring in excited toluene in the collision with H 2 and D 2 have been studied by use of classical trajectory procedures at 300 K. Energy lost by the vibrationally excited toluene to the ground-state H 2 /D 2 is not large, but the amount increases with increasing vibrational excitation from 5000 and 40,000 cm −1 . The principal energy transfer pathway is vibration to translation ( V-T ) in both systems. The vibration to vibration ( V-V ) step is important in toluene + D 2 , but plays a minor role in toluene + H 2 . When the incident molecule is also vibrationally excited, toluene loses energy to D 2 , whereas it gains energy from H 2 instead. The overall extent of energy loss is greater in toluene + D 2 than that in toluene + H 2 . The different efficiency of the energy transfer pathways in two collisions is mainly due to the near-resonant condition between D 2 and C-H vibrations. Collision-induced dissociation of C-H methyl and C-H ring bonds occurs when highly excited toluene (55,000-70,400 cm −1 ) interacts with the ground-state H 2 /D 2 . Dissociation probabilities are low (10 −5 ~10 −2 ) but increase exponentially with rising vibrational excitation. Intramolecular energy flow between the excited C-H bonds occurring on a subpicosecond timescale is responsible for the bond dissociation.
Studies on the rates and mechanisms of energy transfer in molecular collisions are critical to understanding elementary processes which play a vital role in chemical reactions. 1 - 13 The collision-induced relaxation of vibrationally excited polyatomic molecules is of particular interest as it directly affects the fate of bond dissociation and subsequent reaction dynamics. Energy transfer involving polyatomics has been studied experimentally by a variety of techniques. Experimental studies have been commonly carried out using timeresolved infrared flourescence (TR-IRF) 4 - 6 or ultraviolet absorption (UVA) techniques. 7 - 9 Several studies have used techniques based on photothermal processes. 10 - 12 Through these studies the researchers have provided insight into relative collision efficiencies and their dependence on the energy of excited molecules. Recently, Hsu et al. have studied the energy transfer of highly vibrationally excited molecules using crossed-beam techniques. 13 However, theoretical studies of energy transfer processes in polyatomic molecules have lagged behind the experimental studies as models and collision dynamics of many-body interactions are difficult to develop. 14 - 22 When such studies are developed using physically realistic interaction models, the results can complement experimental studies and guide future studies.
In collisions involving a large organic molecule, the average amount of energy transfer per collision is not very large. 23 - 31 One such molecule is toluene, which contains two distinctly different C-H bonds, namely the methyl C-H and ring C-H bonds, along with many C-C bonds. When toluene undergoes a collision with other molecules, both intermolecular energy transfer and intramolecular energy flow from one CH bond to others via C-C bonds can occur. 22 , 31 Hippler and co-workers studied the deactivation of the excited toluene by about 60 different collider gases using the UVA techniques to follow the collisional relaxation of cycloheptatriene. 32 They prepared highly vibrationally excited toluene (52,000 cm −1 = 150 kcal/mol) after irradiating cycloheptatriene and showed the excited molecule dissociates into benzyl radical and H atom, which indicate that the extent of energy transfer can be significant. Toselli and co-workers investigated the collisional loss of vibrational energy from gas-phase toluene by 20 collider gases by monitoring the TR-IRF from the C-H stretch modes near 3.3 mm. 23 Computational studies of collision- induced energy transfer between aromatic polyatomic molecules including toluene have been reported by Bernshtein and Oref. 33
The purpose of this paper is to study energy transfer and bond dissociation in the collision of highly vibrationally excited toluene with H 2 and D 2 using quasiclassical trajectory procedures. We use the model presented in Ref. 22, where the toluene + N 2 /O 2 collision was considered. In the present systems, toluene is in interaction with the light molecules, which create a collision environment greatly different from that of N 2 or O 2 . The vibrational frequencies are now more than twice those of the latter molecules, thus making the contribution of intermolecular vibrational energy transfer more significant. In particular, the vibrational frequency of D 2 is very close to C-H frequency in toluene, creating the near-resonant condition which can lead to increased vibrational energy transfer. The low moments of inertia of H 2 /D 2 lead larger rotational spacings, which in turn lead to a significant contribution of rotational energy transfer. The masses of these molecules are low, so the collision dynamics can be significantly different from that involving heavier masses such as N 2 /O 2 . Further, a comparative study of toluene + H 2 and toluene + D 2 is important as it can reveal the effects of isotopic substitutions on the transfer of energy and bond dissociation in large molecules. To study these aspects, we start with a vibrationally excited toluene (5000-40,000 cm −1 ) and determine energy transfer between toluene and H 2 /D 2 as a function of the toluene vibrational energy content and compare the results with experimental data. Energy transfer processes occurring in the collision are vibration-to-translation ( V-T ), vibration-to-rotation ( V-R ) and vibration-tovibration ( V-V ) pathways. We then take highly excited toluene (55,000-70,400 cm −1 ) to study its collision dynamics with H 2 /D 2 leading to bond dissociation with particular emphases on the importance of intramolecular energy flow between the excited methyl and ring C-H bonds. Dissociation of ring C-H or methyl C-H can occur when the initial energy stored in the molecule redistribute intramolecularly on collision, the process which will be discussed in detail.
Interaction Model and Energies
The model and internal coordinates of toluene have been reported in Ref. 22 for the collision of toluene with N 2 /O 2 and the same interaction model will be adopted here whenever applicable. We briefly recapitulate the essential aspects of the model for clarity using the same notations and conventions used in Ref. 22. We describe the collision of toluene with hydrogen molecule in terms of 39 intramolecular coordinates of the planar ring and the rotational and vibrational modes of the incoming molecule. Toluene is considered to be non-rotating during the collision.
The C-H ring bond is in the clockwise side from C-H methyl. The internal modes of toluene are 6 stretches (χ 1 , χ 2 , χ 3 , χ 4 , χ 14 , χ 15 ) and 12 bends ( φ 1 , φ 2 , φ 3 , φ 4 , φ 5 , φ 6 , φ 7 , φ 1 ', φ 1 ", φ b , φ b ', φ b ") between C-H ring and C-Hm ethyl , the region we consider to be the interaction (or primary) zone, where the C-H methyl and C-H ring bonds are in the direct interaction with H 2 ; see Figure 1 . Other coordinates included are the χ 5 9 (C-C) ring stretches, χ 10 13 (C-H) ring stretches and φ 8 - φ 19 bends around the carbon atoms C 3 , C 4 , C 5 and C 6 , a total of 9 stretches and 12 bends, in the region remote from the primary zone. We refer this region as the secondary zone.
PPT Slide
Lager Image
Collision model. The stretching and bending coordinates of vibrations are defined. All carbon atoms and ring H atoms are coplanar and modes are numbered clockwise. The star denotes the center-of-mass (c.m.) of toluene. The displacement of each vibrational coordinate is identified by χe1 + χ1, etc. (see the text).
The interaction energies needed to describe the collision of H 2 with toluene must contain terms responsible for the coupling of the relative motion with the ring C-H stretch and methyl group C-H stretch as well as the coupling between the stretches and bends. We first identify the two atoms of H 2 as H (1) and H (2) and denote the interatomic distances r 1 - r 4 as shown in Figure 1 , where r represents the distance between the center-of-mass (c.m.) of H 2 and the c.m. of toluene describing the relative motion of the collision system. The interatomic distances can be expressed in terms of the instantaneous coordinates of the C-H methyl bond, C-CH methyl bond, (C-C) ring bond, C-H ring bond, C-C-H methyl bend, C-CC methyl bend, C-C-Hring bend and H-H bond. That is, r i = r i ( r , χ 1 , χ 2 , χ 3 , χ 4 , φ 1 , φ 2 , φ 5 , χ, θ, η, η') for i = 1-4, where χ is the displacement of the H 2 bond from its equilibrium distance d , and θ is the angle of incidence defined in Figure 1 , η and η' are the rotational angles of H 2 . These atom-atom distances are given in Ref. 22. The coupling of these modes with others including those of the inner zone will have to be considered in formulating the overall interaction energy.
The overall interaction energy is the sum of the Morsetype intermolecular terms, Morse-type stretching terms and the harmonic bending terms of toluene, intramolecular coupling terms, the incident molecule’s term and London interaction,
PPT Slide
Lager Image
U ( r β ) = D [ e –(rrβ ) /a –2 e (r– r β )/2a ] , β = 1-4
Us ( xi ) = Σ Di [1– e xi/bi ] 2 , i = 1-1 5
PPT Slide
Lager Image
U int ( YiYj ) = Σ Kij ( Yi –Yei ) ( Yj – Yej ) , i ≠ j
PPT Slide
Lager Image
The values of potential and spectroscopic constants for toluene are listed in Table 1 . 34 - 41 Each bond length will be denoted by (χ ei + χ i ), where χ i is the displacement of the bond length from its equilibrium value χ ei . Similarly, we express each bending coordinate as ( φ ej + φ j ), where φ j is the displacement of the jth bending vibration from the equilibrium angle φ ej . In the London term, Ii is the ionization potential and α i is the polarizability of toluene and H 2 /D 2 . The “inc” in Eq. (1) means the incident molecule, H 2 or D 2 . For the intermolecular interaction, we take D = 116.85 k B and 126.94 k B for toluene + H 2 and toluene + D 2 , respectively, the Lennard- Jones (LJ) parameter for the well depth for the collision pairs calculated by the usual combining rule, 42 where k B is the Boltzmann constant, and a = 0.238 and 0.235 Å for toluene + H 2 and toluene + D 2 , respectively. 41 , 42 For the i th stretch, we use b i = (2 D i i ω ; i 2 ) 1/2 , where D i = D o,i + ½ħ ω ei , to determine the exponential range parameters listed in Table 1 . In the coupling terms, Y i = x i or φ i and the coupling constant K ij are taken from Xie and Boggs' ab initio calculations. 43 Table 2 lists the potential and spectroscopic constants for H 2 and D 2 . 44 - 46
Potential and spectroscopic constants for toluene
PPT Slide
Lager Image
#Number i represents the subscript of the vibrational cooridinate x shown in Figure 1. aReference 34-36. breference 37. creference 38; the dissociation energy Di used throughout the paper including Eq. (1) is Di = Doi + ½ħ ωei. dreference 39, 40. ereference 41. freference 40
The equations of motion which determine the time evolution of the relative motion, 15 stretches, 24 bends and three motions of the incident molecule for the given incident angle can be expressed in the general forms as
PPT Slide
Lager Image
where k = A for the relative motion with the reduced mass M A = μ , k = B for the vibrations χ 1 15 with the corresponding reduced mass M B = μ i , k = C for the bending modes φ 1 - φ 24 and the moment of inertia M C = I j . For the incident molecule, k = χ, η and η ' are associated with the reduced mass μ inc and the moment of inertia I inc . We use the standard numerical routines 47 , 48 to integrate these equations for the initial conditions at t = to and their conjugate quantities d r ( to )/d t , d χi ( to )/d t , d φ j ( to )/d t , d χ ( to )/d t , d η ( to )/d t , and d η '( to )/ d t , where the derivatives are evaluated at t = to . The initial conditions for the relative and vibrational motions in the interaction zone are given in Ref. 29. We sample 40,000 trajectories for each run at 300 K, where the sampling includes determining collision energies ( E ) sampled from the Maxwell distribution and weighting the initial vibrational and rotational energies by the Boltzmann distribution at 300 K.
Potential and spectroscopic constants for toluene, H2and D2
PPT Slide
Lager Image
aReference 44; the ionization energy and polarizability of toluene are shown here for easy comparison with the corresponding values of H2 and D2. breference 45. creference 46. drange parameter is determined from the relation bi = (2Dii)1/2/ωi.
The principal quantity we determine in non-dissociative cases is the difference between the initial and final energies of translation, rotation and vibration of H 2 (or D 2 ). We denote the ensemble-averaged energy transfer by <Δ E > = < Efinal Einitial > such that a positive value represents an energy loss by toluene via the V-T , V-V and V-R pathways. In studying dissociative cases, we determine the time evolution of methyl and ring C-H bonds, and sum all trajectories, which lead to the methyl or ring C-H energy reaching its dissociation threshold. We discuss the latter process in terms of collision-induced intramolecualr energy flow.
Results and Discussion
Energy Loss by Toluene. In Figure 2(a) , we show the dependence of energy loss by toluene on the total vibrational energy of toluene ET in the toluene + H 2 collision at 300 K. Plotted is the ET dependence of the amount of energy transfer from toluene to H 2 via V-T , V-V , V-R pathways, as well as their sum. The lowest value of ET considered in this figure is 6204 cm −1 or 0.769 eV, which is equivalent to the sum of one-quantum vibrational energies of C-H methyl and C-H ring 3060 cm −1 (0.379 eV) and 3144 cm −1 (0.390 eV), respectively. We treat the incident molecules to be in the ground state, so the initial energies are ½(4401) and ½(3115) cm −1 for H 2 and D 2 , respectively. (We give energies in both cm −1 and eV for convenience.) At the lowest value of ET , the energy loss is −6 cm −1 [see “sum” in Figure 2(a) ], which is obtained by adding 13, −26 and 7 cm −1 for V-T , V-V and V-R , respectively. When ET is raised from 6204 cm −1 to 37,221 cm −1 , the range considered in experimental studies, 23 the amount of energy transfer varies from −6 cm −1 to 123 cm −1 , which is less than 1% of ET . At the risk of repetition: positive energy, such as 123 cm −1 , represents the energy loss by toluene. The maximum energy content 37,221 cm −1 considered in Figure 2(a) is the sum of C-H methyl and C-H ring energies, which are taken to be 18,359 cm −1 (2.276 eV) and 18,862 cm −1 (2.339 eV) respectively. The experimental data plotted in Figure 2(a) are the work of Toselli et al. obtained from the measurements of collisional loss of vibrational energy by use of the time-resolved infrared fluorescence from the C-H stretch modes near 3.3 μm. 23 The energy loss gradually increases with increasing ET as found in the computed sum. The latter variation is known to be a general behavior for the relaxation of toluene by rare gases and diatomic molecules in both experimental and theoretical studies. 16 , 17
PPT Slide
Lager Image
Plots of energy loss by toluene vs. total energy of toluene: (a) Toluene + H2 and (b) toluene + D2. In each frame, the transfer of vibrational energy from toluene to various motions of H2/D2 (i.e, V-T, V-R and V-V) is indicated. The sum of all three contributions indicated by “sum” is compared with the experiment data.23
A particularly important aspect of the present result is that we can identify the contributions of individual energy transfer pathways to the vibrational relaxation of toluene. As shown in Figure 2(a) , the V-V values are negative except near the upper end of ET . The V-V curve rises from −26 cm −1 at ET = 6204 cm −1 to −3 cm −1 at the toluene energy content as high as ET = 31,018 cm −1 . As the negative values indicate toluene gains energy, the V-V pathway does not contribute to the relaxation of excited toluene. On the other hand, the amount of V-T energy transfer is positive and significant over the entire ET range. In fact, H 2 removes energy from toluene by 13 cm −1 at ET = 6204 cm −1 to 81 cm −1 at ET = 37,221 cm −1 , the V-T process being the leading pathway for the relaxation process. We note that the V-R value shown in Figure 2(a) is always positive, but its contribution to the sum is minor.
The plots similar to Figure 2(a) are presented in Figure 2(b) for toluene + D 2 . The collision system now involves a heavier incident molecule with a lower force constant. The reduced mass affecting the relative motion of the collision system is now μ μ D2 which is twice the reduced mass μ μ H2 of the toluene + H 2 system. The vibrational frequency of D 2 is 3115 cm −1 , which is very close to the C-H frequencies (C-H methyl = 3144 cm −1 and C-H ring = 3060 cm −1 ). Such near resonant condition enhances V-V energy flow. Thus we find the V-V pathway is more efficient than in toluene + H 2 . Unlike the toluene + H 2 case, the amount of V-V energy transfer is now always positive; see Figure 2(b) . It increases from 13 cm −1 at ET = 6204 cm −1 to 67 cm −1 at ET = 37,221 cm −1 , which are comparable to the V-T values. The contribution of V-R pathway is insignificant as in the toluene + H 2 system. As shown in Figure 2(b) , the magnitude of energy loss indicated by “sum” is larger than the observed values 23 but the ET dependence is in fair agreement. The comparison of Figures 2(a) and 2(b) indicates that the heavier D 2 with a frequency in near resonance with the C-H bonds is more efficient in relaxing the highly vibrationally excited toluene. Thus increased efficiency of V-T and V-V energy transfer in toluene + D 2 is the manifestation of the effects of mass and frequency.
Figure 3 shows the effects of vibrational excitation of H 2 / D 2 on the energy loss of toluene by assigning the energies equivalent to one and two quanta to H 2 /D 2 ; see v = 1 and 2 curves in Figures 3(a) and 3(b) . The v = 0 results are the “sum” curves reproduced from Figures 2(a) and 2(b) . The v = 1 and 2 curves of H 2 clearly represent the incident molecule losing its energy to the target. The extent of H 2 to toluene energy transfer increases with increasing the vibrational excitation of H 2 , which is mainly due to the contribution of the V-V energy pathway working in the direction of H 2 to toluene. A similar trend of energy transfer decreasing with increasing vibrational excitation is seen in Figure 3(b) for toluene + D 2 but the process now transfers energy from toluene to D 2 except at the lower end of ET .
PPT Slide
Lager Image
Dependence of “sum” on the vibrational excitation (v = 0 – 2) of the incident molecules: (a) toluene + H2 and (b) toluene + D2.
Dynamics of Bond Dissociation. When the internal state of a highly vibrationally excited toluene is perturbed by H 2 / D 2 , either C-H methyl or C-H ring bond can dissociate. We discuss important characteristics of the dissociation dynamics in this section. The results presented in Figure 2 suggest collisioninduced bond dissociation is not likely to occur for toluene with low to medium extent of vibrational excitation. Thus we take a sufficiently high vibrational energy state of toluene. For the representative trajectory shown in Figure 4 , the CH methyl and C-H ring bonds start out with energies 3.472 eV and 4.457 eV, respectively, 0.50 eV below the dissociation threshold, while other bonds and bends are initially in the ground state. In Figure 4(a) , we show the time evolution of the toluene-H 2 distance, which is the collision trajectory, the C-H methyl and C-H ring distances.
PPT Slide
Lager Image
The dynamics of a representative trajectory of the CHmethyl bond dissociation in toluene + H2 collision: (a) Time development of toluene-H2, C-Hmethyl and C-Hring distances, (b) Time development of C-Hmethyl and C-Hring vibrational energies. The total energy of toluene above the ground state ET is 63,944 cm−1 (7.928 eV), which is distributed between C-Hmethyl and CHring such that the energy of each bond is 0.50 eV below the dissociation threshold. The collision energy E taken from the Maxwell sampling is 0.070 eV.
The incident H 2 molecule is in the ground state, so it acts primarily as a perturber inducing energy flow from one C-H bond to another. The collision time is scaled such that the instant of the initial impact is zero. For the representative trajectory chosen, H 2 approaches toluene to the closest separation 4.30 Å with the collision energy E = 0.070 eV. The C-H methyl bond begins to be perturbed at t ≈ −0.1 ps, when energy starts to flow from the C-H ring bond; see Figure 4(b) . During this subpicosecond time period, main part of the energy lost by the C-H ring bond passes through the χ 3 and χ 2 C-C bonds and eventually localizes in the C-H methyl bond for its dissociation. As shown in Figure 4(b) , the amount of energy lost by the C-H ring bond is 0.654 eV, whereas the CH methyl bond gains 0.650 eV, thus exceeding the dissociation threshold. The closeness of the latter two values indicates the energy needed for the dissociation of highly excited CH methyl comes entirely from C-H ring through intramolecular flow.
Figure 5(a) shows the dynamics of energy flow in toluene + D 2 . Intramolecular energy flow and bond dissociation occurs in a time-scale of about 0.6 ps, much longer than 0.1 ps of the toluene + H 2 case. As shown in Figure 5 , the excited C-H methyl bond undergoes nearly 20 oscillations before dissociation. It is interesting that such a highly excited CH methyl bond can remain undissociated so long. The bond waits that long until the last trace of energy needed for dissociation arrives from C-H ring . The detail of energy flow is seen in Figure 5(b) . The C-Hmethyl bond distance begins to be perturbed from its initial vibrational motion as the collision trajectory approaches t = 0, when the first impact occurs. The incident molecule remains near toluene for ~0.60 ps, during which the C-H methyl bond first loses a small mount of energy to D 2 and then gains a large amount from C-H ring near the end of collision via intramolecular flow. During the period of ~0.6 ps, the C-H ring bond loses energy as high as ~1.5 eV, which flows through C-C bonds. Only about 0.6 eV of the energy reaches the C-H methyl bond, whereas the rest accumulates in C-C stretches and H-C-C, C-C-C bends. Thus the time evolution suggests that it takes about 0.6 ps for the vibrational energy to travel a distance of 2.8 Å for two C-C bonds for bond dissociation in the toluene + D 2 collision. Comparing the curves shown in Figure 4 with those in Figure 5 , we can see the dissociation of C-H methyl in toluene + D 2 occurs via a long-time, or complex-mode, mechanism. Although it is short, the time scale for bond dissociation in the representative case of toluene + D 2 shows that the dissociation does not occur instantaneously.
As shown in Figures 6 and 7 , the collision dynamics leading to C-H ring dissociation in both toluene + H 2 and toluene + D 2 , bond dissociation proceeds through a complexmode period. In the latter system, the C-H ring bond distance undergoes a large amplitude motion before dissociation; see Figure 7(a) . Even in toluene + H 2 , the period is nearly 0.3 ps, which is significantly longer than the C-H methyl dissociation case shown in Figure 4(a) . As the collision partners interact in the representative trajectory considered in Figure 6(b) , both C-H methyl and C-H ring bonds begin to lose energy mainly to C-C bonds and bends intramoleculary, but near the end of complex-mode period, the C-H ring bond gains back most of its energy from C-C bonds and bends, as well as that from CH methyl . The collision trajectory shown in Figure 6(a) indicates that the benzyl radical and H 2 still undergo another impact as Hring flies away from the primary zone. In the C-H ring dissociation for toluene + D 2 considered in Figure 7(b) , the CH methyl loses its vibrational energy as large as 1.38 eV, part of which flows to C-H ring upon the first impact. The C-H ring bond gains enough energy for dissociation, but the fragmented H ring becomes attracted back to the radical before it has a chance to escape from the primary zone. The large amplitude motion of the C-H ring bond noted above is the manifestation of this recapture period. At the latter period, the rebounding H ring atom gains enough kinetic energy to recede from the reaction zone. We note that in the four representative cases considered in Figures 4 - 7 , the collision energies E are 0.070, 0.035, 0.087 and 0.052 eV, respectively, which lie in the tail of the Maxwell distribution at 300 K.
PPT Slide
Lager Image
Plots are the same as Figure 4, but they are now for the toluene + D2 collision. The collision energy E is 0.035 eV.
PPT Slide
Lager Image
The dynamics of a representative trajectory of the CHring bond dissociation in toluene + H2 collision: (a) Time development of toluene-H2, C-Hring and C-Hmethyl distances. (b) Time development of C-Hring and C-Hmethyl vibrational energies. The total energy of toluene ET is the same as that given in Figure 4, but the collision energy E is 0.087 eV.
We find the collisions taking place at or below the most probable velocity (or energy) of the distribution are inefficient in dissociating even for such highly excited C-H bonds. Thus H 2 /D 2 approaching toluene with collision energy above the most probable value has a chance to induce intramolecular energy flow to a sufficient extent for bond dissociation.
We plot the dissociation probabilities for toluene + H 2 collision at 300 K in Figure 8(a) . The total energy content of highly excited toluene is considered to be localized initially in the C-H bonds as noted above. The initial energy of each bond is maintained below its dissociation threshold value, so that intramolecular energy flow through C-C bonds and bends can lead to dissociation. Dissociation probabilities are found to be low as the light molecules are not efficient in inducing intramolecular energy flow from one C-H bond to another. Bond dissociations begin to occur when the total energy content of toluene is gradually raised to 57,492 cm −1 , which is partitioned 32,719 cm −1 (4.0567 eV) in C-H ring and 24,773 cm −1 (3.0715 eV) in C-H methyl . The latter two energies are ~0.9 eV below the dissociation threshold. The semilogarithmic plot shown in Figure 8(a) indicates the dissociation probability is low (~10 −5 ), but rise very rapidly with increasing energy content. At 70,396 cm −1 , which is partitioned in the two bonds such that their bond energies are 0.1 eV below the threshold, both probabilities are now ~0.01.
PPT Slide
Lager Image
Plots are the same as Figure 6, but they are now for the toluene + D2 collision. The collision energy E is 0.052 eV.
As shown in Figure 8(a) , the probability of C-H methyl dissociation is higher than that of C-H ring dissociation over the entire energy range considered. A similar energy dependence of both probabilities is seen for toluene + D 2 ; see Figure 8(b) . The D 2 probabilities are slightly higher than the H 2 case, but the principal qualitative features of the semilogarithmic energy dependence remain unchanged. It is interesting to note that the probabilities of toluene + N 2 /O 2 are higher than the present values primarily due to the fact that heavier incident molecules lead to larger extent of intramolecular energy flow. 22
PPT Slide
Lager Image
Dissociation probabilities for C-Hmethyl and C-Hring bonds in (a) toluene + H2 and (b) toluene + D2. Both H2 and D2 are in the ground state.
Concluding Comments
We have studied energy loss by vibrationally excited toluene and dissociation of C-H bonds in the toluene + H 2 / D 2 collision systems at 300 K using classical trajectory procedures. The collision system consists of the primary zone of toluene + H 2 /D 2 interaction, where the incident molecule interacts with both C-H methyl and C-H ring bonds, and the secondary zone which includes the stretches and bends of toluene beyond the primary zone. Trajectory calculations are carried out for the target molecule with vibrational energy in its two C-H bonds varying from 5000 to 40,000 cm −1 , while other stretches and bends are initially in the ground state.
The amount of energy lost by vibrationally excited toluene interacting with the ground-state H 2 /D 2 is small, but the energy loss increases when the extent of vibrational excitation increases. The entire energy transfer process taking place on a subpicosecond time scale. The dependence and magnitude of energy loss on the total energy of toluene are in general agreement with experimental data. The amount of energy transfer is significantly larger in toluene + D 2 , where C-H and D 2 vibrations are in near resonance. Further, the heavier D 2 imparts stronger perturbation on the C-H bonds, thus enhancing the relaxation process of toluene. The main contribution to the vibrational relaxation of toluene comes from the V-T energy transfer pathway in both toluene + H 2 and toluene + D 2 collisions. The V-V pathway is inefficient in toluene + H 2 because of the large disparity in their vibrational frequencies. In general, the V-R pathway plays a minor role in the relaxation of toluene. When the incident molecules H 2 /D 2 are vibrationally excited ( v = 1 and 2), the amount of energy lost by toluene is small compared to the v = 0 case as toluene tends to gain energy from the excited collision partners.
When the total energy content ET of toluene is sufficiently high, either C-H bond can dissociate in collisions with H 2 /D 2 molecule. The time evolution of collision events shows that the dissociation occurs when the internal energy of toluene is initially above 55,000 cm −1 (~6.8 eV) and is nearly equally distributed between C-H methyl and C-H ring . Dissociation probabilities are low but rise exponentially with increasing initial energy. Dissociation occurs as a result of the incident molecule (H 2 or D 2 ) inducing intramolecular energy flow from one C-H bond to another via C-C bonds and bends. The dissociation probability of C-H methyl is significantly higher than that of C-H ring as collision-induced energy flow from the high-frequency C-H ring to low-frequency C-H methyl dominates the reverse process. No experimental data are available for C-H bond dissociation of highly excited toluene, but the present values are comparable with those for toluene + N 2 and toluene + O 2 collisions, 22 while they are lower than those of toluene + HF and toluene + Ar. 30 , 31
This work was financially supported by research fund of Chonnam National University in 2012.
Cottrell T. L. , McCoubrey J. C. 1961 Molecular Energy Transfer in Gases Butterworths London
Burnett G. M. , North A. M. 1969 Transfer and Storage of Energy by Molecules Wiley New York
Smith I. W. M. 1977 Gas Kinetics and Energy Transfer; Specialist Periodical Reports 2 1 - 57
Ma J. , Liu P. , Zhang M. , Dai H.-L. 2005 J. Chem. Phys. 123 154306 -    DOI : 10.1063/1.2049271
Poel K. L. , Alwahabi Z. T. , King K. D. 1995 Chem. Phys. 201 263 -    DOI : 10.1016/0301-0104(95)00257-6
Smith I. W. M. , Warr J. F. 1991 J. Chem. Soc. Faraday Trans. 87 807 -    DOI : 10.1039/ft9918700807
Hippler H. , Otto B. , Troe J. 1989 Ber. BunsenGes. Phys. Chem. 93 428 -    DOI : 10.1002/bbpc.19890930404
Nakashima N. , Yoshihara K. 1983 J. Chem. Phys. 79 2727 -    DOI : 10.1063/1.446176
Abel B. , Herzog B. , Hippler H. , Troe J. 1989 J. Chem. Phys. 91 900 -    DOI : 10.1063/1.457141
Tam A. C. , Sontag H. , Hess P. 1985 Chem. Phys. Lett. 120 280 -    DOI : 10.1016/0009-2614(85)87058-5
Wallington T. J. , Scheer M. D. , Braun W. 1987 Chem. Phys. Lett. 138 538 -    DOI : 10.1016/0009-2614(87)80120-3
Beck K. M. , Gordon R. J. 1987 J. Chem. Phys. 87 5681 -    DOI : 10.1063/1.453736
Hsu H. C. , Tsai M.-T. , Dyakov Y. A. , Ni C.-K. 2012 Int. Rev. Phys. Chem. 31 201 -    DOI : 10.1080/0144235X.2012.673282
Gilbert R. G. 1984 J. Chem. Phys. 80 5501 -    DOI : 10.1063/1.446661
Toselli B. M. , Barker J. R. 1990 Chem. Phys. Lett. 174 304 -    DOI : 10.1016/0009-2614(90)85350-L
Lim K. F. 1994 J. Chem. Phys. 101 8756 -    DOI : 10.1063/1.468070
Wright S. M. A. , Sims I. R. , Smith I. W. M. 2000 J. Phys. Chem. A 104 10347 -    DOI : 10.1021/jp0014216
Kable S. H. , Knight A. E. W. 2003 J Phys. Chem. A 107 10813 -    DOI : 10.1021/jp035516u
Zhang J. Z. H. , Li Y. M. , Wang M. L. , Xiang Y. 2004 Adv. Ser. Phys. Chem. 14 209 -
Fu B. , Han Y.-C. , Bowman J. M. , Leonori F. , Balucani N. , Angelucci L. , Occhiogrosso A. , Petrucci R. , Casavecchia P. 2012 J. Chem. Phys. 137 22A532 -
Dagdigian P. J. 2013 Int. Rev. Phys. Chem. 32 229 -    DOI : 10.1080/0144235X.2012.758543
Ree J. , Kim S. H. , Lee S. K. J. 2013 Bull. Kor. Chem. Soc. 34 1494 -    DOI : 10.5012/bkcs.2013.34.5.1494
Toselli B. M. , Brenner J. D. , Yerram M. L. , Chin W. E. , King K. D. , Barker J. R. 1991 J. Chem. Phys. 95 176 -    DOI : 10.1063/1.461473
Damm M. , Deckert F. , Hippler H. , Troe J. 1991 J. Phys. Chem. 95 2005 -    DOI : 10.1021/j100158a022
Lenzer T. , Luther K. , Troe J. , Gilbert R. G. , Lim K. F. 1995 J. Chem. Phys. 103 626 -    DOI : 10.1063/1.470096
Bernshtein V. , Lim K. F. , Oref I. 1995 J. Phys. Chem. 99 4531 -    DOI : 10.1021/j100013a024
Shin H. K. 2000 J. Phys. Chem. A 104 6699 -    DOI : 10.1021/jp0009094
Elioff M. S. , Fang M. , Mullin A. S. 2001 J. Chem. Phys. 115 6990 -    DOI : 10.1063/1.1400782
Ree J. , Kim Y. H. , Shin H. K. 2002 J. Chem. Phys. 116 4858 -    DOI : 10.1063/1.1452726
Ree J. , Kim Y. H. , Shin H. K. 2004 Chem. Phys. Lett. 394 250 -    DOI : 10.1016/j.cplett.2004.06.134
Ree J. , Kim S. H. , Lee T. H. , Kim Y. H. 2006 Bull. Kor. Chem. Soc. 27 495 -    DOI : 10.5012/bkcs.2006.27.4.495
Berkowitz J. , Ellison G. B. , Gutman D. 1994 J. Phys. Chem. 98 2744 -    DOI : 10.1021/j100062a009
Weast R. C. 1986 CRC Handbook of Chemistry and Physics 67th ed CRC Press Boca Raton, FL
Barone V. , Fliszar S. 1996 J. Mol. Struc. (Theochem) 369 29 -    DOI : 10.1016/S0166-1280(96)04581-2
Guan Y. , Thompson D. L. 1990 J. Chem. Phys. 92 313 -    DOI : 10.1063/1.458601
Lapouge C. , Cavagnat D. 1998 J. Phys. Chem. A 102 8393 -    DOI : 10.1021/jp9821230
Kjaergaard H. G. , Turbull D. M. , Henry B. R. 1997 J. Phys. Chem. A 101 2589 -    DOI : 10.1021/jp963100t
Kjaergaard H. G. , Turnbull D. M. , Henry B. R. 1998 J. Phys. Chem. A 102 6095 -    DOI : 10.1021/jp981411t
Lim K. F. 1994 J. Chem. Phys. 100 7385 -    DOI : 10.1063/1.466882
Hirschfelder J. O. , Curtiss C. F. , Bird R. B. 1967 Molecular Theory of Gases and Liquids Wiley New York
Xie Y. , Boggs J. E. 1986 J. Comp. Chem. 7 158 -    DOI : 10.1002/jcc.540070209
Olney T. N. , Cann N. M. , Cooper G. , Brion C. E. 1997 Chem. Phys. 223 59 -    DOI : 10.1016/S0301-0104(97)00145-6
Huber K. P. , Herzberg G. 1979 Constants of Diatomic Molecules Van Nostrand Reinhold New York
Gear C. W. 1971 Numerical Initial Value Problems in Ordinary Differential Equations Prentice-Hall New York
1989 MATH/LIBRARY, Fortran Subroutines for Mathematical Applications IMSL Houston