作业帮如何来调整siesta的收敛速度
来源：学生作业帮助网 编辑：作业帮 时间：2024/11/24 22:30:39
xTMSP+Y31`WMʢ3]#M5SE&"PZV4g^;ؐ!ysϹzW2lL+(O8SsԻv;,{GV#-wS0,XEɼޡ(oTeř&A P㧦AePMq9
\JY] (ަ$r1ސ_}2/6x_H fȳ6Òs ۂ
/,o^7Wy~R ?QO4GDc NW/ r
i5Y^r= yqS[DA(RЧC`ߘxC[~c&еcOLnXgM_B10.TlBMgUP!L6FA_n"i|19El.=+[^Y4B{H D(Ju$ )Qh%"41.B/rړN=yL I'qHQ A
44eZt^Ԭ~zH]lXY$1FE2hJI%jPO>ʺn=SMC:,w@_Gk@ M՞ < bd{C,UqqOvQ ftd|nI,-wP= a2nNv t܅ʶgT&0v
8RVJ {Ps!V3 F 欼kG/~T.A{|^i3BЌB1 C x*?:P
如何来调整siesta的收敛速度
如何来调整siesta的收敛速度
如何来调整siesta的收敛速度
对于金属的计算, 需要调整各种参数来达到或者加速收敛的速度.
在siesta 中, 其中控制计算收敛的参数有OccupationFunction, ElectronicTemperature,
DM.MixingWeights, DM.NumberPulay, DM.NumberKick,
DM.KickMixingWeight 等.
其 中OccupationFunction有两种占据函数可选择, FD/MP 即Fermi-Dirac 分布和
Methfessel-Paxton分布. 对于金属体系来说,MP占据函数可以更精确地对Fermi面取点, 可以达到加速收敛的目的.
对于MP占据函数, 还有关键词OccupationMPOrder来控制MP占据函数的阶数.在手册上提到,
高阶的展开可以得到好的基态收敛结果. 因此在金属体系相关的计算中, 一般来说选择MP占据函数来达到加速收敛的目的.
同时升高ElectronicTemperature也可以有助于加快体系计算的收敛速度, 但是需要了解一下,
升高温度对体系计算结果的影响.
而 剩下的几个参数, DM.MixingWeights 和DM.NumberPulay 同样是可以达到加速收敛的目的,
如果DM.NumberPulay 设置为0或1 则在计算中是采取线性混合, 为了加快收敛,
可以将DM.NumberPulay设置为4-10的样子(这个数目应当是根据中间结果分析来设置, 如果收敛很慢可以增加.) .
有的时候试图增加或者降低DM.MixingWeights, 其默认是0.25,
在收敛慢的情形下可以尝试将其设置在0.50的样子.