Next Article in Journal
Species-Dependent Splice Recognition of a Cryptic Exon Resulting from a Recurrent Intronic CEP290 Mutation that Causes Congenital Blindness
Previous Article in Journal
Development of Small RNA Delivery Systems for Lung Cancer Therapy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mechanistic Studies on the Dibenzofuran Formation from Phenanthrene, Fluorene and 9–Fluorenone

Environment Research Institute, Shandong University, Jinan 250100, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2015, 16(3), 5271-5284; https://doi.org/10.3390/ijms16035271
Submission received: 2 February 2015 / Revised: 2 March 2015 / Accepted: 4 March 2015 / Published: 6 March 2015
(This article belongs to the Section Physical Chemistry, Theoretical and Computational Chemistry)

Abstract

:
We carried out molecular orbital theory calculations for the homogeneous gas‑phase formation of dibenzofuran from phenanthrene, fluorene, 9-methylfluorene and 9-fluorenone. Dibenzofuran will be formed if ∙OH adds to C8a, and the order of reactivity follows as 9-fluorenone > 9-methylfluorene > fluorene > phenanthrene. The oxidations initiated by ClO∙ are more favorable processes, considering that the standard reaction Gibbs energies are at least 21.63 kcal/mol lower than those of the equivalent reactions initiated by ∙OH. The adding of ∙OH and then O2 to phenanthrene is a more favorable route than adding ∙OH to C8a of phenanthrene, when considering the greater reaction extent. The reaction channel from fluorene and O2 to 9-fluorenone and H2O seems very important, not only because it contains only three elementary reactions, but because the standard reaction Gibbs energies are lower than −80.07 kcal/mol.

Graphical Abstract

1. Introduction

Incineration is an efficient means of waste disposal, owing to significant volume reduction and energy recovery [1]. Unfortunately, a large quantity of toxic byproducts, such as polychlorinated dibenzo‑p‑dioxins (PCDDs) and polychlorinated dibenzofurans (PCDFs), are emitted from the incinerators [2]. PCDD/Fs are among the most toxic environmental contaminants due to their potency as aryl hydrocarbon receptor (AHR) ligands [3,4]. PCDD/F emissions have caused serious public concern and sparked the controversy about the acceptance of incineration as a viable solid waste treatment method. Suppressing the PCDD/F formations may be more efficient than end-of-pipe technology for PCDD/Fs in incineration in the future. Despite many studies, there remains much uncertainty about the predominant formation pathways of PCDD/Fs in combustion systems.
PCDD/Fs are formed mainly through three pathways in combustion [5,6,7,8,9]: (I) gas-phase formations from molecular precursors at high temperatures (>600 °C); (II) catalytic reactions on the surface of fly ash across temperature range of 200–600 °C; and (III) chlorination and dechlorination of the existing PCDD/Fs. Chlorophenols are generally considered as the predominant precursors of PCDD/Fs and implicated as key intermediates in the de novo synthesis [10,11]. However, diversify sources of PCDD/Fs have been demonstrated, such as catechol [12], semiquinone radicals [13] and captan pesticide [14] and polycyclic aromatic hydrocarbons (PAHs) [15]. Both catalytic and non‑catalytic condensations of chlorophenols showed that the ratio of PCDDs and PCDFs either in gas phase or on fly ash is more than one, but PCDFs are major congeners in waste incinerators [11,16,17,18]. The catalytic reaction of amorphous 12C- and 13C-labeled carbon revealed that approximately half of PCDDs are formed via the condensation of monocyclic aromatic precursors; however, PCDFs are primarily released from the oxidation of the precursors with two or more aromatic rings, such as PAHs and elemental carbon [19]. As evidence, the multiple regression analysis showed that the coefficient of determination between the toxic equivalent (TEQ) of PCDD/Fs and the concentration of naphthalene, fluorene and phenanthrene was 0.85 [20]. PAHs are recognized as potential carbon sources of PCDFs, because the causes are not only the structural similarities, but also high concentrations in fly ash and tail gas. Total amounts of PAHs in the different ash samples vary from 14 to 77,000 μg/kg [21], and the concentrations in gas phase range from 566 to 8374 μg/Nm3 [22,23].
The chlorination of dibenzofuran (DF) could play an important role in forming PCDFs. On the one hand, the complete distributions for PCDF congeners shows broadly consistent between the chlorination model and the data from incinerators [6,24]; In addition, the chlorination of DF vapor by copper (II) chloride is proved to be a highly efficient process [24]. On the other, the high concentration of DF supports the hypothesis. The concentration of DF in flue gas is 50–700 times that of total PCDFs, and 2600 times during combustion disturbance [16,25].
Cosentino et al. [15] reported that oxidative pathways of pyrene and benzodibenzofuran—a model for the de novo synthesis—are possible routes of dibenzofuran formation. Meanwhile, theoretical studies on the oxidation mechanisms in combustion processes of aromatic and unsaturated hydrocarbons [26,27] suggested that the oxidation of PAHs is a source of DF and dibenzo‑p‑dioxin (DD). The nearly identical congener’s distributions of PCDD/Fs among different emission sources serves as an evidence for the occurrence of heterogeneous mechanism similar to the formation mechanism from PAHs. For the further investigation of the DD/DF formation mechanism from PAHs and seeking after the heterogeneous mechanism of PCDD/Fs in incineration processes, we carried out a density functional theory study on the mechanistic aspects of DF formation from phenanthrene (Phe), fluorene (Flu), 9-menthylfluorene (9MF), as well as 9-fluorenone (9FO).

