Heat fluctuations in Ising models coupled with two different heat baths

Monte Carlo simulations of Ising models coupled to heat baths at two different temperatures are used to study a fluctuation relation for the heat exchanged between the two thermostats in a time tau. Different kinetics (single-spin-flip or spin-exchange Kawasaki dynamics), transition rates (Glauber or Metropolis), and couplings between the system and the thermostats have been considered. In every case the fluctuation relation is verified in the large tau limit, both in the disordered and in the low temperature phase. Finite-tau corrections are shown to obey a scaling behavior.