{"id":"12720839-fd1b-4fa7-940a-c9a2e7727223","arxiv_id":"2411.13254","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"In the thermal hyperbolic Calogero fluid, the Lax-matrix eigenvalue density is computed numerically and matches a generalized log-gas, with weak boundary-condition dependence in the thermodynamic limit.","lead":"Using Monte Carlo simulations, the authors map out the eigenvalue density of the Lax matrix for the hyperbolic Calogero gas in thermal equilibrium. The results support a bridge between random-matrix ensembles and integrable many-body systems, a useful ingredient for generalized hydrodynamics.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"High-density boundary-condition independence is not established: at rho>=5 the cosh-confined system is far from the thermodynamic limit (Fig. 6), and the C-c/C-r agreement in Fig. 5 may be a shared cosh-confinement finite-size artifact rather than the true N->infinity DOS.","rationale":"The reader's weakest assumption identifies the same load-bearing gap I find: the thermodynamic limit is inferred from N<=256 simulations even though the authors' own data (Figs. 5 and 6) show slow convergence at high density. I agree with that assessment. The paper's strongest independent evidence is the finite-N verification of the exact joint law Eq. (2.5) against direct Lax diagonalization at N=64 (Fig. 2), and the low/intermediate-density boundary-condition comparison in Fig. 3. Neither test reaches the high-density thermodynamic limit. At rho=5-11, Fig. 6 shows the cosh trap has not flattened on the particle scale, so the cosh-confined system is not yet equivalent to a box or ring of the same density. The observed agreement between C-c and C-r in Fig. 5 is therefore not independent support for the rational-Calogero limit: both systems share the same cosh confinement and the same finite-size bias. The paper explicitly acknowledges that neither has reached the thermodynamic limit, so there is no internal inconsistency; but the central claim of boundary-condition universality and the high-density approximation remains numerically open in precisely the regime where the paper claims them. A concrete resolution is to solve the variational problem Eq. (2.8) for the hyperbolic Calogero and to push the cosh-confined simulations to larger N at fixed rho. Since the reader already formulated this as the basis for CONDITIONAL, my stress-test does not change the verdict.","tokens_in":19949,"tokens_out":11169,"duration_ms":107681,"concrete_test":"Run the cosh-confined hyperbolic Calogero Monte Carlo at N=1024 (and 2048 if feasible) for rho=5 and rho=11, T=1, with ell=N/rho, and compare the resulting DOS to (i) the box and ring results at the same N and ell, and (ii) the numerical minimizer of the variational functional Eq. (2.8). If the large-N cosh DOS converges to the box/ring curve and to the Eq. (2.8) minimizer, the concern is resolved; if it remains distinct, the claimed boundary-condition independence at high density fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that the thermal Lax DOS is independent of boundary conditions in the thermodynamic limit is supported numerically only where the cosh-confined system has effectively reached a large flat box. At high densities (rho=5, 8, 11, T=1), Fig. 6 shows the cosh-confined particle density profile is still far from flat at N=256, so the N->infinity limit at fixed rho is not reached. Consequently, the comparison in Fig. 5 between cosh-confined Calogero (C-c) and cosh-confined rational Calogero (C-r) cannot validate the high-density approximation: both share the same cosh trap and hence the same finite-size bias. The only converged benchmark, box/ring rational Calogero versus the trigonometric TBA (Eqs. C.4-C.5), is not approached by C-c or C-r. Thus the paper's contribution (iii) — that at high density the Calogero DOS is well approximated by the rational Calogero model — and the general boundary-condition-independence claim are unverified in this regime. This matches the paper's own caveats on p. 15 but still leaves the weakest point of the central claim open.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the eigenvalue density (DOS) of the Lax matrix of the hyperbolic Calogero fluid in thermal equilibrium, using Monte Carlo sampling of positions and momenta followed by direct diagonalization. It compares three confining mechanisms: a cosh trap, a hard box, and a ring (the latter via the elliptic Calogero model), and benchmarks the results against the modified-log-gas joint distribution Eq. (2.5) for the cosh case, the exact low-density Toda-chain DOS, and the TBA solution of the trigonometric Calogero model. The central claims are that the thermal Lax DOS has a deterministic thermodynamic limit, that this limit is independent of boundary conditions, and that the low- and high-density limits are well approximated by the Toda chain and the rational Calogero model respectively.","tokens_in":20175,"tokens_out":6694,"duration_ms":72573,"significance":"If established, the paper would provide a useful bridge between integrable many-body theory and random matrix theory: the thermal Lax DOS is a basic input for generalized hydrodynamics, and the existence of an exact modified-log-gas representation is analytically valuable. The manuscript has genuine strengths: the agreement in Fig. 2 between direct diagonalization and direct Monte Carlo sampling of the eigenvalue joint distribution is an independent numerical check of the mapping in Eq. (2.5), and the low-density comparison with the Toda-chain DOS in Fig. 4 is convincing. The main weakness is that the high-density boundary-independence claim is not closed by the presented simulations, and the paper candidly acknowledges this in the final part of §4. The manuscript is therefore a promising contribution whose central claims are only partially verified.","major_comments":[{"comment":"The high-density boundary-independence claim is not established. In Fig. 5 the comparison is between cosh-confined Calogero (C-c) and cosh-confined rational Calogero (C-r), both of which are subject to the same cosh trap; the paper itself states that neither has reached the thermodynamic limit, and Fig. 6 shows that the cosh-confined particle density profile is still far from flat at N=256 for ρ=5, 8 and 11. The good C-c/C-r agreement at these densities may therefore be a shared finite-trap artifact rather than evidence for the N→∞ DOS, and the only converged benchmark (rational Calogero in a box, which agrees with the trigonometric TBA) is not approached by the displayed finite-size curves. A direct high-density simulation of box-confined Calogero fluid, or an explicit finite-size extrapolation of the C-c DOS at fixed ρ, is needed to support contribution (iii) and the boundary-independence statement.","section":"§4, Figs. 5 and 6"},{"comment":"The numerical verification of the modified-log-gas distribution Eq. (2.5) is performed only at N=64 and is reported without error bars or a discrepancy measure. Since Eq. (2.5) is the central analytical mapping used in the paper (and is taken from Ref. [44]), the claim of excellent agreement should be substantiated by, for example, the N-dependence of the integrated difference between the two methods or a chi-square per bin. Without such a quantitative check, the possibility of a finite-size coincidence is not excluded.","section":"§2 and Fig. 2"},{"comment":"The statement that the Lax DOS becomes independent of boundary conditions with increasing N is supported by only two system sizes in the figure (N=128 and 256), and visible differences remain at ρ=5 (Fig. 3c). The text mentions N=512 but no N=512 data appear in the figure, and no statistical uncertainties are shown. A finite-size scaling analysis, or at least a third and larger system size with reported error bars, would be necessary to close the thermodynamic-limit claim.","section":"§3, Fig. 3"}],"minor_comments":[{"comment":"There are numerous typos, including 'sill define' in §3, 'tignometric Calogeoro' in the Fig. 1 caption, 'po' in the Boltzmann weight definition in §3, and 'Heavyside' in Appendix C. A careful proofread is needed.","section":"Throughout"},{"comment":"The Monte Carlo protocol is not described: no equilibration time, proposal distribution, acceptance rate, or number of effectively independent samples is given. Adding these details, and reporting error bars or confidence bands in all figures, would make the numerical results reproducible.","section":"§3"},{"comment":"The truncation of the Weierstrass potential at n=±2 is introduced without a convergence test. Please quote the size of the omitted terms for the smallest ℓ used in the simulations, or compare the results with the n=±3 truncation.","section":"Appendix A, Eq. (A.28)"},{"comment":"The relation between the pressure P and the fixed density used in the main text is terse; the inversion of Eq. (B.5) should be stated explicitly so that the Toda benchmark is reproducible from the data given.","section":"Appendix B, Eq. (B.5)"}],"recommendation":"major_revision","confidential_remarks":"Ref. [44] is a book by one of the co-authors and is the source of the modified-log-gas formula and the TBA benchmarks. This is not improper, but it means the paper's own numerical verification carries additional weight. The high-density finite-size issue is the main risk; if the authors can add box-confined high-density Calogero data or a convincing finite-size extrapolation, I would support acceptance."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read the paper. Short version: solid, useful numerical study, with the expected caveats. The genuinely new content is the systematic Monte Carlo phase diagram of the Lax DOS for the hyperbolic Calogero fluid across densities, temperatures, and three boundary conditions (cosh trap, hard box, ring), together with an independent finite-N check of the modified log-gas mapping from Spohn's book [44]. The check in Fig. 2, comparing direct diagonalization of L against MC sampling of the eigenvalue joint distribution, is the strongest part of the paper: it confirms the mapping without relying on it. The low-density comparison with the Toda chain and the high-density comparison with the rational/trigonometric Calogero model are sensible benchmarks, and the agreement is convincing where the system has converged.\n\nThe soft spots are real but the authors do not hide them. At the highest densities (rho = 5, 8, 11, T = 1), the cosh-confined systems have not reached the thermodynamic limit at N = 256, as the paper's own Fig. 6 shows: the particle density profile is still not flat. So the boundary-condition independence claim is established at low and intermediate densities, but not at high density. The close match between cosh-confined Calogero and cosh-confined rational Calogero in Fig. 5 may be a shared finite-size artifact of the trap rather than the true thermodynamic-limit DOS. The stress-test note overstates one detail: the paper does show a trend toward the TBA curve as N grows, so the C-c and C-r curves are not static, but they remain far enough at N = 256 that the high-density approximation is not quantitatively closed.\n\nSeparately, the numerics are under-reported: no error bars, no thermalization or autocorrelation diagnostics, and no code. For a numerical paper in this day, that is a legitimate referee request. Novelty is moderate -- the central mapping comes from a coauthor's prior work -- but the phase diagram and boundary-condition comparison are new, and the paper is transparent about its own caveats.\n\nWho this is for: people working on generalized hydrodynamics or on Lax-matrix random-matrix connections who need the thermal Lax DOS of the Calogero fluid as input. It deserves a serious referee. I would send it out, ask for error bars, MC details, code, and either larger-N data at high density or a softened conclusion about boundary-condition independence. My own verdict is conditionally positive: the core numerical content holds together, but the high-density boundary-independence claim should not be taken as established.","headline":"Solid numerical study of the thermal Lax DOS for the hyperbolic Calogero fluid, with an honest but not fully closed boundary-independence claim at high density.","tokens_in":20715,"tokens_out":7935,"would_cite":true,"duration_ms":74639,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["60B20","37J35","82C22"],"pacs":[],"model":"deepseek-v4-flash","headline":"The hyperbolic Calogero fluid's Lax-matrix eigenvalue density has a boundary-independent thermodynamic limit given by a modified log-gas.","keywords":["Lax matrix","Calogero fluid","density of states","log-gas","generalized hydrodynamics","Toda chain","thermodynamic Bethe ansatz","random matrix theory"],"falsifier":"Simulate the high-density $\\cosh$-confined Calogero fluid (for example $\\bar\\rho = 11$, $T=1$) at $N = 512$ and $N = 1024$ and compare the Lax DOS and the scaled particle density profile with the box- and ring-confined results; if the curves continue to shift with $N$ or the density profile fails to flatten, the claimed boundary-condition independence has not been demonstrated.","tokens_in":19714,"feed_emoji":"🧮","tokens_out":11956,"duration_ms":106616,"temperature":0.7,"pith_summary":"This paper tries to establish that the eigenvalue density of the Lax matrix of the classical hyperbolic Calogero fluid, sampled from the thermal Gibbs distribution, becomes a deterministic, boundary-independent function in the thermodynamic limit. For the analytically convenient cosh external potential, it argues the joint distribution of eigenvalues is exactly a modified log-gas with a Bessel-K interaction term, Eq. (2.5), whose large-N free energy functional, Eq. (2.8), has a unique minimizer equal to the Lax density of states. Monte Carlo simulations up to N=256 are presented to show this modified log-gas matches direct diagonalization of the Lax matrix, and that the same DOS is obtained in cosh, box, and ring geometries. The paper further claims the low-density limit recovers the Toda-chain DOS and the high-density limit recovers the rational and trigonometric Calogero DOS, benchmarked against the thermodynamic Bethe ansatz. A sympathetic reader would care because the Lax DOS is the central input for generalized hydrodynamics of these integrable fluids, and a boundary-independent deterministic limit makes it a well-defined hydrodynamic ingredient.","feed_headline":"Lax-matrix spectrum of Calogero fluid is boundary-independent","feed_subtitle":"The same deterministic density appears in cosh, box, and ring traps, matching known integrable limits.","key_machinery":"The central object is the Lax matrix $L_c$ of the hyperbolic Calogero fluid: diagonal entries are the momenta $p_j$ and off-diagonal entries are $i/(2\\sinh((q_i-q_j)/2))$, so that $H_c = \\tfrac12\\operatorname{tr}(L_c^2)$. The argument is carried by the canonical map from $(q,p)$ to scattering coordinates $(\\lambda,\\phi)$, under which the $\\cosh$ trap becomes $\\sum_i e^{-\\ell/2} Y_i \\cosh\\phi_i$ with $Y_i = \\prod_{m\\ne i}(1+(\\lambda_m-\\lambda_i)^{-2})^{1/2}$. The identity $\\int_0^\\infty dt\\, e^{-x\\cosh t} = K_0(x)$ lets one integrate out the scattering shifts one by one, producing the exact modified log-gas joint law (2.5) for the Lax eigenvalues alone. Taking $N\\to\\infty$ in that law yields the free-energy functional $\\mathcal{F}_c[\\varrho]$ of Eq. (2.8), with a quadratic potential, an entropy term, and a two-body scattering-shift term, whose unique minimizer is the claimed DOS. The same scattering machinery produces the Toda limit at low density and the trigonometric Calogero TBA equations (C.4)–(C.5) at high density.","core_discovery":"On the paper's own terms, the discovery is that the thermal Lax DOS of the hyperbolic Calogero fluid, $\\varrho_N(\\lambda) = \\frac{1}{N}\\sum_i \\delta(\\lambda-\\lambda_i)$ for the Lax matrix $L_c$ with entries $[L_c]_{ij} = \\delta_{ij} p_j + i(1-\\delta_{ij})(2\\sinh((q_i-q_j)/2))^{-1}$, converges almost surely to a deterministic $\\varrho(\\lambda)$ as $N,\\ell\\to\\infty$ at fixed density $\\bar\\rho=N/\\ell$. For the $\\cosh$ trap $U_C = \\sum_i e^{-\\ell/2}\\cosh(q_i)$, the scattering-coordinate transformation yields the exact joint eigenvalue distribution (2.5), a modified log-gas whose Boltzmann factor is $\\exp(-\\frac{\\beta}{2}\\sum_i\\lambda_i^2)\\prod_i 2K_0(2 e^{-\\ell/2} Y_i)$ with $Y_i = \\prod_{m\\ne i}(1+(\\lambda_m-\\lambda_i)^{-2})^{1/2}$; at large $N$ the logarithm of this product becomes the free-energy functional (2.8). The authors verify numerically that this distribution reproduces the directly diagonalized Lax DOS, and that the same DOS appears for box and ring confinement, establishing boundary-condition independence. They also identify the Toda and rational or trigonometric Calogero limits.","pith_inferences":["Not in the paper, but if boundary independence holds as claimed, one can initialize generalized-hydrodynamic simulations of the hyperbolic Calogero fluid directly from the box DOS instead of the more expensive cosh-confined Bessel sampling.","Not in the paper, but the reported slow convergence at high density suggests a finite-size scaling study of the central DOS versus $N$; extracting the exponent would tell whether the approach to the flat and TBA profile is algebraic or logarithmic.","Not in the paper, but the Bessel-K interaction in (2.5) is a temperature-dependent deformation of the log-gas, so spectral statistics such as level spacing or spectral form factor may show a crossover from Wigner-like to integrable behavior as density changes; the paper lists such diagnostics as future work.","Not in the paper, but the variational functional (2.8) could in principle predict the full crossover among Gaussian, Toda, flat, and TBA regimes without particle-level simulation, by minimizing $\\mathcal{F}_c[\\varrho]$ numerically."],"forward_implications":["In the thermodynamic limit the thermal Lax DOS is deterministic and independent of confinement, so cosh, box, and ring simulations all target the same $\\varrho(\\lambda)$.","For the cosh trap, the exact joint distribution (2.5) is a valid substitute for diagonalizing $L_c$; the paper verifies that Monte Carlo sampling of the modified log-gas reproduces the directly computed DOS.","The low-density Calogero DOS is the Toda-chain DOS, and the high-density Calogero DOS is the rational or trigonometric Calogero DOS, so the Calogero fluid interpolates between two exactly solvable spectra without fitting parameters.","The limiting DOS is a building block for generalized hydrodynamics of integrable systems, since it is the state variable from which conserved charges and their currents are constructed.","At low density and high temperature the DOS is approximately Gaussian with variance $T$, while at high density it flattens; the paper quantifies deviations with the Binder cumulant."],"supporting_citations":[{"why":"Supplies the exact joint eigenvalue distribution and the large-N free-energy functional that form the paper's analytic core.","marker":"[44]"},{"why":"Provides the TBA equation for the trigonometric Calogero Lax DOS used as the high-density benchmark.","marker":"[46]"},{"why":"Derives the constant-pressure Toda Lax DOS and free-energy functional used for the low-density benchmark.","marker":"[16]"},{"why":"Introduces Lax-matrix random matrix ensembles and the spacing-distribution approach the paper builds on.","marker":"[43]"},{"why":"Provides the elliptic and Weierstrass construction of the periodic Lax pair used for the ring simulations.","marker":"[45]"},{"why":"Gives the exact expression for the Toda DOS via parabolic cylinder functions, used in Appendix B.","marker":"[62]"},{"why":"Supplies the invariant beta ensemble and Gauss-Wigner crossover result used to evaluate the exact Toda DOS.","marker":"[63]"}],"fun_headline_variants":["Calogero Lax spectrum: one density for all traps","Integrable structure fixes random matrix spectrum","Boundary-free Lax DOS from Calogero dynamics","Same eigenvalue density in cosh, box, ring traps","Modified log-gas emerges from Calogero Lax matrix"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The Monte Carlo results at system sizes up to $N=256$ are taken to represent the thermodynamic limit, even though the paper reports that the high-density cosh-confined results converge slowly and the particle density profile has not yet flattened.","fun_headline_variants_meta":{"raw":{"variants":["Calogero Lax spectrum: one density for all traps","Integrable structure fixes random matrix spectrum","Boundary-free Lax DOS from Calogero dynamics","Same eigenvalue density in cosh, box, ring traps","Modified log-gas emerges from Calogero Lax matrix"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000897,"raw_usage":{"total_tokens":3870,"prompt_tokens":953,"completion_tokens":2917,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":569,"completion_tokens_details":{"reasoning_tokens":2839}},"tokens_in":569,"tokens_out":2917,"duration_ms":19921,"temperature":1.0,"reasoning_tokens":2839,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T16:38:30.536639+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Simulate the high-density $\\cosh$-confined Calogero fluid (for example $\\bar\\rho = 11$, $T=1$) at $N = 512$ and $N = 1024$ and compare the Lax DOS and the scaled particle density profile with the box- and ring-confined results; if the curves continue to shift with $N$ or the density profile fails to flatten, the claimed boundary-condition independence has not been demonstrated.","supporting_citations":[{"cited_title":"Choquard, Classical and Quantum Partition Functions of the Calogero-Moser-Sutherland Model","cited_arxiv_id":null,"evidence_quote":"Provides the TBA equation for the trigonometric Calogero Lax DOS used as the high-density benchmark."},{"cited_title":"Bogomolny, O","cited_arxiv_id":null,"evidence_quote":"Introduces Lax-matrix random matrix ensembles and the spacing-distribution approach the paper builds on."},{"cited_title":"Calogero, Classical Many-Body Problems Amenable to Exact Treatments Springer (2014) https://doi.org/10.1007/3-540-44730-X","cited_arxiv_id":null,"evidence_quote":"Provides the elliptic and Weierstrass construction of the periodic Lax pair used for the ring simulations."},{"cited_title":"Opper, Analytical solution of the classical Bethe-ansatz equation for the Toda chain Phys","cited_arxiv_id":null,"evidence_quote":"Gives the exact expression for the Toda DOS via parabolic cylinder functions, used in Appendix B."},{"cited_title":"Allez, J","cited_arxiv_id":null,"evidence_quote":"Supplies the invariant beta ensemble and Gauss-Wigner crossover result used to evaluate the exact Toda DOS."}],"review_version":1}