Modeling the impact of solid surfaces in thermal degradation processes

(1)

DOI: 10.1002/cphc.201200921

Modeling the Impact of Solid Surfaces in Thermal

Degradation Processes

Christian Tuma,

[a]

Teodoro Laino,*

[a]

Elyette Martin,

[b]

Steffen Stolz,

[c]

and

Alessandro Curioni

[a]

Glycerol (GC, propane-1,2,3-triol), propylene glycol (PG, pro-pane-1,2-diol), and triacetin (TA, propane-1,2,3-triyl triacetate) (Figure 1) are commonly used as solvents, food additives, and humectants in the food and tobacco industries. Their use often involves high-temperature conditions, in which thermal de-composition reactions occur, exposing the human body to po-tentially harmful substances formed during pyrolytic processes. For this reason, it is of great importance not only to know their pyrolysis but also to study the underlying reaction mech-anisms. A detailed picture of the energetics and kinetics of the pyrolytic process is the key to making reliable predictions for the products formed under certain external conditions.

Recently we investigated the pyrolysis of GC, PG, and TA in the gas phase. For GC, we provided a novel overall decomposi-tion scheme whose characteristic is the formadecomposi-tion of glycidol as the first intermediate representing the overall rate-limiting step.[1]_{The investigation of PG dehydration led to the}

discov-ery of propylene oxide as a similar intermediate during the

rate-limiting step.[2]_{Regarding TA, we presented a scheme}

de-scribing its complex thermal decomposition pattern for the first time, in which the rate-limiting step is predicted to consist of the initial elimination of acetic acid.[2]_{For both PG and TA,}

we characterized major decomposition products and the corre-sponding energetics, complementing theoretical results with qualitative data from experiments. Unfortunately, for all the aforementioned substances, a comparison with published ex-perimental data, if available, is a complex task. Exex-perimental studies are always performed under specific catalytic condi-tions that are totally different from a clean pyrolytic gas-phase scenario.

Gas-phase simulations of pyrolytic steps offer ideal condi-tions to study a broad range of possible chemical degradation processes in detail. However, it is also clear that carbohydrate and fat pyrolysis chemistry occurs widely wherever biomass is combusted or food is cooked. These rather complex situations can hardly be simplified by investigating such chemical com-pounds under pure gas-phase conditions. In fact, the presence of solid particles and the interactions of organic compounds with a variety of solid surfaces may strongly affect the energet-ics of the degradation processes, leading to pyrolytic patterns that are totally different to what is known for the gas phase. Several experimental studies related to surface catalyzed degradation reactions of GC have been reported,[3–12]

describ-ing optimal conditions to convert GC either into acrolein or acetol. One reason for these extensive investigations is that GC became available widely as a byproduct of the transesterifica-tion of fats into biodiesel. Regarding the reactivity of GC in the presence of catalytically active surfaces, only few computation-al studies exist in the literature : density functioncomputation-al theory simu-lations on the reactivity of GC over zeolites[13] _{and transition}

metals.[14] _{More recently, Vlachos and coworkers}[15] _reported

a QM/MM investigation on the acid catalyzed dehydration of polyols in liquid water. However, to the best of our knowledge, neither experimental data nor computational results for the pyrolysis of PG and TA over solid surfaces/materials are avail-able in the literature.

Nonetheless, the information available from publications in this field is relevant only to conditions extremely different from those found in everyday life. The latter involve the exposure of GC, PG, and TA to carbon and metal-oxide particles, as in ciga-rettes,[16]_{or when pots and pans are used in high-temperature}

food-preparation processes. For this reason, our study focuses on the rate-limiting steps of the gas-phase decomposition of GC, PG, and TA predicted in our recent works[1, 2]_{and simulates}

how barrier heights change when these reactions occur at the surfaces of different solids, characterizing the key features that

Figure 1. Chemical structures of glycerol (GC), propylene glycol (PG), and tri-acetin (TA).

[a] Dr. C. Tuma,+_{Dr. T. Laino, Dr. A. Curioni}

IBM Research—Zurich, Sumerstrasse 4 CH-8803 Rschlikon (Switzerland) E-mail: teo@zurich.ibm.com [b] Dr. E. Martin

