pith. sign in

arxiv: 2509.20132 · v3 · submitted 2025-09-24 · ❄️ cond-mat.soft · cond-mat.dis-nn· math-ph· math.MP

Random close packing fraction of bidisperse discs: Theoretical derivation and exact bounds

Raphael Blumenfeld This is my paper

Pith reviewed 2026-05-18 14:09 UTC · model grok-4.3

classification ❄️ cond-mat.soft cond-mat.dis-nnmath-phmath.MP
keywords random close packingbidisperse discspacking fractiondisordered packingscell order distributionexact boundsgranular materials
0
0 comments X

The pith

A disorder-guaranteeing theory using cell order distributions derives the highest possible random close packing fraction for bidisperse discs along with exact bounds.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a theoretical framework to find the maximum density achievable in random, disordered packings of two differently sized discs. It relies on the concept of cell order distribution to ensure the arrangement stays fully disordered while reaching the highest density possible. This matters because bidisperse systems are commonly used in experiments to avoid ordering, and a prediction of their densest disordered state would help in modeling granular and colloidal materials. The approach also yields precise upper and lower bounds that any such packing fraction must lie between.

Core claim

Using a disorder-guaranteeing theory based on the cell order distribution, the highest mathematically possible value of the random close packing fraction φ_RCP(p, D) is derived for bidisperse discs, where p is the concentration and D the size ratio, together with exact upper and lower bounds on this densest disordered packing fraction.

What carries the argument

The cell order distribution, which is used to mathematically guarantee that the packing remains fully disordered while maximizing the overall density.

If this is right

  • φ_RCP can now be expressed as a function of the size ratio D and the concentration p of each disc type.
  • Exact bounds provide rigorous limits that simulations and experiments must respect for disordered states.
  • Predictions for the densest disordered packings in two dimensions become available for comparison with data.
  • The method offers a way to minimize order while achieving high density in bidisperse systems.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If this framework holds, it could be adapted to predict packing fractions in three-dimensional bidisperse sphere systems.
  • Connections might exist to the jamming transition and mechanical properties of granular materials.
  • Testing the bounds against large-scale simulations could reveal how close real systems get to the theoretical maximum.

Load-bearing premise

That it is possible to define a cell order distribution which guarantees the packing is completely disordered at the mathematically highest density.

What would settle it

Finding a disordered bidisperse disc packing with a density higher than the derived maximum φ_RCP(p,D) or outside the exact bounds, as verified by an order parameter that confirms absence of crystalline order.

Figures

Figures reproduced from arXiv: 2509.20132 by Raphael Blumenfeld.

