Toward Quantitative Predictions in Ffellonics: What Contact-Counting Can and Can't Yet Tell Us

Ffellonics is currently a qualitative framework. It describes a clear hierarchy that begins with the first symmetric touch and proceeds toward higher coordination, driven by symmetric attachment and free-energy minimization. That narrative is geometrically and thermodynamically well-grounded, but it doesn't yet assign numbers, probabilities, or rates to anything. Thermodynamics offers tools to start closing that gap — treating each configuration as a point on a free-energy landscape and each transition as a rate process — but it's worth being precise about how far those tools actually get us today, because the honest answer is: further than nothing, not yet all the way to a validated prediction.

Theoretical Foundation

Consider a system choosing between two configurations of the same N spheres — a higher-coordination configuration and a lower-coordination alternative. For a simple contact potential (energy = −ε per contact, a standard simplified model used in colloidal cluster physics), the energies differ by:

ΔU ≡ U(high-coordination) − U(alternative) = −ε·ΔN_contacts

where ΔN_contacts is the number of additional contacts the higher-coordination configuration makes. This is an energy difference, not yet a free-energy difference — ΔF = ΔU − TΔS, and the entropy term is not automatically small. Two worked examples below show why: one where the entropy term comes from genuine mechanical floppiness, and one where it comes from symmetry alone.

The rate of transition between two configurations, when an additional sphere arrives and must choose an attachment site, follows an Arrhenius-like expression, k = A·exp(−E_a/kT), where A depends on the attempt frequency (particle size, solvent viscosity, concentration) and E_a is the activation barrier. Predicting absolute timescales requires estimating A, which is highly system-dependent; what follows addresses relative probabilities, not absolute rates.

Neither worked example below is a step in Ffellonics' own named hierarchy, which runs dyad (N=2) → triangle (N=3) → tetrahedron (N=4) → octahedron (N=6) → icosahedron. Both are illustrations of the same generative principle — maximize contacts, minimize free energy — applied to cluster sizes at or adjacent to the hierarchy, offered to show what the method can and can't yet deliver.

Worked Example 1: Five Spheres, and Why the Energy Term Isn't the Whole Story

Configuration A — triangular bipyramid (N=5): three equatorial spheres mutually in contact (3 contacts) plus each equatorial sphere touching both poles (6 contacts). Total: 9 contacts, U = −9ε.

Configuration B — square pyramid (N=5): four base spheres in a ring, each touching two neighbors (4 contacts), plus each base sphere touching the apex (4 contacts). Total: 8 contacts, U = −8ε.

Both are real, buildable equal-edge solids. The bipyramid has one more contact, so ΔU = −ε, and the energy-only Boltzmann ratio is:

P_energy(bipyramid) / P_energy(pyramid) = exp(ε/kT)

At ε = 5kT: ≈148. At ε = 10kT: ≈2.2×10⁴. These are correct evaluations of the energy term alone — and that caveat is doing real work, because two separate entropic effects push against the bipyramid, not one.

