| 查看: 827 | 回复: 0 | ||
| 当前主题已经存档。 | ||
[资源]
【转帖】VASP计算NiAl
|
||
|
VASP计算NiAl 发信人: valenhou (AAA), 信区: Gaussian 标 题: Re: VASP计算NiAl 发信站: BBS 大话西游站 (Sat Aug 30 09:17:56 2003), 转信 可以在INCAR文件自己定义NBANDS，但是一般要按下面的规则来确定一个大致的数： 取整数（体系总的价电子/2.0） ＋ 一个整数（至少在5以上，一般用10就好了） 另外的一个办法是，利用所编译出来的makeparam来确定NBANDS的值，这个 makeparam可以检查你所输入是否错误，还可以产生倒格子空间的k点，即IBZKPT文件； 另外可以知道默认的NBANDS、NGX、NGY、NGZ，你可以适当的把这几个数加大些（如果 PREC=High，就没有必要加大了）；还可以知道所要提交的作业将会需要多少的内存。 NBANDS值的改变会导致读入不了WAVECAR和CHGCAR，这个可能与编译器和版本有关， 对4.5.5版本，不会出现此问题，4.4.5版本的会出现，但是用不同的编译器，情况 不同。 如果NGX,NGY,NGZ的值变化了，那就肯定读入不了前次的WAECAR和CHGCAR。 电子收敛很慢，能不能把你的INCAR文件和POSCAR贴出来一下？？？？ 【 在 helpme (dd) 的大作中提到: 】 : ISPIN=1计算以后，得到WAVECAR和CHGCAR， : 但是在ISPIN＝2的计算中，因为NBANDS变化了(原来是72，现在是81）， : 所以实际上并没有用那个WAVECAR。 : 是不是可以在spin＝1的计算中加入：NBANDS＝81， : 然后用这个WAVECAR来计算? : 我现在的计算，在60个电子步内，根本就看不到收敛的影子。 发信人: valenhou (AAA), 信区: Gaussian 标 题: Re: VASP计算NiAl 发信站: BBS 大话西游站 (Mon Sep 1 18:00:17 2003), 转信 1、KPOINTS中取的网格大小不太好，两个数不要相等，因为你的a,b是不一样长的。 2、你的这个slab较大，在运行的时候，提示会说“你可以用实空间的投影算符” ，即设置LREAL=.TRUE.，当你要相当精确的结果时，再用LREAL=.FALSE.计算一下。 你想快点得到大致的结果，可以设置这个，提高计算速度。 3、你用ISIF来设置，让原子移动，这样不太好，因为这种方式，会让所有的原子都移动 起来，没有一个固定点，一般slab的中间层要固定下来。你可以改用Selective dynamics的方式，在POSCAR中的Direct前增加一行"Selective dynamics",然后在 每个原子的坐标的一行后按下面的格式： 0.00000000000000 0.00000000000000 0.50000000000000 F F F 你想要那个原子的坐标在某个分量上移动，这设为T，不移动则为F。 4、你用IBRION=2时,最好改变一下POTIM的默认值，让每一次移动的步骤略大一点也没有 关系，但是不要太夸张。比如设为1.0。 5、在开始试找合适的计算参数时，可以先自己定义一个合适的ENCUT（一般比POTCAR中的 最大的ENMAX大10％，或按PREC=Accurate时，ENCUT所给定的值），再把PREC设为High。 因为ENCUT的值是非常影响计算速度的。 【 在 helpme (dd) 的大作中提到: 】 : POSCAR: : NiAl : 2.895900000000000 : 1.0000000000000000 0.0000000000000000 0.0000000000000000 : 0.0000000000000000 1.4142135623730950 0.0000000000000000 : 0.0000000000000000 0.0000000000000000 12.0208152801700000 : 9 9 : Direct : 0.0000000000000000 0.0000000000000000 0.0025921520508667 : 0.0000000000000000 0.5000000000000000 0.0594327710792201 : 0.0000000000000000 0.0000000000000000 0.1171757892364726 |
回复此楼
