正文:摘要:运用量子化学密度泛函 (DFT)中B3LYP方法,在6-31G*基组水平上,对六种氨基膦酸化合物的缓蚀性能和分子结构进行研究,讨论了量子化学计算结果和缓蚀性能的关系,结果表明:分子的轨道能量与缓蚀性能有很好的相关性,氨基膦酸的缓蚀性能主要与苯环的得电子有关。
关键词:氨基膦酸化合物 缓蚀 量子化学
中图分类号:TQ314
文献标识码:A
Quantum chemistry studies on mechanism of corrosion inhibition for amido phosphonic acids
Abstract:The structure and activity relationships(SAR) of amido phosphonic acids were studied systemically by using Density Functional Theory(DFT) method at B3LYP/6-31G* level. Results of quantum chemical calculations indicated that the orbital energies have great relationships with the corrosion efficiency and the abilities of incepting electron of the benzene ring in amido phosphonic acids has great relationship of the corrosion inhibition efficiency.
Key words: amido phosphonic acids;corrosion inhibition;quantum chemistry
0 引言缓蚀剂是以适当浓度和形式存在于环境介质中,可以防止或减缓腐蚀的一种化学物质或几种化学物质的混合物,其中有机膦酸是一种高效、价廉的缓蚀剂,近年来研究表明
[1-2],将具有大π键的芳环引入有机膦酸中,不但能降低膦含量,抑制细菌的生长,且能提高缓蚀效果。缓蚀剂的缓蚀效率与分子结构有密切的关系,采用量子化学计算方法可以得到表征缓蚀剂分子结构参数,本文选取的六种氨基膦酸化合物中不仅含有膦酸基团,还含有具有大π键的苯环和亲水性好,分散能力强的-SO
3或-SO
2NH
2,运用量子化学密度泛函理论(DFT)中B3LYP方法,在6-31G*基组水平上,对六种氨基膦酸化合物的缓蚀性能和分子结构关系进行研究。
1 计算方法上述六种物质的初始构型由ChemOffice2004
[3]中的Chem3D Ultra 8.0构建,在不施加任何对称性限制的条件下,采用半经验PM3方法对分子构型进行预优化(包括键长、键角以及二面角等)。然后在Gaussian03
[4]的程序中,运用密度泛函理论(DFT)的B3LYP方法,选用6-31G*基组,对预优化后的分子进行几何构型全优化,并进行频率计算和分子轨道的计算,收敛精度取程序内定值,所得结构均为势能面上的极小点(无虚频),全部计算均在Pentum Ⅳ计算机上完成。
2 计算结果与讨论2.1分子构型图1 六种物质的分子结构
Fig 1 The structure of the six compounds
2.2 量子化学参数和缓蚀效率的关系2.2.1 前线轨道能量和缓蚀性能的关系量子化学的前线轨道理论
[5]认为:有机物的前线分子轨道,即最高占据轨道(HOMO)和最低空轨道(LUMO)的性质,决定着其反应活性和反应机理。HOMO轨道的能量
EHOMO代表分子轨道中最高被占轨道的能量,是化合物分子给电子能力的量度。
EHOMO越高,说明原子核对HOMO轨道上电子的吸引力越弱,这些电子越容易供出;
EHOMO越小,该轨道中的电子越稳定,分子的给电子能力越小。
ELUMO与分子的电子亲和能直接相关,其值越小,电子进入该轨道后体系能量降低的越多,该分子接受电子的能力越强。轨道能量差(能隙)Δ
E =
ELUMO-
EHOMO是非常重要的稳定性指标,其值越大稳定性越好,在化学反应中的活性越差。前人对含N化合物的研究
[6-8]均认为物质的缓蚀率和前线轨道能量有很好的相关性。
据文献
[9]介绍,分子的电负性X=-(
EHOMO+
ELUMO)/2,分子的硬度n=(
ELUMO-
EHOMO)/2,计算数据见表1,其中六种物质的缓蚀率来自文献
[10],缓蚀率数值取自药剂量为100mg/L,pH=8的条件下对Q235钢的缓蚀率,将轨道能量参数与缓蚀效率用SPSS软件做线形回归,发现缓蚀率发现缓蚀率
η=0.706
EHOMO-2.903
ELUMO-0.287,相关系数
r=0.950,相关性较好,根据方程中参数前的系数可以发现分子的得电子能力对缓蚀性能比供电子能力影响较大,E
LUMO值越小,物质越容易得电子,缓蚀效果越好。
缓蚀率
η=2.288
X+8.190
n-0.309,相关系数
r=0.950.,相关性较好,根据方程中参数前的系数我们可以发现分子的硬度对缓蚀性能影响较大,说明缓蚀剂在作用时服从软硬酸碱理论。
表1六种物质的缓蚀效率和轨道能量
Tab 1 The corrosion inhibition efficiency and orbital energy of the six compounds
物质 |
η/% |
EHOMO/eV |
ELUMO/eV |
△E/eV |
X/eV |
n/eV |
M1 |
65.73 |
-0.23857 |
-0.09501 |
0.14356 |
0.16679 |
0.07178 |
M2 |
80.54 |
-0.22677 |
-0.03312 |
0.19365 |
0.129945 |
0.096825 |
M3 |
76.59 |
-0.21957 |
-0.03108 |
0.18849 |
0.125325 |
0.094245 |
M4 |
73.91 |
-0.21994 |
-0.02552 |
0.19442 |
0.12273 |
0.097219 |
M5 |
82.75 |
-0.23876 |
-0.03679 |
0.20197 |
0.137775 |
0.10099 |
M6 |
76.14 |
-0.21503 |
-0.01823 |
0.1968 |
0.11663 |
0.09840 |
2.2.2 原子电荷和缓蚀性能的关系这六种物质的结构相似,符合通式:我们把左边的苯环上所有碳原子的电荷之和记为
1/3 1 2 3 下一页 尾页