We study topological crystalline insulators doped with magnetic impurities, in which ferromagnetism at the surface lowers the electronic energy by spontaneous breaking of a crystalline symmetry. The number of energetically equivalent groundstates is sensitive to the crystalline symmetry of the surface, as well as the precise density of electrons at the surface. We show that for a SnTe model in the topological state, magnetic states can have two-fold, six-fold symmetry, or eight-fold degenerate minima. We compute spin stiffnesses within the model to demonstrate the stability of ferromagnetic states, and consider their ramifications for thermal disordering. Possible experimental consequences of the surface magnetism are discussed.PACS numbers: 73.20. At,75.70.Rf,75.30.Gw Introduction -Topological crystalline insulators (TCI's) are a class of materials in which the energy bands can host non-trivial topology protected by a crystalline symmetry [1]. These systems support surface states [2] which remain gapless provided the crystal symmetry is unbroken, and are believed to present themselves in (Sn,Pb)Te and related alloys [3][4][5][6][7][8]. Interesting effects may arise when the symmetry protecting a topological band structure is broken. In topological insulators protected by time-reversal symmetry (TRS), magnetic impurities on a surface break this symmetry and form collective states [9][10][11][12][13], which may be understood in terms of a gap opening in the surface spectrum [14].In contrast, TCI's are not protected by TRS, so the loss of this symmetry does not by itself energetically favor ordering of magnetic moments [15,16]. However, a uniform magnetization can undermine one or more relevant crystalline symmetries [17,18]. Indeed, the most common such symmetry is reflection across a mirror plane, of which there can be several. We show below that spontaneous surface magnetization opens a maximal gap when oriented along axes dictated by the bulk symmetries of the system. For a generic surface with a single mirror plane, there are two surface Dirac points at different momenta and energies [19], and in such cases at low temperature this results in a metallic, Ising-like ferromagnet, with the easy axis determined by the chemical potential µ. Importantly, the number of degenerate low-energy directions is enhanced for surfaces with further symmetries. Rotational symmetries in particular yield multiple mirror planes, and connect distinct surface Dirac cones to one another, yielding a multiplicity of easy axis directions. For sufficiently high symmetry, all the surface Dirac points may be related by symmetry operations, resulting in a fully gapped surface spectrum and a large number of groundstate orientations.To illustrate this physics, we present detailed calculations for the (111) surface of (Sn,Pb)Te [8,[20][21][22], using a known model Hamiltonian [3,23]. The (111) surface states are characterized in this system by four surface Dirac points, one at theΓ point and one at each of three Fig. 1(a).] When the system...