2. Results and Discussion

The structure formulas of Phe, Flu, 9MF, 9FO and DF are shown in Scheme 1. Those structures show high similarity to each other. Therefore, if a generic modeling is used, Phe will be got while “–X–” stands for “–CH=CH–”, etc. The change routes mentioned in this study are mainly focused on breaking the middle ring of these molecules, and then replacing the “–X–” by oxygen atom.
Scheme 1. Illustrative structural diagrams for phenanthrene, fluorene, 9-methylfluorene, 9‑fluorenone, dibenzofuran and a generic modeling for them.
Scheme 1. Illustrative structural diagrams for phenanthrene, fluorene, 9-methylfluorene, 9‑fluorenone, dibenzofuran and a generic modeling for them.
Ijms 16 05271 g004

2.1. Reaction from Phenanthrene to Dibenzofuran

In a flame, aromatic compounds might be attacked by active radicals, and then be disintegrated [28]. Hydroxyl radical (∙OH) turns out to be one of the dominant chain carriers, and plays an important role in various aspects of combustion chemistry [29,30]. In addition to ∙OH, chlorine monoxide radical (ClO∙) becomes important in the flame of chlorinated organic compounds. The mole fraction of ClO∙ reached 4.01 × 10−4, i.e., 1.3 times the concentration of ∙OH, when ammonium perchlorate was used as an oxidizer in a hydrogen-breathing combustion engine [31]. The combustion of CH3Cl, CH2Cl2, and CHCl3 in air showed that the mole fractions of ClO∙ radicals under fuel lean conditions gained at least 100 times higher than those of ∙OH radicals at temperatures lower than 900 K [32].
There are 14 carbon atoms with sp2 hybridization in the three fused aromatic rings of Phe, and ∙OH radical can add to any one of the carbon atoms. Among all the possible additions, C8a (or C10a, they are the equivalent atoms considering symmetry, as shown in Scheme 1) and C9 (or C10) are chosen as the only reaction sites. Radicals could add at other positions in the Phe molecule, but those are not relevant to the formation of DD and DF and hence they were excluded from mechanistic analysis. The bonding between ∙OH and C8a will finally create DF, and Figure 1A illustrates the predicted reaction route. Activation enthalpies ( H ), activation Gibbs energies ( G ), standard enthalpy changes ( r H m o ) and Gibbs energy changes ( r G m o ) for every elementary reaction are listed in Tables S1 and S2 (attached as supplementary materials). We calculated the data every 100 K from 1500 to 298.15 K. The bonding reaction between ∙OH and C8a of Phe is an exothermic process, and the reaction enthalpy amounts to –4.17 kcal/mol at 1000 K. The activation enthalpy of this addition reaches as low as 6.20 kcal/mol. The ring-opening reaction of the addition product—Phe-B1—is an endothermic reaction with a potential barrier of 31.73 kcal/mol. The ring-opening seems to be an unfavorable process for the standard Gibbs energy change is 40.41 kcal/mol at 1000 K, although the value is already significantly lower than those of the C8a–C9 directly dissociation (177.90 kcal/mol). The subsequent H‑migration and C2H3-leaving reactions eventually produce DF. As a whole, the transformation (Phe + ∙OH → DF + ∙C2H3) absorbs 22.28 kcal/mol heat at 1000 K. Affected by this, the r G m o presents positive within the temperature range from 1500 to 298.15 K, as shown in Figure 1B.
Figure 1. Dibenzofuran formation paths from phenanthrene (Phe); (A) Reaction scheme, enthalpies were calculated at 1000 K, and energy unit is kcal/mol; (B) Standard Gibbs energy changes for the overall reactions; (C) Activation Gibbs energies for the overall reactions (We accumulate the Gibbs energy change of each elementary reaction in turn, and the maximum value in a reaction route at given temperature is selected for the G ).
Figure 1. Dibenzofuran formation paths from phenanthrene (Phe); (A) Reaction scheme, enthalpies were calculated at 1000 K, and energy unit is kcal/mol; (B) Standard Gibbs energy changes for the overall reactions; (C) Activation Gibbs energies for the overall reactions (We accumulate the Gibbs energy change of each elementary reaction in turn, and the maximum value in a reaction route at given temperature is selected for the G ).
Ijms 16 05271 g001
The reaction mechanism of ClO∙ adding to C8a of Phe is similar to that of ∙OH. The five elementary reactions are addition, ring-opening, Cl-migration, cyclization and C2ClH2-leaving in turn. The ring‑opening and cyclization reactions have higher potential barriers and more absorbed heat than the equivalent elementary reactions in the OH-Phe addition process. However, the Cl-migration and C2ClH2‑leaving reactions release more heat, and then the enthalpy change for the overall reaction drops into negative, for instance, −20.50 kcal/mol at 1000 K. The r G m o drop correspondingly into negative, i.e., about −25.00 kcal/mol in temperature range 1500–298.15 K.
The removal of acetylene is also a possible reaction of Phe-B2, with the exception of H-migration. The reaction enthalpy for this decomposition reaches up to 43.18 kcal/mol. The r G m o will drop slowly into negative only if the temperature rises above 1100 K, and this suggests that the removal of acetylene occurs mainly in high temperatures. The subsequent cyclization-dehydrogenation reaction has been proved to be an elementary reaction, the intrinsic reaction coordinate (IRC) analysis of the reaction Phe-B5 → DF + ∙H is displayed on the Figure S1 (attached as supplementary materials). This reaction also has a high potential barrier; therefore, the importance of this reaction channel is watered down again. In Figure 1B, the r G m o for the overall reaction (Phe + ∙OH → DF + C2H2 + ∙H) is always greater than 14.39 kcal/mol. Figure 1C presents the variation tendencies of activation Gibbs energies for the overall reactions along with temperature, and we can see that the activation Gibbs energy always greater than 87.40 kcal/mol. Notwithstanding the removal of acetylene from Phe-C2 is exactly alike to that of Phe-B2, the overall reaction (Phe + ∙OCl → DF + C2H2 + ∙Cl) has a negative r G m o owing to the substantial reduce in Gibbs energy during the cyclization and Cl-leaving reaction.
Adding to C8a of Phe has a chance to create DF, but not a superior route. For example, the reaction enthalpy of ∙OH adding to C9 (or C10) of Phe amounts to −22.43 kcal/mol at 1000 K, which is obviously lower than the value of adding to C8a. Phe-D1 will be produced without any transition state. Similarly, H‑adding reaction also has no potential barrier, and Phe-A1 stands for the H-Phe adduct in Figure 1A. The unpaired electron of Phe-A1 mainly locates in the C10, and further oxidation occurs between C10 and O2, the greatest abundant oxidant in flame. The standard enthalpy change of C–O bonding is −22.06 kcal/mol at 1000 K, and this adding process has no transient state. The outer oxygen atom of Phe-A2 may bond with the adjacent unsaturated carbon atom, and an intermediate—Phe‑A3—with a four-membered ring appears. The elementary reaction is an endothermic process, resulting from the tensile force of the ring. On the other hand, Phe-A3 tends to break the O-O and C-C bonds in the four-membered ring simultaneously, and the driving force of bond breaking is also the tensile force. Phe‑A4 is a more stable state comparing with Phe-A2, and then r H m o from Phe-A2 to Phe-A4 reaches −42.15 kcal/mol at 1000 K. DF will be finally formed via cyclization and CH2CHO-leaving, which are very similar to the equivalent reactions in the pathway Phe-B. The cis/trans isomerism will be found in the reaction process, when O2 adds to Phe-D1. The energy data of trans-addition are placed above the arrows, in Figure 1A, and the data of cis-addition are located under the arrows. The subsequent reactions resemble to reaction pathway Phe-A mentioned above, except for the minor numerical difference.
The standard Gibbs free energy changes for the overall reaction from Phe to DF are calculated and presented in Figure 1B. The activation Gibbs energies for the overall reactions are plotted against temperature and shown in Figure 1C. The detailed data for every elementary reaction are attached in Tables S1 and S2. There is little chance that ∙OH oxides Phe solely into DF, but ClO∙ has more chances. There will be opportunity for O2 to take part in and change Phe into DF, once C9 of Phe is saturated. The radicals, which can bond with C9 and facilitate the combination between the radical-Phe adduct and O2, are not restricted to ∙H and ∙OH, i.e., the substituent on C9 has a minimal impact on the transformation process. As shown in Figure 1A, the activation enthalpies for each elementary reactions in channel Phe-A and Phe-D are smaller than 37.63 kcal/mol, and G are smaller than 41.50 kcal/mol within the temperature range 1500–298.15 K. The r G m o for overall reactions (Phe + ∙H + O2 → DF + ∙C2H3O and Phe + ∙OH + O2 → DF + ∙C2H3O2) are below −48.85 kcal/mol. So the oxidation routes—Phe-A and Phe-D—are superior scenarios than Phe-B and Phe-C.

