| 查看: 892 | 回复: 7 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
zxzj05荣誉版主 (著名写手)
|
[交流]
【求助】CG不收敛的问题
|
||
|
我在计算声子谱时，经常碰到CG不收敛的问题，在abinit官方网站上的FAQ里面，也有人提出这个问题： During my calculation for the first derivative of the energy, after many iteration steps, the potential residual begins to fluctuate. It has been found that , due to the 'log file', this potential-based CG line minimization has trouble to converge. 有人提议： Sometimes the CG of WF have trouble converging ; especially for systems requiring very large ecut. In those case since the wave functions do not converge accurately the potential can't converge. To solve this you have to resort to an increase of "nline". nline=15 is a very safe value. Also you may increase nnsclo. 我做了下面修改： nline=15 nnsclo=25（不敢太大，非常浪费时间）， 可是错误依然出现。 请大家给点建设性意见，帮我度过难关！ 不胜感激！ [ Last edited by zxzj05 on 2008-12-5 at 10:33 ] |
回复此楼» 猜你喜欢
考研调剂-材料类-284
已经有3人回复
-
0703化学求调剂
已经有9人回复
-
一志愿深大085601材料工程专业（专硕）300分可以调剂去哪
已经有16人回复
-
085801 总分275 本科新能源 求调剂
已经有4人回复
-
368化学求调剂
已经有9人回复
-
334求调剂
已经有10人回复
-
085404，334分，求调剂
已经有3人回复
-
材料专硕调剂
已经有15人回复
-
材料工程调剂
已经有9人回复
-
材料化工总分334求调剂
已经有3人回复
wuchenwf
荣誉版主 (职业作家)
- 应助: 0 (幼儿园)
- 贵宾: 3.433
- 金币: 19425.8
- 散金: 10
- 红花: 4
- 帖子: 3563
- 在线: 1036.1小时
- 虫号: 398569
- 注册: 2007-06-10
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
- 管辖: 第一性原理
6楼2009-04-20 19:22:42
2楼2008-12-04 22:23:02
zxzj05
荣誉版主 (著名写手)
- 1ST强帖: 1
- 应助: 0 (幼儿园)
- 贵宾: 1.242
- 金币: 5912.5
- 红花: 2
- 帖子: 1528
- 在线: 144.5小时
- 虫号: 406147
- 注册: 2007-06-18
- 性别: GG
- 专业: 金属功能材料
- 管辖: 第一性原理
|
相关输入文件： # Crystalline phonon spectrum ndtset 10 #Set 1 : ground state self-consistency getwfk1 0 kptopt1 1 nqpt1 0 tolvrs1 1.0d-18 rfphon1 0 nline 15 nnsclo3 40 nnsclo4 40 #Q vectors for all datasets nqpt 1 qpt2 0.00000000E+00 0.00000000E+00 0.00000000E+00 ……………………………………………………………………………… qpt10 -2.50000000E-01 5.00000000E-01 2.50000000E-01 #Set 2 : Response function calculation of d/dk wave function iscf2 -3 # Need this non-self-consistent option for d/dk kptopt2 2 # Modify default to use time-reversal symmetry rfphon2 0 # Cancel default rfelfd2 2 # Calculate d/dk wave function only tolvrs2 0.0 # Cancel default for d/dk tolwfr2 1.0d-22 # Use wave function residual criterion instead #Set 3 : Response function calculation of Q=0 phonons and electric field pert. getddk3 2 # d/dk wave functions from last dataset kptopt3 2 # Modify default to use time-reversal symmetry rfelfd3 3 # Electric-field perturbation response only #Sets 4-10 : Finite-wave-vector phonon calculations (defaults for all datasets) getwfk 1 # Use GS wave functions from dataset1 kptopt 3 # Need full k-point set for finite-Q response rfphon 1 # Do phonon response rfatpol 1 3 # Treat displacements of all atoms rfdir 1 1 1 # Do all directions (symmetry will be used) tolvrs 1.0d-8 # This default is active for sets 3-10 nband 25 ecut 40 ngkpt 4 4 4 nshiftk 4 shiftk 0.5 0.5 0.5 0.5 0.0 0.0 0.0 0.5 0.0 0.0 0.0 0.5 iscf 5 错误输出： scfcge : ERROR - Potential-based CG line minimization not converged after 13 restarts. Action : read the eventual warnings about lack of convergence. Some might be relevant. Otherwise, raise nband. 增加nband，ecut只会增加计算量，浪费时间，无法真正解决问题。 |
3楼2008-12-05 10:58:38
baoer3808
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 613.6
- 散金: 15
- 帖子: 172
- 在线: 17.4小时
- 虫号: 725583
- 注册: 2009-03-18
- 专业: 凝聚态物性I:结构、力学和
5楼2009-04-20 09:11:09
