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

Calculate the gradient of the average value of a multicolvar value

This command allows you to calculate the collective variable discussed in [49].

Examples

The input below calculates the gradient of the density of atoms in the manner described in [49] in order to detect whether or not atoms are distributed uniformly along the x-axis of the simulation cell.

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

d1: DENSITY SPECIESthis keyword is used for colvars such as coordination number. 
=1-50 The DENSITY action with label d1
s1: GRADIENT ORIGINwe will use the position of this atom as the origin in our calculation. 
=1 DATAcompulsory keyword 
The multicolvar that calculates the set of base quantities that we are interested
in 
=d1 DIRcompulsory keyword ( default=xyz )
the directions in which we are calculating the gradient. 
=x NBINScompulsory keyword 
number of bins to use in each direction for the calculation of the gradient 
=4 SIGMAcompulsory keyword ( default=1.0 )
the width of the function to be used for kernel density estimation 
=1.0 The GRADIENT action with label s1 calculates a single scalar value
PRINT ARGthe input for this action is the scalar output from one or more other actions. 
=s1 FILEthe name of the file on which to output these quantities 
=colvar The PRINT action with label

The input below calculates the coordination numbers of the 50 atoms in the simulation cell. The gradient of this quantity is then evaluated in the manner described using the equation above to detect whether the average values of the coordination number are uniformly distributed along the x-axis of the simulation cell.

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

d2: COORDINATIONNUMBER SPECIESthis keyword is used for colvars such as coordination number. 
=1-50 SWITCHThis keyword is used if you want to employ an alternative to the continuous switching
function defined above. 
={RATIONAL R_0=2.0}  MORE_THANcalculate the number of variables more than a certain target value. 
={EXP R_0=4.0}  The COORDINATIONNUMBER action with label d2 calculates a single scalar value
s2: GRADIENT ORIGINwe will use the position of this atom as the origin in our calculation. 
=1 DATAcompulsory keyword 
The multicolvar that calculates the set of base quantities that we are interested
in 
=d2 DIRcompulsory keyword ( default=xyz )
the directions in which we are calculating the gradient. 
=x NBINScompulsory keyword 
number of bins to use in each direction for the calculation of the gradient 
=4 SIGMAcompulsory keyword ( default=1.0 )
the width of the function to be used for kernel density estimation 
=1.0 The GRADIENT action with label s2 calculates a single scalar value
PRINT ARGthe input for this action is the scalar output from one or more other actions. 
=s2 FILEthe name of the file on which to output these quantities 
=colvar The PRINT action with label

Glossary of keywords and components

The atoms involved can be specified using

ORIGIN

we will use the position of this atom as the origin in our calculation. For more information on how to specify lists of atoms see Groups and Virtual Atoms

Compulsory keywords

DATA	The multicolvar that calculates the set of base quantities that we are interested in
DIR	( default=xyz ) the directions in which we are calculating the gradient. Should be x, y, z, xy, xz, yz or xyz
NBINS	number of bins to use in each direction for the calculation of the gradient
SIGMA	( default=1.0 ) the width of the function to be used for kernel density estimation
KERNEL	( default=gaussian ) the type of kernel function to be used

Options

NUMERICAL_DERIVATIVES	( default=off ) calculate the derivatives for these quantities numerically
NOPBC	( default=off ) ignore the periodic boundary conditions when calculating distances
SERIAL	( default=off ) do the calculation in serial. Do not use MPI
LOWMEM	( default=off ) lower the memory requirements
TIMINGS	( default=off ) output information on the timings of the various parts of the calculation