A modified embedded-atom method (MEAM) interatomic potential for the Fe−Cu binary system has been developed using previously developed MEAM potentials of Fe and Cu. The Fe−Cu potential was determined by fitting to data on the mixing enthalpy and the composition dependencies of the lattice parameters in terminal solid solutions. The potential gives a value of 0.65eV for the dilute heat of solution and reproduces the increase of lattice parameter of Fe with addition of Cu in good agreement with experiments. The potential was used to investigate the primary irradiation defect formation in pure Fe and Fe−0.5at.%Cu alloy by a molecular dynamics cascade simulation study with a PKA energy of 2keV at 573K. In the paper (Lee et al., Phys. Rev. B, 2005), a tendency for self-interstitial atom-Cu binding, the formation of mixed (Fe−Cu) dumbbells and even Cu−Cu dumbbells was observed.