Crystal building in Atomes
The crystal building process in Atomes will be briefly presented in the next pages. 2 example cases will be used to illustrate the methodology implemented in the program:
The Fd
m space group, N 227, with the first setting: origin 1.The R
c space group, N 167, with the second setting: rhombohedral parameters.
The terminology used in the Atomes program, and therefore used in this manual, is directly taken from the International tables for Crystallography Vol. A. [1]
In the following every single step depends on the considered space group, and corresponds to the way the process is implemented in Atomes:
Determine the origin of the coordinate system.
Calculate the inverse transformation matrix.
Select the space group extra site positions, if any.
Adjust the different symmetrical positions of the first Wyckoff position using [cb:2].
Calculate the crystalline positions using [cb:4] and [cb:3].
Select the unique position(s) in [cb:5] to build the crystal.
Calculate the object position in Cartesian coordinates using [cb:6].
Insert object(s) at the appropriate coordinates based on occupancy.
Origin of the coordinate system
The origin of the coordinate system, if any particular, determine how the object fractional coordinates are modified to match the space group criteria before calculating the Cartesian coordinates. The first step is to build the 4
Ex: (a, b, c) and can be converted in:
(1) Fd | (2) R | |
![]() | ![]() | |
(a- | ( | |
The corresponding P matrix are: | ||
Inverse transformation matrix
To transform the coordinates to the appropriate system we use the corresponding inverse matrix P
(1) Fd | (2) R | |
The P
Space group extra site position(s)
The space group extra site position(s) describe translations to be performed on the initial set of Wyckoff position to get the entire list of coordinates.
The list of extra position(s),
(1) Fd | (2) R | |
Symmetrical positions using the first Wyckoff position
The first step is to convert the list of coordinates from the first Wyckoff position on the form of matrix W:
(1) Fd | (2) R | |
(x, y, z) | (x, y, z) | |
The W' matrix, suitable for the calculation, is calculated using the P
(1) Fd | (2) R | |
The results, converted back to standard coordinates, are presented in the Wyckoff positions table of the space group info window:
(1) Fd | (2) R | |
Object position(s)
The first step in getting the final fractional coordinates
Removing duplicate position(s)
Since one of the goal in Atomes is to build crystal structure of super-structure to work with, all crystalline positions cannot be conserved and duplicate positions (same positions using the periodic boundary conditions) must be removed.
When final positions
In Atomes a duplicate position
and | |
and |
Calculate the object final position(s) in Cartesian coordinates
After the construction/selection process Cartesian coordinates
a, b, c and
Occupancy
The only thing left to do after that it to select which set of coordinates are to be occupied among the list prepared in the previous steps. Depending on the occupancy defined at input Atomes does the following:
If the occupancy is equal to 1.0: all coordinates are to be occupied and no particular selection is performed.
If the occupancy is < 1.0, then the following options are possible to handle occupancy:
"
Random for the initial cell only
": sites are filled randomly in the initial cell only, then the initial cell is simply replicated."
Random cell by cell
": sites are filled randomly for each cell, cell by cell separately."
Completely random
": sites are filled randomly for the entire network, the final crystal is considered as whole."
Successively
": sites are filled successively, all object(s) A are inserted (for the first positions), then all object(s) B are inserted (for the next positions) ..."
Alternatively
": sites are filled alternatively: object A is inserted on the first position, object B is inserted on the second position, object A on the third position, object B on the fourth position ... and so on.
For [i], [ii], [iv] and [v] the number of object(s) by cell is constant, but it can vary for [iii].
For [iv] and [v] the order of the positions to be filled is the order they are found in during the previous steps.
Final check: distance(s) between inserted objects
Finally the last calculation performed by Atomes when creating the crystal is distance check between the object(s) be inserted. Note that this check performed, if and only if, no overlapping is allowed.
If any distance
If
Å then a warning message will pop-up and ask for the user confirmation to continue, or not, to build the crystal.If the user decides to build the crystal anyway and if
then for fail safe reasons object is removed from the list to build.
The "SGL
" "Space group info XML file" format
The data files that contain the information regarding each of the 230 space groups are located in the "bin/library/space_groups
" directory.
In this folder to each space group corresponds a file with a name following the construction:
"space group number
"-"space group name
"."sgl
"
For the Fd227-Fd-3m.sgl
".
These files must have the ".sgl
" extension, follow the "XML
" encoding rules, and the following structure:
<?xml version="1.0" encoding="UTF-8"?>
<!– Space group info XML file –>
<sg-xml>
<space-group>R-3c</space-group>
<sg-num>167</sg-num>
<hm-symbol>R -3 2/c</hm-symbol>
<bravais>Trigonal</bravais>
<settings num="2">
<set name="R_-3_c_:h" x="a" y="b" z="c">
<points num="3">
<pt x="0" y="0" z="0"/>
<pt x="2/3" y="1/3" z="1/3"/>
<pt x="1/3" y="2/3" z="2/3"/>
</points>
</set>
<set name="R_-3_c_:r" x="2/3a+1/3b+1/3c" y="-1/3a+1/3b+1/3c" z="-1/3a-2/3b+1/3c">
<points num="3">
<pt x="0" y="0" z="0"/>
<pt x="0" y="0" z="0"/>
<pt x="0" y="0" z="0"/>
</points>
</set>
</settings>
<wyckoff num="6">
<wyck id="1" mul="12" let="f" site="1">
<pos x="x" y="y" z="z"/>
<pos x="-y" y="x-y" z="z"/>
<pos x="-x+y" y="-x" z="z"/>
<pos x="y" y="x" z="-z+1/2"/>
<pos x="x-y" y="-y" z="-z+1/2"/>
<pos x="-x" y="-x+y" z="-z+1/2"/>
<pos x="-x" y="-y" z="-z"/>
<pos x="y" y="-x+y" z="-z"/>
<pos x="x-y" y="x" z="-z"/>
<pos x="-y" y="-x" z="z+1/2"/>
<pos x="-x+y" y="y" z="z+1/2"/>
<pos x="x" y="x-y" z="z+1/2"/>
</wyck>
<wyck id="2" mul="6" let="e" site=".2">
<pos x="x" y="0" z="1/4"/>
<pos x="0" y="x" z="1/4"/>
<pos x="-x" y="-x" z="1/4"/>
<pos x="-x" y="0" z="3/4"/>
<pos x="0" y="-x" z="3/4"/>
<pos x="x" y="x" z="3/4"/>
</wyck>
<wyck id="3" mul="6" let="d" site="-1">
<pos x="1/2" y="0" z="0"/>
<pos x="0" y="1/2" z="0"/>
<pos x="1/2" y="1/2" z="0"/>
<pos x="0" y="1/2" z="1/2"/>
<pos x="1/2" y="0" z="1/2"/>
<pos x="1/2" y="1/2" z="1/2"/>
</wyck>
<wyck id="4" mul="4" let="c" site="3.">
<pos x="0" y="0" z="z"/>
<pos x="0" y="0" z="-z+1/2"/>
<pos x="0" y="0" z="-z"/>
<pos x="0" y="0" z="z+1/2"/>
</wyck>
<wyck id="5" mul="2" let="b" site="-3.">
<pos x="0" y="0" z="0"/>
<pos x="0" y="0" z="1/2"/>
</wyck>
<wyck id="6" mul="2" let="a" site="32">
<pos x="0" y="0" z="1/4"/>
<pos x="0" y="0" z="3/4"/>
</wyck>
</wyckoff>
</sg-xml>
- "Space-group symmetry," International Tables for Crystallography, vol. A, 2016.