In phospholipid bilayer/monolayer systems the solid like gel (𝐿 𝛽 ) or LC states have essentially frozen lipid tails oriented along the normal direction, and the kinetic energy of the beads is low. Therefore, simulating at reduced temperature 𝑇 𝐷𝑃𝐷 = 1 may be unsuitable. Kranenburg et al simulated gel phase using DPD simulations, setting 𝑇 𝐷𝑃𝐷 as low as 0.3 and obtained an L𝛼 -𝐿 𝛽 transition in the range 𝑇 𝐷𝑃𝐷 =0.3-0.65, depending on the parameter set. Howerver, Equation ( ) is quite simplistic and does not account for the effects of coarse-graining or the use of smooth and soft interaction potentials. It scales too fast and cannot be used as a proper temperature scaling for CG systems. For example, with 𝑇 𝑟𝑒𝑓 = 300 𝐾, 𝑇 𝐷𝑃𝐷 = 0.6 indicates 𝑇 𝑟𝑒𝑎𝑙 = 180𝐾, which is far from a relevant temperature range for colloidal systems. If one sets 𝑇 𝐷𝑃𝐷 to the corresponding real temperature, then 𝑇 𝑟𝑒𝑓 becomes unrealistically high to be parametrized at; 𝑇 𝑟𝑒𝑎𝑙 = 293 𝐾 implies 𝑇 𝑟𝑒𝑓 = 𝑇 𝑟𝑒𝑎𝑙 𝑇 𝐷𝑃𝐷 = 488𝐾, at which practically no parametrization can be performed. Thus, to reproduce the temperature-dependent phase behavior of phospholipid films, we adopt the following strategy: