{"id":"2f435494-8d61-4593-8d0e-0b5c9d347ecb","arxiv_id":"1908.02091","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Even infinitesimal asphericity changes jamming critical exponents: contact number excess grows as (A-1)^{1/4}, shear modulus becomes linear in pressure near the transition, and gap/force power-law singularities are smoothed.","lead":"This paper argues that adding even a tiny amount of asphericity to jammed particles changes the entire critical behavior of the jamming transition, giving new exponents for contact number, shear modulus, and vibrational spectra. The result unifies non-spherical and breathing particle models in a new universality class distinct from spherical jamming.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The universality claim is conditional on a sign assumption (c1>0, c2<0) that the authors themselves state fails for dimers and Platonic solids, so δz∼Δ^{1/2} is not established for non-spherical particles generally.","rationale":"The reader's weakest assumption identifies the same load-bearing concern: Eq. (21) depends on the sign of c1 and c2, and the paper itself concedes that c1 > 0 is violated for dimers and Platonic solids. I agree with this assessment. I considered whether a stronger concern exists, such as whether all terms in Eq. (20) must be proportional to p, or whether the zero-mode counting in Eq. (18) assumes independence of the QN constraints. These are noted in the text and supported for ellipsoids and breathing particles by prior work, so they are less central than the sign assumption. The sign assumption is the weakest link because it is shape-dependent and unrestricted in the abstract. The numerical scaling collapses in Figs. 1, 2, and 4 are real evidence for the c1 > 0 class, and the paper's explicit enumeration of counterexamples is a credit to its honesty. The verdict should remain CONDITIONAL, not REJECT: the paper establishes a well-defined universality class for shapes satisfying the sign condition, but it overstates the generality in the abstract. No change to the reader's conditional verdict is needed, hence UNCHANGED.","tokens_in":11475,"tokens_out":3811,"duration_ms":42678,"concrete_test":"Directly compute c1 and c2 for dimers and Platonic solids. For a jammed configuration at small pressure p and small asphericity Δ, construct the Hessian restricted to the zero modes, project onto the rotational/breathing sector, and extract the coefficients of pΔδz² and pΔ². If c1 ≤ 0 for dimers, Eq. (21) fails for that shape, and the abstract's general statement is false. A complementary test: simulate jammed dimers over at least a decade of Δ and plot z − 2d versus Δ at the jamming point; if the scaling exponent is not 1/2 (e.g., linear in Δ, or identically zero), the universality class is not universal across non-spherical particles.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that infinitesimal asphericity changes the universality class rests on Eq. (21), δz ∼ Δ^{1/2}, which follows from balancing the two terms in Eq. (20), λ_min ∼ p[c1 Δ δz² + c2 Δ²]. This balance requires c1 > 0 (the first-order pre-stress stabilizes the zero modes) and c2 < 0 (the second-order term destabilizes). The authors explicitly state in Sec. 6 that c1 > 0 'is violated for non-spherical particles consisting of spherical particles, such as dimers' and that Platonic solids are isostatic at jamming. For these shapes the predicted exponents do not follow. Yet the abstract and Sec. 6 present the claim for non-spherical particles without this restriction, stating that 'even an infinitesimal asphericity is enough to change the universality class.' The numerical support (Figs. 1, 2, 4, and 5) is limited to ellipsoid-like and breathing-particle models, which are precisely the cases believed to have c1 > 0. No argument is provided that bounds c1 away from zero or determines its sign for general convex shapes. Therefore the strongest form of the claim is not established; the result is proven conditionally for shapes satisfying the sign assumption. This is not an external disagreement with consensus but an internal limitation of the derivation, and it is acknowledged in the text.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the jamming transition of non-spherical particles and argues that even infinitesimal asphericity changes the universality class relative to spherical particles. Using a variational argument combined with marginal stability, the authors derive scaling laws for the excess contact number δz ∼ Δ^{1/2}, the shear modulus G ∼ p/Δ^{1/2} for p ≪ Δ, the truncation of gap and force distributions, and the characteristic frequencies of the density of states. They support these predictions with numerical simulations of the breathing-particle model and by comparisons with literature data for ellipsoids and other shapes. The paper is presented as a longer, more direct version of the authors' earlier PNAS work, with additional numerical data.","tokens_in":11717,"tokens_out":6684,"duration_ms":68709,"significance":"The central claim, if established, is significant: it would imply that the jamming universality class is discontinuous at the spherical limit, and it unifies non-spherical particles with breathing particles into a single class with hypostatic jamming. The paper's strengths are that the exponents are derived rather than fitted to the new data, that the numerical collapses in Figs. 2, 4, 5, and 6 are convincing, and that the predictions are checked against external data in Refs. [7, 15, 16]. However, the derivation of the key scaling δz ∼ Δ^{1/2} relies on a sign assumption about the pre-stress coefficients that the authors themselves state fails for dimers and Platonic solids; the scope of the universality claim therefore needs to be narrowed or supported by additional argument.","major_comments":[{"comment":"The step from Eq. (20) to Eq. (21) requires c1 > 0 and c2 < 0. The paper explicitly states in Sec. 6 that the assumption c1 > 0 is violated for non-spherical particles made of spherical particles such as dimers, and that Platonic solids are isostatic at jamming. Since the abstract and the opening of Sec. 6 present the universality change for non-spherical particles without this restriction, the strongest form of the claim is not supported by the derivation. Please restrict the claim to shapes satisfying the sign assumption, or provide a criterion for determining the sign of c1 and ideally a numerical check of the balance in Eq. (20).","section":"Sec. 6, Eq. (20)"},{"comment":"The text says Fig. 1 compares Eq. (21) with numerical results for 'various shapes of non-spherical particles' from Ref. [15]. It is not stated which shapes are actually included. If dimers or Platonic solids are among them, this would be in tension with the Sec. 6 caveat and needs discussion; if they are excluded, the figure caption should say so explicitly.","section":"Fig. 1"},{"comment":"The expression for g_i(ui) in Eq. (7) appears to contain both f(r_ij, ui) and f(-r_ij, ui), whereas the first-order gap correction in Eq. (3) involves f(r_ij, ui) + f(-r_ij, uj). Summing Eq. (7) over i does not reproduce the first-order expansion of the potential derived from Eq. (3) unless a factor of 1/2 or a reindexing is intended. Please correct or clarify this step, since QN and its stiffness kR are the basis of Eq. (19).","section":"Eqs. (6)-(7)"}],"minor_comments":[{"comment":"In the abstract, 'a sphericity' should be 'asphericity'.","section":"Abstract"},{"comment":"Please clarify the p- and Δ-dependence of the characteristic frequencies: Fig. 5(b) is at fixed p, giving ω0 ∼ Δ^{1/2}, while Fig. 5(c) is at fixed Δ, giving ω0 ∼ p^{1/2}; the text currently mixes these two dependencies without distinguishing the fixed variable.","section":"Sec. 5, Eq. (52)"},{"comment":"The caption says the data are for 'non-spherical particles', but the text explains that the simulations are of the breathing-particle model; this should be reflected in the caption for clarity.","section":"Fig. 4 caption"},{"comment":"The statement that N_c^Q - N0 = N δz/2 > 0 'indeed suffices to stabilize the zero modes' is not immediate, because extra constraints stabilize the modes only if they are independent; a brief justification of the independence assumption would help.","section":"Sec. 2.4, Eq. (18)"}],"recommendation":"major_revision","confidential_remarks":"The main issue is the mismatch between the broad universality claim in the abstract and the acknowledged limitations of the sign assumption in Sec. 6. If the authors qualify the claim to the class of shapes for which c1 > 0 and c2 < 0, and clarify which shapes are included in Fig. 1, the paper will be substantially stronger. The additional data and direct variational derivation are useful, but the central claim as currently worded overreaches the derivation."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Haru—\n\nThis one is worth your time if you care about jamming. It's the extended version of the Brito-Ikeda-Urbani-Wyart-Zamponi PNAS paper, adding two things: a direct variational derivation of the scaling exponents for non-spherical particles, and new numerics on the breathing-particle model (gap distribution collapse, DOS frequencies, band weights). The variational argument is cleaner than the earlier replica-based mapping, and the data collapses in Figs. 2, 4, 5, and 6 are solid for the cases they test. The band structure of the DOS—three bands with ω0∼Δ^{1/2}, ω1∼Δ, ω2∼Δ^{1/2}—is a nice new result, and the weight counting in Fig. 6 checks out.\n\nThe soft spot is the scope of the claim. The derivation hinges on the sign of the first-order pre-stress coefficient, c1>0, and the second-order term c2<0. The authors admit in Sec. 6 that c1>0 fails for dimers and Platonic solids; for those shapes the system is isostatic and the exponents revert to the spherical case. Yet the abstract and the summary state the universality change for non-spherical particles with no caveat. That's an overgeneralization, and it's the main thing a referee should press on. The derivation itself is a scaling argument, not a proof, so it can't settle the sign for general convex shapes. The numerical support is also limited to ellipsoid-like and breathing particles—exactly the cases believed to have c1>0. So the strong form of the claim is not established; the conditional form holds up.\n\nMinor: the paper relies heavily on the authors' own previous work, but the new derivation and numerics do add value. No code or data files are included, though the methods are detailed in the earlier paper.\n\nBottom line: the paper is a solid contribution within its true domain of validity. It should be refereed, not desk-rejected. The revision should either narrow the universality claim explicitly or provide an argument bounding c1 away from zero for a broader class of shapes.","headline":"A solid extension of the authors' PNAS 2018 result with a clearer variational derivation and good breathing-particle numerics, but the universality claim is overstated: the key sign condition fails for dimers and Platonic solids.","tokens_in":12317,"tokens_out":1934,"would_cite":true,"duration_ms":18491,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Even a tiny asphericity changes the universality class of the jamming transition.","keywords":["jamming transition","non-spherical particles","asphericity","universality class","marginal stability","shear modulus","gap distribution","breathing particles"],"falsifier":"Numerically simulate slightly aspherical ellipsoids at fixed small $\\Delta$ and measure the shear modulus as a function of pressure down to $p$ much smaller than $\\Delta$: the predicted linear regime $G \\sim p$ must appear, with crossover to $G \\sim p^{1/2}$ near $p \\sim \\Delta$. If the square-root law persists all the way to $p=0$ at fixed $\\Delta>0$, the central scaling claim is wrong. Equivalently, the excess contact number at jamming should scale as $(A-1)^{1/4}$; a different exponent would falsify Eq. (21).","tokens_in":11176,"feed_emoji":"","tokens_out":10023,"duration_ms":80937,"temperature":0.7,"pith_summary":"Jamming of frictionless particles is usually described by universal power laws set by isostaticity, but this paper argues that those laws are special to perfect spheres. The authors claim that any nonzero asphericity $\\Delta$ puts the transition in a different universality class: particles become hypostatic, the excess contact number scales as $\\delta z \\sim \\Delta^{1/2}$, and the shear modulus becomes linear in pressure sufficiently close to jamming rather than following the spherical square-root law. They further claim the power-law singularities in the gap and force distributions are rounded at scales set by $\\Delta$, and that breathing particles—spheres whose radii can fluctuate—fall in the same class. The practical consequence is that a tiny shape anisotropy should be visible in the critical exponents measured in experiments and simulations.","feed_headline":"A trace of asphericity changes jamming's universality class","feed_subtitle":"For any nonzero asphericity, the shear modulus near jamming turns linear in pressure, not square-root.","key_machinery":"The argument is carried by a variational estimate of the smallest vibrational eigenvalue combined with the assumption of marginal stability. The interaction potential is expanded in a small shape parameter $\\Delta$ around a reference sphere; the leading asphericity introduces a term $Q_N \\sim p\\Delta$ that stabilizes the rotational zero modes, with stiffness $k_R \\sim p\\Delta$. Applying the variational argument of spherical jamming to these zero modes gives $\\lambda_{\\min} \\sim p(c_1 \\Delta \\delta z^2 + c_2 \\Delta^2)$, where the first term is the stabilizing pre-stress and the second (negative, $c_2<0$) is the destabilizing buckling contribution. Requiring $\\lambda_{\\min} \\sim 0$ (marginal stability) balances the two terms and yields $\\delta z \\sim \\Delta^{1/2}$. This single relation then controls the contact number, the shear modulus via $G \\sim p\\delta z/\\Delta$, and the truncation scales of the gap and force distributions through finite-size scaling with $\\delta z$ in place of $1/N$.","core_discovery":"The paper's central claim is that the jamming transition of non-spherical particles belongs to a universality class distinct from that of spherical particles, and that the distinction survives in the limit of infinitesimal asphericity. Concretely, for any fixed $\\Delta > 0$, sufficiently near the transition ($p \\ll \\Delta$) the shear modulus behaves as $G \\sim p/\\Delta^{1/2}$, linear in pressure, whereas for spheres $G \\sim p^{1/2}$; the crossover between the two regimes occurs at $p \\sim \\Delta$. The excess contact number at the jamming point scales as $\\delta z \\sim \\Delta^{1/2} \\sim (A-1)^{1/4}$, and above jamming $z-z_J \\sim p/\\Delta^{1/2}$. The gap distribution $g(h)$ and force distribution $P(f)$, which for spheres diverge as power laws at the transition, remain finite and analytic for $\\Delta > 0$, with the singularities truncated at $h \\sim \\Delta^{\\mu/2}$ and $f \\sim \\Delta^{\\nu}$. The density of states splits into three bands with characteristic frequencies $\\omega_0 \\sim \\Delta^{1/2} p^{1/2}$, $\\omega_1 \\sim \\Delta$, and $\\omega_2 \\sim \\Delta^{1/2}$. The same exponents are found for breathing particles, and the paper presents numerical data supporting the predictions.","pith_inferences":["If the claim holds, any experimental granular material composed of slightly elongated grains should show a linear pressure dependence of the shear modulus in the deeply unjammed regime; this could be checked with photoelastic disks or 3D ellipsoidal particles.","The sign of the first-order stabilizing coefficient becomes a shape classifier: shapes with positive sign should follow the new universality class, while shapes where it is zero or negative (dimers, Platonic solids) remain isostatic—so the universality class may be selected by the geometry of the rotational modes.","The same marginal-stability balance might extend to other internal degrees of freedom, such as particle deformability or internal orientational fields, whenever the first-order stabilizing term is positive.","Because the gap distribution is regularized, force-network criticality and associated avalanches near jamming may be suppressed at scales below $\\Delta$, a testable prediction for mesoscopic simulations."],"forward_implications":["For any nonzero asphericity, the shear modulus of a jammed packing is linear in pressure in the asymptotic unjamming limit, changing the low-frequency mechanics from the spherical square-root form.","The gap and force distributions are regular functions at the jamming point for aspherical particles, so the power-law singularities characteristic of isostatic jamming disappear for any $\\Delta>0$.","The density of states contains three distinct bands with characteristic frequencies $\\omega_0 \\sim \\Delta^{1/2}p^{1/2}$, $\\omega_1 \\sim \\Delta$, and $\\omega_2 \\sim \\Delta^{1/2}$, rather than the single soft mode scale of spheres.","Breathing particles reproduce the same critical exponents, giving a rotationally symmetric model in which the predictions can be tested efficiently.","The contact number at jamming increases with asphericity as $z_J = 2d + O(\\Delta^{1/2})$, and above jamming $z-z_J$ grows linearly in $p/\\Delta^{1/2}$."],"supporting_citations":[{"why":"previous joint work establishing the breathing-particle mapping and the claim of a shared universality class, which this paper extends with a direct variational derivation","marker":"[5]"},{"why":"numerical data on ellipsoids showing the linear-in-pressure shear modulus and its dependence on aspect ratio","marker":"[7]"},{"why":"numerical density of states for ellipsoids whose characteristic frequency scalings this paper reproduces","marker":"[8]"},{"why":"source of the breathing-particle model and of the relation $k \\sim p/\\Delta$ used in the argument","marker":"[10]"},{"why":"the variational marginal-stability argument for spherical jamming that gives $\\delta z \\sim p^{1/2}$; the non-spherical argument is built on it","marker":"[11]"},{"why":"the variational estimate $\\lambda_{\\min} \\sim k\\delta z^2$ for nearly isostatic packings, used as Eq. (13)","marker":"[12]"},{"why":"numerical baseline for spherical jamming ($G \\sim p^{1/2}$, gap distributions) against which the aspherical scaling is compared","marker":"[14]"},{"why":"numerical data for many particle shapes supporting $z-4 \\sim (A-1)^{1/4}$","marker":"[15]"},{"why":"simulations of ellipsoids supporting $z-z_J \\sim p/\\Delta^{1/2}$","marker":"[16]"}],"fun_headline_variants":["Any asphericity changes jamming's universality","Tiny shape deviation alters jamming critical exponents","Jamming transition: infinitesimal asphericity matters","Non-spherical particles jam with different criticality","A sliver of asphericity smooths jamming singularities"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument assumes that the first-order effect of asphericity stabilizes the rotational zero modes; for shapes where that stabilization is absent or negative—the paper names dimers and Platonic solids—the predicted exponents do not follow.","fun_headline_variants_meta":{"raw":{"variants":["Any asphericity changes jamming's universality","Tiny shape deviation alters jamming critical exponents","Jamming transition: infinitesimal asphericity matters","Non-spherical particles jam with different criticality","A sliver of asphericity smooths jamming singularities"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000228,"raw_usage":{"total_tokens":1544,"prompt_tokens":1083,"completion_tokens":461,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":699,"completion_tokens_details":{"reasoning_tokens":384}},"tokens_in":699,"tokens_out":461,"duration_ms":5338,"temperature":1.0,"reasoning_tokens":384,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:54:15.003607+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Numerically simulate slightly aspherical ellipsoids at fixed small $\\Delta$ and measure the shear modulus as a function of pressure down to $p$ much smaller than $\\Delta$: the predicted linear regime $G \\sim p$ must appear, with crossover to $G \\sim p^{1/2}$ near $p \\sim \\Delta$. If the square-root law persists all the way to $p=0$ at fixed $\\Delta>0$, the central scaling claim is wrong. Equivalently, the excess contact number at jamming should scale as $(A-1)^{1/4}$; a different exponent would falsify Eq. (21).","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"previous joint work establishing the breathing-particle mapping and the claim of a shared universality class, which this paper extends with a direct variational derivation"},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"numerical data on ellipsoids showing the linear-in-pressure shear modulus and its dependence on aspect ratio"},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"numerical density of states for ellipsoids whose characteristic frequency scalings this paper reproduces"},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"source of the breathing-particle model and of the relation $k \\sim p/\\Delta$ used in the argument"},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"the variational marginal-stability argument for spherical jamming that gives $\\delta z \\sim p^{1/2}$; the non-spherical argument is built on it"},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"the variational estimate $\\lambda_{\\min} \\sim k\\delta z^2$ for nearly isostatic packings, used as Eq. (13)"},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"numerical data for many particle shapes supporting $z-4 \\sim (A-1)^{1/4}$"},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"simulations of ellipsoids supporting $z-z_J \\sim p/\\Delta^{1/2}$"}],"review_version":1}