Motivated by the presence of a lattice of rotating molecular dipoles in the high temperature phase of methylammonium lead iodide, we investigate the ground state of a simple cubic lattice of dipoles interacting with each other via the dipole-dipole interaction and with an external field via the standard, linear dipole-field interaction. In the absence of an external field, the ground state is infinitely degenerate, and all the configurations in the ground state manifold are periodic along the three lattice axes with a period of two lattice sites. Using a 1000-dipole lattice as a unit cell in numerical simulations of an infinite simple cubic lattice, we determine the ground state dipole configurations in the presence of an external field. We then analyze the polarization, dipole orientation statistics and correlations in these configurations. Our calculations show that for some special directions of the external field the two-site periodicity in the dipole configurations is preserved, while in the general case this periodicity is lost and complex dipole configurations form under the influence of the external field. More specifically, for a general field direction, a sudden transition from two-site-periodic configurations to irregular configurations occurs at a finite value of the applied field strength.Spin models are ubiquitous in theoretical physics in the study of a variety of physical systems, and they can be used to understand a wide range of phenomena such as ordered states of matter, domain boundaries and phase transitions. [1] Recently, they have also served as natural models to investigate entanglement and related properties in the study of condensed matter systems from the viewpoint of quantum information. [2,3] In the field of light harvesting devices, the recently emerging hybrid organic-inorganic lead-halide perovskites have opened a novel application for spin models, because the organic component often has a net dipole moment. The most representative perovskite material in this area is methylammonium lead iodide (CH 3 NH 3 PbI 3 ), which in its high-temperature phase contains a simple cubic lattice of dipoles from the polar molecule CH 3 NH 3 that occupies the A site in the perovskite crystal structure. [4,5] Spin models have been used to explore the possible origin of hysteresis observed in solar cells, [6,7] to understand the dynamics of cation rotation [8] and to explain the enhancement of carrier conductivity at the interface between domains. [9] More recently, spin models have been used to explore the long-range order in perovskites [10] and the structure of the interface between domains with different ordering patterns across the interface. [11] Spin models apply more generally to other systems in nature, some examples being the ordering of the electric dipoles of water molecules in ice [12] and the dipolar ordering of small molecules trapped in regular cage structures. [13][14][15] Artificially fabricated magnetic superstructures can also