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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » pwscf 优化错误 diagonalization (DSYGV*) failed
121/2返回列表12下一页
查看: 1498  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的？ +15 爱看书的可乐 2026-08-26 17/850 2026-08-28 01:23 by xxniao123
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 XLGfUIdOmEAO 2026-08-27 3/150 2026-08-28 00:49 by 3RAFyYXYSwGS
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[教师之家] 导师吐槽：我怎么摊上了这么个极品研究生！ +7 苏东坡二世 2026-08-23 7/350 2026-08-27 18:14 by 瞬息宇宙
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 基金不中，共勉 +10 eulota 2026-08-26 10/500 2026-08-27 17:23 by lqllinqiaoli
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜？ 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合现在查不到了吗？ +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 牛来！米来！面来！ +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 明天应该可查了！？ +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
信息提示
关闭
请填处理意见
关闭