{"id":"724983c4-b108-486a-a6ff-de913a028094","arxiv_id":"2608.03730","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Mutually avoiding directed polymers in a quenched random medium exhibit a density-squared interaction free energy, reproduced by filling single-polymer transfer-matrix eigenvalues.","lead":"A new numerical method fills the largest eigenvalues of a single-polymer transfer matrix to compute the free energy of many mutually avoiding directed polymers in a random medium. The results match the replica-theory prediction that disorder repels the polymers with a density-squared free-energy cost, and reveal a linear spectral edge.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Spectral filling rule (Eq. 3) assumes real positive eigenvalues of non-symmetric W(t); Appendix A's oscillatory proof is incomplete and Appendix B's phase check covers only the simulated windows, so the method's central assumption is not established.","rationale":"The paper proposes a spectral method to compute finite-density non-crossing DPRM free energies by filling logarithmic eigenvalues of a single-polymer transfer-matrix product. The central numerical observation is a linear upper spectral edge that, when filled, yields the replica Bethe ansatz result Δf(ρ)∼ρ². The argument is elegant and the full-filling determinant identity provides an exact anchor. However, the entire construction hinges on Eq. (3), which requires the filled eigenvalues of the non-symmetric W(t) to be real and positive. This is not generic for non-symmetric matrices, and the paper's proof in Appendix A is only an outline: total nonnegativity alone is insufficient, and the oscillatory property is asserted rather than demonstrated. Appendix B offers a numerical phase check, but it covers only the exact simulation parameters (N≤64, t≤512) and deliberately excludes the deep tail. Thus the method's foundational assumption is not established beyond the fitted regime. This is precisely the reader's weakest assumption, and it is load-bearing because a single counterexample would invalidate the spectral route. The appropriate verdict remains CONDITIONAL: the method is plausible and the reported data are consistent with the prediction, but the central spectral-positivity assumption requires either a rigorous proof or a broader numerical certification before the claim can be accepted. No change to the reader's verdict is needed.","tokens_in":10703,"tokens_out":16114,"duration_ms":204542,"concrete_test":"Run the paper's own phase diagnostic at N=64, t=1024 (t/N^{3/2}=2) and at N=128, t/N^{3/2}=1 (t≈1448) with wp=16000; if any nonreal eigenvalue appears in the filled sector, the spectral filling rule is not stable beyond the reported windows. In addition, for three disorder realizations at N=64,t=512, verify that all minors of W(t) up to size 4 are nonnegative (total nonnegativity) and that some small power W(t)^k has all minors strictly positive (oscillatory property), which would give direct evidence for the Gantmacher–Krein hypothesis in the tested regime.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that Δf(ρ)∼ρ² follows directly from the spectral filling rule ln Z_m ≃ Σ_{i=1}^m ln λ_i(t), which is valid only if the filled eigenvalues λ_i of the non-symmetric product W(t) are all real and positive. Perron–Frobenius guarantees only λ_1. Total nonnegativity of W(t) (Appendix A) does not by itself imply real positive eigenvalues for a non-symmetric matrix; one additionally needs the full oscillatory property (some power strictly totally positive). Appendix A's local 3×3 check (Eq. A3) and the assertion of 'oscillatory assumptions' do not constitute a proof for general N, disorder realizations, or times. The only numerical evidence is the phase diagnostic of Appendix B, which covers N≤64, t≤512—exactly the simulation windows—and explicitly excludes the deep spectral tail. It does not certify the long-time/thermodynamic limit in which the ρ² law is claimed. If any filled eigenvalue becomes complex or negative for a realization at larger N or t, Eq. (3) fails and the reported Δf∼ρ² could be an artifact of the high-precision rescaled construction rather than a property of the physical problem.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a numerical method to compute the quenched free energy of a finite density of mutually avoiding directed polymers in a random medium. The method avoids the exponentially large many-polymer transfer matrix by expressing the m-polymer partition function as a sum over products of the m largest eigenvalues of a single-polymer transfer-matrix product W(t), and then using the spectral filling rule (Eq. 3). The authors test the replica Bethe ansatz predictions: at low density, the interaction free energy Δf(ρ) ~ ρ^2, the second disorder cumulant scales as ρ^{1/2}, and a linear upper spectral edge is observed. They also derive an exact identity at full filling (Eq. 7–8). Numerical results for N = 48 and 64 at two scaled times t/N^{3/2} = 3/4 and 1 are presented and reported to be consistent with the predicted exponents.","tokens_in":11046,"tokens_out":3132,"duration_ms":39198,"significance":"If the method is valid, it provides a new numerical route to a longstanding problem in disordered line matter, and it directly connects the many-polymer free energy to the spectral properties of a random transfer-matrix product. The paper contains several strengths: the exact full-filling determinant identity (Eq. 8) is a nontrivial check; the high-precision stabilization procedure (Appendix B) is carefully described; and a data collapse is attempted at two scaled times. The linear spectral edge (Fig. 4) is a clear and testable signature that distinguishes this problem from pure free-fermion and standard random-matrix edges. However, the central claim rests on unproven real-positivity of the filled spectrum, the amplitude fits are not independently benchmarked, and the long-time limit is not established. These issues must be addressed before the result can be considered reliable.","major_comments":[{"comment":"The spectral filling rule ln Z_m(t) ≃ Σ_{i=1}^m ln λ_i(t) requires that the filled eigenvalues of the non-symmetric matrix W(t) are real and positive. Perron–Frobenius gives only λ_1. Appendix A invokes total nonnegativity and oscillatory matrices, but the argument is an outline: Eq. (A3) is a local 3×3 check, and no proof is given for general N, general disorder realizations, or times t→∞. Appendix B verifies realness numerically only up to N=64, t=512, exactly the simulation windows, and excludes the deepest spectral tail. Since the central claim Δf(ρ)~ρ^2 follows directly from this filling construction, a failure of real-positivity in the thermodynamic/long-time limit would invalidate the result. Please provide a rigorous proof of the oscillatory property, or a substantially more extensive numerical certification (e.g., larger N, longer t, and a check that all filled levels are real f","section":"Eq. (3) and Appendix A/B"},{"comment":"The linear fit F_m(t)/(mt) = v_1 − a_1 ρ (Eq. 9) uses a_1 fitted from the same data that are then used to claim Δf(ρ) ~ ρ^2. This tests the exponent but not the amplitude: any data lying on a smooth curve could be locally fit by a line. Likewise, the second cumulant fit a_2 ρ^{1/2} in Fig. 3 is a two-parameter fit. No error bars are shown for any of the plotted quantities. The collapse between N=48 and N=64 at two scaled times is suggestive but does not quantify finite-size or time deviations. Please provide error estimates, show the fit residuals, and if possible compare a_1 and a_2 to independent predictions from the dilute Bose gas mapping.","section":"Eq. (9), Figs. 2–4"},{"comment":"The paper claims an 'approximately stationary long-time regime' based on data at t/N^{3/2} = 3/4 and 1. This is insufficient to establish the t→∞ limit in which the replica predictions apply. The spectral edge in Fig. 4 is also shown only at these two times. Without a systematic extrapolation in t (or a scaling collapse demonstrating convergence), the measured exponents could still be pre-asymptotic. Please show that the results are stable when t is increased at fixed N (or at fixed t/N^{3/2} for several values), and provide a quantitative criterion for convergence.","section":"Figs. 2–4 and the long-time limit"},{"comment":"The text states that 'data (not shown) indicate a third cumulant' and claims consistency with a_3 ≈ 0.15 a_2^2/a_1, but no third-cumulant data are presented. Since the cumulant hierarchy is a central prediction of the replica Bethe ansatz, the omission of third (and higher) cumulants leaves the claim untested. Please include the third-cumulant results or explicitly retract the consistency statement until such data are available.","section":"Third cumulant paragraph"}],"minor_comments":[{"comment":"The paper switches between one-based (Eq. 11) and zero-based (Fig. 4) indexing of eigenvalues. This is confusing; please use a single convention or clearly mark the switch in every equation/figure.","section":"General notation"},{"comment":"The inset shows data only for ρ ≤ 1/4. Please state the fit range explicitly and justify why this is the 'dilute window'.","section":"Fig. 2 inset"},{"comment":"The reconstruction ln|λ_i(s)| = ln|λ̃_i(s)| + C_s assumes that the rescaling does not mix eigenvalues of non-normal W(s). Please comment on the validity of this step for non-normal matrices and on the numerical accuracy of the computed eigenvalues for the filled sector.","section":"Appendix B, Eq. (B1)"},{"comment":"Reference [25] is cited as 'in press (2026)'. If it is not yet published, please provide an arXiv identifier or clarify its availability.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The paper presents an appealing spectral method and an exact full-filling identity, but the central claim is currently supported only by numerical evidence in a narrow parameter window and by a sketch of an oscillatory-matrix proof. The authors should be asked to either supply a rigorous proof of the real-positivity of the filled spectrum or to significantly expand the numerical certification range. The fit-based amplitude tests are not yet convincing without error bars and an independent amplitude prediction. If these issues are addressed, the paper could be a nice contribution. I also note that the third-cumulant claim should not be made without showing the data."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know this before reading: the central result—Delta f ~ rho^2 from filling a linear spectral edge—is a genuine numerical claim with an elegant method, but the paper's own appendices don't fully close the gap they identify. The method is to compute the single-polymer transfer-matrix product W(t) and fill its largest m log-eigenvalues to get the m-polymer free energy. That's a nice trick, building on Ref. 24, and it avoids the combinatorial many-polymer transfer matrix. The paper does a couple of things well. The exact full-filling identity det W(t) = exp(-sum 2E) gives a closed-form anchor at rho=1, and the disorder-averaged cumulative growth per filled level collapses nicely at two scaled times for N=48,64. The linear upper edge in Fig. 4 is a clean observation, distinct from both the pure free-fermion quadratic edge and the GOE square-root edge. If the filling rule holds, the rho^2 law follows by integration, and the second-cumulant sqrt(rho) scaling is at least visually consistent. The soft spot is exactly what the stress-test flags. Eq. 3 assumes all filled eigenvalues of the non-symmetric W(t) are real and positive. Perron-Frobenius gives only the leading one. Appendix A argues total nonnegativity plus oscillatory structure, but it stops at an outline—the Gantmacher-Krein conditions are asserted, not verified, for this family of disorder realizations. Appendix B checks realness only up to N=64, t=512, precisely the simulated range, and it explicitly excludes the deep spectral tail. So the long-time limit in which the rho^2 law is phrased is not actually verified: the filled sector could conceivably develop complex levels at larger t/N^{3/2}, and if it did, the apparent linear edge might be an artifact of the rescaling. The paper is honest about the limitation, but honesty isn't proof. Two more things. The amplitudes a1 and a2 are fitted from the same data used to test the theory, so the test validates exponents, not amplitudes. And the third cumulant—which would discriminate the Bose-gas mapping—is explicitly not shown; the paper says the data are consistent but presents no numbers. There are no error bars anywhere, and no code or data. None of these are fatal; they're addressable. Who is this for? People working on directed polymers, vortex lines, or non-intersecting ensembles will want to know the method and the rho^2 candidate. It deserves a serious referee—conditional acceptance at best, but a referee can push for the missing proof or at least a longer-time phase check, error bars, and the third-cumulant data. My recommendation: engage with it, but don't cite the central law as established yet.","headline":"A clever spectral-filling method for finite-density directed polymers, but the core real-positivity assumption is only partially supported, so the rho^2 law remains a well-argued numerical conjecture rather than a settled benchmark.","tokens_in":11496,"tokens_out":1973,"would_cite":false,"duration_ms":22709,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Filling a single-polymer spectrum yields the ρ² law of directed-polymer glasses.","keywords":["directed polymers in random media","quenched disorder","transfer-matrix spectrum","replica Bethe ansatz","non-crossing paths","spectral edge","free-energy density","disorder cumulants"],"falsifier":"Compute the filled-sector phase diagnostic at working precision beyond 16000 or for t > 512 and check whether any nonreal eigenvalue pair enters the top m levels; a nonreal pair would invalidate the sum-of-log-eigenvalues rule and call the ρ² law into question. A complementary check is to measure F_m(t)/(mt) at densities below 0.03 on larger N: if the linear relation bends, the quadratic free-energy law is a finite-window artifact.","tokens_in":10611,"feed_emoji":"🧵","tokens_out":12525,"duration_ms":131212,"temperature":0.7,"pith_summary":"This paper shows that many mutually avoiding directed polymers in a frozen random environment can be studied without constructing the exponentially large many-polymer transfer matrix. The key equivalence: non-crossing polymers behave like fermions, so at long times the traced many-polymer partition function is the product of the largest eigenvalues of a single-polymer transfer-matrix product. Using this spectral filling rule, the authors compute the free-energy density and its sample-to-sample fluctuations. The central numerical observation is that the mean growth per filled logarithmic level falls linearly with density, which corresponds to a disorder-induced linear upper edge of the logarithmic spectrum; filling this edge gives an interaction free energy Δf(ρ) ~ ρ², matching the replica Bethe ansatz prediction. An exact determinant identity at full filling anchors the unit-density limit, and the second disorder cumulant follows the predicted ρ^{1/2} scaling.","feed_headline":"Filling one polymer spectrum produces the ρ² glass law","feed_subtitle":"Non-crossing polymers act like fermions; filling a spectrum yields the ρ² repulsion law.","key_machinery":"The central object is the single-polymer transfer-matrix product W(t)=T(t)T(t-1)...T(1), whose ordered logarithmic eigenvalues ε_i(t)=t^{-1} ln λ_i(t) define the filled levels. The spectral filling rule ln Z_m(t) ≈ Σ_{i=1}^m ln λ_i(t) reduces the many-polymer trace to a sum over these levels. The mathematical load is carried by total nonnegativity of W(t): by the standard non-intersecting-path determinant identity, every minor is a sum of positive weights of non-crossing path families, giving nonnegative minors. Because W(t) is non-symmetric, the paper adds the oscillatory-matrix structure of the nearest-neighbor transfer product (totally nonnegative matrices whose powers become strictly pos","core_discovery":"The paper's central claim is numerical: at low density ρ = m/N, the disorder-averaged mean cumulative logarithmic growth per filled level is linear in ρ, F_m(t)/(mt) = v_1 - a_1 ρ. Spectrally this means the disorder-averaged logarithmic band has a linear upper edge, ε_0(t) - ε_i(t) ∝ r near the top; filling that edge gives Δf(ρ) ~ ρ², the replica Bethe ansatz prediction. This contrasts with the ρ³ interaction law of the pure non-crossing problem and with the different exponent obtained by filling a standard random-matrix soft edge. The same filled spectrum reproduces the predicted cumulant scalings, and the fully packed determinant identity anchors the ρ = 1 limit exactly.","pith_inferences":["The linear edge may be a generic disorder signature for non-crossing line ensembles; a natural test is to repeat the filling construction on other transfer-matrix ensembles with the same nonnegative-minor path structure and check whether the same ρ² law appears.","Because the method avoids the many-body transfer matrix, it could probe the crossover from the dilute ρ² regime to full filling, where replica scaling is not expected and the exact determinant identity takes over.","The connection to oscillatory-matrix theory suggests an analytic route: the edge shape might be derivable from total nonnegativity plus disorder averaging, turning the ρ² law into a theorem rather than a numerical observation.","The spectral object introduced here likely has its own random-matrix statistics; measuring the distribution of the largest filled log-eigenvalue would test whether its edge fluctuations match standard random-matrix classes or a new universality class."],"forward_implications":["The many-polymer glass free energy, its interaction part, and its disorder cumulants become accessible from a single-polymer spectral computation over a range of densities and times.","The replica Bethe ansatz ρ² law acquires a concrete spectral mechanism: a linear disorder-induced upper edge of the logarithmic band, distinct from the pure and standard-random-matrix cases.","The fully packed determinant identity fixes the spectrum at ρ = 1, so the same construction connects the dilute scaling regime to an exact unit-density limit.","The filled spectrum gives a route to higher disorder cumulants; in particular the third cumulant amplitude can be predicted from the first two and tested with additional statistics.","If the linear edge is robust under changes of disorder and transfer-matrix ensembles, the spectral filling construction carries over to vortex-line arrays, non-crossing disordered interfaces, and other line ensembles."],"supporting_citations":[{"why":"Supplies the dilute disordered-system prediction Δf(ρ)∼ρ² that the numerics target.","marker":"[17]"},{"why":"Provides the replica Bethe ansatz treatment of disordered interfaces whose low-density law is tested.","marker":"[18]"},{"why":"Predicts the disorder cumulant hierarchy F^p_m/(Nt) ∼ a_p ρ^{(5-p)/2} compared with the data.","marker":"[19]"},{"why":"Supplies the dilute Bose-gas mapping giving exact cumulant amplitudes, including the predicted a_3 relation.","marker":"[20]"},{"why":"Supplies the determinant formula identifying non-intersecting path partition sums with minors.","marker":"[21]"},{"why":"Provides the vector-representation path-determinant version used in the positivity argument.","marker":"[22]"},{"why":"Provides the binomial-determinant/path formula completing the identity for nonnegative minors.","marker":"[23]"},{"why":"Early use of spectral filling for vortex lines pinned by columnar defects; precursor of the present method.","marker":"[24]"},{"why":"Defines the specific transfer-matrix ensemble whose product spectrum is filled.","marker":"[25]"},{"why":"Gives the oscillatory-matrix theory converting total nonnegativity into real positive eigenvalues.","marker":"[26]"}],"fun_headline_variants":["Filling polymer spectrum yields ρ² glass law","Edge filling gives ρ² repulsion for directed polymers","Spectrum method reproduces replica ρ² law","Filling spectral edge confirms ρ² interaction","Single-polymer spectrum yields ρ² glass law"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The filling rule requires the top m eigenvalues of the non-symmetric product W(t) to be real and positive for the times, densities, and disorder realizations used; the paper argues this via oscillatory-matrix theory but verifies it numerically only up to working precision 16000 and t = 512, with the deep spectral tail excluded.","fun_headline_variants_meta":{"raw":{"variants":["Filling polymer spectrum yields ρ² glass law","Edge filling gives ρ² repulsion for directed polymers","Spectrum method reproduces replica ρ² law","Filling spectral edge confirms ρ² interaction","Single-polymer spectrum yields ρ² glass law"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00046,"raw_usage":{"total_tokens":2086,"prompt_tokens":638,"completion_tokens":1448,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":382,"completion_tokens_details":{"reasoning_tokens":1391}},"tokens_in":382,"tokens_out":1448,"duration_ms":13876,"temperature":1.0,"reasoning_tokens":1391,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T13:54:42.730564+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the filled-sector phase diagnostic at working precision beyond 16000 or for t > 512 and check whether any nonreal eigenvalue pair enters the top m levels; a nonreal pair would invalidate the sum-of-log-eigenvalues rule and call the ρ² law into question. A complementary check is to measure F_m(t)/(mt) at densities below 0.03 on larger N: if the linear relation bends, the quadratic free-energy law is a finite-window artifact.","supporting_citations":[{"cited_title":"Kardar and D","cited_arxiv_id":null,"evidence_quote":"Supplies the dilute disordered-system prediction Δf(ρ)∼ρ² that the numerics target."},{"cited_title":"Kardar, Replica Bethe ansatz studies of two- dimensional interfaces with quenched random impurities, Nucl","cited_arxiv_id":null,"evidence_quote":"Provides the replica Bethe ansatz treatment of disordered interfaces whose low-density law is tested."},{"cited_title":"Emig and M","cited_arxiv_id":null,"evidence_quote":"Predicts the disorder cumulant hierarchy F^p_m/(Nt) ∼ a_p ρ^{(5-p)/2} compared with the data."},{"cited_title":"Emig and M","cited_arxiv_id":null,"evidence_quote":"Supplies the dilute Bose-gas mapping giving exact cumulant amplitudes, including the predicted a_3 relation."},{"cited_title":"Karlin and J","cited_arxiv_id":null,"evidence_quote":"Supplies the determinant formula identifying non-intersecting path partition sums with minors."},{"cited_title":"Lindstr¨ om, On the vector representations of induced matroids, Bull","cited_arxiv_id":null,"evidence_quote":"Provides the vector-representation path-determinant version used in the positivity argument."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the binomial-determinant/path formula completing the identity for nonnegative minors."},{"cited_title":"Polkovnikov, Y","cited_arxiv_id":null,"evidence_quote":"Early use of spectral filling for vortex lines pinned by columnar defects; precursor of the present method."},{"cited_title":"Directed Polymer Transfer Matrices as a Unified Generator of Distinct One-Point Fluctuation Laws","cited_arxiv_id":"2603.14477","evidence_quote":"Defines the specific transfer-matrix ensemble whose product spectrum is filled."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the oscillatory-matrix theory converting total nonnegativity into real positive eigenvalues."}],"review_version":1}