Spherical Shell Lattice Models

• Spherical shell lattice models are systems where particles, spins, or lattice points are constrained to spherical surfaces, emphasizing finite geometry and inherent topological defects.
• They utilize varied interaction laws—including Lennard-Jones, Hertzian, and Coulomb potentials—to explore energetic landscapes and the resulting order or defect structures.
• Dynamic construction methods and advanced simulation protocols reveal robust kinetic growth, yielding high-symmetry shells analogous to viral capsid and spherical design formations.

Searching arXiv for the cited paper and closely related spherical-shell lattice work. “Spherical shell lattice model” denotes a class of constructions in which discrete degrees of freedom are constrained to a spherical manifold and their organization is determined by curvature, topology, and the interaction law. In the recent literature, the term covers nonequilibrium assembly of identical Lennard-Jones particles on the unit sphere (Golushko et al., 12 Mar 2026), packed soft-core particles on a rigid spherical shell with competing hexagonal and square order (Xie et al., 2 Jan 2025), XY spins on a spherical Fibonacci lattice (Song et al., 2021), ultra-soft cluster crystals on $S^2$ (Franzini et al., 2018), Coulomb-coupled point particles on a sphere in a Thomson-type setting (Bachmann et al., 2024), and normalized shells of the $D_4$ lattice on $S^3$ studied as spherical designs (Hirao et al., 2023). Across these usages, the model family is characterized by a finite spherical geometry, unavoidable topological constraints, and a shell-based notion of order.

1. Terminological scope and geometric setting

In the particle-assembly literature, the shell is a finite cluster whose constituents are constrained to the surface of a sphere. Golushko et al. consider $N$ identical point particles on the surface of a sphere of unit radius $R=1$, with shell growth proceeding sequentially until no new particle can be attached with negative binding energy; the final shell size is denoted $n^\*$, and the energy per particle is $u=E_{\mathrm{tot}}/n^\*$ (Golushko et al., 12 Mar 2026). Xie et al. likewise place particles on a spherical surface $S^2(R)$, but study competition between hexagonal and square patterns under a Hertzian interaction (Xie et al., 2 Jan 2025). In the ultra-soft cluster-crystal setting, the sphere is again the substrate, with distances measured along the great circle and the shell viewed as a curved two-dimensional assembly (Franzini et al., 2018).

In lattice-spin work, the shell may be a discrete spherical sampling rather than a self-assembled particle cluster. The spherical Fibonacci lattice is constructed by placing $N$ points on a sphere of radius $R$ through

$D_4$0

with $D_4$1, followed by the corresponding Cartesian embedding; nearest neighbors are then defined by a cutoff radius $D_4$2 (Song et al., 2021). In the $D_4$3-lattice setting, a shell is instead the set of lattice points of fixed squared norm, normalized onto the unit three-sphere $D_4$4; this is a shell in the arithmetic sense rather than a physical particle shell (Hirao et al., 2023).

This range of usage suggests that “spherical shell lattice model” is not a single standardized formalism. Rather, it names a family of models sharing the constraints of spherical embedding and finite topology, while differing in whether the microscopic objects are particles, spins, or lattice points.

2. Interaction laws and energy functionals

The simplest shell-assembly realization uses an isotropic Lennard-Jones pair potential depending on the geodesic chord distance $D_4$5,

$D_4$6

with equilibrium separation $D_4$7 at which $D_4$8. A weak second minimum can be introduced through the Lennard-Jones-Gauss potential to mimic anisotropy (Golushko et al., 12 Mar 2026). The sole control parameter in the minimal model is $D_4$9, with $S^3$0.

The competing Hex–Sq model employs a purely repulsive Hertzian contact potential of range $S^3$1 and energy scale $S^3$2,

$S^3$3

with the total energy $S^3$4. The spherical metric is

$S^3$5

and distances entering the interaction are computed from the chord distance $S^3$6 (Xie et al., 2 Jan 2025).

The XY model on a spherical Fibonacci lattice replaces particle coordinates by tangent-plane spins and introduces a Gaussian-screened interaction on the neighbor graph,

$S^3$7

with $S^3$8 in the reported calculations (Song et al., 2021). The ultra-soft cluster-crystal model uses the GEM-4 interaction

$S^3$9

