I have a supercell that has 12 layers, each with four atoms. I want to generate an SQS with local ordering at layer numbers 6, 7 and 8. Rest of the layers should be disordered.
I followed the procedure "How to use SQS for disordered materials."
but I get the following error:
Impossible to match point correlations due to incompatible supercell size.
The input code (I am considering only the pairs):
mcsqs -2=4
mcsqs -rc
The rndstr.in
file(A,B,C,D can be any metallic element):
6.92820300 0.00000000 0.00000000
2.30940100 6.53197300 0.00000000
0.00000000 0.00000000 33.94112500
1 0 0
0 1 0
0 0 1
1 0.5 1 A=.25,B=.25,C=.25,D=.25
1 1 1 A=.25,B=.25,C=.25,D=.25
0.5 1 1 A=.25,B=.25,C=.25,D=.25
0.5 0.5 1 A=.25,B=.25,C=.25,D=.25
0.25 0.25 0.916667 A=.25,B=.25,C=.25,D=.25
0.75 0.25 0.916667 A=.25,B=.25,C=.25,D=.25
0.75 0.75 0.916667 A=.25,B=.25,C=.25,D=.25
0.25 0.75 0.916667 A=.25,B=.25,C=.25,D=.25
1 1 0.833333 A=.25,B=.25,C=.25,D=.25
1 0.5 0.833333 A=.25,B=.25,C=.25,D=.25
0.5 1 0.833333 A=.25,B=.25,C=.25,D=.25
0.5 0.5 0.833333 A=.25,B=.25,C=.25,D=.25
0.75 0.25 0.75 A=.25,B=.25,C=.25,D=.25
0.25 0.75 0.75 A=.25,B=.25,C=.25,D=.25
0.25 0.25 0.75 A=.25,B=.25,C=.25,D=.25
0.75 0.75 0.75 A=.25,B=.25,C=.25,D=.25
0.5 0.5 0.666667 A=.25,B=.25,C=.25,D=.25
0.5 1 0.666667 A=.25,B=.25,C=.25,D=.25
1 1 0.666667 A=.25,B=.25,C=.25,D=.25
1 0.5 0.666667 A=.25,B=.25,C=.25,D=.25
0.25 0.25 0.583333 B
0.75 0.25 0.583333 C
0.75 0.75 0.583333 C
0.25 0.75 0.583333 D
1 1 0.5 B
1 0.5 0.5 C
0.5 1 0.5 D
0.5 0.5 0.5 A
0.25 0.75 0.416667 B
0.25 0.25 0.416667 D
0.75 0.75 0.416667 A
0.75 0.25 0.416667 A
0.5 0.5 0.333333 A=.25,B=.25,C=.25,D=.25
0.5 1 0.333333 A=.25,B=.25,C=.25,D=.25
1 1 0.333333 A=.25,B=.25,C=.25,D=.25
1 0.5 0.333333 A=.25,B=.25,C=.25,D=.25
0.25 0.25 0.25 A=.25,B=.25,C=.25,D=.25
0.75 0.75 0.25 A=.25,B=.25,C=.25,D=.25
0.75 0.25 0.25 A=.25,B=.25,C=.25,D=.25
0.25 0.75 0.25 A=.25,B=.25,C=.25,D=.25
0.5 0.5 0.166667 A=.25,B=.25,C=.25,D=.25
0.5 1 0.166667 A=.25,B=.25,C=.25,D=.25
1 0.5 0.166667 A=.25,B=.25,C=.25,D=.25
1 1 0.166667 A=.25,B=.25,C=.25,D=.25
0.25 0.75 0.083333 A=.25,B=.25,C=.25,D=.25
0.75 0.75 0.083333 A=.25,B=.25,C=.25,D=.25
0.75 0.25 0.083333 A=.25,B=.25,C=.25,D=.25
0.25 0.25 0.083333 A=.25,B=.25,C=.25,D=.25
Please suggest if I am doing something wrong, or is there anything else further which I need to specify to the ATAT MCSQS code?