版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

gaods101@163

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 510.3
帖子: 51
在線: 32.3小時
蟲號: 1096198

[交流] 【求助】pwscf錯誤提示diagonalization(ZHEGV*) falied

out輸出：Program PWSCF v.4.2       starts on 23Dec2010 at 10:15:10

   This program is part of the open-source Quantum ESPRESSO suite
   for quantum simulation of materials; please acknowledge
      "P. Giannozzi et al., J. Phys.:Condens. Matter 21 395502 (2009);
      URL http://www.quantum-espresso.org",
   in publications or presentations arising from this work. More details at
   http://www.quantum-espresso.org/ ... ng_Quantum-ESPRESSO

   Parallel version (MPI), running on    1 processors

   Current dimensions of program PWSCF are:
   Max number of different atomic species (ntypx) = 10
   Max number of k-points (npk) =  40000
   Max angular momentum in pseudopotentials (lmaxx) =  3
   Waiting for input...
   file Ni.pbe-nd-rrkjus.UPF: wavefunction(s)  4S renormalized

   Subspace diagonalization in iterative solution of the eigenvalue problem:
   Too few procs for parallel algorithm: we need at least 4 procs per pool
   a serial algorithm will be used

   Planes per process (thick) : nr3 = 320 npp =  320 ncplane = 625
   Planes per process (smooth): nr3s= 216 npps=  216 ncplanes= 324

   Proc/  planes cols    G planes cols G    columns  G
   Pool    (dense grid)    (smooth grid)    (wavefct grid)
      1 320 409 82567  216 187 25245    61    5005

   bravais-lattice index    =          12
   lattice parameter (a_0) =    4.7089  a.u.
   unit-cell volume       = 1134.9842 (a.u.)^3
   number of atoms/cell    =          9
   number of atomic types =          1
   number of electrons    =       90.00
   number of Kohn-Sham states=       285
   kinetic-energy cutoff    =    30.0000  Ry
   charge density cutoff    =    265.0000  Ry
   convergence threshold    =    1.0E-06
   mixing beta             =    0.2000
   number of iterations used =          8  local-TF  mixing
   Exchange-correlation    =  SLA  PW PBE  PBE (1434)
   EXX-fraction             =       0.00

   celldm(1)= 4.708905  celldm(2)= 1.000000  celldm(3)=  12.551612
   celldm(4)=  -0.500000  celldm(5)= 0.000000  celldm(6)= 0.000000

   crystal axes: (cart. coord. in units of a_0)
            a(1) = (  1.000000  0.000000  0.000000 )
            a(2) = ( -0.500000  0.866025  0.000000 )
            a(3) = (  0.000000  0.000000 12.551612 )

   reciprocal axes: (cart. coord. in units 2 pi/a_0)
            b(1) = (  1.000000  0.577350  0.000000 )
            b(2) = (  0.000000  1.154701  0.000000 )
            b(3) = (  0.000000  0.000000  0.079671 )

   PseudoPot. # 1 for Ni read from file Ni.pbe-nd-rrkjus.UPF
   Pseudo is Ultrasoft + core correction, Zval = 10.0
   Generated by new atomic code, or converted to UPF format
   Using radial grid of 1203 points,  6 beta functions with:
            l(1) = 0
            l(2) = 0
            l(3) = 1
            l(4) = 1
            l(5) = 2
            l(6) = 2
   Q(r) pseudized with 0 coefficients

   atomic species valence mass    pseudopotential
      Ni          10.00    1.00000    Ni( 1.00)

   Starting magnetic structure
   atomic species magnetization
      Ni          0.200

   No symmetry found

