[QE-users] Problems with hp.x

dv009200 at fh-muenster.de dv009200 at fh-muenster.de
Mon Jun 15 21:16:53 CEST 2020


Hello everyone,

I'm trying to calculate the hubbard u parameter for Zn in Zinc sulfide
(sphalerite structure) with the help of the hp.x code. The calculations
terminate normally without any errors. The problem is that I get
(presumably) way too high values for U that also won't converge (if I take
the value I got from a one-shot calculation and plug it in the SCF input
and then redo the HP calculation).

For example in the first step I calculate a U = 75.7035 in the second
iteration I get U = 804.2405 and in the third U = 30999.2684.

This seems unreasonable considering that the calculations for the provided
examples in the 'HP' folder work fine and converge fast without such a
massive change to a certain value for U using the above described scheme.

Has someone an idea what is causing this trouble in my system? I already
tried different PPs, functionals, U_projection_type, thresholds and k and
q point grids all without success.

Below is my input for the scf and hp calculation

SCF-input:
 &control
    calculation='scf'
    restart_mode='from_scratch',
    pseudo_dir = '/home/dominik/codes/QE6.5/pseudo/'
    outdir='/home/dominik/codes/QE6.5/tempdir/'
    prefix='zns'
 /
 &SYSTEM
    ibrav = 2
    celldm(1)=10.291937439
    nat   = 2
    ntyp  = 2
    ecutwfc = 60.0
    ecutrho= 720.0
    lda_plus_u = .true.
    lda_plus_u_kind = 0
    U_projection_type = 'atomic'
    Hubbard_U(1) =  1d-8
 /
 &electrons
    mixing_beta=0.7
    conv_thr=1d-15
 /
ATOMIC_SPECIES
 Zn  65.39  Zn.pbe-dn-rrkjus_psl.0.2.2.UPF
 S   32.07  S.pbe-n-rrkjus_psl.0.1.UPF
ATOMIC_POSITIONS {alat}
Zn      0.000000   0.000000   0.000000
S       0.250000   0.250000   0.250000
K_POINTS automatic
12 12 12 0 0 0


HP-input:
 &inputhp
   prefix='zns'
   outdir='/home/dominik/codes/QE6.5/tempdir/'
   nq1 = 2
   nq2 = 2
   nq3 = 2
   conv_thr_chi = 1.0d-10
   iverbosity =2
 /


Best regards


Dominik Voigt

Dominik Voigt
PhD Student University of Applied Sciences Münster
Department of Physical Chemistry



More information about the users mailing list