X. Zhou et al. / Journal of Fluorine Chemistry 153 (2013) 101–106
103
Table 1
coefficient for the BTP pyrolysis, expressed by the Arrhenius
equation as follows:
The fitted equations and thermodynamic data calculated in Fig.5.
Mole
Fitted equation
R2
Apparent
Ea (kJ molꢀ1
)
ꢂ
ꢃ
Ea
fraction
k ¼ A exp
(2)
RT
10%
20%
50%
ln ln[1/(1 ꢀ
ln ln[1/(1 ꢀ
ln ln[1/(1 ꢀ
a
a
a
)] = 36.87–37.37 ꢂ 103/T
)] = 39.36–39.50 ꢂ 103/T
)] = 37.52–37.79 ꢂ 103/T
0.9551
0.9850
0.9908
310.7
328.4
314.2
where Ea (kJ molꢀ1) is the activation energy, R = 8.3145 J molꢀ1 Kꢀ1
is the ideal gas constant, A is the pre-exponential factor and T (K) is
the absolute temperature.
Eq. (3) was obtained by introducing Eq. (2) into Eq. (1). Note
that, for first order chemical reaction, the thermal decomposition
ratio is irrelevant to the initial concentration of the reactant. And in
a narrow temperature range, there must be a linear correlation
Fig. 5. The same goes for the possible transition states. Beyond that,
the bond distances and bond dissociation energies (BDEs) for BTP
molecular associated with energy barriers for three feasible
decomposition pathways, which were listed in Table 2, were also
investigated at the same basis set and theory level. All BDEs and
energy barriers involved in the dissociation and decomposition
reactions of BTP are expressed by thermal enthalpies that
corrected by zero point energy, based upon frequency calculation
(freq = temperature = 873.15) by Gaussian 03.
between ln ln[1/(1 ꢀ
a
)] and 1/T with a fixed residence time (see
Eq. (4)).
ꢂ
ꢃ
ꢀEa
lnð1 ꢀ
a
Þ ¼ ꢀA exp
t
(3)
(4)
RT
ꢀ
ꢁ
1
Ea
From the theoretical calculation results we can see that the
dissociation energy of C-Br bond (308.2 kJ molꢀ1) is the lowest
among all the bonds in the BTP molecule. The long bond distance
and low dissociation energy mean that the bond is loosely attached
and easy to break, and Brꢁ radical can be easily formed at a
relatively lower energy. Besides, bond dissociation energy of C–C
bond (386.3 kJ molꢀ1) is much lower than BDE values of C–F, C–H
and C55C bonds except for C–Br bond. The conclusion draw from
ln ln
¼ ꢀ
þ constant
1 ꢀ
a
RT
Based on the experimental data of the three different initial
mole fraction experiments with a fixed residence time of 5 s, the
calculation results of Eq. (4) are plotted in Fig. 4. From the plot it
can be observed that the relationship between ln ln[1/(1 ꢀ
a)] and
1000/T is linear, and fitted well. So, it is reasonable to assume the
thermal degradation of BTP to be the first order reaction from
above calculation and Fig. 3. The apparent activation energies
calculated are 314.2 kJ molꢀ1, 328.4 kJ molꢀ1, 310.7 kJ molꢀ1 with
the relative error within 6% and the results are categorized in
Table 1.
above is that BTP molecule may provide CF3 and Brꢁ radicals at
ꢁ
high temperature, which are very effective in extinguishing a fire.
In addition, the energy barriers of three possible decomposition
pathways were also calculated using TS method, among which the
pathway 1, CF3CBrCH2–TS1–CF3CCH + HBr, was calculated to have
the lowest energy barrier of 300.1 kJ molꢀ1. Another pathway that
leads to product channel, CF2 + FBrCH2, has an energy barrier of
437.1 kJ molꢀ1. The calculation results indicate that the most
probable dissociation products should be (CF3CCH + HBr) for the
lowest energy barrier. Likewise, (CF3Br + CCH2) with energy barrier
of 502.0 kJ molꢀ1 should be the most difficult products to produce.
8
2.3. Theoretical analysis by Gaussian 03
In order to reveal the thermal decomposition of BTP and its fire
extinguishing mechanisms, theoretical analysis by Gaussian 03
was performed. Density function theory (DFT) method accompa-
nied with 6-311++G(d,p) basis set was employed to do this. During
the calculation, geometrical optimization and frequency analysis
for BTP, including feasible decomposition products or fragments,
and transition states, were probed theoretically with Gaussian 03
at the B3LYP/6-311++G(d,p) theory level. Besides, connections
between reactants, transition structures and products were
confirmed by intrinsic reaction coordinate (IRC) calculation. Part
of the calculation results, such as molecular structure of BTP and
corresponding optimized geometrical parameters, are depicted in
>
>
>
300:1 kJ molꢀ1
ꢀ!
CF3CCH þ HBr
CH2CFBr þ CF2
CF3Br þ CCH2
>
>
TS1
>
>
<
437:1 kJ molꢀ1
CF3CBr¼CH2
ꢀ!
>
>
TS2
>
>
>
502:0 kJ molꢀ1
>
>
:
ꢀ!
TS3
To verify the above theoretical calculation results and investi-
gate the decomposition products of BTP, the gaseous products
were also detected by GC–MS. Analysis of the obtained mass
spectrum reveals two major decomposition products. One with the
characteristic fragments at m/z = 69, 75 and 94 (molecular ion
Table 2
BDEs and energy barriers calculated by Gaussian 03.
Bonds/
Bond
Reaction equations
BDEs/energy
barriers
(kJ molꢀ1
)
pathways
distance
˚
(A)
C–Br
1.906
1.511
1.325
1.350
1.348
1.083
1.081
—
CF3CBrCH2–CF3(Cꢁ)CH2 + Brꢃ
308.2a
386.3a
727.9a
447.7a
479.8a
464.9a
471.5a
300.1b
437.1b
502.0b
C–C
CF3CBrCH2–CF3ꢁ + CH2(Cꢁ)Br
ꢁ
CF3CBrCH2–CF3(Cꢁ)Br + CH2
C5C
(C–F)1
CF3CBrCH2–CH2CBrCF2ꢁ + Fꢁ
CF3CBrCH2–CH2CBrCF2ꢁ + Fꢁ
CF3CBrCH2–CF3CBrCHꢁ + Hꢁ
CF3CBrCH2–CF3CBrCHꢁ + Hꢁ
CF3CBrCH2–TS1–CF3CCH + HBr
CF3CBrCH2–TS2–CF2ꢁ + CFBrCH2
CF3CBrCH2–TS3–CF3Br + CCH2
(C–F)2
(C–H)1
(C–H)2
Pathway 1
Pathway 2
Pathway 3
—
—
a
Bond dissociation energy (BDE) values of different bonds.
b
Fig. 4. Relationship between ln ln[1/(1 ꢀ
a)] and 1000/T.
Reaction energy barriers of three feasible decomposition pathways for BTP.