24小时热门版块排行榜    

CyRhmU.jpeg
南方科技大学公共卫生及应急管理学院2026级博士研究生招生报考通知(长期有效)
查看: 2435  |  回复: 6

奔跑的爷们

银虫 (正式写手)

[求助] 求教各位大侠,关于vasp的运行记录

如题,以下是vasp的运行记录,想就运行记录请教各位大侠几个问题:
1. |      For optimal performance we recommend that you set                      |
|        NPAR = approx SQRT( number of cores)                                 |
|      This will greatly improve the performance of VASP for DFT.  
这段的意思是要设置NPAR吗?我没有设置但是我的八个核仍然还在运转啊
2.WARNING: small aliasing (wrap around) errors must be expected
为什么会出现这个警告?
3. 为什么经过了一段:
      N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1     0.395219316642E+03    0.39522E+03   -0.18815E+04  5424   0.184E+03
……
DAV:  20    -0.742393695060E+02   -0.51618E-04   -0.48508E-06  3728   0.101E-02
之后,紧接着又来了一段
N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.742916604512E+02   -0.52343E-01   -0.17052E+01  5424   0.270E+01    0.110E+00
……
DAV:   7    -0.742797664941E+02    0.15331E-04   -0.27305E-04  4544   0.115E-01
此时电子的能量收敛标准(EDIFF=1E-4)已经达到了,为何还要再从DAV1开始计算,而且重复数次?
   1 F= -.74239370E+02 E0= -.74167302E+02  d E =-.742394E+02
curvature:   0.00 expect dE= 0.000E+00 dE for cont linesearch  0.000E+00
trial: gam= 0.00000 g(F)=  0.308E+00 g(S)=  0.190E-02 ort = 0.000E+00 (trialstep = 0.100E+01)
search vector abs. value=  0.310E+00
bond charge predicted
这一段话包含了那些意义呢?
请各位指教

下面是log文件:

-----------------------------------------------------------------------------
|                                                                             |
|           W    W    AA    RRRRR   N    N  II  N    N   GGGG   !!!           |
|           W    W   A  A   R    R  NN   N  II  NN   N  G    G  !!!           |
|           W    W  A    A  R    R  N N  N  II  N N  N  G       !!!           |
|           W WW W  AAAAAA  RRRRR   N  N N  II  N  N N  G  GGG   !            |
|           WW  WW  A    A  R   R   N   NN  II  N   NN  G    G                |
|           W    W  A    A  R    R  N    N  II  N    N   GGGG   !!!           |
|                                                                             |
|      For optimal performance we recommend that you set                      |
|        NPAR = approx SQRT( number of cores)                                 |
|      This will greatly improve the performance of VASP for DFT.             |
|      The default NPAR=number of cores might be grossly inefficient          |
|      on modern multi-core architectures or massively parallel machines.     |
|      Unfortunately you need to use the default for hybrid, GW and RPA       |
|      calculations.                                                          |
|                                                                             |
-----------------------------------------------------------------------------

LDA part: xc-table for Pade appr. of Perdew
POSCAR, INCAR and KPOINTS ok, starting setup
WARNING: small aliasing (wrap around) errors must be expected
FFT: planning ...(           1 )
WAVECAR not read
entering main loop
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1     0.395219316642E+03    0.39522E+03   -0.18815E+04  5424   0.184E+03
……
DAV:  20    -0.742393695060E+02   -0.51618E-04   -0.48508E-06  3728   0.101E-02
   1 F= -.74239370E+02 E0= -.74167302E+02  d E =-.742394E+02
curvature:   0.00 expect dE= 0.000E+00 dE for cont linesearch  0.000E+00
trial: gam= 0.00000 g(F)=  0.308E+00 g(S)=  0.190E-02 ort = 0.000E+00 (trialstep = 0.100E+01)
search vector abs. value=  0.310E+00
bond charge predicted
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.742916604512E+02   -0.52343E-01   -0.17052E+01  5424   0.270E+01    0.110E+00
……
DAV:   7    -0.742797664941E+02    0.15331E-04   -0.27305E-04  4544   0.115E-01
   2 F= -.74279766E+02 E0= -.74207676E+02  d E =-.403970E-01
trial-energy change:   -0.040397  1 .order   -0.045188   -0.309599    0.219223
step:   0.5721(harm=  0.5854)  dis= 0.01659  next Energy=   -74.327040 (dE=-0.877E-01)
bond charge predicted
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.743293982141E+02   -0.49616E-01   -0.31287E+00  5424   0.115E+01    0.471E-01
……
DAV:   5    -0.743270296764E+02    0.96815E-04   -0.92136E-04  5512   0.179E-01
   3 F= -.74327030E+02 E0= -.74254782E+02  d E =-.876602E-01
