As another property impacted by temperature gradients, we computed the distribution of cholesterol along the membrane normal. case of the original model, or significant differences between the two models (see Figures and). However, the system with the least strict settings (and largest ∆T chol ) tends to be an outlier. Therefore, we discourage the use of lincs order= 4 in combination with a long time step of ∆t = 30 fs in the investigated system. We also note that while there are no systematic differences between the two cholesterol topologies, the quantities determined using the optimized model have smaller variances (Figures and). Moreover, if the