研究论文

一种量子化学组合方法及其应用于汞对硫酸气溶胶形成影响的研究

  • 李程桥 ,
  • 王一波
展开
  • 贵州大学化学与化工学院 贵州省高性能计算化学重点实验室 贵阳 550025

收稿日期: 2021-04-12

  网络出版日期: 2021-06-01

基金资助

国家自然科学基金(21767005); 贵州省科技计划项目(20175788)

A Combination Method of Quantum Chemistry and Its Application to the Study of the Effects of Mercury on the Formation of Sulfuric Acid Aerosol

  • Chengqiao Li ,
  • Yibo Wang
Expand
  • School of Chemistry and Chemical Engineering, Guizhou University, Key Laboratory of High Performance Computational Chemistry, Guiyang 550025, China

Received date: 2021-04-12

  Online published: 2021-06-01

Supported by

National Natural Science Foundation of China(21767005); Science and Technology Plan Project of Guizhou Province(20175788)

摘要

本工作建立了一种量子化学的组合方法并应用于研究在硫和汞复合污染地区, 汞及其化合物是否会参与和驱动硫酸气溶胶形成. 该方法采用Molclus Genmer模块生成数以百万计的巨量团簇初始结构, 用自洽紧束缚GFN2-xTB方法预优化, 再用B3PW91-D3(BJ) DFT方法重新优化, 结合能经ωB97M-D3(BJ)校正后计算出团簇的热力学函数, 以其生成过程Gibbs自由能变大小排序, 获得团簇稳定序列及其对应的结合能和热力学函数. 然后用该方法研究了汞及其化合物、硫酸和水生成的8种团簇模型共计240万个结构, 结果表明HgSO4, HgO等汞的极性化合物以及Hg22+和 Hg2+能够促进硫酸气溶胶的形成, 而Hg以及无极性的HgCl2和Hg2Cl2对硫酸气溶胶成核的影响较小. 这一结论对硫、汞复合污染环境下, 硫酸气溶胶成核机制探讨提供了有价值的预测.

本文引用格式

李程桥 , 王一波 . 一种量子化学组合方法及其应用于汞对硫酸气溶胶形成影响的研究[J]. 化学学报, 2021 , 79(8) : 1065 -1072 . DOI: 10.6023/A21040147

Abstract

Haze weather has been occurring frequently in recent years, and studies have shown that its formation mechanism is related to the nucleation of sulfuric acid aerosols. How mercury and its compounds affect the formation and growth of sulfuric acid aerosols is completely unclear in polluted areas where sulfur and mercury coexist. In order to explore this problem theoretically, we established a new quantum chemistry combination method. Millions of cluster structures are generated by using the Molclus Genmer module. The cluster geometries are preliminarily optimized by the GFN2-xTB method and sorted according to the Gibbs free energy changes of the cluster generation process. The top 100 optimized structures were re-optimized using the B3PW91-D3(BJ) method. In order to accurately calculate the binding energy of clusters, we established a small benchmark set based on the composition and structural characteristics of the clusters studied in the project, and used the CCSD(T)/CBS method to calculate the binding energy of each cluster model in the benchmark set. Using the binding energy as a standard, the best method selected from 42 DFT methods is ωB97M-D3(BJ). The binding energy calculated by B3PW91-D3(BJ) is corrected by the results of the ωB97M-D3(BJ) method to obtain an accurate thermodynamic function value. It is worth mentioning that the combined method was used to study a total of 2.4 million structures in 8 cluster models formed by mercury and its typical compounds, sulfuric acid and water. The results show that polar compounds of mercury such as HgSO4 and HgO, together with Hg22+ and Hg2+, can promote sulfuric acid aerosol formation, while Hg and non-polar HgCl2 and Hg2Cl2 have almost no effect on the nucleation of sulfuric acid aerosol. In addition, we found that the conclusion holds true for the temperature interval of the troposphere at standard atmospheric pressure. This conclusion provides valuable predictions for the studies on sulfuric acid aerosol nucleation mechanisms under the combined sulfur and mercury pollution environment.

参考文献

