24小时热门版块排行榜

返回列表

当前主题已经存档。

zxzj05

荣誉版主 (著名写手)

1ST强帖: 1
应助: 0 (幼儿园)
贵宾: 1.242
金币: 5912.5
红花: 2
帖子: 1528
在线: 144.5小时
虫号: 406147
注册: 2007-06-18
性别: GG
专业: 金属功能材料
管辖: 第一性原理

[交流] 【求助】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 ]

回复此楼

» 猜你喜欢

有时候真觉得大城市人没有县城人甚至个体户幸福已经有9人回复
CSC & MSCA 博洛尼亚大学能源材料课题组博士/博士后招生｜MSCA经费充足、排名优已经有6人回复
天津大学招2026.09的博士生，欢迎大家推荐交流（博导是本人）已经有4人回复
同年申请2项不同项目，第1个项目里不写第2个项目的信息，可以吗已经有3人回复
退学或坚持读已经有28人回复
面上项目申报已经有3人回复
酰胺脱乙酰基已经有9人回复
博士延得我，科研能力直往上蹿已经有7人回复
面上基金申报没有其他的参与者成吗已经有5人回复
遇见不省心的家人很难过已经有22人回复

储氢家族欢迎储氢研究者！

1楼 2008-12-04 22:17:55

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

watermall

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 1981
散金: 18
红花: 4
帖子: 286
在线: 295.6小时
虫号: 118944
注册: 2005-11-27
专业: 固体无机化学

★
wuchenwf(金币+1):3Q！

建议把使用的abinit版本
输入文件
输出文件，尤其是出错的地方给出来
大家才能更清楚你的问题是什么

赞一下(1人)

回复此楼

2楼2008-12-04 22:23:02

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