if (Math.random() < Math.exp(-delta_energy / currT)) {
            acc += 1;
        } else {
            // move back to old coordinates
            lab[i].x = x_old;