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

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

24小時(shí)熱門版塊排行榜

返回列表

下一頁(yè)

hanyanli0475

銅蟲 (小有名氣)

應(yīng)助: 11 (小學(xué)生)
金幣: 173.1
散金: 32
紅花: 3
帖子: 296
在線: 130.5小時(shí)
蟲號(hào): 1160645
注冊(cè): 2010-12-02
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

[求助] LDA+U遇到的一個(gè)奇怪的問(wèn)題？？？

在利用LDA+U計(jì)算InN的能帶結(jié)構(gòu)時(shí)，當(dāng)把 Hubbard_U(1)=4.0,這樣設(shè)置時(shí)，就會(huì)報(bào)錯(cuò)，錯(cuò)誤已經(jīng)在下面給出，當(dāng)設(shè)置&CONTROLHubbard_U(2)=4.0,就是計(jì)算，我就暈了，好像不能給In加U,為什么呢？請(qǐng)高手們指點(diǎn)一下唄，現(xiàn)在很迷惑。
以下是我的輸入文件，報(bào)錯(cuò)信息

calculation = 'scf' ,
restart_mode = 'from_scratch' ,
outdir = '/home/zhao/6.77/tmp/' ,
pseudo_dir = '/opt/espresso-4.2.1/pseudo/' ,
prefix = 'InN' ,
tprnfor = .true. ,
tstress = .true. ,
/
&SYSTEM
ibrav = 4,celldm(1) = 6.796,celldm(3) = 1.618,nat = 4,ntyp = 2,
ecutwfc = 45.0 ,ecutrho = 450.0 ,
occupations = 'smearing' ,degauss = 0.02 ,smearing = 'mp' ,
lda_plus_u =.true. Hubbard_U(1)=4.0,
/
&ELECTRONS
conv_thr = 1.0d-8 ,
mixing_beta = 0.7 ,
diagonalization = 'david',
startingpot = 'atomic',
startingwfc = 'atomic+random',
mixing_mode = 'plain',
/
ATOMIC_SPECIES
In 114.820  In.pz-d-rrkjus.UPF
N 14.007 N.pz-rrkjus.UPF
ATOMIC_POSITIONS crystal
In    0.333333333 0.666666667 0.002089135
In    0.666666667 0.333333333 0.502107613
N       0.333333333 0.666666667 0.380470108
N       0.666666667 0.333333333 0.880533144

K_POINTS automatic
4 4 4    0 0 0

錯(cuò)誤：
file N.pz-rrkjus.UPF: wavefunction(s)  2S renormalized

psd = In

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
   from set_hubbard_l : error #       1
   pseudopotential not yet inserted
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

   stopping ...

回復(fù)此樓

» 猜你喜歡

299求調(diào)劑已經(jīng)有6人回復(fù)
286分人工智能專業(yè)請(qǐng)求調(diào)劑愿意跨考！已經(jīng)有5人回復(fù)
279求調(diào)劑已經(jīng)有3人回復(fù)
336求調(diào)劑已經(jīng)有3人回復(fù)
286求調(diào)劑已經(jīng)有3人回復(fù)
一志愿華中科技大學(xué)071000，求調(diào)劑已經(jīng)有3人回復(fù)
279分求調(diào)劑一志愿211 已經(jīng)有14人回復(fù)
332求調(diào)劑已經(jīng)有3人回復(fù)
0856材料專碩353求調(diào)劑已經(jīng)有3人回復(fù)
材料學(xué)學(xué)碩080502 337求調(diào)劑-一志愿華中科技大學(xué) 已經(jīng)有5人回復(fù)

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

計(jì)算晶體的能帶譜圖發(fā)現(xiàn)了一個(gè)很奇怪的現(xiàn)象已經(jīng)有9人回復(fù)
HPLC定量問(wèn)題－很奇怪已經(jīng)有24人回復(fù)
畫能帶時(shí)遇到奇怪的問(wèn)題---幫忙看下吧已經(jīng)有14人回復(fù)
一個(gè)奇怪的現(xiàn)象已經(jīng)有47人回復(fù)
關(guān)于穩(wěn)定結(jié)構(gòu)的計(jì)算，我遇到一個(gè)奇怪的現(xiàn)象已經(jīng)有3人回復(fù)
關(guān)于液相的一個(gè)奇怪問(wèn)題，請(qǐng)各位高手指教~~~~ 已經(jīng)有30人回復(fù)
做液相遇到一個(gè)奇怪的現(xiàn)象，基線會(huì)忽然規(guī)律波動(dòng)起來(lái)？？？已經(jīng)有3人回復(fù)
發(fā)現(xiàn)一個(gè)特別奇怪的問(wèn)題！已經(jīng)有37人回復(fù)
【求助】Segmentation fault 已經(jīng)有15人回復(fù)
【求助/交流】連接轉(zhuǎn)化時(shí)遇到的奇怪問(wèn)題求解已經(jīng)有12人回復(fù)
【求助】向各位高手請(qǐng)教一個(gè)奇怪的問(wèn)題已經(jīng)有25人回復(fù)

