Pith. sign in

REVIEW 4 major objections 4 minor 37 references

A Spectral Route to Directed-Polymer Glasses

T0 review · 4 major / 4 minor · reviewed 2026-08-05 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Filling a single-polymer spectrum yields the ρ² law of directed-polymer glasses.

desk verdict 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. read the letter →

arxiv 2608.03730 v1 pith:VONE3XXZ submitted 2026-08-04 cond-mat.dis-nn cond-mat.softcond-mat.stat-mechmath-phmath.MP

Sen Mu , Abbas Ali Saberi , Mehran Kardar This is my paper · ORCID
classification cond-mat.dis-nncond-mat.softcond-mat.stat-mechmath-phmath.MP
keywords directedpolymersinrandommediaquencheddisordertransfer-matrixspectrumreplicaBetheansatznon-crossingpathsspectraledgefree-energydensitycumulants
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

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.

What carries the argument

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

What would settle it

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.

Watch

Extended reading notes

Core claim

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.

Load-bearing premise

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.

Editorial extensions

If this is right

  • 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.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • 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.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 4 minor

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.

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 (4)
  1. [Eq. (3) and Appendix A/B] 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
  2. [Eq. (9), Figs. 2–4] 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.
  3. [Figs. 2–4 and the long-time limit] 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.
  4. [Third cumulant paragraph] 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.
minor comments (4)
  1. [General notation] 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.
  2. [Fig. 2 inset] The inset shows data only for ρ ≤ 1/4. Please state the fit range explicitly and justify why this is the 'dilute window'.
  3. [Appendix B, Eq. (B1)] 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.
  4. [References] Reference [25] is cited as 'in press (2026)'. If it is not yet published, please provide an arXiv identifier or clarify its availability.

Circularity Check

0 steps flagged · score 1.0 of 10

No significant circularity: spectral filling is an algebraic construction, the rho^2 comparison uses external replica theory, and self-citations are not load-bearing.

full rationale

The central chain is not circular. Eq. (3) is an algebraic identity given the (separately checked) assumption that the filled eigenvalues of W(t) are real and positive; Eq. (4) then defines the free energy in terms of those eigenvalues. The rho^2 result is obtained by fitting the numerically computed Fm(t)/(mt) to a line in rho (Eq. 9) and substituting into Eq. (4). This is an empirical extraction of an exponent, not a prediction of a quantity already used as input: the exponent rho^2 is the replica Bethe ansatz prediction from Refs. [17-20], which are prior independent theoretical works, even though one author overlaps. The amplitudes a1 and a2 are fitted, and the paper does not claim to predict their values; the cumulant prediction for a3 uses the fitted a1 and a2 to test a relation from [20], again an external consistency check. The full-filling determinant identity (Eqs. 5-8) is exact and anchors the high-density limit without using the low-density fit. Self-citations to [25] for the transfer-matrix ensemble and to [19,20] for the theory are background or external benchmarks, not circular justifications. The unproved real-positive spectrum of non-symmetric W(t) is a correctness/assumption risk, not a circularity, because no equation reduces to itself and the numerical phase check is an independent diagnostic. Overall, the paper is self-contained against external theoretical predictions, so circularity is minimal.

Assumptions & free parameters 4 free parameters · 5 assumptions · 0 invented entities

The central derivation rests on standard determinant combinatorics, on an unproven spectral-positivity assumption for a non-normal product matrix, and on a finite-time stationarity assertion. The quality of the theory test is limited because the amplitudes a_1, a_2 are fitted rather than predicted.

free parameters (4)
  • a_1 = ≈0.9
    Slope of the linear fit F_m(t)/(mt)=v_1 - a_1 rho in Fig. 2 inset; model-dependent amplitude of the rho^2 interaction term.
  • a_2 = ≈2.9
    Amplitude of the second-cumulant scaling F_2^c/(mt)=a_2 rho^{1/2} fitted in Fig. 3 inset.
  • v_1 = ≈4.2
    Intercept of the same fit; the single-polymer growth rate, used in Eq. (10) to subtract the single-polymer free energy.
  • scaled time window = 3/4 and 1 (t/N^{3/2})
    The two scaled times used to define the 'approximate stationary regime'; chosen by hand with no systematic t→infinity extrapolation.
assumptions (5)
  • standard math Karlin-McGregor / Lindström-Gessel-Viennot determinant identity: minors of W(t) count non-crossing path families.
    Used in Eq. (2) and Appendix A; standard combinatorial result, assumed valid for the planar directed graph.
  • domain assumption W(t) is oscillatory (nonsingular totally nonnegative with strict positivity in ordered sectors), so all eigenvalues are real and positive.
    The applicability of Gantmacher-Krein theory to this specific non-symmetric product matrix is argued in Appendix A but not proven; it is verified numerically in Fig. B1 at wp=16000 over the reported window.
  • domain assumption At long times, ln Z_m(t) is dominated by the product of the m largest eigenvalues (Eq. 3).
    Requires real positive levels with a gap between level m and m+1; the paper uses finite times t/N^{3/2}=3/4,1 and asserts stationarity from data collapse rather than a t→infinity extrapolation.
  • domain assumption Quenched free energy is obtained by the disorder average of ln Z_m (self-averaging).
    Standard in DPRM studies; the paper does not check finite-size self-averaging explicitly.
  • ad hoc to paper The alternating-sign bond choice (±E on the two sides of each cell) makes the transfer matrix a positive path matrix.
    The paper states 'This choice is somewhat arbitrary but ensures that the transfer matrix reflects an underlying planar graph.' The entire total-nonnegativity argument depends on this design.

