Dmol3、Castep的基本原理和参数设置.docxVIP

Dmol3、Castep的基本原理和参数设置

materialsstudio

DMol3基本原理和参数设置

materialsstudio

结构制备材料

H=ESCF(DFT)

电子波函数能量一阶导数二阶导数

能带，态密度，电荷密度

力场

应力，弹性力常数，体模量….声子频率，散射谱

表征材料

制备工艺

各种性能参数

怎么使SCF收敛？如何使SCF的结果准确？

materialsstudio

**********NAL*****ohenberg-Kohntheorem

E[(rl,...rN)]E[(r)]EToUExc

WalterKohn

Kohn-Shamtheorem

Kohn-Shamequations

2

原子轨道。DND一双数值轨道基组+d轨道极化函数。与DN基组类似,但是为所有非氢原子加入了d轨道函数极化。为DMol3模块默认基组，在可接受的计算时内，确保精度在可接受的范围内。也是DMol3计算的最低可信精度基组。DNP一双数值轨道基组+轨道极化函数。与DND基组类似，但是为所有氢原子加入了p轨道函数极化。精度更高，计算时间更长。对于氢键计算更为重要。

TNP一三数值轨道基组+轨道极化函数。与DNP基组类似，为所有原子加入极化函数。目前只能够使用与H到CI的元素，且不包括He和Neo精度更高，计算时间更长。

veff[(r)]i,k(r)i,ki,k(r)

ExactonlyforgroundstateNeedsapproximationtoExc

materialsstudio

基于DFT的自洽计算过程nstart(r)生成KS势求解KS方程rveff[n(r

)]

2

veff[n(rr)]i,kr(rr)i,kri,kr(rr)

得到新nout(r)和之前的n(r)比较收敛与否？

rrrnrfiiriri

输出结果

materialsstudio

DFT和磁学性能电荷密度依据电子自旋方向的不同，一分为二：自

旋密度定义为：

总的磁矩定义为：

rrrn(r)n(r)n(r)rrr(r)n(r)n(r)rr3Mn(r)n(r)dr

存在自旋极化的体系中，交换相关势也相应的变为不同自旋方向电荷密度的函数：

rrvxcvxc[n(r),n(r)]

PS:

由于在自洽计算中需要分别考虑电荷密度和自旋密度的收敛，因此计算时间会延长，收敛难度会增加。

materialsstudio

Dmol3的基本原理和参数设置

materialsstudio

DM0L3:原子轨道线形组合法(LCAO)

j(r)Rnl(r)Yim(,)ImRadialportionatomicDFTeqs.numericallyAngularPortion

Rcut周期性和非周期性体系

适合于份子、团簇、份子筛、份子晶体、聚合物等〃开放类结构〃Tips:对于空体积较大的晶体，使用DMol3的效率要高于Castep

materialsstudio

在DMol3模块中，电子密度实际上由各个原子轨道的平方和来确定：

在这里，电子密度实际上是由所有占领的份子轨道＜pi.来决定。份子轨道可能由上自旋电子(Alpha电子)和下自旋电子(Beta电子)占领。当Alpha电子和Beta电子的数目相等的时候，我们可以用单一的份子轨道epi.来进行表述，这种体系称为闭壳层体系(Closed-shell),在DMol3中不需要选中Spinrestricted前面的选项。当Alpha电子和Beta电子的数目不相等的时候，我们将会使用不同的＜pi.来表述Alpha电子和Beta电子，这种体系称之为开壳层体系(Opened-shell)或者者自旋极化。在DMol3计算的时候需要选中Spinrestricted的选项，并指定自旋数目。totalalphabeta

在开壳层体系中，会有两个不同的电子密度：一个是Alpha电子的电

子密度，一个是Beta电子的电子密度。它们的和就是整个体系的总电荷，它们的差就是自

旋密度。

spinalphabeta

materialsstudio

当我们确定了电子密度p之后，传统的薛定谓方程会从对电子波函数的处理转换为对电子密度p进行处理的函数。

其中，动能项的方程为：

注意，动能项实际上是一个常数项，在第一次计算完成后，该数值基本上可以确定，后继计算中，则可以忽略这一步骤。

势能项处理：

materialsstudio

需要注意的是，在势能项中的电子-电子相互作用，指的是两个电子间的相互作