Formation of environmentally-persistent free radicals (EPFR) on α-Al2O3 clusters

Alumina oxides assume prominent catalytic applications in a wide range of industrial processes. However, alumina surfaces also serve as potent promoters in the heterogeneous formation of the notorious environmentally-persistent free radicals (EPFR). Herein, we theoretically examine dissociative adsorption mechanisms of phenol molecules over Al2O3 and hydrated Al2O3$nH2O clusters that mimic dehydrated and hydrated alumina structures, respectively. We show that fission of the phenol's hydroxyl bond over dehydrated alumina systematically incurs lower energy barriers in reference to the hydrated structures. A 1,2-water elimination step marks the most feasible channel in the interaction of phenol with hydrated clusters. The relevance of the acidity sites to the catalytic activity of alumina is clearly supported by the finding that the catalytic activity of the alumina surface in producing the phenoxy/phenolate species reversibly correlates with the degree of hydroxyl coverage. Desorption of adsorbed phenolates requires sizable desorption energies, and thus is expected to facilitate surface-mediated condensation into dioxin-like moieties.


Introduction
Alumina oxide represents one of the most important catalytic and catalyst-support materials, 1-3 for chemical reactions. It also nds direct applications in dielectric materials, substrates for electronics, packing materials and radiation dosimeters. 4,5 Over the last two decades, the surface chemistry of alumina oxides has been a distinct research topic in catalysis chemistry. Most of these studies have focused on the surface/water interface conditions, where the chemistry of the surface is greatly affected and accordingly alters its reactivity and catalytic performance. [6][7][8] The interaction of alumina surfaces with water molecules 6,[9][10][11][12][13][14] results in water-and hydroxyl-covered surfaces, in which the degree of coverage is highly sensitive to temperature. 15 Heating and cooling processes can reversibly either add or remove hydroxyl groups from the surfaces, as conrmed experimentally by IR and NMR measurements. 16,17 For instance, X-ray diffraction data by Dyer et al. 18 revealed the formation of Al(OH) 3 on the g-Al 2 O 3 surface, which disappeared aer heating to 473.15 K. This has also been conrmed by Raman and infrared spectroscopy where several sharp peaks in the region around 3600 cm À1 have been observed aer hydration of the g-Al 2 O 3 powder. X-ray photoemission experiments by Liu et al. 19 pointed out to water dissociation; most likely at surface defect sites.
It has become evident that, 3,15,20,21 the differences of the coordination of Al atoms of the alumina surface lead to form different terminal hydroxyl termination or interactive phases. An experimental study by Pimentel and McClellan 22 revealed that, surfaces encompasses two main categories of surface hydroxyl groups; isolated hydroxyl groups with a sharp infrared bands (>3600 cm À1 ), and self-associated hydroxyl groups with a broad, intense band (<3600 cm À1 ). Apart from this, selfassociated hydroxyl groups differ from the isolated groups in the presence of the hydrogen bond connection between adjacent moieties. Detailed examination of chemical sites on Al 2 O 3 surface 23,24 conrmed the ndings of Pimentel and McClellan. 22 Overall, ve different classes of isolated hydroxyl groups have been detected by IR spectroscopy. 21,24 Knözinger et al. 3 attribute the difference in vibrational frequencies among these classes to the coordination number of alumina (i.e.; either tetrahedral and/or octahedral coordination).
The chemical makeup (i.e.; adsorption and decomposition) of the hydroxyl groups over alumina surface constitutes a key factor in clarifying the reactive/catalytic nature of alumina. As stated earlier, heating can reversibly remove the hydroxyl groups from the alumina surface, 16,17,25 leading to the formation of Lewis acid-base pair sites, and thus increases the catalytic activity of the surface. Along a similar line of enquiry, Ballinger and Yates 26 have performed a detailed investigation into the correlation between Al-OH group removal and the number of Lewis acid sites via CO molecular adsorption. The authors have found that, heating leads to the decomposition of the associated hydroxyl groups, leaving only isolated hydroxyl groups on the surface at a temperature range of 475-800 K. This has been viewed as the driving force of the catalytic activity of alumina, either via hydrogen bonding to the hydroxyl group and/or via lone pair donation to Lewis acid sites. A study by Hendriksen et al. 27 measure the heat of adsorption of water on both aand galumina powders via immersion calorimetry. It has been indicated that, the heat of adsorption on both termination depends primarily on the degree of hydration prior to water adsorption.
Over a wide range of operational pressure and temperatures, alumina exhibits a variety of crystallographic structures such as a phase or corundum, which is the most thermodynamically stable phase, 28 and other alumina phases such as d, b, k, h, q and g. 29,30 Literature density function theory (DFT) investigations have considered water adsorption and dissociation over a-Al 2 O 3 surfaces from the dry to the fully hydroxylated surfaces. A detailed analysis and discussion on the subject was presented by Wang et al. 31 The authors have found that, a fully hydrated alumina a-Al 2 O 3 (0001) surface, when the surface Al atom is triply hydroxylated, exhibits relatively inertness toward further hydroxylation. In an insightful analysis into the interaction of water molecules and a a-Al 2 O 3 (0001) surface, Wirth and Saalfrank 32 found that 1,2-dissociation pathway is the most kinetically feasible mechanism, where Al surface atoms are hydroxylated and the nearby oxygen atoms are protonated.
Recent theoretical investigations addressed the interaction of hydrated and dehydrated a-Al 2 O 3 surfaces towards different molecules. [33][34][35] Shukla and Hill 33 have performed a detailed investigation of the molecular interaction of 2,4,6-trinitrotoluene (TNT) molecule with the a-Al 2 O 3 (0001) surface demonstrating that TNT acquires strong covalent interaction in a parallel orientation with the surface. The same group 34 found that molecular adsorption of TNT over hydroxylated alumina surface is much weaker when compared to the dehydrated surface. Similarly, Rohmann et al. 35 investigated molecular interaction of carbon monoxide with hydroxylated alumina surface, reporting that molecular CO exhibits strong interaction with a 25% degree of hydroxylated surface and this interaction is signicantly stronger when compared with that over the clean dehydrated alumina. In a follow-up study on H 2 S removal via g-Al 2 O 3 surfaces with different hydroxyl coverage, Ren et al. 36 revealed that, the bonding strength of H 2 S on hydrated alumina surface is lower than that the over a dehydrated termination. For example, calculated binding energy for 8.9 OH per nm 2 hydroxyl coverage for g-Al 2 O 3 (110) amounts to À68 kJ mol À1 whereas the corresponding value of the dehydrated surface stands at À114 kJ mol À1 . XPS studies by Kelber et al. 37 investigated the interaction between Cu and the hydrated a-Al 2 O 3 (0001) surface to report an enhanced interaction of copper to the surface of alumina owing to the presence of a surface hydrated layer. There is a rather limited literature account on reactions of hydrocarbons on alumina clusters. Most of relevant studies have considered extended surfaces of alumina. 34,35,38 Nonetheless, it has been shown that alumina clusters encompass effective catalytic active sites present on surfaces; i.e., the surface Al-O bonds. 39,40 Therefore, alumina clusters in many theoretical studies (i.e. (Al 2 O 3 ) n¼4 ) are being applied as representative models for both a-Al 2 O 3 and g-Al 2 O 3 surfaces. Alumina clusters imitate both the experimental measurements and DFT surface properties. Additionally, clusters 41,42 were shown to illustrate an effective stability against relaxation oen encountered in extended alumina surfaces.
Phenol molecule is among the most commonly discussed precursors for the synthesis of poly chlorinated dibenzo-pdioxins and dibenzofurans (PCDD/Fs). 43 50 In our recent study, 51 we theoretically investigate the interaction of phenol molecule with the dehydrated a-Al 2 O 3 (0001) surface. We have found that, alumina surface is catalytically very active in producing phenolate (i.e., surface-bounded phenoxy) with a modest activation energy of 43 kJ mol À1 , thus, signifying as an accessible route for the formation of persistent surface-bound radicals. However, when considering the temperature range of the cooling-down zone of the combustor (i.e.; 746-1146 K), all models of alumina (i.e.; polymorphs, surfaces and\or clusters) assume different hydration coverage. Therefore, the effect of the degree of hydration, on dissociative adsorption of phenol, can play an important role in attaining a detailed understanding of the catalytic activity of the surface. Thus, it is very important to shed light into the generation of the environmentally persistent free radicals (EPFRs) 52,53 on alumina surface with varying hydroxyl coverages.
In summary, the current work is part of our ongoing effort to reveal mechanistic and kinetic factors dictating the catalytic role of alumina surfaces in surface-mediated formation of phenoxy-type EPFR. The study provides a detailed investigation of the interaction of a phenol molecule with different hydration coverage of alumina clusters and therefore offers important insights into the role of the surface adsorbed water in the catalytic activity of alumina toward the generation of adsorbed phenoxy radicals. The present study expands on our recent investigation, 51 in which we considered a dehydrated a-Al 2 O 3 (0001) surface. Overall, this study has a three-fold aim; (i) to report modes of interaction of phenol with different hydration coverages of alumina clusters, (ii) to compare between the catalytic activity of hydrated and dehydrated alumina clusters toward the formation of the adsorbed phenolate via two distinct pathways; namely direct ssion of the phenol' O-H bond and elimination of water molecules and to (iii) to construct a simplied kinetic model for the conversion of phenol into adsorbed phenoxy radicals.

