Finally, after the melt-quenching procedure, the glass model obtained was then equilibrated for 20 ps in the same (N, P, T ) ensemble as during melt-quenching, the duration of one constant temperature plateau. The simulation was then further run for production during 1 ns in the (N, V, T ) ensemble to gather statistics of structural and dynamical properties of the glass. We checked that up to 10 ns of further equilibration in the (N, P, T ) ensemble, even with a flexible cell, does not perceptibly change the measured properties of the glass.