{"id":"84c388c5-9ec7-4c3b-8692-5a773569289f","arxiv_id":"2603.14477","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A single directed-polymer transfer-matrix product reproduces the Tracy-Widom GUE, GOE, GSE and Baik-Rains one-point laws through different boundary contractions, while its leading eigenvalue shows t^{1/3} fluctuations with non-canonical statistics.","lead":"This paper shows that one random transfer-matrix product can generate the four main one-point fluctuation laws of (1+1)-dimensional directed polymers simply by changing how the matrix is contracted at its boundaries. It also finds that the logarithm of the product's largest eigenvalue fluctuates like t^{1/3} over an intermediate window but follows none of the standard Tracy-Widom or Baik-Rains distributions.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Baik–Rains identification rests on a tuned Brownian width σ_B without an independent stationary measure; the match may be a two-cumulant fitting artifact rather than evidence for the unified-contraction claim.","rationale":"The paper honestly discloses the fitting loop in Sec. III C, and its GUE/GOE/GSE reproductions are plausible and supported by standard KPZ expectations. The single load-bearing assumption is that the Brownian-weighted initial vector with σ_B=6.2 is close to the stationary measure; this is not tested independently. If it fails, the Baik–Rains identification is reduced to a two-cumulant fit, hurting the strongest claim that all subclasses are projections of one ensemble. The reader identified exactly this weakness. A two-point stationarity check would settle it. The half-space modification is a caveat to the 'same W' wording but is disclosed and not central. The ln λ1 claim is explicitly framed as intermediate and non-asymptotic, so it does not change the verdict. Thus no verdict adjustment is needed beyond the reader's CONDITIONAL.","tokens_in":10422,"tokens_out":8524,"duration_ms":92178,"concrete_test":"Independent stationary-measure check: with the same disorder ensemble and σ_B = 6.2, compute the spatial structure function S(r,t) = Var[F_Bw(x+r,t) − F_Bw(x,t)] for r ≪ t^{2/3} at t = 512. If S(r) ∼ r over a range of r, the Brownian profile is the stationary measure and the Baik–Rains match is physically grounded. If S(r) saturates or scales differently, run a σ_B scan (5.0, 5.6, 6.2, 6.8, 7.4) and report KS distances to the BR CDF; a sharp minimum at a single σ_B would indicate the cumulant match is a tuning artifact.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central numeric assertion that a Brownian-weighted contraction of W(t) realizes the Baik–Rains subclass is the least secure leg of the unification claim. Section III C admits that the exact stationary measure of the discrete dynamics is unknown and that σ_B is treated as an effective tuning parameter; σ_B = 6.2 is selected by requiring the late-time skewness and excess kurtosis to approach the Baik–Rains values. One free parameter can match two cumulants, so the agreement in Fig. 6 is not an independent test. The subsequent full-PDF comparison is visually suggestive but no quantitative distance is reported, and with 10^6 samples the expected statistical error on the PDF is tiny, making the absence of a KS/χ² statistic conspicuous. A closely related fragility is that the TW-GSE benchmark has excess kurtosis ≈0.04, so Fig. 5 is also weakly discriminating, though this affects a secondary identification. If the Brownian initial vector were truly the stationary measure, the same σ_B should also make two-point observables (e.g., Var[h(x+r)−h(x)]) linear in r; the paper explicitly declines to check this. Thus the load-bearing assumption is not independently verified.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper revisits the transfer-matrix formulation of (1+1)-dimensional directed polymers in random media and argues that a single time-ordered random matrix product W(t) organizes the canonical one-point KPZ fluctuation subclasses. For a fixed bulk disorder realization, the authors compute partition functions as matrix elements or contractions of W(t): point-to-point gives TW-GUE, point-to-line gives TW-GOE, a half-space absorbing-wall implementation gives TW-GSE, and a Brownian-weighted initial vector gives Baik-Rains. They report t^{1/3} free-energy fluctuation growth in all four cases and show standardized distributions and cumulants approaching the corresponding benchmarks using 10^6 disorder realizations with N=128 and sigma=3. The paper then studies the leading eigenvalue ln lambda_1(t), which displays an intermediate t^{1/3} fluctuation regime but whose standardized distribution remains distinct from the canonical TW/BR laws within the simulated range. The central message is that geometry-dependent KPZ subclasses can be viewed as projections of one transfer-matrix ensemble, and that spectral observables reveal fluctuation structures beyond endpoint geometries.","tokens_in":10596,"tokens_out":10299,"duration_ms":97595,"significance":"If the claims are correct, the paper offers a compact numerical demonstration that known one-point KPZ laws arise from different algebraic contractions of a single finite-dimensional random matrix product, and it identifies a spectral observable (ln lambda_1) with novel fluctuation behavior. The numerical effort is substantial, with 10^6 samples for every main figure, and the presentation is generally honest: the effective nature of the Brownian initial condition is stated explicitly, and the paper does not claim an asymptotic law for lambda_1(t). However, the Baik-Rains identification is the least secure leg of the unification claim because sigma_B is tuned to the same cumulants used for verification, with no independent stationarity check. The scaling claims also rely on guide lines rather than fitted exponents, and the GSE benchmark is close to Gaussian, making that identification weakly discriminating. These issues are fixable with additional quantitative analysis, so the paper is promising but not yet fully convincing.","major_comments":[{"comment":"The Baik-Rains identification is not an independent test. The paper states that the exact stationary measure is unknown and treats sigma_B as an effective tuning parameter, fixing it by requiring late-time skewness and excess kurtosis to approach the Baik-Rains values; these same two cumulants are then shown in the inset as evidence of convergence. With one free parameter, matching two cumulants is weak evidence. The full-PDF comparison is in principle independent, but no quantitative distance (KS/chi-square) is reported despite 10^6 samples. The paper also explicitly declines to test two-point stationarity (e.g., Var[h(x+r)-h(x)] linear in r), which would be an independent check of the same sigma_B. This leg of the unification claim needs either an independent determination of the effective stationary initial condition or a quantitative full-distribution comparison plus a sensitivity an","section":"Section III C, Eqs. (16)-(17), Fig. 6"},{"comment":"The t^{1/3} statements are supported by guide lines rather than fits. In Fig. 2, the four standard-deviation curves are compared with a single line of slope 1/3 without fitted exponents or confidence intervals; at late times some curves deviate visibly. In Fig. 7, the intermediate t^{1/3} window is not defined by any objective criterion, and the later crossover to t^{1/2} is indicated by a second guide line. Because the scaling is load-bearing for all four subclasses and for ln lambda_1(t), report fitted exponents with uncertainties over declared windows and justify the windows. Note also that N=128 and t=1024 are close to N^{3/2} approximately 1450, so finite-size contamination at late times should be assessed.","section":"Figs. 2 and 7"},{"comment":"The TW-GSE benchmark has skewness approximately 0.16 and excess kurtosis approximately 0.04, so it is close to Gaussian. Visual agreement with the GSE curve is therefore weakly discriminating; the data may also be consistent with a Gaussian or another near-Gaussian benchmark. Report quantitative distances (KS or Anderson-Darling) between the measured distribution and each candidate (Gaussian, TW-GOE, TW-GUE, TW-GSE) as a function of time. This is needed to support the GSE identification.","section":"Fig. 5, Eqs. (13)-(15)"},{"comment":"The phrase 'different boundary contractions of the same random matrix product W(t)' is too broad. The half-space realization in Eqs. (13)-(15) modifies the transfer matrix at the absorbing boundary (triangular block) and the boundary diagonal entries in Eq. (9); it is not a pure contraction of the unmodified W(t). The abstract qualifies this, but Sec. I and Sec. III do not. Either relax the claim to 'contractions and boundary-modified transfer rules' or show that the half-space observable can be written as a contraction of the same W(t) with a fixed boundary operator.","section":"Title/abstract and Secs. I/III"}],"minor_comments":[{"comment":"Skewness and excess kurtosis are plotted without error bars. With 10^6 samples, bootstrap errors would make the convergence claims quantitative.","section":"Figs. 3-6 insets"},{"comment":"State the numerical method and tolerance used for lambda_1(t), and specify the time window used for the t^{1/3} claim in the text or caption.","section":"Fig. 7"},{"comment":"The sign convention for the standardized free energy is nonstandard; state explicitly that positive skewness corresponds to the cited TW/BR values.","section":"Eq. (11)"},{"comment":"The guide line has no stated vertical offset, so it is hard to judge goodness of fit; consider plotting local effective exponents or residuals.","section":"Fig. 2"},{"comment":"A data/code availability statement would improve reproducibility; the paper currently provides no such statement.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The referee's main concern is the Baik-Rains section; the required additional checks (KS distances, an independent estimate of the stationary initial condition, fitted exponents) are well within the scope of the existing simulations. I do not see a fundamental flaw in the GUE/GOE results, and I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The thing to know: this is a solid numerical paper whose main genuinely new item is the statistics of ln λ1(t) for the DPRM transfer-matrix product. The four canonical subclass reproductions (GUE/GOE/GSE/Baik–Rains) are confirmatory and well executed, but the Baik–Rains leg has a fitted-input issue that the authors disclose but do not resolve.\n\nWhat's good: the transfer-matrix implementation is clean, the 10^6 sample statistics are respectable, and the GUE/GOE identifications are convincing. The spectral observable is new — I don't know of prior work on ln λ1(t) fluctuations for this product — and the paper is careful not to claim an asymptotic law. The explicit caveats in Sec. III C and IV are to their credit.\n\nSoft spots, in order: (1) Baik–Rains: σ_B = 6.2 is chosen so late-time skewness and excess kurtosis approach the BR benchmarks, and then those same cumulants are cited as evidence. That is a two-cumulant fit with one knob, not an independent test. The full-PDF comparison is visually suggestive but no KS or chi-square distance is reported; with 10^6 samples the statistical error bars are tiny, so a quantitative test is cheap and should be added. They also decline to check two-point observables, which means the 'stationary' claim rests entirely on one-point matching. (2) The t^{1/3} statements in Figs. 2 and 7 are supported by guide lines, not fitted exponents. Give me an exponent with uncertainty. (3) The GSE benchmark is nearly Gaussian (excess kurtosis ~0.04), so the half-space agreement is weakly discriminating. That one is secondary.\n\nNone of this breaks the paper. The core claim — that a single W(t) ensemble produces these subclasses through contractions — holds up for GUE, GOE, and arguably GSE; the Baik–Rains realization is plausible but needs a sharper test. The eigenvalue result stands as an interesting new observation.\n\nRecommended for peer review. A good referee would ask for fitted exponents, error bars on cumulants, and a BR identification that either tests σ_B sensitivity or includes a spatial roughness check. I'd cite it for the eigenvalue fluctuation and would take it to reading group.","headline":"Solid numerical paper: the four KPZ subclass reproductions are confirmatory, but the genuinely new item is ln λ1(t) of the transfer-matrix product showing t^{1/3} fluctuations with non-canonical statistics; the Baik–Rains leg rests on a tuned Brownian width and needs a sharper test.","tokens_in":11257,"tokens_out":2169,"would_cite":true,"duration_ms":20995,"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":"A single random transfer-matrix product generates the four canonical one-point KPZ fluctuation laws—Tracy–Widom GUE, GOE, GSE, and Baik–Rains—by different boundary contractions, while its largest eigenvalue defines a fluctuation law outside","keywords":["directed polymers","KPZ universality","Tracy-Widom distributions","Baik-Rains distribution","transfer matrix","random matrix products","one-point fluctuations","half-space boundary"],"falsifier":"Compute the exact stationary measure of the transfer-matrix Markov chain or simulate long enough to measure the stationary spatial covariance of the free-energy profile, then check whether the Brownian-weighted initial vector with fixed σ_B reproduces it; if the one-point Baik–Rains agreement requires retuning σ_B with time or system size, the identification is a fitting artifact rather than a true stationary realization.","tokens_in":10143,"feed_emoji":"🎲","tokens_out":4443,"duration_ms":43822,"temperature":0.7,"pith_summary":"The paper claims that in (1+1)-dimensional directed polymers, the whole family of one-point KPZ fluctuation laws can be generated from a single random matrix product: the time-ordered product of transfer matrices. Point-to-point geometry corresponds to one endpoint contraction and gives Tracy–Widom GUE; point-to-line gives GOE; a half-space absorbing boundary gives GSE; and a Brownian-weighted initial vector gives Baik–Rains. If correct, geometry is not a separate dynamical setting but an algebraic choice of how to contract one universal object. The paper also studies the logarithm of the largest eigenvalue of the product, which shows intermediate t^(1/3) growth but a standardized distribution that does not match any of the canonical benchmarks over the simulated times.","feed_headline":"All four KPZ fluctuation laws come from one transfer matrix","feed_subtitle":"Same product, different contractions: GUE, GOE, GSE, and Baik-Rains.","key_machinery":"The key machinery is the transfer-matrix product W(t) and the observation that boundary conditions are implemented as contractions. Point-to-point free energy is −ln⟨x₀|W(t)|x₀⟩, point-to-line is −ln Σₓ⟨x|W(t)|x₀⟩, the half-space absorbing boundary is encoded by a triangular block at the wall, and the stationary-type construction uses a Brownian-weighted initial vector e^{B(x)}. For the spectral observable, the Perron–Frobenius theorem guarantees that the leading eigenvalue λ₁(t) is real and positive, making F₁(t) = ln λ₁(t) a natural growth measure.","core_discovery":"The central object is W(t) = T(t)T(t−1)···T(1), the time-ordered product of random tridiagonal transfer matrices; every polymer partition function is a matrix element of this product. The paper's numerical evidence shows that different contractions—fixed endpoint, summed endpoint, absorbing boundary, and Brownian-weighted summation—produce standardized free-energy distributions consistent with Tracy–Widom GUE, GOE, GSE, and Baik–Rains, with t^(1/3) fluctuation growth and low-order cumulants approaching the corresponding benchmarks. A separate spectral observable, ln λ₁(t), grows with t^(1/3) over an intermediate window but has skewness and kurtosis that remain distinct from all four benchmar","pith_inferences":["If the unified-contraction picture survives longer simulations, any observable that is a function of W(t)—not just endpoint matrix elements—becomes a candidate for new universality subclasses; the natural next step is to compute correlation functions of the full spectrum.","The need to tune σ_B hints that the Brownian vector is a one-parameter family of approximate stationary measures; a sharper test would be to derive the true stationary measure of the transfer-matrix dynamics and compare, removing the fitting freedom.","The GSE benchmark is nearly Gaussian (excess kurtosis ~0.04), so the half-space agreement in Fig. 5 is weakly discriminating at the accessible times; confirming that identification will require either much longer runs or higher-order cumulants.","If ln λ₁(t) converges to a limiting law distinct from all canonical KPZ benchmarks, then the transfer-matrix product defines a larger 'master' ensemble whose projections include the KPZ subclasses plus a spectral sector with its own universality."],"forward_implications":["If the central claim is correct, the four KPZ subclasses are not separate stochastic evolutions but different linear-algebraic projections of a single random-matrix-product ensemble.","The Brownian-weighted contraction provides a practical finite-dimensional route to the stationary KPZ subclass even without knowing the exact invariant measure, with σ_B as an effective tuning parameter.","The spectral observable ln λ₁(t) exhibits KPZ-like t^(1/3) fluctuation growth while remaining outside the canonical Tracy–Widom and Baik–Rains subclasses in the simulated window, suggesting that matrix-level observables form a new class of fluctuation structures.","The framework invites a systematic study of the full spectrum of W(t), since spectral observables are determined by the internal structure of the product rather than by endpoint geometry."],"fun_headline_variants":["One transfer matrix yields all four KPZ fluctuation laws","Four KPZ classes from one transfer matrix product","Same product, four distinct fluctuation laws","One transfer matrix, four KPZ universality laws","Unified matrix product reproduces all KPZ subclasses"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The identification of the Brownian-weighted contraction with the Baik–Rains stationary subclass rests on the assumption that a random-walk initial vector with a single tuned width σ_B = 6.2 approximates the true stationary measure of the lattice dynamics, which the paper states is unknown; if that proxy is wrong, the agreement reduces to matching two cumulants with one free knob.","fun_headline_variants_meta":{"raw":{"variants":["One transfer matrix yields all four KPZ fluctuation laws","Four KPZ classes from one transfer matrix product","Same product, four distinct fluctuation laws","One transfer matrix, four KPZ universality laws","Unified matrix product reproduces all KPZ subclasses"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000604,"raw_usage":{"total_tokens":2667,"prompt_tokens":766,"completion_tokens":1901,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":510,"completion_tokens_details":{"reasoning_tokens":1838}},"tokens_in":510,"tokens_out":1901,"duration_ms":12496,"temperature":1.0,"reasoning_tokens":1838,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T18:10:15.240911+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the exact stationary measure of the transfer-matrix Markov chain or simulate long enough to measure the stationary spatial covariance of the free-energy profile, then check whether the Brownian-weighted initial vector with fixed σ_B reproduces it; if the one-point Baik–Rains agreement requires retuning σ_B with time or system size, the identification is a fitting artifact rather than a true stationary realization.","supporting_citations":[],"review_version":1}