| 查看: 1001 | 回复: 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 |
» 猜你喜欢
298-一志愿中国农业大学-求调剂
已经有5人回复
-
297求调剂
已经有5人回复
-
302求调剂
已经有10人回复
-
266求调剂
已经有8人回复
-
化学工程321分求调剂
已经有17人回复
-
344求调剂
已经有5人回复
-
312求调剂
已经有9人回复
-
一志愿天津大学化学工艺专业（081702）315分求调剂
已经有9人回复
-
288求调剂，一志愿华南理工大学071005
已经有4人回复
-
304求调剂
已经有13人回复
jinhuan123
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 529.3
- 散金: 21
- 红花: 1
- 帖子: 106
- 在线: 51.7小时
- 虫号: 805643
- 注册: 2009-07-08
- 专业: 催化化学
2楼2009-10-16 08:26:43
wuchenwf
荣誉版主 (职业作家)
- 应助: 0 (幼儿园)
- 贵宾: 3.433
- 金币: 19419.2
- 散金: 10
- 红花: 4
- 帖子: 3560
- 在线: 1035.7小时
- 虫号: 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
- 金币: 7440.4
- 散金: 9654
- 红花: 118
- 沙发: 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
- 金币: 19419.2
- 散金: 10
- 红花: 4
- 帖子: 3560
- 在线: 1035.7小时
- 虫号: 398569
- 注册: 2007-06-10
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
- 管辖: 第一性原理
8楼2009-10-16 20:03:34
