hjcet Hans Journal of Chemical Engineering and Technology 2161-8844 2161-8852 beplay体育官网网页版等您来挑战！ 10.12677/hjcet.2024.146052 hjcet-101090 Articles 工程技术 铁盐除砷中水化作用的微观研究
Microscopic Study on Hydration in Arsenic Removal by Iron Salt
蒋欣欣 李会吉 刘一涵 张东晓 刘真真 马世博 孙海杰 郑州师范学院化学化工学院，河南 郑州 04 11 2024 14 06 482 490 23 9 ：2024 15 9 ：2024 15 11 ：2024 Copyright © 2024 beplay安卓登录 All rights reserved. 2024 This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/ 水溶液中砷的除去与水化作用密切相关。本文借助量子化学方法分析了铁盐除砷时砷酸盐周围水化特征的情况。对于[H mAsO 4(H 2O) n] 3− m (m = 0~2)水合团簇，随着氢质子化，[H mAsO 4(H 2O) n] 3− m的水化减弱，O N-H W距离小于O P-H W，而O N-H W配位数大于O P-H W。Fe 3+与H mAsO 4 m −3缔合改变了H mAsO 4 m −3水化情况。相对于[H mAsO 4] 3− m (m = 0~2)物种，[FeH mAsO 4] m + (m = 0~2)缔合团簇中[H mAsO 4] 3− m的水化能力降低，从而达到沉淀除砷。[FeH mAsO 4] m + (m = 0~2)缔合团簇的水化则随着质子化，[FeH mAsO 4] m + (m = 0~2)的水化增强。
The removal of arsenic from aqueous solution is closely related to hydration. In this paper, the hydration characteristics around arsenate during arsenic removal from iron salt were analyzed by means of quantum chemistry. For the hydrating cluster of [H mAsO 4(H 2O) n] 3− m (m = 0~2), the hydration of [H mAsO 4(H 2O) n] 3− m decreases with hydrogen protonation, and the O N-H W distance is less than the O P-H W, while the O N-H W coordination number is greater than the O P-H W. The association of Fe 3+ with H mAsO 4 m −3 changed the hydration of H mAsO 4 m −3. Compared with [H mAsO 4] 3− m (m = 0~2) species, the hydration capacity of [H mAsO 4] 3− m in [FeH mAsO 4] m + (m = 0~2) association cluster was reduced, so as to achieve precipitation and arsenic removal. The hydration of [FeH mAsO 4] m + (m = 0~2) association clusters was enhanced with the protonation of [FeH mAsO 4] m + (m = 0~2).
铁盐除砷，量子化学，缔合结构，水化作用
Arsenic Removal by Iron Salt
Quantum Chemistry Association Structure Hydration 1. 引言

砷是一种非人体必需的元素，在自然界中以四种氧化形式广泛存在。砷具有很强的毒性和生物利用度，被世界卫生组织(WHO)和国际癌症研究机构列为第一类致癌物。砷的服用量达到0.01~0.052 g就会引起中毒，达到0.06~0.2 g甚至会导致死亡。因此，各国都将出台相关法律，严格限制土壤和水体中砷的浓度。各国地表水砷的最大允许浓度一般为50 μg/L [1] [2] 。研究发现，砷可以作用于神经系统，刺激造血器官，入侵人体后影响红细胞的产生。长期接触高砷污染的水会损害皮肤、肺和消化系统，甚至可能诱发恶性肿瘤 [3] 。由于砷的自然释放和采矿、冶金、煤矿、化工等人类活动的影响，砷污染越来越严重。目前，砷污染在世界范围内受到越来越多的关注，砷的去除已经成为一个亟待解决的问题。

