快捷下单入口 关于 合作 招聘 新人手册 会员中心

    热线:400-152-6858

    测试狗科研服务

    预存 免费试测 登录
    Document
    当前位置:文库百科文章详情
    物理VS化学吸附?RDG让你的文章与众不同
    来源:测试GO 时间:2023-02-13 09:52:25 浏览:8165次


    引言

    约化密度梯度(reduced density gradient, RDG)是密度泛函理论中用来描述电子非均匀性的无量纲参量[注1],其表达式为图片


    2010年,杨伟涛教授提出基于RDG的非共价相互作用(noncovalent interactions, NCI)[1]分析方法,可以用来归属原子间或分子间的相互作用。虽然电子定域化函数(electron localization function,ELF)和分子中的原子(atoms-in-molecules, AIM)理论等能够用来分析化学键,但是对氢键、空间位阻、π-π堆叠等非共价相互作用的讨论有一定的局限性,而NCI分析则可以作为揭示这些非共价相互作用强有力的工具。因为分子的电子密度呈现指数衰减,所以在远离分子处RDG呈现很大值;相反,在相互作用的区域RDG的值则很小,其中非共价相互作用区域的电子密度和RDG都比较低。

    由于氢键和空间位阻区域的电子密度和RDG很接近而难以进行区分,所以需要再借助通过电子密度Hessian矩阵的第二大本征值的正负号(sign(λ2))来进一步判断[注2];电子密度大小则可以反映出相互作用的强度,一般认为vdW作用的区域ρ<0.005,因此可以通过对RDG等值面进行sign(λ2)ρ填色或者做RDG vs sign(λ2)ρ的散点图来讨论所研究体系的相互作用类型。

    在Sobereva的博文中已经介绍了RDG相关的概念和使用Multiwfn对非周期边界系进行NCI分析的详细过程[2]。但固体和表面以及牵扯过渡金属的体系的NCI讨论较少。因此,本文将通过若干实例的计算和相应的分析,对这些体系使用Quantum ESPRESSO(QE)计算RDG的过程和NCI分析方法进行介绍。


    准备工作

    本文所采用的计算程序如下:

    1. Quantum ESPRESSO 6.7 [3] 

    实际上QE在5.1.0版本就已经支持NCI分析,所以即使正在使用旧版本的用户也不必担心无法完成本文的工作。

    2. VESTA 3.1.8 [4] 

    个别新版本有bug,所以本文采用旧版本。

    测试体系为:

    1. 来源于Materials Project [5]的结构

    MoS2(mp-2815)和Ni(OH)2(mp-27912)

    2. H2O2@PMCS来源于文献[6]


    计算参数

    计算参数选择的依据:从文献[7]测试结果来看,vdW-DF3-opt2在S22×5和S66×8的表现都还不错,原文中所使用的是GBRV超软赝势[8],该赝势本身也有很好的可靠性和可移植性[9,10]。建议密度截断(ecutrho)选取充分以保证计算结果可靠以及避免RDG计算过程产生噪声。


    计算过程和数据处理

    A1 使用PW模块

    进行结构优化,并进行自洽计算产生波函数以及电子密度

    A2 通过PW/tools/pwo2xsf将out文件转换为xsf格式[注3],如

    A3 使用PP模块进行后处理,通过plot_num=19产生RDG格点,以及 output_format=5来产生xsf文件

    A4 类似过程A3,使用PP模块通过plot_num=20产生sign(λ2)ρ的xsf文件,但建议提前将文件夹拷贝一份或者通过fileout更改xsf输出格点文件的名称避免上一步的结果被覆盖。

    B1 打开VESTA程序,将A2中out文件产生的xsf文件拖入VESTA窗体

    B2菜单栏中Edit->Edit Data->Volumetric Data,分别点击Isosurface和Surface coloring旁的Import分别导入RDG网格和sign(λ2)ρ的xsf文件。

    B3 菜单栏中Objects->Properties->Isosurfaces,将Isosurfaces中的Isosurface level改为0.5,以及将Surface coloring中的Saturation level的Max改为0.02,Min改为-0.04,这样能够和文献[1]较好符合。

    B4 调整其他显示部分,比如菜单中去除勾选Objects-> Volumetric Data->Show Sections,以及在Objects->Properties->Atoms对原子颜色进行调整,在Edit->Bonds对显示成键的判断进行调整,使得作图更加美观,但具体调整的细节这里不进行赘述。如果对操作过程有问题也欢迎大家讨论。

    案例

    MoS2是层状结构材料,MoS2层间出现了RDG等值面围成的区域,且层间S……S连线的区域呈现淡绿色(ρ<0.005),证明了层间存在相互作用,且通过vdW作用相结合并主要是以层间S……S的相互作用为主导;而层内的三个Mo原子之间显示为红色,表现出位阻效应。图片

    Ni(OH)2与MoS2情况类似,但Ni(OH)2层间作用以O-H……O-H的氢键作用为主导;层内的Ni原子之间也呈现出位阻效应。在Ni-O键中形成的蓝色环形区域是因为Ni-O之间相互作用比较强,而pp.x程序在处理RDG的时候将ρ>0.05的网格直接设置为100[注4],导致共价键或离子键在图中不被显示。

    H2O2@PMCS这一类结构常见于单原子催化过程中,也是我们计算较多的一类结构。文献[6]给出的吸附能为-0.45 eV,仅从吸附能来看可能是物理吸附。尽管本文采用的方法与文献不同,但给出了-0.47 eV的结果与文献[6]十分接近。

    为了进一步确定H2O2与基底相互作用的类型,我们进行了NCI分析,结果表明过氧化氢与Zn-N4的相互作用却是是非共价相互作用,其中H2O2中一个O-H与Zn-N4中的N原子形成了氢键;而另一个O-H的O与Zn原子产生了静电作用,这种静电作用从颜色上看要更强于氢键作用,因此可以认为H2O2在PMCS的吸附作用是由这种静电作用所主导的物理吸附。


    Promolecule近似

    Promolecule近似就是将球对称的自由原子密度直接进行叠加,即

    Promolecule近似下的密度通常可以作为自洽场的初始猜测,也能够用于大体系的定性计算。在QE中也能实现Promolecule近似的计算,这是因为QE中UPF格式的赝势文件中已经存有自由原子的电子密度。在自洽场计算过程中,我们只进行初始猜测而不做自洽过程,因此将electron_maxstep设置为0,且startingwfc设置为’atomic’,就可以得到Promolecule近似下的密度。

    由于Promolecule近似做NCI分析也能给出定性可靠的结果[1],因此在处理更大的体系时可以使用这种方法。这里仍然以H2O2@PMCS体系为例,展示Promolecule近似下NCI分析的结果,下图表明Promolecule近似(等值面取0.035)与自洽计算后的差距并不大。


    散点图

    由于非共价相互作用区域电子密度和RDG都比较低,且λ2的符号为负,通常在散点图中呈现为一个或多个尖峰(spike)。如H2O2@PMCS的散点图中,左下角横坐标位于-0.04~0.02的spike就对应于H2O2中OH与N的氢键。

    而promolecule近似下的散点图与上图的结果较为接近

    那么如何绘制这样的散点图呢?我在这里提供了一个自己写的脚本nci_scatter.py能够直接绘制散点图。需要预先安装python 3.x,调用到的python组件有numpy和matplotlib,只要通过pip install就能够进行安装。只要执行nci_scatter.py脚本按照每一步的提示用鼠标拖入对应的xsf文件就能够直接画出上面的散点图。

    “nci_scatter.py”

    链接:

    https://pan.baidu.com/s/1rtUFsfvlFPtLiqROc2qNHw

    提取码:csgo

    请大家请自行下载~


    总结

    Quantum ESPRESSO是一款功能十分强大且开源免费的第一性原理计算程序,NCI分析也是一种强大的非共价相互作用分析方法。QE已经内置了RDG的计算,但却很少有人写QE做NCI分析的教程,本文以固体和表面中的一些常见体系为案例,对实现过程和分析进行了一些叙述讨论,更多的理论基础和分析细节还是建议大家多研读文献。希望本文能够为大家的科研工作助力添彩,同时也欢迎大家共同参与讨论互相进步。


    注释

    [2] 这些概念来源于AIM理论,电子密度的Hessian矩阵是一个3×3的矩阵,对角化后得到三个本征值λ1>λ2>λ3。若λ2<0意味着λ3也小于0,对应于键临界点(bond critical point,BCP);若λ2>0意味着λ1也大于0,对应于环临界点(ring critical point, RCP)。

    [3] pwo2xsf的作用是将pw计算得到的out文件转化为xsf文件,可选选项有

    -ic 初始结构

    -lc 最终结构

    -oc 最优结构

    -a 将轨迹转化为AXSF

    [4] 严谨起见,我们检查了QE的pp.x程序的源代码,发现在QE 6.x版本分析磁性体系的sign(λ2)ρ和ELF时,电子密度存在bug(这是因为自旋极化的计算中,QE 6.x版本的电子密度不是按照spin up和spin down进行存储的,而是按照total和spin进行存储的,但sign(λ2)ρ和ELF这两段代码忘记进行修改)。这点我们已经向QE官方提交了这个bug并得到了确认,所以先建议大家不要分析磁性体系,后续版本会进行改进。本文中的案例Ni(OH)2的结果展示是我们使用修改后的程序,因此结果可靠。


    参考资料

    [1] Johnson E R, Keinan S, Mori-Sánchez P, et al. Revealing noncovalent interactions[J]. Journal of the American Chemical Society, 2010, 132(18): 6498-6506.

    DOI: 10.1021/ja100936w

    [2] http://sobereva.com/68

    [3] https://github.com/QEF/q-e/releases

    [4] http://www.jp-minerals.org/vesta/en/download.html

    [5] https://www.materialsproject.org

    [6] Xu B, Wang H, Wang W, et al. A Single‐Atom Nanozyme for Wound Disinfection Applications[J]. Angewandte Chemie, 2019, 131(15): 4965-4970.

    DOI: 10.1002/anie.201813994

    [7] Chakraborty D, Berland K, Thonhauser T. Next-Generation Nonlocal van der Waals Density Functional[J]. Journal of Chemical Theory and Computation, 2020, 16(9): 5893-5911.

    DOI: 10.1021/acs.jctc.0c00471

    [8] Garrity K F, Bennett J W, Rabe K M, et al. Pseudopotentials for high-throughput DFT calculations[J]. Computational Materials Science, 2014, 81: 446-452.

    DOI: 10.1016/j.commatsci.2013.08.053

    [9] Marsman M, Marzari N, Nitzsche U, et al. Reproducibility in density functional theory calculations of solids[J]. Science, 2016, 351: 6280.

    DOI: 10.1126/science.aad3000

    [10] https://molmod.ugent.be/deltacodesdft

    评论 / 文明上网理性发言
    12条评论
    全部评论 / 我的评论
    最热 /  最新
    全部 3小时前 四川
    文字是人类用符号记录表达信息以传之久远的方式和工具。现代文字大多是记录语言的工具。人类往往先有口头的语言后产生书面文字,很多小语种,有语言但没有文字。文字的不同体现了国家和民族的书面表达的方式和思维不同。文字使人类进入有历史记录的文明社会。
    点赞12
    回复
    全部
    查看更多评论
    相关文章

    【科研干货】电化学表征:循环伏安法详解(上)

    2019-10-25

    晶体结构可视化软件 VESTA使用教程(下篇)

    2021-01-22

    手把手教你用ChemDraw 画化学结构式:基础篇

    2021-06-19

    【科研干货】电化学表征:循环伏安法详解(下)

    2019-10-25

    基础理论丨一文了解XPS(概念、定性定量分析、分析方法、谱线结构)

    2020-05-03

    Zeta电位的基本理论、测试方法和应用

    2020-08-24

    热门文章/popular

    【科研干货】电化学表征:循环伏安法详解(上)

    晶体结构可视化软件 VESTA使用教程(下篇)

    手把手教你用ChemDraw 画化学结构式:基础篇

    电化学实验基础之电化学工作站篇 (二)三电极和两电极体系的搭建 和测试

    【科研干货】电化学表征:循环伏安法详解(下)

    基础理论丨一文了解XPS(概念、定性定量分析、分析方法、谱线结构)

    微信扫码分享文章

    意见反馈

    有奖举报

    商务合作

    ...

    更多

    公众号

    关注我们 了解更多

    小程序

    随时预约 掌握进度

    举报有奖

    TEL: 191-3608-6524

    如:在网络上恶意使用“测试狗”等相关关键词误导用户点击、恶意盗用测试狗商标、冒称官方工作人员等情形,请您向我们举报,经查实后,我们将给予您奖励。

    举报内容:

    200

    上传附件:
    文件格式不正确,请重新上传文件格式不正确,请重新上传文件格式不正确,请重新上传
    文件格式:jpg、jpeg、png、gif、tif、doc、docx、ppt、pptx、xls、xlsx、pdf、zip、rar
    联系方式
    姓名
    电话
    提交意见

    意见反馈

    Suggestions

    您可以在此留下您宝贵的意见,您的意见或问题反馈将会成为我们不断改进的动力。

    意见类型
    测试服务
    网站功能
    财务报账
    其他类型
    意见内容

    200

    联系方式
    姓名
    电话
    提交意见

    收起

    01

    专属信用额度,先测后付0元下单

    02

    下单享高额积分,万千好礼免费兑

    200
    200元无门槛优惠券
    立即激活 立即下单
    已使用
    已作废

    全流程在线可视化,便捷高效触手可及

    如下单过程中有任何疑问或需要帮助,请随时咨询专属顾问~
    9
    9折无门槛优惠券

    支付一笔订单后可领取

    立即领取 立即下单
    已使用
    已作废

    免费测+惊喜盲盒+高额福利,多重福利大放送

    新人免费测

    双双开盲盒(100%中奖)

    邀请人专享

    受邀人专享

    创建/加入团队,解锁定制化权益

    01

    1500元团队专属优惠券

    02

    万元大额信用额度,享先测后付

    03

    团队成员统一开票报销;

    04

    专业工程师课题专属服务

    领取成功,请下单
    请您支付一笔订单后才可以领取优惠券
    Document
    关于我们 新手帮助 测试干货 商务合作 基金查询 相关资质 模拟计算 现场测试 服务项目 科研绘图 同步辐射 电池行业

    联系方式/contact

    400-152-6858

    工作时间/work time

    09:00-18:00

    测试狗公众号

    关注我们 了解更多

    测试狗小程序

    随时预约 掌握进度

    蜀公网安备51010602000648号

    蜀ICP备17005822号-1

    成都世纪美扬科技有限公司

    Copyright@测试狗·科研服务