Methodology
DMol 3 program 54,55 performs all structural optimization, energy calculations and vibrational frequencies based on the generalized gradient approximation (GGA) of the Perdew-Burke-Ernzerhof (PBE) functional. 56 The computational methodology comprises Grimme dispersion correction, 57 together with the basis set of double numerical plus polarization (DNP). Total energy and forces on each atom converge with tolerance to less than 10 À6 eV and 10 À4 eVÅ À1 , respectively.
We computed the binding energies (E b ) for the adsorption based on the energies of a phenol/Al 2 O 3 system E cluster+phenol , the clean alumina cluster E cluster , and an isolated phenol molecule E phenol : Reaction rate constants were estimated based on the conventional transition state theory (TST), 58 we t reaction rate constants to the two-parameters Arrhenius formula, k(T) ¼ A exp(ÀE a /(RT), over the temperature range of 300 K to 1000 K, according to: where: k B is Boltzmann's constant, h signies Plank's constant, R is the universal gas constant. Calculated vibrational frequencies enable to calculate thermodynamic functions; i.e., entropies and enthalpies. DS s and DH s stand for temperature-depended entropy and enthalpy of activation, correspondingly. s e represents the reaction degeneracy number. DMol 3 computes entropy values based on well-documented statistical thermodynamics equations. 59 The enthalpy (at the desired temperature) has been competed relatively to the zero-point energy at 0 K.

Water adsorption on Al 2 O 3 cluster
Starting with dehydrated alumina clusters, Fig. 1 displays two water-free geometrical monomers of alumina near-cyclic AlO 3 Al that adapt D 3h symmetry, and a linear O]Al-O-Al]O chain structure. As illustrated in Fig. 1, the cyclic dehydrated alumina, AlO 3 Al, contains three Al-O bonds with distances of $1.82Å. This structure was found to be energetically less stable 60 than the linear chain structure due to the repulsion between the two cations (Al 3+ ). In the latter conguration, alumina adapts a linear structure with one bridging and two non-bridging oxygen atoms. Calculated bond angle and bond lengths, in both structures, are in agreement with the literature analogous values. 60,61 Moving on now to consider how the geometries of alumina cluster are compared with those in alumina surface. Fig. 2 displays top and side views of a-Al 2 O 3 (0001) surface (i.e., the thermodynamically most stable alumina surface 62,63 ) and the D 3h structure. As can be seen from the gure, Al-O bonds lengths in both structures are within 0.01Å. This close agreement serves to testify that considered cyclic dehydrated structure contains similar active sites to those present in the surface.
A study by Johnson and Panas 60 examine water adsorption and hydrolysis on molecular Al oxides and hydroxides solvation versus cluster formation. They found that, dehydrated alumina clusters, both cyclic and linear chain, reacts with water molecules through accessible energy barriers in the range of 33-110 kJ mol À1 . Further addition of water molecules leads to a series of Al 2 O 3 $nH 2 O structures. Fig. 3 illustrates mechanism of successive addition of water molecules to the dehydrated cluster of alumina (i.e.; cyclic and linear clusters). Starting with the cyclic structure of alumina, D 3h , the addition of a water molecule to the Al 2 -O 1 bond results in the formation of the Al 2 O 3 $1H 2 O (c) cluster (c) which has been reported to be the most stable structure among the various congurations of Al 2 O 4 H 2 clusters. 60 Typically, this structure assumes cis and trans congurations, with very similar thermodynamic stability. 60 Thus, we have performed our subsequent calculations considering the cis planar molecule congurations with the C 2v symmetry as it is marginally more thermodynamically stable than its trans counterpart; by 3 kJ mol À1 . Calculated Al 2 -O 2 and Al 2 -O w1 bond lengths amount to 1.75Å and 1.70Å, respectively. These values concur very well with the corresponding literature values by Johnson and Panas; 60 i.e., 1.75Å and 1.69Å, respectively.   Fig. 3 has been reported by Johnson and Panas. 60 The authors found that formation of (c) and (e) structures display an exothermic process with calculated enthalpies of À310 kJ mol À1 and À184 kJ mol À1 , respectively. The formation of structure (g) is slightly endothermic with a low activation enthalpy at 33 kJ mol À1 . 60