curvature:  -0.29 expect dE=-0.376E-02 dE for cont linesearch -0.130E-05
trial: gam= 0.03777 g(F)=  0.118E-01 g(S)=  0.104E-02 ort = 0.117E-02 (trialstep = 0.914E+00)
search vector abs. value=  0.134E-01
bond charge predicted
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.743306820981E+02   -0.35556E-02   -0.56505E-01  5424   0.517E+00    0.410E-01
……
DAV:   5    -0.743301754215E+02    0.21081E-04   -0.20817E-04  5384   0.846E-02
   4 F= -.74330175E+02 E0= -.74257952E+02  d E =-.314575E-02
trial-energy change:   -0.003146  1 .order   -0.003170   -0.011806    0.005467
step:   0.6250(harm=  0.6250)  dis= 0.00454  next Energy=   -74.331064 (dE=-0.403E-02)
bond charge predicted
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.743310436880E+02   -0.84719E-03   -0.55433E-02  5424   0.162E+00    0.134E-01
……
DAV:   4    -0.743309981039E+02    0.23216E-04   -0.76480E-05  3560   0.623E-02
   5 F= -.74330998E+02 E0= -.74258764E+02  d E =-.396843E-02
curvature:  -0.32 expect dE=-0.616E-03 dE for cont linesearch -0.130E-06
trial: gam= 0.14599 g(F)=  0.187E-02 g(S)=  0.263E-04 ort =-0.732E-04 (trialstep = 0.857E+00)
search vector abs. value=  0.216E-02
bond charge predicted
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.743322307005E+02   -0.12094E-02   -0.58257E-02  5424   0.178E+00    0.105E-01
……
DAV:   4    -0.743322525693E+02    0.37658E-05   -0.28532E-05  3608   0.353E-02
   6 F= -.74332253E+02 E0= -.74260047E+02  d E =-.125447E-02
trial-energy change:   -0.001254  1 .order   -0.001259   -0.001617   -0.000901
step:   1.9360(harm=  1.9360)  dis= 0.00553  next Energy=   -74.332825 (dE=-0.183E-02)
bond charge predicted
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.743327647452E+02   -0.50841E-03   -0.92389E-02  5424   0.224E+00    0.137E-01
……
DAV:   4    -0.743328070278E+02    0.79076E-05   -0.48127E-05  3608   0.481E-02
   7 F= -.74332807E+02 E0= -.74260611E+02  d E =-.180892E-02
curvature:  -1.11 expect dE=-0.109E-02 dE for cont linesearch -0.649E-07
trial: gam= 0.52253 g(F)=  0.930E-03 g(S)=  0.516E-04 ort =-0.113E-04 (trialstep = 0.107E+01)
search vector abs. value=  0.156E-02
bond charge predicted
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.743332034028E+02   -0.38847E-03   -0.77433E-02  5424   0.204E+00    0.144E-01
……
DAV:   4    -0.743331834775E+02    0.23471E-04   -0.78441E-05  3616   0.614E-02
   8 F= -.74333183E+02 E0= -.74260994E+02  d E =-.376450E-03
trial-energy change:   -0.000376  1 .order   -0.000399   -0.001047    0.000248
step:   0.8668(harm=  0.8668)  dis= 0.00219  next Energy=   -74.333230 (dE=-0.423E-03)
bond charge predicted
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.743332115200E+02   -0.45712E-05   -0.30759E-03  5424   0.403E-01    0.272E-02
DAV:   2    -0.743332168524E+02   -0.53324E-05   -0.72029E-05  3688   0.613E-02
   9 F= -.74333217E+02 E0= -.74260994E+02  d E =-.409825E-03
curvature:  -0.69 expect dE=-0.534E-04 dE for cont linesearch -0.146E-06
trial: gam= 0.11587 g(F)=  0.560E-04 g(S)=  0.210E-04 ort =-0.181E-04 (trialstep = 0.103E+01)
search vector abs. value=  0.938E-04
reached required accuracy - stopping structural energy minimisation
writing wavefunctions

回复此楼

» 收录本帖的淘帖专辑推荐

第一性原理相关文档

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

这家伙太懒了,什么都没有填……
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