什么時(shí)候能算完！�。�！

1樓 2012-06-09 13:24:52

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

ChemiAndy

木蟲 (正式寫手)

1ST強(qiáng)帖: 4
應(yīng)助: 52 (初中生)
金幣: 2573.3
散金: 18
紅花: 104
帖子: 590
在線: 553.8小時(shí)
蟲號(hào): 836907
注冊(cè): 2009-08-30
性別: GG
專業(yè): 理論和計(jì)算化學(xué)

為毛在我這里運(yùn)行的好好的？
我只是修改了路徑設(shè)置，并加上了你復(fù)制時(shí)候漏掉的第一行“&CONTROL“

贊一下

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

hanyanli0475

2樓2012-06-10 06:03:47

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

ChemiAndy

木蟲 (正式寫手)

1ST強(qiáng)帖: 4
應(yīng)助: 52 (初中生)
金幣: 2573.3
散金: 18
紅花: 104
帖子: 590
在線: 553.8小時(shí)
蟲號(hào): 836907
注冊(cè): 2009-08-30
性別: GG
專業(yè): 理論和計(jì)算化學(xué)

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

恩，最后有點(diǎn)錯(cuò)誤，可能是code的問(wèn)題：

CODE:

     Program PWSCF v.5.0        starts on  9Jun2012 at 18: 0: 7 

     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     2 processors

     R & G space division:  proc/pool =    2

     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

               file N.pz-rrkjus.UPF: wavefunction(s)  2S renormalized

     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         708     291     81                35404     9014    1435

     Max         709     292     82                35411     9015    1440

     Sum        1417     583    163                70815    18029    2875

     bravais-lattice index     =            4

     lattice parameter (alat)  =       6.7960  a.u.

     unit-cell volume          =     439.8142 (a.u.)^3

     number of atoms/cell      =            4

     number of atomic types    =            2

     number of electrons       =        36.00

     number of Kohn-Sham states=           22

     kinetic-energy cutoff     =      45.0000  Ry

     charge density cutoff     =     450.0000  Ry

     convergence threshold     =      1.0E-06

     mixing beta               =       0.7000

     number of iterations used =            8  plain     mixing

     Exchange-correlation      =  SLA  PZ   NOGX NOGC ( 1 1 0 0 0)

     EXX-fraction              =        0.00

     celldm(1)=   6.796000  celldm(2)=   0.000000  celldm(3)=   1.618000

     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.618000 )  

     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.618047 )  

     PseudoPot. # 1 for In read from file:

     /home/xijun/espresso/pseudo/In.pz-d-rrkjus.UPF

     MD5 check sum: ec4b7502086f46f0fa2ea5856329c8f0

     Pseudo is Ultrasoft + core correction, Zval = 13.0

     Generated by new atomic code, or converted to UPF format

     Using radial grid of 1241 points,  4 beta functions with: 

                l(1) =   2

                l(2) =   2

                l(3) =   0

                l(4) =   0

     Q(r) pseudized with 0 coefficients 

     PseudoPot. # 2 for N  read from file:

     /home/xijun/espresso/pseudo/N.pz-rrkjus.UPF

     MD5 check sum: 9152467a3ef5f298703596eb0e9fbbb1

     Pseudo is Ultrasoft, Zval =  5.0

     Generated by new atomic code, or converted to UPF format

     Using radial grid of 1257 points,  4 beta functions with: 

                l(1) =   0

                l(2) =   0

                l(3) =   1

                l(4) =   1

     Q(r) pseudized with 0 coefficients 

     atomic species   valence    mass     pseudopotential

        In            13.00   114.82000     In( 1.00)

        N              5.00    14.00700     N ( 1.00)

 -----

     simplified LDA+U calculation, Hubbard_lmax = 2

     atomic species  L   Hubbard U  Hubbard alpha

        In           2    4.000000    0.000000

 -----

      6 Sym. Ops. (no inversion) found

   Cartesian axes

     site n.     atom                  positions (alat units)

         1           In  tau(   1) = (  -0.0000000   0.5773503   0.0033802  )

         2           In  tau(   2) = (   0.5000000   0.2886751   0.8124101  )

         3           N   tau(   3) = (  -0.0000000   0.5773503   0.6156006  )

         4           N   tau(   4) = (   0.5000000   0.2886751   1.4247026  )

     number of k points=    13  Methfessel-Paxton smearing, width (Ry)=  0.0200

                       cart. coord. in units 2pi/alat

        k(    1) = (   0.0000000   0.0000000   0.0000000), wk =   0.0312500

        k(    2) = (   0.0000000   0.0000000   0.1545117), wk =   0.0625000

        k(    3) = (   0.0000000   0.0000000  -0.3090235), wk =   0.0312500

        k(    4) = (   0.0000000   0.2886751   0.0000000), wk =   0.1875000

        k(    5) = (   0.0000000   0.2886751   0.1545117), wk =   0.1875000

        k(    6) = (   0.0000000   0.2886751  -0.3090235), wk =   0.1875000

        k(    7) = (   0.0000000  -0.5773503   0.0000000), wk =   0.0937500

        k(    8) = (   0.0000000  -0.5773503   0.1545117), wk =   0.1875000

        k(    9) = (   0.0000000  -0.5773503  -0.3090235), wk =   0.0937500

        k(   10) = (   0.2500000   0.4330127   0.0000000), wk =   0.1875000

        k(   11) = (   0.2500000   0.4330127   0.1545117), wk =   0.3750000

        k(   12) = (   0.2500000   0.4330127  -0.3090235), wk =   0.1875000

        k(   13) = (   0.0000000  -0.2886751   0.1545117), wk =   0.1875000

     Dense  grid:    70815 G-vectors     FFT dimensions: (  45,  45,  75)

     Smooth grid:    18029 G-vectors     FFT dimensions: (  30,  30,  48)

     Largest allocated arrays     est. size (Mb)     dimensions

        Kohn-Sham Wavefunctions         0.38 Mb     (   1137,   22)

        Atomic wavefunctions            0.45 Mb     (   1137,   26)

        NL pseudopotentials             0.69 Mb     (   1137,   40)

        Each V/rho on FFT grid          1.17 Mb     (  76950)

        Each G-vector array             0.27 Mb     (  35404)

        G-vector shells                 0.02 Mb     (   3159)

     Largest temporary arrays     est. size (Mb)     dimensions

        Auxiliary wavefunctions         1.53 Mb     (   1137,   88)

        Each subspace H/S matrix        0.12 Mb     (  88,  88)

        Each  matrix      0.01 Mb     (     40,   22)

        Arrays for rho mixing           9.39 Mb     (  76950,   8)

     Check: negative/imaginary core charge=   -0.000009    0.000000

     Initial potential from superposition of free atoms

     starting charge   35.99450, renormalised to   36.00000

     Number of +U iterations with fixed ns =  0

     Starting occupations:

 --- enter write_ns ---

 LDA+U parameters:

U( 1)     =  4.00000000

alpha( 1) =  0.00000000

atom    1   Tr[ns(na)] =  10.00000

    eigenvalues: 

  1.000  1.000  1.000  1.000  1.000

    eigenvectors:

  1.000  0.000  0.000  0.000  0.000

  0.000  1.000  0.000  0.000  0.000

  0.000  0.000  1.000  0.000  0.000

  0.000  0.000  0.000  1.000  0.000

  0.000  0.000  0.000  0.000  1.000

    occupations:

  1.000  0.000  0.000  0.000  0.000

  0.000  1.000  0.000  0.000  0.000

  0.000  0.000  1.000  0.000  0.000

  0.000  0.000  0.000  1.000  0.000

  0.000  0.000  0.000  0.000  1.000

atom    2   Tr[ns(na)] =  10.00000

    eigenvalues: 

  1.000  1.000  1.000  1.000  1.000

    eigenvectors:

  1.000  0.000  0.000  0.000  0.000

  0.000  1.000  0.000  0.000  0.000

  0.000  0.000  1.000  0.000  0.000

  0.000  0.000  0.000  1.000  0.000

  0.000  0.000  0.000  0.000  1.000

    occupations:

  1.000  0.000  0.000  0.000  0.000

  0.000  1.000  0.000  0.000  0.000

  0.000  0.000  1.000  0.000  0.000

  0.000  0.000  0.000  1.000  0.000

  0.000  0.000  0.000  0.000  1.000

N of occupied +U levels =  20.0000000

 --- exit write_ns ---

 Atomic wfc used for LDA+U Projector are NOT orthogonalized

     Starting wfc are   26 randomized atomic wfcs

     total cpu time spent up to now is       29.0 secs

     per-process dynamical memory:    21.8 Mb

     Self-consistent Calculation

     iteration #  1     ecut=    45.00 Ry     beta=0.70

     Davidson diagonalization with overlap

     ethr =  1.00E-02,  avg # of iterations =  2.2

 --- enter write_ns ---

 LDA+U parameters:

U( 1)     =  4.00000000

alpha( 1) =  0.00000000