Reaction of phenol with dehydrated alumina cluster
We studied the interaction of a phenol molecule with the two dehydrated alumina clusters (i.e.; cyclic and linear clusters in Fig. 3a and b). Two different reactions have been considered, denoted as pathway a and pathway b. These pathways involve the interaction of phenol molecule with cyclic and linear dehydrated alumina clusters, respectively. The optimized geometries for the reactants, transition states and products are presented in Fig. 4. Table 1 lists energies for physisorbed congurations and dissociative structures as well as prominent interatomic distances for all pathways.
We rst consider the interaction along the cyclic alumina structure, pathway a. Fig. 4(a) shows, phenol molecule physically adsorbed on the cluster through a van der Waals interaction between O phenol atom and Al cluster atom with a binding distance of 1.909Å. The interaction proceeds via ssion of the hydroxyl's bond over Al-O cluster bond where phenoxy moiety attaches to Al cluster and deducted hydrogen atom binds to O cluster atom. This has resulted in the optimised dissociation product P a . Hence, the length of this bond increases by 13.4% with respect to that of the clean cyclic cluster (1.822Å). This value concurs with results obtained in our recent study 51 of the dissociative adsorption of phenol on the dehydrated a-Al 2 O 3 (0001) surface (11.9%). In comparison, Wang et al. 31 showed that the elongation of Al-O bond increased by 6.4% upon the dissociation of a water molecule over this Al-O bond.
The reaction is highly exothermic with an energy barrier of only 41 kJ mol À1 . Calculated binding energy of the reactant, R a , and the product, P a , are À129 kJ mol À1 and À241 kJ mol À1 , respectively (based on eqn (1) in reference to separated reactants). These values also agree with the results obtained in our recent study. 51 We found that, phenol molecule strongly interacts with the dehydrated a-Al 2 O 3 (0001) surface via physisorbed binding energies in the range of À127 kJ mol À1 to À202 kJ mol À1 , leading to the formation of an adsorbed phenolate with a modest activation energy of 48 kJ mol À1 .  (Fig. 4) increases from the equilibrium 0.97Å distance in the gas phase to 1.016Å. The exothermic formation of the product P b requires a trivial activation energy of 10 kJ mol À1 (via direct ssion of the phenol's O-H bond catalysed by the cluster's Al-O linkage). Calculated binding energies of the reactant, R b , and the product, P b , are À152 kJ mol À1 and À345 kJ mol À1 , respectively (in reference to the separated reactants, Fig. 4). P b appears relatively more stable than P a . Activation barrier via pathway b is lower by 31 kJ mol À1 than that along pathway a. This indicates that, the linear chain cluster of alumina is catalytically more active in producing phenolate than the cyclic dehydrated cluster.
Test computations employing the dependence of the reactivity of alumina cluster on the cluster size has been considered. We investigated physisorption and chemisorption interactions of phenol with the larger (Al 2 O 3 ) 4 reported in the literature. 64  Stable optimised structures are shown in Fig. 5 calculated binding energies of the physisorbed and the chemisorbed structures are À150 kJ mol À1 and À315 kJ mol À1 , respectively, both values are in a good agreement with the corresponding values we obtained for the smaller dehydrated clusters (i.e., À129 to À152 kJ mol À1 and À241 to À345 kJ mol À1 , respectively). Furthermore, binding energies obtained for the larger (Al 2 O 3 ) 4 cluster are consistent with the calculated binding energies in our recent study 51 for the interaction of phenol with a seven-layer slab of the 2 Â 2 (0001) surface alumina (À127 kJ mol À1 to À202 kJ mol À1 ). Thus, it can be concluded that energies of the surface interaction of phenol with alumina surface do not display sensitivity toward the cluster size. The same observation has been observed by Pan et al. 65 in their study of the reaction of 2-chlorophenol over the dehydrated and hydroxylated silica clusters with (SiO 2 ) 3 and (SiO 2 ) 8 clusters.
They have found that binding energies of reactants and products using both structures are within (0.4-17.4) kJ mol À1 .

