  We present a hybrid Brownian dynamics / Monte Carlo algorithm for simulating
solutions of highly entangled semiflexible polymers or filaments. The algorithm
combines a Brownian dynamics time-stepping approach with an efficient scheme
for rejecting moves that cause chains to cross or that lead to excluded volume
overlaps. The algorithm allows simulation of the limit of infinitely thin but
uncrossable threads, and is suitable for simulating the conditions obtained in
experiments on solutions of long actin protein filaments.
