Brownian dynamics algorithms integrate Langevin equations numerically and allow to probe long time scales in simulations. A common requirement for such algorithms is that interactions in the system should vary little during an integration time step; therefore, computational efficiency worsens as the interactions become steeper. In the extreme case of hard-body interactions, standard numerical integrators become ill defined. Several approximate schemes have been invented to handle such cases, ...
