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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » QE添加范德瓦尔兹力计算原子受力结果为NaN， 求教应该怎么办
51/1返回列表
查看: 1778  |  回复: 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个字符以上)
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085600材料与化工（高分子）290分求调剂 +4 wengyujian 2026-03-04 4/200 2026-03-04 18:10 by houyaoxu
[考研] 学硕材料275调剂 +9 路三三 2026-03-03 9/450 2026-03-04 17:02 by 梦天888
[考研] 085600材料调剂 总分330 +3 池池丶 2026-03-03 3/150 2026-03-04 16:47 by yuningshan
[考研] 085600 材料与化工 298 +14 小西笑嘻嘻 2026-03-03 14/700 2026-03-04 16:23 by sslc1985
[考研] 环境调剂 +6 柒槿levana 2026-03-01 6/300 2026-03-04 15:55 by hypershenger
[考研] 材料化工调剂 +15 今夏不夏 2026-03-01 18/900 2026-03-04 15:44 by 每天只摆一小会
[考研] 322,求调剂 +3 菜菜爱玩 2026-03-04 3/150 2026-03-04 12:15 by xiongyaxuan
[考研] 291求调剂 +3 咕噜咕噜123123 2026-03-02 4/200 2026-03-04 08:39 by 每天只摆一小会
[考研] 292求调剂 +3 sgbl 2026-03-03 3/150 2026-03-03 17:30 by barlinike
[考研] 0857调剂 +6 一ll半 2026-02-28 8/400 2026-03-03 09:30 by 热情沙漠
[考研] 课题组接收材料类调剂研究生 +6 gaoxiaoniuma 2026-02-28 9/450 2026-03-03 09:27 by gaoxiaoniuma
[考研] 清华大学 材料与化工 353分求调剂 +5 awaystay 2026-03-02 6/300 2026-03-03 09:03 by houyaoxu
[考研] 高分子化学与物理调剂 +6 好好好1233 2026-02-28 16/800 2026-03-03 00:36 by 好好好1233
[考研] 285求调剂 +9 满头大汗的学生 2026-02-28 9/450 2026-03-02 20:29 by hypershenger
[考研] 261求调剂 +3 陆lh 2026-03-01 3/150 2026-03-02 19:32 by zhukairuo
[考研] 一志愿山东大学材料与化工325求调剂 +5 半截的诗0927 2026-03-02 5/250 2026-03-02 18:37 by 明亮9527
[考研] 0856化工专硕求调剂 +15 董boxing 2026-03-01 15/750 2026-03-02 15:06 by 晃晃不许晃
[考研] 322求调剂 +3 熊境喆 2026-03-01 3/150 2026-03-02 08:44 by houyaoxu
[考研] 调剂 +3 简木ChuFront 2026-02-28 3/150 2026-03-01 11:46 by 王伟要上岸啊
[论文投稿] Optics letters投稿被拒求助 30+3 luckyry 2026-02-26 4/200 2026-03-01 09:06 by babero
信息提示
关闭
请填处理意见
关闭