| 查看: 891 | 回复: 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 ] |
回复此楼» 猜你喜欢
328求调剂
已经有15人回复
-
331求调剂
已经有3人回复
-
材料与化工专硕306分找合适调剂
已经有24人回复
-
材料工程调剂
已经有8人回复
-
070300化学279求调剂
已经有7人回复
-
一志愿上海海洋大学083200食品学硕，求调剂，接受其他专业
已经有12人回复
-
材料专硕调剂
已经有14人回复
-
334求调剂
已经有7人回复
-
材料工程085601，270求调剂
已经有10人回复
-
专硕310求调剂
已经有7人回复
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
