应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » pwscf 优化错误 diagonalization (DSYGV*) failed
121/2返回列表12下一页
查看: 1156  |  回复: 11
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

minmin_0082003

金虫 (正式写手)

莫愁莫愁

[求助] pwscf 优化错误 diagonalization (DSYGV*) failed

本来想重复一篇文献的数据，但是计算的时候出现错误
“ from rdiaghg : error #       201
     diagonalization (DSYGV*) failed”
该如何解决，去mail-list 上看了，解决办法也试了，都不行。
现附上输入文件，赝式用的是官网上的PBE超软赝式
“&CONTROL
  calculation    = "relax",
  nstep          = 500,
  prefix       = "c6h5nc-pt"
  pseudo_dir   = "/home/gaomin/PWscf/freq-PP/GGA-uspp",
  outdir       = "/scr/gaomin",
/
&SYSTEM
  ibrav                  = 0,
  nat                    = 22,
  ntyp                   = 4,
  ecutwfc                = 50.0D0,
  ecutrho                = 320.0D0,
  nspin                  = 2,
  starting_magnetization(1) = 0.5,
  starting_magnetization(2) = 0.5,
  starting_magnetization(3) = 0.5,
  starting_magnetization(4) = 0.5,
  occupations              = "smearing"
  degauss                  = 0.01
/
&ELECTRONS
  conv_thr    = 1.0D-5,
  mixing_beta = 0.3D0,
  electron_maxstep = 400,
/
&IONS
  ds                = 1.D0,
/
ATOMIC_SPECIES
H  1.00794   H.pbe.UPF
C  12.0107   C.pbe.UPF
N  14.00674  N.pbe.UPF
Au 196.96655    Au.pbe.UPF
ATOMIC_POSITIONS { angstrom }
Au                 0.00000000    0.00000000    2.26546472
Au                 2.77461630    0.00000000    2.26546472
Au                 5.54923260    0.00000000    2.26546472
Au                -1.38730815    2.40288820    2.26546472
Au                 1.38730815    2.40288820    2.26546472
Au                 4.16192445    2.40288820    2.26546472
Au                -2.77461630    4.80577640    2.26546472
Au                 0.00000000    4.80577640    2.26546472
Au                 2.77461630    4.80577640    2.26546472
   C   0.349789684       0.345710403       0.154967780
   H  0.364068615       0.357588915       0.199276427
   C  0.557854327       0.394622547       0.270754759
   H  0.662499028       0.393210459       0.245646609
   C  0.584252875       0.413015632       0.322598632
   H  0.712994743       0.425952068       0.338539077
   C  0.450123189       0.414694113       0.354832681
   C  0.286864953       0.397757497       0.333230714
   H  0.179802696       0.398486518       0.357573855
   C  0.256933266       0.379098849       0.281544186
   H  0.129306327       0.365454617       0.264814226
   C  0.392824094       0.377111719       0.250101757
   H  0.471654467       0.428447475       0.394606202
CELL_PARAMETERS {cubic}
                 8.32380000    0.00000000    0.00000000
                  -4.16190000    7.2084108    0.00000000
                  0.00000000    0.00000000    18.86460000
K_POINTS gamma
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2012-02-27 11:03:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)

【答案】应助回帖

感谢参与，应助指数 +1
minmin_0082003(金币+5): ★★★很有帮助 找到错误了，结构有问题，哈哈， 2012-02-27 11:24:11
franch(金币+5, 专家考核): 鼓励交流。。呵呵 2012-02-27 13:38:50
ecut比较怪异，怎么是6倍？ecutrho可以取大一些，比如10倍。另外，你用的是什么类型赝势？模守恒吗？那就可以把ecutwfc可以再大一些。这个问题根本上是能量截断，也就是平面波波数不够多引起的。
赞一下 回复此楼
2楼2012-02-27 11:19:48
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

minmin_0082003

金虫 (正式写手)

莫愁莫愁
引用回帖:
2楼: Originally posted by goldenfisher at 2012-02-27 11:19:48:
ecut比较怪异，怎么是6倍？ecutrho可以取大一些，比如10倍。另外，你用的是什么类型赝势？模守恒吗？那就可以把ecutwfc可以再大一些。这个问题根本上是能量截断，也就是平面波波数不够多引起的。

我用的是超软赝式，模守恒估计特别的慢吧
赞一下 回复此楼
3楼2012-02-27 11:28:53
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)
WDD880227(金币+1): 感谢交流提示 2012-02-29 13:57:48
找到错误就好，不过还是应该提示一下，超软的ecutrho至少应该是ecutwfc的8倍
赞一下 回复此楼
4楼2012-02-28 21:05:09
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

minmin_0082003

金虫 (正式写手)

莫愁莫愁
引用回帖:
4楼: Originally posted by goldenfisher at 2012-02-28 21:05:09:
找到错误就好，不过还是应该提示一下，超软的ecutrho至少应该是ecutwfc的8倍

灰常感谢，一般8-10倍。越大越好么
赞一下 回复此楼
5楼2012-02-29 09:33:41
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)
fzx2008(金币+2, 专家考核): 谢谢指教 2012-02-29 10:11:24
引用回帖:
5楼: Originally posted by minmin_0082003 at 2012-02-29 09:33:41:
灰常感谢，一般8-10倍。越大越好么

只增ecutrho不增ecutwfc，超过12倍后一般认为意义不大了。
除了cutoff的因素，偶尔空带数目的多寡也能影响到出现对角化矩阵操作出错。可以适当增加一下nbnd或者降低一下（不要少于不使用该参数时的系统确定的nbnd数就可）。另外，如果用davidson方法求解出错，可以尝试用CG方法求解，CG的成功率是一般认为大于davidson的，不过效率也低很多……
赞一下 回复此楼
思想重于技巧，内涵重于表象
6楼2012-02-29 10:10:05
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)
顺便补一张官方的解释截图：
赞一下 回复此楼
思想重于技巧，内涵重于表象
7楼2012-02-29 10:11:38
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

minmin_0082003

金虫 (正式写手)

莫愁莫愁
引用回帖:
7楼: Originally posted by souledge at 2012-02-29 10:11:38:
顺便补一张官方的解释截图：

图呢
回复此楼
8楼2012-02-29 10:19:49
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)
liliangfang(金币+1): 多谢交流 2012-02-29 14:56:32
引用回帖:
8楼: Originally posted by minmin_0082003 at 2012-02-29 10:19:49:
图呢

俺也不知道，就那么没了…… 就直接到PWGUI里看吧，那里有说明～或者Doc/INPUT_PW也一样～
赞一下 回复此楼
思想重于技巧，内涵重于表象
9楼2012-02-29 13:49:38
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

minmin_0082003

金虫 (正式写手)

莫愁莫愁
引用回帖:
9楼: Originally posted by souledge at 2012-02-29 13:49:38:
俺也不知道，就那么没了…… 就直接到PWGUI里看吧，那里有说明～或者Doc/INPUT_PW也一样～

好吧，谢谢
回复此楼
10楼2012-02-29 14:45:13
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 minmin_0082003 的主题更新
121/2返回列表12下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