written on the sphere through the geodesic distance $N$0 or equivalently $N$1 (Franzini et al., 2018).

The Coulomb spherical-lattice model takes

$N$2

supplemented during equilibration by the viscous drag force $N$3, with all forces projected to the tangent plane of the sphere (Bachmann et al., 2024). In the elastic-shell many-body model, the shell itself is deformable: the discrete per-patch elastic energy is

$N$4

and many-body overlap is handled by the rule $N$5 (Boattini et al., 2020).

3. Dynamical construction and computational protocols

The nonequilibrium Lennard-Jones shell grows in discrete steps. At each step, approximately $N$6 approximately equidistant trial points are placed in a spherical cap of angular radius approximately $N$7 around each already attached particle. For a trial point $N$8, the binding energy of the would-be $N$9th particle is

$R=1$0

The particle is attached at the trial point for which $R=1$1 is most negative, and the entire cluster then undergoes constrained energy minimization on the sphere. Growth stops when no trial point yields a negative binding energy (Golushko et al., 12 Mar 2026).

The same framework admits a finite-temperature generalization. If the local minima of the attachment energy are $R=1$2, then one samples attachment site $R=1$3 with probability

$R=1$4

The limit $R=1$5 reproduces deterministic growth, while $R=1$6 yields a random choice among local minima (Golushko et al., 12 Mar 2026).

Other spherical shell lattice models use standard large-scale simulation protocols. The Hex–Sq study employs overdamped Langevin dynamics in LAMMPS with dimensionless time step $R=1$7, friction $R=1$8, $R=1$9 independent simulated-annealing runs per $n^\*$0 point, and cooling from $n^\*$1 down to $n^\*$2 (Xie et al., 2 Jan 2025). The Fibonacci-sphere XY study uses classical Monte Carlo annealing with local Metropolis updates, with $n^\*$3 MC steps for stable vortex-pattern snapshots and a one-layer graph-convolutional-network classifier trained on low- and high-temperature configurations (Song et al., 2021). The GEM-4 cluster-crystal model combines classical density functional theory, based on minimization of a grand-potential functional on $n^\*$4, with Monte Carlo simulations at fixed $n^\*$5, area $n^\*$6, and temperature (Franzini et al., 2018). The Coulomb “Crystal Ball” framework integrates Newton’s law with adaptive time step, exact reprojection onto the sphere, and a phase-transition protocol in which $n^\*$7 particles are removed around a random seed after equilibration (Bachmann et al., 2024). The elastic-shell many-body model accelerates simulation by fitting the costly shell-deformation energy to $n^\*$8 symmetry functions, reducing evaluation time by $n^\*$9–$u=E_{\mathrm{tot}}/n^\*$0 before Monte Carlo phase-diagram calculations in the two-dimensional $u=E_{\mathrm{tot}}/n^\*$1 ensemble (Boattini et al., 2020).

4. Topology, defects, and order metrics

A defining feature of spherical shell lattice models is that order is topologically frustrated. In the Lennard-Jones shells, most vertices are $u=E_{\mathrm{tot}}/n^\*$2-coordinated, but Euler’s theorem enforces exactly twelve $u=E_{\mathrm{tot}}/n^\*$3-coordinated defects in hexagonal order. When square tiles occur in square–triangular shells, vertices with coordination $u=E_{\mathrm{tot}}/n^\*$4 or $u=E_{\mathrm{tot}}/n^\*$5 appear, corresponding to $u=E_{\mathrm{tot}}/n^\*$6 defect angles (Golushko et al., 12 Mar 2026). In the Hertzian model, the Gauss–Bonnet theorem gives the net disclination charge on $u=E_{\mathrm{tot}}/n^\*$7 as $u=E_{\mathrm{tot}}/n^\*$8, implying twelve $u=E_{\mathrm{tot}}/n^\*$9 disclinations in a purely hexagonal lattice or eight $S^2(R)$0 disclinations in a purely square lattice (Xie et al., 2 Jan 2025).

Xie et al. formulate this competition through bond-orientational order parameters

$S^2(R)$1

with $S^2(R)$2 for Hex and $S^2(R)$3 for Sq. They also introduce the incompatibility angle

$S^2(R)$4

