应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » LDA+U遇到的一个奇怪的问题？？？
151/2返回列表12下一页
查看: 2691  |  回复: 14
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

hanyanli0475

铜虫 (小有名气)

[求助] LDA+U遇到的一个奇怪的问题？？？

在利用LDA+U计算InN的能带结构时，当把 Hubbard_U(1)=4.0,这样设置时，就会报错，错误已经在下面给出，当设置&CONTROLHubbard_U(2)=4.0,就是计算，我就晕了，好像不能给In加U,为什么呢？请高手们指点一下呗，现在很迷惑。
以下是我的输入文件，报错信息

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

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

psd = In


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

     stopping ...
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
什么时候能算完！！！！
1楼 2012-06-09 13:24:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ChemiAndy

木虫 (正式写手)
为毛在我这里运行的好好的？
我只是修改了路径设置，并加上了你复制时候漏掉的第一行“&CONTROL“
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

hanyanli0475
2楼2012-06-10 06:03:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ChemiAndy

木虫 (正式写手)

【答案】应助回帖

感谢参与，应助指数 +1
恩，最后有点错误，可能是code的问题：
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 ...
赞一下 回复此楼
3楼2012-06-10 06:06:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ChemiAndy

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★
hanyanli0475: 金币+5, ★有帮助 2012-06-10 08:04:39
内容已删除
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

hanyanli0475
4楼2012-06-10 06:08:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hanyanli0475

铜虫 (小有名气)
送鲜花一朵
引用回帖:
4楼: Originally posted by ChemiAndy at 2012-06-10 06:08:58
我估计你是pseudopotential文件的问题
我是从QE网站新鲜下载的
wget http://www.quantum-espresso.org/wp-content/uploads/upf_files/In.pz-d-rrkjus.UPF

谢谢您的回答，但是，我不加U时，这个贗势是没有问题，我试过你给我的这个贗势，也相同的报错，我在网上搜这个问题，有人回答说是：
set_hubbard_l.f90设置问题，我照着他说的写法修改，还是不行，您能再帮我想想吗，呵呵！谢谢哈！
赞一下 回复此楼
什么时候能算完！！！！
5楼2012-06-10 08:04:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hanyanli0475

铜虫 (小有名气)
送鲜花一朵
引用回帖:
2楼: Originally posted by ChemiAndy at 2012-06-10 06:03:47
为毛在我这里运行的好好的？
我只是修改了路径设置，并加上了你复制时候漏掉的第一行“&CONTROL“

我查了一下set_hubbard_l.f90，这个里面在过渡金属里面没有In，你看一下，你的set_hubbard_l.f90中，是不是有In啊？我这个没有，怎么加呢？请您帮帮忙吧，现在很纠结啊！
赞一下 回复此楼
什么时候能算完！！！！
6楼2012-06-10 09:16:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ChemiAndy

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★
hanyanli0475: 金币+5, ★★★★★最佳答案 2012-06-10 14:04:22
1. 你修改set_hubbard_l.f90 源文件后是否重新编译了？
2. 你现在用的版本是多少？我用的是最新版5.0.
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

hanyanli0475
7楼2012-06-10 09:51:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ChemiAndy

木虫 (正式写手)
我的set_hubbard_l.f90里面也没有'In'
奇怪，为什么我的能跑呢？
赞一下 回复此楼
8楼2012-06-10 09:54:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ChemiAndy

木虫 (正式写手)

【答案】应助回帖

★
xueht987: 金币+1, 应助指数+1, 谢谢回帖，鼓励交流！ 2012-06-10 10:50:30
问题找到了
5.0版的pwscf中的set_hubbard_l.f90里面是包含In这个元素的。刚才我查的那个文件版本不对。

你下载5.0版重新编译即可。
赞一下 回复此楼
9楼2012-06-10 10:15:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ChemiAndy

木虫 (正式写手)
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
赞一下 回复此楼
10楼2012-06-10 10:16:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 hanyanli0475 的主题更新
151/2返回列表12下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