版塊導航: 正在加載中...

登錄注冊

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

24小時熱門版塊排行榜

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

返回列表

xiaoqiu007

新蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 356.5
散金: 7
帖子: 100
在線: 20.8小時
蟲號: 707167
注冊: 2009-02-23
專業(yè): 金屬功能材料

[求助] Problem - nscf - c_bands: eigenvalues not converged

采用pwscf計算費米面時，前面若干步順利通過，但在band_FS.x步驟時報錯

./bands_FS.x < Bands_NSP.out >bands_fs.out
forrtl: severe (59): list-directed I/O syntax error, unit 5, file stdin
Image             PC       Routine          Line       Source
bands_FS.x       0809DB77  Unknown             Unknown  Unknown
bands_FS.x       0809D197  Unknown             Unknown  Unknown
bands_FS.x       08071A96  Unknown             Unknown  Unknown
bands_FS.x       08050B4C  Unknown             Unknown  Unknown
bands_FS.x       080507E7  Unknown             Unknown  Unknown
bands_FS.x       08061C93  Unknown             Unknown  Unknown

請高手賜教�！菊娌幌Ｍ搯栴}無解】

附1
類似的問題還出現(xiàn)在：
http://qe-forge.org/pipermail/pw_forum/2011-October/096758.html；
http://www.democritos.it/pipermail/pw_forum/2010-July/017605.html；
以及小木蟲上
http://www.gaoyang168.com/html/201207/4667669.html

附件2：band_FS.x程序需要的input_FS如下：
3  7
9.0728
Sc
16 16  8
  1.000000  0.577350  0.000000
  0.000000  1.154701  0.000000
  0.000000  0.000000  0.628094

附3：非自洽計算所用的輸入和輸出結(jié)果文件如下：

3. 設(shè)定nbnd=16
輸入文件
&control
calculation='nscf',
prefix='Sc'
pseudo_dir = '/home/users/spclibtt/home/yexq/shanghai-opt/pw-work/pseudo/',
outdir='/home/users/spclibtt/home/yexq/shanghai-opt/pw-work/environment/tmp/'
/
&system
ibrav = 0,
nat=  2,
ntyp= 1,
ecutwfc = 90,
ecutrho = 400.00,
nbnd=16
/
&electrons
diagonalization='cg'
mixing_beta = 0.7
conv_thr =  1.0d-10
/
#&CELL
#cell_dynamics= 'damp-pr'
# press=0
#/

ATOMIC_SPECIES
Sc 44.9559 Sc.pbe-nsp-van.UPF

CELL_PARAMETERS {bohr}
6.254820382  0.00000000  0.00000000
-3.127410191  5.41683334  0.00000000
0.000000000  0.00000000  9.958412336

ATOMIC_POSITIONS {crystal}
Sc 0.3333333333333286 0.6666666666666714 0.2500000000000000
Sc 0.6666666666666714 0.3333333333333286 0.7500000000000000
K_POINTS
  2601
0.000000 0.000000 0.000000  1.00
0.000000 0.000000 0.078512  1.00
0.000000 0.000000 0.157024  1.00
0.000000 0.000000 0.235535  1.00
0.000000 0.000000 0.314047  1.00
0.000000 0.000000 0.392559  1.00
.....
輸出文件

   Program PWSCF v.5.0.2 (svn rev. 9392) starts on 13Mar2013 at 11:17:57

   This program is part of the open-source Quantum ESPRESSO suite
   for quantum simulation of materials; please cite
      "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/quote.php

   Parallel version (MPI), running on 16 processors
   R & G space division:  proc/nbgrp/npool/nimage =    16

   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...
   Reading input from standard input

   Atomic positions and unit cell read from directory:
   /home/users/spclibtt/home/yexq/shanghai-opt/pw-work/environment/tmp/Sc.save/

   Subspace diagonalization in iterative solution of the eigenvalue problem:
   a serial algorithm will be used

   Parallelization info
   --------------------
   sticks: dense  smooth    PW    G-vecs: dense smooth    PW
   Min       66    59    22                2846    2423    554
   Max       67    60    23                2855    2438    561
   Sum       1069    955 361             45639 38911 8929

   bravais-lattice index    =          0
   lattice parameter (alat)  =    6.2548  a.u.
   unit-cell volume       =    337.4042 (a.u.)^3
   number of atoms/cell    =          2
   number of atomic types =          1
   number of electrons    =       22.00
   number of Kohn-Sham states=          16
   kinetic-energy cutoff    =    90.0000  Ry
   charge density cutoff    =    400.0000  Ry  #應該用到900嗎？
   Exchange-correlation    =  SLA  PW PBE  PBE ( 1 4 3 4 0)
   EXX-fraction             =       0.00

   celldm(1)= 6.254820  celldm(2)= 0.000000  celldm(3)= 0.000000
   celldm(4)= 0.000000  celldm(5)= 0.000000  celldm(6)= 0.000000

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

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

   PseudoPot. # 1 for Sc read from file:
   /home/users/spclibtt/home/yexq/shanghai-opt/pw-work/pseudo/Sc.pbe-nsp-van.UPF
   MD5 check sum: ff19847cdcb3487ca1f25e4b0e9f43ec
   Pseudo is Ultrasoft + core correction, Zval = 11.0
   Generated by new atomic code, or converted to UPF format
   Using radial grid of  849 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  6 coefficients,  rinner = 1.200 1.200 1.200
                                                   1.200 1.200

   atomic species valence mass    pseudopotential
      Sc          11.00 44.95590    Sc( 1.00)

   24 Sym. Ops., with inversion, found (12 have fractional translation)

