My Community

Recently, I try to optimize the structure of Al2O3 by E-V.
My VASP version is 4.6. Related files are as follows:

INCAR file(Relaxation):
System=a_Al2O3
ENCUT=400
ISTART=0
ISMEAR=0
SIGMA=0.05
PREC=Accuracy
IALGO=38
#ICHARG=11
GGA=91;VOSKOWN=1;EDIFF=0.1E-04; EDIFFG=-0.01
IBRION=2
NSW=50
ISIF=4
NPAR=40
LPLANE=.TRUE.
NSIM=4
LREAL=.FALSE.
LWAVE=.FALSE.
LCHG=.FALSE.

INCAR file(Energy):
ISMEAR=-5
ISIF=2

POSCAR file(Relaxation):
a_Al2O3
1.0000
4.122540729635063 -2.380150000000001 0.000000000000000
0.000000000000000 4.760300000000000 0.000000000000000
0.000000000000000 0.000000000000000 12.993300000000000
18 12
Direct
0.3062399923801409 0.0000000000000000 0.2500000000000000
0.9729066590468075 0.3333333333333334 0.5833333333333334
0.6395733257134743 0.6666666666666666 0.9166666666666667
0.0000000000000000 0.3062399923801410 0.2500000000000000
0.6666666666666666 0.6395733257134743 0.5833333333333334
0.3333333333333333 0.9729066590468076 0.9166666666666667
0.6937600076198590 0.6937600076198592 0.2500000000000001
0.3604266742865256 0.0270933409531925 0.5833333333333334
0.0270933409531924 0.3604266742865257 0.9166666666666667
0.6937600076198590 0.0000000000000000 0.7500000000000000
0.3604266742865257 0.3333333333333334 0.0833333333333334
0.0270933409531924 0.6666666666666669 0.4166666666666667
0.0000000000000000 0.6937600076198590 0.7500000000000000
0.6666666666666666 0.0270933409531924 0.0833333333333334
0.3333333333333333 0.3604266742865257 0.4166666666666667
0.3062399923801410 0.3062399923801410 0.7500000000000000
0.9729066590468075 0.6395733257134743 0.0833333333333334
0.6395733257134743 0.9729066590468078 0.4166666666666667
0.0000000000000000 0.0000000000000000 0.3521600067615511
0.6666666666666666 0.3333333333333335 0.6854933400948845
0.3333333333333333 0.6666666666666666 0.0188266734282178
0.0000000000000000 0.0000000000000000 0.1478399932384489
0.6666666666666666 0.3333333333333335 0.4811733265717821
0.3333333333333333 0.6666666666666666 0.8145066599051154
0.0000000000000000 0.0000000000000000 0.6478399932384489
0.6666666666666666 0.3333333333333335 0.9811733265717820
0.3333333333333333 0.6666666666666666 0.3145066599051156
0.0000000000000000 0.0000000000000000 0.8521600067615511
0.6666666666666666 0.3333333333333335 0.1854933400948845
0.3333333333333333 0.6666666666666666 0.5188266734282176

I try several scale coefficient from 0.98 to 1.02. Relation between scale and energy is as follows:
scale energy
0.98 225.31995ev
0.99 226.25659ev
1.00 226.78708ev
1.01 226.94750ev
1.02 226.76866ev
Accordingly, the scale should be 1.0097. But it seems not correct, because the external pressure is not close to ZERO.
scale external pressure
0.98 196.98kB
0.99 94.76kB
1.00 3.35kB
1.01 -68.07kB
1.02 -130.72kB

I try ISIF =3 to relax.

Direct lattice vectors
4.126580148 -2.382482159 0.000000000
0.000000000 4.764964319 0.000000000
0.000000000 0.000000000 12.996271259

It seems reasonable. So i infer there must be some problems in my privious way. Who can give me some sugguestions?

3Q x)

Hi,

We're sorry that we didn’t answer your question. This does not live up to the quality of support that we aim to provide. The team has since expanded. If we can still help with your problem, please ask again in a new post, linking to this one, and we will answer as quickly as possible.

Best wishes,

VASP