Reaction of phenol with Al 2 O 3 $nH 2 O cluster
We now turn our attention to investigate the interaction of a phenol molecule with hydrated alumina clusters entailing different hydroxylation coverages, (i.e., all the hydrated alumina clusters shown in Fig. 3). The detailed mechanism of the interaction is presented in Scheme 1. Fig. 6 depicts optimised geometries of reactants, transition states, and products; along with their energy proles. First, we have studied the interaction of a phenol molecule with the (c) structure, Fig. 3. Three possible pathways have been considered, denoted as c 1 , c 2 and c 3 in Scheme 1. As Scheme 1 demonstrates, these reaction pathways characterise (i) ssion of hydroxyl's bond over the same Al 2 -O 2 bond (bridging bond), (ii) H 2 O elimination mechanism over the same Al-OH bond (1,2 non-bridging) and (iii) H 2 O elimination route over two Al-OH bonds (1,4 non-bridging linkage), respectively (Fig. 3 illustrates atomic numberings). The reaction is initiated by the physisorbed-type interaction between O phenol and Al cluster via binding spacings of 2.007Å, 1.979Å and 1.997Å in the reactants R c 1 , R c 2 and R c 3 , respectively (Fig. 6). Calculated binding energies of the three reactants stand at À143 kJ mol À1 , À147 kJ mol À1 and À129 kJ mol À1 , respectively (in reference to the separated reactants). Regarding to the structural changes induced to the cluster; we observe that Al 2 -O 2 bridging bond in Fig. 6, R c 1 elongates by 2.4% when compared with the corresponding value in the Fig. 3(c) structure while Al 1 -O 1 in R c 2 and Fig. 6, R c 3 increases by 8.1% and 8.2%, respectively. Dissociative desorption of phenol via these pathways leads to the formation of P c 1 , P c 2 and P c 3 (Fig. 6) intermediates with reaction energies of À225 kJ, À167 kJ mol À1 and 17 kJ mol À1 , respectively. The thermodynamic penalty associated with the product P c 3 most likely stem from the low coordinated Al 1 that display a 2-fold coordination. 66 Next, we examine the dissociative adsorption of phenol over the Fig. 3(g) structure. The reaction proceeds via the active Al] O double bond site of the cluster. We found a slight increase in the O-H bond length of the physisorbed phenol in pathway g (1.009Å, Scheme 1), in reference to the gas phase value (0.97Å). The corridor g in Scheme 1 characterises ssion of the O-H bond in the phenol molecule in a noticeably exothermic reaction of 360 kJ mol À1 . The physisorbed state Fig. 6 (R g ) in this channel resides 166 kJ mol À1 below the separated reactants. Compared to Fig. 3(d), the size of hydrated alumina cluster (g) with one water molecule is smaller and it has a planner structure. Therefore, as temperature increases, molecules move faster and they are more likely to collide and hence positively affect the collision frequency, A.
Reaction pathways c 1 , c 2 and g proceeds via modest energy barriers of 71 kJ mol À1 , 104 kJ mol À1 and 23 kJ mol À1 , respectively, whereas pathway c 3 incurs a very sizable energy barrier of 289 kJ mol À1 . Three concluding remarks can be drawn from values in Fig. 6. First, within Al 2 O 3 $1H 2 O system (i.e.; pathways: g, c 1 , c 2 , c 3 ), as expected, reaction pathway g is the most favourable channel, presumably due to the presence of the active Al]O double bond site. Second, phenol dissociation catalysed/accompanied with water elimination (pathways c 2 and c 3 ) is more energy demanding than phenol decomposition through surface-assisted ssion of its hydroxyl O-H bond (pathways c 1 ). For instance, the activation energy in pathway c 1 is lower by 33 kJ mol À1 than that along pathway c 2 . Third, the activation energy in reaction pathway c 3 as well as the calculated binding energy of the product P c 3 are higher when contrasted with the corresponding values in pathway c 2 , this indicates that water elimination preferentially occurs at the same A-OH linkage; i.e., 1.2-elimination step rather than at different A-OH sites (i.e., 1,4-elimination). This has prompted us to exclude this mechanism (i.e.; pathway c 3 ) from further consideration in our subsequent discussion.
Phenol interaction with higher hydrated alumina clusters (Al 2 O 3 $2H 2 O (d) and Al 2 O 3 $3H 2 O Fig. 3(e)) is presented in pathways d and e in Scheme 1, with Fig. 7 illustrates optimised geometries and relative energies for species involved in pathways d and e.
Water elimination via pathways d and e (Scheme 1) results in the formation of a phenolate moiety via activation barriers of 110 kJ mol À1 and 142 kJ mol À1 , respectively. Pathways d and e characterise dissociative adsorption of phenol via water elimination through 1,2-route; i.e., attachment of the phenoxy'O at the Al atom and the simultaneous departure of the OH group attached to the same Al atom with the hydroxyl'H as a water molecule. Higher energy barrier for pathway e in reference to pathways c 2 and d can be rationalised in energy penalty required for the rupture of the Al-O linkage in the former. The degree of hydration of alumina; could be viewed as a factor that partially dictates the catalytic activity of alumina via controlling the surface-acid properties. 67 In this regard, Al atom in structure (e) (the Al 2 O 3 $3H 2 O system), is coordinated by four hydroxyl groups, whereas Al atom in c 2 and d, is less coordinated presenting more acidic structures as illustrated in Scheme 2. From the scheme, we can see that, the addition of water molecule in the structure dblocks the Lewis active site and accordingly decreases the catalytic activity of the cluster. This view is supported by Digne et al. 68 who argued that the lower the Al atom coordination, the strongerthe Lewis acidity. Calculated atomic Mullikan charges on Al atom for (c), (d) and (e) were shown to have more positive values on the four coordinated Al atom of 1.69e (structure (e)) when compared to the three coordinated Al atom (i.e., 1.52e in structures (c) and (d)). A positive correlation is found between the calculated energy barrier and calculated Mulliken atomic charges; i.e., the higher the atomic charges on Al atom, the higher energy required to break Al-OH bond toward generating water molecule (accompanied with the dissociative adsorption of phenol). Thus, the observed increase in barrier for pathway e in reference to pathways c 2 and d correlates with partial atomic charges; i.e., a prominent electronic descriptor that determines the acidity character.
Overall, we can compare between hydrated and dehydrated alumina from two different perspectives. First, within ssion of hydroxyl's bond reactions over the same Al-O bond; contrasting pathway Fig. 6 c 1 with pathway a reveals that, the energy barrier required in hydrated alumina clusters (i.e.; c 1 ¼ 71 kJ mol À1 ) are signicantly higher than that over the dehydrated one (i.e.; a ¼ 41 kJ mol À1 ), respectively. Second, saturation of Al]O double bonds in the course of the hydroxylation of cluster b into d (refer to Fig. 3) signicantly increases the barrier required for water elimination from only 10 kJ mol À1 to 110 kJ mol À1 . Together, these ndings indicate that dehydrated alumina cluster are more active in producing phenolate than hydrated alumina clusters. This observation agrees with the ndings of Pan et al. 65 who found that, the dehydrated silica cluster (with the presence of two ends of S]O double bonds sites) produces chlorophenolate from 2-chlorophenol through an energy barrier of only 21 kJ mol À1 , compared with 70 kJ mol À1 over the hydrated cluster (the same cluster aer the two Si]O double bonds sites become saturated).
The most obvious nding to emerge from inspection of data in Table 1 is that, in both hydrated and dehydrated alumina, the dissociated product via Al]O double bonds site of the cluster is more stable than that leading to Al-O bond cleaving. For example, reaction energies that mark the formation of P g and P b (double bonds sites) amount to À360 kJ mol À1 and À345 kJ mol À1 , respectively, while the analogous energies of P c 1 and P a are À225 kJ mol À1 and À241 kJ mol À1 , (A-O bond ssion), respectively. This nding is in agreement with those obtained by Pan et al. 65 The authors found that dissociative adsorption of 2-chlorophenol through rupture of its hydroxyl's bond over a S]O double bonds site forms more stable product by 136 kJ mol À1 than that leading to Si-O bond cleaving.
Desorption of the phenolate moiety from both dehydrated alumina cluster and hydrated alumina clusters, was found to be highly endothermic by 334-373 kJ mol À1 . This nding is in accord with those reported in our recent study for the interaction of phenol 51 and chlorophenol 69 molecules over alumina surfaces and silica clusters, 394 kJ mol À1 and 379 kJ mol À1 respectively. This profound stability of phenolate enables it to undergo bimolecular reactions via the so called Eley-Rideal (ER) and Langmuir-Hinshelwood (LH) mechanisms 48 to produce dibenzo-p-dioxin and dibenzofuran; respectively. In alternative corridor, decomposition of phenolate on the surface is expected to initiate the formation of soot.