how reviews work

0 comments
Cite this review

Pith. "Pith review of A Spectral Route to Directed-Polymer Glasses." pith.science (2026). https://pith.science/paper/VONE3XXZ

@misc{pith2026260803730,
  author       = {Pith},
  title        = {Pith review of: A Spectral Route to Directed-Polymer Glasses},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/VONE3XXZ}},
  note         = {Machine review of arXiv:2608.03730}
}
abstract

A finite density of mutually avoiding directed polymers in a quenched random medium is a minimal model of glassy line matter. The dilute theory, solved by replica Bethe ansatz, predicts an interaction free energy proportional to $\rho^2$ and disorder cumulants with distinct power-law dependences on the density $\rho$, but direct numerical tests have been hindered by the combinatorially large many-polymer transfer matrix. We recast the problem as filling logarithmic eigenvalues of a single-polymer transfer-matrix product, obtaining the quenched free energy, its cumulants, and a disorder-induced linear spectral edge consistent with the replica prediction.

Figures

Figures reproduced from arXiv: 2608.03730 by the authors.

Figure 1
Figure 1. FIG. 1. Finite-density mutually avoiding directed polymers [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. Mean cumulative growth per filled level, [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3. Second disorder cumulant of the filled free energy, [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

37 extracted references · 37 canonical work pages

  1. [1]

    Kardar and Y.-C

    M. Kardar and Y.-C. Zhang, Scaling of directed polymers in random media, Phys. Rev. Lett.58, 2087 (1987)

    work page

  2. [2]

    Halpin-Healy and Y.-C

    T. Halpin-Healy and Y.-C. Zhang, Kinetic roughening phenomena, stochastic growth, directed polymers and all that, Phys. Rep.254, 215–414 (1995)

    work page 1995

  3. [3]

    Zygouras, Directed polymers in a random environ- ment: A review of the phase transitions, Stochastic Pro- cess

    N. Zygouras, Directed polymers in a random environ- ment: A review of the phase transitions, Stochastic Pro- cess. Appl.177, 104431 (2024)

    work page 2024

  4. [4]

    D. A. Huse and C. L. Henley, Pinning and roughening of domain walls in Ising systems due to random impurities, Phys. Rev. Lett.54, 2708–2711 (1985)

    work page 1985

  5. [5]

    D. S. Fisher, Interface fluctuations in disordered systems: 5−ϵexpansion and failure of dimensional reduction, Phys. Rev. Lett.56, 1964–1967 (1986)

    work page 1964

  6. [6]

    Kardar, G

    M. Kardar, G. Parisi, and Y.-C. Zhang, Dynamic scal- ing of growing interfaces, Phys. Rev. Lett.56, 889–892 (1986)

    work page 1986

  7. [7]

    M. P. A. Fisher, Vortex-glass superconductivity: A pos- sible new phase in bulk high-T c oxides, Phys. Rev. Lett. 62, 1415–1418 (1989)

    work page 1989

  8. [8]

    Blatter, M

    G. Blatter, M. V. Feigel’man, V. B. Geshkenbein, A. I. Larkin, and V. M. Vinokur, Vortices in high-temperature superconductors, Rev. Mod. Phys.66, 1125–1388 (1994)

    work page 1994

Show all 37 references
  1. [9]

    V. L. Pokrovsky and A. L. Talapov, Ground state, spec- trum, and phase diagram of two-dimensional incommen- surate crystals, Phys. Rev. Lett.42, 65–67 (1979)

    1979

  2. [10]

    S. N. Coppersmith, D. S. Fisher, B. I. Halperin, P. A. Lee, and W. F. Brinkman, Dislocations and the commensurate-incommensurate transition in two dimen- sions, Phys. Rev. B25, 349–363 (1982)

    1982

  3. [11]

    M. E. Fisher and D. S. Fisher, Wall wandering and the dimensionality dependence of the commensurate- incommensurate transition, Phys. Rev. B25, 3192–3198 (1982)

    1982

  4. [12]

    Borodin, I

    A. Borodin, I. Corwin, and P. L. Ferrari, Free energy fluctuations for directed polymers in random media in 1 + 1 dimension, Commun. Pure Appl. Math.67, 1129– 1214 (2014)

    2014

  5. [13]

    O’Connell and J

    N. O’Connell and J. Warren, A multi-layer extension of the stochastic heat equation, Commun. Math. Phys.341, 1–33 (2016)

    2016

  6. [14]

    De Luca and P

    A. De Luca and P. Le Doussal, Crossing probability for directed polymers in random media, Phys. Rev. E92, 040102(R) (2015)

    2015

  7. [15]

    De Luca and P

    A. De Luca and P. Le Doussal, Crossing probability for directed polymers in random media. II. Exact tail of the distribution, Phys. Rev. E93, 032118 (2016)

    2016

  8. [16]

    Barraquand and P

    G. Barraquand and P. Le Doussal, A stationary model of non-intersecting directed polymers, J. Phys. A: Math. Theor.56, 045001 (2023)

    2023

  9. [17]

    Kardar and D

    M. Kardar and D. R. Nelson, Commensurate- incommensurate transitions with quenched random impurities, Phys. Rev. Lett.55, 1157–1160 (1985)

    1985

  10. [18]

    Kardar, Replica Bethe ansatz studies of two- dimensional interfaces with quenched random impurities, Nucl

    M. Kardar, Replica Bethe ansatz studies of two- dimensional interfaces with quenched random impurities, Nucl. Phys. B290, 582–602 (1987)

    1987

  11. [19]

    Emig and M

    T. Emig and M. Kardar, Thermodynamic fingerprints of disorder in flux line lattices and other glassy mesoscopic systems, Phys. Rev. Lett.85, 2176–2179 (2000)

    2000

  12. [20]

    Emig and M

    T. Emig and M. Kardar, Probability distributions of line lattices in random media from the 1D Bose gas, Nucl. Phys. B604, 479–510 (2001)

    2001

  13. [21]

    Karlin and J

    S. Karlin and J. McGregor, Coincidence probabilities, Pacific J. Math.9, 1141–1164 (1959)

    1959

  14. [22]

    Lindstr¨ om, On the vector representations of induced matroids, Bull

    B. Lindstr¨ om, On the vector representations of induced matroids, Bull. London Math. Soc.5, 85–90 (1973)

    1973

  15. [23]

    I. M. Gessel and G. Viennot, Binomial determinants, paths, and hook length formulae, Adv. Math.58, 300– 321 (1985)

    1985

  16. [24]

    Polkovnikov, Y

    A. Polkovnikov, Y. Kafri, and D. R. Nelson, Vortex pin- ning by a columnar defect in planar superconductors with point disorder, Phys. Rev. B71, 014511 (2005)

    2005

  17. [25]

    S. Mu, A. A. Saberi, R. Moessner, and M. Kardar, Di- rected polymer transfer matrices as a unified generator of distinct one-point fluctuation laws,Phys. Rev. E, in press (2026); arXiv:2603.14477

    2026 arXiv

  18. [26]

    F. R. Gantmacher and M. G. Krein,Oscillation Matrices 6 and Kernels and Small Vibrations of Mechanical Systems (Gostekhizdat, Moscow, 1950; English translation, AMS Chelsea, Providence, 2002)

    1950

  19. [27]

    Karlin,Total Positivity, Vol

    S. Karlin,Total Positivity, Vol. I (Stanford University Press, Stanford, 1968)

    1968

  20. [28]

    Ando, Totally positive matrices, Linear Algebra Appl

    T. Ando, Totally positive matrices, Linear Algebra Appl. 90, 165–219 (1987)

    1987

  21. [29]

    Brunet and B

    E. Brunet and B. Derrida, Probability distribution of the free energy of a directed polymer in a random medium, Phys. Rev. E61, 6789–6801 (2000)

    2000

  22. [30]

    Brunet and B

    E. Brunet and B. Derrida, Ground state energy of a non- integer number of particles withδattractive interactions, Physica A279, 398–407 (2000)

    2000

  23. [31]

    Barraquand and P

    G. Barraquand and P. Le Doussal, Large time cumulants of the KPZ equation on an interval, J. Stat. Phys.192, 111 (2025)

    2025

  24. [32]

    E. P. Wigner, On the distribution of the roots of certain symmetric matrices, Ann. Math.67, 325–327 (1958)

    1958

  25. [33]

    M. L. Mehta,Random Matrices, 3rd ed. (Else- vier/Academic Press, Amsterdam, 2004), ISBN 978-0- 12-088409-4

    2004

  26. [34]

    C. A. Tracy and H. Widom, Level-spacing distributions and the Airy kernel, Commun. Math. Phys.159, 151–174 (1994)

    1994

  27. [35]

    C. A. Tracy and H. Widom, On orthogonal and sym- plectic matrix ensembles, Commun. Math. Phys.177, 727–754 (1996)

    1996

  28. [36]

    https://www.pks.mpg.de/asg2024

  29. [37]

    Chu, Fluctuating interfaces and paths in disordered and non-equilibrium systems, https://hdl.handle.net/1721.1/123352

    S. Chu, Fluctuating interfaces and paths in disordered and non-equilibrium systems, https://hdl.handle.net/1721.1/123352. End Matter The End Matter collects two technical ingredients used in the main text. Appendix A explains the pos- itivity structure underlying the filled lo...

    2000

Pith tools

Reviewed August 5, 2026 · model on record in the stance chip above.