Maintaining a designated space group for a relaxation

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

Moderator: bguster

Locked
brehmj
Posts: 35
Joined: Thu Jan 20, 2011 3:18 pm

Maintaining a designated space group for a relaxation

Post by brehmj » Wed Oct 26, 2011 5:21 pm

Hi everyone:
maybe I just am not understanding the abinit input variables involving space groups.

Here is my question: How do I perform a relaxation that maintains a particular space group throughout?
say I input starting coords that are for space group 63. How do I keep abinit only searching within space group 63 for a minimum in energy, and not go to solutions outside of this space group?

Thanks,
John

jolafc
Posts: 34
Joined: Wed Aug 26, 2009 5:37 am

Re: Maintaining a designated space group for a relaxation

Post by jolafc » Wed Oct 26, 2011 5:55 pm

Hi,

By symmetry, the forces are supposed to keep your structure in the same space group throughout relaxation. However, the numerical errors can add up to a level where Abinit breaks the symmetry of your structure during relaxation. To help prevent this, it is useful to input coordinates using fractions as much as possible. For instance, with an fcc structure (say, silicon), one could use :

Code: Select all

rprim
1/2 1/2 0
1/2 0 1/2
0 0 1/2
xred
0 0 0
1/4 1/4 1/4

Instead of giving the coordinates using xcart with a lot of digits. Is it what you're currently doing? Posting your input file would help us diagnose your problem...

Jonathan

brehmj
Posts: 35
Joined: Thu Jan 20, 2011 3:18 pm

Re: Maintaining a designated space group for a relaxation

Post by brehmj » Wed Oct 26, 2011 6:43 pm

I am not looking for a ground state though... for example: I would like to see how a NaCl BCC energy compares to a NaCL FCC ground state. Now, this is an easy state as none of the coords are on decimal sites.... but a lot of xray diffraction results show atoms to be on decimal sites. So, for another example, suppose I pull off of ICSD coordinates for a SrMnO3. And the coords are on decimal sites. And now I'd like to see the energy of a BaMnO3 (for this email, let's assume it is not in the ICSD database) in the same space group as SrMnO3. I use as my starting coords the xred for the ICSD SrMnO3. But, I don't want ABINIT to leave that space group.
How can I ensure that it won't?

Locked