The equations describing the distributions and concentrations of trace species are nonlinear and may thus possess more than one solution. Several authors have suggested that the steady‐state equations describing tropospheric and stratospheric chemistry may have multiple solutions, but the existence of such solutions has not been completely demonstrated. This paper develops methods for searching for multiple physical solutions to chemical continuity equations and applies these to subsets of equations describing tropospheric chemistry. The calculations are carried out with a box model and use two basic strategies. The first strategy is a “search” method. This involves fixing model parameters at specified values, choosing a wide range of initial guesses at a solution, and using a Newton‐Raphson technique to determine if different initial points converge to different solutions. The second strategy involves a set of techniques known as homotopy methods. These do not require an initial guess, are globally convergent, and are guaranteed, in principle, to find all solutions of the continuity equations. The first method is efficient but essentially “hit or miss” in the sense that it cannot guarantee that all solutions which may exist will be found. The second method is computationally burdensome but can, in principle, determine all the solutions of a photochemical system. Multiple solutions have been found for models that contain a basic complement of photochemical reactions involving Ox, HOx, NOx, and CH4. In the present calculations, transitions occur between stable branches of a multiple solution set as a control parameter is varied. These transitions are manifestations of hysteresis phenomena in the photochemical system and may be triggered by increasing the NO flux or decreasing the CH4 flux from current mean tropospheric levels.