{"id":"c2abc100-48cb-4a8c-a659-1005bd21938a","arxiv_id":"2504.16924","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Soft-core repulsions convert disordered stealthy hyperuniform ground states into ultradense, size-independent sphere packings, with packing fractions approaching those of jammed disk and MRJ sphere packings in the low-χ limit.","lead":"Adding a short-range soft-core repulsion to the standard stealthy-hyperuniform optimization lets the authors generate disordered sphere packings with packing fractions up to 0.86 in 2D and 0.63 in 3D, independent of system size. Without that repulsion, the maximal achievable packing fraction of such ground states falls to zero as the number of particles grows.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Reported φmax values are not maxima over all SHU ground states: BCC (d=3, φ≈0.680) and square-lattice (d=2, φ≈0.785) configurations satisfy both potentials with Φ=0, contradicting the claimed SAT-UNSAT boundary and the no-soft-core decay.","rationale":"The reader's weakest assumption, the independence approximation in Eq. (9), is a valid numerical concern about the no-soft-core extrapolation rate. However, the more fundamental issue is that φmax is not a maximum over all ground states: explicit lattice configurations satisfy both potentials exactly and have packing fractions above the reported values. This is not a disagreement with the numerical data; it is an interpretation issue in the central claim. The paper's title and abstract emphasize disordered packings, so the authors may intend to restrict to disordered ground states, but the text does not do so consistently: Sec. IV A states flatly that φmax decreases to zero, Sec. V calls Eq. (21) an upper bound on φmax, and Fig. 5 is labeled a SAT-UNSAT phase diagram of Eq. (17), all without the qualifier 'among configurations reached from random initial conditions.' Because BCC is a zero-energy ground state at φ=0.680, the phase diagram as drawn is not the satisfiability boundary of the potential. The numerical ultradense packings remain an interesting result, and the soft-core numbers are reproducible observations, so the appropriate disposition is still CONDITIONAL: the authors must redefine φmax or explicitly restrict claims to typical disordered ground states from their optimization protocol, and they should respond to the independent-identically-distributed approximation in Eq. (9) with a direct test of nearest-neighbor-distance correlations. I disagree with the reader only in that the load-bearing weakness is the conflation of typical and maximal, not the independence approximation, although both should be addressed.","tokens_in":26225,"tokens_out":19848,"duration_ms":191953,"concrete_test":"Take d=3, χ=0.0025 and 0.45, N=2m^3 in a cubic box at density ρ=1, and initialize L-BFGS on the BCC lattice with σ=1.091 (φ≈0.680). Confirm that the final configuration has Φ below the paper's threshold, rmin=σ, and S(k)<10^-20 for |k|≤K. Repeat for d=2 with a square lattice at φ=π/4 and for d=1 with an integer lattice. If these zero-energy ground states are found, the reported φmax values are not the SAT-UNSAT boundary of Eq. (17), and the paper must qualify all 'maximal' statements as maximal only over the disordered branch accessible from random initial conditions.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central results depend on treating φmax as a property of the ground-state manifold of potentials (2) and (17), but the paper actually measures typical values reached from random initial conditions. This distinction is load-bearing. For d=3 and ρ=1, the BCC lattice with nearest-neighbor distance σ=1.091 (packing fraction π√3/8≈0.680) has all pair distances at least σ and shortest reciprocal-lattice vector |G|=2π√2/a_c≈7.05. Since K≤5.43 for every χ≤0.45, S(k)=0 for all |k|≤K and the soft-core term vanishes, so BCC is an exact zero-energy ground state of both (2) and (17). Thus the true SAT-UNSAT threshold of (17) is at least 0.680 in d=3, above the reported φmax=0.63. Analogously, the square lattice (φ=π/4≈0.785) is a valid SHU ground state in d=2 for all χ<0.5, and the integer lattice gives φ=1 in d=1. The no-soft-core statement that φmax decreases to zero as N→∞ is therefore false if 'maximal' is taken literally; the maximum over the ground-state manifold is bounded below by these lattice values for every N. What the simulations actually show is that the packing fraction of typical disordered ground states sampled by L-BFGS from random starts tends to zero; the 'at least 5 successes out of 50 starts' criterion makes φmax an algorithm-dependent threshold, not a phase boundary of the potential. The abstract, Sec. IV A, and Fig. 5 must be rephrased to restrict φmax to the disordered, random-initial-condition ensemble.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a modified collective-coordinate optimization scheme that adds a soft-core repulsion to the standard stealthy pair potential [Eq. (17)] and uses it to generate sphere packings from disordered stealthy hyperuniform (SHU) point patterns. Using the minimum-distance distribution P(rmin;N) and nearest-neighbor distribution HP(r;N), the paper reports that without soft-core repulsions the maximal packing fraction phi_max decreases to zero on average as N→∞ for chi<0.5, whereas with soft-core repulsions phi_max is independent of N and reaches 1.00, 0.86, and 0.63 in the zero-chi limit, decreasing to 1.00, 0.67, and 0.47 at chi=0.45 for d=1,2,3. The paper also provides empirical Padé and Weibull formulas for phi_max, characterizes pair statistics, contact numbers, hyperuniformity order metrics, and spectral densities, and compares the small-chi packings to fast-compression jammed hard-particle packings and MRJ states.","tokens_in":26510,"tokens_out":7784,"duration_ms":70214,"significance":"If interpreted as a characterization of typical disordered SHU ground states obtained from random initial conditions, the numerical results are a useful contribution: they show that a non-compression, fixed-density optimization can produce nearly jammed disordered SHU packings, with the d=3 packings closely matching MRJ packings in pair statistics, contact number, and packing fraction. The tabulated data for phi_max, Z(r=sigma+), and Lambda, together with the spectral-density calculations, are concrete assets for the stealthy-hyperuniform two-phase materials community. However, the paper's literal claims that phi_max is the true maximum over the ground-state manifold and that the boundary in Fig. 5(a) is a SAT-UNSAT phase transition are not supported, because ordered lattices with higher packing fractions are exact ground states of the same potentials. The significance of the work therefore depends on reframing the claims as statements about the random-initial-condition ensemble, not about the potential's ground-state manifold.","major_comments":[{"comment":"The quantity phi_max is defined in Sec. III as the largest target packing fraction for which at least 5 of 50 L-BFGS runs from random initial conditions reach the energy tolerance Phi<7chi*10^-20. This is an algorithm-dependent threshold for a particular optimization protocol, not the satisfiability threshold (SAT-UNSAT boundary) of the potential (17). Exact zero-energy ground states with larger packing fractions exist: in d=3 at rho=1, the BCC lattice with nearest-neighbor distance 1.091 has S(k)=0 for all |k|<=K (the shortest reciprocal vector is about 7.05, while K<=5.43 for chi<=0.45) and rmin>=sigma for the sigma corresponding to every reported phi_max, so it is a ground state of both (2) and (17) with packing fraction pi*sqrt(3)/8=0.680, exceeding the reported phi_max=0.63 (chi->0) and 0.47 (chi=0.45). In d=2, the square lattice with phi=pi/4=0.785 is likewise a ground state for chi<=0.45, exceeding the reported phi_max=0.67 at chi=0.45. Thus the statements that beyond phi_max 'the ground state ceases to exist' and that Fig. 5(a) separates satisfiable and unsatisfiable phases are not correct as statements about the potential; the abstract, Sec. IV.A, Sec. V, and Fig. 5 must be rephrased to refer to typical disordered ground states obtained from random initial conditions.","section":"Sec. III, Sec. II.D, Fig. 5(a)"},{"comment":"The thermodynamic-limit conclusion for the no-soft-core case rests on the approximation that the N nearest-neighbor distances r(i) are independent and identically distributed, with HP(r;N) replaced by HP(r;infinity) and an empirical correction factor gamma=1/2 inserted in Eq. (9). This is an unproven assumption about correlations among nearest-neighbor distances. If the r(i) are positively correlated, the decay of the minimum distance with N could be slower than predicted, or could saturate at a positive value; negative correlations would make the decay faster. Because this approximation directly underlies the probabilistic upper bound in Eq. (21), the estimate in Eq. (22), and the Weibull extrapolation in Eq. (23), the claim that phi_max tends to zero for chi<0.5 needs either a rigorous bound or an explicit restriction to the random-initial-condition ensemble, with the accuracy of Eq. (9) quantified in the main text rather than only asserted via the supplementary material.","section":"Sec. IV.A, Eqs. (8)-(9), (21)-(23)"},{"comment":"There is an internal inconsistency about the chi threshold for the no-soft-core decay: Sec. IV.A states that phi_max decreases to zero in the thermodynamic limit for chi<0.5 in d=1 but only for chi<=0.35 in d=2 and chi<=0.30 in d=3, whereas the abstract and conclusions claim the decay for all chi<0.5. This discrepancy is not cosmetic: Fig. 4 shows a qualitative change near chi=0.4, and the Weibull extrapolation in Eq. (23) is used precisely for the difficult regime 0.35<chi<0.5. The authors should specify the actual threshold as a function of d and reconcile the abstract and main text.","section":"Sec. IV.A and Abstract"}],"minor_comments":[{"comment":"The Padé approximant in Eq. (24) is fitted to the same simulation data shown in Fig. 5(a), so it is an interpolation formula for the random-start threshold, not an independently predictive law or a physical phase boundary; this should be stated explicitly.","section":"Eq. (24) and Table II"},{"comment":"Equation (25) is introduced as a theoretical upper bound on the mean contact number, but the text and Fig. 5(b) treat it as an equality. Please clarify that equality corresponds to isostatic saturation and that the simulation data approach the bound from below.","section":"Eq. (25) and Fig. 5(b)"},{"comment":"The phrase 'phi_max decreases to zero on average' is ambiguous: if phi_max denotes the maximum over all ground states, the statement is false because lattice ground states exist; if it denotes a typical value or an algorithm-dependent threshold, the ensemble and optimization protocol should be part of the definition.","section":"Abstract and Sec. III"},{"comment":"The description of the chi=0.0025 d=2 states as 'highly ordered structures with large triangular coordination domains' is in tension with calling them disordered SHU packings; a quantitative measure of orientational order or a clearer statement about the presence of defects would help avoid confusion.","section":"Sec. VI"}],"recommendation":"major_revision","confidential_remarks":"The BCC and square-lattice counterexamples are decisive against the literal reading of phi_max as a SAT-UNSAT threshold of the potential, and the no-soft-core claim of decay to zero needs to be restricted to the random-initial-condition ensemble. Nevertheless, the numerical data and structural analyses are likely to be of genuine value to the stealthy-hyperuniform community once the claims are reframed. I recommend major revision rather than rejection, provided the authors explicitly narrow their claims, quantify the independence approximation in Eq. (9), and reconcile the chi-threshold inconsistency."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe paper's headline is overstated. The reported φmax values are not maxima over the SHU ground-state manifold. In d=3, the BCC lattice with nearest-neighbor distance σ≈1.091 (φ≈0.680) is an exact zero-energy ground state of both potentials (2) and (17) for every χ≤0.45, since its shortest reciprocal vector is ≈7.05, well above K for those χ. In d=2, the square lattice (φ≈0.785) is an analogous ground state for all χ≤0.45. The integer lattice gives φ=1 in d=1. So the true SAT-UNSAT boundary of the potential is at least 0.680 in 3D and 0.785 in 2D, above the reported 0.63 in 3D and 0.67 in 2D at χ=0.45. The no-soft-core claim that φmax decreases to zero as N→∞ is false if 'maximal' is taken literally. What the simulations actually measure is the largest packing fraction that L-BFGS from random starts can achieve with a 10% success criterion (Sec. III). That is an algorithm-dependent threshold, not a phase boundary of the potential. The abstract, Sec. IV A, and Fig. 5 need to be rephrased to restrict φmax to the disordered, random-initial-condition ensemble.\n\nCredit where due: the numerical study is careful, and the central observation—that soft-core SHU ground states reach φ≈0.86 in 2D and 0.63 in 3D with N-independent statistics—is well supported. The structural comparison to jammed and MRJ states (Fig. 8) is convincing, and the χ-dependence of contact number and the emergence of polymer-like chains is new and interesting. The spectral density results are useful for applications.\n\nSoft spots beyond the main framing: Eq. (9) treats the N nearest-neighbor distances as i.i.d. with an empirical γ=1/2; that uncontrolled approximation underlies the no-soft-core thermodynamic-limit claim. The Padé and Weibull formulas are fits to the same data, so calling them explicit formulas is generous. No code or data are deposited, and the MRJ connection leans on a companion paper.\n\nWho it's for: researchers working on hyperuniform two-phase media, disordered packings, and jamming. The numerical route to ultradense disordered SHU packings is worth knowing about. The paper deserves peer review, but the referee should insist on a major revision: either redefine φmax as a typical threshold for the random-start protocol, or prove that the ordered lattice states are excluded under a precise definition of 'disordered.' With that fix, it would be a solid contribution.\n\nRecommendation: send to peer review with the expectation of heavy revision.","headline":"The φmax values in the abstract are algorithm-dependent thresholds, not true maxima: BCC and square-lattice ground states are exact zero-energy states with higher density.","tokens_in":27196,"tokens_out":10485,"would_cite":false,"duration_ms":86935,"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":"Adding a soft-core repulsion to the collective-coordinate energy turns disordered stealthy hyperuniform point patterns into sphere packings with maximal packing fractions of 1.00, 0.86, and 0.63 in dimensions 1, 2, and 3.","keywords":["stealthy hyperuniform","sphere packings","collective-coordinate optimization","soft-core repulsion","packing fraction","nearest-neighbor statistics","maximally random jammed packings","disordered materials"],"falsifier":"Numerically measure the minimum-distance distribution P(rmin;N) for standard-potential SHU ground states in d = 2 and d = 3 at χ = 0.3 for N up to $10^{5}$ or $10^{6}$ and check whether the mean minimum distance decays according to the Weibull extrapolation (23) or instead saturates at a positive plateau; the prediction fails if a positive plateau appears.","tokens_in":25860,"feed_emoji":"⚪","tokens_out":7144,"duration_ms":59152,"temperature":0.7,"pith_summary":"Disordered stealthy hyperuniform (SHU) point patterns, whose structure factor vanishes over a finite band of wavenumbers, can be converted into sphere packings of unexpectedly high density by adding a short-range soft-core repulsion to the collective-coordinate optimization energy. The paper argues that this modification makes the maximal packing fraction independent of system size, reaching 1.00, 0.86, and 0.63 in the zero-stealthiness limit for dimensions d = 1, 2, 3, and declining to 1.00, 0.67, and 0.47 at χ = 0.45. Without the repulsion, an extreme-value analysis of minimum-pair and nearest-neighbor distances predicts that the maximal packing fraction decreases to zero on average as the particle number N grows for χ < 0.5. The result matters because it opens a route to disordered hyperuniform two-phase materials at densities previously thought inaccessible, with the small-χ packings closely resembling jammed hard-particle packings.","feed_headline":"Soft-core repulsion pushes stealthy hyperuniform packings to 86%","feed_subtitle":"Maximal packing fractions become system-size independent, reaching 1.00, 0.86, 0.63 in 1D, 2D, and 3D.","key_machinery":"The load-bearing object is the modified collective-coordinate potential (17), Φ(r^N) = (ρ/2) Σ_{k≠0} ṽ(k)S(k) + Σ_{i<j} u(r_ij), with ṽ(k)/v0 = Θ(K − |k|) and the soft-core repulsion u(r)/ε0 = (1 − r/σ)^2 Θ(σ − r). Because both sums are non-negative, any zero-energy ground state must simultaneously satisfy S(k) = 0 for |k| < K and rmin ≥ σ, so the point pattern maps directly to a packing of nonoverlapping spheres of diameter σ with packing fraction ρv1(σ/2). The extreme-value connection (8)-(9), which approximates the minimum-distance distribution through the nearest-neighbor distribution HP(r;∞) under an independence assumption with empirical correction γ = 1/2, supplies the no-soft-core upper bounds and the Weibull extrapolation.","core_discovery":"The paper establishes that the maximal packing fraction of a sphere packing derived from a stealthy hyperuniform ground-state point pattern is governed by the minimum pair distance of that pattern. For ground states of the standard stealthy potential with no soft-core repulsion, the nearest-neighbor statistics imply that as N increases the minimum distance shrinks, so that the maximal packing fraction approaches zero on average in the thermodynamic limit for χ below 1/2. With the soft-core repulsion added, every ground state must have all pair separations at least σ, so the packing fraction ρv1(σ/2) can be raised to a well-defined maximum that does not depend on N. The reported maxima are φmax = 1.00, 0.86, 0.63 for d = 1, 2, 3 in the zero-χ limit, decreasing to 1.00, 0.67, and 0.47 at χ = 0.45, and the mean contact number follows the isostaticity bound Z(σ+) = 2d(1 − 2χ).","pith_inferences":["If the independence approximation in Eq. (9) survives scrutiny, the no-soft-core decay of φmax should follow a Weibull extreme-value law; a direct numerical measurement of P(rmin;N) at N = 10^5 for d = 2, 3 would settle this.","The soft-core construction with σ chosen larger than the particle diameter should produce fully connected-matrix microstructures, and the accompanying transport predictions (effective permittivity, permeability, survival time) can be computed directly from the spectral densities the paper reports.","The near-identity of the zero-χ packings with jammed MRJ states suggests that the degenerate SHU ground-state manifold may contain jamming-relevant configurations in all dimensions, which a study of contact-network statistics as N → ∞ could test."],"forward_implications":["In two and three dimensions the soft-core SHU packings achieve φmax = 0.86 and 0.63 at small χ while preserving exact stealthy hyperuniformity, S(k) = 0 for |k| ≤ K, enabling high-density photonic, acoustic, and transport applications.","For χ → 0 the packings become effectively jammed and isostatic, with contact numbers Z(σ+) = 4.42 (2D) and 5.91 (3D), so the optimization provides a non-compression route to such states.","The mean contact number follows Z(σ+) = 2d(1 − 2χ), so increasing χ progressively fragments the packings into polymer-like chains with fewer contacts and lower φmax.","The computed spectral densities of the packings allow quantitative estimates of effective dynamic dielectric response, fluid permeability, and mean survival time of the resulting two-phase dispersions.","The one-dimensional packings are integer lattices with φmax = 1.00 for all χ < 1/2, so the construction saturates the density bound in d = 1."],"supporting_citations":[{"why":"Defines the stealthiness parameter χ and the degree-of-freedom count dc = d(1 − 2χ) that underpin the contact-number upper bound Z(σ+) = 2d(1 − 2χ).","marker":"[13]"},{"why":"Supplies the previous range of achievable packing fractions for SHU packings without soft-core repulsions, providing the baseline that the paper's ultradense packings are contrasted with.","marker":"[21]"},{"why":"Provides the nearest-neighbor distribution functions HP(r;N) for SHU ground states that the extreme-value argument in Eqs. (8)-(9) builds on.","marker":"[56]"},{"why":"Introduced the collective-coordinate optimization scheme and established the disorder and degeneracy of SHU ground states in two dimensions.","marker":"[11]"},{"why":"Introduced the modified collective-coordinate procedure with soft-core repulsion, which the present paper scales to maximal packing fractions.","marker":"[22]"},{"why":"Provides the 3D maximally random jammed packing reference data, including pair statistics and the gap exponent, used to identify the χ → 0 soft-core states as MRJ-like.","marker":"[60]"},{"why":"Companion work that pins down the three-dimensional χ → 0 limit of φmax and the gap exponent γ ≈ 0.44, supporting the paper's claim that the degenerate manifold contains jammed states.","marker":"[61]"},{"why":"Supplies the harmonic contact potential and fast-compression jammed packings used as structural benchmarks for the ultradense soft-core SHU states.","marker":"[68]"}],"fun_headline_variants":["Soft-core repulsion boosts hyperuniform packings to 86% density","Hyperuniform sphere packings reach 86% with soft cores","Disordered hyperuniform packings: soft cores push density to 86%","Maximal packing fraction hits 86% in soft-core hyperuniform systems"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The claim that, without soft-core repulsions, the maximal packing fraction vanishes in the thermodynamic limit rests on treating the N nearest-neighbor distances as independent random variables drawn from the infinite-N distribution, corrected by an empirical factor γ = 1/2; if correlations among these distances become significant at large N, the predicted decay to zero could be wrong in rate or in existence.","fun_headline_variants_meta":{"raw":{"variants":["Soft-core repulsion boosts hyperuniform packings to 86% density","Hyperuniform sphere packings reach 86% with soft cores","Disordered hyperuniform packings: soft cores push density to 86%","Maximal packing fraction hits 86% in soft-core hyperuniform systems"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000248,"raw_usage":{"total_tokens":1671,"prompt_tokens":1197,"completion_tokens":474,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":813,"completion_tokens_details":{"reasoning_tokens":405}},"tokens_in":813,"tokens_out":474,"duration_ms":4905,"temperature":1.0,"reasoning_tokens":405,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T10:52:24.398106+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Numerically measure the minimum-distance distribution P(rmin;N) for standard-potential SHU ground states in d = 2 and d = 3 at χ = 0.3 for N up to $10^{5}$ or $10^{6}$ and check whether the mean minimum distance decays according to the Weibull extrapolation (23) or instead saturates at a positive plateau; the prediction fails if a positive plateau appears.","supporting_citations":[{"cited_title":"Oppenheimer , author D","cited_arxiv_id":null,"evidence_quote":"Defines the stealthiness parameter χ and the degree-of-freedom count dc = d(1 − 2χ) that underpin the contact-number upper bound Z(σ+) = 2d(1 − 2χ)."},{"cited_title":"Aeby , author G","cited_arxiv_id":null,"evidence_quote":"Supplies the previous range of achievable packing fractions for SHU packings without soft-core repulsions, providing the baseline that the paper's ultradense packings are contrasted with."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Companion work that pins down the three-dimensional χ → 0 limit of φmax and the gap exponent γ ≈ 0.44, supporting the paper's claim that the degenerate manifold contains jammed states."}],"review_version":1}