  The computational efficiency of stochastic simulation algorithms is
notoriously limited by the kinetic trapping of the simulated trajectories
within low energy basins. Here we present a new method that overcomes kinetic
trapping while still preserving exact statistics of escape paths from the
trapping basins. The method is based on path factorization of the evolution
operator and requires no prior knowledge of the underlying energy landscape.
The efficiency of the new method is demonstrated in simulations of anomalous
diffusion and phase separation in a binary alloy, two stochastic models
presenting severe kinetic trapping.
