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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » QE添加范德瓦尔兹力计算原子受力结果为NaN， 求教应该怎么办
51/1返回列表
查看: 1780  |  回复: 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个字符以上)
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料与化工，291分，求调剂 +7 咕噜咕噜123123 2026-03-04 7/350 2026-03-04 20:38 by kakakapanpan
[考研] 266求调剂 +7 哇塞王帅 2026-03-03 7/350 2026-03-04 20:30 by 一切OK
[考研] 一志愿武汉理工大学-085602-总分296分-求调剂 +6 紫川葡柚 2026-03-04 6/300 2026-03-04 20:10 by 1124861332
[基金申请] 成果系统访问量大，请15分钟后再尝试。由此给您造成的不便，敬请谅解。 +12 xhuama 2026-03-02 13/650 2026-03-04 16:29 by 69252822697
[考研] 材料328求调剂 +7 一个萝卜02 2026-03-03 7/350 2026-03-04 16:02 by zhukairuo
[考研] 281求调剂 +3 我是小小葱葱 2026-03-03 5/250 2026-03-04 14:23 by kakakapanpan
[考研] 289求调剂 +8 yang婷 2026-03-02 10/500 2026-03-04 09:10 by 无际的草原
[考研] 0856材料工程，初试313调剂 +7 卖个关子吧 2026-03-03 7/350 2026-03-03 21:33 by L135790
[考研] 化工专硕调剂 +4 利好利好. 2026-03-03 7/350 2026-03-03 21:30 by L135790
[考研] 理学，工学，农学调剂，少走弯路，这里欢迎您！ +8 likeihood 2026-03-02 11/550 2026-03-03 19:39 by 一叶凡尘?
[考研] 化工专硕348，一志愿985求调剂 +8 弗格个 2026-02-28 11/550 2026-03-03 14:55 by relum
[考研] 一志愿中科大能动297求调剂，本科川大 +3 邵11 2026-03-03 3/150 2026-03-03 14:50 by 新能源达人
[考研] 求调剂 +4 Guo_yuxuan 2026-03-02 5/250 2026-03-03 14:39 by xiaomc_gzh
[考研] 299求调剂 +5 kkcoco25 2026-03-02 9/450 2026-03-03 12:55 by 公瑾逍遥
[考研] 324求调剂 +4 wxz2 2026-03-03 5/250 2026-03-03 09:25 by 杨杨杨紫
[考研] 085600材料工程一志愿中科大总分312求调剂 +9 吃宵夜1 2026-02-28 11/550 2026-03-02 20:14 by hypershenger
[考研] 材料085601调剂 +5 多多子. 2026-03-02 5/250 2026-03-02 19:15 by zhukairuo
[基金申请] 此成果不能导入原因：元数据必填信息不完整，可 进行补充。 +4 Kittylucky 2026-03-02 5/250 2026-03-02 11:07 by jurkat.1640
[考研] 调剂 +3 13853210211 2026-03-02 4/200 2026-03-02 10:16 by 13853210211
[考研] 化工299分求调剂 一志愿985落榜 +5 嘻嘻(*^ω^*) 2026-03-01 5/250 2026-03-01 19:47 by 无际的草原
信息提示
关闭
请填处理意见
关闭