BRIDGE
This is part of the multicolvar module

Calculate the number of atoms that bridge two parts of a structure

This quantity calculates:

\[ f(x) = \sum_{ijk} s_A(r_{ij})s_B(r_{ik}) \]

where the sum over \(i\) is over all the ``bridging atoms" and \(s_A\) and \(s_B\) are switchingfunction.

The atoms involved can be specified using
ATOMS the atoms involved in each of the collective variables you wish to calculate. Keywords like ATOMS1, ATOMS2, ATOMS3,... should be listed and one CV will be calculated for each ATOM keyword you specify (all ATOM keywords should define the same number of atoms). The eventual number of quantities calculated by this action will depend on what functions of the distribution you choose to calculate. You can use multiple instances of this keyword i.e. ATOMS1, ATOMS2, ATOMS3...
Or alternatively by using
BRIDGING_ATOMS The list of atoms that can form the bridge between the two interesting parts of the structure.
GROUPA The list of atoms that are in the first interesting part of the structure
GROUPB The list of atoms that are in the second interesting part of the structure
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 parallelize
LOWMEM ( default=off ) lower the memory requirements
TIMINGS ( default=off ) output information on the timings of the various parts of the calculation
VERBOSE

( default=off ) write a more detailed output

SWITCH The parameters of the two switchingfunction in the above formula
SWITCHA The switchingfunction on the distance between bridging atoms and the atoms in group A
SWITCHB

The switchingfunction on the distance between the bridging atoms and the atoms in group B

Examples

The following example instructs plumed to calculate the number of water molecules that are bridging betweeen atoms 1-10 and atoms 11-20 and to print the value to a file

BRIDGE BRIDGING_ATOMS=100-200 GROUPA=1-10 GROUPB=11-20 LABEL=w1
PRINT ARG=a1.mean FILE=colvar