To further estimate the balance of bonded and non-bonded interactions in these models, we evaluate a posteriori their ability to describe phase separation in DOPC/DPPC mixtures known to simultaneously exhibit two phase states at 298K (this phase separation is not correctly described using Martini lipid models version 3.0.0). Simulating a bilayer composed of 1152 randomly dispersed lipids at 10/90% mass DOPC/DPPC and starting from a configuration in the liquid phase, the nucleation of the gel phase is observed after ~100 ns of equilibration and the equilibrated system exhibits liquid/gel phase separation after ~1 µs (Fig. ). To characterize the phase separation we employ the LENS descriptor , which allows to evaluate locally the changes in the environment of a molecule across trajectory frames and here enables classifying lipids into phase states. The proportion of lipids in the gel and liquid phases is stable in an additional 10 µs of production simulation, during which the two phases constantly exchange lipids (Fig. : stars, Supplementary Movie 1). Re-iterating this experiment with 20/80% mass DOPC/DPPC in the system, we did not observe gel phase nucleation after 1 µs of equilibration. Lowering the simulation temperature to 293K, stable phase separation could be observed also at 20/80% mass DOPC/DPPC after ~1 µs of equilibration (Fig. ), again with a constant exchange of lipids in between the two phases (Fig. squares).