作业帮如何来调整siesta的收敛速度
来源：学生作业帮助网 编辑：作业帮 时间：2024/05/12 14:50:41
如何来调整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的样子.