趙曉洋， 唐 靜， 劉 偉， 肖 信， 南俊民*

(1. 河南工業(yè)職業(yè)技術學院環(huán)境工程系， 南陽 473009； 2. 華南師范大學化學學院， 廣州510006)

在化學清洗工業(yè)中，通常采用的無機酸具有揮發(fā)性強、氣味濃、對設備腐蝕性強等缺點，而有機酸通過螯合作用溶解污垢，具有酸性小、對金屬腐蝕性小、無毒無污染等優(yōu)點. 乙二胺四乙酸 (EDTA)是一種氨基多羧酸螯合劑，能夠與Fe、Co、Ni、Cu等金屬離子形成穩(wěn)定且溶于水的配合物，已被廣泛應用于工業(yè)、農(nóng)業(yè)、醫(yī)藥、科研等領域[1]. EDTA可作為螯合清洗劑，具有安全、高效、工期短等優(yōu)勢[2]. EDTA與污垢中金屬離子有很強的螯合作用，可以有效清除鈣垢、鐵垢、銅垢等，對金屬基底腐蝕小，可在金屬表面生成保護膜，除垢后無需后續(xù)鈍化保護過程. 乙二胺四亞甲基膦酸(EDTMP)是一種氨基多膦酸螯合劑，已被廣泛應用于水處理. 在水溶液中，EDTMP與Ca2+、Mg2+、Fe3+等可形成松散的螯合物，破壞水垢的正常結(jié)晶，有效抑制碳酸鹽、硫酸鹽、磷酸鹽及水合氧化鐵的沉淀，可阻止硬垢的生成.

垢樣取自平頂山某鹽化有限責任公司的列管式鈦換熱器. 換熱室的殼程材料為16MnR鋼，管程材料為鈦管. 在生產(chǎn)中，鹵水物料進入管程被蒸發(fā)濃縮，過熱蒸汽進入殼程被加熱，運行溫度為143 ℃. 由于NaCl鹽分對蒸汽管道有較強的腐蝕作用，殼程加熱室容易結(jié)生灰色鐵垢，具有鐵磁性，主要成分為Fe3O4. 在高溫下生成的鐵垢非常致密難溶于水，不溶于無機酸溶液. 由于管程使用鈦管，不能使用氫氟酸；而無機酸如鹽酸、硝酸、磷酸、硝酸+鹽酸、硝酸+硫酸均不能溶解鐵垢，只能采用有機酸才能溶解鐵垢. 而在篩選可用于化學清洗的有機酸時發(fā)現(xiàn)：含有膦酸基的羥基乙叉二膦酸(HEDP)效果好[3]. EDTMP含有膦酸基且與EDTA結(jié)構(gòu)相似，通過溶解鐵垢實驗發(fā)現(xiàn)，EDTA可以溶解鐵垢，而EDTMP不能溶解鐵垢. 雖然已有量子化學計算研究[4-7]報道了EDTA和EDTMP與過渡金屬元素生成配合物的相關研究，但是它們與Fe3+配合的分子結(jié)構(gòu)對比研究尚未報道. 使用量子化學計算方法，從分子結(jié)構(gòu)的角度，解釋EDTA清洗鐵垢的效果強于EDTMP的原因.

1 計算方法

采用Gaussview 6.0軟件建立Fe(EDTA)-和Fe(EDTMP)-配合物的分子結(jié)構(gòu)，使用Gaussian 16軟件基于密度泛函理論進行計算[8]. 利用B3LYP方法，對所有原子采用DEF2-TZVP基組. 以水為溶劑，在基于密度的溶劑模型(Solvation Model Based on Density,SMD)中，對2種配合物的幾何結(jié)構(gòu)進行構(gòu)型優(yōu)化和頻率計算，計算收斂閾值取程序默認值. 所得結(jié)構(gòu)振動分析沒有虛頻，對應于勢能面上的穩(wěn)定結(jié)構(gòu). 對所得結(jié)構(gòu)采用Multiwfn 3.6程序[9]進行波函數(shù)分析，包括：幾何構(gòu)型分析、擴展電荷分解分析(Extended Charge Decomposition Analysis,ECDA)、約化密度梯度(Reduced Density Gradient,RDG)分析和拓撲(Atoms in Molecular,AIM)分析. 使用NBO 3.1程序進行自然鍵軌道(Natural Bond Orbital,NBO)分析. 使用(Visual Molecular Dynamics，VMD)軟件[10]對幾何構(gòu)型、RDG、AIM、NBO進行圖形化繪制.

2 結(jié)果與討論

2.1 配合物的結(jié)合能和擴展電荷分解分析

