==== Front Molecules Molecules molecules Molecules 1420-3049 MDPI 33287203 10.3390/molecules25235697 molecules-25-05697 Communication Quantum Chemical Investigation of Polychlorinated Dibenzodioxins, Dibenzofurans and Biphenyls: Relative Stability and Planarity Analysis https://orcid.org/0000-0001-8042-4513Datta Sopanant https://orcid.org/0000-0002-8558-6199Limpanuparb Taweetham * Gatti Carlo Academic Editor Science Division, Mahidol University International College, Mahidol University, Salaya, Phutthamonthon, Nakhon Pathom 73170, Thailand; sopanant.dat@student.mahidol.edu * Correspondence: taweetham.lim@mahidol.edu 03 12 2020 12 2020 25 23 569716 11 2020 01 12 2020 © 2020 by the authors.2020Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).All the possible polychlorinated aromatic compounds in the classes of dibenzodioxins (PCDDs), dibenzofurans (PCDFs), and biphenyls (PCBs) were studied by the quantum chemical methods of HF/6-311++G(d,p), B3LYP/6-311++G(d,p), and MP2/cc-pVTZ. The calculated stabilities and structures of these compounds were compared with the available data on their abundance and toxicity. Prediction models for trends in energy and planarity among these congeners were proposed. The results discussed here can help contribute to the understanding of the role of dioxin-like compounds (DLCs) in the environment. polychlorinated dibenzodioxins (PCDDs)polychlorinated dibenzofurans (PCDFs)polychlorinated biphenyls (PCBs) ==== Body 1. Introduction The Stockholm Convention on Persistent Organic Pollutants (POPs) [1], initially signed in 2001, includes 12 classes of compounds, three of which can be regarded as dioxin-like compounds (DLCs). (see Figure 1). A subset of DLCs was assigned toxic equivalency factors (TEFs) by the World Health Organization (WHO) [2] to indicate their potential health impacts on humans and mammals. These compounds may accumulate and bioconcentrate through the food chain and, as a result, food products in some jurisdictions such as the European Union are strictly monitored for the presence of these compounds [3]. Quantum chemical methods, semi-empirical calculations, and other related computational techniques have previously been applied to polychlorinated dibenzodioxins (PCDDs), polychlorinated dibenzofurans (PCDFs), and polychlorinated biphenyls (PCBs) [4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36]. The goals of these studies were, primarily, to predict their chemical properties and understand their mechanisms of action or degradation in biotic systems. A summary can be found in Supplementary Table S1. This communication describes the semi-automated codes for generating chemical structures, to be used as inputs for quantum chemical calculations in a standard package and to extract information from outputs for further analysis. Three different representative quantum chemical methods were selected for comparison: Hartree-Fock (HF), Density Functional Theory (DFT), and second-order Møller–Plesset perturbation theory (MP2). The lists of 76 PCDDs, 136 PCDFs, and 210 PCBs (inclusive of three non-chlorinated base structures) along with their systematic numberings [37,38,39] are shown in Supplementary Table S2–S4, respectively. The use of numberings herein will be in accordance with those in the lists. (This is not to be confused with PubChem ID in the Supplementary Materials and elution order number which is instrument-specific.) As shown in Supplementary Table S5, the congeners can be categorized into homologue groups by their empirical formulas or degrees of substitution. To the best of our knowledge, in comparison to the existing literature this is the most comprehensive theoretical study in terms of the number of compounds in the study and the level of theory employed for the investigation. The computational results are presented in terms of the relative stability and planarity of compounds. These calculated results are compared to the experimental results in the literature in terms of their isomer distribution and toxicity. A better understanding of the nature of these compounds can lead us to a better management of and more sustainable solutions for POPs in the future. 2. Results 2.1. Relative Stability The association between the substitution patterns of each compound and its relative stability, with respect to its isomers, was assessed. The substitution pattern parameters include the degree of chlorination (using the number of H atoms, or the number of Cl atoms subtracted from 10 for PCBs and 8 for PCDDs and PCDFs) and the intra-ring and cross-ring interactions between two substituents. Intra-ring interactions are quantified by the number of pairs of Cl and non-hydrogen substituents (Cl, O, and C) at o-, m- and p- positions on each phenyl ring. Cross-ring interactions are defined as the number of pairs of Cl substitution over an O-bridge (CR-O) and a C-C bond (CR-C). The possible parameters for intra-ring and cross-ring interactions are unique for each class of compounds, as shown in Table 1. The total number of o-, m- and p- pairs were also calculated from the tally of substituent pairs mentioned earlier. Overall, as listed in the table, there are a total of 11, 15, and 11 potential parameters for PCDDs, PCDFs, and PCBs, respectively. The following simple linear model was proposed for a relative stability prediction based on HF, B3LYP ∆G at 300 K, and MP2 electronic energy S = x0 + 102x1 + 104x2 + 106x3 + 108x4 + … where S is the energy score and x0, x1, x2, … are possible parameters. Coefficients of 102n are used as some variables can take up to two-digit values (0 to 12 for example). The association between energy scores and the calculated energy values was assessed using Spearman’s correlation coefficient (ρ). Inspired by the knapsack problem, parameters that give the highest ρ values were added to the model one by one, until the increase in the ρ value became less than 0.0001. Table 2 shows the resulting parameters, along with the corresponding ρ values. The results from all methodologies gave ρ > 0.99, and were mostly in agreement with each other. The major predictor for all compounds, in addition to the degree of substitution, are the intra-ring interactions. The predictor with the highest priority is the number of substituent pairs at o- positions, regardless of substituent type, followed by substitutions at the m- and p- positions, respectively. This supports the expectations of steric hindrance from the substituents (bulkiness: Cl > O ≈ C > H). Further, a cross-ring interaction over a C-C bond is an important parameter, especially for PCDFs, whereas a cross-ring interaction over an O-bridge is of less priority for PCDDs. With these results, the stability ranks of these compounds, regardless of methodology, can be sufficiently explained by their substituent pair interactions. For example, the energy scores S of PCDDs-13 and 14, when fitted according to the B3LYP results, can be calculated as: SPCDD-13 = x0 + 102x1 + 104x2 + 106x3 + 108x4 + 1010x5 = 1 + 102(0) + 104(1) + 106(2) + 108(3) + 1010(5) = 50,302,010,001 SPCDD-14 = x1 + 102x2 + 104x3+ 106x4+ 108x5 + 1010x6 = 2 + 102(1) + 104(1) + 106(1) + 108(3) + 1010(5) = 50,301,010,102, respectively. PCDD-13 ranks higher than PCDD-14 and is predicted to have a higher energy and less stability. The weights used in this rank modelling are to provide a non-overlapping priority of consideration for each predictor variable. Qualitative interpretation then becomes obvious. Multiple linear regression for the prediction of energy values, following a similar knapsack approach, was also explored and is shown in Supplementary Table S6. In this case, the priorities of the predictors overlap and the already high ρ values can be increased further. Our computational model is able to predict the relative stabilities of all of these classes of compounds with a stronger correlation than the linear models proposed earlier [16,17,18,19], while not increasing the numbers of parameters. 2.2. Isomer Distribution The energetic data in the previous subsection were used to predict the distribution of isomers within each homologue group. For this analysis, energy values were used to evaluate the distribution of isomers using the Boltzmann distribution at 300, 600, and 900 K [40]. The results were compared to available data on the abundance of these compounds in nature [37,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56]. A brief review of these data in the literature is shown in Supplementary Table S7. Our PCDDs and PCDFs results agree with a limited number of reports on experimental isomer distribution from incineration sources, while no clear trends were observed for PCBs (see Supplementary Table S8. The median ρ values calculated for each homologue group are 0.5833, 0.6555, and −0.1044 for PCDDs, PCDFs, and PCBs, respectively). For all three classes, the distribution of compounds found in nature usually differ due to the different accumulation and decomposition mechanisms in biotic and abiotic sources. For PCBs, it was argued that the distribution of isomers from the source is kinetically controlled rather than thermodynamically controlled [26]. 2.3. Planarity Planarity is widely claimed to be an important factor for these compounds to bind with the aryl hydrocarbon receptor (AhR) [57]. Therefore, geometric data were extracted and selected dihedral angles were calculated to represent the coplanarity of two phenyl rings. The designated carbon atoms for dihedral angle calculation are shown in Figure 1 and the list below. PCDD PCDF PCB Ca,Cc,Cd,Ch Cb,Cd,Cc,Cg Ca,Ce,Cf,Ch Cb,Cf,Ce,Cg Ca,Cd,Cf,Cg Cb,Cc,Ce,Ch Ci,Ck,Cl,Cp Cj,Cl,Ck,Co Ci,Cm,Cn,Cp Ci,Cn,Cm,Co Ci,Cl,Cn,Co Cj,Ck,Cm,Cp Cq,Cs,Ct,Cu Cq,Cs,Ct,Cv Cr,Cs,Ct,Cu Cr,Cs,Ct,Cv The lists for PCDDs and PCDFs are separated into two groups representing two orthogonal axes of rotation. Absolute values of the acute angle representations of the angles were used for the mean calculation of each group. Our data show that most PCDDs are planar, and therefore the data are insufficient for further prediction. For PCDFs, with a limited number of compounds exhibiting non-planarity, there appears to be a moderate trend (R values ranging from 0.61 to 0.65) between the substitution pattern and coplanarity. The significant parameters (p < 0.05) involved are the presence or absence of a substituent pair at positions 1 and 9 and at positions 2 and 8. The reasons for this finding may be unclear due to a limited number of data points (n = 13) for analysis and moderate correlation. The prediction models were derived from multiple linear regression and are shown in Table S9. For PCBs, a dihedral angle prediction, based on chlorine substitution position, was also assessed using multiple linear regression. In this model, the dihedral angle A of a compound is predicted by the equation: A = C0 + C1(x2 + x2′ + x6 + x6′) + C2(x2 + x6)(x2′ + x6′) + C3(x3 + x3′ + x5 + x5′) + C4(x4 + x4′) where; C1, C2, C3, and C4 are the coefficients to be determined; and xi is 1 if position i is occupied by a Cl atom, and 0 if otherwise. Table 3 shows the coefficients of significant parameters (p < 0.05) and the corresponding correlation coefficients for the predictions. The results from all methodologies show that the most important parameter is Cl substitution at any one of the four positions ortho to the C-C bond connecting two rings, as mentioned in earlier findings [57]. This is then followed by the same kind of substitution on both rings (product of x2 + x6 and x2′ + x6′). Meta substitutions according to all methods and para substitutions, according to MP2 results, are less important parameters, as C3 and C4 are relatively small. 2.4. Toxicity When comparing our geometric data to the toxic equivalency factor (TEF) values of these compounds [2], as shown in Supplementary Table S10, the following points are observed. The twelve known dioxin-like PCBs with TEF values are among the ones with relatively small dihedral angles or are closer to being planar. According to the B3LYP and MP2 results, there are three non-coplanar PCDDs (PCDDs-46, 65, and 71) from both methods, and two additional (PCDDs-42 and 45) from B3LYP only. All of these are not substituted at all four lateral positions (2,3,7,8), which are the substitution positions of PCDDs known to be toxic [2]. One of these, PCDD-46, is substituted at positions exactly opposite to that of 2,3,7,8-tetrachlorodibenzodioxin (TCDD), the most toxic dioxin-like compound [2]. 3. Discussion, Conclusions and Future Work Overall, this study provides high-quality quantum chemical calculations for the energetic and geometric properties of PCDDs, PCDFs, and PCBs. The comprehensive data provided for all of these three related classes of compounds can be used for further investigations. From our computational results, we have proposed prediction models for the important characteristics of these compounds. The relative stability prediction extends beyond earlier prediction models [16,17,18,19] with a high accuracy and with the allowance of qualitative interpretations. Prediction models have also been proposed for planarity, an important structural parameter for the toxic activities of these compounds. Furthermore, we have made associations between theoretical calculations (relative stability and planarity) and experimental data in the literature (abundance and toxicity). Our work can be used for further analyses on the persistence and activities of the toxic compounds in nature. In terms of methodology, semi-automated codes in the Supplementary Materials can be extended further for the construction and analysis of other related compounds. All three methods, HF, B3LYP, and MP2, yielded similar trends for the classes of compounds in our study. Our results suggest that a computationally expensive treatment for correlation energy may not be needed for this purpose and a low-level method such as HF, in most cases, can sufficiently provide insightful findings on the molecular properties of the studied classes of compounds. The further investigation of different variations of DFT methods may also be useful. For example, range-separated hybrid functionals such as ωB97XD can be used for comparison with the results presented here. 4. Materials and Methods We extended our exhaustive combinatorial investigation approach [58,59,60] to all members of the three classes of compounds. The Becke three-parameter hybrid functional combined with the Lee-Yang-Parr correlation functional (B3LYP), the most popular variation of DFT, was selected for this work. In total, three methods, HF/6-311++G(d,p), B3LYP/6-311++G(d,p), and MP2/cc-pVTZ [61], were used on Q-Chem 5.1 (developer version) [62] for geometry optimization. The HF- and B3LYP-optimised geometries were confirmed to be minima by verifying the absence of imaginary frequencies. All the raw output files and incomplete attempts made with MP2/6-311++G(d,p) are provided in the Supplementary Materials. The process of the structure enumeration and analysis of molecular properties, particularly relative energy and planarity, were completed by semi-automated scripts written in Mathematica for this project. The Mathematica notebook can be found in the Supplementary Materials. These Mathematica codes can be used for the construction and analysis of other compounds and for teaching purposes in classes relating to molecular modelling and cheminformatics. Our computational results were compared to the relative abundance of the compounds in nature and their estimated toxicity in humans and animals. The analyses presented are derived from the energetic and geometric data extracted from the computational results. For energetic information, the electronic energy of each compound was obtained from the results of all methodologies and, additionally, the enthalpy, entropy, and Gibbs energy values were obtained from the frequency jobs of HF and B3LYP. Due to the free rotation about C-C bonds connecting the two phenyl rings in PCBs, 78 potential atropisomers [38] were explored and only the structures with the lower energy values were used for further calculations. Acknowledgments We thank Wanutcha Lorpaiboon for constructive comments on the manuscript and are grateful to Wirote Chansri (UAE-IDEA Advance Analytical Company Limited) for an opportunity to learn about dioxin analysis. Sample Availability: Samples of the compounds are not available from the authors. Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. Supplementary Materials The following are available online, Table S1: Summary of computational data and analyses on PCDDs, PCDFs, and PCBs in the literature; Table S2: List of all PCDDs by substitution positions and numbering, proposed by Ballschmiter et al., revised by Dorofeeva et al.; Table S3: List of all PCDFs by substitution positions and numbering, proposed by Ballschmiter et al.; Table S4: List of all PCBs by substitution positions and numbering, first proposed by Ballschmiter et al., later revised by IUPAC; Table S5: List of congener groups of PCDDs, PCDFs and PCBs by degree of substitution; Table S6: Values for constant C0 and coefficients C1, C2, C3,…, parameters x1, x2, x3,… and correlation coefficients for the energy prediction model E = C0 + C1x1 + C2x2 + C3x3 + …, where E is the predicted HF, B3LYP ∆G at 300 K in Hartree and MP2 electronic energy in Hartree; Table S7: Summary of experimental data on abundance of PCDDs, PCDFs and PCBs in the literature; Table S8: Rank correlation coefficients (ρ values) for association between predicted distribution (using Boltzmann distribution) of PCDDs and PCDFs based on relative stability and observed distribution in fly ash samples from municipal waste incinerators (MWI); Table S9: Values for constant C0, coefficients C1, C2, and correlation coefficients for the PCDFs coplanarity prediction model; Table S10: List of toxic equivalency factor (TEF) values of 7 PCDDs, 10 PCDFs and 12 PCBs adapted from Van den Berg et al. Click here for additional data file. Author Contributions Conceptualization, S.D. and T.L.; methodology, S.D. and T.L.; software, T.L.; validation, S.D. and T.L.; formal analysis, S.D.; investigation, T.L.; resources, T.L.; data curation, S.D. and T.L.; writing—original draft preparation, S.D. and T.L.; writing—review and editing, S.D. and T.L.; visualization, S.D.; supervision, T.L.; project administration, T.L.; funding acquisition, T.L. All authors have read and agreed to the published version of the manuscript. Funding This research was funded by the National Research Council of Thailand (NRCT) and Mahidol University. Conflicts of Interest The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analysis, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results. Figure 1 Substitution position numberings (in red) and selected carbon atoms for dihedral angle calculation (in blue) for PCDDs, PCDFs, and PCBs. molecules-25-05697-t001_Table 1Table 1 Possible numbers of substituent pairs for each compound. Substitution Positions Substituent Pair PCDD PCDF PCB Intra-Ring Interactions o- Cl-Cl 0–6 1 0–6 0–8 Cl-O 0–4 0–2 N/A Cl-C N/A 0–2 0–4 Total 0–10 0–10 0–12 m- Cl-Cl 0–4 0–4 0–8 Cl-O 0–8 0–4 N/A Cl-C N/A 0–4 0–4 Total 0–12 0–12 0–12 p- Cl-Cl 0–2 0–2 0–4 Cl-O 0–4 0–2 N/A Cl-C N/A 0–2 0–2 Total 0–6 0–6 0–6 Cross-ring interactions Over O-bridge Cl-Cl 0–2 2 0–1 3 N/A Over C-C bond N/A 0–1 3 0–4 4 1 Example: A fully substituted PCDD (PCDD-75) has six o-Cl-Cl pairs at positions (1,2), (2,3), (3,4), (6,7), (7,8), and (8,9). Possible pairs are substitutions at positions: 2 (1,9) and (4,6); 3 (4,6) over O-bridge and (1,9) over C-C bond; 4 (2,2′), (2,6′), (2′,6), and (6,6′) due to the free rotation about the C-C bond connecting two rings. N/A indicates a case whereby such a substituent pair is not possible. molecules-25-05697-t002_Table 2Table 2 Parameters (in order of x0, x1, x2, …) for relative stability prediction and the corresponding correlation coefficients (ρ). PCDD HF B3LYP MP2 Parameters CR- O CR- O CR- O p- Cl-Cl p- Cl-Cl p- Cl-Cl m- Cl-Cl m- Cl-Cl m- Cl-Cl o- Cl-Cl o- Cl-Cl o- Cl-Cl o- total o- total o- total total H total H total H ρ 0.9998 0.9997 0.9996 PCDF HF B3LYP MP2 Parameters p- Cl-O m- Cl-Cl m- Cl-Cl m- Cl-Cl p- Cl-O p- Cl-Cl o- Cl-O o- Cl-O p- Cl-O o- Cl-Cl o- Cl-Cl o- Cl-O CR- C CR- C o- Cl-Cl total H o- total CR- C total H total H ρ 0.9997 0.9962 0.9978 PCB HF B3LYP MP2 Parameters p- total m- total p- total m- total o- Cl-Cl m- total o- Cl-Cl o- total o- Cl-Cl o- total total H total H total H ρ 0.9993 0.9990 0.9978 molecules-25-05697-t003_Table 3Table 3 Constants, coefficients, and Pearson’s (R) and Spearman’s (ρ) correlation coefficients for the PCB planarity prediction model. Coefficients R ρ C 0 C 1 C 2 C 3 C 4 HF 50.858 20.325 −10.751 0.724 N/A 0.9354 0.8735 B3LYP 37.926 24.996 −11.419 0.529 N/A 0.9552 0.9599 MP2 33.308 24.202 −9.499 0.668 −0.529 0.9800 0.9560 N/A indicates a case whereby the coefficient is not required and can be regarded as zero. ==== Refs References 1. United Nations Treaties and International Agreements Registered or Filed and Recorded with the Secretariat of the United Nations United Nations New York, NY, USA 2006 2. Van den Berg M. Birnbaum L.S. Denison M. De Vito M. Farland W. Feeley M. Fiedler H. Hakansson H. Hanberg A. Haws L. 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 10.1093/toxsci/kfl055 16829543 3. European Commission Commission Regulation (EC) No 1881/2006 of 19 December 2006 setting maximum levels for certain contaminants in foodstuffs Off. J. Eur. Union 2006 364 5 24 4. Arulmozhiraja S. Fujii T. Tokiwa H. Electron affinity for the most toxic 2,3,7,8-tetrachlorodibenzo-p -dioxin (TCDD):  A density functional theory study J. Phys. Chem. A 2000 104 7068 7072 10.1021/jp994237x 5. Dorofeeva O.V. Moiseeva N.F. Yungman V.S. Novikov V.P. Ideal gas thermodynamic properties of biphenyl Thermochim. Acta 2001 374 7 11 10.1016/S0040-6031(01)00492-0 6. Iino F. Tsuchiya K. Imagawa T. Gullett B.K. An isomer prediction model for PCNs, PCDD/Fs, and PCBs from municipal waste incinerators Environ. Sci. Technol. 2001 35 3175 3181 10.1021/es001857z 11505998 7. Arulmozhiraja S. Fujii T. Morita M. Density functional theory studies on radical ions of selected polychlorinated biphenyls J. Phys. Chem. A 2002 106 10590 10595 10.1021/jp021014p 8. Chana A. Concejero M.A. de Frutos M. González M.J. Herradón B. Computational studies on biphenyl derivatives. Analysis of the conformational mobility, molecular electrostatic potential, and dipole moment of chlorinated biphenyl:  Searching for the rationalization of the selective toxicity of polychlorinated biphenyls (PCBs) Chem. Res. Toxicol. 2002 15 1514 1526 10.1021/tx025596d 12482233 9. León L.A. Notario R. Quijano J. Sánchez C. Structures and enthalpies of formation in the gas phase of the most toxic polychlorinated dibenzo-p -dioxins. A DFT study J. Phys. Chem. A 2002 106 6618 6627 10.1021/jp025602e 10. Dorofeeva O.V. Yungman V.S. Enthalpies of formation of dibenzo-p -dioxin and polychlorinated dibenzo-p -dioxins calculated by density functional theory J. Phys. Chem. A 2003 107 2848 2854 10.1021/jp026973y 11. Lee J.-E. Choi W. DFT calculation on the electron affinity of polychlorinated dibenzo-p -dioxins Bull. Korean Chem. Soc. 2003 24 792 796 10.5012/BKCS.2003.24.6.792 12. Lee J.E. Choi W. Mhin B.J. DFT calculation on the thermodynamic properties of polychlorinated dibenzo-p -dioxins:  intramolecular Cl−Cl repulsion effects and their thermochemical implications J. Phys. Chem. A 2003 107 2693 2699 10.1021/jp027133m 13. Arulmozhiraja S. Morita M. Electron affinities and reductive dechlorination of toxic polychlorinated dibenzofurans: A density functional theory study J. Phys. Chem. A 2004 108 3499 3508 10.1021/jp037620h 14. Dorofeeva O.V. Moiseeva N.F. Yungman V.S. Thermodynamic properties of polychlorinated biphenyls in the gas phase J. Phys. Chem. A 2004 108 8324 8332 10.1021/jp0477710 15. 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 10.1021/es0497227 15506206 16. Wang Z.-Y. Zhai Z.-C. Wang L.-S. Chen J.-L. Kikuchi O. Watanabe T. Prediction of gas phase thermodynamic function of polychlorinated dibenzo-p -dioxins using DFT J. Mol. Struc. Theochem 2004 672 97 104 10.1016/j.theochem.2003.11.030 17. Wang Z.-Y. Zhai Z.-C. Wang L.-S. Prediction of gas phase thermodynamic properties of polychlorinated dibenzo-furans by DFT J. Mol. Struc. Theochem 2005 725 55 62 10.1016/j.theochem.2005.01.037 18. Wang Z.Y. Han X.Y. Zhai Z.C. Wang L.S. Study on the thermodynamic property and relative stability of a series of polychlorinated biphenyls by density functional theory Acta Chim. Sinica 2005 63 964 972 19. Zhai Z.-C. Wang Z.-Y. Chen X.-H. Wang L.-S. DFT calculation on 204 polychlorinated biphenyls: Their thermodynamic function and implication of Cl substitute position J. Mol. Struc. Theochem 2005 714 123 131 10.1016/j.theochem.2004.08.060 20. Dorofeeva O.V. Vishnevskiy Y.V. Moiseeva N.F. Assessment of Gaussian-3X theory for chlorinated organic molecules. Enthalpies of formation of chlorobenzenes and predictions for polychlorinated aromatic compounds Struct. Chem. 2006 17 383 392 10.1007/s11224-006-9047-8 21. Thompson D. Ewan B.C.R. A group additivity algorithm for polychlorinated dibenzofurans derived from selected DFT analyses J. Phys. Chem. A 2007 111 5043 5047 10.1021/jp0709499 17518453 22. Wang L. Heard D.E. Pilling M.J. Seakins P. A Gaussian-3X prediction on the enthalpies of formation of chlorinated phenols and dibenzo-p -dioxins J. Phys. Chem. A 2008 112 1832 1840 10.1021/jp077478a 18232671 23. Jansson S. Antti H. Marklund S. Tysklind M. Multivariate relationships between molecular descriptors and isomer distribution patterns of PCDD/Fs formed during MSW combustion Environ. Sci. Technol. 2009 43 7032 7038 10.1021/es901202w 19806738 24. Wang L. He Y.-L. The enthalpies of formation for polychlorinated dibenzofurans with use of G3XMP2 model chemistry and density functional theory J. Phys. Chem. A 2009 113 238 245 10.1021/jp802091z 19072627 25. Ewan B.C.R. Thompson D. Predicted specific heat behaviour of the polychlorinated dibenzofuran family from DFT analysis Thermochim. Acta 2013 558 46 52 10.1016/j.tca.2013.02.012 26. Jansson S. Grabic R. Multivariate relationships between molecular descriptors and isomer distribution patterns of PCBs formed during household waste incineration Environ. Sci. Pollut. Res. 2014 21 3082 3090 10.1007/s11356-013-2257-x 27. Zhou Q. Su X. Yong Y. Ju W. Fu Z. Li X. Adsorption of 2, 3, 7, 8-tetrachlorodibenzao-p -dioxin (TCDD) on graphane decorated with Ni and Cu: A DFT study Vacuum 2018 149 53 59 10.1016/j.vacuum.2017.12.016 28. Abronin I.A. Volkova L.V. On the relationship between the energy characteristics of the isodesmic reactions of polychlorinated dioxins and their toxicity Russ. Chem. Bull. 2019 68 867 869 10.1007/s11172-019-2498-z 29. Bai N. Wang W. Zhao Y. Feng W. Li P. Theoretical insights into the reaction mechanism between 2,3,7,8-tetrachlorodibenzofuran and hydrogen peroxide: A DFT study ACS Omega 2019 4 358 367 10.1021/acsomega.8b00724 31459335 30. Hou S. Altarawneh M. Kennedy E.M. Mackie J.C. Weber R. Dlugogorski B.Z. Formation of polychlorinated dibenzo-p -dioxins and dibenzofurans (PCDD/F) from oxidation of 4,4′-dichlorobiphenyl (4,4′-DCB) Proc. Combust. Inst. 2019 37 1075 1082 10.1016/j.proci.2018.05.045 31. Behjatmanesh-Ardakani R. Heydari A. Molecular and dissociative adsorption of tetrachlorodibenzodioxin on M-doped graphenes (M = B, Al, N, P): Pure DFT and DFT + VdW calculations J. Mol. Model. 2020 26 164 10.1007/s00894-020-04387-4 32494890 32. Diao J. Li Y. Shi S. Sun Y. Sun Y. QSAR models for predicting toxicity of polychlorinated dibenzo-p -dioxins and dibenzofurans using quantum chemical descriptors Bull. Environ. Contam. Toxicol. 2010 85 109 115 10.1007/s00128-010-0065-2 20628729 33. Larsson M. Kumar Mishra B. Tysklind M. Linusson A. Andersson P.L. On the use of electronic descriptors for QSAR modelling of PCDDs, PCDFs and dioxin-like PCBs SAR QSAR Environ. Res. 2013 24 461 479 10.1080/1062936X.2013.791719 23724952 34. Hirano M. Hwang J.-H. Park H.-J. Bak S.-M. Iwata H. Kim E.-Y. In silico analysis of the interaction of avian aryl hydrocarbon receptors and dioxins to decipher isoform-, ligand-, and species-specific activations Environ. Sci. Technol. 2015 49 3795 3804 10.1021/es505733f 25692546 35. Larsson M. Fraccalvieri D. Andersson C.D. Bonati L. Linusson A. Andersson P.L. Identification of potential aryl hydrocarbon receptor ligands by virtual screening of industrial chemicals Environ. Sci. Pollut. Res. 2018 25 2436 2449 10.1007/s11356-017-0437-9 36. Zeinali N. Oluwoye I. Altarawneh M. Dlugogorski B.Z. Destruction of dioxin and furan pollutants via electrophilic attack of singlet oxygen Ecotoxicol. Environ. Saf. 2019 184 109605 10.1016/j.ecoenv.2019.109605 31505406 37. Ballschmiter K. Buchert H. Niemczyk R. Munder A. Swerev M. Automobile exhausts versus municipal-waste incineration as sources of the polychloro-dibenzodioxins (PCDD) and -furans (PCDF) found in the environment Chemosphere 1986 15 901 915 10.1016/0045-6535(86)90055-X 38. Agudo A. Aronson K.J. Bonefeld-Jorgensen E.C. Cocco P. Cogliano V. Cravedi J.-P. Esch H. Fiedler H. Glauert H.P. Guo Y.-L.L. Polychlorinated Biphenyls and Polybrominated Biphenyls International Agency for Research on Cancer Lyon, France 2016 39. Ballschmiter K. Zell M. Analysis of polychlorinated biphenyls (PCB) by glass capillary gas chromatography Fresenius Z. Anal. Chem. 1980 302 20 31 10.1007/BF00469758 40. Dickson L.C. Lenoir D. Hutzinger O. Quantitative comparison of de novo and precursor formation of polychlorinated dibenzo-p -dioxins under simulated municipal solid waste incinerator postcombustion conditions Environ. Sci. Technol. 1992 26 1822 1828 10.1021/es00033a017 41. Tiernan T.O. Taylor M.L. Garrett J.H. VanNess G.F. Solch J.G. Deis D.A. Wagel D.J. Chlorodibenzodioxins, chlorodibenzofurans and related compounds in the effluents from combustion processes Chemosphere 1983 12 595 606 10.1016/0045-6535(83)90217-5 42. Schwartz T.R. Campbell R.D. Stalling D.L. Little R.L. Petty J.D. Hogan J.W. Kaiser E.M. Laboratory data base for isomer-specific determination of polychlorinated biphenyls Anal. Chem. 1984 56 1303 1308 10.1021/ac00272a025 43. Ballschmiter K. Swerev M. Reaction pathways for the formation of polychlorodibenzodioxins (PCDD) and -furans (PCDF) in combustion processes I Fresenius Z. Anal. Chem. 1987 328 125 127 10.1007/BF00560966 44. Yasuhara A. Ito H. Morita M. Isomer-Specific Determination of polychlorinated dibenzo-p -dioxins and dibenzofurans in incinerator-related environmental samples Environ. Sci. Technol. 1987 21 971 979 10.1021/es50001a010 19994995 45. Faroon O.M. Samuel Keith L. Smith-Simon C. De Rosa C.T. Polychlorinated Biphenyls: Human Health Aspects World Health Organization Geneva, Switzerland 2003 46. Hanari N. Horii Y. Okazawa T. Falandysz J. Bochentin I. Orlikowska A. Puzyn T. Wyrzykowska B. Yamashita N. Dioxin-like compounds in pine needles around Tokyo Bay, Japan in 1999 J. Environ. Monit. 2004 6 305 312 10.1039/b311176h 15054539 47. Kim K.S. Hirai Y. Kato M. Urano K. Masunaga S. Detailed PCB congener patterns in incinerator flue gas and commercial PCB formulations (Kanechlor) Chemosphere 2004 55 539 553 10.1016/j.chemosphere.2003.11.056 15006507 48. Haws L.C. Su S.H. Harris M. DeVito M.J. Walker N.J. Farland W.H. Finley B. Birnbaum L.S. Development of a refined database of mammalian relative potency estimates for dioxin-like compounds Toxicol. Sci. 2006 89 4 30 10.1093/toxsci/kfi294 16120753 49. Ishikawa Y. Noma Y. Yamamoto T. Mori Y. Sakai S.-i. PCB decomposition and formation in thermal treatment plant equipment Chemosphere 2007 67 1383 1393 10.1016/j.chemosphere.2006.10.022 17134732 50. Srogi K. Levels and congener distributions of PCDDs, PCDFs and dioxin-like PCBs in environmental and human samples: A review Environ. Chem. Lett. 2008 6 1 28 10.1007/s10311-007-0105-2 51. Sundqvist K.L. Tysklind M. Geladi P. Hopke P.K. Wiberg K. PCDD/F source apportionment in the Baltic sea using positive matrix factorization Environ. Sci. Technol. 2010 44 1690 1697 10.1021/es9030084 20121084 52. Jansson S. Lundin L. Grabic R. Characterisation and fingerprinting of PCBs in flue gas and ash from waste incineration and in technical mixtures Chemosphere 2011 85 509 515 10.1016/j.chemosphere.2011.08.012 21885088 53. Barahona-Urbina C. Nunez-Gonzalez S. Gomez-Jeria J.S. Model-based quantum-chemical study of the uptake of some polychlorinated pollutant compounds by zucchini subspecies J. Chil. Chem. Soc. 2012 57 1497 1503 10.4067/S0717-97072012000400032 54. Hinwood A.L. Callan A.C. Heyworth J. Rogic D. de Araujo J. Crough R. Mamahit G. Piro N. Yates A. Stevenson G. Polychlorinated biphenyl (PCB) and dioxin concentrations in residential dust of pregnant women Environ. Sci. Process. Impacts 2014 16 2758 2763 10.1039/C4EM00383G 25343442 55. Thompson D. Ooi T.C. Anderson D.R. Fisher R. Ewan B.C.R. The polychlorinated dibenzofuran fingerprint of iron ore sinter plant: Its persistence with suppressant and alternative fuel addition Chemosphere 2016 154 138 147 10.1016/j.chemosphere.2016.03.066 27043380 56. Huang R. Wang P. Zhang J. Chen S. Zhu P. Huo W. Jiang Y. Chen Z. Peng J. The human body burden of polychlorinated dibenzo-p -dioxins/furans (PCDD/Fs) and dioxin-like polychlorinated biphenyls (DL-PCBs) in residents’ human milk from Guangdong Province, China Toxicol. Res. 2019 8 552 559 10.1039/C8TX00337H 57. The Eightieth Meeting of the Joint FAO/WHO Expert Committee on Food Additives (JECFA) Safety Evaluation of Certain Food Additives and Contaminants, Supplement 1: Non-Dioxin-Like Polychlorinated Biphenyls World Health Organization Geneva, Switzerland 2016 58. Datta S. Limpanuparb T. Geometric and energetic data from quantum chemical calculations of halobenzenes and xylenes Data Brief 2020 30 105386 10.1016/j.dib.2020.105386 32258270 59. Limpanuparb T. Datta S. Chinsukserm K. Teeraniramitr P. In silico geometric and energetic data of all possible simple rotamers made of non-metal elements Data Brief 2020 30 105442 10.1016/j.dib.2020.105442 32292807 60. Chinsukserm K. Lorpaiboon W. Teeraniramitr P. Limpanuparb T. Geometric and energetic data from ab initio calculations of haloethene, haloimine, halomethylenephosphine, haloiminophosphine, halodiazene, halodiphosphene and halocyclopropane Data Brief 2019 27 104738 10.1016/j.dib.2019.104738 31763395 61. Riley K.E. Hobza P. Assessment of the MP2 method, along with several basis sets, for the computation of interaction energies of biologically relevant hydrogen bonded and dispersion bound complexes J. Phys. Chem. A 2007 111 8257 8263 10.1021/jp073358r 17649987 62. Shao Y. Gan Z. Epifanovsky E. Gilbert A.T. Wormit M. Kussmann J. Lange A.W. Behn A. Deng J. Feng X. Advances in molecular quantum chemistry contained in the Q-Chem 4 program package Mol. Phys. 2015 113 184 215 10.1080/00268976.2014.952696