 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1784 | 回復(fù): 3 | ||
| 【懸賞金幣】回答本帖問題，作者秋日以往將贈(zèng)送您 5 個(gè)金幣 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||
秋日以往新蟲 (初入文壇)
|
[求助]
QE添加范德瓦爾茲力計(jì)算原子受力結(jié)果為NaN， 求教應(yīng)該怎么辦 已有1人參與
|
|
|
目的是弛豫獲得CsPbI3的立方相晶體結(jié)構(gòu) 材料的范德瓦爾茲力比較強(qiáng)，對晶格影響比較大，在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)閉對稱性（ 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 猜測是由于受力沒有成功計(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ù)此樓» 猜你喜歡
085600材料調(diào)劑 總分330
已經(jīng)有6人回復(fù)
-
歡迎211本科同學(xué)，過A區(qū)國家線，A區(qū)非偏遠(yuǎn)一本，交叉學(xué)科課題組
已經(jīng)有30人回復(fù)
-
高分子化學(xué)與物理調(diào)劑
已經(jīng)有17人回復(fù)
-
化學(xué)專業(yè)調(diào)劑
已經(jīng)有6人回復(fù)
-
0703 東華大學(xué) 理學(xué)化學(xué)均過 a 線，總分 281 求調(diào)劑
已經(jīng)有11人回復(fù)
-
材料工程調(diào)劑
已經(jīng)有6人回復(fù)
-
材料學(xué)碩318求調(diào)劑
已經(jīng)有13人回復(fù)
-
08工科求調(diào)劑
已經(jīng)有4人回復(fù)
-
A區(qū)一本交叉課題組，低分調(diào)劑，招收機(jī)械電子信息通信等交叉方向
已經(jīng)有40人回復(fù)
-
070300調(diào)劑一志愿南京大學(xué)化學(xué)專業(yè)，初試273分，有論文有專利有競賽，求調(diào)劑B區(qū)211
已經(jīng)有11人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
lq8989
至尊木蟲 (著名寫手)
- 應(yīng)助: 4 (幼兒園)
- 金幣: 10104.1
- 散金: 169
- 紅花: 2
- 帖子: 1541
- 在線: 225.1小時(shí)
- 蟲號(hào): 1452335
- 注冊: 2011-10-20
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
3樓2022-09-14 14:10:03
lq8989
至尊木蟲 (著名寫手)
- 應(yīng)助: 4 (幼兒園)
- 金幣: 10104.1
- 散金: 169
- 紅花: 2
- 帖子: 1541
- 在線: 225.1小時(shí)
- 蟲號(hào): 1452335
- 注冊: 2011-10-20
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
【答案】應(yīng)助回帖
感謝參與，應(yīng)助指數(shù) +1
|
QE開發(fā)團(tuán)隊(duì)非常反對在高度對稱的體系里使用ibrav=0（其實(shí)最好永遠(yuǎn)不要用ibrav=0），因?yàn)槟承゜ug容易在ibrav=0時(shí)在對稱性上出現(xiàn)錯(cuò)誤。對于這個(gè)簡單立方體系，應(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à)原子，來完全保證對稱性正確。 |
2樓2022-09-14 14:08:49
秋日以往
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 19.8
- 帖子: 18
- 在線: 3.9小時(shí)
- 蟲號(hào): 8541249
- 注冊: 2018-04-16
- 性別: MM
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
4樓2022-11-28 15:32:49
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研]
 高分子化學(xué)與物理調(diào)劑
+6
|
好好好1233 2026-02-28 | 17/850 |
|
|---|---|---|---|---|
|
[考研] 新疆大學(xué)地質(zhì)與礦業(yè)工程學(xué)院招生 +7 | another12 2026-03-04 | 12/600 |
|
|
|
[考研] 不限學(xué)校專業(yè)的調(diào)劑同學(xué)看過來 +5 | 啊擺啊擺 2026-03-05 | 9/450 |
|
|
|
[考博] 26申博-目前4篇SCI一作 +4 | chen_2024 2026-03-02 | 4/200 |
|
|
|
[考研] 一志愿武理085601專碩347分求調(diào)劑 +4 | 啊歐歐歐 2026-03-04 | 5/250 |
|
|
|
[基金申請]
 本子寫完了，給DS兄弟看了，得了92分
+6
|
Doma 2026-03-01 | 10/500 |
|
|
|
[考研] 一志愿武漢理工大學(xué)-085602-總分296分-求調(diào)劑 +7 | 紫川葡柚 2026-03-04 | 7/350 |
|
|
|
[考研] 學(xué)碩材料275調(diào)劑 +9 | 路三三 2026-03-03 | 9/450 |
|
|
|
[考研] 0857求調(diào)劑 +4 | CDDD1013 2026-03-04 | 4/200 |
|
|
|
[考研] 能動(dòng)297求調(diào)劑，本科川大 +4 | 邵11 2026-03-04 | 4/200 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | G212 2026-03-03 | 5/250 |
|
|
|
[考研] 化工專碩348，一志愿985求調(diào)劑 +8 | 弗格個(gè) 2026-02-28 | 11/550 |
|
|
|
[考研] 278求調(diào)劑 +3 | 滿天星11_22 2026-03-02 | 3/150 |
|
|
|
[考研] 化學(xué)0703求調(diào)劑 學(xué)碩 理/工科均可 總分279 +3 | 1一11 2026-03-03 | 5/250 |
|
|
|
[考研] 0856調(diào)劑 +10 | 劉夢微 2026-02-28 | 10/500 |
|
|
|
[考研] 11408，學(xué)碩276求調(diào)劑 +3 | 崔wj 2026-03-02 | 5/250 |
|
|
|
[考研] 0856材料調(diào)劑 +5 | 沿岸有貝殼OUC 2026-03-02 | 5/250 |
|
|
|
[考研] 材料284求調(diào)劑，一志愿鄭州大學(xué)英一數(shù)二專碩 +15 | 想上岸的土撥鼠 2026-02-28 | 15/750 |
|
|
|
[考研] 一志愿東北大學(xué)材料專碩328，求調(diào)劑 +3 | shs1083 2026-03-02 | 3/150 |
|
|
|
[考研]
尋找調(diào)劑
+4
|
LYidhsjabdj 2026-02-28 | 4/200 |
|