 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 1775 | 回復(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ù)此樓» 猜你喜歡
264求調(diào)劑
已經(jīng)有7人回復(fù)
-
085600 材料與化工 298
已經(jīng)有8人回復(fù)
-
266求調(diào)劑
已經(jīng)有5人回復(fù)
-
0854總分272
已經(jīng)有6人回復(fù)
-
EST拒稿重投
已經(jīng)有3人回復(fù)
-
歡迎采礦、地質(zhì)、巖土、計(jì)算機(jī)、人工智能等專業(yè)的同學(xué)報(bào)考
已經(jīng)有11人回復(fù)
-
0703化學(xué)學(xué)碩 288分
已經(jīng)有4人回復(fù)
-
接收調(diào)劑
已經(jīng)有23人回復(fù)
-
化學(xué) 0703求調(diào)劑 總分293 一志愿211
已經(jīng)有3人回復(fù)
-
304分材料專碩求調(diào)劑
已經(jīng)有5人回復(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ù)。這樣可能問題就解決了。 如果還無(wú)法解決，就嘗試明確空間群，使用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à)原子，來(lái)完全保證對(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
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 266求調(diào)劑 +5 | 哇塞王帥 2026-03-03 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 接收調(diào)劑 +15 | 津萌津萌 2026-03-02 | 23/1150 |
|
|
|
[考研] 化學(xué) 0703求調(diào)劑 總分293 一志愿211 +3 | 土土小蟲 2026-03-03 | 3/150 |
|
|
|
[考研] 291求調(diào)劑 +4 | Afy123456 2026-03-03 | 7/350 |
|
|
|
[考研] 276求調(diào)劑 +8 | 路lyh123 2026-02-28 | 10/500 |
|
|
|
[考研] 292求調(diào)劑 +3 | sgbl 2026-03-03 | 3/150 |
|
|
|
[考研] 289求調(diào)劑 +7 | BrightLL 2026-03-02 | 9/450 |
|
|
|
[考研] 一志愿華中科技大學(xué)，化學(xué)專業(yè)344分，求調(diào)劑 +6 | 邢xing1 2026-03-02 | 6/300 |
|
|
|
[考研] 化工專碩348，一志愿985求調(diào)劑 +8 | 弗格個(gè) 2026-02-28 | 11/550 |
|
|
|
[考研] 化工270求調(diào)劑 +10 | 什么名字qwq 2026-03-02 | 10/500 |
|
|
|
[考研] 080500材料科學(xué)與工程 +4 | 202114020319 2026-03-03 | 4/200 |
|
|
|
[考博] 26申博 +4 | north， 2026-02-28 | 4/200 |
|
|
|
[考研] 285求調(diào)劑 +9 | 滿頭大汗的學(xué)生 2026-02-28 | 9/450 |
|
|
|
[考研] 材料工程274求調(diào)劑 +5 | Lilithan 2026-03-01 | 5/250 |
|
|
|
[考博] 誠(chéng)招農(nóng)業(yè)博士 +3 | 心欣向榮 2026-02-28 | 3/150 |
|
|
|
[考研] 材料與化工328求調(diào)劑 +3 | 。，。，。，。i 2026-03-02 | 3/150 |
|
|
|
[考研] 264求調(diào)劑 +4 | 巴拉巴拉根556 2026-02-28 | 4/200 |
|
|
|
[考研] 275求調(diào)劑 +3 | L-xin? 2026-03-01 | 6/300 |
|
|
|
[考研] 調(diào)劑 +3 | 簡(jiǎn)木ChuFront 2026-02-28 | 3/150 |
|
|
|
[論文投稿]
求助coordination chemistry reviews 的寫作模板
 10+3
|
ljplijiapeng 2026-02-27 | 4/200 |
|