Cartesian axes

   site n.    atom                positions (alat units)
      1          Sc  tau( 1) = ( 0.0000000 0.5773503 0.3980295  )
      2          Sc  tau( 2) = ( 0.5000000 0.2886751 1.1940885  )

   number of k points=  2601

   Number of k-points >= 100: set verbosity='high' to print them.

   Dense  grid: 45639 G-vectors    FFT dimensions: (  40,  40,  64)

   Smooth grid: 38911 G-vectors    FFT dimensions: (  40,  40,  64)

   Largest allocated arrays    est. size (Mb)    dimensions
      Kohn-Sham Wavefunctions       0.08 Mb    ( 318, 16)
      NL pseudopotentials          0.17 Mb    ( 318, 36)
      Each V/rho on FFT grid       0.10 Mb    ( 6400)
      Each G-vector array          0.02 Mb    ( 2855)
      G-vector shells                0.01 Mb    ( 1357)
   Largest temporary arrays    est. size (Mb)    dimensions
      Each subspace H/S matrix       0.00 Mb    (  16,  16)
      Each matrix    0.01 Mb    (    36, 16)

   The potential is recalculated from file :
   /home/users/spclibtt/home/yexq/shanghai-opt/pw-work/environment/tmp/Sc.save/charge-density.dat

   Starting wfc are 20 randomized atomic wfcs

   total cpu time spent up to now is       1.3 secs

   per-process dynamical memory: 13.6 Mb

   Band Structure Calculation #【啟動了Band Structure Calculation ，但特征值不收斂】
   CG style diagonalization
   c_bands:  1 eigenvalues not converged  #這是什么錯誤
   c_bands:  1 eigenvalues not converged
   c_bands:  2 eigenvalues not converged
   c_bands:  3 eigenvalues not converged
  .........#【略去n個...eigenvalues not converged】
   c_bands:  2 eigenvalues not converged
   c_bands:  3 eigenvalues not converged
   c_bands:  2 eigenvalues not converged
   c_bands:  2 eigenvalues not converged
   c_bands:  1 eigenvalues not converged
   c_bands:  1 eigenvalues not converged

   ethr =  4.55E-13,  avg # of iterations = 57.3

   total cpu time spent up to now is    3955.9 secs

   End of band structure calculation

   Number of k-points >= 100: set verbosity='high' to print the bands.

   highest occupied, lowest unoccupied level (ev):    9.9643 8.6137

   Writing output data file Sc.save

   init_run    :    0.50s CPU    0.67s WALL (    1 calls)
   electrons : 3824.17s CPU 3954.63s WALL (    1 calls)

   Called by init_run:
   wfcinit    :    0.00s CPU    0.00s WALL (    1 calls)
   potinit    :    0.11s CPU    0.11s WALL (    1 calls)

   Called by electrons:
   c_bands    : 3824.16s CPU 3954.62s WALL (    1 calls)
   v_of_rho    :    0.08s CPU    0.09s WALL (    1 calls)
   newd       :    0.04s CPU    0.04s WALL (    1 calls)

   Called by c_bands:
   init_us_2 :    1.08s CPU    1.03s WALL ( 2601 calls)
   ccgdiagg    : 3501.84s CPU 3601.94s WALL ( 12559 calls)
   wfcrot    : 317.80s CPU 319.55s WALL ( 12559 calls)

   Called by *cgdiagg:
   h_psi       : 2831.31s CPU 2893.42s WALL ( 2168039 calls)
   s_psi       : 201.42s CPU 204.42s WALL ( 4323519 calls)
   cdiaghg    :    4.95s CPU    5.16s WALL ( 12559 calls)

   Called by h_psi:
   add_vuspsi : 115.68s CPU 118.56s WALL ( 2168039 calls)

   General routines
   calbec    : 507.36s CPU 523.06s WALL ( 4323519 calls)
   fft       :    0.08s CPU    0.09s WALL (    13 calls)
   ffts       :    0.00s CPU    0.00s WALL (    1 calls)
   fftw       : 2282.87s CPU 2296.81s WALL ( 4733656 calls)
   interpolate  :    0.00s CPU    0.00s WALL (    1 calls)
   davcio    :    0.14s CPU    10.04s WALL ( 2601 calls)

   Parallel routines
   fft_scatter  : 1293.12s CPU 1310.33s WALL ( 4733670 calls)

   PWSCF       :    1h 3m CPU       1h 6m WALL

