Abstract Shockwave compression of methane was simulated using the molecular dynamics method with Hugoniostat and MSST approach. The calculation was carried out employing ReaxFFlg potential. We show comparison of simulated Hugoniots for pressure range 0,1-44 GPa and times required for simulation run of two different methods. We also compare two ReaxFF potentials: ReaxFF(2008) and ReaxFF-lg