add a lockstep rate matrix that is treated not so well by matrix exponential and stationary distribution calculation