This run was terminated on:  12:24:32  13Mar2013

=------------------------------------------------------------------------------=
JOB DONE.
=------------------------------------------------------------------------------=

回復此樓

1樓 2013-03-13 14:18:27

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

回帖支持 ( 顯示支持度最高的前 50 名 )

gemucai

木蟲 (正式寫手)

應助: 16 (小學生)
金幣: 1774.6
散金: 157
紅花: 4
帖子: 693
在線: 201.7小時
蟲號: 816478
注冊: 2009-07-27
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

★
WDD880227: 金幣+1, 謝謝指教 2013-03-27 10:21:43

引用回帖:

4樓: Originally posted by xiaoqiu007 at 2013-03-13 22:09:02
做的是非自洽計算，沒找到estimated scf accuracy語句；
奇怪的是如果設(shè)置nbnd=8,則不會出現(xiàn)
c_bands: 1 eigenvalues not converged
的情況。

另外，請問計算費米面之前，是不是要用PWscf對結(jié)構(gòu)進行很 ...

可能是激發(fā)態(tài)確實不太好算吧。
vasp所用的贗勢跟pwscf不一樣，平衡位置肯定也不一樣，vasp優(yōu)化好了，在pwscf里面一般是有受力的。所以，沒辦法，哈哈

贊一下(1人)

回復此樓

5樓2013-03-14 10:53:58

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

gemucai

木蟲 (正式寫手)

應助: 16 (小學生)
金幣: 1774.6
散金: 157
紅花: 4
帖子: 693
在線: 201.7小時
蟲號: 816478
注冊: 2009-07-27
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

★ ★
WDD880227: 金幣+2, 多謝指教 2013-03-27 10:21:56

引用回帖:

6樓: Originally posted by 劉仕晨 at 2013-03-14 10:57:18
不好意思，插到樓主的帖子里問個問題，您好，昨天我的pwscf安裝帖子您也幫助了很多，然后我還有幾個問題，現(xiàn)在我在PC上跑的時候都是按照example在學習使用，但是里面的膺勢都是人家選好的，我想知道我自己算的時候 ...

如何測贗勢，可以參見castep的帖子，過程是一樣的：
http://www.gaoyang168.com/bbs/viewthread.php?tid=5153783

LDA還是GGA，參見此帖：
http://www.gaoyang168.com/bbs/viewthread.php?tid=5464861&fpage=1

贊一下(1人)

回復此樓

7樓2013-03-14 14:40:14

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

goldenfisher

金蟲 (著名寫手)

1ST強帖: 3
應助: 55 (初中生)
貴賓: 0.153
金幣: 5779.1
散金: 2908
紅花: 21
帖子: 1552
在線: 351.5小時
蟲號: 117764
注冊: 2005-11-26
性別: GG
專業(yè): 原子和分子物理

★ ★
WDD880227: 金幣+2, 多謝指教 2013-03-27 10:22:06

OK.
既然你是金屬，在計算nscf的時候就必須設(shè)置nbnd,而且必須設(shè)置為電子占據(jù)能帶數(shù)目的1.2倍以上（程序的說明文件里面有詳細說明）。
為什么在scf的時候不需要呢？因為在scf計算的時候，你設(shè)置了smearing，程序自動為你多設(shè)置了部分能帶。而如果不設(shè)置nbnd又沒有設(shè)置smearing,程序默認的是電子占據(jù)的滿能帶數(shù)，也就是電子數(shù)目的一半，這對金屬而言是不對的。至于為什么不對，我以前的帖子講過很多了，固體物理里面也應該學過。

