{"id":"8d0629ce-b6d5-4d16-a1d8-55822a1bbdf7","arxiv_id":"2507.11179","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"A pressure-tuned transition in 3D amorphous solids shows an intermediate 'anomalous elasticity' phase, where angular correlations of strain-induced displacement diverge with exponent about 1.66.","lead":"This paper reports a pressure-driven transition in how 3D amorphous solids respond to strain, with an intermediate phase where the displacement field is screened and shows anomalous, sign-changing behavior. The authors find that angular correlations of the displacement diverge near the transition, with a critical exponent around 1.66.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The divergence in Eq. (11) may be imposed by the fixed-exponent generalized Lorentzian fit, since for integer angular modes and f† < 1 the fit is constrained mainly by the dc normalization, not by a resolved correlation length.","rationale":"The reader's weakest assumption correctly identifies the fixed-exponent generalized Lorentzian fit as the vulnerable point. My analysis sharpens this: for the angular Fourier variable in Eq. (8), integer wavenumbers are the natural modes, and when f† < 1 the nonzero modes all sample the f^{−2.85} tail. In that regime the fit cannot independently determine f† from the spectral shape; it effectively uses the dc amplitude and the imposed exponent. This makes the reported divergence in Eq. (11) a property of the fitting ansatz unless alternatives are tested. The paper does provide supporting evidence for the qualitative three-phase picture: the exponential decay below jamming, the power-law quasi-elastic response, and the intermediate-phase fit to Eq. (7) with κ ≈ 0.256–0.285 are direct and plausibly robust. The screening theory also has prior support from the group's earlier work. However, the critical claim—the diverging angular correlation length and exponent μ ≈ 1.66—requires the fitting-form robustness check. The discrepancy among pc estimates (2.3, 2.4, 2.9) further weakens the exponent determination. This is an addressable data-analysis concern, not a fatal flaw, so the reader's CONDITIONAL verdict remains appropriate and I do not recommend changing it.","tokens_in":8448,"tokens_out":9277,"duration_ms":132190,"concrete_test":"Refit the S(r,f) data underlying Fig. 4 with a free-exponent generalized Lorentzian S = A/[1 + (f/f†)^η] and with a pure power law f^{−η}, comparing fits over the same f range via AIC or χ². If the best-fit η(p) drifts away from 2.85 near pc, or if a pure power law fits as well as the Lorentzian with f† > 0, then the divergence of θ† in Eq. (11) is a consequence of the fixed exponent and normalization, not a measured critical length.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Eq. (10) fixes the spectral tail exponent at 2.85 and leaves f† as the only free parameter. The divergence θ† = 1/f† → ∞ in Eq. (11) can therefore be an artifact of the fitting form: any growth of the zero-frequency peak relative to the tail can be absorbed into a shrinking f†. On a circle, the angular Fourier variable in Eq. (8) is discrete; once f† < 1, all nonzero integer modes f = 1,2,... lie in the asymptotic tail, where SGL(f) ≈ f^{−2.85} is independent of f†. The value of f† is then controlled only by the normalization A = 1 + f†^{−2.85} (i.e., by the ratio S(0)/S(1)) and by the assumed 2.85 exponent. The paper reports no test against a pure power law, a free exponent, an exponential cutoff, or a dependence on shell radius r and system size rout. The text itself states that pc depends on d0, rin, and rout, and the three estimates pc ≈ 2.3, 2.4, and 2.9 are not tightly pinned; μ ≈ 1.66 is extracted relative to this uncertain pc. Thus Eq. (11) is currently a property of the imposed spectral shape rather than an established critical divergence.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This Letter reports simulations of athermal, three-dimensional packings of Hertzian spheres subjected to inflation of a central sphere, and identifies three pressure regimes in the radial displacement response: exponential decay at pressures below jamming, quasi-elastic r^{-2} decay at high pressure, and an intermediate pressure regime in which the radial displacement changes sign and is fitted by a screened elastic solution with a screening parameter kappa (Eq. 7). The central claim is that the transition between the intermediate and quasi-elastic phases is critical, with an angular correlation length theta-dagger diverging as theta-dagger proportional to (p_c - p)^{-mu} with mu approximately 1.66 and p_c approximately 2.9 (Eq. 11).","tokens_in":8787,"tokens_out":8018,"duration_ms":94635,"significance":"If established, the existence of a dipole-induced transition in the three-dimensional mechanical response of amorphous solids would be conceptually significant, extending the analogy with two-dimensional topological transitions to a new setting. The qualitative observation of three distinct radial response forms, including the sign reversal in the intermediate phase, is a valuable empirical result. However, the central critical-scaling claim is not yet supported by the presented evidence: the divergence is extracted from a fitting form whose single parameter can produce an apparent divergence by construction, and the paper does not provide the statistical or finite-size analysis needed to distinguish a true divergence from a crossover. The paper's strengths are its direct simulation observations and the clear presentation of the three regimes; its weakness is the lack of any test of the assumed spectral shape or of the robustness of the extracted critical exponents.","major_comments":[{"comment":"The divergence in Eq. (11) is obtained by fitting S(r,f) to the generalized Lorentzian SGL(r,f) = A/(1 + (f/f-dagger)^2.85) with A = 1 + (1/f-dagger)^2.85, leaving f-dagger as the only free parameter. Because f is a discrete angular Fourier index, once f-dagger < 1 all nonzero integer modes lie in the tail, where SGL(f) is approximately f^{-2.85} and independent of f-dagger. In that regime f-dagger is determined only by the normalization ratio S(0)/S(1), so any growth of the zero-frequency peak relative to the first nonzero mode is automatically encoded as a shrinking f-dagger and hence a diverging theta-dagger = 1/f-dagger. The paper reports no test against a pure power law, a free spectral exponent, an exponential cutoff, or a dependence on shell radius r and outer radius r_out. Since Eq. (11) is the paper's main result, this is a load-bearing issue that must be addressed with alternative fits and a demonstration that the divergence is not an artifact of the assumed form.","section":""},{"comment":"No error bars, number of independent configurations, or finite-size scaling analysis are reported for the angular correlation length theta-dagger(p). The text itself states that p_c depends on d_0, r_in, and r_out, and that the determination of p_c from the kappa jump is only approximate; the three estimates p_c = 2.3, 2.4, and 2.9 are not tightly pinned. The exponent mu = 1.66 is extracted in a log-log plot against p_c - p with p_c chosen from the same data. Without a quantitative characterization of the statistical uncertainty and a check of whether theta-dagger saturates with system size, the 'apparent critical divergence' is not distinguishable from a smooth crossover or a fitting artifact.","section":""},{"comment":"The estimate p_c = 2.4 is not an independent prediction because it uses the measured value of the screening parameter kappa = 0.256 through the relation Delta Z = 6 kappa, and that relation is imported from Ref. [36], a study of two-dimensional frictionless disks. Applying a two-dimensional relation to three-dimensional Hertzian spheres requires a derivation or at least a justification that the numerical coefficient remains valid in 3D. Moreover, since kappa itself is fitted from the simulation data, the agreement among p_c = 2.3, 2.4, and 2.9 is a consistency check of scaling relations rather than a derivation of the critical pressure. The authors should either derive the Delta Z-kappa relation in three dimensions or present this estimate as a heuristic consistency argument rather than as a theoretical prediction.","section":""},{"comment":"The choice of the second singularity of Eq. (7), kappa = 0.29, over the first, kappa = 0.145, is justified only by the statement that the first value is too small for the system size. This selection rule is not quantified: no criterion is given for what 'too small' means in terms of r_out, and no demonstration is provided that the second singularity is the relevant one for the observed packing. Because this selection underpins the theoretical interpretation of the intermediate phase and the associated p_c estimate, the selection mechanism needs a quantitative formulation or at least a systematic test against system-size variation.","section":""}],"minor_comments":[{"comment":"The caption contains a typo: 'Pnael a:' should be 'Panel a:'.","section":""},{"comment":"The boundary-condition notation is garbled: '(.rin, theta, phi) = d0' should read 'd(r_in, theta, phi) = d0' (and similarly for r_out).","section":""},{"comment":"The text reports kappa = 0.285 for the fit in Fig. 1(d) but later uses kappa = 0.256 as the measured average (Fig. 3 caption and the theoretical estimate section). The discrepancy between these two values should be explained or reconciled.","section":""},{"comment":"The variable f in the Fourier transform is never defined as an integer angular mode. Since the transform is taken along a circle, this should be stated explicitly, as it is relevant for interpreting the fits in Eq. (10).","section":""},{"comment":"The sentence 'this intermediate phase exists both in two and three dimensions' is asserted without a supporting reference or figure in this manuscript; it should be either substantiated or explicitly attributed to prior work.","section":""},{"comment":"Eq. (12) states p ~ (phi - phi_J)^{3/2} ~ Delta Z^3; the relation Delta Z ~ p^{1/3} should be checked against the Hertzian-sphere literature, as the exponent for the excess coordination number in three dimensions may differ from the simple form assumed here.","section":""}],"recommendation":"major_revision","confidential_remarks":"The paper's qualitative phase characterization is interesting and likely worth publishing in some form, but the central critical-scaling claim needs substantially stronger evidence. The authors should be asked to reanalyze their angular spectra with free exponents and alternative functional forms, to report uncertainties and finite-size effects, and to either derive the Delta Z - kappa relation in 3D or downgrade the theoretical p_c estimate to a consistency check. If such analysis cannot be provided within a reasonable revision, the claim of a critical divergence should be removed or reframed as a preliminary observation."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe honest take: this paper has a genuinely interesting qualitative observation—three distinct response regimes in 3D amorphous solids, with an intermediate phase where the radial displacement is screened and even reverses sign. That part is convincing and is a real step beyond the group's earlier work. The new quantitative claims—the exponent mu ≈ 1.66 and the pressure dependence of the angular correlation length—are presented as the central result, but they rest on a fit that may be manufacturing the divergence.\n\nI think the stress-test note is right. Eq. (10) fixes the spectral tail exponent at 2.85 and leaves f-dagger as the only free parameter. Since the angular modes on a circle are discrete integers, once f-dagger drops below 1, all nonzero modes sit in the asymptotic tail, where the Lorentzian is essentially f^{-2.85} independent of f-dagger. The fit is then controlled by the ratio S(0)/S(1) and the assumed exponent. The paper reports no test with a free exponent, a pure power law, an exponential cutoff, or variation with shell radius and system size. So Eq. (11) may currently be a property of the imposed spectral shape. The authors are candid that p_c depends on d0, r_in, r_out and give three estimates (2.3, 2.4, 2.9) that don't tightly pin the critical point; mu ≈ 1.66 is extracted relative to an uncertain p_c, with no error bars. The theoretical estimate of p_c imports Delta Z ≈ 6 kappa from a 2D disk study [36] into 3D Hertzian spheres; that's a stretch, though not fatal.\n\nTo give credit where due: the simulations look clean, the distinction between exponential, power-law, and screened responses is well demonstrated, and the paper does not overstate certainty—it explicitly flags the approximate nature of p_c. The qualitative intermediate phase is already in the group's prior work, so the novelty is incremental rather than a completely new phenomenon.\n\nBottom line: this deserves serious peer review, but a referee should push hard for a robustness analysis of the spectral fit. If the divergence survives a free exponent and system-size checks, it's a nice result. If not, the qualitative phase picture still has value. I'd bring it to a reading group to discuss the fitting issue, but I wouldn't cite the divergence claim in my own work until it's hardened.\n\nRecommendation: send to review, with the expectation of major revisions.","headline":"A plausible three-phase picture with a genuinely interesting qualitative observation, but the headline critical divergence is likely imposed by the fixed-exponent Lorentzian fit; worth refereeing, but the divergence claim needs hardening.","tokens_in":9291,"tokens_out":2860,"would_cite":false,"duration_ms":33278,"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":"This paper reports a dipole-induced critical transition in the three-dimensional mechanical response of amorphous solids, marked by a diverging angular correlation length at a pressure near 2.9.","keywords":["amorphous solids","anomalous elasticity","dipole screening","quadrupolar plastic events","angular correlations","jamming transition","Hertzian spheres","critical scaling"],"falsifier":"Measure $S(r,f)$ for the same packings at pressures just below $p_c$, but fit it without fixing the exponent 2.85; if the best-fit $f^\\dagger$ approaches a positive constant rather than zero, the claimed critical divergence is an artifact of the fitting form. A second check is to repeat the inflation experiment in larger systems and see whether the fitted exponent $\\mu$ stays near 1.66 and whether the divergence sharpens with system size, as a true critical point should.","tokens_in":8219,"feed_emoji":"🧲","tokens_out":9285,"duration_ms":103856,"temperature":0.7,"pith_summary":"By inflating the central sphere in a simulated packing of Hertzian spheres and measuring the resulting displacement field, this paper argues that a three-dimensional amorphous solid has a pressure-tuned, dipole-induced critical transition in its mechanical response. At high pressure the response is quasi-elastic, with the radial displacement decaying as a power law; at zero pressure it is fluid-like, with exponential decay. In between, gradients of quadrupolar plastic events act as effective dipole charges that screen elasticity, producing an intermediate \"anomalous elasticity\" phase in which the radial response remains long-ranged but angular correlations are lost. The central result is that the angular correlation length diverges as $\\theta^\\dagger \\propto (p_c-p)^{-\\mu}$ with $\\mu \\approx 1.66$ and $p_c \\approx 2.9$, so the phase boundary is a genuine critical point. If correct, this provides a three-dimensional analogue of the hexatic and Kosterlitz-Thouless transitions, realized in mechanical response rather than in thermal ordering.","feed_headline":"Angular order diverges at a pressure-tuned 3D solid transition","feed_subtitle":"Strained amorphous packings show a critical intermediate phase with correlation length growing as (pc−p)^−1.66.","key_machinery":"The central object is the angular power spectrum of the radial displacement field in thin spherical shells around the inflated particle, built from the transform $X(r,\\phi,f)=(1/M)\\sum_m d_r(r_m,\\theta_m,\\phi)e^{-if\\theta_m}$ and its ensemble-averaged magnitude squared $S(r,f)=|\\langle X\\rangle|^2$. The physical mechanism is dipole screening: plastic events are quadrupolar Eshelby inclusions forming a field $Q_{\\alpha\\beta}(\\mathbf{r})$, and their gradients define effective dipole charges $P_\\alpha=\\partial_\\beta Q_{\\alpha\\beta}$, which add a $k^2\\mathbf{d}$ term to the elasticity equation and break translational symmetry. The transition is read off by fitting $S(r,f)$ to a generalized Lorentzian $S_{GL}(r,f)=A/[1+(f/f^\\dagger)^{2.85}]$; as the spectrum approaches a pure power law, $f^\\dagger\\to0$ and $\\theta^\\dagger=1/f^\\dagger$ diverges. The observed screening parameter $\\kappa\\approx0.256$ is understood through the requirement that the screened solution sits near one of its discrete singular maxima, the second singularity ($\\kappa\\approx0.29$) being selected because the first would give a screening length too large for the system.","core_discovery":"The paper establishes that the mechanical response of athermal amorphous packings to a central inflation is controlled by a critical pressure $p_c$. Above $p_c$, the classical Lamé equation $\\mu\\Delta\\mathbf{d}+(\\lambda+\\mu)\\nabla(\\nabla\\cdot\\mathbf{d})=0$ describes the displacement field, which decays as $1/r^2$ in the bulk. Below $p_c$, quadrupolar plastic events (Eshelby inclusions) form gradients that act as effective dipoles $P_\\alpha=\\partial_\\beta Q_{\\alpha\\beta}$, adding a screening term $k^2\\mathbf{d}$ to the Lamé equation; the radial solution then takes the screened form involving spherical Bessel functions, with the screening parameter $\\kappa$ jumping from zero to about $0.256$. The critical signature is the divergence of the angular correlation length $\\theta^\\dagger=1/f^\\dagger$, extracted from the angular power spectrum $S(r,f)$, which scales as $\\theta^\\dagger\\propto(p_c-p)^{-\\mu}$ with $\\mu\\approx1.66$ and $p_c\\approx2.9$ near the transition.","pith_inferences":["If the divergence survives larger systems, the exponent $\\mu\\approx1.66$ becomes a new, testable critical exponent of disordered solids, and comparing it across frictional, Hertzian, and differently prepared packings would reveal whether the transition is universal.","Because the avalanche size $\\ell$ depends on the interaction law while $\\kappa^{-1}$ is mostly geometric, the same mechanism predicts that the critical pressure shifts with particle softness and friction even if the divergence itself persists; this is a concrete prediction beyond the paper's explicit claims.","The combination of retained radial range and lost angular correlations suggests the transition may admit an effective lower-dimensional description, in which case $\\mu$ might be connected to known correlation-length exponents of one-dimensional or planar models.","A practical corollary, if the transition is genuine, is that a packing approaching $p_c$ from the elastic side becomes increasingly sensitive to small angular perturbations, which could serve as an early indicator of imminent plastic response in granular or colloidal systems."],"forward_implications":["If the divergence is real, the high-pressure quasi-elastic phase and the intermediate anomalous phase are separated by a true critical point, not a smooth crossover.","The classical Lamé equation fails below $p_c$; the screened equation with $k^2\\mathbf{d}$ and Bessel-function solutions becomes the correct continuum description of the intermediate phase.","The critical pressure can be predicted from the competition between plastic-avalanche size $\\ell$ and screening length $\\kappa^{-1}$, and the scaling argument yields $p_c\\approx2.4$, consistent with the simulation estimates of 2.3 and 2.9.","The framework gives a three-dimensional analogue of the hexatic phase: the intermediate phase keeps radial long-range correlations while losing angular correlations, with both translational and chiral symmetries broken."],"supporting_citations":[{"why":"Supplies the anomalous elasticity theory, including the screened Lamé equation (5), that the paper's intermediate phase is built on.","marker":"[6]"},{"why":"Extends dipole screening to three-dimensional amorphous solids, giving the theoretical basis for the $k^2\\mathbf{d}$ term and the anomalous response.","marker":"[8]"},{"why":"Defines the quadrupolar plastic events (Eshelby inclusions) whose gradients act as the effective dipole charges $P_\\alpha=\\partial_\\beta Q_{\\alpha\\beta}$.","marker":"[15]"},{"why":"Provides the selection principle for the discrete values of the screening parameter $\\kappa$ that determine the observed anomalous solution.","marker":"[17]"},{"why":"Establishes the intermediate phase between jammed and unjammed amorphous solids that the present transition is located within.","marker":"[32]"},{"why":"Provides the generalized Lorentzian form used to fit $S(r,f)$ and to extract the diverging angular correlation length $\\theta^\\dagger=1/f^\\dagger$.","marker":"[33]"},{"why":"Supplies the geometric scaling of the unstable region that the paper uses to estimate the maximal avalanche size $\\ell\\sim p^{-1/3}$.","marker":"[35]"},{"why":"Gives the relation $\\Delta Z \\approx 6\\kappa$ connecting excess coordination to screening, which the paper uses to estimate the critical pressure $p_c\\approx2.4$.","marker":"[36]"}],"fun_headline_variants":["Dipole-induced transition observed in 3D amorphous solids","Pressure-tuned 3D transition with diverging angular order","3D dipole transition: critical scaling in disordered solids","Novel 3D dipole transition found in strained packings","Angular correlations diverge at 3D dipole transition"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"Everything rests on the assumption that the angular fluctuation spectrum is really a generalized Lorentzian with exponent 2.85 and that a two-dimensional disk relation between excess contacts and screening length carries over to three-dimensional spheres; if either is not exact, the extracted correlation length may not genuinely diverge.","fun_headline_variants_meta":{"raw":{"variants":["Dipole-induced transition observed in 3D amorphous solids","Pressure-tuned 3D transition with diverging angular order","3D dipole transition: critical scaling in disordered solids","Novel 3D dipole transition found in strained packings","Angular correlations diverge at 3D dipole transition"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000246,"raw_usage":{"total_tokens":1537,"prompt_tokens":943,"completion_tokens":594,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":559,"completion_tokens_details":{"reasoning_tokens":512}},"tokens_in":559,"tokens_out":594,"duration_ms":7072,"temperature":1.0,"reasoning_tokens":512,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T17:15:04.667304+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure $S(r,f)$ for the same packings at pressures just below $p_c$, but fit it without fixing the exponent 2.85; if the best-fit $f^\\dagger$ approaches a positive constant rather than zero, the claimed critical divergence is an artifact of the fitting form. A second check is to repeat the inflation experiment in larger systems and see whether the fitted exponent $\\mu$ stays near 1.66 and whether the divergence sharpens with system size, as a true critical point should.","supporting_citations":[{"cited_title":"Lema ˆ ıtre, C","cited_arxiv_id":null,"evidence_quote":"Supplies the anomalous elasticity theory, including the screened Lamé equation (5), that the paper's intermediate phase is built on."},{"cited_title":"Charan, M","cited_arxiv_id":null,"evidence_quote":"Extends dipole screening to three-dimensional amorphous solids, giving the theoretical basis for the $k^2\\mathbf{d}$ term and the anomalous response."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the quadrupolar plastic events (Eshelby inclusions) whose gradients act as the effective dipole charges $P_\\alpha=\\partial_\\beta Q_{\\alpha\\beta}$."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the selection principle for the discrete values of the screening parameter $\\kappa$ that determine the observed anomalous solution."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the intermediate phase between jammed and unjammed amorphous solids that the present transition is located within."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the generalized Lorentzian form used to fit $S(r,f)$ and to extract the diverging angular correlation length $\\theta^\\dagger=1/f^\\dagger$."},{"cited_title":"Wyart, S","cited_arxiv_id":null,"evidence_quote":"Supplies the geometric scaling of the unstable region that the paper uses to estimate the maximal avalanche size $\\ell\\sim p^{-1/3}$."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the relation $\\Delta Z \\approx 6\\kappa$ connecting excess coordination to screening, which the paper uses to estimate the critical pressure $p_c\\approx2.4$."}],"review_version":1}