new plot: magnetization squared