©北京泰科博思科技有限公司   /   京ICP备09107432号-1   /   网站建设:中企动力 北京

解决方案

Solution

>
>
>
【MS应用实例】东北大学朱一民教授成果展示:基于CASTEP模拟的氟磷灰石与十二烷基磷酸酯作用机理研究
产品名称

【MS应用实例】东北大学朱一民教授成果展示:基于CASTEP模拟的氟磷灰石与十二烷基磷酸酯作用机理研究

所属分类
产品中心
联系我们
相关资料
方案详情
行业
材料
分类
Materials Studio

摘要:

为探究十二烷基磷酸酯分子结构对氟磷灰石浮选性能的影响,采用量子化学模拟手段,研究了十二烷基磷酸酯分子结构特性及其在氟磷灰石晶体表面的作用机理。采用Materials Studio软件的CASTEP模块,首先对分子态和离子态十二烷基磷酸酯的分子结构进行了优化,得到十二烷基磷酸酯分子及其阴离子的净电荷、Mulliken电荷布居、偶极距及最高占据轨道(HOMO)组成和能量。结果表明,相对于十二烷基磷酸酯分子,十二烷基磷酸酯阴离子有更强的供电子能力、范德华作用能力和反应活性。对十二烷基磷酸酯阴离子在氟磷灰石表面的相互作用模型进行了模拟计算,得到十二烷基磷酸酯与氟磷灰石之间的作用模型和吸附能,十二烷基磷酸酯阴离子与氟磷灰石表面之间的吸附能为负值,二者之间能够自发发生吸附作用,并形成化学键。通过单矿物浮选试验验证了模拟计算结果,即十二烷基磷酸酯可以作为氟磷灰石的浮选捕收剂。

 

引言:

在实际生产中,油酸是目前最常用的磷矿浮选分散性能不好、不耐低温等方面的缺点外,还有选择性差的问题。为了提高捕收剂在矿物表面吸附的选择性,科研人员从在理论和实践2个方面寻求突破。在理论研究方面,加拿大纽芬兰纪念大学的张亚辉通过总结现有浮选理论,提出了浮选药剂与矿物表面作用的“镜像对称规则”。该规则指出,矿物表面具有断裂键的金属离子倾向于与具备矿物晶体阴离子结构的浮选药剂发生作用。氟磷灰石分子式为Ca10(PO4)6F2,其晶体中具有磷氧四面体结构。按镜像对称规则,含有磷酸根或相近极性基团结构的分子易于同氟磷灰石表面的Ca离子发生相互作用。根据以上规则,十二烷基磷酸酯是氟磷灰石的潜在捕收剂之一。本文采用密度泛函理论对十二烷基磷酸酯的分子结构进行分析,并对十二烷基磷酸酯与氟磷灰石晶体表面的相互作用进行了模拟。

 

模拟方法:

1. 氟磷灰石晶胞的结构优化方法

本文的模拟计算使用CASTEP软件模块(MaterialsStudio材料与化学计算软件,BIOVIA公司),基于密度泛函理论,采用广义梯度近似方法(GGA)对氟磷灰石晶胞进行几何优化。氟磷灰石晶体的体相和表面结构优化采用超软赝势处理原子的内层电子。参与计算的原子轨道为H1s1,C2s22p2,O2s22p4,F 2s22p5,P3s23p3,Ca3s23p64s2。采用BFGS算法对模型的结构进行优化,自洽场收敛精度设为2×10-6 eV/atom。几何优化的收敛标准为:原子间相互作用力不大于0.025 eV/nm,晶体内应力不大于0.1GPa,原子最大位移不大于2×10-4 nm。所有计算都在倒易空间中进行。图1 所示为氟磷灰石晶体的XRD 模拟谱图与实际检测值对比。

由图1可以看出,在上述参数设置条件下,氟磷灰石晶体的XRD模拟谱图与实际检测值吻合较好,证明此模拟参数合理,模拟结果可靠。

 

2. 氟磷灰石表面模型优化与性质计算方法

表面能是创造物质表面时对化学键破坏的度量。晶体某晶面的表面能小于其他晶面,表明该晶面具有更好的热力学稳定性,即在矿物破坏过程中存在下来的概率更大。同一晶面,裸露不同原子时,表面的热力学稳定性可能也是不同的,即表面能不同。氟磷灰石中磷与氧以很强的共价键形成阴离子团,而钙作为活泼金属元素,与磷氧阴离子团形成离子键。表面能的计算如式(1)所示:

