SQS for R-3c space group

Alloy Theoretic Automated Toolkit (ATAT)
1

Dear AtAt team,

I would like to create a SQS for a structure of space group 167 (R-3c) with 30 atoms in the unit cell.
Therefore I created an input file (see below) and tried to build a 120 atom SQS super cell. Unfortunately I got the error about incompatible supercell size. In the rndstrgrp.out file the symmetrically equivalent La/Sr sites are not grouped together and in the sym.out file I do not recognize the symmetry operations of this space group ? What am I doing wrong ?

Thank you very much for your help !

5.52728 5.52728 13.42116 90 90 120
1 0 0
0 1 0
0 0 1
0.0000000000E+00 -1.2325951644E-32 2.5000000000E-01 La=0.625,Sr=0.375
0.0000000000E+00 0.0000000000E+00 0.0000000000E+00 Fe
4.5772000000E-01 -1.2325951644E-32 2.5000000000E-01 O
5.4228000000E-01 5.4228000000E-01 2.5000000000E-01 O
1.2325951644E-32 0.0000000000E+00 7.5000000000E-01 La=0.625,Sr=0.375
0.0000000000E+00 0.0000000000E+00 5.0000000000E-01 Fe
1.2325951644E-32 5.4228000000E-01 7.5000000000E-01 O
1.2325951644E-32 4.5772000000E-01 2.5000000000E-01 O
5.4228000000E-01 -1.2325951644E-32 7.5000000000E-01 O
4.5772000000E-01 4.5772000000E-01 7.5000000000E-01 O
6.6666666667E-01 3.3333333333E-01 5.8333333333E-01 La=0.625,Sr=0.375
6.6666666667E-01 3.3333333333E-01 3.3333333333E-01 Fe
1.2438666667E-01 3.3333333333E-01 5.8333333333E-01 O
2.0894666667E-01 8.7561333333E-01 5.8333333333E-01 O
6.6666666667E-01 3.3333333333E-01 8.3333333333E-02 La=0.625,Sr=0.375
6.6666666667E-01 3.3333333333E-01 8.3333333333E-01 Fe
6.6666666667E-01 8.7561333333E-01 8.3333333333E-02 O
6.6666666667E-01 7.9105333333E-01 5.8333333333E-01 O
2.0894666667E-01 3.3333333333E-01 8.3333333333E-02 O
1.2438666667E-01 7.9105333333E-01 8.3333333333E-02 O
3.3333333333E-01 6.6666666667E-01 9.1666666667E-01 La=0.625,Sr=0.375
3.3333333333E-01 6.6666666667E-01 6.6666666667E-01 Fe
7.9105333333E-01 6.6666666667E-01 9.1666666667E-01 O
8.7561333333E-01 2.0894666667E-01 9.1666666667E-01 O
3.3333333333E-01 6.6666666667E-01 4.1666666667E-01 La=0.625,Sr=0.375
3.3333333333E-01 6.6666666667E-01 1.6666666667E-01 Fe
3.3333333333E-01 2.0894666667E-01 4.1666666667E-01 O
3.3333333333E-01 1.2438666667E-01 9.1666666667E-01 O
8.7561333333E-01 6.6666666667E-01 4.1666666667E-01 O
7.9105333333E-01 1.2438666667E-01 4.1666666667E-01 O

2

I tried your input file and ran

mcsqs -2=4
mcsqs -n=120

and got no error messages. Can you tell me exactly what you did and exactly what the error message is?

BTW, your cell is unnecessarily large. The primitive cell has only 10 atoms. You can find it with the command:

cellcvrt.exe -s < rndstr.in
3

Dear Axel,

Thank you very much for your reply. In fact when I copy the input file I posted no error occurs. Maybe my original input file contains invisible characters ? You are right the unit cell I used is not primitive (I created the atomic positions using abinit functionality spgroup, which is very useful if one does not want to type all atomic positions by hand :wink: )

I wish you a nice Sunday,
Kind regard,
Cint