A two-dimensional atomistic realization of Schlögl’s second model for autocatalysis is implemented and studied on a square lattice as a prototypical nonequilibrium model with first-order transition. The model has no explicit symmetry and its phase transition can be viewed as the nonequilibrium counterpart of liquid-vapor phase separations. We show some familiar concepts from study of equilibrium systems need to be modified. Most importantly, phase coexistence can be a generic feature of the model, occurring over a finite region of the parameter space. The first-order transition becomes continuous as a temperature-like variable increases. The associated critical behavior is studied through Monte Carlo simulations and shown to be in the two-dimensional Ising universality class. However, some common expectations regarding finite-size corrections and fractal properties of geometric clusters for equilibrium systems seems to be inapplicable.