钠长石的晶体化学基因特征及其可浮性预测_王燕.pdf

收稿日期2020-03-02 基金项目国家自然科学基金项目编号 51774069， 51974067， 中央高校基本科研业务专项资金资助编号 N2001035。 作者简介王燕 （1996） ， 女， 硕士研究生。 总第 528 期 2020 年第 6 期 金属矿山 METAL MINE 钠长石的晶体化学基因特征及其可浮性预测 王燕 1， 2 朱一民 1， 2 谢瑞琦 1， 2 乘舟越洋 1， 2 李艳军 1， 21 （1. 东北大学资源与土木工程学院， 辽宁 沈阳110819； 2. 难采选铁矿资源高效开发利用技术国家地方联合工程研究 中心， 辽宁 沈阳 110819） 摘要基于密度泛函理论的第一性原理， 采用Materials Studio软件的CASTEP模块对钠长石晶胞进行优化。 计算结果表明 优化采用交换相关泛函为 GGA-PBESOL， k点取样密度为 224， 平面波截断能为 700 eV。Mul- liken电荷布居表明钠长石中， O原子为电子受体， 其余原子均为电子供体； Mulliken键布居表明NaO键最易断裂， Na会溶于水而使钠长石表面荷负电， 易与阳离子捕收剂发生作用； 能带结构表明钠长石为绝缘体； 态密度分析表 明费米能级附近主要是由O原子的2p轨道贡献， 所以O原子在参与化学反应时活性较强， 是矿物浮选的活性位点； 构建钠长石 （010） 表面， 原子层数为 3、 真空层厚度为 1.5 nm， 表面 Mulliken键布居与体相键布居有明显差别， 计算 得此时表面能为1.309 7 J/m2。研究结果可以从微观层面为锂辉石与钠长石的浮选分离提供理论参考。 关键词钠长石密度泛函理论晶体结构 （010） 面可浮性 中图分类号TD923文献标志码A文章编号1001-1250 （2020） -06-081-06 DOI10.19614/ki.jsks.202006012 Crystal Chemical Genes Characteristics of Albite and Prediction of Its Floatability Wang Yan1， 2Zhu Yimin1， 2Xie Ruiqi1， 2Cheng Zhouyueyang1， 2Li Yanjun1， 2 （1. School of Resource and Civil Engineering， Northeastern University， Shenyang 110819， China； 2. National-local Joint Engi- neering Research Center of High-efficient Exploitation Technology for Refractory Iron Ore Resources, Shenyang 110819， China） AbstractBased on the first principle of the density functional theory（DFT） ，the CASTEP module of Materials Studio was used to optimize the albite crystal cell． The parameters of optimization were set as that the exchange correlation function⁃ al was GGA-PBESOL，the Brillouin zone was sampled using Monkhorst and Pack special k points of a 224 grid，and the plane wave cut-off energy was 700 eV． Mulliken charge population showed that O atom was electron acceptor and the rest at⁃ oms were electron donors． Mulliken bond population indicated that Na-O bond was the easiest to break，Nawould dissolve in water and made the surface of albite negatively charged， which was easy to be acted by cation collector． Energy band struc⁃ ture showed that albite was an insulator． Analysis of the density of states showed that the density of state near Fermi energy level was mainly attributed to 2p orbit of O， so the activity of the O atom was strong when it participated in the chemical reac⁃ tion and it was the active site of mineral flotation． The albite（010）surface was selected to studied，and the atomic layers were 3 and the thickness of the vacuum layer was 1．5 nm，there were obvious differences of Mulliken bond population be⁃ tween surface system and bulk system， and the surface energy was calculated to be 1．309 7 J/m2． The results can provide the⁃ oretical reference for flotation separation of spodumene from albite from the crystal chemical characteristics． KeywordsAlbite， Density functional theory， Crystal structure， （010）surface， Floability Series No． 528 June 2020 钠长石 （NaAlSi3O8） 为铝硅酸盐矿物， 常作为脉 石矿物存在于大多数的金属矿及非金属矿中， 特别 是在锂辉石浮选过程中， 钠长石常以锂辉石的脉石 矿物存在。由于钠长石与锂辉石均属于硅酸盐矿 物， 有较低的零电点， 表面特性相似， 因此钠长石与 锂辉石的浮选分离困难， 分离效率低 ［1］。钠长石和锂 辉石分离最常用的方法为浮选法， 矿物的晶体结构 决定其晶体化学特征， 而矿物的晶体化学特征会直 接影响矿物的可浮性， 进而影响浮选效果 ［2］。近年 来， 有很多学者从量子化学计算的微观角度研究矿 物的晶体特征。谭鑫等 ［3］采用Materials Studio软件的 DMol3模块计算了钨锰矿不同解理面的弛豫、 电荷布 81 ChaoXing 金属矿山2020年第6期总第528期 居以及表面能， 分析计算结果发现沿MnO键断裂 形成的 （010） 面是钨锰矿的理想解理面， 表面活性原 子为Mn原子。付博等 ［4］研究了白云石的电子结构， 分析计算结果发现白云石为绝缘体， Mg、 Ca原子都是 电子供体， 在浮选过程中易发生药剂吸附的位点为 Ca、 Mg和O位点等。但目前基于第一性原理对钠长 石矿物晶体特征进行深入研究的报道还较少。 本文基于密度泛函理论， 通过对钠长石体相的 Mulliken电荷布居和键布居、 能带结构、 态密度以及 （010） 面的Mulliken电荷布居和键布居、 态密度、 表面 能进行分析并预测其可浮性， 进而揭示钠长石的晶 体结构及其对浮选效果的影响， 从微观角度为钠长 石与锂辉石的浮选分离研究提供理论参考。 1计算方法与模型 1. 1晶胞模拟与计算 采用基于密度泛函理论的 Materials Studio 软件 的CASTEP模块对钠长石的体相晶胞进行几何优化， 从而确定合理的计算参数以及稳定的钠长石晶体结 构 ［5］。几何优化采用BFGS算法， 自洽迭代过程中体 系能量改变、 原子间作用力、 原子间内应力和原子位 移， 收敛标准分别设为 210-5eV/atom、 0.5 eV/nm、 0.1 GPa、 210-4nm。几何优化的基本参数设为 交 换相关泛函为广义梯度函数GGA中的PBESOL关联 系数； k点取样密度为224， 平面波截断能为700 eV。优化计算所得钠长石的晶胞参数与试验检测参 数对比如表1所示。 由表1可知， 计算所得的晶胞参数与试验检测所 得的晶胞参数较接近， 相对误差仅为0.56， 表明优化 计算所选取的参数和采用的方法具有代表性。 将优化计算后的钠长石晶胞模型的XRD图与检 测所得的钠长石纯矿物的XRD图做对比， 结果见图1。 由图 1分析可知， 模拟计算的钠长石的 XRD图 谱与钠长石纯矿物的实测XRD图谱较吻合。结合表 1和图1可说明， 优化计算后的钠长石晶胞模型与实 际钠长石纯矿物的晶体结构相符， 具有代表性， 可以 采用 Materials Studio软件的模拟计算结果表征钠长 石的晶体结构。优化后的晶胞结构见图2。 1. 2表面模型构建 在优化后的钠长石晶胞基础上， 构建钠长石表 面模型。根据文献 ［6］ 报道， 钠长石 （010） 面为最可 能的解理面。当切割层数为3， 真空层厚度为1.5 nm 时， 表面可以维持稳定的表面性能， 随着切割层数和 真空层厚度增加， 表面能等特性变化不大， 但计算时 间变长， 因此本文沿米勒指数 （010） 方向切割钠长石 晶胞， 切割层数为3， 真空层厚度为1.5 nm。对切割后 的表面进行优化， 表面优化参数选取与晶胞优化参 数一致， 对于优化后的结构再进行性质分析。 1. 3表面能计算 表面能是沿一个晶体表面的方向将矿物晶体分 成2个独立的表面所需要的能量， 它能够表征在形成 新表面时对分子之间各种化学键的破坏程度。表面 能计算公式如式 （1） 所示 Esurf Eslab- Nslab Nbulk 2A Ebulk.（1） 式中，Esurf表示表面能， J/m2；Eslab和Ebulk分别表示模拟 82 ChaoXing 2020年第6期王燕等 钠长石的晶体化学基因特征及其可浮性预测 计算所得到的表面模型和体相模型的能量， J；Nslab和 Nbulk分别表示表面模型和体相模型的原子个数；A表 示所计算表面的面积， m2。 2计算结果及讨论 2. 1钠长石体相分析 2. 1. 1Mulliken布居分析 钠长石原胞中各原子游离态的价电子构型为O 2s22p4、 Na 2s22p63s1、 Al 3s23p1、 Si 3s23p2， 经优化计算后 体相中各体相原子的Mulliken电荷布居如表2所示。 由表2可知 钠长石原胞体相优化后， 氧原子的 价电子构型有 8种， 分别为 O 2s1.852p5.34、 O 2s1.842p5.30、 O 2s1.852p5.31、 O 2s1.812p5.34、 O 2s1.862p5.39、 O 2s1.832p5.30、 O 2s1.852p5.30、 O 2s1.822p5.32， 在氧原子中， 2s轨道失去的 电子少， 2p轨道得到的电子多， 氧原子最后所带的电 荷为负， 是电子受体； 经优化后的钠原子主要是s轨 道失电子， 钠原子的总电子数为 7.86， 失去了 1.14 e 电子， 所以钠原子所带的电荷为1.14 e； 铝原子在体 相优化后的价电子构型为Al 3s0.463p0.78， 其中3s轨道 失去 1.54 e 电子， 3p 轨道失去 1.22 e 电子， 共失去 1.76 e电子， 所以铝原子所带的电荷为1.76 e； 硅原子 在体相优化后的价电子构型有 3 种， 分别为 Si 3s0.623p1.23、 Si 3s0.653p1.25、 Si 3s0.643p1.29， 硅原子的s轨道和 p轨道都存在失电子的现象， 是电子供体。 矿物的Mulliken键布居可以用来表征原子间所 成键的共价性和离子性。首先原子间成键状态的布 居值是大于0的； 若布居值小于0， 则表明原子间状态 是反成键或不成键。在成键原子中， 布居值越大， 键 的共价性越强， 反之布居值越小， 键的离子性越强。 一般情况下， 共价性越强的键其键能越大， 越难断 裂。此外， 键长越长， 键能越小， 越容易发生断裂。 钠长石体相的Mulliken键布居如表3所示。 分析表3可知 在优化后的钠长石晶胞中存在具 有不同键布居值和键长的12种SiO键、 4种AlO 键和6种NaO键； SiO键的各个键布居值相比较 AlO键和NaO键来说是最大的， 这表明在这3种 化学键中， SiO键的共价性最强， 离子性最弱， 最不 容易发生断裂； 4种AlO键的布居值最为接近或一 致； 而NaO键的布居值最小， 键长最长， 这说明Na O键的离子性最强， 而共价性最弱， 最容易发生断 裂。当 NaO 键断裂时， Na会在水中溶解而交换 H， 水中剩余的大量 OH-与矿物表面暴露出来的 Al3、 Si4发生键合作用， 这将导致在较宽的pH值范围 内， 钠长石表面都是荷负电的， 所以选用阴离子捕收 剂难以将钠长石捕收上来 ［7］。 2. 1. 2能带结构分析 禁带宽度 （Band gap） 是表征半导体的重要参量， 矿物的内部晶体结构和原子之间的相互结合性质都 会影响到禁带宽度的大小。一般当矿物的禁带宽度 83 ChaoXing 金属矿山2020年第6期总第528期 大于3 eV时属于绝缘体； 小于3 eV时属于半导体 ［8］。 定义钠长石费米能级Ef为0 eV， 钠长石体相能带结构 如图3所示。 由图3可知， 钠长石的禁带宽度为4.928 eV， 且能 带中还有电子满带， 所以钠长石属于绝缘体。 2. 1. 3态密度分析 图4所示为钠长石体相态密度分布。 由图4可知 钠长石的能带范围在-50.64～20.14 eV， 主要分为4个部分； 在-49.45～-47.83 eV之间的 价带全部由钠原子的s轨道贡献； 在-21.40～-15.71 eV之间的价带， 大多数由钠原子的2p轨道和氧原子 的2s轨道贡献； 在-9.48～0.72 eV之间的顶部价带几 乎全部由氧原子的2p轨道贡献， 而硅原子的3s轨道 和3p轨道、 铝原子的3s轨道和3p轨道只贡献极小一 部分； 在4.40～18.83 eV之间的导带能级主要是由硅 原子的3p轨道和铝原子3p轨道组成， 其余轨道也有 很小一部分的贡献。在费米能级附近的态密度主要 是氧原子的2p轨道贡献， 因为越靠近费米能级的电 子越活跃， 所以在钠长石中氧原子的活性最强。 2. 2钠长石表面分析 钠长石 （010） 面优化后的表面结构如图5所示。 对表面层原子进行Mulliken布居分析， 从而了解表面 层原子的带电及不同原子之间的成键情况， 同时对 表面能进行分析。 2. 2. 1Mulliken布居分析 钠长石 （010） 面表面层原子的Mulliken电荷布居 如表4所示。 由表4可知， 钠长石 （010） 面优化后， 表面层原子 的外电子出现了一定程度的重新排布。对比表2数 据可以发现， 表面层氧原子的2s轨道和2p轨道上的 电子发生了微弱的重新排布， 2s轨道相比于体相失 去的电子少了， 同时2p轨道得到的电子也少了， 最终 氧原子所带电荷量并无太大改变； 表面层硅原子考 察了3种， 与体相硅原子相比， 外电子存在重新排布 但不明显； 表面层钠原子和铝原子均表现出一定的 电子转移能力， 铝原子的s轨道、 p轨道均表现为失电 子， 失电子量分别为0.04 e和0.06 e， 所以总电子数表 现为减少， 最终荷电量由体相的1.76 e增大到1.86 e； 而钠原子的s轨道失去0.02 e电子， p轨道得到0.12 e 84 ChaoXing 2020年第6期王燕等 钠长石的晶体化学基因特征及其可浮性预测 电子， 最终表现为得到 0.10 e 电子， 所以电荷量由 1.14 e减小到1.04 e。 对于 Si3O3、 Si1O3、 AlO1、 AlO7 能成键 是因为本文所考察的是钠长石单晶胞， 而钠长石实 际晶体是很多单晶胞有规律的重复单元， 所以从化 学理论上来看， 以上化学键是可以成键并有其特有 的键布居值和键长。钠长石 （010） 面表面层不同种 类原子之间的Mulliken键布居如表5所示。 分析表5可知 钠长石表面层原子之间的键布居 值和键长与体相原子的布居值与键长相比， 有明显 差别； 对于表面层所考察的 SiO 键， 尤其是 Si1 O4、 Si2O7的键布居值， 相较于体相SiO键的所有 布居值来说， 更接近于1， 表明其共价性更强； 但也有 Si1O2、 Si1O5、 Si1O3、 Si3O8这些键的布居值 低于体相所有SiO键的布居值， 表明它们之间的相 互作用有所减弱； 对于表面层所考察的AlO2键布 居值为0.55， 大于体相所有AlO键的布居值， 同时 AlO2的键长为0.169 251 nm， 小于体相所有AlO 键的键长， 表明两原子间的电子云重叠密度增大， 即 AlO2键的共价性更强， 相较于体相键能增强； 而Al O7键布居值仅为0.31， 表明其共价性有所减弱； 对 于表面层所考察的NaO4键布居值为0.11， 大于体 相中所有 NaO 键的布居值， 而键长为 0.215 143 nm， 小于体相中所有 NaO键的键长， 表明 NaO4 键的离子性有所减弱。键布居结果中， NaO键仍然 为3种化学键中最易断裂键， Na会溶于水， 所以钠长 石表面荷负电， 易与胺类等阳离子捕收剂发生静电 吸附作用或分子吸附作用。 2. 2. 2态密度分析 图 6 所示为钠长石 （010） 面表层原子态密度分 布。 由图6可知， 钠长石 （010） 面表层原子的能带分 布在-53.04 eV～17.82 eV之间， 在费米面以上的导带 底主要是由硅原子的s轨道、 p轨道以及钠原子的p 轨道、 铝原子的p轨道组成， 而在费米面附近的顶部 价带主要是由表面氧原子的2p轨道提供轨道贡献， 表明在浮选过程中， 钠长石 （010） 表面上的氧原子相 对于其他原子有更高的活性， 在参与化学反应时得 失电子的能力更强， 更容易与浮选药剂发生作用， 是 矿物浮选的活性位点。 2. 2. 3表面能分析 优化表面提交的能量计算结果如表6所示。表6 表明数据是收敛的， 根据表面能计算式 （1） ， 最终算 得表面能为1.309 7 J/m2。 3结论 （1） 在钠长石的晶体结构中， O原子为电子受体， 主要是2p轨道得到电子， 其余Na、 Al、 Si原子均为电 子供体， 所带电荷为正； 由于钠长石晶体结构的复杂 性， SiO键、 AlO键和NaO键都有多种不同的布 居和键长， 其中SiO键的布居值最大， 其共价性最 强， AlO 键次之， 而 NaO 键的离子性最强。这 3 种化学键中NaO键最易断裂， Na会溶于水中而使 钠长石表面荷负电， 易与胺类等阳离子捕收剂发生 静电吸附作用或分子吸附作用。 （2） 在钠长石 （010） 面费米能级附近的顶部价带 主要是由氧原子的2p轨道构成， 所以在浮选过程中 氧原子有更高的活性， 更易与浮选药剂发生作用， 是 其浮选的活性位点。 （3） 钠长石 （010） 面的表面层原子的外电子在达 到新平衡时出现了一定程度的电子重构现象， 表面 85 ChaoXing 层原子间的布居值也发生明显改变， 经过计算得构 建表面的表面能为1.309 7 J/m2。 参 考 文 献 谢贞付，王毓华，于福顺，等. 伟晶岩型锂辉石矿浮选研究综 述 ［J］ . 稀有金属， 2013， 37 （4） 641-649. Xie Zhenfu，Wang Yuhua，Yu Fushun，et al. Review on flotation of spodumene of pegmatite type ［J］ .Chinese Journal of Rare Metals， 2013， 37 （4） 641-649. 董栋， 程宏伟， 郭保万， 等. 锂辉石选矿技术现状及展望 ［J］ . 矿产保护与利用， 2018 （4） 130-134. Dong Dong，Chen Hongwei，Guo Baowan，et al. Status and pros- pect of spodumene mineral processing technology ［J］ . Conservation and Utilization of Mineral Resources， 2018 （4） 130-134. 谭鑫， 何发钰， 谢宇. 钨锰矿 （010） 表面电子结构及性质第 一性原理计算 ［J］ . 金属矿山， 2015 （6） 52-58. Tan Xin， He Fayu， Xie Yu. First-principles calculation of electron- ic structure and properties of tungsten-manganese（010）surface ［J］ . Metal Mine， 2015 （6） 52-58. 付博， 李龙江， 张覃， 等. 白云石电子结构及油酸钠在白云 石 （104） 面吸附作用的分子模拟研究 ［J］ . 矿冶工程，2019，39 （5） 69-73. Fu Bo，Li Longjiang，Zhang Qin，et al. Molecular simulation study on electronic structure of dolomite and adsorption of sodium oleate on the（104）surface of dolomite ［J］ . Mining and Metallurgical En- gineering， 2019， 39 （5） 69-73. Clark S J，Segall M D，Pickard C J，et al. First principles s using CASTEP ［J］ . Zeitschrift fuer Kristallograhie，2005，220 （5/ 6） 567-570. Sun Wei， Lu Zhongyuan， Tian Jia， et al. Anisotropic surface phys- icochemical properties of spodumene and albite crystals Implica- tions for flotation separation［J］ . Applied Surface Science，2017， 426 （12） 1005-1022. 徐龙华， 巫侯琴， 田佳， 等. 伟晶岩型铝硅酸盐矿物的晶体化 学特征计算与分析 ［J］ . 有色金属选矿部分， 2017 （6） 22-27. Xu Longhua，Wu Houqin，Tian Jia，et al. Calculation and analysis of crystal chemical characteristics of pegmatite-type aluminosilicate minerals ［J］ . Nonferrous MetalsMineral Processing Section，2017 （6） 22-27. 陈建华. 硫化矿物浮选固体物理研究 ［M］ . 长沙 中南大学出版 社， 2015. Chen Jianhua. Study on Solid Physics of Sulfide Mineral Flotation ［M］ . Changsha Central South University Press， 2015. 李玉琼， 陈建华， 陈晔， 等. 黄铁矿 （100） 表面性质的密度泛 函理论计算及其对浮选的影响 ［J］ . 中国有色金属学报， 2011， 21 （4） 919-926. Li Yuqiong，Chen Jianhua，Chen Ye，et al. Theoretical calculation of density functional theory of pyrite（100）surface properties and its effect on flotation ［J］ . Transactions of Nonferrous Metals Society of China， 2011， 21 （4） 919-926. 王杰， 张覃， 邱跃琴， 等. 方解石晶体结构及表面活性位点 第一性原理 ［J］ . 工程科学学报， 2017， 39 （4） 487-493. Wang Jie， Zhang Qin， Qiu Yueqin， et al. First principles of crystal structure and surface active sites of calcite ［J］ . Chinese Journal of Engineering， 2017， 39 （4） 487-493. 艾晶， 刘崇峻， 吴桂叶， 等. 阴阳离子捕收剂在长石浮选分离 中的研究进展 ［J］ . 中国矿业， 2019， 28（S） 314-316. Ai Jing，Liu Chongjun，Wu Guiye，et al. Research progress of an- ion and cation collector in feldspar flotation separation［J］ . China Mining Magazine， 2019， 28（S） 314-316. ［1］ ［2］ ［3］ ［4］ ［5］ ［6］ ［7］ ［8］ ［9］ ［10］ ［11］ 金属矿山2020年第6期总第528期 86 ChaoXing