2.2. Reaction from Fluorene to Dibenzofuran and 9-Fluorenone

Our calculations show that the bond length for C8a–C9 in a Flu molecule is 1.504 Å, and 1.437 Å is the length in a Phe molecule. The C8a–C9 bond in Flu is a typical single bond, in comparison with the experimental bond lengths for C–C and C=C (that is, 1.536 Å in ethane and 1.339 Å in ethylene [33]), while the C8a–C9 bond in Phe appears a hybrid between single and double bond. This difference in bonding properties makes an impact on the breaking of C8a–C9 bond. The standard enthalpy change for breaking C8a–C9 bond in Flu amounts to 129.89 kcal/mol at 1000 K, and this is significantly smaller than the energy change for Phe (191.84 kcal/mol). In addition, the activation enthalpy of ∙OH adding to C8a of Flu is 6.13 kcal/mol lower than the value of Phe; and the reaction enthalpy is 13.07 kcal/mol lower than that of Phe at 1000 K. In summary, ∙OH becomes much easier to bond with C8a of Flu, and the same goes for ClO∙ addition. As shown in Figure 2A, the ring-opening reactions of Flu-A1 and Flu-B1 occur much easier, for the bond energy of C8a–C9 in Flu seems much lower than that in Phe. The subsequent elementary reactions for Phe–B2 are H-migration, cyclization and CH3-leaving, while the subsequent elementary reactions for Phe–C2 are Cl–migration, cyclization and CClH2-leaving. Throughout the elementary reactions from Flu and ∙OH to DF and ∙CH3, the activation enthalpies are always less than 27.38 kcal/mol. The r G m o for the overall reaction are between −28.88 and −32.68 kcal/mol within the investigated temperatures. Therefore, the reaction (Flu + ∙OH → DF + ∙CH3) is a kinetically and thermodynamically feasible change. The standard Gibbs energy change for the reaction (Flu + ClO∙ → DF + ∙CClH2) is at least 21.63 kcal/mol lower than the r G m o for Flu + ∙OH → DF + ∙CH3, i.e., ClO∙ is much easier to oxidize Flu into DF. Figure 2B,C show the r G m o and G for each overall reaction, respectively. The calculated H , r H m o , G and r G m o for every elementary reaction are listed in Tables S3 and S4.
Besides addition, H-abstraction is also a feasible route to form DF. Flu has a carbon atom with sp3 hybrids unlike most PAHs, and therefore, this atom numbered 9 in Flu cannot participate in delocalized π bonding. If one of the hydrogen atoms bonding with C9 is dislodged, the hybridized type of C9 will be sp2, and the unpaired electron will participate in delocalized π bonding accordingly. The H-abstraction reaction take place easily under the influence of delocalization energy. The comparison has been made among the H-abstractions, in which ∙OH, HOO∙, ∙H, ∙CH3 and O2 are used to eliminate the hydrogen atom bonding with C9 of Phe, Flu and 9MF, and the thermodynamic parameters of these reactions at 1000 K are listed in Table 1. The activation enthalpies for the H-abstractions of Flu are always less than those of Phe, additionally reaction enthalpies, activation Gibbs energies and reaction Gibbs energies under normal conditions show the same regularity. This means that Flu is more likely to lose the hydrogen atom than Phe.
The product of H-abstraction—Flu–C1—can be oxidized by oxygen molecule. The first step of oxidation occurs between Flu–C1 and oxygen molecule, and this addition reaction is an exothermic process without a transition state. This bonding creates a peroxide radical—Flu–C2, and then an intermediate with a four-membered ring like dioxetane will be obtained if the outer oxygen atom is linked with C8a. Though the formation of this dioxetane-like intermediate is an endothermic process, dioxetane-ring-opening reaction releases a lot of energy, and gives rise to Flu-C4. As the sum total of heat effect from Flu–C2 to Flu–C4, 49.95 kcal/mol of heat will be released. Then the cyclization reaction accompanying with the CHO-leaving reaction generates DF finally.
Figure 2. Dibenzofuran and 9-fluorenone formation paths from fluorene (Flu); (A) Reaction scheme, enthalpies were calculated at 1000 K, and energy unit is kcal/mol; (B) Standard Gibbs energy changes for the overall reactions; (C) Activation Gibbs energies for the overall reactions.
Figure 2. Dibenzofuran and 9-fluorenone formation paths from fluorene (Flu); (A) Reaction scheme, enthalpies were calculated at 1000 K, and energy unit is kcal/mol; (B) Standard Gibbs energy changes for the overall reactions; (C) Activation Gibbs energies for the overall reactions.
Ijms 16 05271 g002
In addition to DF, 9FO is also the potential oxidation product of Flu. The loss of hydroxyl radical will change Flu–C2 into 9FO. The activation enthalpy and standard enthalpy change for this reaction at 1000 K are 38.29 and −30.28 kcal/mol, respectively. The extent of this decomposition reaction is great, for the r G m o at 1000 K reaches −64.05 kcal/mol. Another formation channel of 9FO is the reaction between Flu and O2. As mentioned above, one of the hydrogen atoms which bond with C9 of Flu can be abstracted by O2, and give rise to Flu-C1 and ∙OOH. It is worth mentioning that the coupling of these two radicals is accompanied by the spin-flip of one unpaired electron. There are two unpaired electrons with same direction of spin in O2, and when the H-abstraction occurs, the two unpaired electrons dispersed in Flu–C1 and ∙OOH will have the same spin direction. The spin-flip will occur when Flu-C1 and ∙OOH collide with each other, so that the two unpaired electrons can form a covalent bond effectively. The released heat of the spin-flip and bonding process amounts to 49.81 kcal/mol at 1000 K, and this entropy reducing reaction remains to be a spontaneous process below 1200 K, benefiting from reaction heat. The newly generated Flu-D1 has a strong tendency to decompose into 9FO and water, as can be seen from r G m o , e.g., −100.38 kcal/mol at 1000 K. The overall reaction (Flu + O2 → 9FO + H2O) is a bimolecular reaction, thus the activation Gibbs energy is higher than those reactions initiated by radicals, as shown in Figure 2C. The formation of 9FO should not be ignored for some reasons, such as that the concentration of O2 is much higher than radicals, and that there are only three elementary reactions.
Table 1. Comparison of the C9-H abstraction reactions at 1000 K (energy unit is kcal/mol).
Table 1. Comparison of the C9-H abstraction reactions at 1000 K (energy unit is kcal/mol).
H-Abstraction Reactions H r H m o G r G m o
∙OHPhenanthrene3.65−6.2228.79−12.50
Fluorene−0.24−35.3824.01−38.48
9-Methylfluorene−1.38−39.3724.69−48.04
HOO∙Phenanthrene25.8626.2253.2322.81
Fluorene12.12−2.9344.63−3.67
9-Methylfluorene10.39−6.9241.79−12.73
∙HPhenanthrene14.968.5835.541.97
Fluorene4.93−20.5828.27−24.01
9-Methylfluorene3.24−24.5725.71−33.57
∙CH3Phenanthrene17.266.6945.506.58
Fluorene6.27−22.4743.53−19.40
9-Methylfluorene6.75−26.4634.63−28.96
O2Phenanthrene58.3463.5578.9454.91
Fluorene36.5934.4062.9028.93
9-Methylfluorene33.9730.4160.2119.37