Philip Morris International R&D, Quai Jeanrenaud 5 CH-2000 Neuchtel (Switzerland)

[c] Dr. S. Stolz

University of Twente, Faculty EEMCS

PO Box 217, NL-7500 AE Enschede (The Netherlands) [+

] Current address:

Konrad-Zuse-Zentrum fr Informationstechnik Berlin Takustrasse 7, D-14195 Berlin (Germany)

Supporting information for this article is available on the WWW under http://dx.doi.org/10.1002/cphc.201200921.

2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim ChemPhysChem 2013, 14, 88 – 91

88 CHEM

PHYS

CHEM

COMMUNICATIONS

(2)

are relevant for thermal degradation processes. Our candidates for the solid materials are amorphous carbon, amorphous silica, crystalline zirconia, and crystalline alumina. Computation-al details and methods are available in the Supporting Informa-tion.

The pseudo-amorphous carbon surface contains reactive sites (low-valent carbon atoms). They interact with oxygen atoms of the GC hydroxyl groups. With a barrier of 19 kcal mol 1_{, primary hydroxyl groups are split, forming an}

energeti-cally more stable surface oxide intermediate (chemisorption). This intermediate is connected to a barrier height DE° _of

25 kcal mol 1_{for the formation of glycidol, see Figures 2 and 3.}

The corresponding transition structure formed after C O bond splitting is similar to a primary carbocation, and its relative sta-bilization when compared with the situation for the unimolec-ular reaction step in the gas phase[1]_(DE°₌_{59 kcal mol} 1_{) can}

be attributed to the energy released by incorporating the leav-ing oxygen atom into the amorphous carbon framework. This is the main effect, which even leads to a negative apparent barrier height, that is, the energy for the adsorption of GC from the gas phase at the reactive site of the amorphous carbon surface is larger than the barrier height for the conver-sion of adsorbed GC to glycidol.

GC forms hydrogen bonds with surface oxygen atoms of silica and zirconia (Figure 3). The corresponding PBE + D physi-sorption energies are 27 kcal mol 1_{in the case of amorphous}

silica and 20 kcal mol 1 _{for crystalline zirconia. Here,}

desorp-tion competes with the activadesorp-tion process for the formadesorp-tion of glycidol. The intrinsic total energy PBE + D barrier heights are 59 kcal mol 1_{(silica) and 62 kcal mol} 1_{(zirconia), that is, two to}

three times the energy required for the desorption of GC. Fur-thermore, these intrinsic barrier heights are very close to that for the same reaction step in the gas phase (59 kcal mol 1_).

The clean alumina (0001) surface, however, turns out to be more reactive than the other oxide material surfaces. This is

due to the Lewis acidic character of the surface aluminum atoms. As already shown in previous studies related to alumina (0001) water adsorption,[17, 18]_{a barrier of only a few kcal mol} 1

is associated with the proton transfer from physisorbed water molecules to surface oxygen atoms and the simultaneous crea-tion of surface Al OH groups. GC hydroxyl groups exhibit a slightly higher Brønsted acidity than water and, as a conse-quence, lead to a direct (i.e., barrierless) chemisorption of GC via those hydroxyl groups that do not form an internal GC hy-drogen bond. For GC adsorbed via a primary, its secondary, or both a primary and the secondary hydroxyl group, we obtain adsorption energies of 36, 40, and 71 kcal mol 1_,

respec-tively. In these configurations, the secondary hydroxyl group always forms an internal hydrogen bond with a free primary hydroxyl group, resembling the structural entry point for the formation of glycidol. Corresponding barrier heights are 45 kcal mol 1_{if starting from the adsorption via the secondary}

hydrox-yl group only and 50 kcal mol 1_{when GC is adsorbed via both}

the secondary and a primary hydroxyl group. In the latter case, we also arrive at a negative apparent barrier height, a situation similar to the case of amorphous carbon, see above.

Similar to GC, PG forms hydrogen bonds to surface oxygen atoms of silica and zirconia. The intrinsic barrier heights for the formation of propylene oxide are virtually identical to the cor-responding gas-phase value of 59 kcal mol 1_{, see Figure 4.}

Again, this is two to three times the physisorption energy, ren-dering the desorption of PG more likely than its decomposition at the surface.

