 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 1783 | 回復: 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 |
回復此樓» 猜你喜歡
歡迎211本科同學，過A區(qū)國家線，A區(qū)非偏遠一本，交叉學科課題組
已經(jīng)有24人回復
-
新疆大學地質與礦業(yè)工程學院招生
已經(jīng)有10人回復
-
招調(diào)劑
已經(jīng)有11人回復
-
0856材料與化工求調(diào)劑！
已經(jīng)有11人回復
-
材料與化工304求調(diào)劑
已經(jīng)有10人回復
-
A區(qū)一本交叉課題組，低分調(diào)劑，招收機械電子信息通信等交叉方向
已經(jīng)有31人回復
-
中科大材料299求調(diào)劑
已經(jīng)有7人回復
-
301求調(diào)劑
已經(jīng)有3人回復
-
材料考研339求調(diào)劑
已經(jīng)有3人回復
-
材料專碩290求調(diào)劑
已經(jīng)有5人回復
» 本主題相關商家推薦: (我也要在這里推廣)
lq8989
至尊木蟲 (著名寫手)
- 應助: 4 (幼兒園)
- 金幣: 10099.1
- 散金: 169
- 紅花: 2
- 帖子: 1540
- 在線: 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ù)據(jù)。這樣可能問題就解決了。 如果還無法解決，就嘗試明確空間群，使用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 (幼兒園)
- 金幣: 10099.1
- 散金: 169
- 紅花: 2
- 帖子: 1540
- 在線: 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ā)表 | |
|---|---|---|---|---|
|
[考研] 材料與化工304求調(diào)劑 +7 | 邱gl 2026-03-05 | 10/500 |
|
|---|---|---|---|---|
|
|
[考研] 304分材料專碩求調(diào)劑 +11 | qiuzhigril 2026-03-03 | 14/700 |
|
|
[考研] 材料復試調(diào)劑 +7 | 學材料的點 2026-03-01 | 8/400 |
|
|
|
[考研] 一志愿武漢理工大學-085602-總分296分-求調(diào)劑 +7 | 紫川葡柚 2026-03-04 | 7/350 |
|
|
|
[考研] 264求調(diào)劑 +8 | 26調(diào)劑 2026-03-03 | 8/400 |
|
|
|
[考研] 267化工調(diào)劑求助 +5 | 聰少OZ 2026-03-04 | 5/250 |
|
|
|
[考研] 085600 材料與化工 298 +14 | 小西笑嘻嘻 2026-03-03 | 14/700 |
|
|
|
[考研] 一志愿314求調(diào)劑 +7 | 202111120625 2026-03-03 | 7/350 |
|
|
|
[考研] 281求調(diào)劑 +3 | 我是小小蔥蔥 2026-03-03 | 5/250 |
|
|
|
[考研] 0856材料工程，初試313調(diào)劑 +7 | 賣個關子吧 2026-03-03 | 7/350 |
|
|
|
[考研] 292求調(diào)劑 +3 | sgbl 2026-03-03 | 3/150 |
|
|
|
[考研] 一志愿華中科技大學，化學專業(yè)344分，求調(diào)劑 +6 | 邢xing1 2026-03-02 | 6/300 |
|
|
|
[考研] 求調(diào)劑 +4 | Guo_yuxuan 2026-03-02 | 5/250 |
|
|
|
[考研] 材料工程求調(diào)劑 +3 | 1431251 2026-03-03 | 3/150 |
|
|
|
[考研] 288求調(diào)劑 +3 | 少71.8 2026-03-02 | 5/250 |
|
|
|
[考研] 0854復試調(diào)劑 276 +5 | wmm9 2026-03-01 | 7/350 |
|
|
|
[考研] 一志愿東北大學材料專碩328，求調(diào)劑 +3 | shs1083 2026-03-02 | 3/150 |
|
|
|
[考研] 化工專碩342，一志愿大連理工大學，求調(diào)劑 +6 | kyf化工 2026-02-28 | 7/350 |
|
|
|
[考研] 材料學調(diào)劑 +10 | 提神豆沙包 2026-02-28 | 12/600 |
|
|
|
[考研] 311求調(diào)劑 +6 | 亭亭亭01 2026-03-01 | 6/300 |
|