Standard sampling algorithms for lattice QCD suffer from topology freezing (or critical slowing down) when approaching the continuum limit, thus leading to poor sampling of the distinct topological sectors. I will present a modified Hamiltonian Monte Carlo (HMC) algorithm that triggers topological sector jumps during the assembly of Markov chain of lattice configurations. We study its performance in the 2D Schwinger model and compare it to alternative methods, such as fixing topology or master field. We then briefly discuss the difficulties of the algorithm in a SU(2) gauge model in 4D.