贊一下(1人)

回復此樓

12樓2013-03-15 21:12:32

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

gemucai

木蟲 (正式寫手)

應助: 16 (小學生)
金幣: 1774.6
散金: 157
紅花: 4
帖子: 693
在線: 201.7小時
蟲號: 816478
注冊: 2009-07-27
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

11樓: Originally posted by 劉仕晨 at 2013-03-15 14:56:10
您好，謝謝，也即是說測試贗勢主要是看跑出來的結(jié)果和文獻進行比對，相差不大即可；我昨天也找了下資料，總結(jié)了下，也就是說一般LDA包括PAW ,PZ ;GGA包括PBE、PW91、B88-P86、BLYP，然后不能混用，所有的元素要么都 ...

測試贗勢主要指測試其收斂性：ecutrho、degauss、k mesh。后面有一點小錯誤，在這里糾正一下：
PAW是與模守恒和超軟相提并論的，它既可能是LDA也可能是GGA。

贊一下(1人)

回復此樓

16樓2013-03-16 15:31:10

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

普通回帖

gemucai

木蟲 (正式寫手)

應助: 16 (小學生)
金幣: 1774.6
散金: 157
紅花: 4
帖子: 693
在線: 201.7小時
蟲號: 816478
注冊: 2009-07-27
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

這還真不是個錯誤，是因為你計算的結(jié)構(gòu)偏離平衡位置太遠了，導致有的本征值不收斂。你查查，出現(xiàn)警告的那句話后面的estimated scf accuracy很大吧。一開始有能量不收斂沒關(guān)系，只要最后一步是收斂的就行了，放心吧。

贊一下

回復此樓

2樓2013-03-13 15:23:24

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

gemucai

木蟲 (正式寫手)

應助: 16 (小學生)
金幣: 1774.6
散金: 157
紅花: 4
帖子: 693
在線: 201.7小時
蟲號: 816478
注冊: 2009-07-27
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

參見http://www.democritos.it/piperma ... ovember/007832.html
當然，費米面的問題我在那個貼子里也是無能為力的。。。

贊一下

回復此樓

3樓2013-03-13 15:32:32

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

xiaoqiu007

新蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 356.5
散金: 7
帖子: 100
在線: 20.8小時
蟲號: 707167
注冊: 2009-02-23
專業(yè): 金屬功能材料

做的是非自洽計算，沒找到estimated scf accuracy語句；
奇怪的是如果設(shè)置nbnd=8,則不會出現(xiàn)
c_bands:  1 eigenvalues not converged
的情況。

   另外，請問計算費米面之前，是不是要用PWscf對結(jié)構(gòu)進行很好的優(yōu)化，然后再做自洽，非自洽運算？
   可不可以VASP優(yōu)化，(利用VASP優(yōu)化得到的位置形狀)PWscf自洽及其后續(xù)計算?

贊一下

回復此樓

4樓2013-03-13 22:09:02

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

劉仕晨

木蟲 (正式寫手)

應助: 30 (小學生)
金幣: 3276.8
散金: 6
紅花: 11
帖子: 319
在線: 179.2小時
蟲號: 2201333
注冊: 2012-12-23
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

3樓: Originally posted by gemucai at 2013-03-13 15:32:32
參見http://www.democritos.it/pipermail/pw_forum/2007-November/007832.html
當然，費米面的問題我在那個貼子里也是無能為力的。。。

不好意思，插到樓主的帖子里問個問題，您好，昨天我的pwscf安裝帖子您也幫助了很多，然后我還有幾個問題，現(xiàn)在我在PC上跑的時候都是按照example在學習使用，但是里面的膺勢都是人家選好的，我想知道我自己算的時候我該怎么選呢，我在網(wǎng)上找了下，人家都說要測試，要怎么測試呢，是不是從官網(wǎng)下了，然后一個一個測試，算自洽測試么？還有，我怎么看這個膺勢是LDA還是GGA呢，直接打開膺勢文件看么？不好意思，我們組用的是VASP，現(xiàn)在導師突然讓我學習pwscf，我是初學者，很多問題不懂，謝謝了

贊一下

回復此樓

6樓2013-03-14 10:57:18

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

xiaoqiu007

新蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 356.5
散金: 7
帖子: 100
在線: 20.8小時
蟲號: 707167
注冊: 2009-02-23
專業(yè): 金屬功能材料

