| 5 | 1/1 | 返回列表 |
| 查看: 3304 | 回復(fù): 8 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||||
[交流]
【求助】pwscf錯(cuò)誤提示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點(diǎn)省略了些 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 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 那位大俠幫我解決一下,在下感激不盡! |
» 搶金幣啦!回帖就可以得到:
+2/142
+1/129
+1/91
+1/73
+1/57
+1/40
+1/39
+1/16
+1/16
+1/16
+1/11
+1/11
+1/10
+1/9
+1/6
+1/5
+1/5
+3/4
+1/4
+1/3
金蟲 (著名寫手)
|
nbnd=285, 太大了。原因是你的原胞很小,倒空間大,平面波數(shù)目不夠多,而你的nbnd太多,就是說你的導(dǎo)帶數(shù)目太多了。一般你選取電子數(shù)的二分之一加百分之20足夠了。nbnd并不是越大越好的。 還有另外一個(gè)辦法,就是增大ecutwfc和ecutrho,這樣也可以解決整個(gè)問題。你試試,不對(duì)再說 |
專家顧問 (正式寫手)
![]() |
專家經(jīng)驗(yàn): +129 |
|
我對(duì)磁性金屬體系不太熟悉,不過你考慮以下問題: 1. 能帶是不是太多了 2. ecutwfc是不是太小了,相對(duì)于你這么大的體系來說 3. 贗勢(shì)對(duì)么?合適么? 4. 你這個(gè)體系優(yōu)化了么?沒優(yōu)化的話 forc_conv_thr = 1.0D-5 達(dá)得到么? |

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研]
|
Gymno 2026-03-30 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 材料與化工調(diào)劑一志愿大連海事085600,349 +6 | 吃的不少 2026-03-30 | 6/300 |
|
|
[考研] 生物學(xué)學(xué)碩,一志愿湖南大學(xué),初試成績(jī)338 +7 | YYYYYNNNNN 2026-03-26 | 9/450 |
|
|
[考研] 332求調(diào)劑 +14 | 032500 2026-03-25 | 14/700 |
|
|
[考研] 359求調(diào)劑 +5 | 王了個(gè)楠 2026-03-25 | 5/250 |
|
|
[考研] 26考研-291分-廈門大學(xué)(085601)-柔性電子學(xué)院材料工程專業(yè)求調(diào)劑 +5 | min3 2026-03-24 | 6/300 |
|
|
[考研] 322求調(diào)劑:一志愿湖南大學(xué) 材料與化工(085600),已過六級(jí)。 +9 | XX小鄧 2026-03-29 | 9/450 |
|
|
[考研] 化學(xué)0703 調(diào)劑 306分 一志愿211 +8 | 26要上岸 2026-03-28 | 8/400 |
|
|
[考研] 070300一志愿211,312分求調(diào)劑院校 +4 | 小黃鴨寶 2026-03-30 | 4/200 |
|
|
[考研] 材料與化工272求調(diào)劑 +21 | 阿斯蒂芬2004 2026-03-28 | 21/1050 |
|
|
[考研] 085600,專業(yè)課化工原理,321分求調(diào)劑 +5 | 大饞小子 2026-03-28 | 5/250 |
|
|
[考研] 調(diào)劑考研 +3 | 王杰一 2026-03-29 | 3/150 |
|
|
[考研] 11408軟件工程求調(diào)劑 +3 | Qiu學(xué)ing 2026-03-28 | 3/150 |
|
|
[考研] 299求調(diào)劑 +7 | 嗯嗯嗯嗯2 2026-03-27 | 7/350 |
|
|
[考研] 265求調(diào)劑11408 +3 | 劉小鹿lu 2026-03-27 | 3/150 |
|
|
[考博] 26申博 +3 | 加油沖。 2026-03-26 | 3/150 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4 | 崔wj 2026-03-26 | 4/200 |
|
|
[考研] 351求調(diào)劑 +4 | 麥克阿磊 2026-03-24 | 4/200 |
|
|
[考研] 環(huán)境專碩324分求調(diào)劑推薦 +5 | 軒小寧—— 2026-03-26 | 5/250 |
|
|
[考研]
|
13659058978 2026-03-24 | 4/200 |
|