WHAM_WEIGHTS
This is part of the analysis module

Calculate and output weights for configurations using the weighted histogram analysis method.

This shortcut action allows you to calculate and output weights computed using the weighted histogram analysis technique. For more detail on how this technique works see REWEIGHT_WHAM

Examples

The following input can be used to analyze the output from a series of umbrella sampling calculations. The trajectory from each of the simulations run with the different biases should be concatenated into a single trajectory before running the following analysis script on the concatenated trajectory using PLUMED driver. The umbrella sampling simulations that will be analyzed using the script below applied a harmonic restraint that restrained the torsional angle involving atoms 5, 7, 9 and 15 to particular values. The script below calculates the reweighting weights for each of the trajectories and then applies the binless WHAM algorithm to determine a weight for each configuration in the concatenated trajectory.

Click on the labels of the actions for more information on what each action computes
tested on master
#SETTINGS NREPLICAS=4
phi: TORSION 
ATOMS
the four atoms involved in the torsional angle
=5,7,9,15 rp: RESTRAINT
ARG
the arguments on which the bias is acting.
=phi
KAPPA
compulsory keyword ( default=0.0 ) specifies that the restraint is harmonic and what the values of the force constants on each of the variables are
=50.0
AT
compulsory keyword the position of the restraint
=@replicas:{ -3.00000000000000000000 -1.45161290322580645168 .09677419354838709664 1.64516129032258064496 } ... WHAM_WEIGHTS
BIAS
compulsory keyword ( default=*.bias ) the value of the biases to use when performing WHAM
=rp.bias
TEMP
the temperature at which the simulation was run
=300
FILE
compulsory keyword the file on which to output the WHAM weights
=wham-weights

The script above must be run with multiple replicas using the following command:

mpirun -np 4 plumed driver --mf_xtc alltraj.xtc --multi 4
Glossary of keywords and components
Compulsory keywords
BIAS ( default=*.bias ) the value of the biases to use when performing WHAM
STRIDE ( default=1 ) the frequency with which the bias should be stored to perform WHAM
FILE the file on which to output the WHAM weights
Options
TEMP the temperature at which the simulation was run
FMT the format to use for the real numbers in the output file