Structural optimization for monolayer In2Se3

Total energy, geometry optimization, DFT+U, spin....

Moderator: bguster

Locked
suman
Posts: 13
Joined: Fri Jan 24, 2020 5:25 pm

Structural optimization for monolayer In2Se3

Post by suman » Fri Jan 24, 2020 6:11 pm

Hello abinit users,
I am trying to relax structure of monolayer In2Se3.
It turns out that output file shows weird result which means x an y coordinates of ions doesn't change but z coordinates has weird result. But relaxation shouldn't produce such weird result.
The input file is as follows
optcell 0
ionmov 3
ntime 20
dilatmx 1.1
ecutsm 0.5
autoparal 1
ecut 30.0
kptopt 1
ngkpt 6 6 1
nshiftk 1
shiftk 0.0 0.0 0.5
toldfe 1e-10
nstep 300
tolmxf 5.0d-6
natom 5
ntypat 2

rprim
6.5887536421 3.8040186889 0.0000000000000000
-6.5887536421 3.8040186889 0.0000000000000000
0.0000000000 0.0000000000 66.140414361901960
vacuum is along z axis of 35 Angstrom.
typat
1 1 2 2 2

xcart
2.1962512140 3.8040186889 3.9740940401
2.1962512140 3.8040186889 9.3201292467
2.1962512140 3.8040186889 1.5448511069
-2.1962512140 3.8040186889 5.1306064284
0.0000000000 0.0000000000 7.9992106855


znucl 49 34

#== Other ==#
nspinor 1
prtden 1
prtwf 0

The output of log file is

xcart 2.1962512140E+00 3.8040186889E+00 -1.8626829133E+03
2.1962512140E+00 3.8040186889E+00 -1.6728811305E+03
2.1962512140E+00 3.8040186889E+00 -1.8863907632E+03
-2.1962512141E+00 3.8040186889E+00 6.6436517266E+02
-3.3333659690E-11 -7.0388485575E-16 4.7855585259E+03
Comparing input and output xcart coordinates shows that there is no displacement of x and y coordinates of ion but drastic change in the z coordinates.
on top of that i got warning which says
scprqt: WARNING -
nstep= 300 was not enough SCF cycles to converge;
maximum energy difference= 4.022E-06 exceeds toldfe= 1.000E-10


Could you please tell me why such weird result pops out?
Could you please tell me any error in input files?
Any suggestion would be highly appreciated.
Thank you

Locked