2.3. Reaction from 9-Methylfluorene to Dibenzofuran and 9-Fluorenone

Under the influence of methyl, the length of C8a–C9 bond in 9MF is 1.515 Å, which is slightly longer than the length in Flu—1.512 Å. Subsequent oxidation reactions of 9MF take place more easily, as shown in Figure S2. The calculated H , r H m o , G and r G m o for every elementary reaction are listed in Tables S5 and S6. In general, the elementary reactions of 9MF oxidation have slightly lower energy barriers and release a little more heat than those of Flu oxidations. In order to reduce in symmetry, the process of ∙OH and ClO∙ adding to 9MF have cis-trans isomers, as shown in Figure S2A. There are significant differences in forming 9FO: 9MF-C2 cannot eliminate CH3O∙ directly, and CH3OH cannot break away from 9MF-D1. In other words, 9FO will be formed when ∙OH and ∙CH3 depart from 9MF-D1 in turn. The elementary reaction, which ∙OH separates itself from 9MF-D1, breaks the O–O bond, and therefore it is a barrierless process with the gradually increasing enthalpy. The reaction enthalpy of OH-leaving reaches 44.11 kcal/mol at 1000 K. The CH3-leaving is relatively easy to accomplish, for the activation enthalpy is only 10.98 kcal/mol, and the reaction Gibbs energy is −36.84 kcal/mol at 1000 K.

