| 查看: 1805 | 回复: 3 | ||
| 【悬赏金币】回答本帖问题，作者秋日以往将赠送您 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 |
回复此楼» 猜你喜欢
311（085601）求调剂
已经有4人回复
-
305求调剂
已经有4人回复
-
一志愿北化085600材料专硕275|有文章专利｜求调剂
已经有7人回复
-
一志愿华理，数一英一285求A区调剂
已经有12人回复
-
289求调剂
已经有11人回复
-
食品工程专硕一志愿中海洋309求调剂
已经有7人回复
-
生物学学硕，一志愿湖南大学，初试成绩338
已经有6人回复
-
321求调剂
已经有7人回复
-
343求调剂
已经有4人回复
-
本科新能源科学与工程，一志愿华理能动285求调剂
已经有6人回复
» 本主题相关商家推荐: (我也要在这里推广)
lq8989
至尊木虫 (著名写手)
- 应助: 4 (幼儿园)
- 金币: 10139.1
- 散金: 169
- 红花: 2
- 帖子: 1548
- 在线: 225.1小时
- 虫号: 1452335
- 注册: 2011-10-20
- 专业: 凝聚态物性 II ：电子结构
3楼2022-09-14 14:10:03
lq8989
至尊木虫 (著名写手)
- 应助: 4 (幼儿园)
- 金币: 10139.1
- 散金: 169
- 红花: 2
- 帖子: 1548
- 在线: 225.1小时
- 虫号: 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 只输入不等价原子，来完全保证对称性正确。 |
2楼2022-09-14 14:08:49
秋日以往
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 19.8
- 帖子: 18
- 在线: 3.9小时
- 虫号: 8541249
- 注册: 2018-04-16
- 性别: MM
- 专业: 凝聚态物性 II ：电子结构
4楼2022-11-28 15:32:49
