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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 求教各位大侠，关于vasp的运行记录
北京石油化工学院2026年研究生招生接收调剂公告
51/1返回列表
查看: 2623  |  回复: 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的回帖

奔跑的爷们

银虫 (正式写手)
引用回帖:
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的回帖
查看全部 7 个回答

米色空白

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +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的回帖

米色空白

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
奔跑的爷们: 金币+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的回帖
查看全部 7 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 285求调剂 +5 AZMK 2026-04-02 6/300 2026-04-03 08:41 by yulian1987
[考研] 295求调剂 +7 愿旅途永远坦然 2026-04-02 7/350 2026-04-03 08:22 by fangshan711
[考研] 279求调剂 +5 qazplm0852 2026-04-02 5/250 2026-04-02 22:38 by liveagain126
[考研] 312 化工或制药调剂 +7 小小墨123 2026-04-02 8/400 2026-04-02 22:15 by ZXlzxl0425
[考研] 工科 267求调剂 +3 wanwan00 2026-04-02 5/250 2026-04-02 21:42 by wanwan00
[考研] 一志愿复旦材料，英一专硕，总分357调剂 +4 1050389037 2026-04-02 5/250 2026-04-02 21:40 by dongzh2009
[考研] 296求调剂 +4 sdhu 2026-04-02 4/200 2026-04-02 21:29 by baoball
[考研] 环境科学与工程334分求调剂 +7 王一一依依 2026-03-30 9/450 2026-04-02 21:15 by 1104338198
[考研] 070300化学求调剂 +8 起个名咋这么难 2026-03-27 9/450 2026-04-02 20:42 by dongzh2009
[考研] 085602化工求调剂（331分） +9 111@127 2026-03-30 9/450 2026-04-02 20:00 by dick_runner
[考研] 302求调剂一志愿北航070300，本科郑大化学 +8 圣日耳曼条 2026-04-01 11/550 2026-04-02 07:40 by chemdavid
[考研] 安全工程 285 求调剂 +3 Xinyu56 2026-04-01 4/200 2026-04-01 21:50 by 静静静静静静静�
[考研] 0855机械初试280求调剂 +3 kazenotori 2026-03-31 3/150 2026-04-01 10:08 by JourneyLucky
[考研] 070300化学专业279调剂 +10 哈哈哈^_^ 2026-03-31 10/500 2026-03-31 23:13 by liu823948201
[考研] 本科211安全工程，初试290分，求调剂 +3 2719846834 2026-03-28 3/150 2026-03-31 13:52 by 热情沙漠
[考研] 一志愿西电085401数一英一299求调剂 六级521 +4 爱吃大鸭梨 2026-03-31 4/200 2026-03-31 11:51 by 搏击518
[考研] 085701环境工程求调剂 +11 多久上课 2026-03-27 12/600 2026-03-30 21:21 by 研究僧导导
[考研] 071010 323 分求调剂 +3 Baekzhy 2026-03-27 3/150 2026-03-30 14:24 by andresqi
[考研] 085701求调剂初试286分 +5 secret0328 2026-03-28 5/250 2026-03-30 12:54 by fangnagu
[考研] 298调剂 +3 jiyingjie123 2026-03-27 3/150 2026-03-27 11:57 by wxiongid
信息提示
关闭
请填处理意见
关闭