Cartesian axes

   site n.    atom                positions (a_0 units)
      1          Ni  tau(  1) = ( 0.3000000 0.4618802 0.8164967  )
      2          Ni  tau(  2) = ( 0.3000000 0.4618802 3.2659867  )
      3          Ni  tau(  3) = ( 0.3000000 0.4618802 5.7154768  )
      4          Ni  tau(  4) = (  -0.2000000 0.7505554 0.0000000  )
      5          Ni  tau(  5) = ( 0.8000000 0.1732051 1.6329934  )
      6          Ni  tau(  6) = (  -0.7248907 1.6596927 2.4494901  )
      7          Ni  tau(  7) = ( 0.2751093 1.0823424 4.0824834  )
      8          Ni  tau(  8) = (  -0.7248907 1.6596927 4.8989801  )
      9          Ni  tau(  9) = ( 0.2751093 1.0823424 6.5319735  )

   number of k points= 82  gaussian broad. (Ry)=  0.0100    ngauss = 0
                     cart. coord. in units 2pi/a_0
      k( 1) = ( 0.0000000 0.0000000 0.0000000), wk = 0.0123457
      k( 2) = ( 0.0000000 0.1283001 0.0000000), wk = 0.0246914
      k( 3) = ( 0.0000000 0.2566001 0.0000000), wk = 0.0246914
      k( 4) = ( 0.0000000 0.3849002 0.0000000), wk = 0.0246914
      k( 5) = ( 0.0000000 0.5132002 0.0000000), wk = 0.0246914
      k( 6) = ( 0.1111111 0.0641500 0.0000000), wk = 0.0246914
      k( 7) = ( 0.1111111 0.1924501 0.0000000), wk = 0.0246914
.........
   k點省略了些

   G cutoff =  148.8421  (  82567 G-vectors)    FFT grid: ( 25, 25,320)
   G cutoff = 67.4002  (  25245 G-vectors)  smooth grid: ( 18, 18,216)

   Largest allocated arrays    est. size (Mb)    dimensions
      Kohn-Sham Wavefunctions       13.86 Mb    ( 3186, 285)
      NL pseudopotentials          7.88 Mb    ( 3186, 162)
      Each V/rho on FFT grid       6.10 Mb    ( 200000, 2)
      Each G-vector array          0.63 Mb    (  82567)
      G-vector shells                0.03 Mb    ( 4350)
   Largest temporary arrays    est. size (Mb)    dimensions
      Auxiliary wavefunctions       55.42 Mb    ( 3186,1140)
      Each subspace H/S matrix    19.83 Mb    ( 1140,1140)
      Each matrix    0.70 Mb    ( 162, 285)
      Arrays for rho mixing       24.41 Mb    ( 200000, 8)
   writing wfc files to a dedicated directory

   Check: negative/imaginary core charge= -0.000022 0.000000

   Initial potential from superposition of free atoms
   Check: negative starting charge=(component1): -0.000152
   Check: negative starting charge=(component2): -0.000102

   starting charge 89.99587, renormalised to 90.00000

   negative rho (up, down):  0.152E-03 0.102E-03
   Starting wfc are 54 atomic +  231 random wfc

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
   from cdiaghg : error #    298
   diagonalization (ZHEGV*) failed
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

   stopping ...
以下是我的輸出文件：
&CONTROL
               calculation = 'scf' ,
            restart_mode = 'from_scratch' ,
                  outdir = './' ,
                  wfcdir = './' ,
               pseudo_dir = '/home/abc/pwscf/espresso-4.2/pseudo/' ,
                  prefix = 'Ni' ,
            etot_conv_thr = 1.0D-6 ,
            forc_conv_thr = 1.0D-5 ,
                  tstress = .true. ,
/
&SYSTEM
                     ibrav = 12,
               celldm(1) = 4.70890458,
               celldm(2) = 1,
               celldm(3) = 12.5516119789,
               celldm(4) = -0.5,
                     nat = 9,
                     ntyp = 1,
                  ecutwfc = 30 ,
                  ecutrho = 265 ,
                     nbnd = 285,
               occupations = 'smearing' ,
                  degauss = 0.01 ,

                  smearing = 'gaussian' ,
                     nspin=2,
starting_magnetization(1)=0.2,
/
&ELECTRONS
         electron_maxstep = 500,
               mixing_mode = 'local-TF' ,
               mixing_beta = 0.2 ,