2.4. Reaction from 9-Fluorenone to Dibenzofuran

The intermediate 9FO-B1 will be achieved if ∙OH adds to C8a of 9FO, as shown in Figure 3A. This elementary reaction has a low activation enthalpy (1.39 kcal/mol), and the releases heat is 15.44 kcal/mol at 1000 K. The ring-opening reaction of 9FO–B1 is easy to complete in comparison with the equivalent reactions of Flu and 9MF, because the activation enthalpy is significantly decreased (15.85 kcal/mol for 9FO, 26.82 kcal/mol for Flu, and 23.14 kcal/mol for 9MF). There are two possible reaction pathways for 9FO–B2: (I) the H-migration creates Flu–C4 (an intermediate mentioned in the transformation route from Flu to DF); (II) carbon monoxide may slip the leash and leaves Phe–C5 alone. Each of the competing reactions has its advantages and disadvantages. When activation enthalpy and reaction enthalpy are taken into account, the change from 9FO–B2 to Flu–C4 will be favorable; but in terms of standard reaction Gibbs energy, Phe–C5 and CO are the dominant products. The point of view we take is that Phe–C5 is more favorable at higher temperatures and Flu–C4 at lower temperatures. The calculated H , r H m o , G and r G m o for every elementary reaction at different temperatures are listed in Tables S7 and S8.
So take the other case, ∙OH bonds with C9 of 9FO. The first elementary reaction is to form a C–O bond. The bonding process has a low activation enthalpy—8.16 kcal/mol, and the released heat is 13.38 kcal/mol at 1000 K. The subsequent reaction is the breakage of C8a–C9 bond, and the activation enthalpy of this breakage—9.23 kcal/mol—is lowest among all the cases mentioned in this study. Moreover, then there are several routes to change 9FO-A2 into DF, such as O(3P) adding to 9FO–A2 (as shown in Figure 3A), O2 adding to 9FO–A2 and H-migration of 9FO-A2 (Figure S3A). In the case of O(3P) adding, the standard Gibbs energy changes for the overall reaction are significantly lower than those of ∙OH adding to C8a at any temperature between 1500 and 298.15 K, as shown in Figure 3B.
The reaction mechanism of ClO∙ addition is similar to those of ∙OH addition. The bonding between ClO∙ and C8a and the breaking of C8a–C9 bond will be a little more difficult to occur than the reactions of ∙OH, when thought of a little higher activation enthalpies and reaction enthalpies. However, the Cl-migration reaction is much more feasible than H-migration. Under standard conditions and 1000 K, the activation enthalpies for Cl- and H-migration are 9.48 and 13.76 kcal/mol, respectively; the reaction enthalpies are −53.86 and −1.56 kcal/mol; the activation Gibbs energies are 13.96 and 22.25 kcal/mol; and the reaction Gibbs energies are −51.45 and 1.05 kcal/mol. This means that the superiority of the ClO∙ addition exists similarly in the DF formation from 9FO. In the reaction channel of ClO∙ adding to C9, the possible reactions of intermediate 9FO-D2 are O(3P)addition (as shown in Figure 3A), O2 addition and Cl-migration (Figure S3B). The r G m o   of the overall reaction (9FO + ∙OCl + O(3P) → DF + CO2 + ∙Cl) is lower than −171.73 kcal/mol, and such a low Gibbs energy change shows that the reaction extant is great.
Figure 3. Dibenzofuran formation paths from 9-fluorenone (9FO); (A) Reaction scheme, enthalpies were calculated at 1000 K, and energy unit is kcal/mol; (B) Standard Gibbs energy changes for the overall reactions; (C) Activation Gibbs energies for the overall reactions.
Figure 3. Dibenzofuran formation paths from 9-fluorenone (9FO); (A) Reaction scheme, enthalpies were calculated at 1000 K, and energy unit is kcal/mol; (B) Standard Gibbs energy changes for the overall reactions; (C) Activation Gibbs energies for the overall reactions.
Ijms 16 05271 g003

