应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 求教各位大侠，关于vasp的运行记录
北京石油化工学院2026年研究生招生接收调剂公告
71/1返回列表
查看: 2615  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

奔跑的爷们

银虫 (正式写手)

[求助] 求教各位大侠，关于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
回复此楼

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

第一性原理相关文档

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
这家伙太懒了，什么都没有填……
1楼 2013-10-18 16:40:23
已阅   回复此楼   关注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的描述。
赞一下(1人) 回复此楼
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的标准，于是结束了
赞一下(2人) 回复此楼
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里面也能差得到。
赞一下(2人) 回复此楼
6楼2013-10-20 16:41:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

奔跑的爷们

银虫 (正式写手)
内容已删除
回复此楼
这家伙太懒了，什么都没有填……
7楼2013-10-20 18:12:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 奔跑的爷们 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085600 295分求调剂 +17 W55j 2026-03-30 21/1050 2026-04-03 01:10 by BruceLiu320
[考研] 一志愿深大085601材料工程专业（专硕）300分可以调剂去哪 +7 10160315 2026-04-02 7/350 2026-04-02 22:48 by liveagain126
[考研] 372分材料与化工（085600）一志愿湖南大学求调剂 +5 蓝笺片 2026-04-02 6/300 2026-04-02 21:37 by dongzh2009
[考研] 085602 找调剂 +3 逆时针快乐 2026-04-02 3/150 2026-04-02 21:23 by dongzh2009
[考研] 262求调剂 +5 励志一定发文章 2026-04-02 6/300 2026-04-02 12:51 by yulian1987
[考研] 311求调剂 +14 蓝月亮亮 2026-03-30 14/700 2026-04-02 12:18 by 1753564080
[考研] 一志愿北京科技大学085601材料工程英一数二初试总分335求调剂 +9 双马尾痞老板2 2026-04-01 9/450 2026-04-02 12:14 by oooqiao
[考研] 找调剂 +6 楚乔乔 2026-04-01 6/300 2026-04-02 12:06 by xiaoranmu
[考研] 0856，269分求调剂 +8 有学上就行求求� 2026-03-30 11/550 2026-04-01 22:33 by 2026材料调剂
[考研] 材料调剂 +11 一样YWY 2026-03-31 11/550 2026-04-01 22:25 by zhouyuwinner
[考研] 332求调剂 +8 Lyy930824@ 2026-03-29 8/400 2026-04-01 18:40 by 千钧澄玉宇啊
[考研] 301求调剂 +8 axibli 2026-04-01 8/400 2026-04-01 09:51 by 我的船我的海
[考研] 考研材料工程351分调剂 +5 整个好的 2026-03-31 5/250 2026-04-01 09:36 by topgun2009
[考研] 一志愿北京科技大学085601材料工程英一数二初试总分335求调剂 +5 双马尾痞老板2 2026-03-31 5/250 2026-04-01 09:04 by oooqiao
[考研] 求调剂：085600材料与化工，考材科基，总分319 +17 678lucky 2026-03-31 21/1050 2026-04-01 01:40 by 1018329917
[考研] 生物考研337分求调剂 +4 cgxin 2026-03-30 6/300 2026-03-31 14:18 by 记事本2026
[考研] 266分，求材料冶金能源化工等调剂 +8 哇呼哼呼哼 2026-03-27 10/500 2026-03-31 13:35 by Huaxue_Wang
[考研] 调剂求院校招收 +7 鹤鲸鸽 2026-03-28 7/350 2026-03-31 11:21 by oooqiao
[考研] 269求调剂 +4 我想读研11 2026-03-31 4/200 2026-03-31 10:04 by cal0306
[考研] 调剂考研 +3 王杰一 2026-03-29 3/150 2026-03-29 08:09 by fmesaito
信息提示
关闭
请填处理意见
关闭