This is part of the volumes module |
This quantity can be used to calculate functions of the distribution of collective variables for the atoms that lie in a particular, user-specified part of of the cell.
Each of the base quantities calculated by a multicolvar can can be assigned to a particular point in three dimensional space. For example, if we have the coordination numbers for all the atoms in the system each coordination number can be assumed to lie on the position of the central atom. Because each base quantity can be assigned to a particular point in space we can calculate functions of the distribution of base quantities in a particular part of the box by using:
\[ \overline{s}_{\tau} = \frac{ \sum_i f(s_i) w(x_i,y_i,z_i) }{ \sum_i w(x_i,y_i,z_i) } \]
where the sum is over the collective variables, \(s_i\), each of which can be thought to be at \( (x_i,y_i,z_i)\). The function \( w(x_i,y_i,z_i) \) measures whether or not the system is in the subregion of interest. It is equal to:
\[ w(x_i,y_i,z_i) = \int_{xl}^{xu} \int_{yl}^{yu} \int_{zl}^{zu} \textrm{d}x\textrm{d}y\textrm{d}z K\left( \frac{x - x_i}{\sigma} \right)K\left( \frac{y - y_i}{\sigma} \right)K\left( \frac{z - z_i}{\sigma} \right) \]
where \(K\) is one of the kernel functions described on histogrambead and \(\sigma\) is a bandwidth parameter. The function \((s_i)\) can be any of the usual LESS_THAN, MORE_THAN, WITHIN etc that are used in all other multicolvars.
When AROUND is used with the DENSITY action the number of atoms in the specified region is calculated
The following commands tell plumed to calculate the average coordination number for the atoms that have x (in fractional coordinates) within 2.0 nm of the com of mass c1. The final value will be labeled s.mean.
By default the value of the calculated quantity can be referenced elsewhere in the input file by using the label of the action. Alternatively this Action can be used to calculate the following quantities by employing the keywords listed below. These quantities can be referenced elsewhere in the input by using this Action's label followed by a dot and the name of the quantity required from the list below.
Quantity | Keyword | Description |
lessthan | LESS_THAN | the number of cvs in the region of interest that are less than a certain threshold |
morethan | MORE_THAN | the number of cvs in the region of interest that are more than a certain threshold |
between | BETWEEN | the number of cvs in the region of interest that are within a certain range |
sum | SUM | the sum of all the colvars weighted by the function that determines if we are in the region |
mean | MEAN | the average values of the colvar in the region of interest |
ATOMS | the group of atoms that you would like to investigate. For more information on how to specify lists of atoms see Groups and Virtual Atoms |
ORIGIN | the atom whose vicinity we are interested in examining. For more information on how to specify lists of atoms see Groups and Virtual Atoms |
ATOM | an alternative to ORIGIN |
SIGMA | the width of the function to be used for kernel density estimation |
KERNEL | ( default=gaussian ) the type of kernel function to be used |
XLOWER | ( default=0.0 ) the lower boundary in x relative to the x coordinate of the atom (0 indicates use full extent of box). |
XUPPER | ( default=0.0 ) the upper boundary in x relative to the x coordinate of the atom (0 indicates use full extent of box). |
YLOWER | ( default=0.0 ) the lower boundary in y relative to the y coordinate of the atom (0 indicates use full extent of box). |
YUPPER | ( default=0.0 ) the upper boundary in y relative to the y coordinate of the atom (0 indicates use full extent of box). |
ZLOWER | ( default=0.0 ) the lower boundary in z relative to the z coordinate of the atom (0 indicates use full extent of box). |
ZUPPER | ( default=0.0 ) the upper boundary in z relative to the z coordinate of the atom (0 indicates use full extent of box). |
SERIAL | ( default=off ) do the calculation in serial. Do not parallelize |
OUTSIDE | ( default=off ) calculate quantities for colvars that are on atoms outside the region of interest |
SUM | ( default=off ) calculate the sum of all the quantities. |
MEAN | ( default=off ) calculate the average value of the colvar inside the region of interest |
DATA | the label of an action that calculates multicolvars. Weighted sums based on the location of the colvars calculated by this action will be calcualted |
LESS_THAN | calcualte the number of colvars that are inside the region of interest and that are less than a certain threshold |
MORE_THAN | calcualte the number of colvars that are inside the region of interest and that are greater that a certain threshold |
BETWEEN | calculate the number of colvars that are inside the region of interest and that have a CV value that is between a particular set of bounds |