国内外铁盐除砷技术，主要有化学沉淀法、电化学法、吸附法、离子交换法、膜分离法、生物法等 [4] [5] 。化学沉淀可细分为两种类型：一种是中间产物沉淀砷，然后转化为砷产品销售 [6] ，另一种是稳定的化合物沉淀砷。许多常规化学药剂(石灰、氯化铁、聚合硫酸铁、硫酸亚铁、明矾、硫化钠等)在合适的pH值下，可与As(V)反应生成砷沉淀物，其中有些药剂也可与As(Ⅲ)反应形成沉淀物 [7] 。但是，由于As(Ⅲ)的毒性远高于As(V)，铁盐，尤其是三价铁与五价砷形成的砷酸铁(FeAsO4)的溶度积非常小，其除砷效率也较高 [8] 。因此，目前世界上应用最广泛的从废水中沉淀砷的方法主要是基于三价铁与五价砷的反应。比如研究表明，用三氯化铁除砷去除废水中的As(Ⅲ)和As(V)，对于As(Ⅲ)废水，最佳工艺条件为：pH = 9，nFe/As = 2，反应时间1 h；对于As(V)废水，最佳工艺条件为：pH = 8，nFe/As = 2，反应时间1 h。在上述条件下，砷的去除率都可达99.2%以上 [9] 。可见，研究铁盐除砷中水化作用的微观特征具有重要的意义。Bhattacharjee曾对 HAsO 4 2 − 水合作用进行了从头算分子动力学模拟研究 [10] ，得到的结果显示 HAsO 4 2 − 中未质子化的氧原子较好地与附近的2~3个水分子形成相对稳定的氢键，而质子化的氧原子形成1~2个氢键，弱于水–水氢键，由模拟得到的两种As-O距离(未质子氧和质子氧分别为1.68 Å和1.78 Å)，水化数为10，而Tanaka的研究 [11] 却显示 HAsO 4 2 − 周围有15个水分子，与Bhattacharjee研究结果不一致。因此，本文采用量子化学理论计算先研究了[HmAsO4]3−m (m = 0~2)周围的水化特征，进一步对铁盐除砷中水化作用的微观机制进行讨论，为日后更好地除去饮用水中的砷提供参考。

 2. 计算方法

本文使用密度泛函理论计算方法研究了Fe3+与砷物种HmAsO4m−3缔合对HmAsO4m−3水化作用的情况。首先，通过利用B3LYP/6-311G(d,p)方法优化了一系列[HmAsO4(H2O)12]m+ (m = 0~2)稳定的水化团簇。进一步优化了考虑30个水分子时Fe3+与HmAsO4m−3 (m = 0~2)缔合结构，其中，对Fe采用了Stuttgart赝势基组，对于As、H、O原子，则采用了6-311G(d,p)基组。也采用B3LYP/6-311G(d,p)方法优化[FeHmAsO4(H2O)30]m+水合团簇，分析了它们的具体结构参数。对所有优化而得的几何构型，进一步进行freq分析，都是无虚频。条件是温度为298 K，1 atm下，水溶液条件采用的是SMD模型，都是有且只有一个虚频。所有的量子化学计算中都采用了Gaussian 09软件包。

 3. 结果与讨论 3.1. H<sub>m</sub>AsO<sub>4</sub>水化特征

砷酸盐的水化特征直接与水溶液中除砷密切相关。在此节中，我们使用B3LYP/6-311G(d,p)方法首先对[HmAsO4(H2O)12]3−m (m = 0~2)水合团簇进行微观结构分析。 图1 列出了优化后得到的最稳定的[HmAsO4(H2O)12]3−m (m = 0~2)水合团簇的几何构型。

图1. 优化而得的稳定的[HmAsO4(H2O)12]3−m (m = 0~2)团簇构型

为了便于分析，将砷酸中质子化的氢标为H，与氢相连的砷酸中的氧记为OP，不与氢相连的砷酸中的氧记为ON。如 表1 所列出，[AsO4(H2O)12]3−中As-ON键的平均键长是1.73 Å；[HAsO4(H2O)12]2−中有一个OP，其与砷的距离为1.88 Å，As-ON的平均距离为1.69 Å，H-OP的距离为0.97 Å；[H2AsO4(H2O)12]−中有两个O被质子化，As-OP的平均距离1.75 Å，As-ON的平均距离为1.70 Å，H-OP的平均距离为0.99 Å。由此看来，随着砷酸的质子化，As-ON和As-OP的键长在减短，H-OP的键长增长。

<xref></xref>Table 1. Partial bond lengths of [HmAsO4(H2O)12]3−m hydration clusters calculated by B3LYP/6-311G(d,p) methodTable 1. Partial bond lengths of [HmAsO4(H2O)12]3−m hydration clusters calculated by B3LYP/6-311G(d,p) method 表1. B3LYP/6-311G(d,p)方法计算的[HmAsO4(H2O)12]3−m水合团簇的部分键长

