Biased random walks on a lattice: exact numerical method to study the effect of alternating fields in disordered and asymmetric systems of obstacles.

The migration of a particle in a system of obstacles under the action of an external field is often modeled using lattice Monte Carlo algorithms. For example, such simulation methods have been used to study the electrophoresis of charged molecules in sieving gels and the separation of particles using ratchet systems. In the case of constant fields or low-frequency alternating fields, the Monte Carlo simulation method can be mapped onto a numerical or algebraic matrix problem that can be solved exactly. In this Rapid Communication, we generalize this matrix approach to treat periodic time-dependent fields. The evolution of the spatial distribution function during a period is computed using a sequence of transfer matrices, and a steady-state closure relation allows us to calculate the exact mean velocity of the particle during a complete cycle. As an example, we examine the properties of a simple spatially asymmetric ratchet system in the presence of periodic alternating fields (symmetric and asymmetric) as well as random telegraph signals.