EXTERNAL

This is part of the bias module

Calculate a restraint that is defined on a grid that is read during start up

Examples

The following is an input for a calculation with an external potential that is defined in the file bias.dat and that acts on the distance between atoms 3 and 5.

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

d1: DISTANCE 
ATOMS
the pair of atom that we are calculating the distance between. 
=3,5 The DISTANCE action with label d1 calculates a single scalar value
external: EXTERNAL 
ARG
compulsory keyword 
the labels of the scalars on which the bias will act 
=d1 
FILE
compulsory keyword 
the name of the file containing the external potential. 
=bias.grid The EXTERNAL action with label external calculates a single scalar value

The bias.grid will then look something like this:

#! FIELDS d1 external.bias der_d1
#! SET min_d1 1.14
#! SET max_d1 1.32
#! SET nbins_d1 6
#! SET periodic_d1 false
   1.1400   0.0031   0.1101
   1.1700   0.0086   0.2842
   1.2000   0.0222   0.6648
   1.2300   0.0521   1.4068
   1.2600   0.1120   2.6873
   1.2900   0.2199   4.6183
   1.3200   0.3948   7.1055

This should then be followed by the value of the potential and its derivative at 100 equally spaced points along the distance between 0 and 1.

You can also include grids that are a function of more than one collective variable. For instance the following would be the input for an external potential acting on two torsional angles:

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

t1: TORSION 
ATOMS
the four atoms involved in the torsional angle 
=4,5,6,7 The TORSION action with label t1 calculates a single scalar value
t2: TORSION 
ATOMS
the four atoms involved in the torsional angle 
=6,7,8,9 The TORSION action with label t2 calculates a single scalar value
ext: EXTERNAL 
ARG
compulsory keyword 
the labels of the scalars on which the bias will act 
=t1,t2 
FILE
compulsory keyword 
the name of the file containing the external potential. 
=bias2.grid The EXTERNAL action with label ext calculates a single scalar value

The file bias2.grid for this calculation would need to look something like this:

#! FIELDS t1 t2 ext.bias der_t1 der_t2
#! SET min_t1 -pi
#! SET max_t1 pi
#! SET nbins_t1 3
#! SET periodic_t1 true
#! SET min_t2 -pi
#! SET max_t2 pi
#! SET nbins_t2 3
#! SET periodic_t2 true
 -3.141593 -3.141593 0.000000 -0.000000 -0.000000
 -1.047198 -3.141593 0.000000 0.000000 -0.000000
 1.047198 -3.141593 0.000000 -0.000000 -0.000000

 -3.141593 -1.047198 0.000000 -0.000000 0.000000
 -1.047198 -1.047198 0.007922 0.033185 0.033185
 1.047198 -1.047198 0.007922 -0.033185 0.033185

 -3.141593 1.047198 0.000000 -0.000000 -0.000000
 -1.047198 1.047198 0.007922 0.033185 -0.033185
 1.047198 1.047198 0.007922 -0.033185 -0.033185

This would be then followed by 100 blocks of data. In the first block of data the value of t1 (the value in the first column) is kept fixed and the value of the function is given at 100 equally spaced values for t2 between \(-pi\) and \(+pi\). In the second block of data t1 is fixed at \(-pi + \frac{2pi}{100}\) and the value of the function is given at 100 equally spaced values for t2 between \(-pi\) and \(+pi\). In the third block of data the same is done but t1 is fixed at \(-pi + \frac{4pi}{100}\) and so on until you get to the one hundredth block of data where t1 is fixed at \(+pi\).

Please note the order that the order of arguments in the plumed.dat file must be the same as the order of arguments in the header of the grid file.

Glossary of keywords and components

Description of components

By default this Action calculates the following quantities. 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	Description
bias	the instantaneous value of the bias potential

Compulsory keywords

ARG	the labels of the scalars on which the bias will act
FILE	the name of the file containing the external potential.
SCALE	( default=1.0 ) a factor that multiplies the external potential, useful to invert free energies

Options

NUMERICAL_DERIVATIVES	( default=off ) calculate the derivatives for these quantities numerically
NOSPLINE	( default=off ) specifies that no spline interpolation is to be used when calculating the energy and forces due to the external potential
SPARSE	( default=off ) specifies that the external potential uses a sparse grid