物种

As-ON (Å)

As-OP (Å)

H-OP (Å)

[AsO4(H2O)12]3−

1.73

-

-

[HAsO4(H2O)12]2−

1.69

1.88

0.97

[H2AsO4(H2O)12] −

1.70

1.75

0.99

对[HmAsO4(H2O)12]3−m (m = 0~2)水合团簇中[HmAsO4]3−m周围水分子的分布情况进行了概率密度函数(PDF)及配位数分析。如 图2 是[HmAsO4(H2O)12]3−m (m = 0~2)水合团簇中As与水中氧(OW)的概率密度函数及配位数； 图3 为[AsO4(H2O)12]3−中砷酸根中ON与周围水分子中氢(HW)的概率密度函数及配位数； 图4 为[HAsO4(H2O)12]2−中砷酸根中ON与周围水分子中氢(HW)的概率密度函数及配位数； 图5 为[H2AsO4(H2O)12]−中砷酸根中ON与周围水分子中氢(HW)的概率密度函数及配位数；As-OW和O-HW的配位数是通过在PDF的第一配位层的最小值处的距离r作为相互作用的截断距离来确定的。

图2. 优化而得的稳定的[HmAsO4(H2O)12]3−m (m = 0~2)中As-OW键的PDFs

图3. 优化而得的稳定的[AsO4(H2O)12]3−中O-HW键的PDFs

根据 图2(a) ，我们得到[AsO4(H2O)12]3−中As-OW对应的截断距离为4.4 Å，说明在r < 4.4 Å范围内的As-OW均属于第一配位水层，配位数是11，最高峰对应的距离是As-OW在第一配位层内的最佳距离为3.6 Å。根据 图2(b) ，可以得出[HAsO4(H2O)12]2−中As-OW对应的截断距离也是4.4 Å，配位数为8，最佳距离是3.8 Å，与Bhattacharjee等人 [10] 用从头算分子动力学模拟得到[HAsO4]2−的As-OW距离一致；对应Bhattacharjee等人的从头算分子动力学模拟分析中的径向分布函数图，可发现其As-OW第一配位层是将 图2(b) 中前两个峰的拟合，结合 图1 中[HAsO4(H2O)12]2−的键长参数发现，他们将未直接与[HAsO4]2−中O形成氢键，而是通过该水分子的2个H，同时与HAsO42−中O形成氢键的2个水分子的O形成氢键。与As距离为4.72 Å，若将这两个水分子也考虑在内，[HAsO4(H2O)12]2−中As-OW配位数为10，与Bhattacharjee等人的从头算分子动力学模拟结果一致，而与Tanaka等人 [11] 的所得配位数为15不一致。3.4 Å处的小峰为同时与[HAsO4]2−中两个O配位的水分子与As的距离。根据 图2(c) ，能够得到[H2AsO4(H2O)12]−中As-OW对应的截断距离为4.4 Å，最佳距离也是3.8 Å，配位数为8，与Tanaka等人所得配位数基本吻合。随着氢质子化，[HmAsO4(H2O)12]3−m (m = 0~2)水合团簇中As-OW的配位数减小，As-OW键长变长，说明随着氢质子化，[HmAsO4(H2O)12]3−m的水化减弱。

图4. 优化而得的稳定的[HAsO4(H2O)12]2−中O-HW键的PDFs

图5. 优化而得的稳定的[H2AsO4(H2O)12]−中O-HW键的PDFs

结合 图3 ， 图4 和 图5 ，[AsO4(H2O)12]3−中四个未质子氢ON-HW的平均配位数是2.5，截断距离也是2.6 Å，最佳距离都是1.6 Å；[HAsO4(H2O)12]2−中有一个氢，三个ON-HW的平均配位数是2.6，截断距离都是2.4 Å，平均最佳距离是2.0 Å，OP-HW配位数为2，截断距离都是2.4 Å，最佳距离是2.0 Å；[H2AsO4(H2O)12]−中有两个氢，两个ON-HW的平均配位数为2.5，截断距离也是2.4 Å，平均最佳距离是1.8 Å，两个OP-HW的平均配位数为2，截断距离是2.8 Å，平均最佳距离是2.0 Å。随着氢质子化，[HmAsO4(H2O)12]3−m (m = 0~2)水合团簇中ON-HW距离小于OP-HW，而ON-HW配位数大于OP-HW。

 <xref></xref>3.2. Fe<sup>3+</sup>与H<sub>m</sub>AsO<sub>4</sub><sup>m−3</sup>缔合对H<sub>m</sub>AsO<sub>4</sub><sup>m−3</sup>水化的影响