式中,Eslab和Ebulk分别为晶面模型和体相模型的能量;Nslab和Nbulk分别为晶面模型和体相模型的原子数;A为晶面的面积;2指沿着表面z 轴的2个表面。


表面模型由优化后的体相切出。由文献[9-10]可知(001)面和(100)面是氟磷灰石的主要解理面。在此基础上,笔者前期分别计算了氟磷灰石在(001)面和(100)面上不同终端位置时的表面能,选取其中的低能稳定构型作为表面模型。此外,本文还进行了表面厚度和真空层厚度的表面能收敛测试,以保证模拟模型的准确性。通过以上测试,确定矿物表面计算参数为:k点2×2×2,截断动能340eV,真空层厚度1.4nm,表面模型厚度1.028 7nm。

 

3. 捕收剂的结构优化方法

先用CASTEP软件包将捕收剂分子进行结构优化,优化参数与氟磷灰石晶体优化参数保持一致。对优化后的分子结构,采用上述参数进行捕收剂分子的能量计算,最后用DMol3模块计算了捕收剂分子的前线轨道能量。

 

4. 捕收剂分子与氟磷灰石矿物表面作用吸附能计算方法

采用与氟磷灰石表面结构优化相同的计算参数,对捕收剂分子与氟磷灰石矿物表面相互作用模型进行几何优化,然后对相互作用模型的总能量进行计算。捕收剂分子与氟磷灰石矿物表面的吸附能用式(2)进行计算。

 

ΔEads=Ecomplex-(Eadsorbate+Emineral)(2)式中,Eads为吸附能;Ecomplex为捕收剂分子与氟磷灰石矿物表面相互作用模型的总能量;Eadsorbate和Emineral分别为捕收剂分子的能量和矿物表面的能量。

 

计算结果与讨论:

1. 十二烷基磷酸酯的分子结构与性质计算

(1) 十二烷基磷酸酯的分子结构

十二烷基磷酸酯的分子结构如图2所示。十二烷烃基通过O原子与磷酸酯极性基团相连,极性基团中P原子另与3个氧原子成键,形成一个磷氧双键和2个羟基。

(2) 十二烷基磷酸酯的性质计算

十二烷基磷酸酯的净电荷、Mulliken 电荷布居、偶极矩计算结果如表 1所示。

 

从表1可知:分子态和离子态的十二烷基磷酸酯羰基和羧基O原子均荷负电,并且离子态时的O原子荷电值均大于其分子态。由于净电荷数越大,原子的供电子能力越强,理论上离子态十二烷基磷酸酯更有可能与矿物表面发生作用。偶极矩是正、负电荷中心间的距离和电荷中心所带电量的乘积,它可以用来度量分子的极性。偶极矩越大,分子的极性越大,则以取向力为主的范德华力越大。离子型的十二烷基磷酸酯的偶极矩(18.584D)远大于其分子形态的偶极矩(3.852D),增加的偶极矩使离子化后的十二烷基磷酸酯与矿物表面发生范德华作用的能力加强。

 

前线轨道理论认为,参与化学反应的分子其最高占据轨道(HOMO,highest occupied molecular orbital)上的电子最为活泼。HOMO上电子云密度越大的位置,越容易给出电子,化学反应就越容易在该位置发生,图3给出了HOMO的直观形状。从图3可以看出,十二烷基磷酸酯阴离子的HOMO主要由失去H的O原子、羰基O原子和与碳链相连的O原子提供,未失去H 原子的羟基对HOMO贡献不大。上述HOMO主要贡献原子是各自捕收剂与氟磷灰石表面作用的主要位点。

 

根据前线轨道理论,化合物的HOMO与最低未占据轨道(LUMO)之间能级差|ΔE1|越小,则该化合物的稳定性越差,反应活性越强。通过十二烷基磷酸酯的前线轨道能量(FMO,frontier molecular orbital)与带隙计算结果(表2)可知,十二烷基磷酸酯分子的能级差大于其离子态,即分子态十二烷基磷酸酯结构更为稳定,反应活性小,而离子态十二烷基磷酸酯的反应活性较强。

 

前线轨道理论同时提出了物质之间反应难易程度的判定方法:物质之间的HOMO与LUMO的能量差越小,他们之间的反应越容易发生。十二烷基磷酸酯阴离子与氟磷灰石表面发生反应时为供电子方,氟磷灰石为受电子方。故二者之间的能量差计算应为|ΔE2|=|Ereagent,HOMO-Eapatite,LUMO|。通过表2可知,十二烷基磷酸酯阴离子与氟磷灰石之间的能级差较小,十二烷基磷酸酯阴离子与氟磷灰石之间有发生反应的可能性。

 