問題依舊沒解決，頭痛：
1) 優(yōu)化了結(jié)構(gòu)；并在非自洽計算中先后考察了
2）.不同的贗勢：USPP與模守恒贗勢，
3) .修改了mixing_beta--0.2，
4) .增加mixing_ndim=10，
5) . 換了對角化方法'cg'--'david';
6) . 修改了ecutwfc及ecutrho
.......
   問題依然，真的沒轍了。難道真的無解？  PWscf 算費米面就要這樣無果而終.....?

貼出：最近一次的nscf部分的輸入文件及輸出文件如下，期待高手斧正：

輸入：
##################( 3.non-scf calculation is performed to obtain the  file:$name.nscf.out########################
#  to obtain eigenvalues (saving to $name.nscf.out) of the kpoints which are in the file kves_$Sysname

cat > $name.nscf.in << EOF
&control
calculation='nscf',
prefix='$name'
pseudo_dir = '$PSEUDO_DIR/',
outdir='$TMP_DIR/'
/
&system
ibrav = 0,
nat=  2,
ntyp= 1,
ecutwfc = 90,
ecutrho = 360.00,
nbnd=18
/
&electrons
diagonalization='cg'
mixing_beta = 0.2
mixing_ndim = 10
conv_thr =  1.0d-8
/
ATOMIC_SPECIES
Sc 44.9559  Sc.pbe-sp-hgh.UPF

CELL_PARAMETERS {bohr}
6.270069069 0.001932758    -6.25482E-09
-3.137196645 5.431908189    -1.87645E-08
-1.25096E-08 0.000000000    9.7650248520

ATOMIC_POSITIONS {crystal}
Sc 0.333301249 0.666728758 0.250000432
Sc 0.666698751 0.333271242 0.749999568

K_POINTS
EOF

cat kvecs_$Sysname >>$name.nscf.in

$MPIDIR/mpirun -np $NP -machinefile $CURDIR/.nodelist $EXEDIR/pw.x < $name.nscf.in > $name.nscf.out

輸出：
.....

Parallel version (MPI), running on 16 processors
   R & G space division:  proc/nbgrp/npool/nimage =    16

   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...
   Reading input from standard input

   Atomic positions and unit cell read from directory:
   /home/users/spclibtt/home/yexq/shanghai-opt/pw-work/environment/tmp/Sc.save/

   Subspace diagonalization in iterative solution of the eigenvalue problem:
   a serial algorithm will be used

   Parallelization info
   --------------------
   sticks: dense  smooth    PW    G-vecs: dense smooth    PW
   Min       60    60    22                2385    2385    544
   Max       62    62    23                2390    2390    549
   Sum       967    967 367             38201 38201 8755

   bravais-lattice index    =          0
   lattice parameter (alat)  =    6.2701  a.u.
   unit-cell volume       =    332.6407 (a.u.)^3
   number of atoms/cell    =          2
   number of atomic types =          1
   number of electrons    =       22.00
   number of Kohn-Sham states=          18
   kinetic-energy cutoff    =    90.0000  Ry
charge density cutoff    =    360.0000  Ry
   Exchange-correlation    = SLA-PW-PBX-PBC ( 1 4 3 4 0)
   EXX-fraction             =       0.00

   celldm(1)= 6.270069  celldm(2)= 0.000000  celldm(3)= 0.000000
   celldm(4)= 0.000000  celldm(5)= 0.000000  celldm(6)= 0.000000

   crystal axes: (cart. coord. in units of alat)
            a(1) = ( 1.000000 0.000308 0.000000 )
            a(2) = (  -0.500345 0.866323 0.000000 )
            a(3) = ( 0.000000 0.000000 1.557403 )

   reciprocal axes: (cart. coord. in units 2 pi/alat)
            b(1) = (  0.999822  0.577447  0.000000 )
            b(2) = ( -0.000356  1.154098  0.000000 )
            b(3) = (  0.000000  0.000000  0.642095 )

   PseudoPot. # 1 for Sc read from file:
   /home/users/spclibtt/home/yexq/shanghai-opt/pw-work/pseudo/Sc.pbe-sp-hgh.UPF
   MD5 check sum: 97710d3f2ecb2a27396effe433afd68f
   Pseudo is Norm-conserving, Zval = 11.0
   Generated in analytical, separable form
   Using radial grid of 1173 points,  5 beta functions with:
            l(1) = 0
            l(2) = 0
            l(3) = 1
            l(4) = 1
            l(5) = 2

   atomic species valence mass    pseudopotential
      Sc          11.00 44.95590    Sc( 1.00)

   4 Sym. Ops., with inversion, found ( 2 have fractional translation)

