  We present a new computational methodology for the investigation of gel
electrophoresis of polyelectrolytes. We have developed the method initially to
incorporate sliding motion of tight parts of a polymer pulled by an electric
field into the bond fluctuation method (BFM). Such motion due to tensile force
over distances much larger than the persistent length is realized by non-local
movement of a slack monomer at an either end of the tight part. The latter
movement is introduced stochastically. This new BFM overcomes the well-known
difficulty in the conventional BFM that polymers are trapped by gel fibers in
relatively large fields. At the same time it also reproduces properly
equilibrium properties of a polymer in a vanishing filed limit. The new BFM
thus turns out an efficient computational method to study gel electrophoresis
in a wide range of the electric field strength.
