【入门基础教程】单点计算与BAND的基本参数设置、计算效率
ADFinput窗口,用于建模、参数设置等。File - Import Coordinates,可以导入cif文件,也可以在ADFinput右上角的搜索按钮搜索一些简单常见的晶体结构。
说明:
单点计算一般只需要简单设置Main窗口即可
Task:Single Point,单点的含义是指计算的体系几何结构确定、不变,不做优化,不扫描势能面
Followed by:指计算完毕之后,是否需要计算其他
Periodicity:Bulk表示体系沿着ABC三个方向无限重复;Slab则是二维体系,单胞也是二维的,在AB方向无限重复,上下为半无限大的真空;Chain是一维材料,一维材料(例如聚合物)周围是无限大真空;None表示体系没有周期性
Total Charge:体系单包的带电量,2019版不允许周期性体系带电,2020版以后将会允许带电
Spin Polarization:单包内有多少个电子没有配对?如果此处设置为2,则表示单包是三重态;如果勾选Unrestricted,此处留白,则计算结果会自动收敛到能量最低的自旋多重度,例如O2分子,将得到三重态的结果。但是一般只建议在磁性材料基态计算时这样做,作为基态的初始猜测,然后降低、升高自旋极化,去确认是否基态正确。还有一种情况也可以用到这种功能,即,存在自旋翻转的分子动力学模拟,例如有O2参与化学反应的分子动力学模拟,其自旋是随时可能变化的,因此不设置更准确。
如果任何情况下,都盲目地不设置自旋极化,同时勾选Unrestricted,一方面导致每一步SCF迭代的计算量至少增大一倍,同时有可能导致收敛困难,最后浪费自己的大量的时间,实际效率下降上百倍。
Unrestricted:如果勾选,则α电子与β电子分别求解电子态,
每一步SCF迭代计算量都会增大一倍以上;如果不勾选,则β电子的能级、波函数沿用α电子的结果。因此,如果是
闭壳层体系,请尽量不要勾选。只有对开壳层体系做pEDA-NOCV分析时,才可能对闭壳层片段也使用Unrestricted(但此时,也请务必设定正确Spin Polarization值)。
XC functional:选择泛函
Relativity:是否考虑相对论效应,默认不考虑。如果考虑相对论对动能的修正,则选择Scalar,如果选择Spin-Orbit则不仅考虑相对论对动能的修正,还考虑相对论效应带来的自旋轨道耦合。没有重元素的情况下,相对论实际上可以设置为None,这样计算量会降低约60%。Scalar一般用于金、银、汞、镧系、锕系这种级别的重元素。
Basis Set:选择基组。基本原则是更重的元素选择更大的基组,具体参考:
BAND:如何设置基组。对于吸附问题的计算,包括氢键、范德华力的研究,都需要增大基组至少一个级别。对NMR、ESR研究,往往使用最大基组。
Frozen Core:该设置是为了减小计算量,内层电子不参与自洽迭代,直接沿用真空中单原子的结果,冻芯越大,节省的计算量越大,具体冻结到什么级别,可参考:
BAND:如何设置基组
-
Calculate:如果要计算DOS、Band Structure就勾选一下对应的选项即可
Details → SCF → Cryterion:这是SCF收敛阈阈值,默认10-6这个标准是很严格的,如果体系较大,这个总能量就会更大,这个标准可以按比例降低标准,例如有100原子的体系,这个标准就设置为10-5就可以了,更大的体系,甚至可以降低标准到10-4
计算完毕之后,体系的电子云密度空间分布,可以通过SCM - View - Add isosurface - Select Field - Density - Total density查看。
注意
BAND的能级是以无限远处真空能级为参考0点的,所以得到的费米能级(一般小于0)、能带能量某种程度上说,是绝对量,不是相对量。