I will present a numerical method that we are currently developing in the NAMD software and which aim at simulation reactive trajectories. Numerical and theoretical results on this technique will be described.