工业水处理, 2021, 41(7): 100-104 doi: 10.19965/j.cnki.iwt.2020-0976

试验研究

聚天冬氨酸阻垢缓蚀性能量子化学研究

程玉山,, 吴晋英, 黄长山

Quantum chemistry study on scale-corrosion inhibition performance of polyaspartic acid

Cheng Yushan,, Wu Jinying, Huang Changshan

收稿日期: 2021-05-1  

基金资助: 河南省重点科技攻关项目.  182102210075
河南省科学院科研开发专项.  200606064

Received: 2021-05-1  

作者简介 About authors

程玉山(1980-),硕士,副研究员电话:0371-60103605,E-mail:chengyushan09@163.com , E-mail:chengyushan09@163.com

Abstract

The structure-activity relationship between molecule structure of polyaspartic acid(PASP) and its scale-corrosion inhibition performance was studied by quantum chemical calculation method. The results showed that PASP had good scale inhibition performances, the nitrogen atom on the imine group and the oxygen atom on the carboxyl group had strong negative electricity, which made the oxygen atom on the carboxyl group very easy to interact with calcium ion on the crystal surface. The relative distance between the two nitrogen and oxygen ions in the PASP molecule was very close to the calcium ion spacing on the calcite crystal surface, the excellent matching could significantly enhance the adsorption force between PASP molecule and specific crystal face, therefore, the aim that inhibiting CaCO3 scale was achieved. The static scale inhibition experiment conclusion was completely consistent with the calculation results. PASP molecule in the corrosion inhibition performance was obtained by quantitative calculation. The maximum occupied orbital energy, the lowest space orbit energy and the difference between the lowest space orbit energy and the highest occupied orbit energy of PASP molecule were obtained by quantitative calculation. The calculation results were basically consistent with the mechanism description, and the experimental results were completely consistent with the calculated results.

Keywords: quantum chemistry ; polyaspartic acid ; scale inhibition ; corrosion inhibition

