castep参数优质精品文档.docxVIP

CASTEP参数设置

1、castep在做几何优化时，结构收敛判据，也就是Energy,maxforce,maxstress,maxdisplacement,这4个判据，给你分了差中好超好4个级别，这4个级别对应的判据数值是非常合理的。

如果你的体系是一个比较小而且对称性超高的原胞，那么几何判据用Ultrafine级别；如果是复合体系，或体系较大，用fine级别。如果体系超大，用中等级别。当然，所有情况下，你都可以先用中等级别做优化，然后再提升到好的级别。在写论文时，最终结果至少要用到好的级别，但你论文中不能直接说我用了finelevel，你要把每个判据，按照castep的fine设定写出来。

评析：先用中等级别优化，再用fine级别优化，特殊情况采用ultrafine。

cutoff给多少eV，这个你要看看文献一般都用多少，优化的时候少用，有个经验就是，你把catoff那里点开，然后回头点优化质量，差等，中等，好，超好，你每个都点一下，你会发现cutoff不一样，用fine级别的cutoff值去优化，就一定没有问题。根据体系不一样，castep会自己?判定一个在fine级别下cutoff应该是多少，你就用这个值是多少。对于很大的体系castep即使在fine级别判定的cutoff也可能很小这时候是不行的cutoff一般不要低于400eV。

评析：Ecut根据文献值，不低于400eV。优化时可以取小些，性质计算可以大些，比如优化400eV性质计算500eV。

做几何优化的时候，scf的收敛标准也不用那么高，用5e-6eV/atom足以，这个数值比fine差一些，比中等好一些。几何优化做完，做能量和性质计算的时候，scf至少要用1e-6,如果体系很小，你可以用5e-7.

评析：scf标准：几何优化用5.0-6eV/atom或者fine值;计算性质时使用至少1.0-6，如果体系很小，你可以用5e-7.

k-point设置多少也要看文献经验办法是你自己定义kpoint比如nxnxn他下面有一个对应的实空间间隔当你吧kpoint调整到0.3左右的时候就用这套k就可以，这是几何优化一般大体系333中小体系555

优化好结构算能量和性质的时候对应实空间间隔调整到0.1附近就用这套一般做能量性质计算大体系777或999中小体系999或111111或131313视情况而定

评析：k-point设置看文献值，几何优化实空间间隔0.3左右；性质计算实空间间隔0.1附近。大体系设置小的值，比如333,；小体系设置打的值，比如151515。当然一般人懒的做收敛性测试直接用大些151515等。

慕：看你是什么材料晶格有多大

实空间越大倒空间越小

人：晶格3.64A

慕：那你优化的时候用551基本就行了

慕：文献用的很大是因为他懒得做收敛测试了

人：我的真空层大15A

慕：真空无论多少c都用1

二维材料k点都是nn1，最后一个1表示在c方向上电子不流通，如果你是111的话电子就只能在晶格内流通晶格之间不流通这就是gama点。

人：这个551和3.64A有什么关系计算出来吗

慕：我经验判断出来的，优化结构不用那么纠结

但是3.64的实空间晶格如果做能量k点一般15151足矣，没必要搞那么高，文献里弄那么高绝对是他懒得做收敛测试。

2、Densitymixing默认就行，99%不用管，除非scf不收敛才调整。

3、性质计算

空带数目设定为50，保证有足够的虚轨道，不论是计算band，还是dos（castep中光学性质和gaussion等不一样，gaussian只能在gama点计算，相当于111，这个计算的是光学激发，castep计算的光学性质物理意义不明确，许多中国人喜欢计算这个，外国人很少用castep计算光学性质），后面的k-point不用设置，直接设置seperation为0.005，这样画出的能带就很完美了。

注意：kpoint是表示每隔多少做一个能量点，那个点做的越多能量点在能带上分布就越密集，能带看起来就越平滑。因为kpoint即使是fine，seperation才0.015，不够用。

能带计算more中有UseseparateXCfunctionalforbandstructurecalculation，只有在体系中涉及过渡金属,并且晶胞很小（不超过10个原子）才考虑勾选前面的复选框，这是基于GGA的电子密度波函数用非局域密度泛函计算能带，有过渡金属的体系，GGA会低估gap，所以才有HSE去做能带，但是能量迭代时还是GGA,只是做能带时用