| 查看: 2436 | 回复: 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  |
回复此楼» 收录本帖的淘帖专辑推荐
 第一性原理相关文档 |
» 猜你喜欢
三甲基碘化亚砜的氧化反应
已经有4人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有5人回复
-
孩子确诊有中度注意力缺陷
已经有12人回复
-
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
-
论文投稿，期刊推荐
已经有4人回复
-
硕士和导师闹得不愉快
已经有13人回复
» 本主题相关价值贴推荐，对您同样有帮助:
VASP建模，超胞模型有几百个原子，如何写程序来建模啊？
已经有5人回复- 求教，关于使用MS建石墨烯超胞模型，然后导入vasp中运算
已经有4人回复
- 用VASP计算单点能的时候总是出错，怎么办啊？大侠来给知道啊？
已经有8人回复
- vasp运行，第一离子步结束后自动停了
已经有7人回复
- 各位大侠，vasp下计算完能带如何作图
已经有12人回复
- 饲养实验过程中鸡死亡了饲料转化效率怎么计算啊？
已经有7人回复
- 想配个台式机算VASP，求各位大侠推荐个配置
已经有10人回复
- 求 河南省新乡县DEM图，不知哪位大侠有，请大侠指条名路，感激不尽
已经有8人回复
- 各位大侠，目前，植物类油脂研究的热点是什么，
已经有15人回复
- 请问大侠，我初学Vasp，对于linux系统应该怎样学？
已经有16人回复
- 制备活性干酵母干燥过程酵母变褐色怎么预防
已经有5人回复
- 求助！关于FM和AFM如何设置的问题
已经有3人回复
- 【求助完结】J1需要导师的CV，能直接问导师要吗？
已经有4人回复
- 【求助】原子的电荷密度分布怎么求？
已经有3人回复
- 【求助】p4vasp出现的问题 哪位大侠帮忙解决
已经有5人回复
米色空白
金虫 (小有名气)
- 应助: 15 (小学生)
- 金币: 1925.4
- 红花: 3
- 帖子: 73
- 在线: 30小时
- 虫号: 2132467
- 注册: 2012-11-17
- 性别: MM
- 专业: 高分子材料的加工与成型
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
奔跑的爷们: 金币+10, ★★★★★最佳答案 2013-10-20 10:54:40
奔跑的爷们: 金币+10, ★★★★★最佳答案 2013-10-20 10:54:40
|
应该是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
米色空白
金虫 (小有名气)
- 应助: 15 (小学生)
- 金币: 1925.4
- 红花: 3
- 帖子: 73
- 在线: 30小时
- 虫号: 2132467
- 注册: 2012-11-17
- 性别: MM
- 专业: 高分子材料的加工与成型
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
奔跑的爷们: 金币+10, ★★★很有帮助 2013-10-20 10:27:55
感谢参与，应助指数 +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
3楼2013-10-20 10:29:23
5楼2013-10-20 10:55:34