24小时热门版块排行榜

返回列表

【悬赏金币】回答本帖问题，作者秋日以往将赠送您 5 个金币

秋日以往

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 19.8
帖子: 18
在线: 3.9小时
虫号: 8541249
注册: 2018-04-16
性别: MM
专业: 凝聚态物性 II ：电子结构

[求助] 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

回复此楼

» 猜你喜欢

请问哪里可以有青B申请的本子可以借鉴一下。已经有4人回复
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗已经有6人回复
孩子确诊有中度注意力缺陷已经有14人回复
三甲基碘化亚砜的氧化反应已经有4人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有5人回复
2025冷门绝学什么时候出结果已经有3人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复
AI论文写作工具：是科研加速器还是学术作弊器？已经有3人回复

» 本主题相关商家推荐: (我也要在这里推广)

1楼2022-09-14 12:20:51

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

lq8989

铁杆木虫 (著名写手)

应助: 4 (幼儿园)
金币: 9894.1
散金: 169
红花: 2
帖子: 1499
在线: 224.9小时
虫号: 1452335
注册: 2011-10-20
专业: 凝聚态物性 II ：电子结构

【答案】应助回帖

感谢参与，应助指数 +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

铁杆木虫 (著名写手)

应助: 4 (幼儿园)
金币: 9894.1
散金: 169
红花: 2
帖子: 1499
在线: 224.9小时
虫号: 1452335
注册: 2011-10-20
专业: 凝聚态物性 II ：电子结构

引用回帖:

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