Abstract: An interatomic potential for pure phosphorus, an element that has van der Waals, covalent and metallic bonding character, simultaneously, has been developed for the purpose of application to metal-phosphorus systems. As a simplification, the van der Waals interaction, which is less important in metal-phosphorus systems, was omitted in the parameterization process and potential formulation. On the basis of the second-nearest-neighbor modified embedded-atom method (2NN MEAM) interatomic potential formalism applicable to both covalent and metallic materials, a potential that can describe various fundamental physical properties of a wide range of allotropic or transformed crystalline structures of pure phosphorus could be developed. The potential was then extended to the Fe-P binary system describing various physical properties of intermetallic compounds, bcc and liquid alloys, and also the segregation tendency of phosphorus on grain boundaries of bcc iron, in good agreement with experimental information. The suitability of the present potential and the parameterization process for atomic scale investigations about the effects of various non-metallic impurity elements on metal properties is demonstrated.