During the minimization and solvation, we maintain position restraints with force constant 1000 kJ/nm 2 on the head group and tail ends. After solvation, we slowly release the restraints, decreasing the force constant of restraint to 500, 250, 100, 50, 10, and 0, while minimizing and equilibrating the system for 1 ns at each step. Once the restraints are fully removed we run a long (up to 1 µs) NPT equilibration to see if the system stays in the gyroid phase for an extended period of time. Figure shows the gyroid structure obtained after the 1 µs of equilibration.