introduce reaction operator BindingReaction to simulate stationary binding sites

The operator BindingReaction has a single state variable that accounts for the number of bound molecules. Molecules are bound with the given rate (using the rho-sigma method) and can dissociate from the site with rate diss_rate. One can also specify an optional initial state.