2. 十二烷基磷酸酯与氟磷灰石表面相互作用模

型与吸附能计算图4为十二烷基磷酸酯阴离子与氟磷灰石表面Ca原子作用的吸附模型。值得注意的是,十二烷基磷酸酯阴离子的吸附结构考虑了氢键的作用,该结构能使十二烷基磷酸酯阴离子与氟磷灰石表面之间的吸附作用最大化。

 

模拟计算中,吸附能是体现捕收剂与矿物表面相互作用能否发生的最可靠的证据之一。若吸附能为负,表明捕收剂与氟磷灰石表面之间可以自发发生吸附作用。若吸附能为正值,则吸附作用不会自发发生。十二烷基磷酸酯阴离子与氟磷灰石表面的吸附能为负值,这表明十二烷基磷酸酯阴离子与氟磷灰石表面之间能够自发发生吸附。此外,十二烷基磷酸酯阴离子的O原子与氟磷灰石表面的Ca原子的距离(表3)小于Ca和O原子共价半径的1.15倍,即O原子与Ca原子之间形成了化学键。吸附能为负值和Ca—O成键这2方面的计算结果均说明十二烷基磷酸酯阴离子有作为氟磷灰石捕收剂的潜质。

 

 

3. 磷灰石单矿物浮选试验

为了验证模拟计算结果,以十二烷基磷酸酯为捕收剂(用量200 mg/L),在矿浆pH=7.0条件下,进行了氟磷灰石单矿物浮选条件试验(图5)。结果显示,氟磷灰石回收率为88.62%。浮选试验结果验证了模拟结果,即十二烷基磷酸酯可以作为氟磷灰石的浮选捕收剂。

讨论:

氟磷灰石的晶体结构中,Ca—O多面体呈三方柱状,附加阴离子F-填充于Ca相间排列构成的平行于晶体c 轴的通道中。晶体中,Ca与O以离子键键合,P—O键以很强的共价键结合。由于两种键的键强差异性,磷灰石在形成解理面和断裂面时,主要沿Ca—O键断裂,在矿物表面产生带正电荷的Ca2+和负电荷的O2-离子。在进行磷灰石捕收剂的结构设计时,应主要考虑捕收剂与磷灰石晶体表面,尤其是与Ca离子的相互作用强度。基于本文理论计算工作分析,磷灰石捕收剂的结构应符合以下几个特点:

(1) 亲固基团原子带较多负电荷,以期具有较强的供电子能力。

 

(2) 分子结构的极性较强,具有较大偶极矩,以期提高捕收剂的范德华作用能力。

 

(3) 捕收剂分子自身的FMO差较小,以使其具有较强的反应活性。

 

(4) 捕收剂HOMO与磷灰石LUMO的能量差较小,以期二者之间的反应易于发生。

 

结论:

以密度泛函理论为基础,采用CASTEP和DMol3两个模块,对十二烷基磷酸酯的结构性质及其在氟磷灰石表面的相互作用进行了模拟计算。结论如下:

(1) 离子态十二烷基磷酸酯的O原子荷电数大于其分子态,说明离子态十二烷基磷酸酯更有可能与矿物表面发生作用;离子态十二烷基磷酸酯的偶极矩远大于其分子态,离子化后十二烷基磷酸酯与矿物表面发生范德华作用的能力加强。

 

(2) 前线轨道计算结果表明,离子态十二烷基磷酸酯的反应活性强于其分子态;十二烷基磷酸酯阴离子与氟磷灰石之间的能量差较大,二者之间有发生反应的可能性。

 

(3) 根据相互作用模型计算结果,十二烷基磷酸酯阴离子与氟磷灰石表面之间能够自发发生吸附,并且十二烷基磷酸酯阴离子的O原子与磷灰石表面的Ca原子之间形成了化学键,十二烷基磷酸酯阴离子有作为氟磷灰石捕收剂的潜质。

 

(4) 单矿物浮选试验验证了模拟计算结果,十二烷基磷酸酯可以作为氟磷灰石的浮选捕收剂。

 

文章详情:10.19614/j.cnki.jsks.202006015

 

北京泰科博思科技有限公司是BIOVIA Materials Studio官方指定代理商,有关软件详情或者技术支持请咨询北京泰科。


电话:010-64951848

邮箱:sales@tech-box.com.cn