PDF (0KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

程玉山, 吴晋英, 黄长山. 聚天冬氨酸阻垢缓蚀性能量子化学研究. 工业水处理[J], 2021, 41(7): 100-104 doi:10.19965/j.cnki.iwt.2020-0976

Cheng Yushan. Quantum chemistry study on scale-corrosion inhibition performance of polyaspartic acid. Industrial Water Treatment[J], 2021, 41(7): 100-104 doi:10.19965/j.cnki.iwt.2020-0976

阻垢缓蚀剂的理论与应用研究一直以来都是水处理领域和腐蚀科学领域的热点问题,尤其是对有机缓蚀阻垢剂分子结构与其性能之间关系的研究,已经由量子化学进入到“分子设计”的水平。近年来出现许多方法用于研究性能与结构间的关系1,并在某些方面取得了成功,例如线性自由能关系(LFER)或Hemmet方程、软硬酸碱理论(HSAB)和自由Wilsin和Hansch关系。量子化学方法(QM)于1971年被Vosta和Eliasek首次使用。目前为止,QM在研究分子性能、结构与反应活性等方面已取得很大成功2

根据国内外有关文献报道3-5,阻垢剂的阻垢机理大多可用沉淀物的结晶过程来解释。迄今为止,对于阻垢剂机理的研究一般采用实验方法,主要围绕晶格畸变理论、螯合增溶作用以及分散作用等三种假说进行实验解释。笔者通过对聚天冬氨酸(PASP)分子量子化学计算,得到其优化后的空间构型,PASP分子各原子的净电荷、键长、原子间的空间距离、二面角等相关参数,并利用相关数据分析PASP分子在方解石晶面的吸附过程及对阻垢性能的影响,采用国标法碳酸钙静态阻垢试验与量子化学计算结果进行对比。

有机物缓蚀性能与其最高占据轨道能量HOMO值有关,有些研究者将HOMO值与有机物缓蚀效率作了“定量结构-活性相关(QSAR)”分析,并绘制出半经验线性回归曲线,用以作为有机缓蚀剂合成及筛选的重要依据6-8。笔者通过量子化学计算出PASP分子的最高占据轨道能量EHOMO和最低空轨道能量ELUMO,并判断其亲核反应的强弱,根据轨道能量差ΔE=ELUMO-EHOMO判断PASP分子的稳定性,从而研究PASP分子的缓蚀性能。同时利用传统的试验方法对PASP缓蚀性能进行研究,与量子化学计算结果进行对比。

PASP相对分子质量大约为1 000~5 000,是典型的绿色环保型水处理剂,兼有阻垢和缓蚀双重功能,结构式如图 1所示9

图1

图1   PASP分子结构

Fig.1   Molecular structure of PASP


利用Gaussian 16和GaussView 6软件对PASP分子进行量子化学计算,将微观电子性能与宏观阻垢缓蚀性能对比分析,进一步了解阻垢缓蚀的微观机理,进而为PASP设计高效复配水处理剂提供可靠的理论依据。

1 PASP阻垢性能研究

1.1 计算方法

首先利用GaussViev6获得PASP(m=5,n=5)分子的初始构型,接着在不施加任何对称性限制的条件下,在RHF(Restricted Hartree-Fock)水平上,采用半经验(semi-empirical)PM6-SCF-MO(Self Consistent Field-Molecular Orbitals)方法(Opti+Freq)优化分子构型(包括键长、键角以及二面角等),然后在Gaussian16程序包DFT-B3LYP/6-31G水平上对预优化后的分子进行几何构型全优化,最后计算相关参数及轨道计算,收敛精度取程序内定值,所得结构均为势能面上的极小点(无虚频)10,优化后的PASP分子空间结构如图 2所示。

图2

图2   PASP的结构优化

Fig.2   Structure optimization of PASP


1.2 PASP部分原子的键长

用Gaussian16软件计算得出PASP部分原子的键长,如表 1所示。

表1   PASP部分原子的键长

Table 1  Partial atomic bond lengths of PASP

原子键C1—N2C1—C5C1—C10C5—O6C5—O7C10—C13
键长/10-10 m1.473 661.535 731.547 441.210 271.356 181.515 20
原子键C13—O14C13—N16N16—C15C15—C18C18—O19C18—O20
键长/10-10 m1.213 891.409 931.470 561.542 651.208 001.357 46
原子键C15—C23C23—C26C26—O27C26—N29N29—C28C28—C31
键长/10-10 m1.536 271.515 081.223 351.394 531.477 121.559 06
原子键C31—N40N40—C39C42—O44C39—C47C47—C50C50—O54
键长/10-10 m1.396 181.476 841.360 961.531 361.495 001.373 19

新窗口打开| 下载CSV


1.3 PASP原子上的净电荷与阻垢性能分析

用Gaussian16软件计算得出PASP分子主要原子的电荷分布,如表 2所示。

表2   PASP(m=5, n=5)部分原子的电荷分布

Table 2  Atomic charges distribution on PASP

C1N2C5O6O7C10C13N16
-0.016-0.5080.614-0.528-0.505-0.4940.628-0.537
C15C18O19O20C23C26O27N29
0.0330.588-0.522-0.537-0.4490.622-0.565-0.516
C28C31O32C34C37O38C39N40
0.0020.549-0.523-0.4200.674-0.5340.048-0.510
C42O43O44C47C50O51O52O54
0.596-0.523-0.498-0.4380.648-0.510-0.505-0.521

新窗口打开| 下载CSV


表 2可知,PASP分子总体呈电中性,羧基和亚胺基是其主要的功能基团,羧基基团上的氧原子和亚胺基基团上的氮原子是带负电荷最多的原子,所带净电荷如表 3所示。

表3   PASP(m=5,n=5)官能团上的原子所带净电荷

Table 3  Net charge of the atom on the functional group of PASP

羧基基团O6O19O38O43O51
(5个)-0.528-0.522-0.534-0.523-0.510
亚胺基N2N16N29N40
(4个)-0.508-0.537-0.516-0.510

新窗口打开| 下载CSV


PASP分子中的亚胺基基团(—NH—)上的氮原子和羧基基团(—COO—)上的双键氧原子易于与碳酸钙表面(生长面)的钙离子形成强静电相互作用,作用强度与空间匹配密度有关11。如果羧基基团(—COO—)上双键氧原子的间距与碳酸钙表面(生长面)上钙离子的间距相当,二者则具有很好的匹配性,PASP分子会在静电作用下吸附在碳酸钙晶面上,占据碳酸钙表面的活性生长点,诱导晶体晶格发生严重的畸变,进而表现出较强的抑制碳酸钙晶面生长性能,即阻垢性能。显而易见,亚胺基基团和羧基基团属于阻垢功能基团,在阻垢剂分子中引入亚胺基和羧基,有利于阻垢性能的提高。

1.4 PASP分子上的空间尺度分析

用Gaussian16计算得出PASP分子上各原子间的空间距离,如表 4所示。

表4   PASP分子上原子间的空间距离

Table 4  Space distance between atoms on PASP

原子相对位置N29—O19N29—O29N40—O19N40—O38O6—O14O7—O14
间距/10-10 m3.8613.8894.8924.0914.0954.021
原子相对位置O14—O19O20—O43O38—O44O43—O52O43—O54O44—O54
间距/10-10 m4.9393.8274.8924.9364.9123.972

新窗口打开| 下载CSV


表 5为方解石三类典型的低Miller指数晶面上两个最近邻钙离子之间的空间距离12

表5   方解石晶体典型晶面上最近邻钙离子之间的空间距离

Table 5  Space distance between the nearest calcium ions on the typical face of calcite crystal

晶面名称(104)(102)(001)
相邻空间距离/10-10 m4.0484.0484.990
4.9904.0484.990

新窗口打开| 下载CSV


表 4表 5可以看出:方解石(104)、(102)活性生长面的间距为4.048×10-10 m,能与PASP中的N40-O38(4.091×10-10 m)、O6-O14(4.095×10-10 m)和O7-O14(4.021×10-10 m)很好的匹配。同样,方解石(104)、(001)活性生长面的间距为4.990×10-10 m,能与PASP中的O14-O19(4.939×10-10 m)、O43-O52(4.936×10-10 m)和O43-O54(4.912×10-10 m)很好的匹配。根据以上数据关系可以推断出PASP分子容易吸附在方解石晶体表面以阻止晶体的成长,对碳酸钙晶体具有较强的阻垢效果;从能量方面考虑,N40、O14、O27和O32四个原子非常容易成为阻垢剂分子与晶面作用的吸附位,能够有效地抑制这些晶面的持续成长,或使晶格发生畸变,导致生成的水垢比较松软,很容易被流水冲走,从而表现出优良的阻垢性能。

1.5 PASP阻垢性能试验研究

阻碳酸钙垢性能试验条件13:Ca2+质量浓度为600 mg/L(以CaCO3计),HCO3-质量浓度600 mg/L(以CaCO3计),pH 8.0,水温55 ℃,测试时间为10 h。PASP用量与阻垢率的实验结果如图 3所示。

图3

图3   PASP阻碳酸钙垢率关系

Fig.3   Scale inhibition rates of PASP


图 3可以得出,当PASP投加质量浓度为8 mg/L时,阻垢率达到100%,对碳酸钙垢具有显著的阻垢效果。PASP分子属于柔性分子,可以通过自身优势构型的改变达到阻止多种垢晶体的生长功效,对硫酸钡和硫酸锶同样具有较好的阻垢性能14

2 PASP缓蚀性能研究

2.1 PASP分子的偶极矩和前线轨道能量

表 6为量子化学计算得到的PASP分子前线轨道能量和能隙,其中EHOMO为PASP分子最高占有轨道(HOMO)能量,ELUMO为PASP分子最低空轨道(LUMO)能量。EHOMO的负值与PASP分子第一电离能比较接近,ELUMO的负值为PASP分子的电子亲核能。因此,EHOMO值越低,表明PASP分子越难失去电子;ELUMO值越高,表明PASP分子越难得到电子。ΔEEHOMOELUMO两者之差,决定着化合物的相对稳定性。ΔE越小,稳定性越好;反之,稳定性越差。

表6   PASP的分子轨道能量

Table 6  Molecular orbital energies of PASP

偶极矩ELUMO/eVEHOMO/eVΔE/eV
5.2367-0.380 92-0.012 15-0.368 77

新窗口打开| 下载CSV


PASP分子的EHOMO为-0.012 15 eV,数值较大,表明它的第一电离能力比较强,水溶性较好,极易参与亲核反应,容易与金属表面的空d轨道形成配位键,增强金属表面化学吸附膜的稳定性15。同时,PASP分子的ΔE较小,说明分子稳定性非常好,由此可以推测出PASP分子的缓蚀性能较好。

2.2 PASP缓蚀性能研究

2.2.1 扫描电镜分析

扫描电子显微镜(SEM)测试条件:把待测碳钢试样用导电胶固定在JSM-6510扫描电子显微镜(SEM)样品台上,加速电压设定为15 kV,放大倍数设置为1 000倍,观察碳钢试样的腐蚀状况,并获得碳钢试样的扫描电镜图片。

取两个新的碳钢试片用乙醇、丙酮处理并干燥称重后,分别放入空白溶液(自来水)和PASP溶液(8 mg/L)中进行旋转挂片试验(参照GB/T 18175—2000《水处理剂缓蚀性能的测定—旋转挂片法》),之后用扫描电镜进行扫描分析,腐蚀行为结果如图 4所示。

图4

图4   碳钢试片的扫描电镜分析

Fig.4   Scanning electron microscope analysis of carbon steel specimens


空白溶液中碳钢试片表面存在明显的坑状腐蚀现象,既不平整也不光滑,腐蚀非常严重;PASP溶液中碳钢试片表面光滑平整,几乎没有腐蚀。显然PASP可以保护碳钢样品不受腐蚀,具有良好的缓蚀效果。

2.2.2 电化学分析

电化学试验条件:极化曲线测量是在室温下采用经典的三电极系统,即辅助电极、参比电极和工作电极,分别为铂电极、饱和甘汞电极和A3碳钢(总表面积为1.0 cm2)。极化扫描范围Ecorr =±0.3 V,扫描速度0.5 mV/s,测试温度为室温,用水为自来水,pH 6.4,总硬度(以CaCO3计)161 mg/L,总碱度130 mg/L。在测试前,工作电极用2~20 μm的砂纸抛光,经蒸馏水、无水乙醇和丙酮清洗,然后自然干燥,以8 mg/L PASP水溶液为工作电极的电解液测定极化曲线。

利用电化学工作站,碳钢试片在空白溶液和PASP溶液中的极化曲线如图 5所示。

图5

图5   碳钢试片在两种不同溶液中的极化曲线

Fig.5   Polarization curves of carbon steel specimens in two different solutions


不添加PASP的溶液中,腐蚀电位和腐蚀电流密度分别为E0=-0.338 V,I0=1.91×10-8 A/cm2,添加PASP的溶液中,腐蚀电位和腐蚀电流密度分别为E0=-0.307 V,I0=1.22×10-8 A/cm2。加入PASP后,腐蚀电位发生正向移动,随着自身腐蚀电位增加,腐蚀电流密度减小,碳钢的极化曲线也得到了明显的改变,由此可见,PASP是一种阳极控制为主的混合型缓蚀剂。

3 结论

用量子化学计算方法,研究了聚天冬氨酸(PASP)的分子结构和阻垢性能的构效关系。计算结果显示,在阻垢性能方面,聚天冬氨酸分子中羧基和亚胺基是其主要的功能基团,羧基基团上的氧原子和亚胺基基团上的氮原子为带负电荷最多的原子,呈较强的负电性,羧基上的氧原子极易与晶面上的钙离子发生静电交互作用,且PASP分子中2个氮氧原子间距与方解石有关晶面上钙离子间距匹配性较好,因而显著加强阻垢剂分子与特定晶面之间的吸附力,从而能有效地抑制这些晶面的持续生长,或使晶格发生畸变,导致生成的水垢比较松软,很容易被流水冲走,从而表现出优良的阻垢性能。通过静态阻垢实验验证PASP对碳酸钙垢具有较好的阻垢性能,与量子化学计算结果完全相符。

PASP分子的EHOMO数值较大,第一电离能较强且水溶性较好,非常容易参与亲核反应,容易与金属表面的空d轨道形成配位键,增强金属表面化学吸附膜的稳定性。同时,PASP分子的ΔE数值较小,说明稳定性较好,因此可以说明PASP分子的缓蚀性能比较好。扫描电镜分析结果验证了PASP分子具有良好的缓蚀性能,电化学试验验证了PASP属于阳极型缓蚀剂,与量子化学计算结果完全相符。

参考文献

刘小辉, 周集义.

亚甲基膦酸阻垢剂构效关系研究进展

[J]. 化学推进剂与高分子材料, 2004, 2 (4): 16- 18.

URL     [本文引用: 1]

毛海英, 徐章法.

聚天冬氨酸阻垢缓蚀剂的量子化学研究

[J]. 工业用水与废水, 2007, 38 (5): 12- 15.

URL     [本文引用: 1]

Ter H J H , Wong F S K E , Vreugd C H , et a1 .

Adsorption behaviour of polyelectrolytes on calcium fluoride Part Ⅱ: Molecularmodeling of the adsorption behaviour

[J]. Colloids and Surfaces A, 1999, 154, 273- 284.

DOI:10.1016/S0927-7757(98)00883-8      [本文引用: 1]

石文艳, 张曙光, 夏明珠, .

甲叉膦酸型化合物阻垢机理的量子化学研究

[J]. 水处理技术, 2006, 1 (32): 37- 40.

URL    

郭稚弧. 缓蚀剂及其应用[M]. 武汉: 华中工学院出版社, 1987: 89- 163.

[本文引用: 1]

周公度, 段连云. 结构化学基础[M]. 北京: 北京大学出版社, 1994: 252- 256.

[本文引用: 1]

Casta J M , Luch J M .

The use of quantum mechanics calculations for the study of corrosion inhibitors

[J]. Corrosion Science, 1984, 24, 929- 933.

DOI:10.1016/0010-938X(84)90113-6     

方健, 李杰.

有机膦酸化合物阻垢缓蚀性能的量子化学研究

[J]. 同济大学学报, 2002, 30 (4): 52- 53.

URL     [本文引用: 1]

王娴.

聚天冬氨酸缓蚀剂缓蚀机理分析研究

[J]. 精细化工原料及中间体, 2012, (12): 31- 34.

URL     [本文引用: 1]

王钢, 邵波, 王溯. 两种绿色水处理剂的量子化学研究[C]//江西省电机工程学会年会论文集. 南昌: 江西省电机工程学会2014: 159-164.

[本文引用: 1]

霍宇凝, 陆柱.

环保型阻垢剂聚天冬氨酸对碳酸钙结晶过程的影响

[J]. 精细化工, 2004, 21 (3): 206- 208.

URL     [本文引用: 1]

廖沐真, 吴国是, 刘洪霖. 量子化学从头计算方法[M]. 北京: 清华大学出版社, 1984: 26- 32.

[本文引用: 1]

程玉山, 孙寅聪, 师新广, .

环境友好型水处理剂阻垢缓蚀性能研究

[J]. 河南科学, 2018, 36 (11): 1715- 1722.

URL     [本文引用: 1]

陈琼婵, 张冰如, 李风亭.

聚天冬氨酸和其他阻垢分散剂的协同效应

[J]. 工业用水和废水, 2006, 37 (2): 69- 72.

URL     [本文引用: 1]

刘靖疆. 基础量子化学与应用[M]. 北京: 高等教育出版社, 2004: 313- 352.

[本文引用: 1]

/