Part 4: Building a stoichiometric cluster
The remaining tasks of Parts 4, 5, and 6 of this tutorial consist of geometry setup tasks and total energy calculations, similar in spirit to Tutorial Basics of Running FHIaims, except for larger and more complex structures.
As mentioned elsewhere, one critical step in defining our scientific question (effectiveness of a particular set of catalytic reaction steps at Fe\(_2\)O\(_3\) nanoparticles) involves defining one or more hypotheses for the structure that such nanoparticles would assume.
We here pursue the hypothesis that such a cluster structure can be found by creating a reasonablysized, maximally compact, stoichiometric cutout from a bulk phase (hematite Fe\(_2\)O\(_3\)). This is a crude hypothesis an certainly not the only one, or the best motivated one. However, it will provide us with cluster models with well defined facets and edges, suitable for further testing of chemical reaction steps later.
We call the definition of these clusters "WignerSeitz" clusters, borrowed from the concept of a compact WignerSeitz cell in a periodic crystal lattice. The DFTPBE relaxed structure and spin state of the primive unit cell of hematite Fe\(_2\)O\(_3\) will be our starting point.
We will also need to ensure a good charge and spin initialization of the electron density for our calculations.
As laid out in detail in Tutorial Basics of Running FHIaims, this step is an inherent step of setting up the input of an electronic structure calculation. Without considering at least the target spin state, a calculation may not converge at all or it may converge to an entirely unwanted, unphysical selfconsistent solution of the spin state that is not the ground state or the experimentally desired state.
The physical principle for our chosen charge and spin initialization is simple, based on chemical expectations:

Assume the formal charge state of an ionic solid:
initial_charge 3.0
for Fe andinitial_charge 2.0
for O. For a nonperiodic system, the current version of FHIaims allows such a charged initialization. This initialization can greatly help improve the initial charge density, moving it closer to the expected final charge distribution of an ionic solid. 
For the Fe(III) ions, assume a highspin state  5 d electrons, spin up:
initial_moment 5
orinitial_moment 5
, with up and down spins initially distributed according to the bulk spin ground state of hematite Fe\(_2\)O\(_3\). 
For the O anions, assume zero initial moments.
Note that these choices only define the initial charge density  the actual charge and spin densities will sort themselves out during the s.c.f. cycle of an electronic structure calculation.
To create these initial geometry, charge, and spin states in geometry.in
, we use a script (see below). The script takes the initial moments and charges specified in the primitive structure and constructs a larger supercell with identically initialized atoms. After the construction of this supercell, all atoms are folded back into the WignerSeitz cell and finally the lattice vectors are removed, to leave behind a nonperiodic cutout of the bulk structure. Et Voila: We have our cluster.

Initialize the atoms with proper spin moments and charges in the
geometry.in
file representing the PBErelaxed primitive unit cell. 
Execute the following script. To reproduce the solution
geometrywscluster.in
we used the following clims command:climswignerseitzcluster translatebasis 1 0 0 2 2 2
and the relaxed structure (
geometry.in.next_step
) from theTutorial2/Part1/uudd_relaxation_light
folder. The clims command shifts the initial basis atoms by 1 A in the Cartesian x direction and creates a 2x2x2 supercell. Obviously, you could vary these parameters to obtain different compact initial cluster representations cut out of the hematite bulk lattice. 
Spend some time playing with the translation of the basis atoms. Ideally, attempyt to find a cluster with only a few singly bonded outermost atoms. Atoms that are completely removed from their stable bonding environment will not be represented well in DFT. They will also typically not exist in lowenergy, chemically reasonable atomic arrangements of a material. But worst of all, such unhappy atoms will often significantly impede the convergence of the s.c.f. cycle of DFT. The wave function associated with such atoms is often not well represented by a single Slater determinant (or by two single determinants in the case of a spinpolarized calculation), creating a significant difficulty for practical density functional approximations.
The cluster we have built here contains only 80 atoms. We here restrict ourselves to a small cluster size since this is already a large system for a tutorial, given that spin and charge density convergence of hematite clusters already pose significant computational challenges of their own.
However, we expect that 80atom clusters will generally be too small to maintain a bulk hematitelike structure that might be found in an experiment. For publishable simulations of clusters with a stable hematitelike structure, larger structures could be created using the script infrastructure above.