Cartesian axes

   site n.    atom                positions (alat units)
      1          Sc  tau( 1) = (  -0.0002930 0.5777055 0.3893514  )
      2          Sc  tau( 2) = ( 0.4999482 0.2889262 1.1680516  )

  number of k points=  2601

   Number of k-points >= 100: set verbosity='high' to print them.

   Dense  grid: 38201 G-vectors    FFT dimensions: (  40,  40,  60)

   Largest allocated arrays    est. size (Mb)    dimensions
      Kohn-Sham Wavefunctions       0.08 Mb    ( 306, 18)
      NL pseudopotentials          0.12 Mb    ( 306, 26)
      Each V/rho on FFT grid       0.10 Mb    ( 6400)
      Each G-vector array          0.02 Mb    ( 2390)
      G-vector shells                0.02 Mb    ( 2214)
   Largest temporary arrays    est. size (Mb)    dimensions
      Each subspace H/S matrix       0.00 Mb    (  18,  18)
      Each matrix    0.01 Mb    (    26, 18)

   The potential is recalculated from file :
   /home/users/spclibtt/home/yexq/shanghai-opt/pw-work/environment/tmp/Sc.save/charge-density.dat

   Starting wfc are 0 randomized atomic wfcs

   total cpu time spent up to now is       1.8 secs

   per-process dynamical memory: 21.1 Mb

   Band Structure Calculation
   CG style diagonalization
   c_bands:  1 eigenvalues not converged
   c_bands:  1 eigenvalues not converged
   c_bands:  1 eigenvalues not converged
   c_bands:  1 eigenvalues not converged
   c_bands:  1 eigenvalues not converged
   c_bands:  1 eigenvalues not converged
   c_bands:  1 eigenvalues not converged
.............

贊一下

回復此樓

8樓2013-03-15 06:42:21

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

goldenfisher

金蟲 (著名寫手)

1ST強帖: 3
應助: 55 (初中生)
貴賓: 0.153
金幣: 5779.1
散金: 2908
紅花: 21
帖子: 1552
在線: 351.5小時
蟲號: 117764
注冊: 2005-11-26
性別: GG
專業(yè): 原子和分子物理

【答案】應助回帖

★
感謝參與，應助指數(shù) +1
WDD880227: 金幣+1, 謝謝指教 2013-03-27 10:22:40

你用的是超軟贗勢，ecutwfc = 90, ecutrho = 400.00,的設(shè)置是不對的。ecutrho至少要為ecutwfc的8-12倍。
Sc是金屬嗎？如果是金屬，必須要有smearing相關(guān)設(shè)置。

贊一下

回復此樓

9樓2013-03-15 09:50:40

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

xiaoqiu007

新蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 356.5
散金: 7
帖子: 100
在線: 20.8小時
蟲號: 707167
注冊: 2009-02-23
專業(yè): 金屬功能材料

引用回帖:

9樓: Originally posted by goldenfisher at 2013-03-15 09:50:40
你用的是超軟贗勢，ecutwfc = 90, ecutrho = 400.00,的設(shè)置是不對的。ecutrho至少要為ecutwfc的8-12倍。
Sc是金屬嗎？如果是金屬，必須要有smearing相關(guān)設(shè)置。

這是我仿照侯博的Cu的計算例子文件寫的腳步，包括ESPRESSO 5.0.2 pp/examples/example02或老版的example08計算費米面的例子文件中非自洽計算部分均無需設(shè)定smearing。Sc是金屬。
   Sc.pbe-sp-hgh.UPF是模守恒贗勢，之前用的是超軟贗勢Sc.pbe-nsp-van.UPF：ecutwfc = 90, ecutrho = 960.00。結(jié)果同上，仍然是不收斂。
   贗勢網(wǎng)址如下：
  http://www.quantum-espresso.org/ ... m_k=&origin_id=


#######################【自洽計算，產(chǎn)生$name.scf.out_$a文件，即Sc.scf.out_0文件】################################ self-consistent calculation -Non-Spin-Polarised case  #【為什么要進行自洽計算，目的是.....】

for a in 0
do
cat > $name.scf.in_$a << EOF
&control
calculation = 'scf'
restart_mode='from_scratch',
prefix='$name',
pseudo_dir = '$PSEUDO_DIR/',
outdir='$TMP_DIR/'
tstress=.t.,
tprnfor=.t.
/
&system
  ibrav=0,
  nat=2,
  ntyp=1,
  ecutwfc=90.0, ecutrho = 360.0
  occupations='smearing',
  smearing='methfessel-paxton',
  degauss=0.02
