m = IMP.kernel.Model()# Brownian Dynamicsbd=IMP.atom.BrownianDynamics(m)bd.set_maximum_time_step(1000)for i in range(1, nrounds): e = bd.optimize(10)