  We study the stability of cosmological scaling solutions within the class of
spatially homogeneous cosmological models with a perfect fluid subject to the
equation of state p_gamma=(gamma-1) rho_gamma (where gamma is a constant
satisfying 0 < gamma < 2) and a scalar field with an exponential potential. The
scaling solutions, which are spatially flat isotropic models in which the
scalar field energy density tracks that of the perfect fluid, are of physical
interest. For example, in these models a significant fraction of the current
energy density of the Universe may be contained in the scalar field whose
dynamical effects mimic cold dark matter. It is known that the scaling
solutions are late-time attractors (i.e., stable) in the subclass of flat
isotropic models. We find that the scaling solutions are stable (to shear and
curvature perturbations) in generic anisotropic Bianchi models when gamma <
2/3. However, when gamma > 2/3, and particularly for realistic matter with
gamma >= 1, the scaling solutions are unstable; essentially they are unstable
to curvature perturbations, although they are stable to shear perturbations. We
briefly discuss the physical consequences of these results.
