Skip to Content.
Sympa Menu

forum - Re: [abinit-forum] Problem with symmetry finder

forum@abinit.org

Subject: The ABINIT Users Mailing List ( CLOSED )

List archive

Re: [abinit-forum] Problem with symmetry finder


Chronological Thread 
  • From: Xavier Gonze <xavier.gonze@uclouvain.be>
  • To: forum@abinit.org
  • Subject: Re: [abinit-forum] Problem with symmetry finder
  • Date: Sat, 20 Mar 2010 18:27:26 +0100

Dear Xenophon,

It is a numerical precision issue.

On 19 Mar 2010, at 12:27, Xenophon Krokidis wrote:

Dear All,

I am attaching here the input and output files of an Abinit calculation on the primitive cell of faujasite (space group FD-3, origin 2) with Abinit 5.8.4 and 6.0.2. In both cases the space group found is not correct. Furthermore, both Abinit versions differ: 5.8.4 finds at least a space group with less operations, while 6.0.2 is not able to find a space group at all. In the second attachment there is the input for the same system, but this time using the conventional cell. Here Abinit fails with both versions to generate the full list of atoms.

I had a look at the primitive cell case.
If one uses

rprim 1.0000000000 0.0000000000 0.0000000000
0.5000000000 sqrt(0.75d0) 0.0000000000
0.5000000000 sqrt(0.083333333333d0) sqrt(0.6666666666666d0)

or (even better)

rprim 0 sqrt(.5) sqrt(.5)
sqrt(.5) 0 sqrt(.5)
sqrt(.5) sqrt(.5) 0

The space group is correctly recognized.
Interestingly, with your original rprim,

rprim 1.0000000000 0.0000000000 0.0000000000
0.5000000000 0.8660254038 0.0000000000
0.5000000000 0.2886751346 0.8164965809
(that is accurate up to nine digits)
I would have expected the use of tolsym to allow ABINIT
to recognize the Bravais lattice. But I tested tolsym 1.e-3,
without success.
I will try to see what is problematic.

Best wishes,
Xavier


Best regards,

Xenophon

<FAU_conventional.zip><FAU_primitive.zip><Xenophon_Krokidis.vcf>




Archive powered by MHonArc 2.6.16.

Top of Page