ORIGINAL RESEARCH article

Front. Chem., 01 September 2021

Sec. Green and Sustainable Chemistry

Volume 9 - 2021 | https://doi.org/10.3389/fchem.2021.580643

Efficiency and Quantitative Structure-Activity Relationship of Monoaromatics Oxidation by Quinone-Activated Persulfate

  • 1. State Environmental Protection Key Laboratory of Soil Environmental Management and Pollution Control, Nanjing Institute of Environmental Sciences, Ministry of Ecology and Environment, Nanjing, China

  • 2. College of Environment, Hohai University, Nanjing, China

Abstract

Quinones and quinone-containing organics have potential of activating persulfate to produce sulfate radical. In this work, the optimal condition for quinone activation of persulfate was investigated. It was found representative monoaromatics were degraded fastest in alkaline environment (pH 10.0), but excessive alkalinity restrained the reaction instead. The mechanisms to explain this phenomenon were speculated. The effect of initial quinone concentration on persulfate oxidation was also investigated at pH 10.0. In addition, a quantitative structure-activity relationship model was established with 15 kinds of monoaromatics, which revealed the most negative atomic net charges on carbon atom played an important role on degradation rates. Chemicals with a smaller were easier oxidized in quinone-activate system. This finding helps further exploration of effective activator in alkaline environment.

Introduction

Persulfate-based oxidation process is widely used in water and soil remediation, in which persulfate is activated to produce sulfate radical anion () as the main active material to degrade recalcitrant organics (Waldemer et al., ; Xian et al., ). The common activation methods include thermal energy, transition metal, alkali and UV radiation, etc. Recently, new activating methods which utilize naturally occurring materials such as metal ores and organics have attracted much attention (Lin et al., ; Wu et al., ; Zhang et al., ). Quinones and organic matters containing quinone-like groups are ubiquitous in natural environment (Siqueira et al., ; Cory and McKnight, ), and they have been found to be effective as the accelerator or activator in advanced oxidation process (Zhou et al., ; Zhang et al., ). For instance, quinones increased the oxidation rate of contaminants in Fenton-like system, which was ascribed to their role as an electron shuttle (Duesterberg and Waite, ). The presence of quinone contributed to build up semiquinone/quinone and Fe(III)/Fe(II) cycles (Ma et al., ). Zhu et al. (, ) observed the production of hydroxyl radical (OH•) by halogenated quinones and hydrogen peroxide. Fang et al. () found that was generated from persulfate in the presence of 1,4-benzoquinone (BQ), 2-methyl-1,4-benzoquinone (MBQ), or 2-chloro-1,4-benzoquinone (MCQ). These reactions were all semiquinone radical dependent. Zhou et al. () revealed that 1,4-benzoquinone could efficiently activate peroxymonosulfate for the degradation of sulfamethoxazole, and the catalytic mechanism was that a dioxirane intermediate formed and was decomposed into 1O2 subsequently.

Despite those advantages, there are still some unsolved confusions, including the most suitable target chemicals and the optimal working condition of quinone-activated persulfate system. As a result, the influences of pH and quinone concentration on the degradation efficiency of representative volatile organic compound were investigated in this study, and a quantitative structure-activity relationship (QSAR) model for the degradation rate constant (k) of monoaromatics was established. This comprehensive work would further reveal the quinone activation mechanisms and provide data support for the application of this oxidation method.

Methodology

Chemicals

Sodium persulfate (Na2S2O8), target contaminants (listed in Table 1), sodium thiosulfate (Na2S2O3), sulfuric acid (H2SO4), sodium hydroxide (NaOH), sodium hydrogen phosphate (Na2HPO4), sodium dihydrogen phosphate (NaH2PO4), and sodium borate (Na2B4O7·10H2O) were obtained from Sinopharm Chemical Reagent Co., Ltd. The BQ, MBQ, 2,5-dimethyl-1,4-benzoquinone (DMBQ) and the spin-trapping agent 5,5-dimethyl-1-pyrroline-N-oxide (DMPO) were purchased from Aladdin Chemistry Co., Ltd. All reagents were of high purity analytic grade. Solutions in this study were prepared with ultrapure water purified by Millipore deionized water purification system.

Table 1

Microwave frequency9.77 GHz
Microwave power19.92 mW
Modulation frequency100 kHz
Modulation amplitude1.00 G
Time constant40.96 ms
Sweep time83.87 s

EPR determination condition.

