As the high entropy alloy has become a more and more popular research field, numerous experimental research work has been published with different applications, such as tension, fatigue, and fracture behavior, irradiation damage, oxidation, corrosion…etc. Reliable molecular dynamics simulations are thus needed to benchmark with experimental results for larger-scale calculations, such as the slip of dislocations, grain-boundary and irradiation cascade. However, none of the current mixing method from the embedded atom method (EAM) or Lennard-Jones (LJ) potential can accurately estimate the lattice structure and stacking fault energy. Therefore an accurate potential function is urgently required. We choose the modified embedded atom method (MEAM) as the fitting form, which can be expressed as the summation of the embedding function and pair potential between any two atoms. We will construct a database for different alloy compositions and calculate the total energy by the density function theory (DFT). The actual fitting process would be completed by the General Utility Lattice Program (GULP) code package.