We consider Schloegl models where particles on a square grid annihilate at rate , and are created at rate that related to neighboring particles in . Simulation reveals a discontinuous transition between populated and vacuum states, but equistability determined by stationarity of planar interfaces between states depends on interface orientation at least for smaller . Behavior for large follows from continuum equations. These depend on interface orientation and -shape, but a unique equistability point emerges  imposing a Gibbs phase rule.