| 查看: 950 | 回复: 2 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
charles0116金虫 (初入文坛)
|
[求助]
高手看过来：内部原子弛豫时自洽计算max grad (force/stress）无法满足收敛条件tolmxf
|
|
|
优化结构之前是要先进行内部原子弛豫吧？ 但在做内部原子弛豫时（optcell=0 ionmov=2 or 3) 单次scf在开始的几轮满足条件退出，后面就时而converged,时而not converged,这是为什么呢？ e.g. scprqt: WARNING - nstep= 200 was not enough SCF cycles to converge; potential residual= 1.475E-05 exceeds tolvrs= 1.000E-10 nstep设为200己经差不多了吧？再调大也只是浪费时间？ 最让人不解的是，在前十几轮成功converged scf计算或是后面一直not enough SCF cycles to converge 的scf计算每轮次末尾都会有 max grad (force/stress) = ××××× > tolmxf= 5.0000E-05 ha/bohr (free atoms) 类似字样，而开始下一轮scf计算，直到in文件中设置ntime的100次scf计算完成，整个计算才退出完成。 为此，我调整了occopt=1 or 7 ;iscf=7 or 5; 还调整了各种收敛参数 tolmxf; tolvrs or toldff or tolrff or tolfe ; diemac 等值，也曾尝试过直接进行几何优化（optcell=/0 ionmov=2 or 3)但都无法在ntime之内收敛。 请问各位高手我这是什么问题呢？影响max grad (force/stress)值的因素有什么呢？ 我的认为不管是内部原子弛豫或是后来的几何优化，都应该是满足自洽计算,即 max grad (force/stress) < tolmxf的。不知我这样的理解对吗？ in文件内容如下： #Structural optimisation #************************************ nshiftk 1 shiftk 0 0 0 #********************************* ecut 26 ecutsm 0.5 ngkpt 3 3 4 occopt 7 tsmear 0.01 #********************************* #Common input data chksymbreak 0 chkprim 0 brvltt -1 #Definition of the unit cell #********************************** acell 2*11.135000 4.376000 Angstroms angdeg 90 90 90 #Definition of the atoms #********************* spgroup 82 ntypat 5 znucl 48 30 16 7 6 natom 14 natrd 5 typat 1 2 3 4 5 xred 0.000000 0.000000 0.000000 0.000000 0.500000 0.250000 -0.124500 0.152800 -0.302900 -0.048000 0.359600 0.006400 -0.078300 0.273400 -0.117700 #Definition of the scf procedure #****************************** iscf 7 nstep 300 ixc 1 toldfe 5.0d-6 diemac 10.0 #Structural relaxation #************************************ ionmov 2 optcell 0 tolmxf 5.0d-5 ntime 100 #************************************ out文件我附在了帖子中 图片是我摘出的每轮次scf结束的收敛判定条件，感觉开始max grad (force/stress)还貌似有些下降的趋势，可是后面就毫无规律。  2013-08-20 21:06:34的屏幕截图.png |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : zctc.out
2013-08-20 21:35:14, 250.75 K
» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有52人回复
-
时间戳又变了8-15
已经有7人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有9人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有6人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有7人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
filecode
已经有11人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
charles0116
金虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 1062.5
- 红花: 1
- 帖子: 43
- 在线: 112.5小时
- 虫号: 1684499
- 注册: 2012-03-12
- 专业: 凝聚态物性 II ：电子结构
3楼2013-08-23 01:59:25
charles0116
金虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 1062.5
- 红花: 1
- 帖子: 43
- 在线: 112.5小时
- 虫号: 1684499
- 注册: 2012-03-12
- 专业: 凝聚态物性 II ：电子结构
2楼2013-08-21 11:27:09
