{"id":"e541f311-9664-4d9f-abdb-b4409d710ee0","arxiv_id":"1908.02053","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"For infinite-horizon Lorentz gas and stadium channel, the spreading density is a Lambert-corrected Gaussian core with power-law corridor tails; the renewal Lévy walk works for the Lorentz gas but not for the stadium channel.","lead":"This paper derives analytical formulas for how a cloud of particles spreads between bouncing off circular or stadium-shaped walls in two billiard models with open straight corridors. It shows the spreading develops fat power-law tails along corridors, so the usual Gaussian picture only appears after absurdly long times.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 'no fitting' claim for the Lorentz gas rests on the empirically chosen Lambert split η=1, which controls the O(1/Ω) Kummer correction; its sensitivity is untested.","rationale":"I read the paper as claiming that the geometrically derived inter-collision CDF, inserted into a Lambert-scaled Lévy walk, gives a parameter-free description of the Lorentz gas packet. The CDF derivation itself is strongly supported: it matches the simulated CDF, reproduces the known large-τ tail, and yields ⟨τ⟩ close to the rigorous value. The weakest step is the closure of the Lévy-walk approximation, where the paper must choose a velocity model F_d(v) and a logarithmic split η. The paper explicitly acknowledges in Eq. (60) that η cannot be determined uniquely and that η=1 is chosen because it gives good results. Since the O(1/Ω) Kummer correction is the paper's main non-Gaussian signature, and Ω is only logarithmic in t, the no-fit comparison has an unexamined degree of freedom at the same order as the effect being highlighted. The η-sensitivity test I propose would settle whether this is a real vulnerability or a harmless convention. I do not find a critical internal inconsistency, and the stadium-channel section is honestly presented as an effective-parameter description rather than a parameter-free prediction. The reader's CONDITIONAL verdict remains appropriate, so I recommend no change.","tokens_in":30085,"tokens_out":16822,"duration_ms":189012,"concrete_test":"Take the R=0.4 and R=0.3 cases of Figs. 9 and 10 and recompute the generalized version of Eq. (26) obtained from the split in Eq. (60) with η=1/2 and η=2, keeping all billiard-derived constants (τ0, ⟨τ⟩, Cψ, q) fixed. Then compute the fractional difference max_x |P2_η(x,0) − P2_η=1(x,0)| / P2_η=1(x,0) in the corridor tail region 10^2 ≤ x ≤ 10^3, and similarly for y=x. If the simulated histogram lies inside the η∈{1/2,1,2} band, the concern is resolved and the choice is not load-bearing. If the data track the theory only at η=1, the no-fit claim is conditional on this unconstrained constant.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central Lorentz-gas claim is that Eq. (26) reproduces the simulated density \"without any fitting.\" The weakest load-bearing point is the Lambert split parameter η introduced in Eq. (60). There, the logarithm in the characteristic function is split at an arbitrary point: ln[2C_f^2 k^2/(NΩ)] is written as ln[2C_f^2 η/(NΩ)] + ln(k^2/η), and the first term is used to fix Ω while the second is kept as the O(1/Ω) correction. The paper chooses η=1 because it \"produces good results\"; this is not derived from the billiard dynamics. The same O(1/Ω) correction is precisely the Kummer-type term that distinguishes the theory from Bleher's Gaussian and produces the visible corridor tails. For t=10^4, Ω is only of order 20, so a change of η by a factor of e changes the correction by roughly 1/Ω, i.e. a few percent, comparable in size to the effect being displayed. Thus the agreement in Figs. 2, 3, 9, and 10 could depend on this unconstrained convention even though no billiard parameter is fitted. This is not an accusation of curve fitting; it is an untested sensitivity of the parameter-free claim.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper studies finite-time transport in two deterministic billiard models with infinite horizons: the Lorentz gas on a square lattice with circular scatterers (two or four open corridors) and the stadium channel. The authors derive, from geometry and the assumption of uniform collision parameters, the cumulative distribution function of inter-collision times for each model (Eqs. (28), (33)-(47), (49), (55)-(58)), and the CDFs match single-trajectory simulations. They then model the spatial spreading with a Lévy walk whose velocity distribution is discrete (Eqs. (19), (24)), using a Lambert-W scaling to handle the logarithmic divergence at the border between Gaussian and Lévy statistics. The main claim is that for the Lorentz gas the position PDF (Eq. (26)) reproduces simulations without fitting, with constants τ0, ⟨τ⟩, Cψ, and q computed from the analytic ψ(τ), while for the stadium channel Eq. (23) describes the data only with two effective fitted parameters because the renewal assumption fails.","tokens_in":30367,"tokens_out":8108,"duration_ms":83130,"significance":"If the parameter-free claim held in full, the paper would provide a practically relevant mesoscopic propagator for the infinite-horizon Lorentz gas, bridging the regime where Bleher's Gaussian limit is not yet valid, and would explain the cross-like and British-flag shapes and the power-law corridor tails. The paper has real strengths: the geometric CDF formulas are tested against numerical histograms; the Lambert scaling is a useful resummation of the slow logarithmic factors; the correlation plots in Fig. 12 provide a direct diagnostic of the renewal assumption; and the authors are explicit that the stadium-channel case requires effective parameters and that renewal fails there. The main caveat is that the 'without any fitting' claim for the Lorentz gas is qualified by an unconstrained Lambert split parameter η and by the modeling ansatz for the velocity distribution.","major_comments":[{"comment":"The Lambert split parameter η introduced in Eq. (60) is a free parameter, and the paper chooses η=1 because it 'produces good results' rather than from a dynamical or asymptotic condition. This is load-bearing because the O(1/Ω) term in Eqs. (10), (23), and (25), namely the Kummer correction, is precisely the term that produces the non-Gaussian tails emphasized in the paper. For the R=0.3 Lorentz gas and t=10^4, Eq. (21) gives Ω(t)≈8.2, so replacing η by eη changes the subleading term by roughly 12%; this is the same order as the difference between the theory and Bleher's Gaussian shown in Figs. 9 and 10. Consequently, the statement that Figs. 2, 3, 9, and 10 reproduce the simulation 'without any fitting' is not literally supported. Please either derive η from the asymptotic expansion, demonstrate insensitivity of the conclusions to η, or explicitly reclassify η as an empirically fixed resummation constant.","section":"Section V, Eq. (60)"},{"comment":"The velocity distribution F_2(v) used for the Lorentz gas is an effective ansatz rather than a consequence of the collision dynamics: the actual post-collision velocity is continuous and correlated with the incoming flight, while Eq. (24) replaces it by an IID draw restricted to the corridor directions, with the diagonal weight q taken from Eq. (47). The excellent agreement with particle simulations provides empirical support, but the paper does not test how much of the match depends on this particular ansatz, nor why the continuous components should be irrelevant. Because this choice is part of the claimed parameter-free mapping to the Lévy walk, the text should state explicitly that Eq. (24) is an effective distribution and should include at least one robustness check, such as an alternative q or a velocity distribution with continuous angular support.","section":"Section III, Eq. (24)"},{"comment":"The central claim of parameter-free reproduction of the Lorentz-gas density is supported only by visual log-density maps and two cross-sections. Since the claim is quantitative and the η ambiguity affects exactly the displayed non-Gaussian correction, the paper needs a quantitative comparison metric, such as relative L1 error or residual plots as a function of r at fixed t, to substantiate 'reproduces well' and 'without any fitting'. This is particularly important because the two cross-sections shown are along symmetry directions, while the full 2D comparison in Figs. 2 and 3 is only qualitative.","section":"Section V / Figs. 2, 3, 9, 10"}],"minor_comments":[{"comment":"The headings 'LAMBER T SCALING' and 'LÉVY W ALK' contain spacing errors; they should read 'Lambert scaling' and 'Lévy walk'.","section":"Section II/III headings"},{"comment":"The text uses 'Lorenz gas' in the discussion of correlation patterns; this should be 'Lorentz gas'.","section":"Section V"},{"comment":"The condition for four open horizons is stated as '√20≤ R < 1/√8'; from the context and the later use of 1/√20, this should be 1/√20 ≤ R < 1/√8.","section":"Section IV A"},{"comment":"Reference [17] has a typographical error: 'H. K. ZhangCommun. Math. Phys.' should read 'H. K. Zhang, Commun. Math. Phys.'","section":"Reference [17]"},{"comment":"The uniformity of a and α in Eq. (49) is assumed rather than derived; the excellent CDF match in Fig. 11 supports the assumption, but a brief justification or citation would strengthen the presentation.","section":"Section IV B"},{"comment":"The statement that Eq. (23) 'can indeed describe the stadium channel model' is stronger than the evidence, since the constants in Fig. 4 are obtained by a two-parameter fit; 'can be fitted to' or 'can effectively describe' would be more accurate.","section":"Section V / Fig. 4"}],"recommendation":"major_revision","confidential_remarks":"The main obstacle is the role of η in the Lorentz-gas 'no fitting' claim; if the authors can either derive η, show robustness of the displayed results to η, or clearly label it as an effective resummation parameter, the paper would be publishable. The requested quantitative comparison is standard for a no-fitting claim and would significantly strengthen the paper."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know about this one: it is a serious analytical attempt at the infinite-horizon Lorentz gas and stadium channel, and the Lorentz gas part largely works. The authors derive the geometric CDF of inter-collision times, plug it into a Lévy-walk propagator with Lambert scaling, and reproduce the full 2D density without fitting the final distribution. That is the real news. The CDF construction via shadowing inequalities is original and checked against numerics; the tau0 expressions agree with known limits, and the mean collision times match the rigorous formula to within 0.02%. Credit where due.\n\nThe stadium channel is weaker by design. The renewal assumption fails, the paper says so plainly, and the authors fall back on a two-parameter fit to the t = 10^4 PDF. That is honest, but it means the stadium claim is an effective description, not a prediction.\n\nThe soft spot I would flag is the Lambert split parameter eta in Eq. (60). The logarithm is split at an arbitrary point, and eta = 1 is chosen because it produces good results. Since Omega(t) is only about 20 at t = 10^4, changing eta by a factor of e shifts the O(1/Omega) Kummer correction by a few percent, which is comparable to the effect being displayed. So the \"without any fitting\" claim for the Lorentz gas has a small asterisk: the outcome is not fully parameter-free. I do not think this is curve fitting—the same eta works for the IID sum and both billiards, which points to a convention rather than a fudge—but the sensitivity is untested. A referee should ask for a robustness check, varying eta and showing the cross-sections stay within the numerical error.\n\nThe stadium-channel effective parameters are stable across time, which is nice, but they were not derived from the CDF. That is a second, clearly acknowledged limitation.\n\nOverall, I buy the Lorentz gas result. The math is consistent with rigorous limits, the numerical match is convincing, and the candid reporting of the stadium failure raises my confidence in the rest. This deserves peer review, not desk rejection. I would tell the editor to send it out.","headline":"A serious analytical push on the infinite-horizon Lorentz gas that mostly lands: geometric CDFs plus a Lambert-scaled Lévy walk reproduce the 2D density without fitting, though an arbitrary split constant and the stadium-channel effective fit keep it short of airtight.","tokens_in":30935,"tokens_out":1727,"would_cite":true,"duration_ms":18332,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["05.60.-k","05.40.Fb"],"model":"deepseek-v4-flash","headline":"The paper claims that for infinite-horizon Lorentz gases the full two-dimensional particle density follows from the geometric inter-collision time distribution, without fitting.","keywords":["infinite-horizon Lorentz gas","Lévy walk","Lambert scaling","inter-collision time distribution","billiard transport","stadium channel","power-law tails","Gaussian-Lévy border"],"falsifier":"A direct test: compute the particle density for the Lorentz gas at other scatterer radii within the two-corridor and four-corridor regimes, for example $R=0.35$ or $R=0.28$, and compare with Eq. (26) using only the geometry-derived constants; a mismatch in the corridor plume tails would show that the parameter-free reproduction is specific to the two sampled radii. A second check is a correlation plot of consecutive inter-collision times for $\\tau$ values above the plateau region, where the paper's renewal diagnosis is read directly from the scatter.","tokens_in":29849,"feed_emoji":"🎱","tokens_out":7703,"duration_ms":74236,"temperature":0.7,"pith_summary":"The paper sets out to predict how a packet of particles spreads in two deterministic billiard systems—the infinite-horizon Lorentz gas and the stadium channel—using only the geometry of the scatterers. It derives exact cumulative distributions of the time between collisions, whose tails decay as $\\psi(\\tau)\\sim \\tau_0^2/\\tau^3$, and inserts them into a Lévy-walk model rescaled with the Lambert $W$-function. For the Lorentz gas with two or four open corridors, the resulting position density reproduces numerical simulations without fitting, including the cross-like and British-flag-like plumes along the corridors. The same construction fails for the stadium channel, where consecutive flights are strongly correlated, although an effective-parameter version of the formula still tracks the density.","feed_headline":"Geometry alone reproduces the Lorentz gas particle density","feed_subtitle":"Inter-collision times drawn from the scatterer layout feed a Lévy-walk density that matches simulation without fitting.","key_machinery":"The central object is the inter-collision time distribution $\\psi(\\tau)$, computed as a sum over all reachable scatterers of integrals over the collision impact parameter $b$ and recoil angle $\\beta$, with integration boundaries fixed by no-blocking conditions. Its tail $\\psi(\\tau)\\sim \\tau_0^2/\\tau^3$ places the model at the border between Gaussian and Lévy statistics; Lambert scaling via the secondary branch $W_{-1}$ absorbs the logarithmically growing scale $\\Omega_d(t)$ so finite-time results can be compared with theory. A discrete velocity distribution $F_2(\\mathbf{v})$ with diagonal weight $q$ converts the geometry into the Lévy-walk (Montroll-Weiss) propagator.","core_discovery":"The central claim is that the geometry of an infinite-horizon Lorentz gas determines its finite-time transport statistics: the inter-collision time distribution $\\psi(\\tau)$ computed from the scatterer layout, combined with a discrete corridor velocity distribution and Lambert scaling, yields the full two-dimensional position density $P_2(\\mathbf{r},t)$ that matches simulation 'without any fitting' for $R=0.3$ and $R=0.4$. The power-law corridor tails are part of this density, not a separate asymptotic correction, while the rigorous Gaussian limit exists only after an astronomically large number of collisions. For the stadium channel the paper claims a weaker statement: the Lévy-walk formula can still describe the density if effective waiting-time parameters are fitted, even though the renewal assumption itself fails there.","pith_inferences":["If the parameter-free match persists at other scatterer radii, the same geometry-to-CDF recipe should transfer to other periodic lattices, where corridor directions and blocking conditions change but the $\\tau^{-3}$ tail index suggests a common universality class.","The paper's correlation plots of consecutive inter-collision times could be turned into a quantitative diagnostic: a measure of dependence in $(\\tau_n,\\tau_{n+1})$ would predict in advance which billiards admit a renewal description.","The Lambert split parameter $\\eta$ is the one empirical element in an otherwise derived chain; matching higher-order terms in the asymptotic expansion could fix it and remove the remaining freedom."],"forward_implications":["For the infinite-horizon Lorentz gas, finite-time transport is describable by a Lévy walk with geometry-derived $\\psi(\\tau)$, and the Gaussian propagator is not observable on practical timescales.","The Lambert-scaled formula supplies explicit finite-time corrections, including Kummer function terms, so simulations at times like $t=10^4$ can be compared with theory instead of waiting for $\\sqrt{N\\ln N}$ convergence.","The packet's shape—cross-like for two corridors, British-flag-like for four—is encoded in the velocity distribution $F_2(\\mathbf{v})$ and the parameter $q$ that measures diagonal corridor weight.","For the stadium channel, an effective-parameter Lévy walk reproduces the density after fitting two constants even when the renewal assumption is invalid.","The non-analytic, plateau-rich structure of the inter-collision time CDF reflects the periodic scatterer array and is inherited by the spreading density."],"supporting_citations":[{"why":"Supplies the large-$\\tau$ tail behavior and the uniform collision-parameter assumptions used in the geometry calculation.","marker":"[2]"},{"why":"Establishes the Gaussian limiting distribution that the finite-time theory is compared against.","marker":"[12]"},{"why":"Previous exact result for $\\tau_0$ and the cross-like packet, extended here to four corridors and the channel.","marker":"[23]"},{"why":"Provides the Lévy-walk model framework used to turn $\\psi(\\tau)$ into a spatial density.","marker":"[25]"},{"why":"Supplies the velocity-model expansion of the Montroll-Weiss equation used in the Lambert scaling derivation.","marker":"[27]"},{"why":"Gives the exact Fourier-Laplace propagator (Montroll-Weiss equation) from which the approximate densities are derived.","marker":"[31]"},{"why":"Supports renewal for the Lorentz gas through stretched-exponential decay of temporal correlations.","marker":"[16]"},{"why":"Provides the polynomial correlation decay bound that explains the stadium channel's failure of the renewal assumption.","marker":"[17]"}],"fun_headline_variants":["Geometry alone sets Lorentz gas density, no fitting","Scatterer layout yields Lévy density in Lorentz gas","Finite-time Lorentz gas: geometry beats Gaussian limit","Lévy walks from geometry, no fit for Lorentz gas","From corridors to plumes: exact density without fitting"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is the renewal step: after each collision the next flight duration is drawn independently from the geometry-derived distribution, so deterministic correlations between flights can be dropped; the stadium-channel results show this premise is model-dependent and can fail.","fun_headline_variants_meta":{"raw":{"variants":["Geometry alone sets Lorentz gas density, no fitting","Scatterer layout yields Lévy density in Lorentz gas","Finite-time Lorentz gas: geometry beats Gaussian limit","Lévy walks from geometry, no fit for Lorentz gas","From corridors to plumes: exact density without fitting"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000214,"raw_usage":{"total_tokens":1389,"prompt_tokens":870,"completion_tokens":519,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":486,"completion_tokens_details":{"reasoning_tokens":441}},"tokens_in":486,"tokens_out":519,"duration_ms":5946,"temperature":1.0,"reasoning_tokens":441,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:56:19.430311+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct test: compute the particle density for the Lorentz gas at other scatterer radii within the two-corridor and four-corridor regimes, for example $R=0.35$ or $R=0.28$, and compare with Eq. (26) using only the geometry-derived constants; a mismatch in the corridor plume tails would show that the parameter-free reproduction is specific to the two sampled radii. A second check is a correlation plot of consecutive inter-collision times for $\\tau$ values above the plateau region, where the paper's renewal diagnosis is read directly from the scatter.","supporting_citations":[{"cited_title":"In order to find the leading behavior of the second term of Eq","cited_arxiv_id":null,"evidence_quote":"Supplies the large-$\\tau$ tail behavior and the uniform collision-parameter assumptions used in the geometry calculation."},{"cited_title":"The pair{a,α} and the trio{x0,y 0,α} are related by a simple transformation","cited_arxiv_id":null,"evidence_quote":"Establishes the Gaussian limiting distribution that the finite-time theory is compared against."},{"cited_title":"Boldrighini, L","cited_arxiv_id":null,"evidence_quote":"Provides the Lévy-walk model framework used to turn $\\psi(\\tau)$ into a spatial density."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the velocity-model expansion of the Montroll-Weiss equation used in the Lambert scaling derivation."},{"cited_title":"Chernov,J","cited_arxiv_id":null,"evidence_quote":"Gives the exact Fourier-Laplace propagator (Montroll-Weiss equation) from which the approximate densities are derived."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supports renewal for the Lorentz gas through stretched-exponential decay of temporal correlations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the polynomial correlation decay bound that explains the stadium channel's failure of the renewal assumption."}],"review_version":1}