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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » QE添加范德瓦尔兹力计算原子受力结果为NaN， 求教应该怎么办
51/1返回列表
查看: 1779  |  回复: 3
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【悬赏金币】回答本帖问题，作者秋日以往将赠送您 5 个金币
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

秋日以往

新虫 (初入文坛)

[求助] QE添加范德瓦尔兹力计算原子受力结果为NaN， 求教应该怎么办 已有1人参与

目的是弛豫获得CsPbI3的立方相晶体结构
材料的范德瓦尔兹力比较强，对晶格影响比较大，在qe中选的是DFT-d3的方法
计算结果报错，CRASH文件里这样显示：
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
     task #         0
     from checkallsym : error #         1
     some of the original symmetry operations not satisfied
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

尝试关闭对称性（ nosym = .true.）之后仍然有相同的报错，检查out文件发现：

     atom    1 type  1   force =            NaN           NaN           NaN
     atom    2 type  2   force =            NaN           NaN           NaN
     atom    3 type  3   force =            NaN           NaN           NaN
     atom    4 type  3   force =            NaN           NaN           NaN
     atom    5 type  3   force =            NaN           NaN           NaN
     atom    1 type  1   force =            NaN           NaN           NaN
     atom    2 type  2   force =            NaN           NaN           NaN
     atom    3 type  3   force =            NaN           NaN           NaN
     atom    4 type  3   force =            NaN           NaN           NaN
     atom    5 type  3   force =            NaN           NaN           NaN
     Total force =          NaN     Total SCF correction =     0.000002

猜测是由于受力没有成功计算才导致CRASH中的报错，但不知为何会有这样的错误（去掉范氏力参数则没有报错，但不能得到真实的晶格常数）

求各位大佬给出宝贵指导！

以下是in文件：
&CONTROL
    calculation='vc-relax', prefix='CsPbI3',
    outdir='./', verbosity='high'
    pseudo_dir = '/home/zhaizh/2-CsPbI3/qe/CON_PBE_rel/pseudo',
    tprnfor=.true., tstress=.true., forc_conv_thr=1.0d-4
/
&SYSTEM
    ibrav= 0,
    nat= 5, ntyp= 3,
    ecutwfc= 60, ecutrho = 300,
!   noncolin = .true.
!   lspinorb=.true.
    vdw_corr = 'DFT-D3'
    dftd3_version = 4   !D3 with BJ damping
    nosym= .true.

/
&ELECTRONS
    electron_maxstep = 500
    conv_thr = 1.0d-12
/
&IONS
    ion_dynamics='bfgs'
/
&CELL
    press_conv_thr=0.1
/


&IONS
    ion_dynamics='bfgs'
/
&CELL
    press_conv_thr=0.1
/



ATOMIC_SPECIES
    Cs  132.90545   Cs_ONCV_PBE_FR-1.0.upf
    Pb  207.2       Pb_ONCV_PBE_FR-1.0.upf
    I   126.90447   I_ONCV_PBE_FR-1.1.upf

ATOMIC_POSITIONS {crystal}
    Cs      0.5000000000   0.500000   0.500000
    Pb      0.0000000000   0.000000   0.000000
    I       0.5000000000   0.000000   0.000000
    I       0.0000000000   0.500000   0.000000
    I       0.0000000000   0.000000   0.500000

CELL_PARAMETERS {angstrom}
   6.273000000   0.000000000   0.000000000
   0.000000000   6.273000000   0.000000000
   0.000000000   0.000000000   6.273000000

K_POINTS {automatic}
    11 11 11 0 0 0
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)
1楼 2022-09-14 12:20:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lq8989

至尊木虫 (著名写手)
引用回帖:
2楼: Originally posted by lq8989 at 2022-09-14 14:08:49
QE开发团队非常反对在高度对称的体系里使用ibrav=0（其实最好永远不要用ibrav=0），因为某些bug容易在ibrav=0时在对称性上出现错误。对于这个简单立方体系，应当在&system里使用 ibrav=0 和 A= 6.273，去掉 ...

笔误，应当是ibrav=1 和 A= 6.273
赞一下 回复此楼
3楼2022-09-14 14:10:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

lq8989

至尊木虫 (著名写手)

【答案】应助回帖

