[Wien] Question of bad E(TOP) with no error

夏宇阳 harriron at sjtu.edu.cn
Wed Apr 10 04:52:52 CEST 2024


Dear all,
I am doing a project about KFe2As2 surface.
After relaxing the surface slab in prec 1, i did a prec 2 calculation with cc 0.0001.
There is no error in the process, but when i check the scf1up and scf1dn, i found that E(TOP) is not good.

ATOMIC SPHERE DEPENDENT PARAMETERS FOR ATOM  As1       
:e__0004: OVERALL ENERGY PARAMETER IS    0.0158
          OVERALL BASIS SET ON ATOM IS LAPW
:E2_0004: E( 2)=    0.0158
             APW+lo
:E2_0004: E( 2)=   -2.3962   E(BOTTOM)=   -2.455   E(TOP)=   -2.337  0  1   120
             LOCAL ORBITAL
:E0_0004: E( 0)=    0.4158
             APW+lo
:E0_0004: E( 0)=   -0.8200   E(BOTTOM)=   -1.495   E(TOP)= -520.000  3 -1   175
             LOCAL ORBITAL
:E1_0004: E( 1)=    0.0158
             APW+lo
:E1_0004: E( 1)=    0.0158
             LOCAL ORBITAL(SECDER)

          ATOMIC SPHERE DEPENDENT PARAMETERS FOR ATOM  As2       
:e__0005: OVERALL ENERGY PARAMETER IS    0.0158
          OVERALL BASIS SET ON ATOM IS LAPW
:E2_0005: E( 2)=    0.0158
             APW+lo
:E2_0005: E( 2)=   -2.4026   E(BOTTOM)=   -2.462   E(TOP)=   -2.344  0  1   124
             LOCAL ORBITAL
:E0_0005: E( 0)=    0.4158
             APW+lo
:E0_0005: E( 0)=   -0.8200   E(BOTTOM)=   -1.506   E(TOP)= -520.000  3 -1   175
             LOCAL ORBITAL
:E1_0005: E( 1)=    0.0158
             APW+lo
:E1_0005: E( 1)=    0.0158
             LOCAL ORBITAL(SECDER)

Is there any influence to the results? if is,what should i do?
The structure is 1X1X1 supercell with 25 vacuum in z(for less time),rmt reduced 3%.

init_lapw -perc 1 -ecut -7.0 -sp 
runsp_lapw -i 1000 -fc 3 -p
save_lapw unrelaxed
runsp_lapw -i 1000 -min -fc 1 -p
save_lapw relaxed
init_lapw -prec 2 -ecut -7.0 -sp 
runsp_lapw -i 1000 -cc 0.0001 -p

in outputst,the As is

                    As                                  RHFS

 FOR SPIN  1   R-MT=   2.12000   R*V=  -3.20524  V=  -1.51191
 FOR SPIN  2   R-MT=   2.12000   R*V=  -2.68025  V=  -1.26427

 TOTAL CHARGE FOR SPIN            1 :                  18.000000000108162     
 TOTAL CHARGE FOR SPIN            1  INSIDE SPHERE:    15.842689296044156     
 TOTAL CHARGE FOR SPIN            2 :                  15.000000000108740     
 TOTAL CHARGE FOR SPIN            2  INSIDE SPHERE:    14.639384914097263     
 TOTAL CHARGE in sigma FOR SPIN            1 :                 0.0000000000000000     
 TOTAL CHARGE in sigma FOR SPIN            1 INSIDE SPHERE:    0.0000000000000000     
 TOTAL CHARGE in sigma FOR SPIN            2 :                 0.0000000000000000     
 TOTAL CHARGE in sigma FOR SPIN            2 INSIDE SPHERE:    0.0000000000000000     

          E-up(Ry)      E-dn(Ry)   Occupancy   q/sphere  core-state
  1S    -861.070611   -861.062271  1.00  1.00    1.0000  T
  2S    -109.198830   -109.195361  1.00  1.00    1.0000  T
  2P*    -97.557529    -97.550205  1.00  1.00    1.0000  T
  2P     -94.870103    -94.863015  2.00  2.00    1.0000  T
  3S     -13.958022    -13.947526  1.00  1.00    1.0000  T
  3P*    -10.047612    -10.035593  1.00  1.00    0.9999  T
  3P      -9.678540     -9.666373  2.00  2.00    0.9999  T
  3D*     -2.998617     -2.979436  2.00  2.00    0.9965  F
  3D      -2.945248     -2.925727  3.00  3.00    0.9963  F
  4S      -1.117588     -0.948715  1.00  1.00    0.6527  F
  4P*     -0.448722     -0.278655  1.00  0.00    0.4165  F
  4P      -0.428475     -0.258969  2.00  0.00    0.3984  F

 TOTAL CORE-CHARGE:                   18.000000
 TOTAL CORE-CHARGE INSIDE SPHERE:     17.999236
 TOTAL CORE-CHARGE OUTSIDE SPHERE:     0.000764

 TOTAL ENERGY (RYD):            -4521.980893
     SUM OF EI:-2.6347206E+03     NUC:-1.0916132E+04     COUL:-7.0407530E+03
     V-XC SPIN 1:-1.1047898E+02     E-XC SPIN 1:-8.4909182E+01
     V-XC SPIN 2:-1.0750097E+02     E-XC SPIN 2:-8.2641601E+01

In addition, in the bulk calulation, the As still has a bad E(TOP).

Looking forward to your reply.

With regards!

Yuyang Xia


More information about the Wien mailing list