atom    1   Tr[ns(na)] =   9.99375

    eigenvalues: 

  0.999  0.999  0.999  1.000  1.000

    eigenvectors:

  1.000  0.000  0.000  0.000  0.000

  0.000  0.296  0.033  0.145  0.526

  0.000  0.033  0.296  0.526  0.145

  0.000  0.068  0.603  0.258  0.071

  0.000  0.603  0.068  0.071  0.258

    occupations:

  0.999 -0.000 -0.000 -0.000 -0.000

 -0.000  1.000 -0.000  0.000  0.000

 -0.000 -0.000  1.000  0.000 -0.000

 -0.000  0.000  0.000  0.999  0.000

 -0.000  0.000 -0.000  0.000  0.999

atom    2   Tr[ns(na)] =   9.99375

    eigenvalues: 

  0.999  0.999  0.999  1.000  1.000

    eigenvectors:

  1.000  0.000  0.000  0.000  0.000

  0.000  0.218  0.115  0.667  0.000

  0.000  0.115  0.218  0.000  0.667

  0.000  0.231  0.437  0.000  0.333

  0.000  0.437  0.231  0.333  0.000

    occupations:

  0.999 -0.000 -0.000 -0.000 -0.000

 -0.000  1.000 -0.000  0.000 -0.000

 -0.000 -0.000  1.000 -0.000 -0.000

 -0.000  0.000 -0.000  0.999  0.000

 -0.000 -0.000 -0.000  0.000  0.999

N of occupied +U levels =  19.9875049

 --- exit write_ns ---

     total cpu time spent up to now is       34.0 secs

     total energy              =    -311.96202638 Ry

     Harris-Foulkes estimate   =    -312.89328229 Ry

     estimated scf accuracy    <       1.21206838 Ry

     iteration #  2     ecut=    45.00 Ry     beta=0.70

     Davidson diagonalization with overlap

     ethr =  3.37E-03,  avg # of iterations =  3.2

     total cpu time spent up to now is       38.9 secs

     total energy              =    -312.21161943 Ry

     Harris-Foulkes estimate   =    -313.12944003 Ry

     estimated scf accuracy    <       2.09038321 Ry

     iteration #  3     ecut=    45.00 Ry     beta=0.70

     Davidson diagonalization with overlap

     ethr =  3.37E-03,  avg # of iterations =  3.0

     total cpu time spent up to now is       51.8 secs

     total energy              =    -312.57861633 Ry

     Harris-Foulkes estimate   =    -312.57990562 Ry

     estimated scf accuracy    <       0.00351443 Ry

     iteration #  4     ecut=    45.00 Ry     beta=0.70

     Davidson diagonalization with overlap

     ethr =  9.76E-06,  avg # of iterations =  3.6

     total cpu time spent up to now is       57.2 secs

     total energy              =    -312.57917939 Ry

     Harris-Foulkes estimate   =    -312.57945786 Ry

     estimated scf accuracy    <       0.00068670 Ry

     iteration #  5     ecut=    45.00 Ry     beta=0.70

     Davidson diagonalization with overlap

     ethr =  1.91E-06,  avg # of iterations =  2.0

     total cpu time spent up to now is       60.9 secs

     total energy              =    -312.57921637 Ry

     Harris-Foulkes estimate   =    -312.57923000 Ry

     estimated scf accuracy    <       0.00003665 Ry

     iteration #  6     ecut=    45.00 Ry     beta=0.70

     Davidson diagonalization with overlap

     ethr =  1.02E-07,  avg # of iterations =  2.8

     total cpu time spent up to now is       65.2 secs

     total energy              =    -312.57922451 Ry

     Harris-Foulkes estimate   =    -312.57922516 Ry

     estimated scf accuracy    <       0.00000158 Ry

     iteration #  7     ecut=    45.00 Ry     beta=0.70

     Davidson diagonalization with overlap

     ethr =  4.40E-09,  avg # of iterations =  3.0

     total cpu time spent up to now is       85.9 secs

     total energy              =    -312.57922467 Ry

     Harris-Foulkes estimate   =    -312.57922533 Ry

     estimated scf accuracy    <       0.00000204 Ry

     iteration #  8     ecut=    45.00 Ry     beta=0.70

     Davidson diagonalization with overlap

     ethr =  4.40E-09,  avg # of iterations =  2.5

     total cpu time spent up to now is       90.2 secs

     End of self-consistent calculation

 --- enter write_ns ---

 LDA+U parameters:

U( 1)     =  4.00000000

alpha( 1) =  0.00000000

