| 查看: 2331 | 回復(fù): 43 | |||
[求助]
緊急求助:運(yùn)行pwscf出現(xiàn)對(duì)稱性操作警告
|
|
我在運(yùn)行pwscf進(jìn)行結(jié)構(gòu)優(yōu)化時(shí),出現(xiàn)下面的錯(cuò)誤,麻煩大家?guī)兔χ更c(diǎn)一下啊。 Subspace diagonalization in iterative solution of the eigenvalue problem: scalapack distributed-memory algorithm (size of sub-group: 2* 2 procs) warning: symmetry operation # 3 not compatible with FFT grid. -1 0 0 1 1 0 0 0 -1 warning: symmetry operation # 4 not compatible with FFT grid. 1 0 0 -1 -1 0 0 0 -1 warning: symmetry operation # 7 not compatible with FFT grid. 1 0 0 -1 -1 0 0 0 1 warning: symmetry operation # 8 not compatible with FFT grid. -1 0 0 1 1 0 0 0 1 輸入文件: &CONTROL calculation = 'relax', restart_mode = 'from_scratch', prefix = "101-2", pseudo_dir = "/pkg/suse11/quantum/5.0/pseudo", outdir = "/home/n7388683/QE/tmp", tstress = .true tprnfor = .true / &SYSTEM ibrav = 12, a = 3.287, b = 6.237385, c = 23.868112, cosab = -0.263491891704172, nat = 8, ntyp = 2, occupations = 'smearing', degauss = 0.02D0, ecutwfc = 30, ecutrho = 300, / &ELECTRONS conv_thr = 1.D-8, mixing_beta = 0.4, / &IONS ion_dynamics = 'bfgs', / ATOMIC_SPECIES Zn 65.3799972534 Zn.pbe-van.UPF O 15.9989995956 O.pbe-van_ak.UPF ATOMIC_POSITIONS crystal O 0.2199834158725568 0.3844112496986048 0.0175122228436065 O 0.8318951398867058 0.6082346977269033 0.1225855599052453 O 0.5632431858843476 0.0709307897221868 0.1050733370616389 O 0.9513314618701982 0.8471073416938884 0.0000000000000000 Zn 0.0741775052813360 0.0927994285161632 0.0569885856639280 Zn 0.6860892292954852 0.3166228765444618 0.1620619227255668 Zn 0.8055255512789776 0.5554955205114469 0.0394763628203215 Zn 0.4174372752931270 0.7793189685397455 0.1445496998819603 K_POINTS automatic 8 8 1 0 0 0 計(jì)算的是氧化鋅的一個(gè)極性面,會(huì)不會(huì)因?yàn)檫@個(gè)面不穩(wěn)定出現(xiàn)上面情況啊,計(jì)算其他晶面的時(shí)候都沒出現(xiàn)過,謝謝了。! |
木蟲 (正式寫手)
|
QE5.0改成這樣的警告了?挺有意思的 是這樣的,你所使用的實(shí)空間FFT網(wǎng)格的格點(diǎn)并沒有剛好取到平移對(duì)稱性所在的位置,所以它只能忽略你體系的平移對(duì)稱性。 解決方法也好辦: 1.如果你對(duì)對(duì)稱群沒要求,或者不用進(jìn)行聲子模式分析,可忽略 2.結(jié)構(gòu)優(yōu)化時(shí)不用管,真正自洽計(jì)算時(shí)如果出現(xiàn)該警告再考慮消除。 3.如果想消除本警告,可設(shè)置nr1,nr2,nr3的值,這三個(gè)值決定了你FFT網(wǎng)格在三個(gè)方向上有多少格點(diǎn)。如果不清楚要取多少,請(qǐng)參考輸出文件里的“FFT dimensions”一行,然后酌情稍作修改,使得你取的FFT網(wǎng)格中包含平移對(duì)稱性。 |
木蟲 (正式寫手)
木蟲 (正式寫手)
|
可以參考這個(gè)http://www.gaoyang168.com/html/201010/2409072.html 雖然里面沒有提出來(lái)如何解決,但是如果你能像帖子里那樣改變一下晶格類型,將警告 warning: symmetry operation # 3 not compatible with FFT grid. -1 0 0 1 1 0 0 0 -1 變成如下形式: warning: symmetry operation # 3 not allowed. fractional translation: a b c in crystal coordinates 出現(xiàn)平移坐標(biāo)(a,b,c)之后,具體解決方案就是如下: 設(shè)置的nr1一定要能整除a, nr2要能整除b, nr3要能整除c,然后就沒有警告了。 |
木蟲 (正式寫手)
|
你好,我在輸出文件里找到兩個(gè)FFT dimensions,應(yīng)該看哪一個(gè)? Dense grid: 279535 G-vectors FFT dimensions: ( 36, 72, 250) Smooth grid: 70809 G-vectors FFT dimensions: ( 24, 45, 160) 另外你說(shuō)的貼子我之前看過了,怎么改晶格類型,讓警告變成下面形式呢?我已經(jīng)按照那貼子的建議用了a,b,c,cosab的。 warning: symmetry operation # 3 not allowed. fractional translation: a b c in crystal coordinates 還有改ecutrho又該怎么改呢,麻煩您拿我這例子講下啊,謝謝了! |
木蟲 (正式寫手)
木蟲 (正式寫手)
|
結(jié)果出來(lái)了,并沒有出現(xiàn)任何的警告。然后我用優(yōu)化好的結(jié)構(gòu)去做scf計(jì)算也沒有出現(xiàn)警告。你用的QE是什么版本? 另外,nr1, nr2, nr3的值是Dense grid,不是smooth grid。 比方說(shuō)你的Dense grid是36, 72, 250,你可以設(shè)置nr1=40, nr2=80, nr3=300,然后運(yùn)行pw.x,再看輸出文件就會(huì)發(fā)現(xiàn)Dense grid變成了40, 80 ,300。 但我重復(fù)你的計(jì)算并沒有任何這個(gè)警告啊。。。 |
木蟲 (小有名氣)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿南昌大學(xué),327分,材料與化工085600 +6 | Ncdx123456 2026-03-19 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 工科材料085601 279求調(diào)劑 +7 | 困于星晨 2026-03-17 | 9/450 |
|
|
[考研]
|
不想起名字112 2026-03-19 | 3/150 |
|
|
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +5 | 楓橋ZL 2026-03-18 | 7/350 |
|
|
[考研] 286求調(diào)劑 +6 | lemonzzn 2026-03-16 | 10/500 |
|
|
[考研] 085600材料與化工求調(diào)劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|
|
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +7 | 小材化本科 2026-03-18 | 7/350 |
|
|
[考研] 一志愿天大材料與化工(085600)總分338 +5 | 蔡大美女 2026-03-13 | 5/250 |
|
|
[考研] 一志愿985,本科211,0817化學(xué)工程與技術(shù)319求調(diào)劑 +10 | Liwangman 2026-03-15 | 10/500 |
|
|
[考研] 332求調(diào)劑 +3 | ydfyh 2026-03-17 | 3/150 |
|
|
[考研] 311求調(diào)劑 +11 | 冬十三 2026-03-15 | 12/600 |
|
|
[考研] 一志愿西南交大,求調(diào)劑 +4 | 材化逐夢(mèng)人 2026-03-18 | 4/200 |
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
[考研] 材料與化工專碩調(diào)劑 +5 | heming3743 2026-03-16 | 5/250 |
|
|
[考研]
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
[基金申請(qǐng)]
今年的國(guó)基金是打分制嗎?
50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
[考研] 070303 總分349求調(diào)劑 +3 | LJY9966 2026-03-15 | 5/250 |
|
|
[考研] 求老師收留調(diào)劑 +4 | jiang姜66 2026-03-14 | 5/250 |
|
|
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3 | 麥田耕者 2026-03-14 | 3/150 |
|