A much stronger attraction between PG and the material surface is predicted in the case of alumina ( 44 kcal mol 1_).

Adopting a configuration suitable as precursor for propylene oxide formation, the PG secondary hydroxyl group in the ad-sorption structure both acts as electron donor for an aluminum surface atom and forms an internal hydrogen bond with the primary hydroxyl group. A structure very similar in nature, as-sociated with a barrier height of 50 kcal mol 1 _{for the epoxide}

formation, was already found for GC chemisorbed at the alumi-na surface (see previous subsection). This is consistent with the value of 48 kcal mol 1 _{found for the intrinsic barrier height of}

PG to propylene oxide at the alumina surface.

As for GC, the smallest physisorption energy of PG at the material surface is obtained for amorphous carbon ( 8 kcal mol 1_{). However, overcoming a modest barrier of 15 kcal mol} 1_,

a much more stable surface oxide can be formed (chemisorp-tion), such as in the case of GC. Accordingly, this intermediate is the precursor for the formation of the epoxide structure con-nected to a barrier height of 24 kcal mol 1_{. There is, however,}

a lower lying transition path, with a barrier height of only 19 kcal mol 1_{, leading to acetone which is also energetically}

more stable than propylene oxide. Here, the proton of the sec-ondary hydroxyl group is transferred to the surface framework, and a hydride transfer from the secondary to the primary carbon atom, which previously splits off from the surface oxygen bridge, completes the formation of acetone. Note that this reaction step represents a shortcut compared with the de-composition in the gas phase, where the formation of acetone involves the formation of propylene oxide as its precursor[2] to-Figure 2. DFT (PBE + D) total energy reaction profiles for the decomposition

of glycerol in the gas phase and at the surface of various materials (numbers in kcal mol1

).

89 CHEM

PHYS

CHEM

(3)

gether with an intrinsic activation energy about four times higher than the one observed here starting from the surface intermediate.

TA contains ester groups which are less reactive than the hy-droxyl groups of GC or PG. This reduces its interaction with the amorphous carbon surface to a level as weak as dispersion. Moreover, the carbon surface does not even exhibit partial ionic character, such as that of oxide materials, and therefore, it can be expected that the decomposition of TA is hardly af-fected by the presence of amorphous carbon. However, attrac-tive forces exist between carbon atoms of the ester groups

and the surface oxygen atoms, as well as between oxygen atoms of the ester groups and the Lewis acidic aluminum atoms of the alumina surface. Including van der Waals interac-tions, the TA physisorption energies for silica, zirconia, and alu-mina sum up to 10, 19, and 28 kcal mol 1_{, respectively,}

see Figure 5. The corresponding intrinsic barrier heights for the decomposition of TA into propene-1,3-diol diacetate and acetic acid are slightly below the value of 43 kcal mol 1_{obtained for}

the gas phase.[2]_{This reduction is due to a transition structure}

stabilization effect resulting from an interaction with the polar oxide surface. This effect is less pronounced for alumina than

Figure 3. Decomposition of glycerol (column 1) to glycidol and water (column 3) over amorphous carbon (row 1), pseudo-amorphous silica (row 2), crystalline zirconia (row 3), and crystalline alumina (row 4) via the transition structures shown in column 2 (atom colors : red—O, white—H, blue—C, yellow—Si, green— Zr, purple—Al).

90 CHEM

PHYS

CHEM

(4)

for silica and zirconia because here the structural flexibility of the TA molecule at the surface is smaller owing to stronger binding effects.

In summary, we performed first-principles simulations for the initial steps of the thermal decomposition of GC, PG, and TA over the surfaces of pseudo-amorphous carbon and silica, crystalline zirconia, and crystalline alumina. Our results show that surfaces that do not contain reactive sites hardly influence the decomposition behavior. We have seen, for example, that in the case of silica and zirconia, the intrinsic barrier heights virtually do not change with respect to the values obtained for the gas phase. Furthermore, the adsorption energies are much smaller than the intrinsic activation energies, which renders de-sorption more likely than a decomposition event at the materi-al surface. The picture changes when the adsorption is con-trolled by active surface sites. In the case of the clean alumina surface, we observed a strong interaction between Lewis acid

