Quantitative quality evaluation of the eect of CHOL on membrane properties: To ease the evaluation of simulations against experimental data with non-matching CHOL concentrations, we interpolated the eect of CHOL in simulations and experiments for CHOL concentrations ranging from 0% to 46%. 2D matrices were created by interpolation for the oleate and palmitate chain order parameters (as a function of carbon atoms in the acyl chains and CHOL concentration) using the interp2 function in Matlab based on linear interpolation. Similar matrices were also generated for the electron density proles (as a function of normal distance from bilayer center and CHOL concentration) but with the scatteredInterpolant function in Matlab based on Delaunay triangulation. 66 Linear 1D interpolation with the interp1 function in Matlab was used for the CHOL-dependence of the 1 st and 2 nd form factor minima locations and diusion coecients. These interpola- tions were then used to calculate deviations (in %) from experimental values across CHOL concentrations to quantify the quality of the lipid force elds. For the 2D matrices, the absolute values of the dierences between matrices from simulations and experiments were rst calculated. The averages of dierences over the carbon atoms in the acyl chain (order parameters) or across the entire z coordinate in the simulation (density proles) were then calculated. This resulted in 1D deviation vectors as a function of CHOL concentration.