[1]
Sipilä, M.; Berndt, T.; Petäjä, T.; Brus, D.; Vanhanen, J.; Stratmann, F.; Patokoski, J.; Mauldin, R. L.; Hyvärinen, A. P.; Lihavainen, H.; Kulmala, M. Science 2010, 327, 1243.
[2]
Asaduzzaman, A.; Riccardi, D.; Afaneh, A. T.; Cooper, S. J.; Smith, J. C.; Wang, F.; Parks, J. M.; Schreckenbach, G. Acc. Chem. Res. 2019, 52, 379.
[3]
Hong, Q.; Xie, Z.; Liu, C.; Wang, F.; Xie, P.; Kang, H.; Xu, J.; Wang, J.; Wu, F.; He, P.; Mou, F.; Fan, S.; Dong, Y.; Zhan, H.; Yu, X.; Chi, X.; Liu, J. Atmos. Chem. Phys. 2016, 16, 13807.
[4]
Li, L.; Zhang, Y. Y.; Jiao, C. Y.; Yao, Y. W.; Zhang, H.; Tian, Y. M. Environ. Monit. China 2019, 35, 40 (in Chinese.)
[4]
李亮, 张艳艳, 焦聪颖, 姚雅伟, 张辉, 田英明, 中国环境监测, 2019, 35, 40.)
[5]
Doyle, G. J. J. Chem. Phys. 1961, 35, 795.
[6]
Jaecker, V. A.; Mirabel, P. Atmos. Environ. 1989, 23, 2053.
[7]
Sucarrat, M. T.; Francisco, J. S.; Anglada, J. M. J. Am. Chem. Soc. 2012, 134, 20632.
[8]
Kildgaard, J. V.; Mikkelsen, K. V.; Bilde, M.; Elm, J. J. Phys. Chem. A 2018, 122, 5026.
[9]
Humphries, R. S.; Schofield, R.; Keywood, M. D.; Ward, J.; Pierce, J. R.; Gionfriddo, C. M.; Tate, M. T.; Krabbenhoft, D. P.; Galbally, I. E.; Molloy, S. B.; Klekociuk, A. R.; Johnston, P. V.; Kreher, K.; Thomas, A. J.; Robinson, A. D.; Harris, N. R.P.; Johnson, R.; Wilson, S. R. Atmos. Chem. Phys. 2015, 15,13339.
[10]
Lu, T. Molclus program, Version 1.8.7, Beijing Kein Research Center for Natural Science, China, 2018, http://www.keinsci.com/research/molclus.html
[11]
Bannwarth, C.; Ehlert, S.; Grimme, S. J. Chem. Theory Comput. 2019, 15, 1652.
[12]
Tsuzuki, S.; Honda, K.; Uchimaru, T.; Mikami, M.; Tanabe, K. J. Am. Chem. Soc. 2002, 124, 104.
[13]
Hu, M.; Shang, D. J.; Guo, S.; Wu, Z. J. Acta Chim. Sinica 2016, 74, 385 (in Chinese.)
[13]
胡敏, 尚冬杰, 郭松, 吴志军, 化学学报, 2016, 74, 385.)
[14]
Elm, J.; Passananti, M.; Kurtén, T.; Vehkamäki, H. J. Phys. Chem. A 2017, 121, 6155.
[15]
Yang, P.; Ye, Z. L.; Jiang, G. Y.; Li, Z.; Ding, C. F.; Hou, H. Q. Acta Chim. Sinica 2009, 67, 2031 (in Chinese.)
[15]
杨鹏, 叶招莲, 蒋公羽, 李周, 丁传凡, 侯惠奇, 化学学报, 2009, 67, 2031.)
[16]
Lu, Q.; Luo, Q.; Huang, S.; Li, Y.; Wan, J. J. Chem. Phys. A 2016, 120, 4560.
[17]
Lu, Q.; Luo, Q.; Huang, S.; Li, Y. Phys. Chem. Chem. Phys. 2017, 19, 28434.
[18]
Rasmussen, F. R.; Besel, V.; Mikkelsen, K. V.; Bilde, M.; Elm, J. J. Phys. Chem. A 2020, 124, 5253.
[19]
Dohm, S.; Bursch, M.; Hansen, A.; Grimme, S. J. Chem. Theory Comput. 2020, 16, 2002.
[20]
Myllys, N.; Olenius, T.; Kurtén, T.; Vehkamäki, H.; Riipinen, I.; Elm, J. J. Phys. Chem. A 2017, 121, 4812.
[21]
Becke, A. D. J. Chem. Phys. 1993, 98, 5648.
[22]
Perdew, J. P.; Wang, Y. Phys. Rev. 1992, 45, 13244.
[23]
Amaro-Estrada, J. I.; Maron, L.; Ramirez-Solis, A. Phys. Chem. Chem. Phys. 2014, 16, 8455.
[24]
Castro, L.; Dommergue, A.; Renard, A.; Ferrari, C.; Ramirez-Solis, A.; Maron, L. Phys. Chem. Chem. Phys. 2011, 13, 16772.
[25]
Yoo, S.; Apra, E.; Zeng, X. C.; Xantheas, S. S. J. Phys. Chem. Lett. 2010, 1, 3122.
[26]
Neese, F.; Hansen, A.; Liakos, D. G. J. Chem. Phys. 2009, 131,064103.
[27]
He, Y.; Wang, Y. B. Acta Phys.-Chim. Sin. 2017, 33, 1149 (in Chinese.)
[27]
何禹, 王一波, 物理化学学报, 2017, 33, 1149).
[28]
Najibi, A.; Goerigk, L. J. Chem. Theory Comput. 2018, 14, 5725.
[29]
Mardirossian, N.; Pestana, L. R.; Womack, J. C.; Skylaris, C. K.; Head-Gordon, T.; Head-Gordon, M. J. Phys. Chem. Lett. 2017, 8, 35.
[30]
Mardirossian, N.; Gordon, M. H. J. Chem. Phys. 2016, 144,214110.
[31]
Boys, S. F.; Bernardi, F. Mol. Phys. 1970, 19, 553.
[32]
Kruse, H.; Grimme, S. J. Chem. Phys. 2012, 136,154101.
[33]
Grimme, S. Growing String Method, xTB Documentation, 2019, https://xtb-docs.readthedocs.io/en/latest/gsm.html
[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.; Nakatsuji, H.; Caricato, M.; Li, X.; Hratchian, H. P.; Izmaylov, A. F.; Bloino, J.; Zheng, G.; Sonnenberg, J. L.; Hada, M.; Ehara, M.; Toyota, K.; Fukuda, R.; Hasegawa, J.; Ishida, M.; Nakajima, T.; Honda, Y.; Kitao, O.; Nakai, H.; Vreven, T.; Montgomery Jr., J. A.; Peralta, J. E.; Ogliaro, F.; Bearpark, M.; Heyd, J. J.; Brothers, E.; Kudin, K. N.; Staroverov, V. N.; Keith, T.; Kobayashi, R.; Normand, J.; Raghavachari, K.; Rendell, A.; Burant, J. C.; Iyengar, S. S.; Tomasi, J.; Cossi, M.; Rega, N.; Millam, J. M.; Klene, M.; Knox, J. E.; Cross, J. B.; Bakken, V.; Adamo, C.; Jaramillo, J.; Gomperts, R.; Stratmann, R. E.; Yazyev, O.; Austin, A. J.; Cammi, R.; Pomelli, C.; Ochterski, J. W.; Martin, R. L.; Morokuma, K.; Zakrzewski, V. G.; Voth, G. A.; Salvador, P.; Dannenberg, J. J.; Dapprich, S.; Daniels, A. D.; Farkas, O.; Foresman, J. B.; Ortiz, J. V.; Cioslowski, J.; Fox, D. J. Gaussian 16, Revision B. 01, Gaussian, Inc., Wallingford CT, 2016.
[35]
Neese, F., Software update: the ORCA program system, version 4.2 WIREs: Comput. Mol. Sci. 2018, 8,e1327.
[36]
Werner, H. J.; Knowles, P. J.; Knizia, G.; Manby, F. R.; Schütz, M.; Celani, P.; Györffy, W.; Kats, D.; Korona, T.; Lindh, R.; Mitrushenkov, A.; Rauhut, G.; Shamasundar, K. R.; Adler, T. B.; Amos, R. D.; Bennie, S. J.; Bernhardsson, A.; Berning, A.; Cooper, D. L.; Deegan, M. J.O.; Dobbyn, A. J.; Eckert, F.; Goll, E.; Hampel, C.; Hesselmann, A.; Hetzer, G.; Hrenar, T.; Jansen, G.; Köppl, C.; Lee, S. J.R.; Liu, Y.; Lloyd, A. W.; Ma, Q.; Mata, R. A.; May, A. J.; McNicholas, S. J.; Meyer, W.; Miller III, T. F.; Mura, M. E.; Nicklass, A.; O’Neill, D. P.; Palmieri, P.; Peng, D.; Pflüger, K.; Pitzer, R.; Reiher, M.; Shiozaki, T.; Stoll, H.; Stone, A. J.; Tarroni, R.; Thorsteinsson, T.; Wang, M.; Welborn, M. MOLPRO, version 2018, a package of ab initio programs, http://www.molpro.net
[37]
Yang, Z. Z.; Meng, X. F.; Zhao, D. X.; Gong, L. D. Acta Chim. Sinica 2009, 67, 2074 (in Chinese.)
[37]
杨忠志, 孟祥凤, 赵东霞, 宫利东 化学学报 2009, 67, 2074).
[38]
Temelso, B.; Morrell, T. E.; Shields, R. M.; Allodi, M. A.; Wood, E. K.; Kirschner, K. N.; Castonguay, T. C.; Archer, K. A.; Shields, G. C. J. Phys. Chem. A 2012, 116, 2209.
文章导航

/