This work is the first application of classical atomistic theory to a comprehensive treatment of γ-Fe2O3 surfaces. The surface energy and attachment energy of several low-index surfaces of γ-Fe2O3 have been calculated by the classical atomistic simulation methods. A mean-field approximation involving scaling the short-range potentials, charges, and force constants of the iron ions by a fraction corresponding to the partial occupation of octahedral iron sites in the spinel structure was employed. Reconstruction of the polar surfaces was required to remove the dipole perpendicular to the surface and to achieve structural stability. This was accomplished through the addition of vacancies. The two-dimensional periodic calculations were carried out with the MARVIN code. The (112) surface consisting of iron and oxygen ions had the smallest relaxed surface energy of 1.86 J/m2, while several others including (001), (011), and (012) were within 0.1 J/m2 of this value. The calculated surface energies of these plane...