The structure of pyrrhotite (Fe(1 - x)S with 0.05 < or = x < or = 0.125) has been reinvestigated in the framework of the superspace formalism. A common model with a centrosymmetric superspace group is proposed for the whole family. The atomic domains in the internal space representing the Fe atoms are parametrized as crenel functions that fulfil the closeness condition. The proposed model explains the x-dependent space groups observed and the basic features of the structures reported up to now. Our model yields for any x value a well defined ordered distribution of Fe vacancies in contrast to some of the structural models proposed in the literature. A new (3 + 1)-dimensional refinement of Fe(0.91)S using the deposited dataset [Yamamoto & Nakazawa (1982). Acta Cryst. A38, 79-86] has been performed as a benchmark of the model. The consistency of the proposed superspace symmetry and its validity for other compositions has been further checked by means of ab initio calculations of both atomic forces and equilibrium atomic positions in non-relaxed and relaxed structures, respectively.