\epsilon_s \left[ \left(\frac{\sigma_s}{r_{s}}
\right)^{\nu}-\left( \frac{\sigma_s}{h_{0}}
\right)^{\nu}\left(1-\frac{(r_s-h_0)\nu}{h_0} \right) \right] & 0<r_{s}\leq h_0\\
0 & otherwise\,,\\
\end{cases}
\end{split}\end{aligned}$$

where $\epsilon_s=1.5$ is the energy scale, $\nu=2$, $r_s \equiv r_{ij}-\sigma$ is the inter-particle distance (with $r_s > 0$ to ensure that particles never overlap), $\sigma_s=0.5$ is the characteristic separation length between particles, and $h_0$ the distance at which the potential is shifted to ensure its continuity at $r_s=h_0$. [^1]

The attractive potential consists of a radial $V_{att}(r_{ij})$ and an angular $\phi(\theta_i,\theta_j$) contribution. Concerning the radial part of the potential, we have considered two different truncated Lennard-Jones: one long ($r_c =2.5 \sigma$) and another short ($r_c =1.5 \sigma$) range. While the following equation: $$\begin{aligned}
\label{eq:VLJmidrange}
\begin{split}
&V_{att} (r_{ij}) = \\
&\begin{cases}
4 \epsilon \left[ \left( \frac{\sigma}{r_{ij}} \right)^{12} - \left( \frac{\sigma}{r_{ij}} \right)^6 \right] & (2^{1/6}\sigma)<r_{ij}<r_{c}\,,\\
-\epsilon & r_{ij}< (2^{1/6}\sigma)\,, \\
0 & otherwise\,, \
\end{cases}
\end{split}\end{aligned}$$

represents the former, it has been modified for the latter by adding a smooth transition function, so that the potential $V_{att}$ and its derivative vanish continuously as $r_{ij}$ approaches $r_c$ [@pra_smoothfunct] (see appendix \[Appendix:attra\]). Therefore, whenever $(2.2^{1/6}\sigma)<r_{ij}<r_{c}$

$$\begin{aligned}
\label{eq:VLJshortrange}
V_{att} (r_{ij}) = 4\epsilon \left[ \left( \frac{\sigma}{r_{ij}}
\right)^{12} - \left( \frac{\sigma}{r_{ij}} \right)^{6} \right] 
 + 4\epsilon \left[ 6 \left( \frac{\sigma}{r_{c}}
\right)^{12} - 