Here we have implemented three (3) versions of bond boost acceleration method:
1. Compute_Bond_Boost_Force (): Apply boost to only one atom.
    a. type 1 in control file
    b. First, we find the bond with largest \epsilon ( most possible to break ),
       and record the atom as atom i;
    c. If \epsilon < q (usually 0.3), turn on boost, otherwise turn off boost;
    d. Select atom i to boost
2. Compute_Bond_Boost_Force_All (): Apply boost to all the atoms simutaneously,
   but the envelop function is applied on each atomic according to the tagged bonds
   of this specified atom.
   a. type 2 in control file;
   b. First, calculate the number (n) of radicals in the simulation box. Apply boost
      only if n = 0 (no radical in the simulation).
   c. Select atom sequentially. Calculate the \epsilon max. If \epsilon < q, apply
      boost to this atom. Go to next atom.
3. Coumpute_Bond_Boost_Force_All_Couple (): Apply boost to all the atoms simutaneously, and
   apply envelop function to all the atoms as if they are all coupled.
   a. type 3 in control file;
   b. as 2
   c. Loop all the atom and find the atom with max \epsilon; If \epsilon < q, apply
      boost to the atoms, and store the accumulation of boost energy.
   d. Apply the envelop function to the atom with max \epsilon.