atom    1   Tr[ns(na)] =   9.99457

    eigenvalues: 

  0.999  0.999  0.999  1.000  1.000

    eigenvectors:

  1.000  0.000  0.000  0.000  0.000

  0.000  0.291  0.033  0.138  0.538

  0.000  0.033  0.291  0.538  0.138

  0.000  0.069  0.607  0.258  0.066

  0.000  0.607  0.069  0.066  0.258

    occupations:

  0.999 -0.000 -0.000 -0.000 -0.000

 -0.000  1.000 -0.000  0.000  0.000

 -0.000 -0.000  1.000  0.000 -0.000

 -0.000  0.000  0.000  0.999  0.000

 -0.000  0.000 -0.000  0.000  0.999

atom    2   Tr[ns(na)] =   9.99457

    eigenvalues: 

  0.999  0.999  0.999  1.000  1.000

    eigenvectors:

  1.000  0.000  0.000  0.000  0.000

  0.000  0.214  0.111  0.673  0.003

  0.000  0.111  0.214  0.003  0.673

  0.000  0.230  0.445  0.001  0.323

  0.000  0.445  0.230  0.323  0.001

    occupations:

  0.999 -0.000 -0.000 -0.000 -0.000

 -0.000  1.000 -0.000  0.000 -0.000

 -0.000 -0.000  1.000 -0.000 -0.000

 -0.000  0.000 -0.000  0.999  0.000

 -0.000 -0.000 -0.000  0.000  0.999

N of occupied +U levels =  19.9891406

 --- exit write_ns ---

          k = 0.0000 0.0000 0.0000 (  2273 PWs)   bands (ev):

    -7.3732  -6.7012  -6.6830  -6.6830  -6.6253  -6.6253  -6.4556  -6.4556

    -6.3059  -6.3059  -5.5012  -3.0986   2.9317   7.7145   7.7145   8.1880

     8.5229   8.5453   8.5453  11.0606  15.1946  17.1656

          k = 0.0000 0.0000 0.1545 (  2261 PWs)   bands (ev):

    -7.2968  -6.8424  -6.6834  -6.6834  -6.6410  -6.6410  -6.4249  -6.4249

    -6.3208  -6.3208  -5.1904  -3.5603   3.6377   6.8119   7.8324   7.8324

     8.4202   8.4202   9.7463  11.0565  15.0711  16.6442

          k = 0.0000 0.0000-0.3090 (  2232 PWs)   bands (ev):

    -7.0949  -7.0941  -6.6697  -6.6697  -6.6694  -6.6694  -6.3655  -6.3655

    -6.3654  -6.3654  -4.4460  -4.4448   5.1679   5.1684   8.1222   8.1222

     8.1232   8.1232  10.7362  10.7370  15.3032  15.3059

          k = 0.0000 0.2887 0.0000 (  2239 PWs)   bands (ev):

    -7.4260  -7.0717  -6.6339  -6.6160  -6.5391  -6.5002  -6.4932  -6.4287

    -6.3928  -6.2927  -4.5562  -2.7425   3.5096   5.3049   5.9152   6.9778

     7.9983   8.1394  10.7783  12.8514  14.3961  15.3742

          k = 0.0000 0.2887 0.1545 (  2232 PWs)   bands (ev):

    -7.3721  -7.1185  -6.6380  -6.5998  -6.5714  -6.5169  -6.4582  -6.4097

    -6.3984  -6.3129  -4.3536  -3.0978   3.8704   5.2443   5.9952   7.1358

     7.8518   7.9554  11.0119  12.6997  14.3241  15.1471

          k = 0.0000 0.2887-0.3090 (  2234 PWs)   bands (ev):

    -7.2431  -7.2425  -6.6152  -6.6150  -6.5591  -6.5589  -6.4208  -6.4207

    -6.3616  -6.3615  -3.7961  -3.7950   4.5772   4.5781   6.9760   6.9767

     7.5321   7.5333  11.8831  11.8843  14.6551  14.6566

          k = 0.0000-0.5774 0.0000 (  2256 PWs)   bands (ev):

    -7.4939  -7.3415  -6.6290  -6.6185  -6.5159  -6.4657  -6.3863  -6.3835

    -6.3158  -6.2757  -3.2163  -2.6520   3.2819   4.1209   5.2722   6.3573

     6.9147   7.7285  12.1377  12.7394  13.6133  14.5958

          k = 0.0000-0.5774 0.1545 (  2234 PWs)   bands (ev):

    -7.4612  -7.3538  -6.6489  -6.5807  -6.5572  -6.4539  -6.4104  -6.3367

    -6.3152  -6.2930  -3.1548  -2.7550   3.1327   3.7326   5.9496   6.5404

     7.2695   7.5061  11.5713  12.5582  14.2917  14.4878

          k = 0.0000-0.5774-0.3090 (  2252 PWs)   bands (ev):

    -7.3969  -7.3964  -6.6415  -6.6414  -6.4908  -6.4906  -6.3592  -6.3591

    -6.3107  -6.3106  -2.9765  -2.9758   3.2280   3.2291   7.0032   7.0047

     7.0519   7.0521  11.7454  11.7459  14.6830  14.6852

          k = 0.2500 0.4330 0.0000 (  2247 PWs)   bands (ev):

    -7.4466  -7.3108  -6.6840  -6.5127  -6.5014  -6.4746  -6.4421  -6.4344

    -6.3448  -6.3097  -3.3372  -2.4923   3.7952   4.2515   5.3367   5.4586

     6.9663   7.0976  12.7954  13.5943  14.3580  14.6997

          k = 0.2500 0.4330 0.1545 (  2238 PWs)   bands (ev):

    -7.4165  -7.3183  -6.6382  -6.5709  -6.5201  -6.4797  -6.4357  -6.4043

    -6.3434  -6.3243  -3.2408  -2.6445   3.6077   4.0148   5.4488   6.1609

     6.6530   7.4107  12.6295  13.4633  14.5491  14.6991

          k = 0.2500 0.4330-0.3090 (  2232 PWs)   bands (ev):

    -7.3564  -7.3559  -6.5788  -6.5788  -6.5297  -6.5295  -6.4321  -6.4320

    -6.3247  -6.3246  -2.9708  -2.9701   3.6431   3.6441   5.9635   5.9648

     7.2488   7.2491  12.8298  12.8318  14.6667  14.6671

          k = 0.0000-0.2887 0.1545 (  2232 PWs)   bands (ev):

    -7.3721  -7.1185  -6.6381  -6.5997  -6.5713  -6.5170  -6.4583  -6.4096

    -6.3984  -6.3129  -4.3536  -3.0977   3.8704   5.2441   5.9957   7.1358

     7.8517   7.9554  11.0116  12.6997  14.3254  15.1460

     the Fermi energy is     8.5915 ev