3. Experimental Section

The high-accuracy quantum chemical calculations were performed with Gaussian 09 programs [34], using the method based on the density functional theory including dispersion correction (DFT-D). DFT-D is an atomic pairwise correction that is added to the result of a converged standard DFT calculation [35]. The D3 version of Grimme’s dispersion (DFT-D3) has been refined regarding higher accuracy, broader range of applicability, and less empiricism [36,37]. As one of the DFT-D3 methods, the M06‑2X‑GD3 provides anisotropy of van der Waals interactions that agrees best with both higher level of theory and experiments [38,39]._ENREF_24 All the geometrical conformations and vibration parameters were calculated on the M06‑2X‑GD3/6-311+G(d,p) level using Gaussian 09. In addition, the corrected factors of enthalpies and Gibbs energies (HCF, GCF) were also acquired through the Freq (vibrational frequencies) calculations on the M06‑2X‑GD3/6-311+G(d,p) level, by means of specifying the temperature at which substances inhabit. These calculations carried out at every 100 K from 1500 to 298.15 K. The M06‑2X‑GD3/aug-cc-pVQZ was employed to calculate the single point energies (ESP) of various species. The final enthalpies (H) at 1000 K were worked out by the expression—H1000 K = ESP + HCF, 1000 K. Analogously, all the H and G at different temperature were figured out. Each transition state was verified to connect the related reactants and products by performing an intrinsic reaction coordinate (IRC) analysis.

4. Conclusions

Calculations show that Phe, Flu, 9MF and 9FO can be oxidized into DF by ∙OH, ClO∙, O2 or O(3P) in flame or flue gas. The C8a–C9 bond could be broken when ∙OH or ClO∙ adds to C8a of these molecules. The middle ring is broken most easily by 9FO, followed by 9MF, Flu, and, lastly, Phe. In general, the addition and ring-opening reactions of ClO∙ are a little more difficult to accomplish compared with those of ∙OH, but the Gibbs energy changes of the overall reactions show that the ClO∙ addition are dominant channels for the r G m o are at least 21.63 kcal/mol lower than the value for the equivalent reactions initiated by ∙OH.
Adding O2 to C10 of the radical-Phe adduct, whose C9 has been saturated by ∙H or ∙OH, is a more feasible oxidation route than the channel Phe-B (Phe + ∙OH → DF + ∙C2H3). Formation paths of 9FO are explored during the oxidation of Flu and 9MF, and the results show that the extents of producing 9FO are great, and yet the energy barriers are higher than those of creating DF. DF may also stem from the fact that ∙OH or ClO∙ adds to the C9 of 9FO, and it is worth mentioning that (I) an additional oxidant—molecular or atomic oxygen—is needed and (II) more heat is released in this reaction route.

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/1422-0067/16/03/5271/s1.

Acknowledgments

This work was supported by NNSFC (National Natural Science Foundation of China, Project Nos. 21337001, 21177077).

Author Contributions

