{"id":"b1283d94-a68e-4be4-8caa-ad5ce2366685","arxiv_id":"1908.06776","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"For one-dimensional systems with nearest-neighbor interactions, the pair connectedness can be computed exactly from the pair correlation function by solving a Volterra integral equation.","lead":"This paper derives an exact integral equation that links the probability that two particles belong to the same connected cluster to the pair correlation function in one-dimensional systems. It offers a practical route from structural data to connectivity for these systems, with a sketch of how the idea extends to three dimensions.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The printed kernel equations are internally inconsistent: Eq. (18) and Eq. (27) use a different normalization for omega-prime than Eqs. (22)-(25) and Eq. (33), so the central derivation cannot be followed as written until the convention is fixed.","rationale":"The reader's weakest assumption, the factorization Eq. (19), is defensible: for one-dimensional nearest-neighbor systems the density hierarchy is Markov, and Eq. (19) is the standard renewal factorization, so I do not regard that as a live objection. My concern is instead that the manuscript as printed does not define omega' consistently. Eq. (18) treats it as a density, while Eqs. (22)-(25) and (33) treat it as the dimensionless ratio h; Eq. (27) belongs to neither convention. This is a concrete, checkable inconsistency in the central derivation. It does not show the final integral equation is wrong, since the paper reproduces the Domb and Drory solutions and Monte Carlo data, so I would keep the reader's CONDITIONAL verdict rather than reject. A revision should fix the normalization convention and correct Eq. (27); after that the central claim appears sound.","tokens_in":17375,"tokens_out":30197,"duration_ms":319650,"concrete_test":"Set g(r)=1 (ideal gas) in the printed formulas and compare: Eq. (18) gives omega'=rho e^{-rho r}; Eq. (22) has solution e^{-rho r}; Eq. (27) with K=omega' gives K=rho - integral_0^r K, whose solution is rho e^{-r}; Eq. (33) states e^{-rho r}. Recompute the pair connectedness for rho=2, d=1 from each candidate kernel in Eq. (24) and compare with Domb's result and Eq. (36). Only one convention reproduces the known solution; adopt that convention consistently in Eqs. (18), (22), (27), and (28), then rerun Figs. 5-8 to verify the numerical claims.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim rests on extracting a nearest-neighbor kernel from g and using it in Eq. (24). The manuscript's own equations for this kernel do not have one consistent meaning for omega'. Eq. (18) defines omega' as a number density, omega' = rho g * P(empty), which for the ideal gas would be rho e^{-rho r}. Eq. (22) is only consistent with the later solution Eq. (33), omega'=e^{-rho r}, and with the renewal identity g = omega' + rho omega' * g, if omega' is the dimensionless ratio h = (nearest-neighbor probability density)/rho. Eq. (27) then prints K(r)=rho g(r)-integral K(r-x)g(x) dx, which for g=1 gives K'=-K and K=rho e^{-r} (or rho e^{-rho r} under a different convention), not the e^{-rho r} used in Eq. (33). The factor rho is on the wrong term. Because Eq. (27) is the stated bridge from g to the kernel in Eqs. (25)-(28), a reader implementing the method from the printed formulas cannot reproduce the Domb or Drory results; the internal inconsistency must be resolved before the exactness claim is certified. This is not a question of the underlying physics, which is likely correct, but of the derivation as written.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript develops an integral-equation formalism for continuum percolation in which the pair connectedness P is expressed through the nearest-neighbor distribution. For one-dimensional systems with nearest-neighbor interactions, the authors further relate the nearest-neighbor distribution to the pair correlation function g via a Volterra equation (Eqs. (22)–(23)), yielding a closed equation (Eq. (24)) for P in terms of g alone. The approach is shown to reproduce the two known exact one-dimensional solutions, for ideal fully penetrable rods (Domb) and for hard rods (Drory), and the paper compares favorably with Monte Carlo simulations for a repulsive nearest-neighbor potential and for a Lennard-Jones fluid. Extensions to external fields, long-ranged interactions, and a three-dimensional ideal gas with a geometric closure are outlined.","tokens_in":17718,"tokens_out":18975,"duration_ms":178280,"significance":"If the one-dimensional results are correct, the paper provides a unified and computationally simple route from the pair correlation function to the pair connectedness for arbitrary nearest-neighbor interactions, avoiding system-specific derivations. The reproduction of both known exact solutions is a strong check, and the Monte Carlo comparisons add credible numerical support. A notable strength is that the kernel in the generalized integral equation is built from a physically observable conditional probability, which in principle can be sampled in simulations—an advantage over the direct connectivity in the Ornstein–Zernike approach. However, the central exactness claim is currently obscured by an internal normalization inconsistency in the printed equations, which must be resolved before readers can implement the method.","major_comments":[{"comment":"The normalization of ω' is inconsistent between Eq. (18) and the equations used to obtain the kernel. Eq. (18) defines ω'(τ)=ω(τ)P(empty|particle at τ)=ρg(τ)P(...), giving ω' units of inverse length, whereas Eqs. (22)–(25) and the ideal-gas solution Eq. (33) require ω' to be dimensionless, namely the ratio of the nearest-neighbor probability density to ρ. As a consequence, Eq. (27), which reads K(r)=ρg(r)-∫ K(r-x)g(x)dx, is dimensionally inconsistent and does not follow from Eq. (22). With ω' dimensionless, the correct relation is K(r)=g(r)-ρ∫_0^r K(r-x)g(x)dx. The authors should adopt one convention throughout (for instance ω' = (nearest-neighbor probability density)/ρ) and correct Eqs. (18) and (27). As printed, a reader implementing the method cannot reproduce the Domb or Drory solutions.","section":"Section III, Eqs. (18)–(27)"},{"comment":"The claimed equivalence of Eq. (25) to Coniglio's expansion (Eq. (8)) is established through an informal cancellation-diagram argument. In particular, the treatment of configurations with multiple outlying particles is summarized by 'we can repeat this procedure,' and the final conclusion that 'both expansions can indeed be brought in perfect unison' is asserted rather than demonstrated. Since exactness of the one-dimensional integral equation is the central claim, the authors should either supply a complete proof, for example a probabilistic derivation of Eq. (22) from the factorization Eq. (19), or explicitly present the diagrammatic equivalence as a consistency check supported by the two exact benchmarks rather than as a proof.","section":"Section III, after Eq. (28)"}],"minor_comments":[{"comment":"The phrase 'how the formalism can applied' should be 'can be applied'.","section":"Abstract"},{"comment":"The integral is printed as ρ∫_r^d ... but for r>d the limits should be ρ∫_d^r ...; the lower limit appears to be a typo.","section":"Eq. (35)"},{"comment":"The sentence 'at larger, a slight discrepancy' is incomplete; it should read 'at larger r' or similar.","section":"Figure 10 and surrounding text"},{"comment":"The three-dimensional closure is presented as a 'purely geometrical treatment,' but the proportionality constant is fixed only by the limiting condition Eq. (61); the approximation is uncontrolled. The manuscript is appropriately cautious elsewhere, but this phrase overstates the systematic character of the choice.","section":"Section IV C, Eq. (63)"},{"comment":"The notation ρ(1) is used both for the density profile in the presence of the external field and for the bulk density; this should be clarified to avoid confusion.","section":"Section IV A, Eq. (45)"}],"recommendation":"major_revision","confidential_remarks":"The normalization inconsistency in Eqs. (18) and (27) is a technical error that appears fixable without changing the underlying physics, which is strongly supported by the reproduction of known exact solutions and by the Monte Carlo comparisons. The informal diagrammatic equivalence argument is a further gap that needs tightening. With these addressed, the manuscript would be a worthwhile contribution to the journal."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear [Colleague],\n\nThe paper has a solid idea: in one dimension with nearest-neighbor interactions, the pair connectedness P can be obtained from the pair correlation function g alone through a Volterra equation. That claim is real and it unifies the two known exact solutions (Domb's ideal gas, Drory's hard rods), which is a strong check. The Monte Carlo comparisons in figures 6-9 look good, and the paper is honest about where exactness stops.\n\nBut the manuscript as printed has an internal inconsistency that has to be fixed before the derivation is usable. Equation (18) defines ω' as a number density, while equations (22)-(25) and the worked examples require ω' to be the dimensionless reduced nearest-neighbor probability density. Concretely, for the ideal gas, Eq. (22) gives ω' = e^{-ρr}, but Eq. (27) as printed, with K set to ω', gives K = ρ e^{-r}. The paper even attributes the correct exponential to Eq. (27) when it actually follows from Eq. (22). This is not a cosmetic issue: a reader implementing the method from the printed formulas cannot reproduce Domb's or Drory's results. The underlying physics is almost certainly correct, because the derivation from Eq. (19) to Eq. (24) is coherent once the normalization is fixed. But the exactness claim does not hold for the equations as written.\n\nTwo smaller issues. The equivalence to Coniglio's diagrammatic expansion is argued through representative diagrams and cancellation arguments, not a full proof for all orders. And the simulation comparisons have no error bars, so the agreement is only visual.\n\nWho is this for? Researchers in continuum percolation and liquid-state theory who want a benchmark or a way to turn g(r) data into connectivity predictions in one dimension. It deserves a serious referee after the normalization issue is resolved; the idea is useful enough that I would cite it once the equations are consistent.\n\nRecommendation: send to peer review, with a clear request to fix the ω' convention and either add error bars or soften the visual-agreement claim.","headline":"The one-dimensional relation is real and useful, but the printed equations mix two normalizations of ω', so the derivation cannot be followed as written.","tokens_in":18155,"tokens_out":9093,"would_cite":false,"duration_ms":84000,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["82B43","82B21","45D05"],"pacs":["64.60.ah","05.20.-y"],"model":"deepseek-v4-flash","headline":"The paper claims that in one-dimensional nearest-neighbor fluids, the pair connectedness follows exactly from the pair correlation function through a Volterra integral equation.","keywords":["continuum percolation","pair connectedness","Volterra integral equation","nearest-neighbor distribution","pair correlation function","one-dimensional fluids","hard rods","ideal gas"],"falsifier":"Take a one-dimensional fluid with a strictly nearest-neighbor pair potential such as the inverse-square potential $V_\\varepsilon$ used in the paper, measure $g(r)$ in simulation, build $\\omega'$ from Eq. (22), solve Eq. (25), and compare the predicted pair connectedness with the simulated $P(r)$. If the two differ by more than statistical error at any density, the claimed exactness fails.","tokens_in":17161,"feed_emoji":"🔗","tokens_out":6854,"duration_ms":67192,"temperature":0.7,"pith_summary":"This paper tries to show that the connectivity problem of continuum percolation can be solved from thermal pair structure, not from the full many-body correlation hierarchy. For one-dimensional fluids with nearest-neighbor interactions, it derives an exact Volterra integral equation that gives the pair connectedness $P$ from only the pair correlation function $g$ and the number density. If the claim is right, cluster connectivity for this class of systems is no harder than solving a one-dimensional integral equation, and the known exact solutions become special cases of one unified procedure. The authors demonstrate this by reproducing the only two analytic results available, the one-dimensional ideal gas and hard rods, and by matching Monte Carlo simulations for nearest-neighbor and short-ranged pair potentials. They also outline how the same equation extends to external fields, long-range forces, and higher dimensions at the cost of a closure for a conditional probability.","feed_headline":"One integral equation gives exact connectivity in 1D fluids","feed_subtitle":"Nearest-neighbor statistics plus a Volterra equation yields exact cluster connectivity in one dimension.","key_machinery":"The load-bearing object is a Volterra integral equation of the second kind whose kernel is the nearest-neighbor distribution $\\omega'$. In one dimension, ordered $n$-point correlation functions factorize into products of pair functions, so $\\omega'$ is a functional of $g$ alone and the connectivity diagrams reduce to ordered chains that the Volterra resolvent sums exactly. In the generalized setting, the same chain carcass remains but the kernel becomes the constrained ratio $p(0|\\tau,r)/p(0)$, replacing exactness by a closure problem for a geometrically interpretable observable.","core_discovery":"The central discovery is an exact integral equation for the pair connectedness $P(0,r)$, the density of particles connected to a tagged particle at the origin by overlapping connectivity shells of diameter $d$. For one-dimensional nearest-neighbor systems the equation is $$P(0,r)=\\Theta(d-r)g(0,r)+\\Theta(r-d)\\rho\\int_0^d d\\tau\\,\\omega'(0,\\tau)P(\\tau,r),$$ where $g$ is the pair correlation function and $\\omega'$ is the nearest-neighbor distribution. Because ordered correlations factorize for these systems, $\\omega'$ is itself determined by $g$ through another Volterra equation, so $P$ follows from $g$ alone without any closure. The paper shows the chain diagrams generated by this equation reproduce the standard connectivity diagrammatic expansion, and it uses the equation to recover the previously known analytic solutions for the one-dimensional ideal gas and hard rods. In higher dimensions and for long-ranged potentials the same equation holds with a constrained conditional probability in the kernel, which must then be approximated or sampled.","pith_inferences":["Because $P$ is a functional of $g$ in one dimension, approximate liquid-state closures could be benchmarked by how well their predicted $g$ reproduces a simulated pair connectedness, giving a stricter test than structure alone.","The geometric closure used for the three-dimensional ideal gas suggests a testable route: replace the conditional probability $c(r,t,u)$ with a directly measured quantity, converting the Fredholm equation into an exact numerical procedure for arbitrary three-dimensional fluids.","The one-dimensional invertibility implies that percolation observables such as cluster-size distributions should also be derivable from $g$ by iterating the same Volterra chain, although the paper does not carry out that construction."],"forward_implications":["For one-dimensional nearest-neighbor fluids, the pair connectedness can be computed exactly from $g$ without closures, so any analytic or simulated pair structure yields connectivity predictions.","The known analytic pair connectedness for the one-dimensional ideal gas and hard rods is recovered as a special case, unifying two previously unrelated derivations.","External fields do not break the exact scheme; they only make the kernel position-dependent through the one-particle density.","For long-range interactions and higher dimensions, the equation still holds but with a conditional probability in the kernel that can be sampled in simulation or approximated geometrically, as illustrated for the three-dimensional ideal gas.","The relation is invertible: pair connectedness determines the nearest-neighbor kernel, so in these systems $g$ and $P$ carry the same information."],"supporting_citations":[{"why":"It defines the pair connectedness and gives the connectivity Ornstein-Zernike equation that the paper's diagrammatic comparison starts from.","marker":"[9]"},{"why":"It supplies the exact one-dimensional ideal-gas pair connectedness used as a benchmark.","marker":"[30]"},{"why":"They provide the exact impenetrable hard-rod solutions that Eq. (24) must recover.","marker":"[31, 32]"},{"why":"It supplies the Volterra resolvent-kernel machinery used to solve the integral equations.","marker":"[35]"},{"why":"It establishes how nearest-neighbor distributions relate to correlation functions and why they generally cannot be inferred from $g$ alone.","marker":"[37]"},{"why":"It provides the factorization of ordered correlations in one dimension that makes $\\omega'$ a functional of $g$.","marker":"[38]"},{"why":"It gives the exact nearest-neighbor gap distribution for hard rods, used as the kernel input.","marker":"[39]"},{"why":"It provides the analytic hard-rod pair distribution from which the paper constructs $g$.","marker":"[41]"}],"fun_headline_variants":["Exact connectivity in 1D from one integral equation","Single equation yields exact cluster statistics in 1D","Nearest-neighbor distributions unlock exact percolation","Integral equation recovers exact 1D fluid connectivity","Percolation from density: exact 1D solution"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that, for ordered one-dimensional configurations, higher-order correlations factor into products of two-particle correlations; if this factorization fails, as it does for interactions beyond nearest neighbors, the nearest-neighbor distribution is not determined by $g$ alone and Eq. (24) is no longer exact.","fun_headline_variants_meta":{"raw":{"variants":["Exact connectivity in 1D from one integral equation","Single equation yields exact cluster statistics in 1D","Nearest-neighbor distributions unlock exact percolation","Integral equation recovers exact 1D fluid connectivity","Percolation from density: exact 1D solution"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000195,"raw_usage":{"total_tokens":1344,"prompt_tokens":922,"completion_tokens":422,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":538,"completion_tokens_details":{"reasoning_tokens":343}},"tokens_in":538,"tokens_out":422,"duration_ms":4252,"temperature":1.0,"reasoning_tokens":343,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:35:18.027002+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a one-dimensional fluid with a strictly nearest-neighbor pair potential such as the inverse-square potential $V_\\varepsilon$ used in the paper, measure $g(r)$ in simulation, build $\\omega'$ from Eq. (22), solve Eq. (25), and compare the predicted pair connectedness with the simulated $P(r)$. If the two differ by more than statistical error at any density, the claimed exactness fails.","supporting_citations":[{"cited_title":"Coniglio, U","cited_arxiv_id":null,"evidence_quote":"It defines the pair connectedness and gives the connectivity Ornstein-Zernike equation that the paper's diagrammatic comparison starts from."},{"cited_title":"Domb, in Mathematical Proceedings of the Cambridge Philosophical Society , Vol","cited_arxiv_id":null,"evidence_quote":"It supplies the exact one-dimensional ideal-gas pair connectedness used as a benchmark."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It supplies the Volterra resolvent-kernel machinery used to solve the integral equations."},{"cited_title":"Torquato, B","cited_arxiv_id":null,"evidence_quote":"It establishes how nearest-neighbor distributions relate to correlation functions and why they generally cannot be inferred from $g$ alone."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It provides the factorization of ordered correlations in one dimension that makes $\\omega'$ a functional of $g$."},{"cited_title":"Zernike and J","cited_arxiv_id":null,"evidence_quote":"It gives the exact nearest-neighbor gap distribution for hard rods, used as the kernel input."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It provides the analytic hard-rod pair distribution from which the paper constructs $g$."}],"review_version":1}