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.