{"id":"07a2c55c-db94-493b-a088-9da5d9d60062","arxiv_id":"2602.17873","paper_version":2,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":4,"one_line_summary":"Including higher-order correction terms and analyticity classes in the long-time spreadability fit recovers the spectral-density exponent α with errors down to ~10^-6 on exact model data.","lead":"This paper improves the standard way to extract the scaling exponent α of a two-phase material's structure from long-time diffusion spreadability data by adding higher-order corrections to the fit. The method yields much smaller errors than the previous leading-order fit and works for hyperuniform, typical nonhyperuniform, and antihyperuniform models, with relevance to NMR-based material characterization.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Type I fit assumes correction exponents β_i = i, while the general expansion (9) allows arbitrary β_i; non-integer corrections are untested and may bias α.","rationale":"The reader's weakest assumption concerns the general power-law expansion premise and the exclusion of non-expandable media. My concern is more specific: even within the power-expandable class, the fitting functions do not implement the general expansion (9) with arbitrary β_i. They only handle commensurate corrections (half-integer or integer t exponents, plus logarithms for integer α). This matters because the central claim is not merely that higher-order terms help for the three chosen models; it is that the Type I procedure is the most broadly applicable and yields high-precision α. Non-integer correction exponents are physically plausible (e.g., corrections to scaling at critical points, or class III hyperuniform media with non-trivial small-k structure), and the paper provides no test or explicit limitation covering them. The proposed concrete test is direct and inexpensive: since Eq. (1) is exact and any nonnegative integrable spectral density defines a valid autocovariance via Bochner's theorem, one can generate exact spreadability data with a non-integer correction and see whether Type I recovers α. If the test fails, the method's generality is overstated, though its value for the tested commensurate cases remains intact. I therefore recommend a conditional acceptance: either the authors add a validation for non-integer β_i or they explicitly narrow the scope claim. This is not a rejection; the exact-data demonstrations for α = 0, 2, -1, the noise robustness, and the two-point Padé construction are genuine strengths.","tokens_in":24539,"tokens_out":12308,"duration_ms":140615,"concrete_test":"Take a nonnegative, integrable spectral density in d = 3, \\tildeχ_V(k) = k^{1/2}(1 + c k^{0.7})/(1 + k^2)^4 with c = 0.5 (normalized so that χ_V(0) = φ1φ2). Generate exact s_ex(t) via Eq. (1) over the same 1000 log-spaced points t ∈ [100, 10000] used in Sec. V, with no noise, and apply the Type I fit (Eq. 24) for n = 0, ..., 8. If the optimal-order estimate of α deviates from 0.5 by more than ~10^-5, the restriction β_i = i is a genuine source of bias; if it recovers α to that precision, the concern is resolved. Repeat for a second non-integer β (e.g., β = 1.3) to confirm.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"Eq. (9) is the paper's general long-time expansion with arbitrary exponents β_i. The three fitting functions, however, all impose commensurate corrections: Type I (Eq. 24) assumes β_i = i, giving t^{-i/2} terms; Type II (Eq. 28) adds logarithms only for integer α; Type III (Eq. 30) assumes β_i = 2i. Thus a spectral density with, e.g., α = 1/2 and a correction k^{β} = k^{0.7} yields an s_ex(t) correction t^{-0.35}, which Type I cannot represent regardless of order n. The validations cover α = 0 (integer/even corrections), α = 2 (even corrections), and α = -1 (integer corrections plus logs); no test includes non-integer α or non-integer β_i. The paper states Type I is the 'weakest assumption' and applies to the broadest range (Sec. IV.A), and the central claim of high-precision α depends on that generality. For a real medium with non-integer correction exponents, the fitted α may be biased by the attempt to approximate t^{-β/2} with the t^{-i/2} basis; the claimed |δα| ~ 10^-5–10^-6 is not established for this allowed class. Section VI excludes only non-power-expandable spectral densities (stealthy hyperuniform, quasicrystals), so this gap is not acknowledged.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents an improved method for extracting the small-wavenumber exponent α from the long-time asymptotics of the normalized excess diffusion spreadability s_ex(t) of two-phase media. Starting from the exact Fourier representation, the authors derive higher-order asymptotic expansions (Eqs. 8–10) for s_ex(t) in terms of the small-k expansion of the spectral density, and propose three types of fitting functions: Type I with half-integer correction powers, Type II with logarithmic corrections for integer α, and Type III for analytic spectral densities with even integer α. They also construct a two-point Padé approximant for approximating s_ex(t) at all times. The fitting scheme is benchmarked on exact spreadability data for Debye random media (α=0), disordered hyperuniform media in d=2 and d=3 (α=2), and an antihyperuniform model (α=-1), including Gaussian-noise robustness tests. The reported noise-free exponent errors reach ~10^-6–10^-5, and noisy-data errors scale roughly as |δα|~η (Fig. 9).","tokens_in":24932,"tokens_out":9297,"duration_ms":94842,"significance":"If the procedure is accepted, it provides a concrete improvement over the earlier Wang–Torquato fitting algorithm, enabling more accurate microstructure characterization from spreadability data obtained by NMR or simulation. The derivation from the exact Fourier representation is clean, and the coefficient formula (10) is a direct, parameter-free consequence of that representation. The classification Table I is a useful synthesis, and the noise tests are carefully designed and reported. The central caveat is that the proposed fitting bases are highly commensurate; the claimed high precision is demonstrated only for media whose small-k expansions contain integer or half-integer correction powers. The paper would be strengthened by explicitly scoping this limitation or by adding tests for non-commensurate corrections.","major_comments":[{"comment":"Type I is described as making the weakest assumption, but Eq. (24) fixes β_i=i, whereas the general expansion (9) allows arbitrary β_i. Thus Type I can represent only corrections of the form t^{-i/2}; it cannot represent, for example, a spectral density k^{1/2}(B + B_{0.7} k^{0.7} + ...), whose correction enters s_ex(t) as t^{-0.35}. All validations in Sec. V use spectral densities with integer/even correction powers (DRM, DHM) or integer powers plus logarithms (AHM). No test involves non-integer β_i or non-integer α. The claimed |δα|≲10^-5–10^-6 is therefore not established for the broad class that Type I is asserted to cover. Please add at least one benchmark with non-commensurate correction exponents, or explicitly restrict the applicability claim in Sec. VI.","section":"§IV.A, Eq. (24) vs. Eq. (9)"},{"comment":"The optimal fitting order n_o is selected by a qualitative comparison of fits across orders, described as separating 'underfitting' from 'overfitting' regions. The headline accuracy values (Figs. 4–8) are reported at the n_o chosen by this heuristic. No objective selection criterion is given, so the procedure cannot be reproduced by a user without the authors' judgment. I recommend defining a quantitative rule (e.g., stability of α across consecutive n, an information criterion, or cross-validation) and reporting the sensitivity of α̂ to n in a neighborhood of n_o.","section":"§IV.B"},{"comment":"The synthetic benchmarks are closed-loop: exact s_ex(t) is computed from the same analytic χ_V(r) or spectral density whose small-k expansion is used as the fitting model. This is a valid test of the truncation procedure but not of its behavior when the true spectral density has a correction structure outside the fitted basis. This compounds the issue raised in the first major comment. A test on data generated from a model not of the fitted form (e.g., a numerically simulated microstructure or a spectral density with an ad hoc non-analytic term) would make the 'real data' claim in Sec. VI more convincing.","section":"§V"}],"minor_comments":[{"comment":"The two-point Padé parameters t0 and ν are not specified. Figure 2(c) states that the approximant involves only three independent coefficients (A1, C0, C1), but t0 and ν are also free parameters. Please clarify how these are chosen.","section":"§III.B, Eq. (20)"},{"comment":"The notation |kσ|^3 and |kσ|^4 is confusing; presumably |k|σ is meant. Please define σ and use consistent notation.","section":"§II.D, Eq. (14)"},{"comment":"There are several typographical errors, including 'synethetic' in the Introduction, 'efficient' in Sec. II, and 'as as' before Eq. (18). A careful proofreading pass would be helpful.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The paper is technically sound in its derivations and the examples are well chosen, but the scope of the 'weakest assumption' claim for Type I is overstated relative to what is tested. The fix is local: either add a non-commensurate correction benchmark or explicitly limit the central claim. The optimal-order selection issue is also likely to be raised by other readers and should be addressed in revision. I would support publication after these points are resolved."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read this one if you care about extracting α from spreadability data. The paper does what it says: it extends the Wang–Torquato fit by adding higher-order corrections, classifies the long-time expansion types, and gives a two-point Padé approximant that works across all t. The exact-model benchmarks are well chosen—α = 0, 2, and −1—and the error numbers (|δα| down to 10⁻⁵ or better, |δα| ~ noise in the noisy tests) are convincing for those models. The derivation from Eq. (1) and the Fourier transforms in Appendix A is standard and checkable. This is a genuine, useful extension of the spreadability toolkit, not a repackaging.\n\nSoft spots: the generality of Type I is overstated. Eq. (9) allows arbitrary correction exponents β_i, but Type I fixes β_i = i, i.e. t^{−i/2} corrections. The paper calls Type I the weakest assumption and says it applies to the broadest range, but that is only true within the subclass of spectral densities whose corrections occur at integer steps in k. A medium with, say, α = 1/2 and a k^{0.7} correction yields a t^{−0.35} term that no order of Type I can represent, and the fitted α can absorb the mis-specified correction. None of the benchmark models tests a non-integer α or non-integer β_i; the antihyperuniform model has α = −1 plus logs, but still integer corrections. The stress-test note has this right. The paper would be stronger with a class III hyperuniform example (non-integer α) and at least one synthetic spectral density with a non-integer correction exponent, to see whether the bias actually appears and whether the optimal-order heuristic catches it.\n\nOther smaller points: the optimal fitting order n_o is chosen by visual under/overfitting inspection; that should be acknowledged more candidly as a heuristic, not a criterion. There is no released code, and validation is synthetic-only. These are minor relative to the core claim.\n\nBottom line: for the tested classes the method works and the math holds up. The central claim—higher-order fits reduce α bias—is established for those classes; the scope of that claim needs tightening. This deserves a serious referee; I'd send it out and ask for a non-integer-α test case and a more careful statement of Type I's domain of validity.","headline":"Solid incremental methods paper: higher-order spreadability fits extract α more accurately on exact benchmarks, but the Type I generality claim is broader than what is actually tested.","tokens_in":25373,"tokens_out":3495,"would_cite":true,"duration_ms":38598,"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":"By fitting next-order terms in the long-time spreadability expansion, the paper recovers the microstructural scaling exponent α with errors near the noise level.","keywords":["diffusion spreadability","two-phase media","hyperuniformity","spectral density","long-time asymptotics","Padé approximant","scaling exponent","microstructure characterization"],"falsifier":"Take a disordered hyperuniform model with a known noninteger exponent in the range 0<α<1 (class III), generate its spreadability exactly from Eq. (1) via its spectral density, and apply the Type I fit with the order optimized; the fitted α should converge to the known value as the time window and fitting order are optimized — if it does not, the claimed bias removal is false.","tokens_in":24453,"feed_emoji":"🧪","tokens_out":8842,"duration_ms":74603,"temperature":0.7,"pith_summary":"The paper refines the way the long-time tail of the diffusion spreadability is fitted, so that the scaling exponent α that encodes how a two-phase material suppresses density fluctuations can be measured far more precisely than before. By adding higher-order correction terms to the asymptotic power-law fit, and by choosing those correction terms based on whether the material's spectral density is analytic at the origin, the scheme recovers α to within about 10^-5 on clean data and to within the noise level on noisy data. This matters because the spreadability is an experimentally accessible probe tied to NMR relaxation measurements, and α is the single number that separates hyperuniform from typical or antihyperuniform media. The paper also constructs a two-point Padé approximant that reproduces the full time decay of the spreadability with only a few parameters.","feed_headline":"Corrected fits recover the spreadability exponent to 10^-5 or better","feed_subtitle":"Adding correction terms reduces the bias in fitted alpha, making it a reliable microstructure classifier.","key_machinery":"The central object is the long-time asymptotic expansion of the normalized excess spreadability, s_ex(t) = t^{-(d+α)/2} Σ_{i=0}^∞ C_{β_i/2} t^{-β_i/2}, obtained by substituting a power-law expansion of the spectral density about k=0 into the exact integral representation of s_ex(t). The expansion's exponents β_i encode the analyticity of the spectral density at the origin, so choosing the right correction terms is what removes the bias in the fitted α. The three proposed fitting functions (Type I, II, III) are the practical embodiment of this machinery; the two-point Padé approximant blends the short- and long-time expansions to cover all times.","core_discovery":"The paper claims that the standard single-power-law fit to the long-time excess spreadability is systematically biased, and that fitting the full asymptotic expansion s_ex(t) = t^{-(d+α)/2} Σ C_{β_i/2} t^{-β_i/2} removes most of that bias. The correction exponents mirror the powers of k in the spectral density near the origin, so three fitting functions are needed: a general Type I, a Type II for integer α with logarithmic terms, and a Type III for analytic spectral densities. On Debye random media (α=0), disordered hyperuniform media in 2D and 3D (α=2), and an antihyperuniform medium (α=-1), the optimal-order fit recovers the exponent to within 10^-6, 10^-5, and 10^-2 respectively; with Gau","pith_inferences":["The same correction-term strategy should transfer to any Gaussian-smoothed transform of a two-point statistic, so analogous bias-corrected fits could sharpen exponent extraction from other scattering-derived quantities.","The three-type hierarchy suggests a model-selection protocol: run Type I first, then test whether coefficient-ratio criteria justify switching to Type II or III; this order-restricted logic could inform other asymptotic-fitting pipelines.","Because only two-point statistics enter, the high-precision α recovered here may serve as a sensitive, cheap discriminator for simulated microstructures, potentially flagging finite-size or noise artifacts that direct spectral-density fits would miss.","For noninteger class-III hyperuniform media (α = 1/2, say), Type I should in principle work; testing the procedure there would close a gap in the paper's numerical demonstrations."],"forward_implications":["For typical nonhyperuniform, hyperuniform, and antihyperuniform disordered media, the exponent α can now be read off from spreadability data with errors near 10^-5–10^-6 in the noise-free case and errors of order the noise level when noise is added.","The fitting output also reveals structural constraints: whether the spectral density is analytic at the origin, which moments of the autocovariance function exist, and how fast the autocovariance decays at large distances.","The two-point Padé approximant provides a compact all-time parameterization of s_ex(t) with a few coefficients, which the paper argues can support inverse design of microstructures with targeted spreadability.","For antihyperuniform media with integer negative α, including logarithmic correction terms (Type II) is necessary to capture the exact long-time expansion; the paper demonstrates this for α = -1."],"fun_headline_variants":["Corrected fits recover spreadability exponent to 1e-5","Bias-free spreadability exponent from improved fits","Higher-order terms sharpen spreadability exponent extraction","Refined fits nail spreadability exponent in two-phase media"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The method assumes the spectral density admits a small-k expansion in powers of k (possibly with logarithmic factors) leading with k^α; if a real material's spectral density instead vanishes in a finite window (stealthy hyperuniform) or is dense in Bragg peaks near the origin (quasicrystals), the fitted α no longer represents the infinite-wavelength scaling exponent.","fun_headline_variants_meta":{"raw":{"variants":["Corrected fits recover spreadability exponent to 1e-5","Bias-free spreadability exponent from improved fits","Higher-order terms sharpen spreadability exponent extraction","Refined fits nail spreadability exponent in two-phase media"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000595,"raw_usage":{"total_tokens":2739,"prompt_tokens":979,"completion_tokens":1760,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":723,"completion_tokens_details":{"reasoning_tokens":1705}},"tokens_in":723,"tokens_out":1760,"duration_ms":13211,"temperature":1.0,"reasoning_tokens":1705,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T22:06:45.668883+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a disordered hyperuniform model with a known noninteger exponent in the range 0<α<1 (class III), generate its spreadability exactly from Eq. (1) via its spectral density, and apply the Type I fit with the order optimized; the fitted α should converge to the known value as the time window and fitting order are optimized — if it does not, the claimed bias removal is false.","supporting_citations":[],"review_version":1}