第一性原理计算软件的使用
主输出文件
Hirshfeld Analysis(原子电荷分析)
Species Ion Hirshfeld Charge (e)
================================
O
1
-0.31
O
2
-0.31
O
3
-0.31
Ti
1
0.53
Ba
1
0.40
==================================
基态能量与电子结构的计算
在Electrtonic选项卡上,大部分的参数由前面的选 择自动给出,可以不改。具体的可以按“More”按钮 进行进一步设定。
电子结构计算的设定共有四个选项卡,分别为Basis, SCF,K-point和Potential。
精品课件
基态能量与电子结构的计算
电子结构计算设定的Basis选项卡 上,可以设定基组的能量截断值, 也即基组的大小。
在Potential选项卡上,可以选 择赝势参数。赝势的类型有二 种:Norm-Conserving和 Ultrasoft。势的表象可以在实 空间(real 空间),也可以在 倒易空间(Reciprocal space)。
精品课件
基态能量与电子结构的计算
Properties选项卡用来选取在第一性原理的计算中要计 算的物理量。每几物理量还有对应的参数可以选择。
点击“More”按钮后,出现几何 结构优化参数的具体设定对话框, 共有三个选项卡Minimizer, Options和Stress。
精品课件
晶体结构驰豫 在“Minimizer”选项卡上,可 以设定几何优化的收敛条件, 可以按等级自动给出,也可以 指定具体数值。
缺省的最大迭代数是100,可以 根据需要进行修改。
精品课件
点击“Density of States”之 后,出现态密度计算的选项,可 以选到计算PDOS。
点击“More”按钮后,可以进一 步设定K点的取值,或者指定交 换关联函数。
点击“Optical Properties”后, 可以计算光学性质,各种选项的 与态密度计算相同。
基态能量与电子结构的计算
可以设定FFT变换的网格的密度。 设定基组修正的方法。
在SCF选项卡上,可以设定SCF计 算的收敛判据,如能量变化量, 循环次数;
还可以设定电荷密度混合的参数 和方法。点击“More”按钮有更 多的选项。
精品课件
基态能量与电子结构的计算
在K-points选项卡上,可以设 定K点的选取参数,只计算 Gamma点,由Quality等级自动 给出,指定间隙,或直接指定 网格的数量。还可以设定原点 的漂移。
精品课件
点击“Band Structure”后, 下面出现能带计算的参数,如 空带数量,K点网格。
点击“More”按钮后出现更细 的选项,如指定交换关联函数, 和设定K点路径。
点击“Path”按钮,就可以得 到在布里渊区路径的对话框。 可以增加或减少能带结构的路 径。
基态能量与电子结构的计算
Material Studio中的 Visualizer提供了处理 大部分第一性原理计算 的方法。
右图为BaTiO3结构优化 完成后的界面。
精品课件
主输出文件
CASTEP模块计算的主输出文件是<名称>.castep,其中 的<名称>是原子结构模型的名称,记录了计算的全部过 程,包括出错信息。
1.84 4.90 0.00 0.00
Ti 11.04
1 0.96
2.35 6.58 2.12 0.00
Ba 8.74
1 1.26
2.07 5.98 0.69 0.00
精=品=课=件====================================== 在晶体中,因为没有明确的边界,每个原子的电荷数是
精品课件
Material Studio的界面
精品课件
主窗口中三种窗口:
原子结构窗口
K空间窗口,显示 了与原子结构对应 的K空间的单胞和 计算能带结构的路 径。
CASTEP计算的主 输出文件,列出了 计算参数,计算过 程,大部分的算结 果。
基态能量与电子结构的计算
下面以ZnO为例,显示基态能量计算的方法,以及得 到的电子结构和物理性质的结果。计算步骤如下: 1. 导入模型:从CASTEP提供的原子结构中导入ZnO的 原子结构文件,如图所示。
Species Ion s
p
d
Charge (e)
f Total
=========================================
O 6.74
1 -0.74
1.84 4.90 0.00 0.00
O 6.74
2 -0.74
1.84 4.90 0.00 0.00
O 6.74
3 -0.74
根据得到的能量变化和应力张量,调整晶体中原子的位 置和晶格常数。
当满足收敛条件时,结构优化结束,得到优化的晶格常 数,原子位置。
对于BaTiO3,晶格常数从4.01000变到4.00636 ,原子 在晶胞中的相对位置受到对称性的限制没有变化。
精品课件
主输出文件 Atomic Populations (Mulliken)
File菜单用来进行项目管理,包括新建,打开和保存项 目和窗口,导入、导出原子结构模型和其他文件。
View菜单用来调整界面上工具条和窗口的显示和隐藏, 可以根据需要进行调整。
Modify可以用来调整个原子结构,如修改原子,成键, 化学键杂化状态,调整电荷和电子态,设定原子运动的 约束条件。
Build菜单可用来构建原子结构模型,调整晶体结构的 对称性等。
可以设定是否要优化晶格参数。 在晶胞可变时有二种基组大小
设定的方法:固定基组大小或 固定基组品质。
在“Options”选项卡上,可以 选取算法,CASTEP中只有二种 算法可选:BFGS和DampedDM。
精品课件
晶体结构驰豫
精品课件
在“Stress”选项卡上,可 以设定晶体受到的应力,可 以计算晶体在高压下的性能。
“Run in parallel on”是指要用几个进程来进行计算。 双核超线程的计算机上可以选4个进程。
精品课件
自治计算收敛情况
精品课件
晶体结构驰豫
晶体结构豫是第一性原理计算中 经常需要进行的计算,目的在于 找到能量最低的基态。
要进行晶体结构的驰豫,在 “CASTEP Calculation”对话框 的“Setup”选项卡上的“Task” 下拉框中选取“Geometry Optimization”。
除了主窗口外的三个窗口也可以通过View菜单来显示和 隐藏。
项目窗口管理第一性原理计算的所有文件,包括原子模 型,计算条件设定,计算结果等。
属性窗口可用来显示和修改当前窗口中对象的属性。 作业调度窗口显示正在进行的计算作业的状态,可以用 精品来课件中断提交的计算。
Material Studio的界面
属性菜单显激活窗口中选中元素的属性。可以显示属性 值,也可以进行修改
可以用Window菜单对主窗口进行操作,如各种排列方式, 图标化显示,关闭窗口等。
在原子结构窗口,按住鼠标左键用于选取原子结构图中 的单元,如原子,键等,按住鼠标右键拖动可以对模型 进行旋转,以便更好地观察。
可以用测量工具测量原子结构中的原子间距,夹角,扭 角等。
偏振光需要给出偏振方向,非偏振 光需要给出入射方向。
介电常数
精品课件
反射率
精品课件
吸收系数
精品课件
折射率
精品课件
光导率
精品课件
Loss Function
精品课件
热力学性质-德拜温度
精品课件
作业
选取一种金属材料(如铜,铝等),一种绝缘体材料, 分别用利用第一性原理计算其DOS分布图和能带结构图, 说明金属和绝缘体的区别。
<名称>.castep首先列出了用于计算的所有参数,包括 人工设定的和系统默认的,然后列出了体系的晶格参数, 原子种类,采用的赝势,原子坐标等。
列出计算中用到的K点,然后是对称性分析有约束条件。 通过三个不同的截断能量计算得到的基态能量值,估算
基组引起的修正值。
精品课件
主输出文件
开始结构优化的大循环,每个循环步聚中都要计算电子 结构的自治计算,得到基态能量和应力张量。调整晶体 中原子的位置和晶格常数。
精品课件
Total DOS
精品课件
能带结构
精品课件
可以得到能带结构图。带隙 1.88eV。直接带隙。
比实验值小,是DFT的固有缺 点。
光学性质
精品课件Байду номын сангаас
通过第一性原理计算得到电子结构 后,可以计算物质的光学性质,包 括介电常数,反射率,吸收率,折 射率,光导率,和Loss Function。
可以计算偏振光,非偏振光,和多 晶体的光学性质。
调整Display Style对话框中的isovalue选项卡中的isovalue值,就可以看到 等高面的变化,也可以看剖面上的密度分布,详见实际演示。
精品课件
DOS/PDOS分布图
Ti原子
精品课件
DOS/PDOS分布图
Ba原子
精品课件
DOS/PDOS分布图
氧原子
精品课件
DOS/PDOS分布图
Tool菜单提供了一些常用的小工具,帮助处理计算中的 一些问题。
Module列出了Material Studio中所有的模块,以及对 精品他课件们的操作。
Material Studio的界面
主窗口是一个多文档窗口,可以容纳其他各种窗口,如 原子结构窗口,各种图表窗口,文本窗口等。