为了比较Fe3+与HmAsO4m−3缔合对HmAsO4m−3水化的影响，在这里我们使用B3LYP/6-311G(d,p)方法优化Fe3+和HmAsO4m−3缔合结构[FeHmAsO4(H2O)30]m+ (m = 0~2)水合团簇进行微观结构分析(见 图6 )。

图6. 优化而得的稳定的[FeHmAsO4(H2O)30]m+ (m = 0~2)团簇几何构型

图7. 优化而得的稳定的[FeHmAsO4(H2O)30]m+(m = 0~2)水合团簇中As-OW键的PDFs

为了有利于分析Fe3+与HmAsO4m−3缔合对HmAsO4m−3水化的影响，我们比较了As-OW的配位数在缔合前后的变化。在上节中我们已对[HmAsO4(H2O)12]3−m (m = 0~2)中的As-OW键的配位数进行统计。按照同样的方法对缔合后的[FeHmAsO4(H2O)30]m+ (m = 0~2)水合团簇(如 图7 所示)进行分析。 图7(a) 中[FeAsO4(H2O)30]0的As-OW的PDF中第一配位层出现了2个峰，总的配位数为13，第一峰(3.6 Å)对应着通过氢键直接与 AsO 4 3 − 中氧配位的水中OW与As的距离，配位数为6，而4.4 Å的峰对应是与Fe配位且与ON配位的水成氢键桥的水分子与As的距离。 图7(b) 中[FeH2AsO4(H2O)30]2+的As-OW的PDF中有同样的情况。[FeHAsO4(H2O)30]+的As-OW的PDF中第一峰(3.6 Å)配位数为7，第一水化层总配位数为13。[FeH2AsO4(H2O)30]2+的As-OW的PDF中第一峰(3.8 Å)配位数为8，第一水化层总配位数为15。可见，不管是只考虑与HmAsO4m−3中氧配位的水还是将As-OW距离小于5.2 Å(根据As-OW的PDF)的水分子都考虑在内，随着质子化，[FeHmAsO4(H2O)30]m+ (m = 0~2)中HmAsO4m−3物种水化数都增加，即随着质子化，[FeHmAsO4(H2O)30]m+ (m = 0~2)中HmAsO4m−3物种水化增强，与[HmAsO4(H2O)12]3−m (m = 0~2)出现相反的情况。可见，Fe3+与HmAsO4m−3缔合改变了HmAsO4m−3水化情况。而相对于[HmAsO4]3−m (m = 0~2)，[FeHmAsO4]m+ (m = 0~2)中[HmAsO4]3−m的水化能力降低，从而达到沉淀除砷。

图8. 优化而得的稳定的[FeHmAsO4(H2O)30]m+ (m = 0~2)水合团簇中Fe-Ow键的PDFs

 <xref></xref>3.3. [FeH<sub>m</sub>AsO<sub>4</sub>]<sup>m+</sup>团簇的水化特征

