The properties of (Mg,Fe)SiO₃ perovskite at lower mantle conditions are still not well understood, and particular attention has recently been given to determining the Fe spin state. A major challenge in spin states studies is interpretation of Mössbauer spectra to determine the electronic structure of iron under extreme conditions. In this paper ab initio methods are used to predict quadrupole splitting values of high-, intermediate- and low-spin Fe²⁺ and Fe³⁺ in perovskite, as a function of pressure and composition. The calculations in (Mg_0.75Fe_0.25)SiO₃ yield quadrupole splitting values in the range of 0.7–1.7 mm/s for all spin and valence states except high-spin Fe²⁺, which has two possible quadrupole splittings, 2.3 and 3.3 mm/s. The unexpected multiple quadrupole splitting values for high-spin Fe²⁺ are explained in terms of small changes in local structure and d-orbital occupations. The computed results are applied to interpret existing perovskite Mössbauer data for iron's spin state.