We use the Wang-Landau algorithm to calculate a density of states for an asymmetric Ising model on a triangular lattice with two- and three-body interactions in an external field. An accurate density of states allows us to determine the phase diagram and to study the critical behavior of this model at and near the critical endpoint. We observe a divergence of the curvature of the spectator phase boundary at the critical endpoint in accordance with theoretical predictions.

We review recently developed decomposition algorithms for molecular dynamics and spin dynamics simulations of many-body systems. These methods are time reversible, symplectic, and the error in the total energy thus generated is bounded. In general, these techniques are accurate for much larger time steps than more standard integration methods. Illustrations of decomposition algorithms performance are shown for spin dynamics simulations of a Heisenberg ferromagnet.