The equilibrium areas per lipid at which the 𝛾 𝐵𝐿 = 0 are found to be ∼ 52 𝐴 𝑜 2 at 293 K, ∼ 55 𝐴 𝑜 2 at 304 K and 62 𝐴 𝑜 2 at 324K. The corresponding experimental values vary within a certain range. In the gel phase, values reported by different groups include 47.9 Å 2 (293K) , 50Å 2 (294 K) , 48.6 Å 2 (298K) , 48.6Å 2 (300K) and 52.3 Å 2 (298K). Similarly, the range of values reported for liquid crystalline phase at 323K spans 57 -71 Å 2 . Our value of 62 Å 2 at 324 K is in good agreement with the value 62.9 Å 2 reported by Nagle et al and some of the atomistic and coarse grained force fields. Partial number density profiles 𝜌 𝑖 (𝑧) along the normal z-direction, of different bead types in bilayer systems of equilibrium area per lipid at various temperatures are given in Figure . The bilayer thickness, estimated as the peak-to peak distance of the choline groups, 𝑡 𝑁𝑁 = 4.73 nm in the gel phase, is in agreement with the experiments. At 324 K, the thickness of the equilibrium liquid crystalline bilayer is 4.1 nm, also close to experimental result, while the thickness of the bilayer of 𝑎 𝐿 = 69.9 𝐴 𝑜 2 is 3.8 nm. In conclusion, our DPPC model reproduces most of the experimentally observed characteristics, with parameters transferable over the temperature range that involves the main phase transition behavior of the lipid.