EDTA具有較強的鰲合能力，其結(jié)構(gòu)中2個sp3雜化的N原子和4個羧基O原子都參與配位，可與大多數(shù)金屬離子形成物質(zhì)的量之比為1∶1的配合物. EDTMP與EDTA結(jié)構(gòu)相似，僅有羧基和膦酸基的區(qū)別. EDTMP與EDTA都可與Fe3+形成配合物，但EDTA對Fe3+的鰲合能力比EDTMP的強. 采用量子化學方法，有望從2種配合物的分子結(jié)構(gòu)差異來揭示這一現(xiàn)象的本質(zhì)原因.

對2種配合物在低自旋和高自旋狀態(tài)下的結(jié)構(gòu)分別進行了優(yōu)化，配合物的能量和結(jié)合能見表1. 從配合物的結(jié)合能可以看出2種配合物結(jié)合的緊密程度，結(jié)合能越大，配體與中心離子結(jié)合得越緊密.

表1 配合物的結(jié)合能和轉(zhuǎn)移電荷Table 1 The binding energies and charge transfer of the complexes

注:Ec為配合物的能量；El為配體的能量；EFe3+為Fe3+的能量；Eb為配合物的結(jié)合能.

分別計算配合物、配體和Fe3+的能量，由公式

Eb=Ec-El-EFe3+

可得配合物的結(jié)合能. 可以看出，同種配合物的高自旋態(tài)能量比低自旋態(tài)能量略低，但同種配合物的低自旋態(tài)結(jié)合能比高自旋態(tài)的大，且結(jié)合能的能量差遠大于高、低自旋態(tài)的能量差. Fe(EDTA)-的高自旋態(tài)能量比低自旋態(tài)能量低0.816 eV，但低自旋態(tài)的結(jié)合能比高自旋態(tài)的低5.388 eV，F(xiàn)e(EDTMP)-的高自旋態(tài)能量比低自旋態(tài)能量低1.061 eV，但低自旋態(tài)的結(jié)合能比高自旋態(tài)的低5.143 eV. 結(jié)合能數(shù)據(jù)說明低自旋態(tài)有利于配體與Fe3+的結(jié)合. 在高溫狀態(tài)下，F(xiàn)e3+由高自旋轉(zhuǎn)態(tài)變?yōu)榈妥孕龖B(tài)，發(fā)生了自旋交叉的現(xiàn)象，并與配體形成配合物[11].

擴展電荷分解分析可以定量描述配體向金屬離子的電荷轉(zhuǎn)移數(shù)量[12]. 從表1中的轉(zhuǎn)移電荷數(shù)據(jù)可以看出，同類配合物的低自旋態(tài)配體向Fe3+轉(zhuǎn)移的電子數(shù)均大于高自旋態(tài)的，表明低自旋態(tài)更有利于配體向Fe3+轉(zhuǎn)移電子. 在低自旋態(tài)下，從EDTA向Fe3+轉(zhuǎn)移的電荷數(shù)為2.708 e，從EDTMP向Fe3+轉(zhuǎn)移的電荷數(shù)為2.632 e，說明EDTA的轉(zhuǎn)移電荷數(shù)比EDTMP的多.

2.2 配合物的分子結(jié)構(gòu)

優(yōu)化后的低自旋配合物分子結(jié)構(gòu)如圖1所示. Fe(EDTA)-的結(jié)構(gòu)由1個EDTA分子的4個羧基和2個N原子參與配位，每個羧基提供1個O原子和2個N原子，與中心Fe3+形成1個六配位的八面體結(jié)構(gòu). 對于Fe(EDTMP)-，其結(jié)構(gòu)與Fe(EDTA)-的相似，不同的是羧基變成了膦酸基. 盡管EDTA和EDTMP分子與Fe3+形成的配合物結(jié)構(gòu)相似，但二者溶解鐵垢的性能存在差異，說明二者分子的差異性起了重要作用. 從表2可以看出，F(xiàn)e(EDTA)-比Fe(EDTMP)-的在八面體平面位置的對應鍵長略短. Fuzzy鍵級反映了2個原子間共享的電子對數(shù). 從鍵級上可以看出，F(xiàn)e(EDTA)-的鍵級比Fe(EDTMP)-同類型鍵的大. 通過鍵長和鍵級的比較可以說明：EDTA與Fe3+結(jié)合得更加緊密，這與結(jié)合能的結(jié)果一致.

圖1 配合物的優(yōu)化構(gòu)型

表2 低自旋態(tài)配合物的幾何結(jié)構(gòu)參數(shù)和Fuzzy鍵級Table 2 The geometry parameters and the Fuzzy bond orders ofthe low-spin complexes

EDTA和EDTMP的分子結(jié)構(gòu)相似，差異在于EDTA有4個羧基，而EDTMP有4個膦酸基. 從表3中可以看出，由于C原子的電負性(2.55)比P原子的電負性(2.19)大，通過對電荷進行自然布局分析(NPA)，EDTMP的膦酸基配位O原子的負電荷比EDTA的羧基配位O原子的負電荷多. 從單個基團的電荷分析，EDTMP上配位O原子的給電子能力應該比EDTA上配位O原子的更強，但實際效果并非如此，如果分析導致EDTMP溶解鐵垢能力不如EDTA的強的因素，單方面用電荷多少判斷結(jié)合力大小來分析并不合理，還需要考慮其他因素，例如空間位阻、拓撲結(jié)構(gòu)、軌道相互作用等.

表3 配合物中部分原子的NPA電荷Table 3 The partial NPA charges of the complexes e

2.3 配合物的空間位阻分析

約化密度梯度是一種可視化研究弱相互作用的方法[13]，不僅可反映范德華力、氫鍵和靜電等弱相互作用，還能顯示空間位阻作用. 弱相互作用臨界點的電子密度ρ(r)是衡量相互作用強度的重要指標之一，其數(shù)值大小和鍵的強度存在正相關性，可用sign(2)函數(shù)表示鍵的類型. 鍵臨界點的sign(2)=-1，環(huán)、籠臨界點的sign(2)=+1. 將電子密度ρ(r)和sign(2)的乘積投影到RDG等值面上，可以可視化研究弱相互作用的強度與類型. 2種配合物分子的RDG等值面和RDG散點圖如圖2所示.

圖2 配合物的RDG等值面和散點圖

RDG等值面和散點圖的顏色一致，圖中藍色區(qū)域代表起吸引作用的弱相互作用，綠色區(qū)域代表相互作用強度弱的范德華力；由于該區(qū)域電子密度很小，所以可正可負，其中綠色區(qū)域代表弱吸引作用、棕色區(qū)域代表弱互斥作用. 紅色區(qū)域代表環(huán)、籠中原子間的互斥作用.

從圖2中可以看出，2種配合物五元環(huán)內(nèi)都有紅色區(qū)域，表現(xiàn)為互斥作用. 2種分子的亞甲基上H原子和配位O原子間有綠色的弱相互吸引作用. 但2種配合物的配位基團之間相互作用并不相同，F(xiàn)e(EDTMP)-的棕色區(qū)域明顯大于Fe(EDTA)-的，表明Fe(EDTMP)-中相鄰的膦酸基之間表現(xiàn)為弱互斥作用. 這是由于EDTA中羧基C原子為sp2雜化，EDTMP中膦酸基P原子為sp3雜化，sp3雜化比sp2雜化具有更大的空間伸展，導致Fe(EDTMP)-的膦酸基之間表現(xiàn)為弱位阻互斥作用，而Fe(EDTA)-的羧基之間位阻互斥作用更弱.

散點圖中X、Y軸分別為sign(2)ρ(r)和RDG函數(shù)，將sign(2)ρ(r)范圍定義為[-0.050～0.050],其中，sign(2)ρ(r)在[-0.025～-0.015]范圍內(nèi)表示強吸引，在[-0.015～0.015]范圍內(nèi)表示范德華力作用，在[0.015～0.025]范圍內(nèi)表示空間位阻作用. 從圖2可以看出，在五元環(huán)中間有紅色梭形區(qū)域，體現(xiàn)較強的位阻效應，對應于散點圖最右邊的紅色峰. Fe(EDTMP)-在散點圖[-0.015～0.015]范圍內(nèi)存在多個峰值，棕色區(qū)域體現(xiàn)膦酸基之間的弱位阻互斥作用，對應于[0.000～0.015]之間的峰；綠色區(qū)域體現(xiàn)配位O原子和亞甲基H原子之間存在微弱的范德華力，對應于[-0.015～0]之間的峰.

2.4 配合物的拓撲分析

BADER等[14]提出分子中的原子(Atoms in Molecules,AIM)理論，將化學結(jié)構(gòu)、化學鍵等與電子密度分布函數(shù)的拓撲性質(zhì)聯(lián)系起來，定量描述了分子中原子及原子間化學鍵. 圖3給出了2種配合物的拓撲分析圖. 圖中紫色點是(3,-3)，對應于核臨界點，與原子核的位置基本一致；桔色點是(3,-1)，對應于鍵臨界點(Bond Critical Point,BCP)，位于2個有作用力的原子之間；黃色是環(huán)臨界點. AIM理論認為，鍵臨界點處的電子密度和能量密度與化學鍵性質(zhì)存在密切聯(lián)系. 表4列出2種配合物的配位鍵、1個羧基和1個膦酸基的鍵臨界點性質(zhì).

圖3 配合物的拓撲分析圖

表4 配合物的鍵臨界點的電子密度拓撲性質(zhì)Table 4 The topological properties of bond critical points of the complexes

