Page 1 of 1

berry phase-BUG

Posted: Tue Nov 16, 2010 12:53 pm
by mina
Hello,
i'm calculating polarization of bifeo3 by abinit.when i run my job i have a BUG for my berry phase.i have sent my input file and my bug.
please guide me what part of my input is rong?what i have to do with this BUG?
thankyou.Best regards.
input file:
berryopt -1
rfdir 1 1 1
kptopt 1
ngkpt 2 2 2
nshiftk 1
shiftk 0.5 0.5 0.5
ntypat 3
znucl 83 26 8
nstep 100
toldfe 1.0d-10
diemac 9.0
natom 10
natrd 3
typat 6*1 6*2 18*3
acell 5.58132 5.58132 13.87698
angdeg 90 90 120
chkprim 0
spgroup 161
brvltt -1
xred 0.000000000 0.000000000 0.000000000
0.000000000 0.000000000 0.221200004
0.442999989 0.012000000 0.954299986
ecut 30
ixc 1
iscf 3

berryphase_new : BUG -
For k-point # 3,
the determinant of the overlap matrix is found to be 0.
Action : contact ABINIT group.

Re: berry phase-BUG

Posted: Tue Nov 16, 2010 7:52 pm
by jzwanzig
I'm not sure your kpt grid is consistent with your lattice, I think you should be using shiftk 0 0 1/2 in this case. Use prtkpt 1 to find a list of useable kpt grids.

Also, your natom is not consistent with typat: you have natom 10, but define 30 atoms in your typat statement. I take it you are trying to construct the primitive unit cell from the conventional cell (because you have brvltt -1) right?

Re: berry phase-BUG

Posted: Sat Nov 20, 2010 7:13 am
by mina
Thankyou professor jzwanzig,i think my problem is this i have problem with primitive and conventinal cell.could you please explain them?
best regards.

Re: berry phase-BUG

Posted: Mon Nov 22, 2010 6:07 pm
by jzwanzig
Just review any basic book on crystallography. The idea is this: for some space groups, the unit cell that displays all symmetries (called the conventional unit cell) is bigger than the minimal unit cell; for other space groups, the two cells are the same. If the space group name starts with a P, then the primitive cell coincides with the conventional cell. But for example group 225, Fm-3m, which is face-centered cubic, has a conventional unit cell that is a cube but a smaller primitive cell. You will find simple examples like this in solid state physics texts like Ashcroft and Mermin, for more complicated examples (body centered tetragonal, face centered monoclinic etc) you have to go to a real crystallography book.