/
ATOMIC_SPECIES
Ni 1.00000  Ni.pbe-nd-rrkjus.UPF
ATOMIC_POSITIONS crystal
Ni 0.566666678 0.533333369 0.065051141
Ni 0.566666678 0.533333369 0.260204565
Ni 0.566666678 0.533333369 0.455357989
Ni 0.233333344 0.866666702 0.000000000
Ni 0.900000011 0.200000035 0.130102282
Ni 0.233333344 1.916448002 0.195153424
Ni 0.900000011 1.249781335 0.325255706
Ni 0.233333344 1.916448002 0.390306847
Ni 0.900000011 1.249781335 0.52040913
K_POINTS automatic
  9 9 1 0 0 0
那位大俠幫我解決一下，在下感激不盡！

回復(fù)此樓

» 猜你喜歡

309求調(diào)劑已經(jīng)有8人回復(fù)
一志愿南航 335分 | 0856 | GPA 4.07 | 有科研經(jīng)歷已經(jīng)有8人回復(fù)
求調(diào)劑已經(jīng)有5人回復(fù)
求化學(xué)調(diào)劑已經(jīng)有5人回復(fù)
085600，材料與化工321分求調(diào)劑已經(jīng)有10人回復(fù)
0703 化學(xué) 求調(diào)劑，一志愿山東大學(xué) 342 分已經(jīng)有5人回復(fù)
一志愿南開大學(xué)0710生物學(xué)359求調(diào)劑已經(jīng)有3人回復(fù)
085600，專業(yè)課化工原理，320分求調(diào)劑已經(jīng)有4人回復(fù)
材料與化工272求調(diào)劑已經(jīng)有16人回復(fù)
290求調(diào)劑已經(jīng)有3人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

PWSCF在安裝完之后用pwgui運行，馬上就stop了，然后出現(xiàn)了錯誤已經(jīng)有5人回復(fù)
【已解決】【求助】為何PWscf不能得到類似于CASTEP那樣的收斂測試曲線呢？已經(jīng)有14人回復(fù)
【求助】pwscf并行已經(jīng)有7人回復(fù)
【求助】利用abinit或pwscf軟件在單機上可以計算半導(dǎo)體材料的哪些性質(zhì)? 已經(jīng)有8人回復(fù)
【原創(chuàng)】PWSCF中的nbnd和degauss 已經(jīng)有11人回復(fù)
【求助】關(guān)于PWSCF介電常數(shù)和聲子譜的計算已經(jīng)有7人回復(fù)
【求助】pwscf中的問題已經(jīng)有4人回復(fù)
【求助】PWSCF計算介電極化時為何出現(xiàn)如下提示？已經(jīng)有19人回復(fù)

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

1樓 2010-12-23 11:29:22

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

lfhuang

木蟲 (著名寫手)

樓主能否說明一下您輸入文件里面：
                  ecutwfc = 30 ,
                  ecutrho = 265 ,
                     nbnd = 285,
三個參數(shù)為什么這么設(shè)置？還有您的編譯器，數(shù)學(xué)庫用的是什么？

贊一下

回復(fù)此樓

2樓2010-12-23 13:15:13

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

gaods101@163

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 510.3
帖子: 51
在線: 32.3小時
蟲號: 1096198

引用回帖:

Originally posted by lfhuang at 2010-12-23 13:15:13:
樓主能否說明一下您輸入文件里面：
                  ecutwfc = 30 ,
                  ecutrho = 265 ,
                     nbnd = 285,
三個參數(shù)為什么這么設(shè)置？還有您的編譯器，數(shù)學(xué)庫用的是 ...

你好:ecutwfc和ecutrho是我經(jīng)過測試后決定的，一般來說ecutrho是ecutwfc的4倍，但我設(shè)定之后收斂不是很好，增加ecutwfc導(dǎo)致運算速度很慢，所以我增加了ecutrho,nbnd是隨便設(shè)定的，多少我感覺無所謂，多了能帶多一些好看些，少了能帶稀疏一些。我用的是fortran編輯器，莫非我的錯誤跟這些有關(guān)？\
感謝回帖！�。�

贊一下

回復(fù)此樓

3樓2010-12-23 16:10:20

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

oxox6085

專家顧問 (正式寫手)

★ ★
zhang668(金幣+2):多謝交流 2010-12-23 19:02:24
gaods101@163(金幣+5):謝謝你的意見，我剛學(xué)不是太久，很多參數(shù)還在摸索階段，希望多多指導(dǎo)！ 2010-12-24 10:33:12
gaods101@163(金幣+5): 2011-03-04 20:20:23