Degradation Procedures

All degradation experiments were performed at 20°C in 40 mL glass vial reactors with Teflon inner membrane lids. The concentration of Na2S2O8 and monoaromatics were 5 mM and 100 μg/L, respectively, and the total volume was 40 mL. Samples were withdrawn and quenched with Na2S2O3 periodically before analysis (Huber et al., ; Zhou et al., ; Zhang et al., ). A preliminary study demonstrated that Na2S2O3 did not react with typical organics (1,2-dichloroethane, bromobenzene, and 1,2,4-trichlorobenzene) and had no interference with the analysis of the organics. All the experiments were conducted in triplicate.

In order to explore the optimal conditions of the reaction, different initial pH (pH 3, 7, 10, and 12) and quinone concentrations (0.01, 0.05, 0.1, 0.5, and 1 mM) were set respectively to investigate their effects on the degradation efficiency.

The reaction rate constants used to build the QSAR model were obtained from experiments under the optimal conditions acquired in the above investigation.

Analysis

The concentration of the monoaromatic was measured by purge & trap autosampler (Eclipse 4552&4660) and GC-MS system (Agilent 7890A/5975C) equipped with a capillary column (DB-624 60 m × 0.25mm × 1.4μm). The column temperature program on GC-MS was set as follows: the beginning temperature of the column oven was 40°C, held for 5 min, increased up to 200°C with a heating rate of 8°C/min and increased up to 230°C with a heating rate of 10°C/min, remained for 2 min. The injector temperature was 250°C and the split ratio was 50:1. The carrier gas was helium.

Electron paramagnetic resonance (EPR) is a spectral method for detecting unpaired electron and can be used to identify free radicals during the degradation process by the technique of spin trapping, which is conducted to produce spin adduct with the diamagnetic compound named spin trapping agent (Johnson et al., ; Kevan, ). EPR determination in this work was performed as follows: 450 μL samples were withdrawn and mixed with 50 μL DMPO (1 M), and then the solution was injected into a capillary tube, which was inserted into the EPR spectrometer (Bruker, EMX10/12) with a resonance frequency of 9.77 GHz, microwave power of 19.97 mW, modulation frequency of 100 kHz, modulation amplitude of 1.0 G, time constant of 40.96 ms, and sweep time of 83.87 s.

Model Derivation and Validation

The QSAR model helps understand the oxidation mechanism of quinone-activated persulfate based on the relation between degradation efficiency and molecular structure property. In this work, the structural descriptors (shown in Table 2) came from the manuscript by Zhu et al. (), in which the structure optimization and frequency calculation were conducted at B3LYP/6-31+G(d,p) level using density functional theory (DFT) method with Gaussian 09. The molecular descriptors included energy of the highest occupied molecular orbital (EHOMO), energy of the lowest unoccupied molecular orbital (ELUMO), the largest negative partial net charge on a carbon atom (), the largest positive net charge on a hydrogen atom (), dipole moment (μ), total energy (TE), surface area (SA), molecular volume (Vm), hydration energy (HE), refractivity (R), partition coefficient (log P) and polarizability (P). The total data set (n = 20) was divided into training set (n = 15) and test set (n = 5) at intervals. The chemicals were ordered according to their structures. The ones with the number of 4n (n = 1, 2, 3, 4, and 5) were classified as the test set for external validation, and the rest were set as the training set to derive the model. The test set chemicals were marked with asterisk (*) in Table 2.

Table 2