!    total energy              =    -312.57922491 Ry

     Harris-Foulkes estimate   =    -312.57922500 Ry

     estimated scf accuracy    <       0.00000013 Ry

     The total energy is the sum of the following terms:

     one-electron contribution =     -79.18311959 Ry

     hartree contribution      =      68.63097517 Ry

     xc contribution           =     -94.84178571 Ry

     ewald contribution        =    -207.18697592 Ry

     Hubbard energy            =       0.00159519 Ry

     smearing contrib. (-TS)   =       0.00008596 Ry

     convergence has been achieved in   8 iterations

     Forces acting on atoms (Ry/au):

     atom    1 type  1   force =     0.00000000    0.00000000    0.00047308

     atom    2 type  1   force =     0.00000000    0.00000000    0.00021429

     atom    3 type  2   force =     0.00000000    0.00000000   -0.00024660

     atom    4 type  2   force =     0.00000000   -0.00000000   -0.00044077

     Total force =     0.000724     Total SCF correction =     0.000647

     SCF correction compared to forces is large: reduce conv_thr to get better values

     entering subroutine stress ...

贊一下

回復(fù)此樓

3樓2012-06-10 06:06:04

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

ChemiAndy

木蟲 (正式寫手)

1ST強(qiáng)帖: 4
應(yīng)助: 52 (初中生)
金幣: 2573.3
散金: 18
紅花: 104
帖子: 590
在線: 553.8小時(shí)
蟲號(hào): 836907
注冊(cè): 2009-08-30
性別: GG
專業(yè): 理論和計(jì)算化學(xué)

【答案】應(yīng)助回帖

★ ★ ★ ★ ★
hanyanli0475: 金幣+5, ★有幫助 2012-06-10 08:04:39

內(nèi)容已刪除

贊一下

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

hanyanli0475

4樓2012-06-10 06:08:58

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

hanyanli0475

銅蟲 (小有名氣)

應(yīng)助: 11 (小學(xué)生)
金幣: 173.1
散金: 32
紅花: 3
帖子: 296
在線: 130.5小時(shí)
蟲號(hào): 1160645
注冊(cè): 2010-12-02
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

送鮮花一朵

引用回帖:

4樓: Originally posted by ChemiAndy at 2012-06-10 06:08:58
我估計(jì)你是pseudopotential文件的問(wèn)題
我是從QE網(wǎng)站新鮮下載的
wget http://www.quantum-espresso.org/wp-content/uploads/upf_files/In.pz-d-rrkjus.UPF

謝謝您的回答，但是，我不加U時(shí)，這個(gè)贗勢(shì)是沒(méi)有問(wèn)題，我試過(guò)你給我的這個(gè)贗勢(shì)，也相同的報(bào)錯(cuò)，我在網(wǎng)上搜這個(gè)問(wèn)題，有人回答說(shuō)是：
set_hubbard_l.f90設(shè)置問(wèn)題，我照著他說(shuō)的寫法修改，還是不行，您能再幫我想想嗎，呵呵！謝謝哈！

