如何来调整siesta的收敛速度
来源:学生作业帮 编辑:神马作文网作业帮 分类:物理作业 时间:2024/11/12 01:51:41
如何来调整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的样子.
在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的样子.