相对于HmAsO4，[FeHmAsO4]m+不仅考虑As-OW的水化数，还有考虑Fe-Ow的水化数，所以我们统计的水化数是As-OW和Fe-Ow综合考虑的结果。根据 图8(a) ，我们得到[FeAsO4(H2O)12]0中Fe-Ow第一配位水层对应的截断距离为2.4 Å，第一配位水层的配位数是5，最高峰对应的距离是Fe-Ow在第一配位层内的最佳距离为2.0 Å；根据 图8(b) ，可以得出[FeHAsO4(H2O)12]+中Fe-Ow对应的截断距离为2.4 Å，配位数为5，最佳距离是2.0 Å；根据 图8(c) ，能够得到[FeH2AsO4(H2O)12]2+中Fe-Ow对应的截断距离为2.4 Å，配位数为5，最佳距离也是2.0 Å。如 图7 所示，[FeAsO4(H2O)30]0、[FeHAsO4(H2O)30]+、[FeH2AsO4(H2O)30]2+水合团簇中As-OW的配位数前面也已经分析过，按[FeAsO4(H2O)30]0、[FeHAsO4(H2O)30]+、[FeH2AsO4(H2O)30]2+中As-OW的配位数依次为6、7、8。通过VESTA直观观察发现，则共用水个数依次是4、4、4，不受质子化影响。那么[FeAsO4(H2O)30]0、[FeHAsO4(H2O)30]+、[FeH2AsO4(H2O)30]2+中[FeHmAsO4]m+ (m = 0~2)第一配位层水分子数目依次为10、11、12，即随着质子化，[FeHmAsO4]m+ (m = 0~2)的水化增强，归因于[FeHmAsO4]m+ (m = 0~2)的水化能顺序依次是[FeH2AsO4]2+ > [FeHAsO4]+ > [FeAsO4]0 [12] 。

 3.4. Fe<sup>3+</sup>与 <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> <msub> <mi> H </mi> <mi> m </mi> </msub> <mi> A </mi> <mi> s </mi> <msubsup> <mi> O </mi> <mn> 4 </mn> <mrow> <mstyle mathvariant="bold" mathsize="normal"> <mi> m </mi> </mstyle> <mo> − </mo> <mn> 3 </mn> </mrow> </msubsup> </mrow> </math>缔合对 <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> <msub> <mi> H </mi> <mi> m </mi> </msub> <mi> A </mi> <mi> s </mi> <msubsup> <mi> O </mi> <mn> 4 </mn> <mrow> <mstyle mathvariant="bold" mathsize="normal"> <mi> m </mi> </mstyle> <mo> − </mo> <mn> 3 </mn> </mrow> </msubsup> </mrow> </math>结构的影响

根据前人颜丹 [13] 所研究的[FeHmAsO4(H2O)n]m+水合团簇都倾向于以六配位的形式存在。所以在n = 30时，Fe3+与[HmAsO4]m−3缔合形成的水合团簇中Fe3+倾向于六配位，且Fe3+与[HmAsO4(H2O)n]m−3配位形式倾向于单齿结构。根据 表2 中所示[FeAsO4(H2O)30]0、[FeHAsO4(H2O)30]+、[FeH2AsO4(H2O)30]2+水合团簇中不同类As-O进行分析。[FeAsO4(H2O)30]0、[FeHAsO4(H2O)30]+、[FeH2AsO4(H2O)30]2+水合团簇中的As-ON分别有3、2、1个，对应的平均键长依次是1.68 Å、1.67 Å、1.66 Å，与 表1 中[HmAsO4(H2O)12]m+ (m = 0~2)的As-ON进行对比，Fe3+与[HmAsO4]m+ (m = 0~2)的缔合使As-ON键长略变短；[FeHAsO4(H2O)30]+、[FeH2AsO4(H2O)30]2+水合团簇中的As-OP分别有1、2个，平均键长依次是1.81 Å、1.74 Å，与[HmAsO4(H2O)12]m+ (m = 0~2)对比，As-OP键长也略变短。[FeHmAsO4(H2O)30]m+水合团簇都是As-OFe均只有1个，键长依次是1.81 Å、1.73 Å、1.71 Å。即随着质子化，[FeHmAsO4(H2O)30]m+水合团簇中As-OFe距离变短。[FeHAsO4(H2O)30]+、[FeH2AsO4(H2O)30]2+水合团簇中H-OP键分别是0.97 Å、0.99 Å，与 表1 中[HmAsO4(H2O)12]m+ (m = 0~2)的H-OP进行对比，[HAsO4(H2O)12]2−、[H2AsO4(H2O)12]−水合团簇中H-OP键也分别是0.97 Å、0.99 Å，即Fe3+与[HmAsO4]m+ (m = 1~2)的缔合不影响H-OP键。

<xref></xref>Table 2. The As-O bond length of the [FeH<sub>m</sub>AsO<sub>4</sub>(H<sub>2</sub>O)<sub>30</sub>]<sup>m</sup><sup>+</sup> hydrated cluster is calculated by the B3LYP/6-311G(d,p) methodTable 2. The As-O bond length of the [FeHmAsO4(H2O)30]m+ hydrated cluster is calculated by the B3LYP/6-311G(d,p) method 表2. B3LYP/6-311G(d,p)方法计算的[FeHmAsO4(H2O)30]m+水合团簇的As-O的键长