aluminum atoms and the hydroxyl groups of GC and PG. On the one hand, this leads to a strong binding between the alco-hol molecule and the surface (chemisorption) ; on the other hand, it can increase the acidity of the hydroxyl groups such that the barrier heights of the epoxide formation decrease. Overall, this results in a significant speedup of the decomposi-tion of GC and PG with respect to the gas phase. Of all the oxide material surfaces investigated in this work, the alumina surface exhibits the largest TA adsorption energy. Although the intrinsic barrier for the decomposition into propene-1,3-diol di-acetate and acetic acid remains virtually constant at all oxide surfaces, in the case of alumina the energy required for the de-sorption of TA is high enough for the decomposition step to become competitive. As a consequence, an acceleration effect for the decomposition of TA can be expected at the clean alu-mina surface. The effects induced by the surface of amorphous carbon differ from those of the oxide materials: Low valent carbon atoms induce the chemisorption of GC and PG, forming surface alkoxides. They not only exhibit a significantly lower barrier for the formation of the epoxides, but also open the door to decomposition mechanisms not predicted for the gas phase.

Keywords: computer chemistry · first-principles calculations · reaction profiles · surfaces · thermal degradation

[1] T. Laino, C. Tuma, A. Curioni, E. Jochnowitz, S. Stolz, J. Phys. Chem. A 2011, 115, 3592.

[2] T. Laino, C. Tuma, P. Moor, E. Martin, S. Stolz, A. Curioni, J. Phys. Chem. A 2012, 116, 4602.

[3] H. Atia, U. Armbruster, A. Martin, J. Catal. 2008, 258, 71.

[4] A. Alhanash, E. F. Kozhevnikova, I. V. Kozhevnikov, Appl. Catal. A 2010, 378, 11.

[5] E. Tsukuda, S. Sato, R. Takahashi, T. Sodesawa, Catal. Commun. 2007, 8, 1349.

[6] C.-W. Chiu, M. A. Dasari, G. J. Suppes, W. R. Sutterlin, AIChE J. 2006, 52, 3543.

[7] M. Pagliaro, R. Ciriminna, H. Kimura, M. Rossi, C. Della Pina, Angew. Chem. 2007, 119, 4516; Angew. Chem. Int. Ed. 2007, 46, 4434.

[8] C.-W. Chiu, A. Tekeei, W. R. Sutterlin, J. M. Ronco, G. J. Suppes, AIChE J. 2008, 54, 2456.

[9] L. Cheng, X. P. Ye, Catal. Lett. 2009, 130, 100.

[10] S.-H. Chai, H.-P. Wang, Y. Liang, B.-Q. Xu, Green Chem. 2008, 10, 1087. [11] S.-H. Chai, H.-P. Wang, Y. Liang, B.-Q. Xu, Appl. Catal. A 2009, 353, 213. [12] S. Sato, M. Akiyama, R. Takahashi, T. Hara, K. Inui, M. Yokota, Appl. Catal.

A 2008, 347, 186.

[13] K. Kongpatpanich, T. Nanok, B. Boekfa, M. Probst, J. Limtrakul, Phys. Chem. Chem. Phys. 2011, 13, 6462.

[14] D. Coll, F. Delbecq, Y. Aray, P. Sautet, Phys. Chem. Chem. Phys. 2011, 13, 1448.

[15] S. Caratzoulas, T. Courtney, D. G. Vlachos, J. Phys. Chem. A 2011, 115, 8816.

[16] A. Rossmann, D. Bernhard, G. Wick, IUBMB Life 2005, 57, 805.

[17] K. C. Hass, W. F. Schneider, A. Curioni, W. Andreoni, Science 1998, 282, 265.

[18] K. C. Hass, W. F. Schneider, A. Curioni, W. Andreoni, J. Phys. Chem. B 2000, 104, 5527.

Received: November 8, 2012

Published online on November 23, 2012 Figure 4. DFT (PBE + D) total energy reaction profiles for the decomposition

of propylene glycol in the gas phase and at the surface of various materials (numbers in kcal mol 1

).

Figure 5. DFT (PBE + D) total energy reaction profiles for the decomposition of triacetin in the gas phase and at the surface of various materials (numbers in kcal mol1_).