The positions of the quenched impurities in the simulation box are
fixed at the beginning of the simulation run following a uniform
random distribution. We also require that to impurities do not
overlap. To achieve this we throw the impurities at random within
the simulation box and then rethrow all overlapping impurities
until the required configuration is obtained. In
Fig. 4.1 we sow the radial distribution function
of a typical configuration of impurities.
Figure 4.1:
Radial distribution function of the impurities for a given
one configuration of disorder