SIMPLECUBIC
This is part of the symfunc module
It is only available if you configure PLUMED with ./configure –enable-modules=symfunc . Furthermore, this feature is still being developed so take care when using it and report any problems on the mailing list.

Calculate whether or not the coordination spheres of atoms are arranged as they would be in a simple cubic structure.

We can measure how similar the environment around atom \(i\) is to a simple cubic structure is by evaluating the following quantity:

\[ s_i = \frac{ \sum_{i \ne j} \sigma(r_{ij}) \left[ \frac{ x_{ij}^4 + y_{ij}^4 + z_{ij}^4 }{r_{ij}^4} \right] }{ \sum_{i \ne j} \sigma(r_{ij}) } \]

In this expression \(x_{ij}\), \(y_{ij}\) and \(z_{ij}\) are the \(x\), \(y\) and \(z\) components of the vector connecting atom \(i\) to atom \(j\) and \(r_{ij}\) is the magnitude of this vector. \(\sigma(r_{ij})\) is a switchingfunction that acts on the distance between atom \(i\) and atom \(j\) and its inclusion in the numerator and the denominator of the above expression as well as the fact that we are summing over all of the other atoms in the system ensures that we are calculating an average of the function of \(x_{ij}\), \(y_{ij}\) and \(z_{ij}\) for the atoms in the first coordination sphere around atom \(i\). This quantity is once again a multicolvar so you can compute it for multiple atoms using a single PLUMED action and then compute the average value for the atoms in your system, the number of atoms that have an \(s_i\) value that is more that some target and so on. Notice also that you can rotate the reference frame if you are using a non-standard unit cell.

Examples

The following input tells plumed to calculate the simple cubic parameter for the atoms 1-100 with themselves. The mean value is then calculated.

Click on the labels of the actions for more information on what each action computes
tested on master