Hi Chol-Jun,
The Mishin1-Potential out of the paper cited before is provided in the LAMMPS code. This potential is fitted to the properties you are interested in. In order to extract the properties from the potential using LAMMPS you have to write some scripts and make some postprocessing. When I implemented the potential I did this and was able to reproduce the bulk properties to which to potential has been fitted.
More easily you can calculate bulk properties using GULP. For the GULP syntax I refer to its manual.
/Gerolf