米色空白

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
奔跑的爷们: 金币+10, ★★★很有帮助 2013-10-20 10:27:55
1.NPAR是并行设置,多个核一起计算的时候最好设置成总核数开根号,这样计算效率更高,但是不设也无妨,而且如果遇到开根号不是整数的时候也可以选择使用别的数,比如32个核设成4也是可以算的。
2.起因恐怕是很深层次的量子力学或者固体物理方面的问题,我们也不是很清楚,但是可以通过INCAR参数调整得更精密来避免,不过一般情况下可以忽略这点提示,对结果影响不会太大。
3.我不是太明白是什么情况。看底下的LOG文件似乎挺正常的。DAV1…20是表示完成了20个电子迭代步,然后收敛了,于是有了1F=…也就是说完成了第1个离子迭代步,然后进入下一个循环,也就是第2个离子步,这里面有包括了7次电子步。具体的其他信息可以参照VASP手册里面关于OSZICAR的描述。
2楼2013-10-20 10:26:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

奔跑的爷们

银虫 (正式写手)

引用回帖:
2楼: Originally posted by 米色空白 at 2013-10-20 10:26:01
1.NPAR是并行设置,多个核一起计算的时候最好设置成总核数开根号,这样计算效率更高,但是不设也无妨,而且如果遇到开根号不是整数的时候也可以选择使用别的数,比如32个核设成4也是可以算的。
2.起因恐怕是很深层 ...

3.我不是太明白是什么情况。看底下的LOG文件似乎挺正常的。DAV1…20是表示完成了20个电子迭代步,然后收敛了,于是有了1F=…也就是说完成了第1个离子迭代步,然后进入下一个循环,也就是第2个离子步,这里面有包括了7次电子步。具体的其他信息可以参照VASP手册里面关于OSZICAR的描述。
前辈,是不是算完第一个离子迭代步循环,这时候还没有达到EDIFFG的收敛标准,所以要继续进行计算?换句话说,如果EDIFFG的标准越高,那么这种循环的次数也就越多?不知道您是否知道我要表达的意思
这家伙太懒了,什么都没有填……
3楼2013-10-20 10:29:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

米色空白

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
奔跑的爷们: 金币+10, ★★★★★最佳答案 2013-10-20 10:54:40
引用回帖:
3楼: Originally posted by 奔跑的爷们 at 2013-10-20 10:29:23
3.我不是太明白是什么情况。看底下的LOG文件似乎挺正常的。DAV1…20是表示完成了20个电子迭代步,然后收敛了,于是有了1F=…也就是说完成了第1个离子迭代步,然后进入下一个循环,也就是第2个离子步,这里面有包括 ...

应该是EDIFF和EDIFFG两个参数吧,默认的话前者是后者的10倍。当电子迭代到能量差小于EDIFFG的时候电子迭代结束,完成第一个离子步,但是此时你的体系的离子能量波动还很大( d E =-.742394E+02),还远没有达到EDIFF要求的标准,于是乎再将这个结构进行微扰,再一个个电子步地迭代,直到这次又迭代了7次终于达到了EDIFFG的标准,第2个离子步完成,但是离子步还没完(d E =-.403970E-01),所以还要这样周而复始地算,直到9 F= -.74333217E+02 E0= -.74260994E+02  d E =-.409825E-03终于达到了EDIFF的标准,于是结束了
4楼2013-10-20 10:41:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

奔跑的爷们

银虫 (正式写手)

引用回帖:
4楼: Originally posted by 米色空白 at 2013-10-20 10:41:52
应该是EDIFF和EDIFFG两个参数吧,默认的话前者是后者的10倍。当电子迭代到能量差小于EDIFFG的时候电子迭代结束,完成第一个离子步,但是此时你的体系的离子能量波动还很大( d E =-.742394E+02),还远没有达到EDI ...

多谢前辈!您的回答正是我所想知道的,另外问一句,9 F= -.74333217E+02 E0= -.74260994E+02分别代表什么?应该是总能一类的东西吗?为什么有两个呢?
这家伙太懒了,什么都没有填……
5楼2013-10-20 10:55:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

米色空白

金虫 (小有名气)

引用回帖:
5楼: Originally posted by 奔跑的爷们 at 2013-10-20 10:55:34
多谢前辈!您的回答正是我所想知道的,另外问一句,9 F= -.74333217E+02 E0= -.74260994E+02分别代表什么?应该是总能一类的东西吗?为什么有两个呢?...

9表示离子步数,F= -.74333217E+02表示收敛到这个程度的体系的总自由能是-74.333217eV,其中参考原子的能量已经被扣除;E0=-.74260994E+02代表的是把INCAR里面那个sigma外推到等于0的时候的能量是-74.260994。往往这个值在OUTCAR里面也能差得到。
6楼2013-10-20 16:41:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

奔跑的爷们

银虫 (正式写手)

内容已删除
这家伙太懒了,什么都没有填……
7楼2013-10-20 18:12:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 奔跑的爷们 的主题更新
信息提示
请填处理意见