which measures the mismatch between the interior angle of a spherical polygon and its Euclidean counterpart, and use it as a geometric criterion for the nucleation of counter-domains inside larger domains (Xie et al., 2 Jan 2025).

In the XY model, the analogous topological objects are vortices rather than disclinations. The vortex winding number is

$S^2(R)$5

and the Poincaré–Hopf theorem enforces total topological charge $S^2(R)$6 for a continuous tangent-vector field on $S^2(R)$7, forcing at least two vortices in the ground state (Song et al., 2021). In the ultra-soft cluster-crystal model, topological frustration appears at the level of cluster occupancy and local stability: fivefold disclination clusters have smaller occupancy $S^2(R)$8, smaller nearest-neighbor spacing $S^2(R)$9, and higher local grand potential $N$0 than non-defective clusters (Franzini et al., 2018).

The arithmetic shell model of $N$1 encodes an algebraic analogue of spherical order. The normalized shell

$N$2

is an antipodal spherical $N$3-design for every $N$4, and the root shell $N$5 is a tight $N$6-design of size $N$7 (Hirao et al., 2023). Here the “defect” language is replaced by harmonic vanishing conditions and linear-programming bounds, but the shell is still organized by spherical symmetry constraints.

5. Nonequilibrium Lennard-Jones shell assembly

Scanning the single parameter $N$8 with $N$9, in steps down to $R$0, yields $R$1 distinct shells with $R$2 under the deterministic sequential protocol. These shells fall into three classes: icosahedral shells, intermediate-symmetry shells with a single $R$3-fold axis $R$4, and low-symmetry shells with no axis $R$5 (Golushko et al., 12 Mar 2026).

The icosahedral class contains $R$6 and $R$7 shells with symmetry $R$8, corresponding to the classic $R$9 and $D_4$00 capsids. A deformed $D_4$01 shell at $D_4$02 appears with $D_4$03 symmetry and, over a small $D_4$04 range, with $D_4$05 rather than perfect $D_4$06. The classic $D_4$07 shell at $D_4$08 and the perfect $D_4$09 icosahedron at $D_4$10 do not assemble in the identical-LJ model (Golushko et al., 12 Mar 2026).

A major result is the spontaneous formation of square–triangular surface order. The reported shells include $D_4$11 with symmetry $D_4$12, identified as the snub-cube Archimedean net; $D_4$13 with symmetry $D_4$14, described as a cubic S–T net with vacancies on the four-fold axes; $D_4$15 with symmetry $D_4$16; $D_4$17 with symmetry $D_4$18; and $D_4$19 with symmetry $D_4$20 (Golushko et al., 12 Mar 2026). These spherical polyhedral graphs correspond closely to experimental protein nanocages, including alpha-tocopherol transfer protein cages, ferritin mutants, and allophycocyanin assemblies, after superposition on PDB coordinates.

The energy landscape is organized into $D_4$21-bands. For each shell size one obtains a band in $D_4$22; between bands, no shell of that size is stable, and band overlap can produce multiple isomers with the same $D_4$23 but different symmetries. Deterministic growth yields energy-minimized shells for $D_4$24. In stochastic assembly, $D_4$25 runs at $D_4$26, $D_4$27, and $D_4$28 show more than $D_4$29 perfect yield for high-symmetry shells such as $D_4$30 even at $D_4$31; most intermediate-symmetry shells with $D_4$32 remain above $D_4$33 yield for $D_4$34; and low-symmetry shells with $D_4$35 drop in yield but remain at least $D_4$36 even at $D_4$37 (Golushko et al., 12 Mar 2026).

6. Other spherical-shell lattice realizations

The Hertzian Hex–Sq model identifies six distinct regimes as the reduced density $D_4$38 is varied: pure Hex with point disclinations or scars, Hex background with Sq domains, nested domain/counter-domain patterns, interwoven bands of both tilings, Sq background with Hex domains, and pure Sq with eight $D_4$39 defects. The rich defect morphologies occur in the narrow window $D_4$40, where Hex and Sq free energies are nearly equal, and the excess disclination number follows the scaling $D_4$41 with $D_4$42 (Xie et al., 2 Jan 2025).

