However, it has also been shown that the expansion series do not necessarily always converge rapidly enough to warrant a cuto↵ after the first non-additive term . The contribution of many-body dispersion interactions were included in the benchmark ↵-phase crystal geometry via the MBD@rsSCS dispersion correction, which includes atomic many-body e↵ects up to infinite order 63,64 , yet it still exhibits the aforementioned increased density along the dispersion-governed a-axis. Thus, it cannot be guaranteed that including atomic manybody e↵ects via. e.g. an additional ATM term in our model would necessarily decrease the resulting liquid density. In any case, the correct form of the short-range damping function is unknown for these types of dispersion corrections. Therefore, another possible, and perhaps more direct way to improve the liquid density in the future could also simply be through further refining the parametrization of our potential based on benchmarks from higher levels of theory that reproduces electronic correlation with greater accuracy and generality, without the need for dispersion .