No.CompoundEHOMOELUMOμTESAVmHERlogPP
(eV)(eV)(e)(e)(debye)(Å2)(Å3)(10−30esu)
1benzene−6.7190.073−0.0840.0840.000−232.3237.4329.9−2.11026.062.05010.43
2toluene−6.4110.081−0.3810.1220.344−271.6268.4384.2−1.39031.102.51012.27
3ethylbenzene−6.3720.098−0.3280.1120.459−310.9290.3432.70.19035.702.91014.10
4*n-propylbenzene−6.4080.123−0.3170.1030.325−350.2323.7490.01.71040.303.31015.94
5isopropylbenzene−6.4250.098−0.3040.1070.284−350.2320.8480.81.77040.253.24015.94
6n-butylbenzene−6.3970.125−0.3180.1030.361−389.5360.3545.73.19044.903.70017.77
7tert-butylbenzene−6.4120.101−0.3100.1050.325−389.5332.3518.81.87044.723.67017.77
8*para-xylene−6.1420.173−0.3810.1200.096−310.9297.3436.9−0.72036.142.98014.10
9ortho-xylene−6.2480.195−0.3810.1180.585−310.9286.3427.8−0.77036.142.98014.10
101,3,5-trimethylbenzene−6.1480.154−0.3810.1200.088−350.2322.6489.8−0.12041.183.45015.94
111,2,4-trimethylbenzene−6.0140.203−0.3820.1190.290−350.2315.0480.6−0.15041.183.45015.94
12*para-isopropyltoluene−6.1560.100−0.3810.1210.092−389.5349.0533.52.40045.293.71017.77
13chlorobenzene−6.714−0.365−0.0930.1101.918−691.9264.5376.1−2.90030.862.56012.36
14bromobenzene−6.590−0.366−0.0860.1091.810−2803.4274.7394.7−2.89033.682.84013.06
15para-chlorotoluene−6.455−0.333−0.3820.1262.350−731.2292.6429.4−2.22035.903.03014.20
16*ortho-chlorotoluene−6.568−0.265−0.3780.1361.665−731.2285.2422.3−2.27035.903.03014.20
171,4-dichlorobenzene−6.750−0.764−0.0940.1200.000−1151.4288.2421.1−3.72035.673.08014.29
181,2-dichlorobenzene−6.848−0.673−0.0880.1182.761−1151.4285.2416.4−3.76035.673.08014.29
191,2,4-trichlorobenzene−6.926−1.037−0.0930.1391.401−1611.0310.0460.7−4.63040.473.60016.22
20*styrene−6.047−0.853−0.2520.1010.187−309.7283.5416.0−3.73035.742.70013.91

The molecular descriptors of 20 kinds of monoaromatics.

*

Test set chemicals.

The reaction kinetics were explored and the rate constants of 20 kinds of monoaromatics were calculated. The multiple stepwise linear regression (MSLR) analytical method was carried out by SPSS version 17.0 to build the QSAR model.

The quality of derived QSAR was evaluated with the squared regression coefficient (R2), the standard error of the estimate (SEE) and the Fisher's test value (F). The model internal predictivity, which reflected the stability of model, was measured by the leave-one-out cross validation coefficient () and the root mean square of the cross-validation (RMSEcv). The external validation capability was measured by external validation coefficient () and root mean square error of prediction (RMSEp). It is generally believed that RMSEcv < 1 and > 0.7 indicates the model has a great internal prediction ability, and that >0.4 and RMSEp < 1 indicates the model has a great external prediction ability (Gramatica, ; Puzyn et al., ).

Results and Discussion

Effect of pH

It has been reported that pH value had a significant effect on persulfate oxidation (Liu et al., ; Zhou et al., ). We chose an acid (pH 3), a neutral (pH 7), a basic (pH 10.0) and a strongly alkaline environment (pH 12) to investigate the effect of pH on 1,4-benzoquinone-activated persulfate oxidation efficiency. A total of 10 kinds of volatile organic compounds, including 4 kinds of monoaromatics (benzene, p-xylene, chlorobenzene, and p-chlorotoluene), 3 kinds of chloroalkanes (chloroform, 1,2-dichloroethane, and 1,2-dichloropropane), and 3 kinds of chloroalkenes (trans-1,2-dichloroethylene, trans-1,3-dichloropropylene, and trichloroethylene) were chosen as the target contaminants to investigate the effects of pH. The initial concentrations of Na2S2O8, 1,4-benzoquinone and target contaminant were 5 mM, 0.1 mM, and 100 μg/L, respectively. The pH of acid solution was adjusted with sulfuric acid, and reactions in neutral or alkaline environment occurred in the 10 mM phosphate or borate buffer solutions which had negligible influences on persulfate oxidation efficiency (Fang and Shang, ; Fang et al., ). The sample were taken and analyzed at 0, 10 min, 30 min, 1 h, 1.5 h, 2 h, 3 h, and 4 h, respectively to observe the variation of concentration with time. Controls at different pH with and without persulfate were set up respectively.

