The periodic boundary conditions

Taking into account the finite size of model/simulation box is crucial to computing correctly many of the structural characteristics (e.g. ring statistics) of the system being studied.
The importance of the finite size of model box can be illustrated using a 1 dm3 edged cube of water (1 L) at room temperature. This cube contains approximately 3.3×1025 water molecules, each of them can be considered as a sphere having a diameter of 2.8 Å. Following this scheme surface interactions can affect up to 10 layers of spheres (water molecules) far from the surface of the model cubic box. In this case the number of water molecules exposed to the surface is about 2×1019, which is a small fraction of the total numebr of molecules in the model.
Currently structure models often contain somewhere from 1 thousand to several thousands of molecules/atoms. As a result a very substantial fraction of them will be influenced by the finite size of the simulation/model box. The problems is solved by applying the so-called Periodic Boundary Conditions "PBC" which means surrounding the simulation box with its translational images in the 3 directions of space, as illustrated below. Users of I.S.A.A.C.S. should take special care that their model boxes are inhearently periodic so that when the periodic boundary conditions are applied the structural characteristics computed are not compromized.


Figure 5.1: Schematic representation of the idea of periodic boundary conditions.
Image pbc-seb



Figure [Fig. 5.1] illustrates the principle of the periodic boundary conditions that can be used5.1 in I.S.A.A.C.S.: a particle which goes out from the simulation box by one side is reintroduced in the box by the opposite side (in the 3 dimensions of space). The maximum inter-atomic distance rcut which is taken into account in the calculations is therefore equal to the half of the edge of the simulation box:

rcut = L/2 (5.1)

The surface/finite model size effects would therefore be small, if any. In general, the larger the simulation box and the number of molecules/atoms in it, the smaller the surface/size effects will be.
Sébastien Le Roux 2011-02-26