On the spherical Fibonacci lattice, the near-uniformity of the graph discretization is central. For $D_4$43 and $D_4$44, the lattice has $D_4$45 neighbor pairs; $D_4$46 sites have $D_4$47 neighbors, $D_4$48 have $D_4$49, and $D_4$50 have $D_4$51. At $D_4$52, $D_4$53, and $D_4$54, two $D_4$55 vortices remain after full annealing and sit approximately $D_4$56 apart; for $D_4$57, they become nearly antipodal at $D_4$58. The graph-convolutional-network predicts $D_4$59 for $D_4$60, $D_4$61 for $D_4$62, and $D_4$63 for $D_4$64 (Song et al., 2021).

The GEM-4 cluster-crystal model generalizes the bulk clustering criterion to the sphere by expanding the potential in Legendre modes and identifying a $D_4$65-line through

$D_4$66

At fixed $D_4$67 and $D_4$68, one finds $D_4$69–$D_4$70, above which a $D_4$71-cluster crystal is stable. At fixed $D_4$72, varying $D_4$73 yields cluster phases with $D_4$74 clusters. The occupancy per site grows linearly with $D_4$75, while $D_4$76 remains nearly constant, which is the reported hallmark of cluster crystals (Franzini et al., 2018).

The Coulomb spherical-lattice model studies rearrangement rather than static pattern formation. For $D_4$77 and $D_4$78, with $D_4$79 Monte Carlo trials per pair, the analytic fit

$D_4$80

gives $D_4$81 at $D_4$82 confidence and reduced $D_4$83. At $D_4$84, the reported average peak kinetic energy rises from approximately $D_4$85 eV for $D_4$86 to approximately $D_4$87 eV for $D_4$88 (Bachmann et al., 2024).

The elastic-shell many-body model yields a simpler phase diagram. For all $D_4$89, only two stable phases are reported: a low-density disordered fluid with $D_4$90 and a high-density hexagonal crystal with $D_4$91. At $D_4$92, coexistence occurs at $D_4$93–$D_4$94, corresponding to $D_4$95–$D_4$96 (Boattini et al., 2020).

7. Significance, limitations, and broader implications

The nonequilibrium Lennard-Jones study shows that a minimal model with only one geometric parameter, $D_4$97, reproduces both small icosahedral viral capsids and a broad set of octahedral, tetrahedral, and square–triangular shells. It further shows that nonequilibrium sequential attachment into the local deepest binding site often leads directly to global or near-global energy minima, especially for high symmetry, and that finite-temperature stochastic growth remains robust for the most symmetric shells (Golushko et al., 12 Mar 2026). This suggests that kinetic growth rules alone can be structurally selective on curved manifolds.

The competing-lattice and cluster-crystal studies emphasize a complementary principle: curvature does not merely perturb flat-space order, but changes the admissible defect content, domain morphology, and stability criteria. In the Hex–Sq system, the incompatibility angle provides a simple geometric predictor for when a domain must nucleate a counter-domain (Xie et al., 2 Jan 2025). In the ultra-soft cluster-crystal system, the spherical harmonic decomposition of the interaction identifies the $D_4$98-line and links preferred cluster number to the first negative harmonic mode $D_4$99 (Franzini et al., 2018).

The arithmetic shell perspective extends the topic beyond soft matter. The $S^3$00 results show that shells normalized onto $S^3$01 can be classified through spherical design theory, linear-programming bounds, and modular forms, with the root shell uniquely realizing the antipodal spherical $S^3$02-design of size $S^3$03 (Hirao et al., 2023). A plausible implication is that “shell lattice” can refer either to self-assembled physical shells or to highly symmetric norm shells of a Euclidean lattice.

Limitations are model-specific and explicit in the cited work. The Coulomb model is purely classical and ignores quantum tunneling, electron screening, magnetic fields, and relativistic effects (Bachmann et al., 2024). The elastic-shell many-body model replaces a direct shell-deformation calculation by a symmetry-function regression surrogate (Boattini et al., 2020). The Lennard-Jones assembly model uses identical particles and isotropic pair interactions unless the Lennard-Jones-Gauss extension is invoked (Golushko et al., 12 Mar 2026). Even with these simplifications, the models collectively provide a blueprint for rational shell design in protein nanocages, colloidosomes, micelles, spherical spin systems, and mathematically defined spherical designs.