R. Wang et al. / Catalysis Today 149 (2010) 117–121
121
Further combination of the above two equations leads to
in the presence of the H3PW6Mo6O40 catalyst. Hence, the macro-
kinetic results derived from this study should lay a solid foundation
for better understanding and optimum design of the overall process.
’
1 þ
nA;2
¼
nA0
(8)
(9)
’
where
4. Conclusions
nac
’
¼ K
nn
The catalytic performance of H3PWxMo12ꢀxO40 (x = 1, 3, 6) was
found to follow the order of H3PWMo11O40 < H3PW3Mo9O40 <
-
o
The constant K is slight dependant on temperature within the
temperature range of the experiment conditions and can therefore
be regarded as fixed. As for nac and nn-o, the amounts of acetonitrile
H3PW6Mo6O40. Moreover, H3PW6Mo6O40 was also found to be
superior to H3PW12O40. The transformation of phosphotungstic acid
to its cesium salt sharply weakens the catalytic effect, as verified by
and n-octane are identical for all the experiments. Hence
w can be
Cs2.5H0.5PW12O40 < H3PW12O40. Under the simulated system with
regarded as an identical constant for all the experiments.
an initial sulfur concentration of 500 ppmw, the ODS process
catalyzed by H3PW6Mo6O40 can proceed very efficiently at 60 8C
and an O/S molar ratio of 15, achieving a desulfurization efficiency of
99.79% or higher after 60 min. Hence, H3PW6Mo6O40 is promising for
thedeepdesulfurizationofactualODSprocess.Theprocessefficiency
can be improved by an increase in any of the factors including
temperature, O/S molar ratio, and catalyst dosage, whereby the
favourable operating conditions were recommended as T = 60 8C, O/
S = 15, and a catalyst dosage of 6.93 g (H3PW6Mo6O40)/L (simulated
diesel). The effects of both the O/S molar ratio and the initial sulfur
concentration conform to the law of mass action. With the aid of GC–
MS analysis, sulfone species was confirmed to be the only product
after reaction for 150 min. Furthermore, macro-kinetics of the
process catalyzed by H3PW6Mo6O40 was studied, from which the
reactionorderswerefoundtobe1.02toDBTand0.38toH2O2,andthe
activation energy of the reaction was found to be 43.3 kJ/mol.
Moreover, the used catalyst can berecoveredand reused withalmost
the same activity as the fresh.
Owing to the equal volume of n-octane and acetonitrile used in
all the experiments, according to Eq. (8), the concentration of DBT
in acetonitrile phase can be written as
’
1 þ
CA;2
¼
CA0
(10)
’
where CA0 is the initial concentration of DBT in n-octane phase.
From the experimental test, the value of was determined as 1.17.
w
In the presence of H3PW6Mo6O40 and H2O2 (B) in acetonitrile
phase, the irreversible oxidation of DBT will proceed at the initial
rate of
ri ¼ kCAm;2CBn
(11)
where the kinetic factor k can be expressed in terms of activation
energy E and temperature T as the well-known Arrhenius equation
ꢀ
ꢁ
E
k ¼ A exp
ꢀ
(12)
(13)
(14)
RT
Acknowledgement
Combining Eqs. (10)–(12) yields
ꢀ
ꢁ
E
ri ¼ A0 exp ꢀ
CAm0CBn
Financial support from the Ministry of Education of China under
New Century Excellent Talents Project (NCET-05-0584) is grate-
fully acknowledged.
RT
where
ꢀ
ꢁ
m
References
’
1 þ
A0 ¼ A
’
[1] E.W. Qian, J. Jpn. Petrol. Inst. 51 (2008) 14.
[2] C. Song, Catal. Today 86 (2003) 211.
[3] A. Ro¨thlishberger, R. Prins, J. Catal. 235 (2005) 229.
[4] I.V. Babich, J.A. Moulijn, Fuel 82 (2003) 607.
In the light of above assumptions, since the steady state of the
overall process can be achieved instantaneously, the initial rate of
the overall process can be represented by Eq. (13).
The value of ri in Eq. (13) can be calculated as the mean slope of
the initial stage of the concentration–time curve of DBT in n-
octane. It should be noted that, owing to the instant attainment of
equilibrium through mass transfer, the starting point at time zero
for determining the mean slope of the above curve corresponds to a
[5] R.T. Yang, A.J. Herna´ndez-Maldonado, F.H. Yang, Science 301 (2003) 79.
[6] S. Mondal, Y. Hangun-Balkir, L. Alexandrova, D. Link, B. Howard, P. Zandhuis, A.
Cugini, C.P. Horwitz, T.J. Collins, Catal. Today 116 (2006) 554.
[7] S. Liu, B. Wang, B. Cui, L. Sun, Fuel 87 (2008) 422.
[8] L. Cedeno-Caero, H. Gomez-Bernal, A. Fraustro-Cuevas, H.D. Guerra-Gomez, R.
Cuevas-Garcia, Catal. Today 133–135 (2008) 244.
[9] A. Chica, A. Corma, M.E. Domine, J. Catal. 242 (2006) 299.
[10] J.L. Garcia-Gutierrez, G.A. Fuentes, M.E. Hernandez-Teran, P. Garcia, F. Murrieta-
Guevara, F. Jimenez-Cruz, Appl. Catal. A: Gen. 334 (2008) 366.
[11] H. Lu¨, J. Gao, Z. Jiang, F. Jing, Y. Yang, G. Wang, C. Li, J. Catal. 239 (2006) 369.
[12] G. Yu, S. Lu, H. Chen, Z. Zhu, Carbon 43 (2005) 2285.
[13] K. Yazu, Y. Yamamoto, T. Furuya, K. Miki, K. Ukegawa, Energy Fuels 15 (2001) 1535.
[14] J.L. Garcia-Gutierrez, G.A. Fuentes, M.E. Hernandez-Teran, F. Murrieta, J. Navar-
rete, F. Jimenez-Cruz, Appl. Catal. A: Gen. 305 (2006) 15.
DBT concentration of CA0–CA,2, rather than CA0
.
In the above experimental results, by means of variable
isolation, i.e. leaving only one changeable variable while fixing
the rest, the pattern of the effects of each variable has been
obtained. Based on these results, setting Eq. (13) as objective
function, through successive least-squares regressions of ln ri
against CA,2, ln ri against CA,2, and further ln ri against 1/T, the
values of m, n, E in Eq. (13) were determined as follows:
[15] F. Al-Shahrani, T. Xiao, S.A. Liewellyn, S. Barri, Z. Jiang, H. Shi, G. Martinie, M.L.H.
Green, Appl. Catal. B: Environ. 73 (2007) 311.
[16] B. Zapata, F. Pedraza, M.A. Valenzuela, Catal. Today 106 (2005) 219.
[17] F.M. Collins, A.R. Lucy, C. Sharp, J. Mol. Catal. A: Chem. 117 (1997) 397.
[18] M.T.C. Fairbridge, Z. Ring, Appl. Catal. A: Gen. 219 (2001) 267.
[19] C. Li, J. Gao, Z. Jiang, S. Wang, H. Lu, Y. Yang, F. Jing, Top. Catal. 35 (2005) 169.
[20] K. Yazu, T. Furuya, K. Miki, J. Jpn. Petrol. Inst. 46 (2003) 379.
[21] C. Li, Z. Jiang, J. Gao, Y. Yang, S. Wang, F. Tian, F. Sun, X. Sun, P. Ying, C. Han, Chem.
Eur. J. 10 (2004) 2277.
[22] I.V. Kozhevnikov, Catalysts for Fine Chemical Synthesis, vol. 2: Catalysis by
Polyoxometalates, John Wiley & Sons., Chichester, 2002.
[23] M.T. Pope, Heteropoly and Isopoly Oxometalates, Springer, Berlin, 1983.
[24] K. Nakamoto, Infrared Spectra of Inorganic and Coordination Compounds, John
Wiley & Sons, London, 1963.
m ¼ 1:02 ðR ¼ 0:996Þ
n ¼ 0:38 ðR ¼ 0:989Þ
E ¼ 43:3 kJ=mol ðR ¼ 0:992Þ
Meanwhile, the value of A0 was determined as
(1.01 ꢁ 0.06) ꢂ 107 (mol/L)ꢀ0.40 minꢀ1
.
Hence, the H3PW6Mo6O40 catalyzed oxidation of DBT by H2O2 is
1st order to DBT and 0.4th order to H2O2. Therefore, it can be drawn
fromthisstudythat,whenhighO/Smolarratioareused, thereaction
can simply be taken as a pseudo-1st order reaction. The activation
energy of 43.3 kJ/mol indicates that the reaction process is quite fast
[25] O.A. Kholdeeva, L.A. Kovaleva, R.I. Maksimovskaya, G.M. Maksimov, J. Mol. Catal.
A: Chem. 158 (2000) 223.
[26] C. Aubry, G. Chottard, N. Platzer, J.M. Bregeault, R. Thouvenot, F. Chauveau, C.
Huet, H. Ledon, Inorg. Chem. 30 (1991) 4409.