我對磁性金屬體系不太熟悉，不過你考慮以下問題：
1. 能帶是不是太多了
2. ecutwfc是不是太小了，相對于你這么大的體系來說
3. 贗勢對么？合適么？
4. 你這個體系優(yōu)化了么？沒優(yōu)化的話 forc_conv_thr = 1.0D-5 達得到么？

贊一下(1人)

回復(fù)此樓

4樓2010-12-23 17:51:12

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

goldenfisher

金蟲 (著名寫手)

★ ★
zhang668(金幣+2):多謝交流 2010-12-23 19:02:07
gaods101@163(金幣+5):謝謝，第一次發(fā)帖，沒想到這么多人幫忙，太感謝了，我剛學(xué)沒多久，很多參數(shù)的設(shè)定還在摸索之中！多多照顧！ 2010-12-24 10:31:15

nbnd=285,
太大了。原因是你的原胞很小，倒空間大，平面波數(shù)目不夠多，而你的nbnd太多，就是說你的導(dǎo)帶數(shù)目太多了。一般你選取電子數(shù)的二分之一加百分之20足夠了。nbnd并不是越大越好的。
還有另外一個辦法，就是增大ecutwfc和ecutrho，這樣也可以解決整個問題。你試試，不對再說

贊一下(1人)

回復(fù)此樓

5樓2010-12-23 17:55:58

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

lfhuang

木蟲 (著名寫手)

★ ★
zhang668(金幣+2):多謝交流 2010-12-23 19:01:56
gaods101@163(金幣+5):感謝你多次的為我提意見，我學(xué)這個不是很久，很多參數(shù)還在探索之中，希望多多給予指導(dǎo) 2010-12-24 10:34:48
gaods101@163(金幣+5): 2011-03-04 20:21:05

引用回帖:

Originally posted by gaods101@163 at 2010-12-23 16:10:20:
一般來說ecutrho是ecutwfc的4倍

一般來說ecutrho是ecutwfc的10倍(to ultrasoft pseudopotential)

引用回帖:

Originally posted by gaods101@163 at 2010-12-23 16:10:20:
一般來說ecutrho是ecutwfc的4倍nbnd是隨便設(shè)定的，多少我感覺無所謂，多了能帶多 ...

您的問題很可能是nbnd太大了。

PS: ecutwfc=30還是可以的，很常用的一個數(shù)值。不過具體如何還得測試。

贊一下(1人)

回復(fù)此樓

6樓2010-12-23 18:44:20

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

gaods101@163

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 510.3
帖子: 51
在線: 32.3小時
蟲號: 1096198

引用回帖:

Originally posted by oxox6085 at 2010-12-23 17:51:12:
我對磁性金屬體系不太熟悉，不過你考慮以下問題：
1. 能帶是不是太多了
2. ecutwfc是不是太小了，相對于你這么大的體系來說
3. 贗勢對么？合適么？
4. 你這個體系優(yōu)化了么？沒優(yōu)化的話 forc_conv_thr = 1.0D- ...

ecuwfc感覺還是可以的，這個體系并不是很大，我試過ecuwfc增大，但速度不是很樂觀。我是沒做過優(yōu)化的，我一直在想一個問題，是不是在做自洽之前必須進行一次優(yōu)化？最后我感覺應(yīng)該是能帶錯了，但還沒試，而我是做了另一個工作，做了一個超晶胞，原子個數(shù)增大了四倍，結(jié)果計算很合理！可能是能帶問題吧，我再試試，完了再看結(jié)果吧

贊一下

回復(fù)此樓

8樓2010-12-24 10:45:42

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

gaods101@163

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 510.3
帖子: 51
在線: 32.3小時
蟲號: 1096198

引用回帖:

Originally posted by goldenfisher at 2010-12-23 17:55:58:
nbnd=285,
太大了。原因是你的原胞很小，倒空間大，平面波數(shù)目不夠多，而你的nbnd太多，就是說你的導(dǎo)帶數(shù)目太多了。一般你選取電子數(shù)的二分之一加百分之20足夠了。nbnd并不是越大越好的。
還有另外一個辦法，就 ...