物种

As-ON (Å)

As-OP (Å)

As-OFe (Å)

H-OP (Å)

[FeAsO4(H2O)30]0

1.68

−

1.81

-

[FeHAsO4(H2O)30]+

1.67

1.81

1.73

0.97

[FeH2AsO4(H2O)30]2+

1.66

1.74

1.71

0.99
4. 结论

铁盐除水溶液中砷时，不可避免与阴阳离子、缔合团簇的水化作用相关。本文利用量子化学方法优化了[HmAsO4(H2O)n]m−3 (m = 0~2)水合团簇和[FeHmAsO4]m+ (m = 0~2)水合团簇的优势物种，分析它们水化特征得出：

(1) 受质子化的影响，随着[HmAsO4(H2O)n]m−3 (m = 0~2)水合团簇中氢离子数量的增加，ON-HW距离小于OP-HW，而ON-HW配位数大于OP-Hw；HmAsO4m−3水化作用减弱。

(2) 相对于[HmAsO4]3−m (m = 0~2)物种，[FeHmAsO4]m+ (m = 0~2)缔合团簇中[HmAsO4]3−m的水化能力降低，从而达到沉淀除砷。随着质子化，[FeHmAsO4]m+ (m = 0~2)整个缔合团簇的水化则增强。

 基金项目

郑州师范学院科研启动专项经费资助(2018)；河南省高等学校大学生创新创业训练计划项目(202412949004)；郑州师范学院大学生创新创业训练计划项目(DCY2023002)。

 NOTES

*通讯作者。

 References 张德超, 张耀阳, 罗劲松, 等. 含砷酸性废水除砷工艺现状与展望[J]. 云南冶金, 2023, 52(2): 49-58. Bennett, G.F. (2002) Environmental Chemistry of Arsenic. Journal of Hazardous Materials, 92, 213-215. >https://doi.org/10.1016/s0304-3894(02)00012-2 杨富糈, 狄军贞, 王明玉. 多源数据的水土污染风险评价管控研究进展[J]. 环境科学与技术, 2024, 47(8): 131-145. 陈芳芳, 陈辉. 用双氧水氧化-铁盐沉淀法从含砷废水中去除砷试验研究[J]. 湿法冶金, 2022, 41(4): 372-375. 邵谱生. 电凝法在处理湿法钨冶炼含砷废水中的研究与应用[D]: [硕士学位论文]. 南昌: 南昌航空大学, 2014. 邓志华, 李嘉艳, 伏鸭生, 等. 含砷工业废水净化技术研究进展[J]. 环境科学导刊, 2017, 36(S2): 152-155. 李文辉, 唐业梅. 酸性高砷废水处理方法的探讨[J]. 湖南有色金属, 2001(6): 40-41. 方兆珩, 石伟, 韩宝玲, 等. 高砷溶液中和脱砷过程[J]. 化工冶金, 2000(4): 359-362. 许根福. 处理高砷浓度工业废水的化学沉淀法[J]. 湿法冶金, 2009, 28(1): 12-17. Bhattacharjee, A., Pribil, A.B., Lim, L.H.V., Hofer, T.S., Randolf, B.R. and Rode, B.M. (2010) Structural and Dynamic Aspects of Hydration of Haso 4 −2: An ab Initio QMCF MD Simulation. The Journal of Physical Chemistry B, 114, 3921-3926. >https://doi.org/10.1021/jp911860y Tanaka, M., Takahashi, Y., Yamaguchi, N., Kim, K., Zheng, G. and Sakamitsu, M. (2013) The Difference of Diffusion Coefficients in Water for Arsenic Compounds at Various Ph and Its Dominant Factors Implied by Molecular Simulations. Geochimica et Cosmochimica Acta, 105, 360-371. >https://doi.org/10.1016/j.gca.2012.12.004 Yan, D., Li, H., Cai, H., Wang, M., Wang, C., Yi, H., et al. (2018) Microscopic Insight into Precipitation and Adsorption of As(V) Species by Fe-Based Materials in Aqueous Phase. Chemosphere, 194, 117-124. >https://doi.org/10.1016/j.chemosphere.2017.11.150 颜丹. 溶液中离子缔合与吸附的理论计算研究[D]: [硕士学位论文]. 长沙: 湖南大学, 2017.
Baidu
map