$ if $r < 2^{\frac{1}{6}} \sigma$, otherwise it is zero. The potential has been truncated and shifted. The bonded interactions between the beads, both inter and intra interactions, are finitely extensible nonlinear elastic (FENE) potential to keep the equilibrium distance between the beads (monomers) around $r_{0}= \sigma$, $V \left( r \right)= -\frac{1}{2} K_{\mathrm F} \Delta r_{\mathrm {max}}^{2} \ln \left[ 1-\left( \left(r-r_{0}\right)/\Delta r_{\mathrm {max}} \right)^{2} \right]$, where $\Delta r_{\mathrm {max}}$ is the maximal bond stretching and $K_{\mathrm F}$ is the spring constant. Besides that, to make rod-like stiff segments, we used the cosine squared bond angle potential between united monomers in each segment which increases the intra-stiffness, the bending potential, $  V \left( \phi \right)= K_{\mathrm a}/2 \left[\cos\left(\phi\right)-\cos\left(\phi_{0}\right) \right]^{2} $, $K_{\mathrm a}$ is the bond angle bending constant and $\phi_{0}=\pi$ is the equilibrium bond angle.

The total number of monomers per each polymer chain is $N_{m}=1365$ where we have $4$ chains in each simulation. Each chain is composed of serial connection of similar segments (united monomers). By changing the number of monomers per segment $n_{m}$, Fig. \[fig:seg\_poly\](a), and changing the intra-stiffness (persistence length) of the segments, we create flexible, normal-stiff short, extra-stiff short, and extra-stiff long segments where the flexible, short, 