For phase separating lipid bilayers, Martini 2 simulations with typical parameter settings have recently found to suffer from substantial artificial temperature gradients across the phase boundaries that impacted other physical properties of the system such as the ratio of diffusion coefficients between the saturated and unsaturated lipids, the degree of phase separation, or the distribution of cholesterol along the membrane normal direction (Fig- ). The origin of the artifact was traced back to insufficient convergence of the highly coupled bond constraints in cholesterol, one of the major components in such bilayers.