| 查看: 1107 | 回复: 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 |
» 猜你喜欢
基础研究怎么拉横向，学校到款任务越来越多，难以完成 拉横向，都有哪些途径啊
已经有10人回复
-
2027年申博
已经有3人回复
-
一个有机合成实验室都需要哪些设备？
已经有7人回复
-
2026年国自然面上资助率
已经有20人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有22人回复
-
微信指数没变化，科研之友没阅读
已经有18人回复
-
HXDI做水性聚氨酯乳液，是不是特别容易出渣
已经有3人回复
-
系统今天又提示维护了，估计离放榜不远了
已经有15人回复
-
你们的时间戳变了吗
已经有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
- 金币: 19427.4
- 散金: 10
- 红花: 4
- 帖子: 3566
- 在线: 1037小时
- 虫号: 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
- 金币: 19427.4
- 散金: 10
- 红花: 4
- 帖子: 3566
- 在线: 1037小时
- 虫号: 398569
- 注册: 2007-06-10
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
- 管辖: 第一性原理
8楼2009-10-16 20:03:34