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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » QE添加范德瓦尔兹力计算原子受力结果为NaN， 求教应该怎么办
北京石油化工学院2026年研究生招生接收调剂公告
51/1返回列表
查看: 1808  |  回复: 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的回帖

秋日以往

新虫 (初入文坛)
引用回帖:
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 个回答

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

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 个回答
不应助 确定回帖应助 (注意：应助才可能被奖励，但不允许灌水，必须填写15个字符以上)
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0856求调剂 +5 楒桉 2026-03-28 5/250 2026-03-28 23:09 by 无际的草原
[考研] 321求调剂 +6 璞玉~~ 2026-03-25 7/350 2026-03-28 17:48 by 璞玉~~
[考研] 求调剂推荐 材料 304 +15 荷包蛋hyj 2026-03-26 15/750 2026-03-28 04:13 by fmesaito
[考研] 一志愿211院校 344分 东北农业大学生物学学硕，求调剂 +5 丶风雪夜归人丶 2026-03-26 8/400 2026-03-27 19:22 by 丶风雪夜归人丶
[考研] 一志愿上海理工能源动力（085800）310分求调剂 +3 zhangmingc 2026-03-27 4/200 2026-03-27 19:01 by 给你你注意休息
[考研] 一志愿 西北大学 总分282 英语一62 求调剂 +7 18419759900 2026-03-25 8/400 2026-03-27 16:38 by 18419759900
[考研] 085601 材料工程 313分 求调剂 +5 Ong3 2026-03-27 5/250 2026-03-27 12:24 by goldfish51
[考研] 314求调剂 +3 溪云珂 2026-03-26 3/150 2026-03-27 11:20 by sanrepian
[考研] 一志愿郑大085600，310分求调剂 +5 李潇可 2026-03-26 5/250 2026-03-27 11:14 by 不吃魚的貓
[考研] 324求调剂 +5 hanamiko 2026-03-26 5/250 2026-03-27 10:33 by wangjy2002
[考研] 349求调剂 +4 李木子啊哈哈 2026-03-25 4/200 2026-03-26 22:49 by fmesaito
[考研] 0703化学求调剂 +3 丹青奶盖 2026-03-26 5/250 2026-03-26 20:11 by macy2011
[考研] 281求调剂 +6 Koxui 2026-03-24 7/350 2026-03-26 15:37 by 无际的草原
[考研] 303求调剂 +6 蓝山月 2026-03-25 6/300 2026-03-25 22:47 by 418490947
[考研] 考研一志愿苏州大学初始315（英一）求调剂 +3 sbdksD 2026-03-24 4/200 2026-03-25 18:16 by xcjcqu
[考研] 296求调剂 +4 汪！？！ 2026-03-25 7/350 2026-03-25 16:41 by 汪！？！
[考研] 材料调剂 +3 iwinso 2026-03-23 3/150 2026-03-25 11:29 by greychen00
[考研] 材料考研调剂生 +3 黄粱一梦千年 2026-03-24 3/150 2026-03-24 17:00 by barlinike
[考研] 341求调剂(一志愿湖南大学070300) +5 番茄头--- 2026-03-22 6/300 2026-03-23 23:45 by Txy@872106
[考研] 336化工调剂 +4 王大坦1 2026-03-23 5/250 2026-03-23 18:32 by allen-yin
信息提示
关闭
请填处理意见
关闭