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 ]

回复此楼

» 猜你喜欢

面上可以超过30页吧？已经有13人回复
网上报道青年教师午睡中猝死、熬夜猝死的越来越多，主要哪些原因引起的？已经有10人回复
为什么中国大学工科教授们水了那么多所谓的顶会顶刊，但还是做不出宇树机器人？已经有13人回复
什么是人一生最重要的？已经有8人回复
版面费该交吗已经有17人回复
体制内长辈说体制内绝大部分一辈子在底层，如同你们一样大部分普通教师忙且收入低已经有19人回复
【博士招生】太原理工大学2026化工博士已经有8人回复
280求调剂已经有4人回复

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

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