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 ]

回复此楼

» 猜你喜欢

博士读完未来一定会好吗已经有6人回复
小论文投稿已经有3人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有9人回复
心脉受损已经有3人回复
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有8人回复
申请2026年博士已经有6人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有5人回复

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

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的回帖