中国有色金属学报 2004,(08),1437-1441 DOI:10.19476/j.ysxb.1004.0609.2004.08.032
苯氧乙酸类浮选抑制剂性能的量子化学计算
华东理工大学生物化学研究所,中南大学资源加工与生物工程学院,中南大学资源加工与生物工程学院,中南大学资源加工与生物工程学院 上海200237中南大学资源加工与生物工程学院长沙410083 ,长沙410083 ,长沙410083 ,长沙410083
摘 要:
运用密度泛函理论(DFT)B3LYP方法,取6 31G 基组,对苯氧乙酸类浮选抑制剂分子进行了量子化学计算。通过官能团羧基的键电荷密度及其比值ρ(O—H)/ρ(C—O)的计算,比较了各药剂分子离解成一价负离子的难易。根据各分子的前线轨道能量差ΔEHOMO LUMO和一价负离子的HOMO能级,对分子的反应活性、一价负离子与矿物表面不同金属离子作用的强弱进行了比较。从药剂与矿物表面作用前后的能量变化,比较了各药剂的一价负离子与不同矿物表面发生作用的稳定性,并通过对比矿物与捕收剂作用后的能量变化,证明苯氧乙酸类抑制剂对方解石的抑制作用是靠抑制剂在矿物表面强烈的化学吸附而形成的。
关键词:
中图分类号: TD923
作者简介:张剑锋(1968),男,讲师,博士.电话:02164253156;E mail:zjf_sh@163.com;
收稿日期:2003-12-26
基金:国家重点基础研究发展规划资助项目(G19990649011);
Quantum chemical calculation on properties of phenoxy acetic acids depressants
Abstract:
Quantum chemistry of phenoxy acetic acids depressants was studied at B3LYP/6-31G* level by density functional theory(DFT). The ratio of bond charge density(ρ(O—H)/ρ(C—O)) for dissociating carboxylic group in each reagent was put out as the judgment for dissociation extent. Frontier orbital energy calculation of the reagents and their corresponding univalence anions approved the traditional theory that molecules in lower ΔEHOMO-LUMO has greater reactivity, and that univalence anion in lower HOMO energy causes greater reactivity with active metal ion on the mineral surface. Energy variation between the univalence anion and the mineral surface reveals that such reagents make intensive adsorption on the surface of calcite and diaspore.
Keyword:
phenoxy acetic acid; density functional theory; bond charge density; quantum chemistry calculation;
Received: 2003-12-26
苯氧乙酸类化合物具有多种用途, 已越来越引起人们的重视
图1 苯氧乙酸类化合物的结构图
Fig.1 Structures of phenoxy acetic acids
1 方法
量子化学计算和分析在美国MSI公司的Cerius2软件包上进行。 量子化学计算顺序如下: 先按有关标准结构数据
2 结果与讨论
2.1药剂的离解性能与羧基键电荷密度比ρ(O—H)/ρ(C—O)的关系
苯氧乙酸类化合物属于有机羧酸或有机多元羧酸, 羧基是它们的活性官能团, 羧基的离解情况直接关系到酸性的强弱, 也关系到与其他物质作用的难易。
对于无机含氧酸的强度规律曾有许多研究工作
f=ρ(O—H)/ρ(X—O)
f值愈小, 酸性愈强。 为了找到一条评判苯氧乙酸类化合物酸性强弱的标准, 计算了化合物1~7分子中各个羧基的键电荷密度ρ值, 并将其中最易离解羧基的C—O单键和O—H键的键电荷密度及其比值(亦即各个分子中最小的羧基ρ(O—H)/ρ(C—O)值)和一级酸离解常数列于表1。
表1 苯乙酸类分子中离解羧基的键电荷密度比
Table 1 Ratio of bond charge density fordissociating carboxyl group inphenoxy acetic acid1
Compound |
ρ(O—H) | ρ(C—O) | f2 | pKa |
1 |
0.967 097 | 1.138 82 | 0.849 21 | 3.15 |
2 |
0.951 368 | 1.146 42 | 0.829 86 | 2.85 |
3 |
0.967 822 | 1.132 87 | 0.854 31 | 3.49 |
4 |
0.967 269 | 1.134 68 | 0.852 46 | 3.18 |
5 |
0.967 650 | 1.133 65 | 0.853 57 | 3.32 |
6 |
0.968 589 | 1.133 94 | 0.854 18 | 3.47 |
7 |
0.968 063 | 1.132 86 | 0.854 53 | 3.53 |
1—Listed value of bond charge density and f in 2-7 is the smallest one in each compound's carboxylic group(exclusive for compound 1); 2—f=ρ(O—H)/ρ(C—O);3—pKa1 is negative logarithm of first order acidic dissociation constant for compounds 2-7, while negative logarithm of acidic dissociation constant for compound 1.
从表1可以看出, 化合物1~7中各离解羧基的f值与它们的pKa1值的大小顺序是完全对应的, 其中f值最小的化合物2具有最小的pKa1值和最强的酸性, 化合物7的f值最大, 其对应的pKa1最大, 具有最弱的酸性。 因此由量子化学计算得到的羧基键电荷密度比ρ(O—H)/ρ(C—O) 值完全可以作为比较苯氧乙酸类化合物酸性强弱的依据。
由此可知, 苯氧乙酸类抑制剂在水溶液中容易离解出氢离子, 变成一价负离子, 从而与矿物表面的金属离子作用并吸附于矿物表面, 如果抑制剂分子本身的亲水性足够大, 则将会对矿物产生有效的抑制。
2.2 前线轨道能量与反应活性
按照化学反应性的前线分子轨道理论, 过渡态的形成是由于反应物的前线轨道的相互作用导致的。 对于化合物本身来说, HOMO-LUMO能级差ΔEHOMO-LUMO又是一个重要的稳定性指标(又叫稳定化能), ΔEHOMO-LUMO值大意味着具有高稳定性, 在化学反应中有低的反应活性, 而ΔEHOMO-LUMO值小时, 则易给出电子, 具有高的反应活性
表2的结果表明, 药剂分子6和7的HOMO与LUMO之间能量差要小于其它的药剂分子, 故它们的反应活性在苯氧乙酸类分子中相对来说要大些, 而药剂分子1的前线轨道能量相差最大, 故它的反应活性最低, 其它几种药剂分子的前线轨道能级差和反应活性则介于这两者之间。
在浮选过程中, 苯氧乙酸类抑制剂在弱碱性条件下对矿物的抑制能力强, 这时的抑制剂分子大多是以离子形式存在, 离子的状态参数对药剂与矿物的作用有着密切关系。 各药剂一价负离子的前线轨道能级对应地列于表2中。
根据软硬酸碱理论, 药剂负离子与矿物阳离子作用时属于碱的范畴, 其HOMO能级越低, 其硬度也就越大。 由表2可以看出: 所有苯氧乙酸类抑制剂负离子的HOMO能级均较低, 它们都是较硬的碱, 但相对说来药剂7的负离子最硬, 药剂6的负离子次之, 其余的5种药剂负离子相对要软一些。 由于方解石、 一水硬铝石和黄铁矿表面的阳离子作为酸的硬度大小顺序为
由此可知: 尽管药剂分子6和7在水溶液里的离解相对要稍难一点, 但是在碱性体系中变成一价负离子以后, 则成了较硬的碱, 故与硬酸Ca2+的作用变得容易, 从而产生对方解石很强的抑制作用。
2.3药剂与矿物表面作用前后的能量变化
一水硬铝石的晶体模型按其结构数据在工作站上自建, 方解石和黄铁矿的晶体模型来自工作站本身的数据库。 分别选取每一种矿物的2个晶胞作为计算对象。 矿物与药剂离子作用的表面为矿物的解理面, 分别为
由表3可知, 所有苯氧乙酸类化合物与方解石
表2 苯氧乙酸类化合物分子及其一价阴离子的前线轨道能量
Table 2 Calculation of HOMO and LUMO energy for molecule and anion of reagents(eV)
Compound(anion) |
EHOMO | ELUMO | ΔEHOMO-LUMO |
1(1-1) |
-9.048 06(-4.415 59) | -1.121 91(1.671 34) | -7.926 15(-6.086 93) |
2(2-1) |
-8.912 28(-4.310 83) | -1.204 23(1.811 87) | -7.708 05(-6.122 70) |
3(3-1) |
-8.571 72(-4.273 28) | -1.187 35(1.522 22) | -7.384 37(-5.7955 0) |
4(4-1) |
-8.694 73(-4.484 98) | -1.281 50(1.676 37) | -7.413 23(-6.161 35) |
5(4-1) |
-8.577 84(-4.376 68) | -1.116 33(1.632 56) | -7.461 51(-6.009 24) |
6(6-1) |
-8.574 04(-4.583 20) | -1.194 56(1.689 16) | -7.379 48(-6.272 36) |
7(7-1) |
-8.337 58(-4.655 05) | -1.286 14(1.719 76) | -7.051 44(-6.374 81) |
Data in parentheses associate with corresponding univalence anions.
表3 苯氧乙酸类化合物一价阴离子与矿物作用后的能量变化
Table 3 Energy variation between mineralsurface and univalence anion(eV)
Reagents |
Calcite
site |
Diaspore {010} site |
Pyrite {100} site |
1-1 |
-21.049 7 | -26.783 6 | 25.912 1 |
2-1 |
-21.696 7 | -29.531 0 | 7.000 31 |
3-1 |
-21.429 0 | -29.586 8 | 4.221 6 |
4-1 |
-21.597 4 | -29.501 4 | 4.519 8 |
5-1 |
-21.406 4 | -29.672 8 | 5.007 7 |
6-1 |
-24.690 4 | -30.255 6 | 2.392 8 |
7-1 |
-28.277 4 | -29.383 0 | 1.787 7 |
Collector |
-21.784 1 | -28.425 8 | -55.949 7 |
Collector for calcite, diaspore and pyrite is oleate, lauryl amine and n-butyl-xanthate, respectively.
表3的结果还说明:1) 当用油酸钠作捕收剂浮选方解石时, 方解石
3 结论
1) 药剂分子的羧基C—O—H键电荷密度比值ρ(O—H)/ρ(C—O)可以作为判断苯氧乙酸类分子离解难易的定量依据。
2) 药剂分子6和7的前线轨道能量差ΔEHOMO-LUMO较小, 它们的反应活性较大, 与矿物表面的作用相对较强。 药剂离解后产生的一价负离子HOMO能级均较低, 它们都是较硬的碱, 与方解石、 一水硬铝石和黄铁矿表面的金属离子产生吸附程度上的差异, 从而产生对不同矿物抑制性能上的差别。
3) 药剂离子与方解石和一水硬铝石表面作用后的能量变化均为负值, 有利于吸附作用的产生, 从而产生抑制作用。 药剂与黄铁矿作用后的能量变化为正值, 不利于吸附, 故抑制作用不明显。
参考文献