The use of optimized Monte Carlo methods for studying spin glasses

Abstract
We start from recently published numerical data by Hatano and Gubernatis to discuss properties of convergence to equilibrium of optimized Monte Carlo methods (bivariate multi-canonical and parallel tempering). We show that these data are not thermalized, and they lead to an erroneous physical picture. We shed some light on why the bivariate multi-canonical Monte Carlo method can fail.