The minimized bilayer was then duplicated in z using GROMACS and minimized again using the same approach as before. The stacked bilayer systems were then relaxed The Berendsen 187 and v-rescale 198 thermostat were used during relaxation and production run, respectively. During both steps τT = 1 ps and lipids were coupled separately from the solvent.