We study magnetism in simple models for rare earth quasicrystals, by considering Ising spins on a quasiperiodic tiling, coupled via RKKY interactions. Computing these interactions from a tight-binding model on the tiling, we find that they are frustrated and strongly dependent on the local environment. Although such features are often associated with spin glass behaviour, we show using Monte Carlo simulations that the spin system has a phase transition to a low-temperature state with long-range quasiperiodic magnetic order.PACS numbers: 75.50. Kj, Understanding how the magnetic properties of complex materials arise from their atomic structure represents one of the fundamental challenges of condensed matter physics. In the context of quasicrystals [1] this leads to the question of how their magnetism [2,3] is influenced by their unusual electronic properties, which include a pseudogap at the Fermi energy [4,5], multifractal wave functions [6] and anomalous electronic transport [7,8]. Magnetic rare earth (R) quasicrystals in particular serve to motivate simple models, as they contain welldefined local moments at concentrations of 5-10% interacting via long-range RKKY interactions [9,10]. Examples of these quasicrystals In this paper we take a two-step theoretical approach to studying such systems. First, we compute RKKY interactions from a tight-binding Hamiltonian defined on a two-dimensional quasiperiodic tiling. Second, we examine the statistical mechanics of Ising spins with these interactions, using extensive Monte Carlo simulations. While there has been much previous work studying spin models for magnetic quasicrystals [13][14][15][16][17][18][19][20], it has been based on choices of interaction (nearest neighbour exchange [13,14], dipolar and further neighbour exchange [15][16][17], or RKKY interactions with a form taken from periodic systems [18][19][20]) that do not incorporate the unique coupling of electronic to magnetic properties that is to be expected in quasicrystals. Our model is designed to address this aspect of the physics in a simple way.We find substantial differences between RKKY interactions on a quasiperiodic tiling and those in simple metals [21]. While both ferromagnetic and antiferromagnetic interactions arise in the tiling, they have no well-defined spatial period since there is no Fermi wavevector. And as there is no translational invariance, the RKKY coupling between pairs of sites is a function not only of their separation, but is also strongly dependent on the local environment. One might expect that the combination of frustration and aperiodicity would lead to spin freezing [22], and in other settings quasiperiodic systems are known to behave like random ones [23]. Our results, however, exclude canonical spin glass behaviour, via the scaling of domain wall energy with length. We demonstrate that instead there is quasiperiodic spin order at low temperature, by examining the Fourier transform of spin configurations.We model conduction electrons by a pure hopping Hamiltonian H el = l,...