我先試試能帶減少吧，也許你說的對，因為我把原子個數(shù)增加到四倍的時候出結(jié)果了，呵呵，謝謝！

贊一下

回復(fù)此樓

9樓2010-12-24 10:47:08

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

簡單回復(fù)

dashu-973327樓

2010-12-23 21:45 回復(fù)

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 gaods101@163 的主題更新

返回列表

普通表情龍兔虎貓高級回復(fù) (可上傳附件)

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 一志愿南開大學(xué)0710生物學(xué)359求調(diào)劑 +3	兔兔兔111223314 2026-03-29	3/150	2026-03-29 23:14 by Evan_Liu
[碩博家園] 求調(diào)劑有機化學(xué)考研356分 +4	Nadiums 2026-03-25	5/250	2026-03-29 21:00 by 唐沐兒
[考研] 299求調(diào)劑 +10	15188958825 2026-03-25	10/500	2026-03-29 17:51 by 王亮_大連醫(yī)科大
[考研] 一志愿211，335分，0856，求調(diào)劑院校和導(dǎo)師 +5	傾____蕭 2026-03-27	6/300	2026-03-29 16:35 by 唐沐兒
[考研] 材料學(xué)碩333求調(diào)劑 +11	北道巷 2026-03-24	11/550	2026-03-29 08:31 by fmesaito
[考研] 317求調(diào)劑 +6	十閑wx 2026-03-24	6/300	2026-03-28 13:27 by Iveryant
[考研] 070300求調(diào)劑306分 +4	26要上岸 2026-03-27	4/200	2026-03-28 13:06 by 唐沐兒
[考研] 藥學(xué)105500求調(diào)劑 +3	Ssun。。 2026-03-28	3/150	2026-03-28 11:24 by lxf170613
[考研] 085602 307分求調(diào)劑 +7	不知道叫什么！ 2026-03-26	7/350	2026-03-28 09:57 by 神馬都不懂
[考研] 292求調(diào)劑 +14	鵝鵝鵝額額額額?/a> 2026-03-25	15/750	2026-03-28 08:45 by WYUMater
[考研] 求調(diào)劑 +4	零八# 2026-03-27	4/200	2026-03-27 18:07 by yu221
[考研] 材料求調(diào)劑 +5	.m.. 2026-03-25	5/250	2026-03-27 11:08 by 不吃魚的貓
[考研] 276求調(diào)劑。有半年電池和半年高分子實習(xí)經(jīng)歷 +10	材料學(xué)257求調(diào)劑 2026-03-23	11/550	2026-03-27 10:13 by YCIT- LHL
[考研] 求調(diào)劑一志愿本科北科大化學(xué) 343 +6	13831862839 2026-03-24	7/350	2026-03-26 22:57 by 不吃魚的貓
[考研] 321求調(diào)劑 +6	Ymlll 2026-03-24	6/300	2026-03-26 20:50 by 不吃魚的貓
[考研] 0703化學(xué)求調(diào)劑 +3	丹青奶蓋 2026-03-26	5/250	2026-03-26 20:11 by macy2011
[考研] 生物學(xué)學(xué)碩求調(diào)劑 +7	小羊睡著了? 2026-03-23	10/500	2026-03-25 02:24 by 清風(fēng)拂揚。 m
[考研] 化工專碩求調(diào)劑 +3	question挽風(fēng) 2026-03-24	3/150	2026-03-24 18:48 by jhhcooi
[考研] 080500求調(diào)劑 +3	zzzzfan 2026-03-24	3/150	2026-03-24 16:38 by barlinike
[考研] 344求調(diào)劑 +3	desto 2026-03-24	3/150	2026-03-24 10:09 by 搏擊518

24小時熱門版塊排行榜

[交流] 【求助】pwscf錯誤提示diagonalization(ZHEGV*) falied

» 猜你喜歡

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

» 搶金幣啦！回帖就可以得到: 查看全部散金貼

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

» 搶金幣啦！回帖就可以得到:

查看全部散金貼