Pd-Cu binary systems based on the second nearest-neighbor modified embedded-atom method formalism. This potential reproduces various fundamental properties of the alloys (the structural, elastic and thermodynamic properties of compound and solution phases, and order-disorder transition temperature) in reasonable agreements with experimental data, first-principles calculations and CALPHAD assessments.