| 查看: 1061 | 回复: 7 | ||||
| 当前主题已经存档。 | ||||
爱惜花草木虫 (著名写手)
|
[交流]
【讨论】我的这个任务结构优化时间过长是不是体系设置有问题
|
|||
|
我INCAR里面设置的NSW=600（当然一般运动不了这么多步数） 做的体系有点大，大概一个cell里有近60个离子 开始做结构优化，现在已经跑了200多圈了还没有达到平衡位置 不过输出文件里仔细检查了能量值，发现总体趋势是在下降的，只是变化不大，大概在0.0001这个数量级左右，也没有出现那种振荡的现象 这种情况是否应该将任务kill？初始结构给的有问题？ PS：精度设置也不是很高 EDIFF = 1E-3 EDIFFG = -1E-2 IBRION=2 POTIM=0.2 [ Last edited by 爱惜花草 on 2009-10-16 at 10:12 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 VASP |
» 猜你喜欢
论文撤稿了
已经有8人回复
-
博士申请
已经有7人回复
-
青B发送上会通知了吗
已经有8人回复
-
化学专业申博
已经有4人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
-
评审有感
已经有27人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有5人回复
jinhuan123
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 529.3
- 散金: 21
- 红花: 1
- 帖子: 106
- 在线: 51.7小时
- 虫号: 805643
- 注册: 2009-07-08
- 专业: 催化化学
2楼2009-10-16 08:26:43
wuchenwf
荣誉版主 (职业作家)
- 应助: 0 (幼儿园)
- 贵宾: 3.433
- 金币: 19426.6
- 散金: 10
- 红花: 4
- 帖子: 3564
- 在线: 1036.8小时
- 虫号: 398569
- 注册: 2007-06-10
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
- 管辖: 第一性原理
3楼2009-10-16 08:46:31
linggang87
银虫 (小有名气)
- 应助: 4 (幼儿园)
- 金币: 48.4
- 散金: 60
- 红花: 1
- 帖子: 198
- 在线: 108.8小时
- 虫号: 635564
- 注册: 2008-10-24
- 专业: 金属材料的制备科学与跨学
4楼2009-10-16 08:51:47
cenwanglai
荣誉版主 (知名作家)
老和山猥琐派九段
- 1ST强帖: 5
- 应助: 46 (小学生)
- 贵宾: 8.842
- 金币: 7441.4
- 散金: 9654
- 红花: 119
- 沙发: 23
- 帖子: 5306
- 在线: 1961.4小时
- 虫号: 537452
- 注册: 2008-04-01
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 仿真模拟
5楼2009-10-16 09:14:44
爱惜花草
木虫 (著名写手)
- 应助: 7 (幼儿园)
- 金币: 3380.2
- 散金: 762
- 红花: 9
- 沙发: 1
- 帖子: 1173
- 在线: 196.5小时
- 虫号: 805026
- 注册: 2009-07-07
- 专业: 半导体材料
6楼2009-10-16 10:11:18
nkleof
银虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 2093.2
- 散金: 649
- 红花: 5
- 帖子: 813
- 在线: 165.9小时
- 虫号: 152581
- 注册: 2006-01-02
- 专业: 凝聚态物性 II ：电子结构
★ ★ ★ ★
aylayl08(金币+3,VIP+0):感谢你的解答，欢迎常来 10-16 17:35
爱惜花草(金币+1,VIP+0): 10-17 10:14
aylayl08(金币+3,VIP+0):感谢你的解答，欢迎常来 10-16 17:35
爱惜花草(金币+1,VIP+0): 10-17 10:14
|
设置IBRION=2的时候，可以不用设置POTIM，用vasp的默认值，它会在优化结束后的任务报告中给出一个合适的trialstep（在任务报告的最末），下一次用IBRION=1优化的时候就可以采用这个trialstep*0.5作为POTIM的值。在vasp的说明中有介绍。 http://cms.mpi.univie.ac.at/vasp/vasp/node115.html Within all relaxation algorithms ( IBRION=1,2 and 3) the parameter POTIM should be supplied in the INCAR file. For IBRION$>$0, the forces are scaled internally before calling the minimization routine. Therefore for relaxations, POTIM has no physical meaning and serves only as a scaling factor. For many systems, the optimal POTIM is around 0.5. Because the Quasi-Newton algorithm and the damped algorithms are sensitive to the choice of this parameter, use IBRION=2, if you are not sure how large the optimal POTIM is. In this case, the OUTCAR file and stdout will contain a line indicating a reliable POTIM. For IBRION=2, the following lines will be written to stdout after each corrector step (usually each odd step): trial: gam= .00000 g(F)= .152E+01 g(S)= .000E+00 ort = .000E+00 (trialstep = .82) The quantity gam is the conjugation parameter to the previous step, g(F) and g(S) are the norm of the force respectively the norm of the stress tensor. The quantity ort is an indicator whether this search direction is orthogonal to the last search direction (for an optimal step this quantity should be much smaller than ( g(F) + g(S)). The quantity trialstep is the size of the current trialstep. This value is the average step size leading to a line minimization in the previous ionic step. An optimal POTIM can be determined, by multiplying the current POTIM with the quantity trialstep. [ Last edited by nkleof on 2009-10-16 at 10:33 ] |
7楼2009-10-16 10:31:20
wuchenwf
荣誉版主 (职业作家)
- 应助: 0 (幼儿园)
- 贵宾: 3.433
- 金币: 19426.6
- 散金: 10
- 红花: 4
- 帖子: 3564
- 在线: 1036.8小时
- 虫号: 398569
- 注册: 2007-06-10
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
- 管辖: 第一性原理
8楼2009-10-16 20:03:34
10