当前位置:
文档之家› 第六讲第一原理计算方法简介及MaterialsStudio中Castep使用方案
第六讲第一原理计算方法简介及MaterialsStudio中Castep使用方案
Visualizer:图形化建模模块
可构建计算的模型:晶胞,分子,晶体表面, 纳米结构,聚合物等
锐钛矿TiO2
TiO2(111)
Pt(110)-CO(2x1)
碳纳米管
TiO2纳米棒
Material Studio 晶体结构模型建立
建立方铅矿PbS晶体结构模型(实验15) 步骤 (1)查晶体结构数据(如ICSD、PDF数据
在CASTAP计算中有很多运行步骤,可分为如下几组:
结构定义:必须规定包含所感兴趣结构的周期性的3D模型 文件,有大量方法规定一种结构:可使用构建晶体 (Build Crystal)或构建真空板(Build Vacuum Stab)来 构建,也可从已经存在的结构文档中引入,还可修正已存 在的结构。
注意: CASTEP仅能在3D周期模型文件基础上进行计算, 必须构建超单胞,以便研究分子体系。
晶体、结晶与X射线 衍射
MS.Polymorph Predictor MS.Morphology MS.X-Cell MS.Reflex MS.Reflex Plus MS.Reflex QPA
分子力学与分子动力学 MS.DISCOVER PASS MS.Amorphous Cell MS.Forcite MS.Forcite Plus MS.GULP MS.Equilibria MS.Sorption
Hartree Fock方程
薛定谔方程简化为:
将总Hamilton分解成单电子贡献H0和电子-电子相互作 用U。应用变分法计算多电子波函数方程,可得HartreeFock方程。
量子化学分子轨道方法
分子轨道方法:在Hartree-Fock框架下,将单电子波函数 用原子轨道(Slater型-STO,Gaussian型-GTO)的 线性叠加表示来求解。
√
√
√
√
Reflex-Powder Indexing
√√
Reflex-Powder Refinement
√√
Reflex Plus
√
√ √ √ Reflex QPA
√
√ √ √ Sorption
√ √ √ Synthia
√ √ √ VAMP
√√
X-Cell
√
√ √ √ Mesotek
√
√ √ √ Morphology
密度泛函理论
赝势(pseudo potential) 赝势就是把离子实的 内部势能用假想的势能 取代真实的势能,但在 求解波动方程时,不改 变能量本征值和离子实 之间区域的波函数。模 守恒赝势NCP (Norm Conserving Pseudopotential) 和 超软赝势 USPP(Ultrasoft Pseudoptential)
精 组态相互作用方法(采用多个Slater行列式考虑电子关联)
度 Mφller-Plesset(MP)修正(将关联作用作为微扰修正)
, 计
Hartree-Fock方法(忽略交换作用,严格计算电子积分)
算 半经验方法,如CNDO,MNDO,MINDO,AM1,PM3
量 等(同样忽略交换作用,近似计算电子积分)
理常数,如光速c、Planck常数h、电子电 量质量e、,电因子此质处量理m不e以同及体原系子时的候各具种有同较位好素的的可
移植性(transferability)。但是,在具 体实行时,仍依赖于具体近似方法的选取, 从而带来系统误差。
多粒子体系(电子+核)的薛定谔方程
三个近似
a. 非相对论近似(忽略了电子运动的相对论效应) ve<<c,ve~108cm/s<3×1010cm/s),me=m0 求解非相对论的薛定谔方程,而不是相对论的狄拉克方程
Module
parallel Windows Linux Linux IA32 IA64
Module
parallel Windows Linux Linux IA32 IA64
Materials Visualizer Adsorption Locator
Amorphous Cell
Blends
√
ONETEP
密度泛函理论
基组(basis set) 求解Kohn-Sham方程,选取适当的基组, 将波函数对其展开,将方程求解转化为线 性代数问题。 一般选用如下基组展开:
(Linearized) augmented plane waves (L)APW’s
(Linearized) muffin-tin orbitals - (L)MTO’s Projector augmented waves -PAW’s
第一原理计算软件
Code Basis Name Set
ABINIT Plane wave
CASTEP Plane wave
PWscf Plane wave
VASP Plane wave
LAPW WIEN2K
Potentials Plane Wave Pseudopotential Codes
Pseudo, PAW
密度函数
电子与原子核间的库仑势 电子间的库仑势 交换关联势 (未知)
密度泛函理论
LDA和GGA近似 Kohn-Sham方程原则是精确的,但遗憾的 是交换关联势是未知的。要进行具体计算, 就必须使用近似方法求出交换关联势。常 用的近似方法有局域密度近似(Local Density Approximation)和广义梯度近似 (Generalized Gradient Approximation), 在某些情况下,广义梯度近似改善了局域 密度近似的计算结果,但它并不总是优于 局域密度近似。
√
√
√√
√√
Polymorph
√ √ √ QMERA
√√
QSAR and QSAR Plus
√
√
√
√
√
√
CASTEP and NMR CASTEP
Reflex-Pattern
√
√ √ √ Processing and
Powder Diffraction
√
COMPASS
CCDC
Conformers Discover DMol3 DPD Equilibria Forcite Gaussian GULP MesoDyn
第一原理常用计算软件
根据对势函数及内层电子的处理方法不同 主要分为两大类,一种是波函数中包含了 高能态和内层电子,而势函数只是原子核 的贡献,这称为全电子(all electron calculation)法,另一种处理方法是势函 数为原子核和内层电子联合产生的势,称 为离子赝势,波函数只是高能态电子的函 数,这称为赝势(pseudo-potential)法。
Pseudo
Pseudo
Pseudo, PAW all-electron
操作系统
Linux
Web Site
www.abinit. org
Windows Linux
Linux
www.tcm.ph / castep/
www.pwscf.o rg/
Linux Linux
cms.mpi.un ivie.ac.at/v asp
密度泛函理论
Hartree-Fock方法的主要缺限:(1)完全忽略电子 关联效应;(2)计算量偏大,随系统尺度4次方关系 增长。
20世纪60年代,Hohenberg,Kohn和Sham(沈吕九) 提出了密度泛函理论(DFT)。DFT理论奠定了将多电 子问题转化为单电子方程的理论基础,给出了单电 子有效势计算的可行方法,DFT在计算物理、计算化 学、计算材料学等领域取得巨大成功。1998年,W. Kohn与分子轨道方法的奠基人Pople分享了诺贝尔化 学奖。
库) (2)通过软件建模(如Material Studio中
模块 Visualizer 、Diamond)
Castep使用
CASTEP模块 Cambridge Serial Total Energy Package)
CASTEP是特别为固体材料学而设计的一个现代的量子力学基本程序,其 使用了密度泛函(DFT)平面波赝势方法,进行第一原理量子力学计算,以 探索如半导体,陶瓷,金属,矿物和沸石等材料的晶体和表面性质。
b. Born-Oppenheimer近似,核固定近似 中子/质子的质量是电子质量的约1835倍,即电子的运 动速率比核的运动速率要高3个数量级,因此可以实现 电子运动方程和核运动方程的近似脱耦。这样,电子可 以看作是在一组准静态原子核的平均势场下运动。
c.单电子近似 把体系中的电子运动看成是每个电子在其余电子的平均 势场作用中运动,从而把多电子的薛定谔方程简化单电 子方程。
提示: CASTAP计算所需时间随原子数平方的增加而增加。 因此,建议是用最小的初晶胞来描述体系,可使用 Build\Symmetry\Primitive Cell菜单选项来转换成初晶
CASTEP的任务
计算设置:合适的3D模型文件一旦确定,必须选择计算类型 和相关参数,例如,对于动力学计算必须确定系综和参数, 包括温度,时间步长和步数。选择运行计算的磁盘并开始 CASTEP作业。 结果分析:计算完成后,相关的CASTEP作业的文档返回用户, 在项目面板适当位置显示。这些文档进一步处理能获得所需 的观察量如光学性质。
5. 可模拟的内容:催化剂、聚合物、固体化学、 结晶学、晶粉衍射以及材料特性等。
主要模块:
建模模块
VisualiBlends CASTEP Conformers DMol3 DPD Discover Equilibria Forcite
GULP MesoDyn Morphology Onetep Polymorph QMERA Reflex Synthia VAMP Gaussian