Modeling Reactive Rarefied Systems Using A Novel Quasi-Particle Boltzmann Solver