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.
berry phase-BUG
Moderator: bguster
Re: berry phase-BUG
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?
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?
Josef W. Zwanziger
Professor, Department of Chemistry
Canada Research Chair in NMR Studies of Materials
Dalhousie University
Halifax, NS B3H 4J3 Canada
jzwanzig@gmail.com
Professor, Department of Chemistry
Canada Research Chair in NMR Studies of Materials
Dalhousie University
Halifax, NS B3H 4J3 Canada
jzwanzig@gmail.com
Re: berry phase-BUG
Thankyou professor jzwanzig,i think my problem is this i have problem with primitive and conventinal cell.could you please explain them?
best regards.
best regards.
Re: berry phase-BUG
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.
Josef W. Zwanziger
Professor, Department of Chemistry
Canada Research Chair in NMR Studies of Materials
Dalhousie University
Halifax, NS B3H 4J3 Canada
jzwanzig@gmail.com
Professor, Department of Chemistry
Canada Research Chair in NMR Studies of Materials
Dalhousie University
Halifax, NS B3H 4J3 Canada
jzwanzig@gmail.com