For each system, 50'000 steps of steepest descent minimization were performed with a step size of 0.01 nm, stopping if the maximum force dropped below 1000 kJ mol -1 nm -1 . Then, the system temperature and pressure were equilibrated using 100 ps of NVT simulation followed by 1 ns of NpT simulation.