Materials-Studio软件CASTEP模块知识分享.pptxVIP

Materials Studio软件CASTEP模块——架构原理及主要参数意义主要内容 1 矿物晶体结构搭建 2 CASTEP整体架构、其中假设 3 CASTEP参数设置及其意义矿物晶体结构搭建（1）从MS数据库中引入结构 MS数据库中为用户提供了较为常见的结构模型：Catalysts 催化剂Organics 有机物Ceramics 陶瓷Polymers 聚合物Glasses玻璃Repeat-units 重复单元Metal-oxides 金属氧化物 Semiconductors 半导体Metals金属Zeolites沸石族Minerals 矿物Molecular-crystals 分子晶体Nanotubes 纳米管矿物晶体结构搭建（2）自助搭建结构 对于结构库中没有的晶体结构，要从文献、XRD软件标准卡片等资料中查找参数自行搭建。如：AlAs晶体结构 空间群 F-43m（代号216） 晶格参数 a=b=c=5.6622? α=β=γ=90o 原子占位 Al（0 0 0） As（0.25 0.25 0.25）矿物晶体结构搭建（3）外部数据库导入结构无机晶体学数据库ICSDAmerican Mineralogist Crystal Structure Database可输入矿物名称、作者、化学元素、晶胞参数及对称性等关键词查找，可直接获得.cif结构文件和XRD检测数据，还可以预览3D模型。矿物晶体结构搭建注意：引入结构或网上查找参数及结构文件时，注意识别需要的结构构型或数据，有的结构是人工合成的或者经过高温、高压处理的，不能采用，如MS库中的ZnOCASTEP整体架构、其中假设CASTEP的整体构架就是基于密度泛函理论求解K-S方程的一种方法决定微观的电子结构基本物理性质获取求解多粒子系统的薛定谔方程密度泛函近似求解多粒子系统薛定谔方程密度泛函理论内容Kohn-Sham方程得到简化哈密顿算符定理1：多电子体系的基态总能量 是电子密度的唯一函数定理2：电子密度确定，体系的基 态性质就唯一确定多电子运动状态→电子密度分布状态K-S方程自洽迭代求解过程nstart(r) 生成KS势Density Mixing:推荐使用All Bands/EDFT:速度慢，占内存。当使用DM计算金属体系SCF不收敛时推荐使用求解KS方程得到新nout(r)和之前的n(r)比较收敛与否 ?输出结果K-S求解近似处理方法-交换相关泛函电-电库伦作用能电 动， 常数一次迭代确定电-电交换相关能核-电赝势处理局域密度近似 LDA广义梯度近似 GGA交换相关泛函-局域密度近似LDA假设原子核外电子云均匀分布LDA适用的情况：（1）电荷密度变化缓慢的体系（如金属）（2）电荷密度较高的体系（如过渡金属）（3）适用于大多数晶体结构（对晶胞参数描述准确）真实情况下原子核外电子并非均匀分布LDA不适用的情况：（1）电子分布定域性较强，电荷密度分布不均匀（化学反应中的过度态）（2）体系束缚能绝对值估计不准确（3）低估禁带宽度的绝对值（固有缺陷）交换相关泛函-广义梯度近似GGA GGA克服了LDA在描述真实体系在电子密度变化剧烈的情况下的缺陷，提高了交换相关能计算结果的精度，也提高了密度泛函方法计算的精度。非定域泛函：HF HF-LDA sX sX-LDA PBE0 B3LYP……一般研究半导体材料采用，精度有所提高，但运算量是LDA和GGA的几十到几百倍，在计算精度与计算时间二者中做出取舍，一般没有必要采用交换相关泛函的选择 虽然GGA方法弥补了LDA的缺陷，提高了计算精度，但也很难说GGA一定优于LDA，不能说对于某一体系哪一种方法一定适用。 实际使用过程中交换相关泛函的选择方法有两种： （1）阅读文献，根据他人已发表的成功经验直接选择，这种方法最省时省力 （2）以交换相关泛函为单一变量做收敛性测试，与试验结果比较，权衡计算精度与计算成本择优选用K-S求解近似处理方法-赝势电-电库伦作用能电 动， 常数一次迭代确定电-电交换相关能核-电赝势处理局域密度近似 LDA广义梯度近似 GGApseudo wave functionpseudopotentialK-S求解近似处理方法-赝势CASTEP将外层电子波函数通过平面波函数展开，将原子实近似为新的“核”有效的减少平面波数目赝势其合理性在于，采用赝势前后：能量本征值不发生变化价电子波函数在Rc外的分布不变在Rc处的电子波函数对数的倒数不变化只适用于周期性体系RcK-S求解近似处理方法-赝势赝势设置Ultrasoft (USP)超软推荐用，精度好，效率高Norm-conserving (NCP) 模守恒 拉曼光谱计算不支持USP，选用NCP为了与已发表文献比较