Figure 1
Figure 1. Figure 1: FIG. 1. The four possible 3-cell configurations and their areas [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. The packing fraction, Φ [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. A typical example of Φ [PITH_FULL_IMAGE:figures/full_fig_p003_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. As in Fig. 5 for several size ratios 1 [PITH_FULL_IMAGE:figures/full_fig_p004_5.png] view at source ↗
read the original abstract

A long-standing problem has been a theoretical prediction of the densest packing fraction of random packings, $\phi_{RCP}$, of same-size discs in $d=2$ and spheres in $3$. However, to minimize order, experiments and numerical simulations often use two-size discs and a prediction of the highest possible packing fraction, $\phi_{RCP}$, for these packings could be very useful. In such bidisperse packings, $\phi_{RCP}$ is a function of the sizes ratio, $D$, and concentrations, $p$, of the disc types. A disorder-guaranteeing theory is formulated here to derive the highest mathematically possible value of $\phi_{RCP}(p,D)$, using the concept of the cell order distribution. I also derive exact upper and lower bounds on this densest disordered packing fraction.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript formulates a disorder-guaranteeing theory for bidisperse disc packings that derives the highest mathematically attainable random close packing fraction φ_RCP(p,D) as a function of size ratio D and concentration p. The derivation relies on the cell order distribution to enforce disorder while maximizing density, and the paper also supplies exact upper and lower bounds on this densest disordered packing fraction.

Significance. A rigorous, parameter-free derivation of φ_RCP(p,D) together with exact bounds would constitute a substantial advance for the long-standing problem of predicting random close packing in two dimensions. If the cell-order construction truly guarantees the absence of crystalline domains or long-range order at the reported maximum density, the result would be directly useful for interpreting simulations and experiments on bidisperse systems and could serve as a benchmark for future numerical work.

major comments (2)
  1. [Cell order distribution definition and optimization] The central claim that the cell order distribution mathematically guarantees a fully disordered state at maximum density is load-bearing for the entire derivation. The manuscript must demonstrate explicitly (e.g., in the section defining the distribution and the subsequent optimization) that no configuration satisfying the distribution can contain small crystalline domains or emergent long-range correlations that would raise the density above the accepted RCP threshold while still obeying the local cell statistics.
  2. [Bounds derivation] The exact upper and lower bounds on φ_RCP(p,D) are presented as rigorous; however, the proof that these bounds remain valid under the disorder constraint (rather than relaxing to ordered packings) is not shown in sufficient detail. A concrete counter-example or a step-by-step argument that any packing achieving the upper bound must violate the cell-order condition if it contains crystalline order would strengthen the result.
minor comments (2)
  1. [Introduction] Notation for the cell order distribution and the parameters p and D should be introduced with a clear table or equation early in the text to aid readability.
  2. [Results] The manuscript would benefit from a short comparison, even qualitative, between the derived φ_RCP(p,D) and existing simulation data for a few representative (p,D) pairs.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their thorough review and valuable suggestions. We address the major comments point by point below and will revise the manuscript accordingly to strengthen the presentation of our results.

read point-by-point responses

  1. Referee: [Cell order distribution definition and optimization] The central claim that the cell order distribution mathematically guarantees a fully disordered state at maximum density is load-bearing for the entire derivation. The manuscript must demonstrate explicitly (e.g., in the section defining the distribution and the subsequent optimization) that no configuration satisfying the distribution can contain small crystalline domains or emergent long-range correlations that would raise the density above the accepted RCP threshold while still obeying the local cell statistics.

    Authors: We agree that an explicit demonstration is crucial for the validity of our claims. In our manuscript, the cell order distribution is defined such that it only includes local configurations consistent with random, disordered packings, excluding those that would allow for crystalline ordering. The optimization maximizes the packing fraction subject to these constraints. To clarify this, we will add a dedicated subsection in the revised manuscript that proves by contradiction: suppose a configuration satisfying the cell order distribution contains a crystalline domain; then the local cell orders in that domain would have probabilities outside the support of our distribution, leading to a contradiction. This ensures that no such domains or correlations can form while obeying the distribution, thereby guaranteeing the disorder at the maximum density. revision: yes

  2. Referee: [Bounds derivation] The exact upper and lower bounds on φ_RCP(p,D) are presented as rigorous; however, the proof that these bounds remain valid under the disorder constraint (rather than relaxing to ordered packings) is not shown in sufficient detail. A concrete counter-example or a step-by-step argument that any packing achieving the upper bound must violate the cell-order condition if it contains crystalline order would strengthen the result.

    Authors: We appreciate this suggestion for strengthening the rigor of our bounds. Our derivation of the bounds incorporates the disorder constraint directly through the cell order distribution, meaning that the upper bound is the maximum achievable only within the space of disordered configurations. We will include in the revision a step-by-step argument showing that achieving the upper bound with any crystalline order would necessitate a cell order distribution that includes ordered local structures, which are excluded by our disorder-guaranteeing condition. Additionally, we will provide a concrete counter-example of a bidisperse packing with crystalline domains that exceeds the bound but violates the cell order statistics, thereby confirming that the bound holds specifically for disordered states. revision: yes

Circularity Check

0 steps flagged

No circularity detected; derivation relies on independent mathematical construction of cell order distribution

full rationale

The paper formulates a disorder-guaranteeing theory based on the cell order distribution to derive φ_RCP(p,D) and exact bounds. No load-bearing steps reduce by construction to fitted parameters, self-citations, or ansatzes imported from prior work by the same author. The central derivation is presented as a self-contained theoretical framework that mathematically enforces disorder while maximizing density, with no visible renaming of known results or uniqueness theorems invoked from overlapping authorship. The approach is self-contained against external benchmarks as a first-principles construction.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central construction rests on the existence and utility of a cell order distribution that enforces disorder while maximizing density; no free parameters or new entities are mentioned in the abstract.

axioms (1)
  • domain assumption A cell order distribution exists that can be used to guarantee disorder while deriving the maximum packing fraction.
    This is the explicit foundation of the disorder-guaranteeing theory stated in the abstract.

pith-pipeline@v0.9.0 · 5678 in / 1192 out tokens · 48772 ms · 2026-05-18T14:09:15.000811+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

16 extracted references · 16 canonical work pages · 1 internal anchor

  1. [1]

    Random close packing fraction of bidisperse discs: Theoretical derivation and exact bounds

    ensures no ‘rattlers’ within enclosures of three large discs. In the following, the COD is first defined and the advantages of using it are briefly reviewed. An exact upper bound on Φ RCP is then derived as a function ofpandD. Next, the disorder criterion, developed in [1], is extended and used to identify the range ofp(D), within which packings are assur...

    work page internal anchor Pith review Pith/arXiv arXiv 2025

  2. [2]

    Blumenfeld, Phys

    R. Blumenfeld, Phys. Rev. Lett.127, 118002 (2021)

    work page 2021

  3. [3]

    Torquato, F.H

    S. Torquato, F.H. Stillinger, Rev. Mod. Phys.82, 2633 (2010) and references therein

    work page 2010

  4. [4]

    Bernal, Nature185, 68 (1960)

    J.D. Bernal, Nature185, 68 (1960)

    work page 1960

  5. [5]

    Torquato, T.M

    S. Torquato, T.M. Truskett, P.G. Debenedetti, Phys. 5 Rev. Lett.84, 2064 (2000)

    work page 2064

  6. [6]

    Visscher, M

    W.M. Visscher, M. Bolsterli, Nature239. 504 (1972)

    work page 1972

  7. [7]

    Bideau, A

    D. Bideau, A. Gervois, L. Oger, J. P. Troadec, J. Physique47, 1697 (1986)

    work page 1986

  8. [8]

    Okubo, T

    T. Okubo, T. Odagaki, J. Phys.: Condens. Matter16, 6651 (2004)

    work page 2004

  9. [9]

    N. Xu, J. Blawzdziewicz, C.S. O’Hern, Phys. Rev. E71, 061306 (2005)

    work page 2005

  10. [10]

    D. J. Koeze, D. V˚ agberg, B.B.T. Tjoa, B.P. Tighe, EPL 113, 54001 (2016)

    work page 2016

  11. [11]

    Brouwers, arxiv: 2501.16005 (2025) and refer- ences therein

    H.J.H. Brouwers, arxiv: 2501.16005 (2025) and refer- ences therein

    work page arXiv 2025

  12. [12]

    Matsushima, R

    T. Matsushima, R. Blumenfeld, Phys. Rev. Lett.112, 098003 (2014)

    work page 2014

  13. [13]

    Wanjura, P.A

    C.C. Wanjura, P.A. Gago, T. Matsushima, R. Blumen- feld, Granular Matter22, 91 (2020)

    work page 2020

  14. [14]

    Supplemental material containing detailed calculations and more plots

    work page

  15. [15]

    Fernique, Comp

    T. Fernique, Comp. Geom. Volumes124–125, 102134 (2025)

    work page 2025

  16. [16]

    Matsushima, R

    T. Matsushima, R. Blumenfeld, Phys. Rev. E95, 032905 (2017)

    work page 2017