We benchmarked our simulations against the published results of ref. [1]}.
All simulations were perfomed using the LAMMPS software package [2]}. Particles interact via a 12-6 Lennard-Jones potential (LJ) \(V_{LJ}(r)=4\epsilon ^{LJ}_{ij}[(\sigma ^{LJ}_{ij}/r)^{12}-(\sigma ^{LJ}_{ij}/r)^6]\)  shifted and truncated at \(r= r_{\text{cut}}\) , such that
\(V_{TS}(r)={\left\lbrace \begin{array}{ll}V_{LJ}(r)-V_{LJ}(r_{\text{cut}}), & \text{if}\ r \le r_{\text{cut}} \\0, & \text{otherwise}.\end{array}\right.}\) 
