Model reduction of fast-slow chemical reaction networks based on the quasi-steady state
approximation fails when the fast subsystem has first integrals. We call these first integrals …