We investigate the applicability of these classical FF adaptations by computing finite temperature bulk moduli, with the same methods we used on ReaxFF, starting with the finite strain difference method (illustrated on Figure ) to allow direct comparison with AIMD results. The parameters used for each method are detailed in the supporting information. We note, however, that there is an important variance in the bulk moduli of the glasses, with potentially divergent values with respect to the FF adaptation for the same glass model, as shown on Figure . It highlights the importance of using a large enough number of glass configurations to obtain average mechanical properties for models with a small system size.