The number of atoms, T max , r c , and t m were variables that depended on the target system used (described in each section). All ensembles for these steps were NVT. An equilibration run at 300 K was carried out for 100 ps under NVT, and a production run was subsequently performed for 100 ps under the isothermal-isobaric ensemble (NPT) with a time step of 1 fs at 300 K and 1 atm. Three simulation runs on each system were independently performed using three different initial structures. This was to evaluate the deviation of each simulation run from the final glass structure. The standard errors in structural data such as density and bond distances were calculated from three independent runs.