text: dependence of diﬀerent samples during the measurement 
  stage. For each temperature we have usually performed 
  three independent simulations, each one containing at 
  least 2×105 measurements, taken after discarding up to 
  5×104 Monte Carlo steps in order to assure thermal equi- 
  libration. 
  
  In the proximity of the critical region the multiple his- 
  togram (MH) technique was also employed21, as it allows 
  us to estimate the physical observables of interest over a 
  whole temperature range in a substantially continuous 
  way by interpolating results obtained from sets of simu- 
  lations performed at some diﬀerent temperatures. 
  
  For all the quantities of interest, the average value and 
  the error estimate were obtained by the bootstrap re- 
  sampling method22 given that, as pointed out in Ref. 23, 
  for a large enough number of measurements, this method 
  turns out to be more accurate than the usual blocking 
  technique. In our implementation, we pick out randomly 
  a sizable number of measurements (typically, between 1 
  and 1×103 for the single simulation, and between 1 and 
  5×104 for the MH technique), and iterate the re-sampling 
  at least one hundred times. 