Self-assembled clusters of mutually repelling particles in confinement

Mutually repelling particles form spontaneously ordered clusters when forced into confinement. The clusters may adopt similar spatial arrangements even if the underlying particle interactions are contrastingly different. Here we demonstrate with both simulations and experiments that it is possible to induce particles of very different types to self-assemble into the same ordered geometric structure by simply regulating the ratio between the repulsion and confining forces. This is the case for both long- and short-ranged potentials. This property is initially explored in systems with two-dimensional (2D) circular symmetry and subsequently demonstrated to be valid throughout the transition to one-dimensional (1D) structures through continuous elliptical deformations of the confining field. We argue that this feature can be utilized to manipulate the spatial structure of confined particles, thereby paving the way for the design of clusters with specific functionalities.

💡 Research Summary

**
The paper introduces a unifying framework for the self‑assembly of mutually repelling particles confined within a finite region. By defining a single dimensionless parameter, the repulsion‑to‑confinement ratio λ, the authors demonstrate that clusters of vastly different physical nature—charged particles interacting via Coulomb forces, magnetic dipoles, hard spheres, and soap bubbles—can be driven to adopt identical equilibrium geometries. The interaction between any two particles is modeled by a generalized pair potential

 u(r) = u₀ (r/R)^{‑α},

where R is the characteristic size of the confining region, u₀ sets the energy scale, and the exponent α encodes the range of the interaction (α = 1 for Coulomb, α ≈ 3 for magnetic dipoles, α → ∞ for short‑range contact forces).

Confinement is represented by M auxiliary particles placed uniformly along the boundary of a circle (or ellipse) of radius R. These boundary particles carry a reduced “charge” λ < 1 relative to the N particles that form the cluster. The total potential energy of a single cluster particle is the sum of the repulsive contribution from the other N − 1 particles (U_R) and the attractive (confining) contribution from the M boundary particles (U_C). In dimensionless form the equilibrium condition reads

 ∂U_C/∂ρ = –λ ∂U_R/∂ρ, with ρ = r/R.

Because ∂U_R/∂ρ diverges as ρ → 0 and decays monotonically toward the boundary, while ∂U_C/∂ρ behaves oppositely, there always exists a solution 0 ≤ ρ ≤ 1 for any positive λ and α. The key insight is that changing α (i.e., changing the physical interaction) can be compensated by an appropriate adjustment of λ, leaving the equilibrium radius ρ—and therefore the spatial arrangement of the particles—unchanged.

The authors first explore the circular case (ε = 0, where ε denotes the eccentricity of an ellipse). They show analytically and numerically that for a given α there is a one‑to‑one mapping λ(α) that preserves the same shell radius. Figure 1 illustrates this balance for Coulomb (α = 1, λ = 1) and for a steeper potential (α = 3, λ ≈ 2.2). The same principle extends to short‑range interactions: for hard spheres the authors impose a maximal packing fraction criterion to determine λ, while for soap bubbles they minimize the variance of Voronoi cell areas, which effectively enforces equal “pressure” among neighboring bubbles. Both procedures recover the experimentally observed configurations.

The study then generalizes to elliptical confinement, where a second geometric parameter ε is introduced. In this case a single λ no longer suffices; instead a surface in the (λ, α, ε) space defines the set of parameters that yield an identical equilibrium configuration. Figure 2 demonstrates that a cluster of five particles confined in ellipses of different eccentricities can be made to occupy the same positions when the interaction exponent α is changed from 1 (Coulomb) to 3 (magnetic dipole) provided λ and ε are co‑adjusted accordingly.

Experimental validation is performed on three model systems: (i) monodisperse hard spheres, (ii) floating permanent‑magnet dipoles, and (iii) monodisperse soap bubbles. For each system the authors prepare clusters of N = 5 and N = 10 particles inside ellipses with ε = 0, 0.44, 0.66, 0.87, 0.95. Photographs (Table I) reveal that despite the completely different interaction laws, the clusters adopt the same topologies—ring‑like arrangements for low ε and linear chains as ε → 1. The experimental configurations match the minima obtained from the energy‑minimization simulations, confirming the universality of the λ‑controlled design rule.

Beyond reproducing known “classical Wigner clusters,” the work expands the concept to a broad class of potentials and confinement geometries. It shows that the competition between a monotonically decreasing repulsive force and a monotonically increasing confining force is the essential ingredient for ordered self‑assembly, irrespective of the microscopic origin of the forces.

The authors discuss potential applications. By tuning λ (through, e.g., the strength of an external electric or magnetic field, the number density of boundary charges, or the surface tension in foams) and the shape of the confinement, one can deliberately engineer clusters with prescribed symmetry and connectivity. This could be exploited in the fabrication of metamaterials with tailored electromagnetic response, in microfluidic devices where particle positioning dictates flow patterns, or in colloidal crystal engineering where defect‑free arrangements are required.

In summary, the paper provides a concise yet powerful theoretical model, corroborated by extensive simulations and three distinct experimental platforms, that establishes a universal design principle for confined repulsive particle clusters. The repulsion‑to‑confinement ratio λ, together with the eccentricity ε of the confining field, constitutes a minimal set of control parameters capable of dictating the equilibrium geometry across a wide spectrum of interaction ranges. This insight opens new avenues for the purposeful design of functional particle assemblies in both hard‑matter and soft‑matter contexts.