贊一下

回復(fù)此樓

什么時(shí)候能算完�。。�！

5樓2012-06-10 08:04:26

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

hanyanli0475

銅蟲 (小有名氣)

應(yīng)助: 11 (小學(xué)生)
金幣: 173.1
散金: 32
紅花: 3
帖子: 296
在線: 130.5小時(shí)
蟲號(hào): 1160645
注冊(cè): 2010-12-02
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

送鮮花一朵

引用回帖:

2樓: Originally posted by ChemiAndy at 2012-06-10 06:03:47
為毛在我這里運(yùn)行的好好的？
我只是修改了路徑設(shè)置，并加上了你復(fù)制時(shí)候漏掉的第一行“&CONTROL“

我查了一下set_hubbard_l.f90，這個(gè)里面在過(guò)渡金屬里面沒(méi)有In，你看一下，你的set_hubbard_l.f90中，是不是有In�。课疫@個(gè)沒(méi)有，怎么加呢？請(qǐng)您幫幫忙吧，現(xiàn)在很糾結(jié)啊！

贊一下

回復(fù)此樓

什么時(shí)候能算完�。。�！

6樓2012-06-10 09:16:12

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

ChemiAndy

木蟲 (正式寫手)

1ST強(qiáng)帖: 4
應(yīng)助: 52 (初中生)
金幣: 2573.3
散金: 18
紅花: 104
帖子: 590
在線: 553.8小時(shí)
蟲號(hào): 836907
注冊(cè): 2009-08-30
性別: GG
專業(yè): 理論和計(jì)算化學(xué)

【答案】應(yīng)助回帖

★ ★ ★ ★ ★
hanyanli0475: 金幣+5, ★★★★★最佳答案 2012-06-10 14:04:22

1. 你修改set_hubbard_l.f90 源文件后是否重新編譯了？
2. 你現(xiàn)在用的版本是多少？我用的是最新版5.0.

贊一下

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

hanyanli0475

7樓2012-06-10 09:51:25

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

ChemiAndy

木蟲 (正式寫手)

1ST強(qiáng)帖: 4
應(yīng)助: 52 (初中生)
金幣: 2573.3
散金: 18
紅花: 104
帖子: 590
在線: 553.8小時(shí)
蟲號(hào): 836907
注冊(cè): 2009-08-30
性別: GG
專業(yè): 理論和計(jì)算化學(xué)

我的set_hubbard_l.f90里面也沒(méi)有'In'
奇怪，為什么我的能跑呢？

贊一下

回復(fù)此樓

8樓2012-06-10 09:54:28

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

ChemiAndy

木蟲 (正式寫手)

1ST強(qiáng)帖: 4
應(yīng)助: 52 (初中生)
金幣: 2573.3
散金: 18
紅花: 104
帖子: 590
在線: 553.8小時(shí)
蟲號(hào): 836907
注冊(cè): 2009-08-30
性別: GG
專業(yè): 理論和計(jì)算化學(xué)

【答案】應(yīng)助回帖

★
xueht987: 金幣+1, 應(yīng)助指數(shù)+1, 謝謝回帖，鼓勵(lì)交流！ 2012-06-10 10:50:30

問(wèn)題找到了
5.0版的pwscf中的set_hubbard_l.f90里面是包含In這個(gè)元素的。剛才我查的那個(gè)文件版本不對(duì)。

你下載5.0版重新編譯即可。

贊一下

回復(fù)此樓

9樓2012-06-10 10:15:12

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

ChemiAndy

木蟲 (正式寫手)

1ST強(qiáng)帖: 4
應(yīng)助: 52 (初中生)
金幣: 2573.3
散金: 18
紅花: 104
帖子: 590
在線: 553.8小時(shí)
蟲號(hào): 836907
注冊(cè): 2009-08-30
性別: GG
專業(yè): 理論和計(jì)算化學(xué)

5.0版的set_hubbard_l.f90

CODE:

!

! Copyright (C) 2001-2010 Quantum ESPRESSO group

! This file is distributed under the terms of the