Kinetic consideration
We conclude the study by implementing a microkinetic analysis of the abovementioned reactions in this study. Fig. 8 represents Arrhenius plots while Table 2 displays the calculated reaction rate parameters. Within the considered temperature limit, reactions pathway c 3 displays the highest activation energy and relatively high dependence on temperature. For all reactions, the tted Arrhenius energy of activation (E a ) largely reects corresponding energy barriers presented in Fig. 3, 6 and 7. Calculations of the conversion-temperature proles were carried out based on a simplied plug ow reactor (PFR) model. The material balance equations were solved by Polymath soware: 70 In the considered reaction, (A / B), A and B represents the molecularly adsorbed phenol and the dissociative structures; respectively. r denotes the reaction rate, F symbolizes the molar ow rate of phenol in mol s À1 , and W stands for alumina catalyst weight in kg. The reaction is considered to be rst order with respect to phenol. The model utilises reaction rate constants given in Table 2. The catalytic destruction of phenol molecule over Al 2 O 3 $nH 2 O catalysts as a function of temperature is presented in Fig. 9. It can be seen that all hydrated coverage of alumina clusters exhibits a highly catalytic destruction activity for phenol molecule, evident by high conversion at low temperatures. Calculated temperature of 90% destruction of phenol molecule (T 90 %) occurs at 350 K, 925 K, 425 K and 425 K for pathways c 2 , c 3 , d and e, respectively. To the best of our knowledge, literature provides no conversion values pertinent to interaction of phenol with a-Al 2 O 3 . Considering the very high desorption energy for phenolate, formation of gas  phase phenoxy should be hindered. The LR and LH mechanisms operated by the highly stable surface-bounded phenolate species, remains largely speculative.

Conclusion
This study sets out to assess the effect of surface acidity on the interaction of alumina with phenol molecules. Interaction of phenol molecule with dehydrated alumina clusters proceeds via ssion of the hydroxyl's bond over either Al-O single bond or Al]O double bond. Clusters with the active Al]O double bond (i.e., b) are catalytically more active in producing phenolate when compared to structures where all bonds are saturated (i.e., a, Al-O single bonds). Overall, c channels that proceeds through ssion of the hydroxyl's bond over Al-O (i.e.; c 1 ) are more feasible than those proceeding via H 2 O elimination routes (i.e.; c 2 and c 3 ). The most obvious nding to emerge from this study is that adding water molecules increases the coordination number of Al atoms and acts to black available Lewis acid sites for uptake of phenol molecules. Simplied kinetics model suggests dissociation of phenol molecules at low temperatures; however, subsequent desorption of adsorbed phenolate into gas phase phenoxy radicals is hindered by very sizable desorption energies.

Conflicts of interest
There are no conicts of interests to declare.