[QE-users] Convergence issue of a 2x2x2 supercell with Ni (111) surface
Suraj P
surajp at iitkgp.ac.in
Sat Jul 6 12:04:40 CEST 2024
Convergence issue of a 2x2x2 supercell with Ni (111) surface
Dear all,
My objective of using QE is to calculate the adsorption energy of diphenyl ether over Ni (111) surface.
Initially I have taken CIF file of FCC Ni and created a 2x2x2 supercell (7 layer) of Ni in VESTA.
Next I have transformed the (111) surface of this supercell, so as to expose the (111) surface to top of this supercell.
Next I have added 10A vacuum in z direction.(After this the supercell was periodic in x and y direction but not in z direction).
Next I have done a relaxation calculation where in Im facing the convergence problem.
I have tried various pseudopotentials (PAW and ultrasoft), also I have tried with Cg and Davidson diagonilization.
I have tried with ibrav =8 (Orthorhombic) and 14 (triclinic).
Ecutwfc and Ecutrho were chosen from the pseudopotentiafile. However the convergence is not achieved in any case.
Kindly have a look into my input file and help me in troubleshooting this issue.
______________________________________________________________________________________________________________________________________________________
&CONTROL
calculation = "relax"
prefix= 'Ni222'
tstress = .true.
tprnfor = .true.
verbosity = 'high'
etot_conv_thr = 1e-6
forc_conv_thr = 1e-5
pseudo_dir='./pseudo'
outdir = './outdir'
/
&SYSTEM
a = 8.51233e+00
b = 4.91460e+00
c = 2.20383e+01
ibrav = 8
nat = 56
ntyp = 1
ecutwfc =45
ecutrho =250
noncolin =.true.
lspinorb =.true.
!nspin =2
occupations ='smearing'
smearing ='gaussian'
degauss =0.002
starting_magnetization = 0.5
/
&ELECTRONS
conv_thr = 1.00000e-06
mixing_beta = 7.00000e-01
startingpot = "atomic"
startingwfc = "atomic+random"
/
&IONS
ion_positions='from_input'
ion_dynamics = "bfgs"
/
&CELL
/
ATOMIC_SPECIES
Ni 58.69340 Ni.pbe-spn-kjpaw_psl.1.0.0.UPF
ATOMIC_POSITIONS {angstrom}
Ni 0.000000 0.000000 0.000000
Ni 0.000000 0.000000 12.038267
Ni 5.674892 0.000000 8.025519
Ni 2.837442 0.000000 4.012748
Ni 2.837442 0.000000 10.031882
Ni 0.000000 0.000000 6.019134
Ni 5.674892 0.000000 2.006385
Ni 7.093609 2.457299 10.031882
Ni 4.256167 2.457299 6.019134
Ni 1.418725 2.457299 2.006385
Ni 1.418725 2.457299 8.025519
Ni 7.093609 2.457299 4.012748
Ni 4.256167 2.457299 0.000000
Ni 4.256167 2.457299 12.038267
Ni 0.709358 1.228650 10.031882
Ni 6.384250 1.228650 6.019134
Ni 3.546809 1.228650 2.006385
Ni 3.546809 1.228650 8.025519
Ni 0.709358 1.228650 4.012748
Ni 6.384250 1.228650 0.000000
Ni 6.384250 1.228650 12.038267
Ni 7.802976 3.685949 8.025519
Ni 4.965525 3.685949 4.012748
Ni 2.128084 3.685949 0.000000
Ni 2.128084 3.685949 12.038267
Ni 2.128084 3.685949 6.019134
Ni 7.802976 3.685949 2.006385
Ni 4.965525 3.685949 10.031882
Ni 0.709358 3.685949 10.031882
Ni 6.384250 3.685949 6.019134
Ni 3.546809 3.685949 2.006385
Ni 3.546809 3.685949 8.025519
Ni 0.709358 3.685949 4.012748
Ni 6.384250 3.685949 0.000000
Ni 6.384250 3.685949 12.038267
Ni 7.802976 1.228650 8.025519
Ni 4.965525 1.228650 4.012748
Ni 2.128084 1.228650 0.000000
Ni 2.128084 1.228650 12.038267
Ni 2.128084 1.228650 6.019134
Ni 7.802976 1.228650 2.006385
Ni 4.965525 1.228650 10.031882
Ni 7.093609 0.000000 10.031882
Ni 4.256167 0.000000 6.019134
Ni 1.418725 0.000000 2.006385
Ni 1.418725 0.000000 8.025519
Ni 7.093609 0.000000 4.012748
Ni 4.256167 0.000000 0.000000
Ni 4.256167 0.000000 12.038267
Ni 5.674892 2.457299 8.025519
Ni 2.837442 2.457299 4.012748
Ni 0.000000 2.457299 0.000000
Ni 0.000000 2.457299 12.038267
Ni 0.000000 2.457299 6.019134
Ni 5.674892 2.457299 2.006385
Ni 2.837442 2.457299 10.031882
K_POINTS {automatic}
1 1 1 0 0 0
______________________________________________________________________________________________________________________________________________________
Thanking you
Suraj
Research student, IIT Kharagpur
More information about the users
mailing list