Despite the many advantages (e.g., suitable band gap, exceptional optical absorptivity, earth abundance) of pyrite as a photovoltaic material, its low open circuit voltage (OCV) has remained the biggest challenge preventing its use in practical devices. Two of the most widely-accepted reasons for the cause of the low OCV are: (i) Fermi level pinning due to intrinsic surface states that appear as gap states; (ii) the presence of the metastable polymorph, marcasite. In this paper, we investigate these claims, via density-functional theory (DFT), by examining the electronic structure, bulk, surface, and interfacial energies of pyrite and marcasite. Regardless of whether the Hubbard U correction is applied, the intrinsic {100} surface states are found to be of d z 2 character, as expected from ligand field theory. However, they are not gap states, but located at the conduction band edge. Thus, ligand field splitting at the symmetry-broken surface cannot be the sole cause of the low OCV. We also investigate epitaxial growth of marcasite on pyrite. Based on the surface, interfacial, and strain energies of pyrite and marcasite, we find from our model that only one layer of epitaxial growth of marcasite is thermodynamically favorable. Within all methods used (LDA, GGA-PBE, GGA-PBE+U , GGA-AM05, GGA-AM05+U , HSE06, ∆-sol), the marcasite Kohn-Sham gap is not less than the pyrite Kohn-Sham gap, and is even larger than the experimental marcasite gap. Moreover, gap states are not observed at the pyrite-marcasite interface. We conclude that intrinsic surface states or presence of marcasite is unlikely to undermine the photovoltaic performance of pyrite.