In a singularly perturbed limit of small diffusivity ε of one of the two chemical species, equilibrium spike solutions to the Gray-Scott (GS) model on a bounded one-dimensional domain are constructed asymptotically using the method of matched asymptotic expansions. The equilibria that are constructed are symmetric k-spike patterns where the spikes have equal heights. Two distinguished limits in terms of a dimensionless parameter in the reaction-diffusion system are considered: the low feed-rate regime and the intermediate regime. In the low feed-rate regime, the solution branches of k-spike equilibria are found to have a saddle-node bifurcation structure. The stability properties of these branches of solutions are analyzed with respect to the large eigenvalues λ in the spectrum of the linearization. These eigenvalues, which have the property that λ = O(1) as ε → 0, govern the stability of the solution on an O(1) time scale. Precise conditions, in terms of the nondimensional parameters, for the stability of symmetric k-spike equilibrium solutions with respect to this class of eigenvalues are obtained. In the low feed-rate regime, it is shown that a large eigenvalue instability leads either to a competition instability, whereby certain spikes in a sequence are annihilated, or to an oscillatory instability (typically synchronous) of the spike amplitudes as a result of a Hopf bifurcation. In the intermediate regime, it is shown that only oscillatory instabilities are possible, and a scaling-law determining the onset and 9600 Garsington Road, Oxford, OX4 2DQ, UK.
22T. Kolokolnikov et al. of such instabilities is derived. Detailed numerical simulations are performed to confirm the results of the stability theory. It is also shown that there is an equivalence principle between spectral properties of the GS model in the low feed-rate regime and the Gierer-Meinhardt model of morphogenesis. Finally, our results are compared with previous analytical work on the GS model.