Probability distributions of the work in the 2D-Ising model

Probability distributions of the magnetic work are computed for the 2D Ising model by means of Monte Carlo simulations. The system is first prepared at equilibrium for three temperatures below, at and above the critical point. A magnetic field is then applied and grown linearly at different rates. Probability distributions of the work are stored and free energy differences computed using the Jarzynski equality. Consistency is checked and the dynamics of the system is analyzed. Free energies and dissipated works are reproduced with simple models. The critical exponent $\delta$ is estimated in an usual manner.