 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 1782 | 回復(fù): 3 | |||
| 【懸賞金幣】回答本帖問題，作者秋日以往將贈(zèng)送您 5 個(gè)金幣 | |||
秋日以往新蟲 (初入文壇)
|
[求助]
QE添加范德瓦爾茲力計(jì)算原子受力結(jié)果為NaN， 求教應(yīng)該怎么辦 已有1人參與
|
||
|
目的是弛豫獲得CsPbI3的立方相晶體結(jié)構(gòu) 材料的范德瓦爾茲力比較強(qiáng)，對(duì)晶格影響比較大，在qe中選的是DFT-d3的方法 計(jì)算結(jié)果報(bào)錯(cuò)，CRASH文件里這樣顯示： %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% task # 0 from checkallsym : error # 1 some of the original symmetry operations not satisfied %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 嘗試關(guān)閉對(duì)稱性（ nosym = .true.）之后仍然有相同的報(bào)錯(cuò)，檢查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 猜測(cè)是由于受力沒有成功計(jì)算才導(dǎo)致CRASH中的報(bào)錯(cuò)，但不知為何會(huì)有這樣的錯(cuò)誤（去掉范氏力參數(shù)則沒有報(bào)錯(cuò)，但不能得到真實(shí)的晶格常數(shù)） 求各位大佬給出寶貴指導(dǎo)！ 以下是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 |
回復(fù)此樓» 猜你喜歡
歡迎211本科同學(xué)，過A區(qū)國家線，A區(qū)非偏遠(yuǎn)一本，交叉學(xué)科課題組
已經(jīng)有12人回復(fù)
-
成果系統(tǒng)訪問量大，請(qǐng)15分鐘后再嘗試。由此給您造成的不便，敬請(qǐng)諒解。
已經(jīng)有14人回復(fù)
-
281電子信息求調(diào)劑
已經(jīng)有8人回復(fù)
-
材料277分求調(diào)劑
已經(jīng)有7人回復(fù)
-
求調(diào)劑院校
已經(jīng)有12人回復(fù)
-
293求調(diào)劑
已經(jīng)有3人回復(fù)
-
化學(xué) 0703求調(diào)劑 總分293 一志愿211
已經(jīng)有10人回復(fù)
-
320材料一志愿西工大專碩英二數(shù)二 有兩年光伏工作經(jīng)驗(yàn)
已經(jīng)有7人回復(fù)
-
材料328求調(diào)劑
已經(jīng)有8人回復(fù)
-
304分材料專碩求調(diào)劑
已經(jīng)有12人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
lq8989
至尊木蟲 (著名寫手)
- 應(yīng)助: 4 (幼兒園)
- 金幣: 10099.1
- 散金: 169
- 紅花: 2
- 帖子: 1540
- 在線: 225.1小時(shí)
- 蟲號(hào): 1452335
- 注冊(cè): 2011-10-20
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
【答案】應(yīng)助回帖
感謝參與，應(yīng)助指數(shù) +1
|
QE開發(fā)團(tuán)隊(duì)非常反對(duì)在高度對(duì)稱的體系里使用ibrav=0（其實(shí)最好永遠(yuǎn)不要用ibrav=0），因?yàn)槟承゜ug容易在ibrav=0時(shí)在對(duì)稱性上出現(xiàn)錯(cuò)誤。對(duì)于這個(gè)簡(jiǎn)單立方體系，應(yīng)當(dāng)在&system里使用 ibrav=0 和 A= 6.273，去掉nosym， 并刪除 CELL_PARAMETERS 這一塊兒數(shù)據(jù)。這樣可能問題就解決了。 如果還無法解決，就嘗試明確空間群，使用space_group = 221， 結(jié)合 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 只輸入不等價(jià)原子，來完全保證對(duì)稱性正確。 |
2樓2022-09-14 14:08:49
lq8989
至尊木蟲 (著名寫手)
- 應(yīng)助: 4 (幼兒園)
- 金幣: 10099.1
- 散金: 169
- 紅花: 2
- 帖子: 1540
- 在線: 225.1小時(shí)
- 蟲號(hào): 1452335
- 注冊(cè): 2011-10-20
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
3樓2022-09-14 14:10:03
秋日以往
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 19.8
- 帖子: 18
- 在線: 3.9小時(shí)
- 蟲號(hào): 8541249
- 注冊(cè): 2018-04-16
- 性別: MM
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
4樓2022-11-28 15:32:49
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 秋日以往 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑院校 +6 | 云朵452 2026-03-02 | 12/600 |
|
|---|---|---|---|---|
|
|
[考研] 320材料一志愿西工大專碩英二數(shù)二 有兩年光伏工作經(jīng)驗(yàn) +4 | 錘某人 2026-03-04 | 7/350 |
|
|
|
[考研] 材料328求調(diào)劑 +8 | 一個(gè)蘿卜02 2026-03-03 | 8/400 |
|
|
[考研] 304分材料專碩求調(diào)劑 +9 | qiuzhigril 2026-03-03 | 12/600 |
|
|
|
[考研] 286 +6 | ksncj 2026-03-04 | 6/300 |
|
|
|
[考研] 302材料工程求調(diào)劑 +6 | Doleres 2026-03-01 | 7/350 |
|
|
|
[考研]
 一志愿中南大學(xué)理學(xué)化學(xué)
+6
|
15779376950 2026-03-01 | 8/400 |
|
|
|
[考研] 歡迎采礦、地質(zhì)、巖土、計(jì)算機(jī)、人工智能等專業(yè)的同學(xué)報(bào)考 +9 | pin8023 2026-02-28 | 12/600 |
|
|
|
[考研] 347求調(diào)劑 +5 | 啊歐歐歐 2026-03-03 | 5/250 |
|
|
|
[碩博家園] 2025屆雙非化工碩士畢業(yè)，申博 +4 | 更多的是 2026-02-27 | 5/250 |
|
|
|
[考研] 299求調(diào)劑 +5 | kkcoco25 2026-03-02 | 9/450 |
|
|
|
[考研] 26考研報(bào)考西工大材料308分求調(diào)劑 +4 | weizhong123 2026-03-01 | 5/250 |
|
|
|
[考研] 考研復(fù)試調(diào)劑，過國家線的同學(xué)都可報(bào)名 +7 | 黑！在干嘛 2026-02-28 | 8/400 |
|
|
|
[考研] 321求調(diào)劑一志愿東北林業(yè)大學(xué)材料與化工英二數(shù)二 +5 | 蟲蟲蟲蟲蟲7 2026-03-01 | 9/450 |
|
|
|
[考研] 化學(xué)，材料，環(huán)境類求調(diào)劑 +7 | 考研版棒棒 2026-03-02 | 7/350 |
|
|
|
[考研] 材料085601調(diào)劑 +5 | 多多子. 2026-03-02 | 5/250 |
|
|
|
[考研] 材料調(diào)劑 +6 | 愛擦汗的可樂冰 2026-02-28 | 7/350 |
|
|
|
[考研] 279求調(diào)劑 +3 | dua1 2026-03-01 | 4/200 |
|
|
|
[考研] 化工299分求調(diào)劑 一志愿985落榜 +5 | 嘻嘻(*^ω^*) 2026-03-01 | 5/250 |
|
|
|
[考研]
尋找調(diào)劑
+4
|
LYidhsjabdj 2026-02-28 | 4/200 |
|