CyRhmU.jpeg
南方科技大学公共卫生及应急管理学院2025级博士研究生招生报考通知
查看: 777  |  回复: 0
当前主题已经存档。

zdhlover

荣誉版主 (职业作家)


[资源] 【转帖】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
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zdhlover 的主题更新
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复(可上传附件)
信息提示
请填处理意见