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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » pwscf 优化错误 diagonalization (DSYGV*) failed
121/2返回列表12下一页
查看: 1468  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 0XLacIJUOj8D 2026-08-24 8/400 2026-08-25 05:12 by BZKMTicpDhFj
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 0XLacIJUOj8D 2026-08-24 3/150 2026-08-25 01:47 by BZKMTicpDhFj
[基金申请] 明天应该可查了！？ +5 chengyan1220 2026-08-23 5/250 2026-08-24 23:28 by 我4大白菜
[基金申请] 没有任何消息-是不是就凉了 +7 图啦图啦 2026-08-24 8/400 2026-08-24 22:00 by maomao_da
[论文投稿] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 0XLacIJUOj8D 2026-08-24 3/150 2026-08-24 21:59 by BZKMTicpDhFj
[基金申请] 能否退出参与的面上项目解除限项 +21 koalala 2026-08-24 24/1200 2026-08-24 19:25 by 家与远方
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 今日不放榜？网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 21/1050 2026-08-24 14:21 by refreshing11
[基金申请] 范进中举一文的中心思想 +6 炎黄贵胄 2026-08-22 7/350 2026-08-24 11:58 by 6543yes
[基金申请] 让我中一个面上吧！ +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 17/850 2026-08-24 10:20 by echo8914667
[基金申请] 科研孤儿太难了 +18 我4大白菜 2026-08-20 19/950 2026-08-24 09:47 by 凯恩广盛催化分�
[基金申请] 今天放榜吗？ +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] 时间戳今天，20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
关闭
请填处理意见
关闭