We developed a python script using the MDAnalysis 19 package to compute λ max from a single molecular configuration. The use of a single configuration is justified only if the constrained particles do not undergo significant fluctuations during their motion. This condition holds for the Martini 2 cholesterol model, as the two coupled triangles never substantially deviate from coplanarity. Based on the framework for the generation of equimomental systems introduced in the previous section, we minimized λ max of A n computed by the script with respect to equimomental configurations of a fixed number of rigid sites. In this way, we aimed at reducing the lincs order required for properly constraining Martini 2 cholesterol.