  In this paper, we perform molecular dynamics (MD) simulations to study the
random packing of spheres with different particle size distributions. In
particular, we deal with non-Gaussian distributions by means of the L\'evy
distributions. The initial positions, as well as the radii of five thousand
non-overlapping particles, are assigned inside a confining rectangular box.
After that, the system is allowed to settle under gravity towards the bottom of
the box. Both the translational and rotational movements of each particle are
considered in the simulations. In order to deal with interacting particles, we
take into account both the contact and long-range cohesive forces. The normal
viscoelastic force is calculated according to the nonlinear Hertz model,
whereas the tangential force is calculated through an accurate nonlinear-spring
model. Assuming a molecular approach, we account for the long-range cohesive
forces using a Lennard-Jones(LJ)-like potential. The packing processes are
studied assuming different long-range interaction strengths.