The 4-h degradation percentages of the target contaminants at different pH are shown in Figure 1A. The average data and standard deviations of parallel samples are presented. All the control samples showed much lower degradation ratios than test groups. It could be seen that monoaromatics exhibited the highest degradation level compared with chloroalkanes and chloroalkenes, which might be due to the selectivity of . It is generally accepted that reacts through three main pathways, including hydrogen atom abstraction, addition on unsaturated bonds, and electron transfer, among which electron transfer is of the most importance. It has been revealed that aromatic compounds were easily oxidized by , and the presence of chlorines was expected to slow oxidation rate (Ahmed et al., ; Xiao et al., ). Chlorinated saturated hydrocarbon was resistant to oxidation according to previous studies (Huang et al., ; Tsitonaki et al., ; Zhu et al., ). As for all the monoaromatics and most of the chlorohydrocarbon, the highest degradation ratio was obtained at pH 10.0. For example, chlorobenzene was degraded by about 96% at pH 10.0, while the degradation ratios were 64, 53, and 45% at pH 3.0, 7.0, and 12.0, respectively. In order to explore the reason, the EPR responses of the oxidation system at different pH were determined, and the results are shown in Figure 1B. The strongest signal of was observed at pH 10.0. The reason might be that on the one hand, the base acted as the activator of persulfate to produce , as shown Equations (1–3) (Furman et al., ). On the other hand, quinone reacted with HO2 to generate semiquinone radical, which accelerated the production of , as presented in Equations (4, 5). In addition, anionic form of phenols at alkaline pH was the activating species in persulfate systems (Ocampo, ; Ahmad et al., ), as shown in Equations (6, 7). But as the alkalinity became stronger, might react with the reductive radicals, thus weakened the oxidation effect, as shown in Equation (8). An additional experiment showed that carbon tetrachloride was degraded by 10.62 ± 1.03% during 4 h at pH 10.0 while the degradation ratio was negligible in the control without quinone. It verified the occurrence of reductive radicals in this oxidation system.

Figure 1

Effect of Initial Quinone Concentration

In this work, we compared the effect of BQ, MBQ, and DMBQ on the 4-h degradation rate of representative monoaromatics (benzene, p-xylene, chlorobenzene, and p-chlorotoluene) at pH 10.0. The initial quinone concentrations were set as 0.01, 0.05, 0.1, 0.5, and 1.0 mM, respectively, while the concentration of Na2S2O8 and monoaromatics were the same as above experiments. The results are shown in Figure 2, which indicated that BQ activated persulfate showed the strongest oxidation effect, followed by MBQ and DMBQ. The degradation rate of the monoaromatics increased with the quinone concentration, and reached a plateau at 0.1 or 0.5 mM of quinone, and then dropped gradually. When the amount of quinone was further increased, MBQ and DMBQ were difficult to be dissolved. The results were a little different from those observed by Fang et al. (). The amount of oxidative radicals reached the maximum when MBQ increased to 0.2 mM, and then decreased with a higher MBQ concentration in their observation. The difference might be attributed to the fact that the effective quinone concentration range was enlarged in alkaline condition.

Figure 2

Calculation of Kinetic Constant

As shown in Figure 1C, a rapid decrease of the monoaromatic concentration was observed in the first 2 h, followed by a slow reduction process, which might be related with the decreased amount of (Fang et al., ).

The degradation kinetics of 20 kinds of monoaromatics (100 μg/L) by 1,4-benzoquinone (0.1 mM) activated persulfate at pH 10.0 was investigated. The samples were collected at 0, 2, 5, 10, 15, and 20 min, respectively, and it was found that each of the oxidation process followed pseudo-first-order kinetics in the first 20 min. The pseudo-first-order kinetic model is shown in Equation (9), in which k is the pseudo-first-order reaction rate constant, t is the reactive time, and Ct is the concentration at t time (Shi et al., ; Xu et al., ).

According to the pseudo-first-order kinetic model, the reaction rate constant and R2 of 20 kinds of monoaromatics were calculated and shown in Table 3. The R2 showed that the values of k all fitted well with the pseudo-first-order kinetic model.

Table 3

No.Compoundk–logkR2
1benzene0.03671.4350.9551
2toluene0.08321.0800.9366
3ethylbenzene0.08891.0510.9735
4*n-propylbenzene0.10990.9590.9931
5isopropylbenzene0.08411.0750.9911
6n-butylbenzene0.07281.1380.8979
7tert-butylbenzene0.05521.2580.9866
8*para-xylene0.07781.1090.9854
9ortho-xylene0.11020.9580.9984
101,3,5-trimethylbenzene0.11610.9350.9999
111,2,4-trimethylbenzene0.09291.0320.9737
12*para-isopropyltoluene0.10170.9930.8631
13chlorobenzene0.02401.6200.9378
14bromobenzene0.04091.3880.9750
15para-chlorotoluene0.07191.1430.9202
16*ortho-chlorotoluene0.07521.1240.9481
171,4-dichlorobenzene0.01711.7670.9211
181,2-dichlorobenzene0.02751.5610.9673
191,2,4-trichlorobenzene0.03711.4310.9972
20*styrene0.08931.0490.9945

