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:DISTANCEATOMS=3,5the pair of atom that we are calculating the distance between.external:EXTERNALARG=compulsory keywordthe labels of the scalars on which the bias will actd1FILE=bias.gridcompulsory keywordthe name of the file containing the external potential.

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:TORSIONATOMS=4,5,6,7the four atoms involved in the torsional anglet2:TORSIONATOMS=6,7,8,9the four atoms involved in the torsional angleext:EXTERNALARG=compulsory keywordthe labels of the scalars on which the bias will actt1,t2FILE=bias2.gridcompulsory keywordthe name of the file containing the external potential.

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 |