where 0 is the center of the bilayer and Dz is a point in the water region where P(W)z = 1. In our calculations, DB is integrated over both leaflets simultaneously to avoid defining a bilayer center. The water density was smoothened using a convolution window of 5 ns (5 frames) before integration.