| 查看: 2476 | 回复: 12 | ||
[求助]
ENCUT测试
|
|
体系为280个原子，选用PAW_PBE势。2种原子A,B的 ENMIN ，ENMAX分别为 A B 184:245 300:400 为了测试ENCUT，之前我首先在ENCUT=700时测试了KPOINTS，INCAR如下 SYSTERM ISTART=0 #ICHARG= LWAVE=.False. PREC=high ALGO=Fast LREAL=Auto LELF=.TRUE. #LCHARG=.FALSE. ################electronic relaxation############ NELM=80 ENCUT=700 EDIFF=1E-05 #ISPIN=2 ###################ion relaxtion################ #NSW=400 #IBRION=1 #ISIF=2 #EDIFFG=-0.001 ISMEAR=-5 分别测试了1x1x1，2x2x2，3x3x3，4x4x4的K点，KPOINTS格式如下： Automatic 0 Gamma x x x 0 0 0 其能量分别为-2290.2356 -2290.7178 -2290.7170 -2290.7176，所以测试ENCUT时我采用了2x2x2的KPOINTS。 分别测试ENCUT=400 500 600 700 800 900 1000，能量如下： 400 -2289.7561 500 -2290.1378 600 -2290.1329 700 -2290.7178 800 -2290.7593 900 -2290.8174 1000 -2290.8880 我的问题是为什么都达到1000了，能量还没有收敛。哪里出现了问题吗？ |
回复此楼
emilyoyang
木虫 (正式写手)
- 应助: 232 (大学生)
- 金币: 4340.4
- 红花: 42
- 帖子: 667
- 在线: 189.9小时
- 虫号: 1630291
- 注册: 2012-02-20
- 性别: MM
- 专业: 理论和计算化学
2楼2013-01-09 13:43:16
3楼2013-01-09 14:54:23
emilyoyang
木虫 (正式写手)
- 应助: 232 (大学生)
- 金币: 4340.4
- 红花: 42
- 帖子: 667
- 在线: 189.9小时
- 虫号: 1630291
- 注册: 2012-02-20
- 性别: MM
- 专业: 理论和计算化学
4楼2013-01-09 15:54:10
colinchen2
金虫 (小有名气)
- 应助: 15 (小学生)
- 金币: 1652
- 帖子: 59
- 在线: 33.7小时
- 虫号: 1536228
- 注册: 2011-12-13
- 性别: GG
- 专业: 半导体物理
5楼2013-01-09 16:00:43
6楼2013-01-09 16:03:03
7楼2013-01-09 16:40:14
emilyoyang
木虫 (正式写手)
- 应助: 232 (大学生)
- 金币: 4340.4
- 红花: 42
- 帖子: 667
- 在线: 189.9小时
- 虫号: 1630291
- 注册: 2012-02-20
- 性别: MM
- 专业: 理论和计算化学
8楼2013-01-09 17:30:34
|
running on 11 nodes distr: one band on 1 nodes, 11 groups vasp.5.2.2 15Apr09 complex POSCAR found type information on POSCAR N Si POSCAR found : 2 types and 280 ions scaLAPACK will be used LDA part: xc-table for Pade appr. of Perdew POSCAR, INCAR and KPOINTS ok, starting setup FFT: planning ...( 1 ) reading WAVECAR entering main loop N E dE d eps ncg rms rms(c) DAV: 1 0.648341706010E+04 0.64834E+04 -0.62749E+05 6292 0.148E+03 DAV: 2 -0.198633704846E+04 -0.84698E+04 -0.82035E+04 7590 0.316E+02 DAV: 3 -0.260881746868E+04 -0.62248E+03 -0.61848E+03 7282 0.898E+01 DAV: 4 -0.262567178229E+04 -0.16854E+02 -0.16800E+02 7876 0.148E+01 DAV: 5 -0.262604758145E+04 -0.37580E+00 -0.37522E+00 7381 0.226E+00 0.103E+02 RMM: 6 -0.229243796515E+04 0.33361E+03 -0.13557E+03 6700 0.332E+01 0.288E+01 RMM: 7 -0.229130484895E+04 0.11331E+01 -0.13174E+02 7342 0.996E+00 0.187E+01 RMM: 8 -0.229118137113E+04 0.12348E+00 -0.16175E+01 7178 0.440E+00 0.578E+00 RMM: 9 -0.229007908022E+04 0.11023E+01 -0.38694E+00 7377 0.183E+00 0.131E+00 RMM: 10 -0.229033387579E+04 -0.25480E+00 -0.28771E+00 7375 0.147E+00 0.247E+00 RMM: 11 -0.229013411261E+04 0.19976E+00 -0.59105E-01 7374 0.741E-01 0.318E-01 RMM: 12 -0.229013715473E+04 -0.30421E-02 -0.43355E-02 7378 0.204E-01 0.451E-02 RMM: 13 -0.229013781928E+04 -0.66455E-03 -0.17970E-03 7418 0.521E-02 0.856E-02 RMM: 14 -0.229013780946E+04 0.98257E-05 -0.24151E-03 7377 0.437E-02 0.276E-02 RMM: 15 -0.229013780307E+04 0.63870E-05 -0.86009E-05 6487 0.860E-03 1 F= -.22901378E+04 E0= -.22901378E+04 d E =0.000000E+00 ----------------------------------------------- 这是ENCUT=500时，自洽计算的输出，INCAR，KPOINTS（2x2x2）如1楼所示。体系没有磁性，“虽然你自洽计算结束了 但计算的能量可能不对”什么意思，为什么自洽计算结束了，能量会不正确？ |
9楼2013-01-09 20:52:11
10楼2013-01-09 20:56:06
50