Total energy of BiFeO3 for two space groups

MULTIBINIT, aTDEP, ANADDB, optics, cut3d, mrg* postprocessors

Moderators: MMNSchmitt, gonze

Locked
SAIDI
Posts: 4
Joined: Tue Dec 17, 2013 11:37 am

Total energy of BiFeO3 for two space groups

Post by SAIDI » Thu Jan 09, 2014 4:24 pm

Hello,
I use lda+u ( U =4 ev) for calculate the total energy of BiFeO3-R3c
I'm trying to have a total energy convergence as fonction of ecut, but I have a problem with RHOIJ
I have in log file:
(print_ij: WARNING -
The matrix seems to have high value(s) !
( 7 components have a value greater than 10.0).
It can cause instabilities during SCF convergence.
Action: you should check your atomic dataset (psp file)
and look for "high" projector functions...)
I don't know how I can to resolve this problem?

thank you for your help
Best regards.

mverstra
Posts: 655
Joined: Wed Aug 19, 2009 12:01 pm

Re: Total energy of BiFeO3 for two space groups

Post by mverstra » Mon May 26, 2014 10:50 pm

This question is much too vague - we need an input and some more log to see what happens - please read the nettiquette.

Further, this has nothing to do with post-processing tools. Please post in the appropriate sections.

One main possibility is that your PAW datasets are not good or appropriate. Try using them separately for the elemental solids to see if they work well. If one of them gives the same large values of rho_ij, it is probably the culprit, and you need a better dataset. Have you tried the JTH set on the download page of PAW for abinit (link at the top of the page)?

Matthieu
Matthieu Verstraete
University of Liege, Belgium

SAIDI
Posts: 4
Joined: Tue Dec 17, 2013 11:37 am

Re: Total energy of BiFeO3 for two space groups

Post by SAIDI » Tue Jul 08, 2014 3:47 pm

Dear professor Verstraete,
Thank you for your help, Sorry i was in the process of writing my thesis that's why i don't see your response .
Please find attached the input files.


Best regards.
Attachments
BFO-R3c.in
(3.52 KiB) Downloaded 537 times
Bi22.in
(518 Bytes) Downloaded 422 times
Fe-a.in
(945 Bytes) Downloaded 431 times
aO.in
(650 Bytes) Downloaded 400 times

ebousquet
Posts: 469
Joined: Tue Apr 19, 2011 11:13 am
Location: University of Liege, Belgium

Re: Total energy of BiFeO3 for two space groups

Post by ebousquet » Wed Jul 09, 2014 4:29 pm

Dear Saidi,

a few questions:
- why did you use pawovlp=63% precisely?
If you need 63% of overlap between the spheres means the atomic distances in your crystal are probably wrong at some point.

- Did you make your own PAW atomic data or did you get them from somewhere?

A question not related to the crash:
You are using a 40 atoms cell while the unit cell contain 10 atoms, do you really need to run a 40 atoms cell?

All the best,
Eric

SAIDI
Posts: 4
Joined: Tue Dec 17, 2013 11:37 am

Re: Total energy of BiFeO3 for two space groups

Post by SAIDI » Thu Jul 10, 2014 12:21 pm

Dear Eric,
Thank you for your kindness
The input file is already written since January, I forget why i put precisely pawovlp=63%, normally it must be lower to 5 %.
I have the overlap ratio of 10.83 % of the PAW spheres, and in the spherical density matrix (augmentation wave occupancies Rhoij) 7 components have a value greater than 10.0.
I tried to decrease cutoff radius of paw dataset but the problem is still there.

I use 40 atoms cell because the BFO-R3c phases have the G-type AFM
Best regards.

ebousquet
Posts: 469
Joined: Tue Apr 19, 2011 11:13 am
Location: University of Liege, Belgium

Re: Total energy of BiFeO3 for two space groups

Post by ebousquet » Thu Jul 10, 2014 5:01 pm

SAIDI wrote:Dear Eric,
Thank you for your kindness
The input file is already written since January, I forget why i put precisely pawovlp=63%, normally it must be lower to 5 %.
I have the overlap ratio of 10.83 % of the PAW spheres, and in the spherical density matrix (augmentation wave occupancies Rhoij) 7 components have a value greater than 10.0.
I tried to decrease cutoff radius of paw dataset but the problem is still there.


OK you can put pawovlp to 12% this is still fine, don't worry.

SAIDI wrote:I use 40 atoms cell because the BFO-R3c phases have the G-type AFM
Best regards.

Well, the 10 atoms cell can fit the G-type as well but you have to double it along the [111] direction.

SAIDI
Posts: 4
Joined: Tue Dec 17, 2013 11:37 am

Re: Total energy of BiFeO3 for two space groups

Post by SAIDI » Fri Jul 11, 2014 6:02 pm

Dear ebousquet,
Thank you for your help
I'll try to change the number of atoms per unit cell to 10 atoms by unit cell and the pawovlp to 12 % .
Best regards.

cgs12
Posts: 16
Joined: Thu Apr 27, 2017 6:27 pm

Re: Total energy of BiFeO3 for two space groups

Post by cgs12 » Fri Mar 09, 2018 12:38 am

Hi,

Thanks for all the posts and replies. I am also interested in performing calculations for rhombohedral phase (R3c) for PZT (PbZrO3/PbTiO3) system. I just want to ask how we can generate the positions of all the ions in a unit cell for this phase? I understand that it is similar, if not identical to the one posted here for BiFeO3 but is there a logical / more intuitive way of doing so?

Thanks in advance for your kind attention!

ebousquet
Posts: 469
Joined: Tue Apr 19, 2011 11:13 am
Location: University of Liege, Belgium

Re: Total energy of BiFeO3 for two space groups

Post by ebousquet » Fri Mar 09, 2018 10:23 am

Dear cgs12,
PZT can have many different possible unit celles depending on the fraction of Zr and Ti you have and on how you arrange these atoms in the cell. Either you find some papers or thesis, etc that give you the cells and positions or you have to build them yourself. I'm not aware of a program doing so, usually people do their own stuff. A way to go is to build Cartesian supercells (like 2x2x2, 4x4x4, etc), depending on your needs, and either keeping these supercells in your calculations or reducing to their asymmetric unit cell through a program like FINDSYM from the ISOTROPY suite (http://stokes.byu.edu/iso/isotropy.php).
Best wishes,
Eric

cgs12
Posts: 16
Joined: Thu Apr 27, 2017 6:27 pm

Re: Total energy of BiFeO3 for two space groups

Post by cgs12 » Fri Mar 09, 2018 7:28 pm

Dear Eric,

Thank you so much for your response.

I just want to confirm that when you mention unit cell, you are referring to the usual cubic-like unit cell involving 5 atoms. It is also possible to construct a rhombohedral unit cell which in this case would be a primitive unit cell involving 10 atoms?

Best

ebousquet
Posts: 469
Joined: Tue Apr 19, 2011 11:13 am
Location: University of Liege, Belgium

Re: Total energy of BiFeO3 for two space groups

Post by ebousquet » Tue Mar 13, 2018 8:35 am

Dear cgs12,
By unit cell I mean the minimal cell you need to build your crystal from period boundary conditions. If you have P(Zr,Ti)O3 50% Ti and 50% Zr you'll have 10 atoms in the unit cell (two 5-atoms perovskite cell), though several unit cells can be built depending if you order the Zr/Ti along the 100, 110, 111 directions, etc.
All the best,
Eric

Locked