  We report a numerical analysis of the unforced break-up of free cylindrical
threads of viscous Newtonian liquid whose interface is coated with insoluble
surfactants, focusing on the formation of satellite droplets. The initial
conditions are harmonic disturbances of the cylindrical shape with a small
amplitude $\epsilon$, and whose wavelength is the most unstable one deduced
from linear stability theory. We demonstrate that, in the limit $\epsilon \to
0$, the problem depends on two dimensionless parameters, namely the Laplace
number, $La=\rho \sigma_0 \bar{R}/\mu^2$, and the elasticity parameter,
$\beta=E/\sigma_0$, where $\rho$, $\mu$ and $\sigma_0$ are the liquid density,
viscosity and initial surface tension, respectively, $E$ is the Gibbs
elasticity and $\bar{R}$ is the unperturbed thread radius. A parametric study
is presented to quantify the influence of $La$ and $\beta$ on two key
quantities: the satellite droplet volume and the mass of surfactant trapped at
the satellite's surface just prior to pinch-off, $V_{sat}$ and $\Sigma_{sat}$,
respectively. We identify a weak-elasticity regime, $\beta \lesssim 0.05$, in
which the satellite volume and the associated mass of surfactant obey the
scaling law $V_{sat} = \Sigma_{sat} = 0.0042 La^{1.64}$ for $La \lesssim 2$.
For $La \gtrsim 10$, $V_{sat}$ and $\Sigma_{sat}$ reach a plateau of about $3
\%$ and $2.9 \%$ respectively, $V_{sat}$ being in close agreement with previous
experiments of low-viscosity threads with clean interfaces. For $La<7.5$, we
reveal the existence of a discontinuous transition at a critical elasticity
$\beta_c (La)$, with $\beta_c \to 0.98$ for $La \lesssim 0.2$, such that
$V_{sat}$ and $\Sigma_{sat}$ abruptly increase. The jumps experienced by both
quantities reach a plateau when $La \lesssim 0.2$, while they decrease
monotonically as $La$ increases up to $La = 7.5$, where both become zero.
