Magnetization is measured in experiments for both ferromagnetic and antiferromagnetic materials to investigate the magnetic properties of materials, and the susceptibility of total magnetization as a function of external field is used to determine Neel or Cued temperatures. In the Monte Carlo simulation, it is important to define the proper order parameters to describe the spin model, where the magnetization is used as an order parameter for ferromagnetic spin model and the staggered magnetization is used as an order parameter for the antiferromagnetic spin model without geometrical frustration. However, it is difficult to define an order parameter for frustrated spin models. We perform the Monte Carlo simulation for the antiferromagnetic Heisenberg spin model using the damage spreading as an order parameter, and also perform simulation using both the magnetization and the staggered magnetization as order parameters. Then, we measure the critical temperatures and the critical exponents on three-dimensional regular lattice estimated by different order parameters, and then study the dependency of critical behaviors on different order parameters for the antiferromagnetic Heisenberg spin models.