北京泰科为广大学习分子模拟科研人员提供了交流讨论平台,泰科建立了BIOVIA Pipeline Pilot和Materials Studio交流群,群里有专业老师解答问题,如有兴趣一起交流,欢迎来电/邮申请入群,作者期待您的参与!

 

北京泰科涉及行业

材料研发

基于BIOVIA Materials Studio材料设计平台,提供涉及电池、航空航天、国防军工、建筑、涂料涂层等多领域材料研发软件及综合解决方案

 

药物研发

针对药物设计、药物研发等提供基于Discovery Studio、COSMOLOGIC等软件的ADME、构象比对、溶剂筛选、结晶、成盐、共晶筛选、稳定性、溶解度pKa、分配系数等性质的模拟预测软件及方案

 

化工设计

面向精细化工、新能源、石油化工等领域提供精馏萃取催化剂设计、热力学性质(溶解度、粘度等)、提纯表面处理吸附等性质模拟软件平台及解决方案

 

数据挖掘

基于Pipeline Pilot提供数据搜集、数据清洗、特征工程、机器学习、流程设计等多种数据挖掘综合解决方案

 

一体化实验室

• 实验室信息管理

• 电子实验记录本/SOP执行

• 试剂耗材管理

• 仪器管理

• 数据管理

 

部分产品

量化材料类

• Crystal:固体化学和物理性质计算软件

• Diamond:晶体结构数据可视化分析

• Endeavour:强大的求解晶体结构的软件

• Molpro:高精度量化软件

• Molcas:多参考态量软件                   

• Turbomole:快速稳定量化软件

• TeraChem:GPU上运行的量化计算软件

• Spartan:分子计算建模软件

 

数据分析类

• GelComparll:凝胶电泳图谱分析软件

• SimaPro:生命周期评估软件             

• Unscrambler:完整多变量数据分析和实验设计软件      

• CSDS:剑桥晶体结构数据库

• lCDD:国际衍射数据中心数据库                

• ICSD:无机晶体结构数据库

• Pearson’s CD:晶体数据库

 

公司简介:

北京泰科博思科技有限公司(Beijing Tech-Box S&T Co. Ltd.)成立于2007年,是国内领先的分子模拟及虚拟仿真综合解决方案提供商。

 

北京泰科博思科技有限公司与国际领先的模拟软件厂商、开发团队深入合作,为高校、科研院所和企业在材料、化工、药物、生命科学、环境、人工智能及数据挖掘、虚拟仿真教学等领域提供专业的整体解决方案。用户根据需要在我们的平台上高效的进行各种模拟实验,指导实际的生产设计。

 

北京泰科博思科技有限公司拥有一支一流的技术服务团队和资深的专家咨询团队,以客户真正需求出发,服务客户,为客户创造价值。我们秉承“职业、敬业、担当、拼搏、合作”的企业精神,致力于用国际领先的软件产品和专业全面的技术支持服务,成为客户可信赖的合作伙伴。 

未找到相应参数组,请于后台属性模板中添加
暂未实现,敬请期待
暂未实现,敬请期待

更多解决方案

——

【Materials Studio应用实例】上海理工大学成果展示:氯化铁改性碳气凝胶除汞机理的实验与理论研究
【Materials Studio应用实例】上海理工大学成果展示:氯化铁改性碳气凝胶除汞机理的实验与理论研究
BIOVIA Materials Studio
【COSMOlogic 应用实例】通过分子模拟了解四丙基多功能深共晶溶剂的CO2捕获潜力
【COSMOlogic 应用实例】通过分子模拟了解四丙基多功能深共晶溶剂的CO2捕获潜力
COSMOLOGIC
【Materials Studio应用实例】卡马西平/奥卡西平双模板磁性聚合物的制备及其去除性能
【Materials Studio应用实例】卡马西平/奥卡西平双模板磁性聚合物的制备及其去除性能
Materials Studio
【COSMOlogic 应用实例】通过分子模拟了解四丙基多功能深共晶溶剂的CO2捕获潜力
【COSMOlogic 应用实例】通过分子模拟了解四丙基多功能深共晶溶剂的CO2捕获潜力
COSMOLOGIC
【COSMOlogic 应用实例】探讨不同含水量对萃取精馏分离三元共沸混合物性能的影响
【COSMOlogic 应用实例】探讨不同含水量对萃取精馏分离三元共沸混合物性能的影响
COSMOlogic
高能材料的力场优化
高能材料的力场优化
Materials Studio