! GNU General Public License. See the file `License'

! in the root directory of the present distribution,

! or http://www.gnu.org/copyleft/gpl.txt .

!

!---------------------------------------------------------------------------

FUNCTION set_hubbard_l( psd ) RESULT( hubbard_l )

  !---------------------------------------------------------------------------

  !

  USE io_global, ONLY : stdout

  !

  IMPLICIT NONE

  !

  INTEGER                      :: hubbard_l

  CHARACTER(LEN=2), INTENT(IN) :: psd

  !

  !

  SELECT CASE( TRIM(ADJUSTL(psd)) )

     !

     ! ... transition metals

     !

     CASE( 'Ti', 'V',  'Cr', 'Mn', 'Fe', 'Co', 'Ni', 'Cu', 'Zn', &

           'Zr', 'Nb', 'Mo', 'Tc', 'Ru', 'Rh', 'Pd', 'Ag', 'Cd', &

           'Hf', 'Ta', 'W',  'Re', 'Os', 'Ir', 'Pt', 'Au', 'Hg'  )

        !

        hubbard_l = 2  

        !

     !

     ! ... rare earths

     !

     CASE('Ce','Pr','Nd','Pm','Sm','Eu','Gd','Tb','Dy','Ho','Er','Tm','Yb','Lu', & 

          'Th','Pa','U', 'Np','Pu','Am','Cm','Bk','Cf','Es','Fm','Md','No','Lr' )

        !

        hubbard_l = 3

        !

     !

     ! ... other elements

     !

     CASE( 'H' )

        !

        hubbard_l =  0

        !

     CASE( 'C', 'N', 'O' )

        !

        hubbard_l =  1

        !

     CASE( 'Ga', 'In' )

        !

        hubbard_l =  2

        !

     CASE DEFAULT

        !

        hubbard_l = -1

        !

        WRITE( stdout, '(/,"psd = ",A,/)' ) psd

        !

        CALL errore( 'set_hubbard_l', 'pseudopotential not yet inserted', 1 )

        !

  END SELECT

  !

  RETURN  

  !

END FUNCTION set_Hubbard_l

贊一下

回復(fù)此樓

10樓2012-06-10 10:16:15

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

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

返回列表

下一頁(yè)

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

[考研] 0856材料專碩353求調(diào)劑 +3	NIFFFfff 2026-03-20	3/150	2026-03-21 10:23 by luoyongfeng
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +5	脫穎而出 2026-03-16	15/750	2026-03-21 10:16 by 脫穎而出
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +7	楓橋ZL 2026-03-18	9/450	2026-03-21 09:56 by JineShine
[考研] 材料工程（專）一志愿985 初試335求調(diào)劑 +3	hiloiy 2026-03-17	4/200	2026-03-21 03:04 by JourneyLucky
[考研] 265求調(diào)劑 +9	梁梁校校 2026-03-17	9/450	2026-03-21 02:17 by JourneyLucky
[考研] 278求調(diào)劑 +6	煙火先于春 2026-03-17	6/300	2026-03-21 01:57 by JourneyLucky
[考研] 求調(diào)劑一志愿南京航空航天大學(xué)289分 +3	@taotao 2026-03-19	3/150	2026-03-20 21:34 by JourneyLucky
[考研] 材料學(xué)碩318求調(diào)劑 +5	February_Feb 2026-03-19	5/250	2026-03-19 23:51 by 23Postgrad
[考研] 0703化學(xué)調(diào)劑 +4	18889395102 2026-03-18	4/200	2026-03-19 16:13 by 30660438
[考研] 286求調(diào)劑 +6	lemonzzn 2026-03-16	10/500	2026-03-19 14:31 by lemonzzn
[考研] 344求調(diào)劑 +6	knight344 2026-03-16	7/350	2026-03-18 20:13 by walc
[考研] 085601專碩，總分342求調(diào)劑，地區(qū)不限 +5	share_joy 2026-03-16	5/250	2026-03-18 14:48 by haxia
[考研] 0854，計(jì)算機(jī)類招收調(diào)劑 +3	胡辣湯放糖 2026-03-15	6/300	2026-03-18 12:09 by 上岸上岸……..
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3	Eember. 2026-03-17	9/450	2026-03-18 10:28 by Eember.
[考博] 26申博 +4	八6八68 2026-03-16	4/200	2026-03-17 13:00 by 輕松不少隨
[考研] 302求調(diào)劑 +4	小賈同學(xué)123 2026-03-15	8/400	2026-03-17 10:33 by 小賈同學(xué)123
[考研] 11408 一志愿西電，277分求調(diào)劑 +3	zhouzhen654 2026-03-16	3/150	2026-03-17 07:03 by laoshidan
[考研] 326求調(diào)劑 +4	諾貝爾化學(xué)獎(jiǎng)覬?/a> 2026-03-15	7/350	2026-03-16 17:11 by 諾貝爾化學(xué)獎(jiǎng)覬?/a>
[考研] 326求調(diào)劑 +3	mlpqaz03 2026-03-15	3/150	2026-03-16 07:33 by Iveryant
[考研] 289求調(diào)劑 +4	這么名字咋樣 2026-03-14	6/300	2026-03-14 18:58 by userper