感谢参与，应助指数 +1
QE开发团队非常反对在高度对称的体系里使用ibrav=0（其实最好永远不要用ibrav=0），因为某些bug容易在ibrav=0时在对称性上出现错误。对于这个简单立方体系，应当在&system里使用 ibrav=0 和 A= 6.273，去掉nosym， 并删除 CELL_PARAMETERS 这一块儿数据。这样可能问题就解决了。

如果还无法解决，就尝试明确空间群，使用space_group = 221， 结合 nat=3 以及

ATOMIC_POSITIONS crystal_sg
    Cs      0.5000000000   0.500000   0.500000
    Pb      0.0000000000   0.000000   0.000000
    I       0.5000000000   0.000000   0.000000

只输入不等价原子，来完全保证对称性正确。
赞一下(1人) 回复此楼
2楼2022-09-14 14:08:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

秋日以往

新虫 (初入文坛)
引用回帖:
2楼: Originally posted by lq8989 at 2022-09-14 14:08:49
QE开发团队非常反对在高度对称的体系里使用ibrav=0（其实最好永远不要用ibrav=0），因为某些bug容易在ibrav=0时在对称性上出现错误。对于这个简单立方体系，应当在&system里使用 ibrav=0 和 A= 6.273，去掉 ...

谢谢！
回复此楼
4楼2022-11-28 15:32:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答
不应助 确定回帖应助 (注意：应助才可能被奖励，但不允许灌水，必须填写15个字符以上)
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 环境工程专硕307 一志愿211 四六级已过 求调剂 +5 ccc! 2026-03-03 5/250 2026-03-04 17:39 by ZYNG
[考研] 求调剂 +4 熬夜的猫头鹰 2026-03-02 4/200 2026-03-04 17:22 by 梦天888
[考研] 材料复试调剂 +6 学材料的点 2026-03-01 7/350 2026-03-04 17:21 by 梦天888
[考研] 085701环境工程 求调剂 +6 xiiiia 2026-03-04 6/300 2026-03-04 16:45 by yuningshan
[考研] 化学 0703求调剂 总分293 一志愿211 +8 土土小虫 2026-03-03 8/400 2026-03-04 16:45 by xwxstudy
[考研] 295求调剂 +4 小赛不吃香菜 2026-03-04 4/200 2026-03-04 16:13 by sslc1985
[考研] 295求调剂 +6 等春来， 2026-03-04 6/300 2026-03-04 16:11 by zhukairuo
[考研] 085601 材料305分求助 +4 泡泡邮件 2026-03-03 6/300 2026-03-04 13:00 by 无际的草原
[考研] 322,求调剂 +3 菜菜爱玩 2026-03-04 3/150 2026-03-04 12:15 by xiongyaxuan
[考研] 材料工程269求调剂 +7 白刺玫 2026-03-02 7/350 2026-03-04 08:47 by zhyzzh
[基金申请] 请问大家，研究风险与应对措施那里， 大家都怎么写呢 ？ +3 cauasen 2026-03-02 3/150 2026-03-03 19:43 by Tide man
[考研] 085602化学工程350，调剂，有没有211的 +5 利好利好. 2026-03-02 9/450 2026-03-03 17:06 by 利好利好.
[考研] 材料类考研调剂 +6 gemmgemm 2026-03-01 7/350 2026-03-03 11:16 by EBSD
[论文投稿] 通讯作者写谁，问题是你意想不到的问题 15+3 阿尔法啊 2026-03-01 3/150 2026-03-03 09:13 by 北京莱茵润色
[考研] 一志愿中石油（华东）本科齐鲁工业大学 +3 石能伟 2026-03-02 3/150 2026-03-02 18:54 by caszguilin
[基金申请] 此成果不能导入原因：元数据必填信息不完整，可 进行补充。 +4 Kittylucky 2026-03-02 5/250 2026-03-02 11:07 by jurkat.1640
[考研] 275求调剂 +3 L-xin? 2026-03-01 6/300 2026-03-02 10:22 by 热情沙漠
[基金申请] 本子写完了，给DS兄弟看了，得了92分 +3 Doma 2026-03-01 7/350 2026-03-02 00:00 by jnzsy
[考研] 化工299分求调剂 一志愿985落榜 +5 嘻嘻(*^ω^*) 2026-03-01 5/250 2026-03-01 19:47 by 无际的草原
[考研] 307求调剂 +4 73372112 2026-02-28 6/300 2026-03-01 00:04 by ll247
信息提示
关闭
请填处理意见
关闭