We propose a method for quantifying charge-driven instabilities in clusters, based on equilibrium simulations under confinement at constant external pressure. This approach makes no assumptions about the mode of decay and allows different clusters to
be compared on an equal footing. A comprehensive survey of stability in model clusters of 309 Lennard-Jones particles augmented with Coulomb interactions is presented. We proceed to examine dynamic signatures of instability, finding that rate constants for ejection of charged particles increase smoothly as a function of total charge with no sudden changes. For clusters where many particles carry charge, ejection of individual charges competes with a fission process that leads to more symmetric division of the cluster into large fragments. The rate constants for fission depend much more sensitively on total charge than those for ejection of individual particles.
When a cluster or nanodroplet bears charge, its structure and thermodynamics are altered and, if the charge exceeds a certain limit, the system becomes unstable with respect to fragmentation. Some of the key results in this area were derived by Rayle
igh in the nineteenth century using a continuum model of liquid droplets. Here we revisit the topic using a simple particle-based description, presenting a systematic case study of how charge affects the physical properties of a Lennard-Jones cluster composed of 309 particles. We find that the ability of the cluster to sustain charge depends on the number of particles over which the charge is distributed---a parameter not included in Rayleighs analysis. Furthermore, the cluster may fragment before the charge is strong enough to drive all charged particles to the surface. The charged particles in stable clusters are therefore likely to reside in the clusters interior even without considering solvation effects.
The influence of quadrupolar interactions on the structure of small clusters is investigated by adding a point quadrupole of variable strength to the Lennard-Jones potential. Competition arises between sheet-like arrangements of the particles, favour
ed by the quadrupoles, and compact structures, favoured by the isotropic Lennard-Jones attraction. Putative global potential energy minima are obtained for clusters of up to 25 particles using the basin-hopping algorithm. A number of structural motifs and growth sequences emerge, including star-like structures, tubes, shells and sheets. The results are discussed in the context of colloidal self-assembly.
