I am using v-rescale and conserved energy drifts linearly, which seems consistent with the discussion here: "Conserved En." is increasing linearly with simulation time - Is this normal? though I did not check the amount of drift…
Not that many people want to see it in their simulations, but, at least in principle, at 1M spontaneous crystallization is actually reasonably physical, so, according to some of my colleagues, “maybe your Gromacs finally started working correctly.” :)
In any case, I am moving back to a 2022 version to check. Will report.