/
&electrons
conv_thr = 1.0d-8
mixing_beta = 0.7
/
#&CELL
#cell_dynamics= 'damp-pr'
# press=$a
#/

ATOMIC_SPECIES
Sc 44.9559 Sc.pbe-nsp-van.UPF

CELL_PARAMETERS {bohr}
6.254820382  0.00000000  0.00000000
-3.127410191  5.41683334  0.00000000
0.000000000  0.00000000  9.958412336

ATOMIC_POSITIONS {crystal}
Sc 0.3333333333333286 0.6666666666666714 0.2500000000000000
Sc 0.6666666666666714 0.3333333333333286 0.7500000000000000

K_POINTS {automatic}
8 8 4 0 0 0
EOF
$PW_ROOT/pw.x  < $name.scf.in_$a >$name.scf.out_$a
done

#########【準備并用kvecs_FS.x程序產(chǎn)生計算費米面要用到的密集網(wǎng)格k點，各k點坐標存儲于文件kves_$Sysname中】#############
#kves_$Sysname的定義見keves_FS.f源程序

# prepare input file $name.fs.in
#
Sysname='$name'
Calc_Type='FS'
nabc=' 16 16  8'    #【k點分割數(shù)與一般的計算相同嗎】
n_start=3
n_last=6
#
E_Fermi=`grep Fermi $name.scf.out_$a| cut-c 26-36` #【$name.scf.out_$a，由前一步自洽計算產(chǎn)生】
a1=`grep 'b(1)' $name.scf.out_$a| cut -c 24-54`
a2=`grep 'b(2)' $name.scf.out_$a| cut -c 24-54`
a3=`grep 'b(3)' $name.scf.out_$a| cut -c 24-54`

cat > kvecs_FS.in < $a1                         #$a1具有xx,xx,xx的數(shù)字格式，相當于一維向量，與nabc格式類似，見上
$a2
$a3
$nabc
$Sysname
EOF

kvecs_Fs.x < kvecs_FS.in >kvecs_FS.out

#########################【非自洽計算，產(chǎn)生文件$name.nscf.out】########################
# 進行一次非自洽計算，計算上一步產(chǎn)生的各k點的本征值；

cat > $name.nscf.in << EOF
&control
calculation='nscf',
prefix='$name'
pseudo_dir = '$PSEUDO_DIR/',
outdir='$TMP_DIR/'
/
&system
ibrav = 0,
nat=  2,
ntyp= 1,
ecutwfc = 90,
ecutrho = 400.00,
#  occupations='smearing', #【非自洽計算這些參數(shù)均無必要】
  # smearing='methfessel-paxton',
  # degauss=0.02
nbnd=8  #【一定要定義nbnd嗎？默認值是.....】
/
&electrons
  diagonalization='cg'
mixing_beta = 0.7
conv_thr =  1.0d-10       #【1.0d-10與1.0e-10含義一樣嗎】
/
#【diagonalization='cg' 英文注釋：conjugate-gradient-like band-by-band diagonalization. Typically slower than 'david' but it uses less memory    and is more robust (it seldom fails)】

#&CELL
#cell_dynamics= 'damp-pr'
# press=$a
#/

ATOMIC_SPECIES
Sc 44.9559 Sc.pbe-nsp-van.UPF

CELL_PARAMETERS {bohr}
6.254820382  0.00000000  0.00000000
-3.127410191  5.41683334  0.00000000
0.000000000  0.00000000  9.958412336

ATOMIC_POSITIONS {crystal}
Sc 0.3333333333333286 0.6666666666666714 0.2500000000000000
Sc 0.6666666666666714 0.3333333333333286 0.7500000000000000
K_POINTS #【K_POINTS已在上一步產(chǎn)生】
EOF

cat kvecs_$Sysname >>$name.nscf.in 【將上一步產(chǎn)生的各k點坐標，附加到$name.nscf.in 文件后面】

$PW_ROOT/pw.x< $name.nscf.in >$name.nscf.out

##########【采用bands_FS.x程序?qū)⑸弦徊降玫降膋點及相應本征值轉(zhuǎn)化成xcrysden能讀取的bxsf格式】##########################

#【下面是即將見證奇跡的時刻-采用bands_FS.x將計算的k點以及相應的本征值，轉(zhuǎn)換成xcrysden軟件的bxsf格式，以便采用xcrysden來畫圖。】