鍵臨界點位置的電子密度ρ(r)和勢能密度V(r)與化學鍵強度有密切關系. 對于同類化學鍵，通常ρ(r)越大、V(r)越負，則化學鍵強度越大. 由表4可知，F(xiàn)e(EDTA)-中配位的4個O原子和2個N原子的ρ(r)均比Fe(EDTMP)-中的大，其V(r)則均比Fe(EDTMP)-的更負，說明Fe(EDTA)-的同類配位鍵強度均大于Fe(EDTMP)-的情況. 羧基鍵臨界點的ρ(r)明顯比膦酸基鍵臨界點的大，這是由于羧基的OCO形成了離域π鍵，而膦酸基的P—O鍵為σ鍵. 羧基的離域程度大，影響了配位鍵臨界點的電子密度，F(xiàn)e(EDTA)-中鍵臨界點的ρ(r)比Fe(EDTMP)-同類型鍵的大. 能量密度H(r)體現(xiàn)了某個點的電子能量，是動能密度G(r)與勢能密度V(r)之和. 當H(r)<0時，化學鍵是共價作用，而當H(r)>0時，化學鍵是非共價作用[15].H(r)/ρ(r)的物理意義是BCP位置上單位電子的能量密度，H(r)/ρ(r)越負則共價作用越強，而H(r)/ρ(r)越正則非共價作用越弱[16]. 羧基比膦酸基的配位能力更強(表4)，2種配合物的配位鍵臨界點的H(r)和H(r)/ρ(r)均為負值，說明配位鍵為共價作用，且Fe(EDTA)-的配位鍵臨界點的H(r)/ρ(r)均比Fe(EDTMP)-的更負，說明Fe(EDTA)-配位鍵的共價作用更強. 在2種配合物中，F(xiàn)e—N配位鍵的H(r)/ρ(r)均比Fe—O配位鍵的更負，說明Fe—N配位鍵的共價作用強于Fe—O配位鍵. 這是因為N原子的電負性小于O原子，N原子束縛電子的能力小于O原子，其電子云具有更大的延伸性，與Fe3+形成共價鍵的共價作用更強.

2.5 配合物的自然鍵軌道分析

使用NBO程序?qū)ε浜衔镞M行自然鍵軌道(Natural Bond Orbital,NBO)分析. 表5列出部分電子供體(Donor)軌道i、電子受體(Acceptor)軌道j以及由二階微擾理論得到軌道之間相互作用的穩(wěn)定化能E(2).E(2)越大表明它們之間的相互作用越強，即軌道i供電子給軌道j的傾向越大，電子的離域化程度越大[17].

表5 配合物的主要電子供體和受體軌道的E(2) Table 5 The E(2) values of the main electron donor and acceptor orbits of the complexes

2種配合物分子的Fe—N反鍵軌道與其對面Fe—O反鍵軌道有最大離域作用，與其它Fe—O配位鍵的反鍵軌道也有離域作用，且Fe(EDTA)-的同類型E(2)均比Fe(EDTMP)-的大(表5). 羧基上O孤對電子與C—O反鍵軌道離域作用，明顯強于膦酸基上O孤對電子與P—O反鍵軌道的相互作用，這與拓撲分析的結(jié)果一致. 圖4展示了上述有強相互作用的軌道相互作用圖，NBO軌道之間發(fā)生了重疊，使其軌道能量降低，這是其軌道穩(wěn)定化能E(2)較大的原因. 拓撲分析和NBO分析均表明，EDTA中羧基形成離域π鍵，使得金屬離子與配位原子間電子離域化程度變大，金屬離子與配體共價作用增強，配位鍵強度增大. EDTMP由于膦酸基沒有離域π鍵，不利于配位鍵共價作用的增強.

圖4 配合物NBO軌道的相互作用圖

3 結(jié)論

探討了在鐵垢的化學清洗過程中，只有EDTA可以溶解鐵垢，而EDTMP不能夠溶解鐵垢的原因. 通過量子化學計算研究了兩種配合物Fe(EDTA)-和Fe(EDTMP)-的分子結(jié)構(gòu)，配合物為低自旋態(tài)的六配位八面體結(jié)構(gòu). 通過配合物分子的幾何構(gòu)型、電荷分布、擴展電荷分解分析、約化密度梯度分析、拓撲分析、自然鍵軌道分析發(fā)現(xiàn)，EDTA的羧基形成的離域π鍵，加強了配位鍵的相互作用，而EDTMP的膦酸基只有σ鍵. Fe(EDTA)-中配位鍵的共價作用強于Fe(EDTMP)-. EDTMP膦酸基之間的空間位阻不利于EDTMP與Fe3+結(jié)合. 結(jié)構(gòu)上的差異導致EDTA與Fe3+的結(jié)合能力強于EDTMP. 本文從分子結(jié)構(gòu)層面解釋了EDTA常用于化學清洗，而EDTMP不能用于化學清洗的原因，為下一步的化學清洗分子設計奠定了基礎.