As seen in Fig. , GLST achieves O (N) scaling as a result of its cubature and domain decomposition scheme. The RMSE of the energy and 3 components of the force all meet the target error threshold of 10 -8 (Fig. ). The slight increase in RMSE with system size is a consequence of the cubature being used, as more energies and forces are being approximated by several different cubatures as described in Sec. III B, each with different errors on the order of the target error threshold value. Despite the errors from multiple cubatures, the RMSE is still on the order of 10 -8 kcal/mol even for the largest system used containing 1.5 million atoms. We also computed histograms of the per-atom error for representative systems to understand our distribution of errors (Figs. and). While maximum errors for the energies and forces slightly exceed the target error threshold, they are clustered around 0 resulting in a desirable RMSE.