 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 1788 | 回復: 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文件發(fā)現(xiàn)： 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中的報錯，但不知為何會有這樣的錯誤（去掉范氏力參數(shù)則沒有報錯，但不能得到真實的晶格常數(shù)） 求各位大佬給出寶貴指導！ 以下是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 |
回復此樓» 猜你喜歡
一志愿中科大材料與化工。353分能調劑到哪？
已經有7人回復
-
337一志愿華南理工0805材料求調劑
已經有5人回復
-
伙伴們，祝我生日快樂吧
已經有4人回復
-
301求調劑
已經有3人回復
-
材料與化工304求B區(qū)調劑
已經有4人回復
-
讀博申請
已經有7人回復
-
一志愿江南大學085701環(huán)境工程專碩總分287求調劑
已經有5人回復
-
材料工程330分求調劑，一志愿985
已經有5人回復
-
0856化工原理
已經有4人回復
-
一志愿安徽大學材料工程專碩313分，求調劑的學校
已經有9人回復
» 本主題相關商家推薦: (我也要在這里推廣)
lq8989
至尊木蟲 (著名寫手)
- 應助: 4 (幼兒園)
- 金幣: 10104.1
- 散金: 169
- 紅花: 2
- 帖子: 1541
- 在線: 225.1小時
- 蟲號: 1452335
- 注冊: 2011-10-20
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
【答案】應助回帖
感謝參與，應助指數(shù) +1
|
QE開發(fā)團隊非常反對在高度對稱的體系里使用ibrav=0（其實最好永遠不要用ibrav=0），因為某些bug容易在ibrav=0時在對稱性上出現(xiàn)錯誤。對于這個簡單立方體系，應當在&system里使用 ibrav=0 和 A= 6.273，去掉nosym， 并刪除 CELL_PARAMETERS 這一塊兒數(shù)據。這樣可能問題就解決了。 如果還無法解決，就嘗試明確空間群，使用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
lq8989
至尊木蟲 (著名寫手)
- 應助: 4 (幼兒園)
- 金幣: 10104.1
- 散金: 169
- 紅花: 2
- 帖子: 1541
- 在線: 225.1小時
- 蟲號: 1452335
- 注冊: 2011-10-20
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
3樓2022-09-14 14:10:03
秋日以往
新蟲 (初入文壇)
- 應助: 0 (幼兒園)
- 金幣: 19.8
- 帖子: 18
- 在線: 3.9小時
- 蟲號: 8541249
- 注冊: 2018-04-16
- 性別: MM
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
4樓2022-11-28 15:32:49
相關版塊跳轉
我要訂閱樓主 秋日以往 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0856化工原理 +4 | z2839474511 2026-03-10 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 哈工大材料324求調劑 +6 | 閆旭東 2026-03-10 | 8/400 |
|
|
|
[考研] 歡迎211本科同學，過A區(qū)國家線，A區(qū)非偏遠一本，交叉學科課題組 +23 | lisimayy 2026-03-04 | 37/1850 |
|
|
|
[考研] 收調劑 +7 | 調劑的考研學生 2026-03-10 | 7/350 |
|
|
|
[考研] 0856材料與化工309分求調劑 +4 | ZyZy…… 2026-03-10 | 4/200 |
|
|
|
[考研] 一志愿：武漢理工，材料工程，英二數(shù)二 總分314 +3 | 2202020125 2026-03-10 | 4/200 |
|
|
[基金申請]
PDF在線壓縮軟件
 10+3
|
haxia 2026-03-08 | 4/200 |
|
|
|
[考博] 求材料讀博院校 +9 | yanglei131 2026-03-08 | 9/450 |
|
|
|
[考研] 復試調劑 +6 | 呼呼？~+123456 2026-03-08 | 8/400 |
|
|
|
[考研] 物理學求調劑 +4 | sx宋 2026-03-05 | 4/200 |
|
|
|
[碩博家園] 木蟲好像不熱鬧了，是不是？ +4 | 偏振片 2026-03-10 | 4/200 |
|
|
|
[考研] 085600材料與化工，一志愿廣州985，求調劑 +15 | qqyyaill 2026-03-05 | 15/750 |
|
|
|
[考研] 085701環(huán)境工程專碩求調劑 +5 | 幼稚園大班打手 2026-03-08 | 5/250 |
|
|
|
[考研] 288求調劑085600材料與化工 +13 | Daunrin 2026-03-07 | 15/750 |
|
|
|
[考研] 286求調劑 +12 | Faune 2026-03-06 | 14/700 |
|
|
|
[考研] 322分 085600求調劑，有互聯(lián)網+國金及主持省級大創(chuàng)經歷 +4 | 熊境喆 2026-03-05 | 4/200 |
|
|
|
[考研] 第一志愿上海大學，專業(yè)化學工程與技術，總分288，求調劑 +3 | 1829197082 2026-03-07 | 3/150 |
|
|
|
[考研] 誠求調劑，323分 有專利/科研/實習經歷 +6 | 橙子cyx 2026-03-06 | 8/400 |
|
|
|
[考研] 材料與化工304求B區(qū)調劑 +4 | 邱gl 2026-03-06 | 4/200 |
|
|
|
[考研] 304求調劑 +3 | 曼殊2266 2026-03-04 | 3/150 |
|