The degradation rate constant (k) and R2 of 20 kinds of monoaromatics.

*

Test set chemicals.

Construction of QSAR Model

The QSAR model was established by the training set data, as shown in Equation (10). The model revealed that the reaction rate constant had a significant correlation with . The R2 of the model was >0.80, and the SEE and high F values was low, indicating a good fit of the model. The and RMSECV values meant that the internal predictability of the model was satisfactory. The and RMSEp-values suggested that the established models had a good external predictive ability.

R2 = 0.825, SEE = 0.105, F = 61.215, = 0.766, RMSECV = 0.121, = 0.636, RMSEp = 0.133.

The is an electrostatic descriptor which reflects the non-uniformity of the charge distribution of the molecular, and it has been found to be correlated to the kinetic constant in several previous researches (Zhu et al., , ; Huang et al., ). A stronger negative net charge of carbon atoms indicated the ease of some valence-bond breakage of organic molecules, so would have a negative influence on degradation rate. A molecule with a smaller often gets a higher reaction rate.

The scatter plots of experimental vs. predicted values of –logk for the training and test sets are shown in Figure 3. In general, the values predicted by QSAR model were close to the experimental values. The relative residuals between the observed and predicted values were all below 20%. Figure 3 also showed that the chemicals were generally divided into two groups according to their degradation rates. Chemicals with alkyl groups (especially methyl groups) seemed to have a higher degradation rate, while those with halogen substitution usually suffered a low degradation rate.

Figure 3

Conclusion

This manuscript revealed the influencing factors of monoaromatic oxidation process by quinone-activated persulfate in aqueous solution. The pH, quinone species and concentration all affected the degradation efficiency. Monoaromatics were more easily oxidized in alkaline condition. The reasons were speculated as that base acted as the activator of persulfate to produce , and the reductive radicals accelerated the semiquinone/quinone cycle, facilitating the production of . Anionic form of phenols at alkaline pH was also the activating species in persulfate systems. However, excessive reductive radicals in strong alkaline environment inhibited the reaction. The degradation rate increased with quinone concentration, but excess quinone would compete for the oxidant and suppress the degradation of target contaminant. The initial degradation rates of monoaromatics fitted well with the pseudo-first-order equation, and a QSAR model relating the reaction rate with the chemical structure showed good correlation and predictability. The model revealed that the degradation efficiencies of monoaromatics were negatively correlated with . A smaller meant the non-uniform charge distribution of the molecular, resulting in easy breakage of the bond. This finding would contribute to further explore the effective activating method of persulfate in alkaline environment.

Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Statements

Data availability statement

The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.

Author contributions

JS designed the experiments and wrote the manuscript. CJ, DP, and YZ performed the experiments. LW analyzed the data. TL and XZ reviewed this paper and gave some advices. All authors contributed to the article and approved the submitted version.

Funding

The research was financially supported by the National Key Research and Development Project of China (2018YFC1801101), National Natural Science Foundation of China (Grants No. 21707041), and the Natural Science Foundation of Jiangsu Province (Grants No. BK20160101).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

Summary

Keywords

persulfate, quinone, degradation kinetics, quantitative structure-activity relationship, monoaromatics

Citation

Shi J, Long T, Zhou Y, Wang L, Jiang C, Pan D and Zhu X (2021) Efficiency and Quantitative Structure-Activity Relationship of Monoaromatics Oxidation by Quinone-Activated Persulfate. Front. Chem. 9:580643. doi: 10.3389/fchem.2021.580643

Received

06 July 2020

Accepted

03 March 2021

Published

01 September 2021

Volume

9 - 2021

Edited by

Guangshan Zhang, Qingdao Agricultural University, China

Reviewed by

Matteo Guidotti, National Research Council (CNR), Italy; Junhe Lu, Nanjing Agricultural University, China

Updates

Copyright

*Correspondence: Xin Zhu

This article was submitted to Green and Sustainable Chemistry, a section of the journal Frontiers in Chemistry

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics