如图,我只改了下收敛精度,这收敛效果也差太多了吧 S`)@%ZFX{Z_G(MFUU`I%XS2.png 1ASSL{7}W(@~GU}O9H8D86A.png 返回小木虫查看更多
你的.castep里面 general parameters 会有一行 random number generator seed = xxxx 之类。 每次相同设置这个随机数会不同。但是可以在 .param 文件指定这个随机数种子来达到完全相同的计算 这个随机数种子似乎跟第一个电子步的初猜有关。 这个随机数种子不同,初猜和计算过程都会不一样。但是由于是数值求解,所以你最终的算出的总能不会有很大偏差。 如果计算在第一个电子步不收敛的话跟这个随机数种子没有关系: 首先确定是否考虑自旋极化,是否是正确的自旋设置; 如果DFT+U,建议先用DFT优化一边的结果再做DFT+U; 如果有重元素,建议用模守恒赝势试一下; 如果还不行,把所有的正价元素的 formal charge 设置为化合价,把所有的负价元素的 formal charge 设置为 0 (对CASTEP在某些情况下有用,没有依据,不确保有用。) 再不行,在计算资源允许的情况下做 EDFT 计算。(这个被认为是一定可以收敛的。)
不同的精度,收敛步骤肯定不一样啊,一般发文章都需要FINE 的精度,建议一些细微参数,还是要做下收敛性测试
你的.castep里面 general parameters 会有一行
random number generator seed = xxxx 之类。
每次相同设置这个随机数会不同。但是可以在 .param 文件指定这个随机数种子来达到完全相同的计算
这个随机数种子似乎跟第一个电子步的初猜有关。
这个随机数种子不同,初猜和计算过程都会不一样。但是由于是数值求解,所以你最终的算出的总能不会有很大偏差。
如果计算在第一个电子步不收敛的话跟这个随机数种子没有关系:
首先确定是否考虑自旋极化,是否是正确的自旋设置;
如果DFT+U,建议先用DFT优化一边的结果再做DFT+U;
如果有重元素,建议用模守恒赝势试一下;
如果还不行,把所有的正价元素的 formal charge 设置为化合价,把所有的负价元素的 formal charge 设置为 0
(对CASTEP在某些情况下有用,没有依据,不确保有用。)
再不行,在计算资源允许的情况下做 EDFT 计算。(这个被认为是一定可以收敛的。)
不同的精度,收敛步骤肯定不一样啊,一般发文章都需要FINE 的精度,建议一些细微参数,还是要做下收敛性测试
这个精度指的是scf tolerance的精度
非常感谢楼主的回复,我确实在做一个过渡金属磁性体系的加u计算==0好难收敛啊,无论用模守恒还是超软,我有用EDFT进行计算,大约为7到10小时跑一次计算,因为确实很慢而且后续能带要更多的k点,我跑了一周看还没收敛就放弃了,但若EDFT一定能够收敛(问一下大佬这方面相关的文献)’--0那就用好了(收敛性测试可以不加u做吗)另外就是在castep下我只能对某种原子进行加u不知道可否对某些原子进行加u。还有就是楼主您提到的形式电荷的设置,我一般都用默认值了,(会用您的方法试试)很感谢您的回复和所提供的帮助
收敛性测试主要是用来考察计算效率和计算结果是否一致来确定最经济的参数的,收敛与否跟其他的设置有关。k点和截断能一般对收敛与否没有很大影响
模守恒的话需要很高的截断能,一般在7,800eV以上。
如果你确定你各个原子的自旋设置没错的话,(不清楚的话就设成 undefined)
你最开始要试试先不加U算一遍得到一个结果,之后再在这个结果之上再做+U计算看能不能收敛(而且也节省时间一些)。
EDFT是最后实在没办法的手段,他需要的时间差不多是4,5倍以上,所以基本很少用这种的。
CASTEP文档里面有教你修改 density mixing 的参数的,个人感觉这个效果不大。
一般有很多用VASP的文献用超软赝势在比较低的截断能就收敛的,CASTEP就一定有某种设置(不用EDFT)可以收敛(这两个软件是同源的)。
另外CASTEP在新一些的版本(2015以后)里有 TPSD 和 BFGS 两种算法(damped MD看到的不多)。
TPSD更快
BFGS慢一些(感觉是慢2倍以上)但是有些TPSD比较难收敛的计算BFGS更适合,
我对中间具体的算法不是特别了解,自己也是各种试的一些经验,试试形式电荷吧,我的CeO2表面死活不收敛,但是设置形式电荷感觉还可以。我怀疑这个可能跟初猜有关,但是没有依据
,
TPSD在ms7.0时就已经写进程序了。是可以直接通过在输入文件写入指令使用了,另外这2种都是用于结构优化使用的啊。我目前是已经用GGA优化了结构,仅是在计算能带时选用加U(因为=。=貌似很多文献说加u不会对结构有太大扰动)不加u计算能带已经能收敛,而且速度很快。所以就是加u尤其是那种需要在p轨道加u的(==0额这个加u方法貌似争议很大,但确实有很多文献在氧的p轨道加u)=。=就一直震荡,那我先试试形式电荷的方法,谢谢你了