基于密度泛函理论的聚丙烯酸钠对钨矿的吸附性能研究

基于密度泛函理论的聚丙烯酸钠对钨矿的吸附性能研究

高海拔地区是一个重要的战略资源。含钨矿物主要有黑钨矿、白钨矿和黑白混合钨矿,易选的黑钨矿已基本枯竭,但由于白钨矿多属细粒嵌布的矿石,适于采用浮选进行分离回收,广大选矿科技工作者在这方面已进行了大量的研究。白钨矿中往往含有与其可浮性相似的含钙脉石矿物,如萤石和方解石等,从而使得白钨矿与萤石、方解石难于分离。矿物的浮选行为取决于矿物的性质,而矿物的性质取决于矿物的电子结构,因此从白钨矿、萤石和方解石3种矿物电子结构方面找到它们的内在差异,并指导浮选实践就显得非常重要。近年来随着计算机技术以及量子化学理论的发展,通过计算机模拟矿物晶体结构及其与药剂的相互作用,从而揭示脉石矿物与目的矿物结构的差异以及与药剂间的作用机制,并指导浮选实践变为了现实,基于密度泛函的第一性原理是研究矿物晶体结构及界面性质的有力工具,在硫化矿浮选理论研究中应用比较广泛,但采用第一性原理研究白钨矿、萤石和方解石浮选性质差异的报道较少。本文采用密度泛函第一原理研究了白钨矿、萤石和方解石晶体及表面电子结构性质,并模拟了水体系中聚丙烯酸钠在3种矿物表面吸附行为的差异,研究结果对于从本质上认清3种矿物的可浮性差异具有重要的意义,从而为研究白钨矿与萤石、方解石浮选分离提供理论参考。1计算参数结果对比基于密度泛函的第一性原理方法,采用软件MaterialStudio4.4中的CASTEP模块,对白钨矿、萤石和方解石3种含钙矿物的能带结构、电子态密度和Mulliken布居进行计算。为了使计算晶格常数和禁带宽度接近实验值,首先对3种含钙矿物的原胞模型进行优化处理,以选取较佳的交换关联函数和平面波截断能。计算所选取的其他参数均采用平面波截断能和所选取的函数所对应的默认值,3种矿物的原胞模型如图1所示。通过对几种函数计算结果的对比,局域密度泛函近似(LDA)下的CA-PZ交换关联函数计算出的晶格常数与实验值最为接近,因而3种含钙矿物的交换关联函数均采用LDA下的CA-PZ交换关联函数。截断能的计算结果如表1所示。表1的平面波截断能测试表明:白钨矿的截断能为277eV时较为合理,其晶格参数的计算值a=b=0.52372nm与实验值a=b=0.52429nm的误差为0.1%;萤石的截断能为270eV时,其晶格参数的计算值a=b=c=0.54624nm与实验值a=b=c=0.54631nm的误差仅为0.01%;方解石的截断能为275eV时,模拟计算得出的晶格参数值为a=b=0.49868nm,与实验值a=b=0.4988nm的误差仅为0.02%。3种矿物的计算结果与实验结果的误差均很小,表明计算所采用的方法以及选取的参数是可靠的。2结果与讨论2.1禁带宽度的计算白钨矿、萤石和方解石的能带结构如图2所示,取费米能级(EF)作为能量零点。计算结果表明:白钨矿的禁带宽度为4.184eV,与MinoruItoh等计算出的4.25eV比较接近;萤石的禁带宽度为6.902eV;方解石的禁带宽度为4.835eV,与AndrewJSkinner等计算出的(4.4±0.2)eV比较接近。模拟计算出的禁带宽度高于或低于实验值主要是由于GGA近似下的DFT,对电子与电子之间的交换关联作用处理不足引起的。由于半导体的禁带宽度一般在2eV以下,绝缘体的禁带宽度则更大,因而从能带结构计算结果可知,白钨矿、萤石和方解石均属于绝缘体。2.2种矿物的表面ca白钨矿、萤石和方解石3种矿物的态密度如图3所示。从图3可以看出,白钨矿、方解石费米能级附近的态密度主要由氧的2p轨道构成,而萤石主要由氟的2p轨道构成。白钨矿、方解石、萤石晶体中Ca的态密度组成非常相似,3种矿物在-40eV附近的态密度主要由Ca3s轨道组成,而处于-20eV附近的态密度主要由Ca3p轨道贡献。这就使得矿物表面Ca的物化性质相似,因此在浮选过程中矿物表面的Ca表现出相似的化学活性,难以分离。3种矿物态密度的微小差别主要在于:(1)在-75eV附近白钨矿W5s轨道对其态密度有小的贡献,而其他两种矿物则没有,但由于其所处位置的能量太低,一般认为其难以参加化学反应;(2)白钨矿和方解石费米能级附近的态密度主要由氧的2p轨道构成,萤石费米能级附近的价带主要由F2p轨道组成,由于处于费米能级附近的能态组成具有较高的化学活性,因此,白钨矿和方解石在参与化学反应时是O的活性较强,而萤石则是F的活性较强;(3)萤石和方解石的导带底部的态密度主要是Ca3d轨道贡献,而白钨矿导带底部态密度主要由W5d轨道贡献。以上这些3种矿物态密度的差别均有可能成为它们选择性分离的依据。2.3氟原子所带负荷的计算Mulliken重叠布居是Mulliken提出的一种表示电荷在各组成原子之间分布情况的方法。通过Mulliken布居分析可以考察模拟体系的电荷分布、转移和所形成的键的性质等情况。白钨矿的O原子、Ca原子和W原子在优化前的价电子构型为O2s22p4,Ca3s23p64s2,W5s25p65d4,优化后原子的Milliken布居值如表2所示。从表2中可知:白钨矿优化后的价电子构型为O2s1.862p4.84,Ca3s2.133p6.003d0.53,W5s2.305p6.545d3.71。Ca原子和W原子为电子供体,钙原子所带电荷为+1.35e,钨原子所带电荷为+1.46e。O原子为电子受体,所带电荷为-0.70e。萤石的F原子和Ca原子在优化前的价电子构型为F2s22p5,Ca3s23p64s2,优化后的价电子构型为F2s1.962p5.71,Ca3s2.163p6.004d0.51,Ca原子所带电荷为+1.32e,为电子供体,F原子所带电荷为-0.66e,为电子受体。方解石的C原子、O原子和Ca原子在优化前的价电子构型为C2s22p2,O2s22p4,Ca3s23p64s2,优化后的价电子构型为C2s0.832p2.41,O2s1.802p4.914s0.42,Ca3s2.143p6.003d0.49,Ca原子和C原子为电子供体,钙原子所带电荷为+1.37e,碳原子所带电荷为+0.76e,O原子为电子受体,所带电荷为-0.71e。综合表2的数据分析可知:白钨矿、萤石和方解石中的钙原子均带正电,所带电荷大小顺序为方解石>白钨矿>萤石;白钨矿和方解石中的氧原子带负电,钨和碳原子带正电;萤石中的氟原子带负电。键的Mulliken布居值能体现出键的离子性和共价性的强弱,布居值越大表明键的共价性越强,越小则表明键的离子性越强。白钨矿、萤石和方解石键的Mulliken布居值列于表3中,由表3中数据分析可知:白钨矿和方解石中的O-Ca键及萤石中的F-F键均呈现出较大的离子性,易产生离子键断裂;白钨矿中的O-W键和方解石中的C-O键主要以共价性为主,但由于处于基团结构中,在破碎磨矿过程中难以断裂,因此白钨矿和方解石在破碎过程中仍主要为离子键断裂;从键的布居值可以看出,白钨矿中Ca-O键和方解石中CaO键的布居相差很小,只有0.001nm,从而它们表面键的断裂情况相似,白钨矿和方解石表面都有价键未饱和的钙离子,萤石表面的F-F键最容易断裂,但氟离子的水化自由能比表面钙离子的小,容易优先水化进入溶液,使得表面也产生价键不饱和钙离子,3种矿物表面的化学活性非常相似,因而,在浮选过程中表现出相似的浮选性能。2.4真空层厚度的表面结构的收敛效果表4列出3种矿物表面能最小时的原子层数和真空层厚度,单胞模型如图4所示。白钨矿(111)表面含有6层原子及1.2nm的真空层厚度的表面结构能够给出较满意的收敛结果,此时的表面能为0.7202J·m-2;萤石(111)表面含有3层原子及1.0nm的真空层厚度的表面结构能够给出较满意的收敛结果,表面能为0.5517J·m-2;方解石(104)表面含有21层原子及1.2nm的真空层厚度的表面结构能够给出较满意的收敛结果,表面能为0.4350J·m-2。2.5单矿物浮选试验采用密度泛函理论模拟水分子和聚丙烯酸钠分子在白钨矿(111)、萤石(111)和方解石(104)表面的吸附,将水分子和聚丙烯酸钠分别置于3种矿物表面不同位置上,通过优化构型和计算吸附能,吸附能最小的吸附位置为最稳定的吸附模式。吸附能为吸附前后各物质中能量的变化,其大小可以表示吸附体系的稳定性。药剂在矿物表面的吸附能根据下式计算:式中,Eads为吸附能(kJ·mol-1);E(矿物+药剂)为吸附后体系的总能量(kJ·mol-1);E矿物为矿物表面的能量(kJ·mol-1);E药剂为药剂的能量(kJ·mol-1)。吸附能为负,代表吸附为放热反应,且其绝对值越大,说明药剂在表面的作用越强,吸附过程越容易发生。由于聚丙烯酸钠分子链太长在计算中难以收敛,因此选取它的部分官能团结构单元进行模拟。水分子和聚丙烯酸钠分子在矿物表面最稳定吸附模式的吸附能分别如表5,6所示。由表5和6的结果可以看出,水在白钨矿、萤石和方解石表面的吸附均较弱,而聚丙烯酸钠在3种矿物表面的吸附更强,且更容易发生。图5和6表明,水分子与白钨矿作用时吸附在同一WO3基团的两个氧原子上的吸附能最小为-225.70kJ·mol-1,与萤石作用时吸附在不同电子层的两个氟原子上的吸附能最小为-112.03kJ·mol-1,与方解石作用时吸附在同一CO3基团的两个氧原子上的吸附能最小为-110.52kJ·mol-1。聚丙烯酸钠分子与白钨矿作用时吸附在同一电子层的钙原子上的吸附能最小为-359.27kJ·mol-1,吸附在萤石(111)面同一电子层的钙原子上的吸附能最小为-405.32kJ·mol-1,聚丙烯酸钠吸附在方解石(104)面同x轴的钙原子上(CP1)的吸附能最小为-289.06kJ·mol-1。水分子和聚丙烯酸钠分子在3种矿物表面吸附后的模型分别如图5,6所示。由于聚丙烯酸钠与3种矿物是在水体系中发生吸附,吸附前必然要排挤开矿物表面已经吸附的水化膜,因此聚丙烯酸钠在水溶液中与3种矿物表面吸附能力大小可以定义为:根据式(2)进行吸附能的计算,结果见表7。从计算结果看出3种矿物在水体系中与聚丙烯酸钠发生反应的吸附能均为负值,其绝对值大小顺序为萤石>方解石>白钨矿。说明聚丙烯酸钠在3种矿物表面的吸附使它们受到抑制作用,因此聚丙烯酸钠在自然pH(未添加pH调整剂)条件下对3种矿物抑制能力顺序为萤石>方解石>白钨矿。为了进一步验证模拟计算结果,以白钨矿、萤石和方解石纯矿物为对象,以聚丙烯酸钠(25mg·L-1)为抑制剂、碳酸钠/盐酸为pH调整剂、731(75mg·L-1)为捕收剂进行了单矿物浮选试验,矿浆pH值对3种矿物可浮性的影响如图7所示。图7的浮选结果表明:采用聚丙烯酸钠做抑制剂、731为捕收剂时,随着pH值的增加,3种矿物的回收率均先增加后降低;pH=8.7~10范围内,白钨矿的回收率明显高于萤石、方解石的回收率,此时白钨矿的回收率均保持在80%左右,萤石的回收率在15%~10%之间,方解石的回收率低于40%。从而可以说明,聚丙烯酸钠是萤石和白钨矿的有效抑制剂,有可能实现白钨矿与萤石、方解石的浮选分离。在自然pH(未添加pH调整剂)条件下,白钨矿、萤石和方解石浮选时的pH值分别为8.78,8.20和9.25,对应的回收率分别为76.03%,18.59%和33.68%,聚丙烯酸钠对3种矿物的抑制能力强弱顺序为萤石>方解石>白钨矿,与模拟计算结果相符。3聚丙烯酸钠在水体系中的吸附能力和作用机理1.白钨矿、萤石和方解石最优的模拟函数均为LDA-CA-PZ函数,截断能分别为277,270和275eV;能带结构表明3种矿物均属于绝缘体。2.白钨矿、萤石和方解石3种矿物的Ca原子的态密度组成很相似,白钨矿和方解石费米能级附近的态密度主要由氧的2p轨道构成,萤石费米能级附近的价带主要由F2p轨道组成,因此,白钨矿和方解石在参与化学反应时是O的活性较强,而萤石则是F的活性较强。3.白钨矿和方解石表面都有价键未饱和的钙离子,萤石表面的F-F键最容易断裂,但氟离子的水化自由能比表面钙离子的小,容易优先水化进入溶液,使得表面也产生价键不饱和钙离子,在浮选过程中表现出相似的浮选性能。4.在水体系中,聚丙烯酸钠对3种矿物吸附能绝对值的大小顺序为白钨矿>方解石>萤石,分别为-133.57,-178.54和-293.19kJ·mol-1,说明聚丙烯酸钠在3种矿物表面的吸附使它们受到抑制作用,且抑制力强弱顺序为萤石>方解石>白钨矿。5.在纯矿物浮选试验中,当pH=8.7~10时,聚丙烯酸钠对萤石和方解石