Shanqing Li conceived the article, performed the calculations. Shanqing Li and Qingzhu Zhang wrote the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Arena, U. Process and technological aspects of municipal solid waste gasification. A review. Waste Manag. 2012, 32, 625–639. [Google Scholar] [CrossRef] [PubMed]
  2. McKay, G. Dioxin characterisation, formation and minimisation during municipal solid waste (MSW) incineration: Review. Chem. Eng. J. 2002, 86, 343–368. [Google Scholar] [CrossRef]
  3. Fernandez-Salguero, P.; Pineau, T.; Hilbert, D.M.; McPhail, T.; Lee, S.S.T.; Kimura, S.; Nebert, D.W.; Rudikoff, S.; Ward, J.M.; Gonzalez, F.J. Immune system impairment and hepatic fibrosis in mice lacking the dioxin-binding Ah receptor. Science 1995, 268, 722–726. [Google Scholar] [CrossRef] [PubMed]
  4. Van den Berg, M.; Birnbaum, L.S.; Denison, M.; De Vito, M.; Farland, W.; Feeley, M.; Fiedler, H.; Hakansson, H.; Hanberg, A.; Haws, L.; et al. The 2005 world health organization reevaluation of human and mammalian toxic equivalency factors for dioxins and dioxin-like compounds. Toxicol. Sci. 2006, 93, 223–241. [Google Scholar] [CrossRef] [PubMed]
  5. Altarawneh, M.; Dlugogorski, B.Z.; Kennedy, E.M.; Mackie, J.C. Mechanisms for formation, chlorination, dechlorination and destruction of polychlorinated dibenzo-p-dioxins and dibenzofurans (PCDD/Fs). Prog. Energ. Combust. 2009, 35, 245–274. [Google Scholar] [CrossRef]
  6. Ryu, J.Y.; Mulholland, J.A.; Dunn, J.E.; Iino, F.; Gullett, B.K. Potential role of chlorination pathways in PCDD/F formation in a municipal waste incinerator. Environ. Sci. Technol. 2004, 38, 5112–5119. [Google Scholar] [CrossRef] [PubMed]
  7. Lomnicki, S.; Dellinger, B. A detailed mechanism of the surface-mediated formation of PCDD/F from the oxidation of 2-chlorophenol on a CuO/silica surface. J. Phys. Chem. A 2003, 107, 4387–4395. [Google Scholar] [CrossRef]
  8. Tuppurainen, K.A.; Ruokojärvi, P.H.; Asikainen, A.H.; Aatamila, M.; Ruuskanen, J. Chlorophenols as precursors of PCDD/Fs in incineration processes:  Correlations, PLS modeling, and reaction mechanisms. Environ. Sci. Technol. 2000, 34, 4958–4962. [Google Scholar] [CrossRef]
  9. Ryu, J.-Y.; Choi, K.-C.; Mulholland, J.A. Polychlorinated dibenzo-p-dioxin (PCDD) and dibenzofurna (PCDF) isomer patterns from municipal waste combustion: Formation mechanism fingerprints. Chemosphere 2006, 65, 1526–1536. [Google Scholar] [CrossRef] [PubMed]
  10. Altarawneh, M.; Dlugogorski, B.Z.; Kennedy, E.M.; Mackie, J.C. Quantum chemical investigation of formation of polychlorodibenzo-p-dioxins and dibenzofurans from oxidation and pyrolysis of 2‑chlorophenol. J. Phys. Chem. A 2007, 111, 2563–2573. [Google Scholar] [CrossRef] [PubMed]
  11. Briois, C.; Visez, N.; Baillet, C.; Sawerysyn, J.P. Experimental study on the thermal oxidation of 2-chlorophenol in air over the temperature range 450–900 degrees C. Chemosphere 2006, 62, 1806–1816. [Google Scholar] [CrossRef] [PubMed]
  12. Altarawneh, M.; Dlugogorski, B.Z.; Kennedy, E.M.; Mackie, J.C. Theoretical study of unimolecular decomposition of catechol. J. Phys. Chem. A 2010, 114, 1060–1067. [Google Scholar] [CrossRef] [PubMed]
  13. Alsoufi, A.; Altarawneh, M.; Dlugogorski, B.Z.; Kennedy, E.M.; Mackie, J.C. A DFT study on the self-coupling reactions of the three isomeric semiquinone radicals. J. Mol. Struct. 2010, 958, 106–115. [Google Scholar] [CrossRef]
  14. Chen, K.; Wojtalewicz, D.; Altarawneh, M.; Mackie, J.C.; Kennedy, E.M.; Dlugogorski, B.Z. Formation of polychlorinated dibenzo-p-dioxins and dibenzofurans (PCDD/F) in oxidation of captan pesticide. Proc. Combust. Inst. 2011, 33, 701–708. [Google Scholar] [CrossRef]
  15. Cosentino, U.; Pitea, D.; Moro, G. Computational modelling of de novo synthesis of dibenzofuran: Oxidative pathways of pyrene and benzodibenzofuran. Theor. Chem. Acc. 2012, 131, 1–12. [Google Scholar] [CrossRef]
  16. Zimmermann, R.; Blumenstock, M.; Heger, H.J.; Schramm, K.W.; Kettrup, A. Emission of nonchlorinated and chlorinated aromatics in the flue gas of incineration plants during and after transient disturbances of combustion conditions:  Delayed emission effects. Environ. Sci. Technol. 2001, 35, 1019–1030. [Google Scholar] [CrossRef] [PubMed]
  17. Rivera-Austrui, J.; Borrajo, M.A.; Martinez, K.; Adrados, M.A.; Abalos, M.; Van Bavel, B.; Rivera, J.; Abad, E. Assessment of polychlorinated dibenzo-p-dioxin and dibenzofuran emissions from a hazardous waste incineration plant using long-term sampling equipment. Chemosphere 2011, 82, 1343–1349. [Google Scholar] [CrossRef] [PubMed]
  18. Onofrio, M.; Spataro, R.; Botta, S. The role of a steel plant in north-west Italy to the local air concentrations of PCDD/Fs. Chemosphere 2011, 82, 708–717. [Google Scholar] [CrossRef] [PubMed]
  19. Hell, K.; Stieglitz, L.; Dinjus, E. Mechanistic aspects of the de-novo synthesis of PCDD/PCDF on model mixtures and MSWI fly ashes using amorphous 12C- and 13C-labeled carbon. Environ. Sci. Technol. 2001, 35, 3892–3898. [Google Scholar] [CrossRef] [PubMed]
  20. Yan, M.; Li, X.; Zhang, X.; Liu, K.; Yan, J.; Yan, K. Correlation between PAHS and PCDD/Fs in municipal solid waste incinerators. J. Zhejiang Univ. (Eng. Sci.) 2010, 44, 1118–1121. [Google Scholar]
  21. Johansson, I.; van Bavel, B. Levels and patterns of polycyclic aromatic hydrocarbons in incineration ashes. Sci. Total Environ. 2003, 311, 221–231. [Google Scholar] [CrossRef] [PubMed]
  22. Weber, R.; Iino, F.; Imagawa, T.; Takeuchi, M.; Sakurai, T.; Sadakata, M. Formation of PCDF, PCDD, PCB, and PCN in de novo synthesis from PAH: Mechanistic aspects and correlation to fluidized bed incinerators. Chemosphere 2001, 44, 1429–1438. [Google Scholar] [CrossRef] [PubMed]
  23. Lee, W.J.; Liow, M.C.; Tsai, P.J.; Hsieh, L.T. Emission of polycyclic aromatic hydrocarbons from medical waste incinerators. Atmos. Environ. 2002, 36, 781–790. [Google Scholar] [CrossRef]
  24. Ryu, J.-Y.; Mulholland, J.A.; Chu, B. Chlorination of dibenzofuran and dibenzo-p-dioxin vapor by copper (II) chloride. Chemosphere 2003, 51, 1031–1039. [Google Scholar] [CrossRef] [PubMed]
  25. Blumenstock, M.; Zimmermann, R.; Schramm, K.W.; Kettrup, A. Influence of combustion conditions on the PCDD/F-, PCB-, PCBz- and PAH-concentrations in the post-combustion chamber of a waste incineration pilot plant. Chemosphere 2000, 40, 987–993. [Google Scholar] [CrossRef] [PubMed]
  26. Tokmakov, I.V.; Kim, G.-S.; Kislov, V.V.; Mebel, A.M.; Lin, M.C. The reaction of phenyl radical with molecular oxygen:  A G2M study of the potential energy surface. J. Phys. Chem. A 2005, 109, 6114–6127. [Google Scholar] [CrossRef] [PubMed]
  27. Sebbar, N.; Bockhorn, H.; Bozzelli, J.W. Thermochemical similarities among three reaction systems: Vinyl + O2–phenyl + O2–dibenzofuranyl + O2. Combust. Sci. Technol. 2008, 180, 959–974. [Google Scholar] [CrossRef]
  28. Li, Y.; Zhang, L.; Tian, Z.; Yuan, T.; Wang, J.; Yang, B.; Qi, F. Experimental study of a fuel-rich premixed toluene flame at low pressure. Energ. Fuel 2009, 23, 1473–1485. [Google Scholar] [CrossRef]
  29. Khachatryan, L.; Burcat, A.; Dellinger, B. The role of chlorine atoms and hydroxyl radicals in the formation of PCDDs from the oxidative pyrolysis of 2,4,6-trichlorophenol. Int. J. Chem. Kinet. 2010, 42, 90–97. [Google Scholar] [CrossRef]
  30. Aizawa, T. Diode-laser wavelength-modulation absorption spectroscopy for quantitative in situ measurements of temperature and OH radical concentration in combustion gases. Appl. Opt. 2001, 40, 4894–4903. [Google Scholar] [CrossRef] [PubMed]
  31. Yuasa, S.; Yushina, S.; Uchida, T.; Shiraishi, N. Ignition and combustion of ammonium perchlorate in a hydrogen atmosphere. Proc. Combust. Inst. 2000, 28, 863–870. [Google Scholar] [CrossRef]
  32. Canneaux, S.; Hammaecher, C.; Cours, T.; Louis, F.; Ribaucour, M. Theoretical study of H‑abstraction reactions from CH3Cl and CH3Br molecules by ClO and BrO radicals. J. Phys. Chem. A 2012, 116, 4396–4408. [Google Scholar] [CrossRef] [PubMed]
  33. Johnson III, R.D. NIST computational chemistry comparison and benchmark database. Release 16 August 2013. Available online: http://cccbdb.nist.gov/ (accessed on 20 December 2014).
  34. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G.A.; et al. Gaussian 09; Gaussian, Inc.: Wallingford, CT, USA, 2009. [Google Scholar]
  35. Goerigk, L. How do DFT-DCP, DFT-NL, and DFT-D3 compare for the description of London-dispersion effects in conformers and general thermochemistry? J. Chem. Theory Comput. 2014, 10, 968–980. [Google Scholar] [CrossRef]
  36. Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu. J. Chem. Phys. 2010, 132. [Google Scholar] [CrossRef]
  37. Grimme, S.; Ehrlich, S.; Goerigk, L. Effect of the damping function in dispersion corrected density functional theory. J. Comput. Chem. 2011, 32, 1456–1465. [Google Scholar] [CrossRef] [PubMed]
  38. Steinmann, S.N.; Corminboeuf, C. Comprehensive benchmarking of a density-dependent dispersion correction. J. Chem. Theory Comput. 2011, 7, 3567–3577. [Google Scholar] [CrossRef]
  39. Goerigk, L.; Grimme, S. A thorough benchmark of density functional methods for general main group thermochemistry, kinetics, and noncovalent interactions. Phys. Chem. Chem. Phys. 2011, 13, 6670–6688. [Google Scholar] [CrossRef] [PubMed]

Share and Cite

MDPI and ACS Style

Li, S.; Zhang, Q. Mechanistic Studies on the Dibenzofuran Formation from Phenanthrene, Fluorene and 9–Fluorenone. Int. J. Mol. Sci. 2015, 16, 5271-5284. https://doi.org/10.3390/ijms16035271

AMA Style

Li S, Zhang Q. Mechanistic Studies on the Dibenzofuran Formation from Phenanthrene, Fluorene and 9–Fluorenone. International Journal of Molecular Sciences. 2015; 16(3):5271-5284. https://doi.org/10.3390/ijms16035271

Chicago/Turabian Style

Li, Shanqing, and Qingzhu Zhang. 2015. "Mechanistic Studies on the Dibenzofuran Formation from Phenanthrene, Fluorene and 9–Fluorenone" International Journal of Molecular Sciences 16, no. 3: 5271-5284. https://doi.org/10.3390/ijms16035271

APA Style

Li, S., & Zhang, Q. (2015). Mechanistic Studies on the Dibenzofuran Formation from Phenanthrene, Fluorene and 9–Fluorenone. International Journal of Molecular Sciences, 16(3), 5271-5284. https://doi.org/10.3390/ijms16035271

Article Metrics

Back to TopTop