My Community

Dear VASp' users,

Im trying to do an optimization of As on GaSb, but never finish it I don't know why... Maybe something in my system (POSCAR, KPOINT or INCAR) is wrong... Someone could you help me? please help me because I spent a lot of time trying to do this calculation.

Thank you very much in advance for your help...

INCAR:


PREC = Accurate
ENCUT = 415
IBRION = 2
NSW = 100
ISIF = 2
ALGO = Normal (blocked Davidson)
NELM = 60
NELMIN = 2
EDIFF = 1.0e-05
EDIFFG = -0.02
VOSKOWN = 1
NBLOCK = 1
ISPIN = 1
INIWAV = 1
ISTART = 0
ICHARG = 2
LWAVE = .FALSE.
LCHARG = .FALSE.
ADDGRID = .FALSE.
ISMEAR = 1
SIGMA = 0.2
LREAL = Auto

KPOINT:

Automatic mesh
0
Monkhorst Pack
2 2 1
0. 0. 0.

POSCAR:


Ga36 Sb36 As (P1) ~ CONTCAR (VASP)
1.0
8.62020000 0.00000000 0.00000000
-4.31010000 7.46531219 0.00000000
0.00000000 0.00000000 41.67250000
36 36 1
Direct
0.33333333 0.16666667 0.80611112
0.33333333 0.16666667 0.54485970
0.33333333 0.16666667 0.28887956
0.00000000 0.00000000 0.89259977
0.00000000 0.00000000 0.63175657
0.00000000 0.00000000 0.37332773
0.16666667 0.33333333 0.71894097
0.16666667 0.33333333 0.45777591
0.16666667 0.33333333 0.20443138
0.33333333 0.66666667 0.80611112
0.33333333 0.66666667 0.54485970
0.33333333 0.66666667 0.28887956
0.00000000 0.50000000 0.89259977
0.00000000 0.50000000 0.63175657
0.00000000 0.50000000 0.37332773
0.16666667 0.83333333 0.71894097
0.16666667 0.83333333 0.45777591
0.16666667 0.83333333 0.20443138
0.83333333 0.16666667 0.80611112
0.83333333 0.16666667 0.54485970
0.83333333 0.16666667 0.28887956
0.50000000 0.00000000 0.89259977
0.50000000 0.00000000 0.63175657
0.50000000 0.00000000 0.37332773
0.66666667 0.33333333 0.71894097
0.66666667 0.33333333 0.45777591
0.66666667 0.33333333 0.20443138
0.83333333 0.66666667 0.80611112
0.83333333 0.66666667 0.54485970
0.83333333 0.66666667 0.28887956
0.50000000 0.50000000 0.89259977
0.50000000 0.50000000 0.63175657
0.50000000 0.50000000 0.37332773
0.66666667 0.83333333 0.71894097
0.66666667 0.83333333 0.45777591
0.66666667 0.83333333 0.20443138
0.33333333 0.16666667 0.87149992
0.33333333 0.16666667 0.60952201
0.33333333 0.16666667 0.35221569
0.00000000 0.00000000 0.69649803
0.00000000 0.00000000 0.43666387
0.00000000 0.00000000 0.18331934
0.16666667 0.33333333 0.78393379
0.16666667 0.33333333 0.52255941
0.16666667 0.33333333 0.26776751
0.33333333 0.66666667 0.87149992
0.33333333 0.66666667 0.60952201
0.33333333 0.66666667 0.35221569
0.00000000 0.50000000 0.69649803
0.00000000 0.50000000 0.43666387
0.00000000 0.50000000 0.18331934
0.16666667 0.83333333 0.78393379
0.16666667 0.83333333 0.52255941
0.16666667 0.83333333 0.26776751
0.83333333 0.16666667 0.87149992
0.83333333 0.16666667 0.60952201
0.83333333 0.16666667 0.35221569
0.50000000 0.00000000 0.69649803
0.50000000 0.00000000 0.43666387
0.50000000 0.00000000 0.18331934
0.66666667 0.33333333 0.78393379
0.66666667 0.33333333 0.52255941
0.66666667 0.33333333 0.26776751
0.83333333 0.66666667 0.87149992
0.83333333 0.66666667 0.60952201
0.83333333 0.66666667 0.35221569
0.50000000 0.50000000 0.69649803
0.50000000 0.50000000 0.43666387
0.50000000 0.50000000 0.18331934
0.66666667 0.83333333 0.78393379
0.66666667 0.83333333 0.52255941
0.66666667 0.83333333 0.26776751
0.52000000 0.52000000 0.94200000

These are some results from my OSZICAR:

N E dE d eps ncg rms rms(c)
DAV: 1 -0.263583914908E+03 -0.26358E+03 -0.13693E+01 5064 0.500E+00 0.242E-01
DAV: 2 -0.263599298881E+03 -0.15384E-01 -0.16500E-01 5520 0.642E-01 0.218E-01
DAV: 3 -0.263602228428E+03 -0.29295E-02 -0.55579E-02 5916 0.408E-01 0.304E-01
DAV: 4 -0.263600736938E+03 0.14915E-02 -0.39464E-02 5736 0.229E-01 0.422E-01
DAV: 5 -0.263597314228E+03 0.34227E-02 -0.74471E-03 5244 0.136E-01 0.614E-02
DAV: 6 -0.263597489764E+03 -0.17554E-03 -0.13234E-03 5928 0.138E-01 0.570E-02
DAV: 7 -0.263597287784E+03 0.20198E-03 -0.72931E-04 5040 0.976E-02 0.196E-02
DAV: 8 -0.263597296952E+03 -0.91683E-05 -0.24401E-04 4644 0.389E-02 0.107E-02
DAV: 9 -0.263597297147E+03 -0.19447E-06 -0.50769E-05 4500 0.161E-02
22 F= -.26359730E+03 E0= -.26357106E+03 d E =-.168973E-01

Hi,
I would start imposing less tight conditions.
Prec=Accurate
EDIFF = 1.0e-05
EDIFFG = -0.02
And also a KPOINT 2x2x1
can be quite time consuming.

Start with something of "milder"...
Gamma point,
Prec=Medium,
EDIFFG=-0.08,
EDIFF=-1.0e-03.
Then you get a rough converged structure you can restart (ISTART=1) just refining the initial values up to those you need.
Mind that you are studying an interstitial defect without relaxing the lattice parameter. Thus there could be some residual stress that could not facilitate the achievement of your (quite tight) convergence parameters.

I hope it helps,
Giacomo