# prepare input data (input_FS, Bands.out)for bands_FS

mv $name.nscf.out  Bands_NSP.out
#【注意這利用了非自洽計算產(chǎn)生的$name.nscf.out文件，該步的意義可能是bands_FS.x陳旭默認讀取文件名為 Bands_NSP.out?】

cat > input_FS < $n_start $n_last
$E_Fermi
$Sysname
$nabc
$a1
$a2
$a3
EOF

$PW_ROOT/bands_FS.x < Bands_NSP.out >bands_fs.out #【這一操作將產(chǎn)生Bands_FS.bxsf文件】
mv Bands_FS.bxsf  $name.fs_NSP.bxsf #【這一操作意義不大，是為了記憶方便】

$ECHO "  Fermi surface plot: use 'xcrysden --bxsf $name.fs_NSP.bxsf' to plot ...\c"
$ECHO" done"

#######################################################

ESPRESSO 5.0.2 pp/examples/example02解讀

贊一下

回復此樓

10樓2013-03-15 14:29:34

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

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

返回列表

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

[考研] 310求調(diào)劑 +13	爭取九點睡 2026-03-30	13/650	2026-04-02 14:04 by asdfzly
[考研] 找調(diào)劑 +6	楚喬喬 2026-04-01	6/300	2026-04-02 12:06 by xiaoranmu
[考研] 319求調(diào)劑 +13	太容易1018 2026-04-01	13/650	2026-04-02 11:39 by Sammy2
[考研] 266分，一志愿電氣工程，本科材料，求材料專業(yè)調(diào)劑 +10	哇呼哼呼哼 2026-04-01	11/550	2026-04-02 11:31 by lnilvy
[考研] 考研調(diào)劑0702 +3	梅子味晚霞 2026-03-28	3/150	2026-04-02 10:10 by jp9609
[考研] 重慶大學材料與化工085600，初試370+，求求調(diào)劑建議 +5	shzhou_ 2026-04-01	6/300	2026-04-02 10:05 by olim
[考研] 302求調(diào)劑一志愿北航070300，本科鄭大化學 +8	圣日耳曼條 2026-04-01	11/550	2026-04-02 07:40 by chemdavid
[考研] 266分，一志愿電氣工程，本科材料，求材料專業(yè)調(diào)劑 +10	哇呼哼呼哼 2026-04-01	11/550	2026-04-01 21:48 by chyhaha
[考研] 353求調(diào)劑 +4	拉鉤不許變 2026-04-01	4/200	2026-04-01 18:10 by 記事本2026
[考研] 350求調(diào)劑 +7	阿佳～ 2026-03-31	7/350	2026-04-01 16:12 by yanflower7133
[考研] 291求調(diào)劑 +3	迷蒙木木 2026-04-01	4/200	2026-04-01 11:07 by 逆水乘風
[考研] 0855機械初試280求調(diào)劑 +3	kazenotori 2026-03-31	3/150	2026-04-01 10:08 by JourneyLucky
[考研] 346求調(diào)劑一志愿070303有機化學 +11	蘿卜燉青菜 2026-03-28	12/600	2026-03-31 17:54 by xhai2011
[考研] 085600材料與化工調(diào)劑 +16	kikiki7 2026-03-30	16/800	2026-03-31 10:03 by 氯化亞硝酰
[考研] 一志愿中海洋320化學工程與技術(shù)學碩求調(diào)劑 +8	披星河 2026-03-30	8/400	2026-03-31 08:53 by lbsjt
[有機交流] 考研調(diào)劑 +8	watb 2026-03-26	8/400	2026-03-30 18:40 by 544594351
[考研] 一志愿南開大學0710生物學359求調(diào)劑 +5	兔兔兔111223314 2026-03-29	7/350	2026-03-30 18:29 by 兔兔兔111223314
[考研] 085600，材料與化工321分求調(diào)劑 +10	大饞小子 2026-03-28	10/500	2026-03-29 23:35 by 飛行日記西
[考研] 305求調(diào)劑 +8	RuiFairyrui 2026-03-28	8/400	2026-03-29 08:22 by fmesaito
[考研] 331環(huán)境科學與工程求調(diào)劑 +3	熠然好運氣 2026-03-27	3/150	2026-03-28 04:11 by fmesaito

24小時熱門版塊排行榜

[求助] Problem - nscf - c_bands: eigenvalues not converged

» 收錄本帖的淘帖專輯推薦

» 猜你喜歡

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

【答案】應助回帖