The floppy-mode effect. A cluster of N spheres in three dimensions is mechanically rigid, generically, only when it has at least 3N−6 contacts — for N=5, that's 9. The bipyramid has exactly 9: it's isostatic, with no internal freedom to flex. The square pyramid, with only 8, is one short of that threshold: its square base can flex like a four-bar linkage while every existing contact stays intact, with the apex adjusting to hold its own four contacts as the base deforms. This isn't a point competing with the bipyramid on equal footing — it's a positive-dimensional set of configurations, and the sticky-sphere literature (Arkus, Manoharan & Brenner's enumeration; Meng, Arkus, Brenner & Manoharan, Science 327, 560, 2010; Holmes-Cerfon et al., "Sticky-sphere clusters") treats exactly this kind of comparison by integrating the Boltzmann weight over that floppy manifold rather than treating it as a second rigid state. That integral favors the pyramid and hasn't been computed here.

The symmetry-number effect — a separate penalty, same direction. This is a second, independent correction, and it's worth stating explicitly rather than folding into "entropy" as a single unexamined term. A rigid cluster's rotational partition function scales as 1/σ, where σ is the order of its point-symmetry group — a more symmetric structure is suppressed relative to a less symmetric one, because more of its rotations map the cluster onto an indistinguishable copy of itself. The triangular bipyramid has symmetry group D₃h, order 12; the square pyramid has symmetry group C₄v, order 8. The bipyramid is the more symmetric structure, so this term also favors the pyramid — in the same direction as the floppy-mode effect, not offsetting it.

So the honest statement is: the bipyramid is favored by a known, calculable energy term of exp(ε/kT), and at least two distinct, independently-documented entropic effects — floppy-mode configurational entropy and rotational symmetry-number suppression — work against that advantage by amounts this piece does not compute. The reported figures of 148 and 2.2×10⁴ are energy-only ratios, not free-energy ratios, and should be read as upper bounds on how strongly the bipyramid is favored, not as the final answer. (Meng et al. report observing only the bipyramid at N=5 under their experimental depletion conditions, which is qualitatively consistent with the energy term winning — but that observation isn't a measurement of the specific ratio computed here, and shouldn't be cited as if it were.)

Worked Example 2: Six Spheres — A Cleaner Test, Because Both Structures Are Rigid

N=5 has the complication of comparing a rigid structure to a floppy one. N=6 offers a cleaner comparison, because both candidates are genuinely rigid, which removes the floppy-mode unknown and isolates the symmetry-number effect on its own.

At N=6, the maximum-contact (12-contact) packings are the octahedron and the polytetrahedron (three face-sharing tetrahedra). Both are isostatic. Their symmetry groups, however, are starkly different: the octahedron's rotation group has order 24 (it's the same group, S₄, that appears at Ffellonics' own Level 4 milestone), while the polytetrahedron's point group has order 2. Because both structures are rigid, the "entropy is small if both clusters are rigid" shortcut is at least available in principle here — but the symmetry-number ratio alone is already a factor of 12 before any vibrational calculation is done, which is large enough that it cannot be waved away as a minor correction even in the rigid-rigid case. Whether the octahedron's energy advantage (if it has one at equal contact count, the comparison is purely entropic) survives a symmetry-number penalty of that size is exactly the kind of calculation this method should do next — and the sticky-sphere literature above provides the machinery to do it. It has not been carried out here.

General Method, Honestly Scoped

The same approach extends to any transition where competing configurations can be identified and enumerated:

For the same N, enumerate the candidate configuration and its most likely alternative.
Count contacts under the specified potential, and compute ΔU.
Check whether each configuration is isostatic (3N−6 contacts, generically rigid) or under-constrained (fewer contacts, generically floppy). A floppy competitor requires integrating over its soft manifold, not treating it as a second point-like state.
Compute the symmetry-number ratio for both configurations, whether or not either is floppy — this term exists even when both competitors are rigid, as the N=6 case shows, and can be large.
Only once both of these are accounted for does ΔU become ΔF. Determine the activation barrier from the minimum energy path, and combine with detailed balance for the rate ratio.

The transitions between Ffellonics' own named milestones — tetrahedron (6 contacts) to octahedron (12 contacts), for instance — involve adding spheres rather than rearranging a fixed-N cluster, so they require a grand-canonical or kinetic treatment rather than the canonical comparisons used above. That calculation hasn't been done here; the qualitative expectation (an incoming sphere gains stability by attaching in the Ffellonic-preferred position) is reasonable, but it remains qualitative until it's carried out with the same rigor — rigidity-checking and symmetry-number accounting included — as the two examples above.

Significance and Limitations

What this approach can honestly offer:

A real, checkable method for comparing candidate configurations at fixed N, directly comparable to data from colloidal self-assembly experiments, molecular dynamics, and cluster databases such as the Cambridge Cluster Database.
A concrete demonstration that the contact-maximizing configuration is favored by a substantial, calculable energy term at experimentally relevant interaction strengths — with that term now correctly labeled as an energy-only bound rather than a free-energy result.
A clear statement of exactly which calculations remain: the floppy-manifold integral at N=5, the symmetry-number-corrected partition functions at N=6, and a grand-canonical treatment of the actual milestone-to-milestone transitions.

What it can't yet offer: a validated numerical prediction for any named Ffellonic transition. Both worked examples here sit adjacent to the hierarchy rather than inside it, and neither number (148, 2.2×10⁴, or the factor-of-12 symmetry ratio at N=6) is a settled free-energy ratio yet. The contact-potential model itself is also a simplification — real colloidal systems have Lennard-Jones, DLVO, or depletion potentials that modify contact energies and introduce range-dependence, so ε, E_a, and the entropy terms all require system-specific calculation or experiment before this becomes more than semi-quantitative.

Conclusion

Thermodynamics doesn't yet give Ffellonics a settled set of numbers — it gives the hierarchy a precise method for generating them, and an honest accounting of which pieces of that method are still open. The five-sphere example shows the contact-maximizing configuration favored by a known energy term, exp(ε/kT), bounded above by roughly 150-to-22,000-fold at realistic interaction strengths, with two distinct, independently-documented entropic effects — floppy-mode configurational entropy and rotational symmetry-number suppression — working against that margin by amounts not yet computed. The six-sphere example shows that even comparing two fully rigid structures doesn't remove the symmetry-number effect, which can run to a factor of 12 or more before vibrations are even considered. That's a narrower, more defensible claim than a settled prediction, and it points directly at the calculations — the floppy-manifold integral, the symmetry-corrected partition functions, then a grand-canonical treatment of the real milestone transitions — that would turn this from an illustration of method into a tested prediction.