{"id":"81082cdf-9787-4b4f-b2d4-d7105b098158","arxiv_id":"2505.00073","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Sequential random MPS ensembles are biased relative to the natural Fubini-Study measure, and the paper derives the correcting weight plus a sampling algorithm.","lead":"This paper shows that the usual way of generating random matrix product states, by placing independent random unitaries at each site, samples the state space unevenly and favors one end of the chain. It derives the missing weight that makes the sampling uniform and gives a Metropolis algorithm for drawing from the corrected ensemble.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The measure-correction rests on an unproven local metric determinant; a direct numerical check of det g versus ∏|Γ_i|^{D(d-1)} is required.","rationale":"The reader's strongest claim is accurate, and the reference-point extension idea is a valid way to obtain a global Radon-Nikodym density from local data. The weakest assumption is not so much the surjectivity of the parameterization, since a local exponential chart exists at every point and the cut locus is measure zero; rather it is the unreported local metric determinant formula. I verified the dimensional logic (per-site tangent dimension (d-1)D², metric Tr[δx†δ'x Γ_i] giving det Γ_i^{(d-1)D}) and found it plausible, but the paper does not provide the calculation, and the SM contains a separate Jacobian inconsistency. The proposed test directly computes the Gram determinant and settles the identity. If the test passes, the central claim is supported and CONDITIONAL would become ACCEPT after minor clarifications; if it fails, the central claim is wrong. Thus the reader's CONDITIONAL verdict is appropriate and unchanged.","tokens_in":10753,"tokens_out":51107,"duration_ms":536409,"concrete_test":"For N=3, d=2, D=2, sample 10 random MPS from the RMPS ensemble. For each, compute the right environments Γ_i from the recurrence Γ_{i-1} = ∑_n A[i]_n Γ_i A[i]_n† with Γ_N = I_D. Then construct the 12×12 complex Gram matrix G of the tangent vectors ∂_{α,a,i}|Ψ⟩ obtained by finite differences of the exponential parameterization centered at that reference point, and compute det G. Verify det G = ∏_{i=1}^3 det(Γ_i)^{(d-1)D} = det(Γ_1)^2 det(Γ_2)^2 det(Γ_3)^2 to relative precision 1e-6. Repeat for one point with a degenerate Γ_i eigenvalue near zero; if the identity fails at any point, the measure correction is disproven. A second check: recompute the Jacobian J in Eq. (8) using the SM's derivation to resolve the Δ²(x̂) versus Δ²(x̂²) discrepancy.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central result dμ_FS = dμ_RMPS ∏|Γ_i|^{D(d-1)} (Eq. 12) is derived from the local Fubini-Study volume element at a reference point. That volume element, Eq. (11), is obtained from the asserted determinant formula |g| = ∏|Γ_i|^{D(d-1)} (Eq. 10), which is stated as 'standard MPS contraction calculations' without a full derivation. If this local formula is off by a power, a missing boundary term, or a normalization of Γ_i, the density correction fails globally even though the reference-point argument itself is sound. The paper's own Supplementary Material raises doubts: Eq. (26) contains a Vandermonde factor Δ²(sin²x̂)/Δ²(x̂²) while the main text Eq. (8) has Δ²(x̂), an inconsistency not addressed. Also, the right environments Γ_i are defined with an implicit boundary condition (Γ_N appears to be I_D from the recurrence, but this is never stated), and the determinant power depends on that boundary. The extension from a single reference point to the whole manifold is logically valid only if the local metric formula holds at the reference point; it is not independently verified anywhere in the manuscript. The numerical figures test properties of the Metropolis-sampled distribution, not the determinant identity itself.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript considers the standard random matrix product state (RMPS) ensemble, formed by applying independent Haar-random unitaries at each site of a left-canonical MPS, and argues that this ensemble is biased with respect to the Fubini-Study volume on the MPS manifold. The central claim is Eq. (12): the Fubini-Study measure is related to the RMPS measure by dμ_FS = dμ_RMPS ∏_i |Γ_i|^{D(d-1)}, where Γ_i is the right environment at site i. The authors derive this from a local metric determinant, propose a Metropolis-Hastings sampler for the corrected measure, and study properties of the new ensemble, including a resolution of identity, entanglement asymmetry under spatial inversion, and an approximate Marchenko-Pastur spectrum with aspect ratio 1/(2d−1).","tokens_in":11027,"tokens_out":29058,"duration_ms":302500,"significance":"If fully established, the paper provides a parameter-free, gauge-invariant correction to a commonly used random-MPS ensemble. The predicted spatial-inversion asymmetry is a sharp and falsifiable claim, and the correction factor depends only on the entanglement spectrum, which makes it easy to apply. The paper also gives a concrete sampling algorithm and makes an analytic prediction for the typical entanglement spectrum. The central formula is elegant and the numerical investigations are appropriate. However, several load-bearing derivations are either missing or contain inconsistencies, so the manuscript requires substantial revision before the claims are fully supported.","major_comments":[{"comment":"The Jacobian printed in Eq. (8) has Δ²(ˆx) in the denominator, whereas the supplement's Eq. (26) has Δ²(ˆx²). Since the text defines ˆx² = x†x, the main-text expression gives J(x,x)→0 as x→0 for D>1, which would make the coordinate comparison at the origin singular; the supplement version is the correct one. This inconsistency affects the derivation of Eq. (12) as written and must be fixed by correcting Eq. (8) and any statement that depends on it.","section":"MPS Measure, Eq. (8) and Supplement Eq. (26)"},{"comment":"The central determinant identity |g| = ∏_i |Γ_i|^{D(d−1)} is asserted after \"standard MPS contraction calculations\" with no derivation. Because Eq. (12) is the main result, this is a load-bearing gap. Please provide a full derivation in the supplement or an explicit reference, and ideally a direct numerical check of the local metric determinant against the formula at several reference points. The skeptic's requested numerical verification is appropriate and would substantially strengthen the paper.","section":"MPS Measure, Eq. (10)"},{"comment":"The resolution-of-identity calculation contains an algebraic error. For a fixed state, Haar averaging over all local unitaries gives I_{d^N}/d^N, so the equality in Eq. (15) should read (1/d^N)∫dμ_FS(Ψ) Tr(|Ψ⟩⟨Ψ|) I_{d^N} = (1/d^N) I_{d^N} ∫dμ_FS(Ψ), and the coefficient in Eq. (16) is Z/d^N rather than Z. As written, the statement that the volume is proportional to the identity with coefficient Z is dimensionally inconsistent.","section":"Properties, Eqs. (14)–(16)"},{"comment":"The pointwise argument \"since we could have chosen any point as the reference\" is sound only if the exponential parameterization Eq. (6) is a local chart at every point of the gauge-fixed manifold. The paper asserts this via Ref. [43] but does not state the gauge condition under which the chart is valid, nor the boundary condition on the right environment (Γ_N, which appears to be I_D from the recurrence but is never declared). Please make these assumptions explicit and state that the formula extends to degenerate points by continuity; otherwise the global statement of Eq. (12) is not fully supported.","section":"MPS Measure, Eq. (6) and global extension"}],"minor_comments":[{"comment":"The proof of free independence of A and B is sketched rather than completed; the crucial step \"this guarantees...\" is an assertion. If this is to stand as a proof, more details are needed. Since the FS spectrum is described as approximate, the heuristic character should also be made explicit for the RMPS fixed point, which is currently presented as proven.","section":"Supplement Eq. (44)–(46)"},{"comment":"Equation (2) is incomplete: A[i] = U[i] 0 should define A[i] as the first D columns of U[i], and Eq. (3) is garbled and should read Γ_{i−1} = ∑_σ A^{[i]}_σ Γ_i (A^{[i]}_σ)†. Please fix the notation.","section":"Matrix product states, Eqs. (2)–(3)"},{"comment":"The acceptance ratio in Eq. (17) should specify that proposals are drawn from a symmetric distribution on the unitary U_j; if proposals are instead made in the x coordinates, the Jacobian J(x,x) must be included. Please clarify.","section":"Sampling, Eq. (17)"}],"recommendation":"major_revision","confidential_remarks":"The paper is likely of interest to the tensor-network and quantum-chaos communities. The central formula is plausible and potentially useful, but the derivation contains a concrete algebraic error in the resolution-of-identity step and an internal inconsistency in the Jacobian formula, and the key determinant identity is unproven. These issues are fixable within the scope of the manuscript, so I do not see grounds for rejection. I recommend major revision and ask the authors to provide the missing derivation and to correct the identified inconsistencies."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The central claim here is real and worth taking seriously: the usual sequentially generated RMPS ensemble is biased toward one end of the chain, and the unbiased Fubini-Study measure differs by a product over sites of |Γ_i|^{D(d-1)}. I don't know of this correction appearing in the earlier literature, so as a matter of novelty the paper is on solid ground. It also ships a Metropolis sampler, a resolution-of-identity argument, and a free-probability estimate that the typical entanglement spectrum shifts from Marchenko-Pastur with c = 1/d to roughly c' = 1/(2d−1). There are no fitted parameters in the central derivation, and the self-citations are motivational rather than load-bearing. That is genuine credit.\n\nThe soft spots are concentrated in one load-bearing place: Eq. (10), |g| = ∏|Γ_i|^{D(d−1)}, is stated after 'standard MPS contraction calculations' but never shown. The entire measure correction rests on that determinant formula, so leaving it as an assertion is a real gap. The supplement makes it worse by writing Δ²(sin²x̂)/Δ²(x̂²) in Eq. (26) while the main text Eq. (8) has Δ²(x̂); those are not the same object, and the paper does not address the mismatch. The boundary condition for the right environment (Γ_N is implicitly the identity) is also never stated, and the determinant power depends on that choice. The extension from a reference point to the whole manifold relies on the claimed surjectivity of the exponential parameterization; that is plausible and supported by a citation, but degenerate points deserve a comment.\n\nI want to push back on one reader complaint: 'the numerics lack error bars' is overstated. Fig. 1 reports means and standard deviations over 200 samples; Figs. 2 and 3 are CDFs and histograms where error bars would be nice but are not essential. The sampler is lightly specified — proposal distribution, burn-in scaling unknown — but the authors say so themselves. The free-independence proof is a sketch ('it can be easily checked' for d = 2, then 'can be easily extended'); that is a weakness but only in the auxiliary RMPS spectrum derivation, not in the main measure correction.\n\nNone of this kills the result. The cheapest fix is a direct numerical check of det g against ∏|Γ_i|^{D(d-1)} at a sample of random points; if that passes, the paper is largely right. I would send it to a serious referee and ask for that check plus a cleaner supplement. The paper is for tensor network practitioners and people studying typical entanglement in many-body systems; it should not be desk-rejected.","headline":"A genuinely new and plausible correction to the standard random MPS ensemble, but the load-bearing metric determinant is asserted rather than proven; deserves serious refereeing with a request for a full derivation or a direct numerical check.","tokens_in":11508,"tokens_out":1766,"would_cite":true,"duration_ms":21103,"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":"Sequential random matrix product states are not uniform on the MPS manifold; this paper derives and samples the unbiased Fubini-Study measure.","keywords":["matrix product states","Fubini-Study measure","random tensor networks","entanglement asymmetry","Haar random unitaries","Metropolis sampling","free probability","Marchenko-Pastur distribution"],"falsifier":"Compute the Fubini-Study metric determinant numerically at many points of the exponential coordinates for a small chain, including points where some Schmidt coefficient vanishes, and compare the local volume element with $\\prod_i |\\Gamma_i|^{D(d-1)}$; a mismatch on any set of positive measure would falsify Eq. (12). A complementary check is to generate Metropolis samples and test for exact left-right symmetry of the entanglement entropy profile.","tokens_in":10579,"feed_emoji":"🎲","tokens_out":15835,"duration_ms":147074,"temperature":0.7,"pith_summary":"This paper establishes that the standard ensemble of random matrix product states, built by placing independent Haar-random unitaries on every site, is not uniform with respect to the Fubini-Study measure that the manifold inherits from the full Hilbert space. The bias is visible as an asymmetry in entanglement entropy when the chain is reflected, an artifact of the left-canonical construction. The paper derives the exact conversion factor between the two measures, $d\\mu_{\\mathrm{FS}} = d\\mu_{\\mathrm{RMPS}}\\prod_i |\\Gamma_i|^{D(d-1)}$, where $D$ is the bond dimension, $d$ the physical dimension, and $\\Gamma_i$ the right-environment matrices that encode the Schmidt spectrum across each cut. It then constructs a Metropolis-Hastings sampler for the corrected measure and studies its properties, including a resolution of identity over MPS and a typical entanglement spectrum that differs from the sequentially generated ensemble. A sympathetic reader should care because statements about typical entanglement, chaos, and thermalization in tensor networks are only meaningful if the sampling measure is the natural unbiased one.","feed_headline":"Random MPS are biased; corrected measure restores uniformity","feed_subtitle":"Sequential Haar-random tensors skew entanglement; the correction reweights by right-environment determinants.","key_machinery":"The load-bearing object is the right-environment matrix $\\Gamma_i$, defined recursively from the right edge of the chain; its determinant appears in the Fubini-Study volume element as $|g| = \\prod_i |\\Gamma_i|^{D(d-1)}$, so it directly controls the measure correction. The second ingredient is the exponential parameterization $U(x) = U_0 \\exp\\!\\begin{pmatrix} 0 & -x^\\dagger \\\\ x & 0 \\end{pmatrix}$, which fixes the gauge and, as the paper asserts with a cited reference, covers the whole MPS manifold; the Haar measure over unitaries is then transformed into Lebesgue measure over the complex matrices $x$ with a Jacobian computed from random matrix theory. Comparing that Jacobian with the metric determinant gives the measure ratio, and makes the correction computable and samplable.","core_discovery":"The central claim is the measure relation (12): the sequentially generated RMPS measure and the unbiased Fubini-Study measure on the matrix product state manifold differ by a state-dependent factor, $$d\\mu_{\\mathrm{FS}} = d\\mu_{\\mathrm{RMPS}}\\prod_{i=1}^{N} |\\Gamma_i|^{D(d-1)},$$ with $\\Gamma_i$ the right environment at cut $i$, whose eigenvalues are the squared Schmidt coefficients of the bipartition across that cut. The derivation computes the metric determinant of the MPS manifold in a gauge-fixed exponential parameterization and compares it with the Jacobian of the Haar measure on the same coordinates; because the reference point can be chosen anywhere, the paper concludes the ratio holds globally. The direct consequence is that RMPS samples over-represent states with small right environments, producing a left-right asymmetry in entanglement entropy and a typical entanglement spectrum that differs from the unbiased ensemble. The paper also shows the corrected measure yields a resolution of identity proportional to the many-body identity, and provides an explicit Metropolis-Hastings algorithm with numerical verification.","pith_inferences":["Because the correction depends only on the Schmidt spectrum, any unbiased MPS sampling scheme in a different gauge must reproduce the same state-space volume; comparing central-gauge samples with FS samples would provide a direct numerical check.","The orientation bias found here should afflict other tensor-network ensembles built from independent local unitaries, and the same determinant-ratio strategy may extend naturally to loopless networks such as tree tensor networks, which the paper names as a plausible next step.","The paper leaves open three stated soft spots: its Marchenko-Pastur prediction for the FS spectrum is a first-order, non-self-consistent approximation; the Metropolis burn-in scaling with $N$ and $D$ is unknown; and an efficient sequential (rather than Markov-chain) sampler does not yet exist.","If the resolution of identity holds, the FS measure gives a bona fide classical phase space for the time-dependent variational principle, which would allow projected Lyapunov spectra and statistical-mechanics treatments of entanglement dynamics to be defined; this is an application the paper anticipates but does not carry out."],"forward_implications":["Entanglement profiles drawn from the unbiased FS ensemble are symmetric under spatial inversion, while RMPS profiles are not; shifting the orthogonality center to the middle does not fix this but introduces a discontinuity.","The determinant factor repels the eigenvalues of the right environment away from zero, so typical FS samples carry more entanglement entropy than RMPS samples at the same bond dimension.","The FS measure admits a resolution of identity, $\\int d\\mu_{\\mathrm{FS}}(\\Psi)|\\Psi\\rangle\\langle\\Psi| = Z\\,\\mathbb{I}_{d^N}$ with $Z = \\int \\prod_i dU_i\\, e^{-D(d-1)\\sum_i \\log|\\Gamma_i|^{-1}}$, turning the MPS manifold into a usable semiclassical phase space.","The typical right-environment spectrum in the FS ensemble is approximately Marchenko-Pastur with aspect ratio $1/(2d-1)$ instead of $1/d$, with deviations visible near the spectral edge.","A Metropolis-Hastings sweep over the chain mixes rapidly in the tested cases, making the unbiased ensemble practically samplable."],"supporting_citations":[{"why":"Supplies the differential-geometric description of the MPS manifold and its metric, used to obtain the Fubini-Study volume element.","marker":"[40]"},{"why":"Cited for the claim that the exponential parameterization globally covers the gauge-fixed MPS manifold, the premise for extending the local measure ratio to all points.","marker":"[43]"},{"why":"Provides the random-matrix Jacobian technology that converts the Haar measure on unitaries into Lebesgue measure on the exponential coordinates.","marker":"[44]"},{"why":"Defines the Marchenko-Pastur law used as the RMPS baseline and as the starting approximation for the FS entanglement spectrum.","marker":"[45]"}],"fun_headline_variants":["RMPS bias revealed; unbiased MPS measure via Metropolis","Unbiased MPS ensemble: fix for Haar-random bias","Sequential RMPS bias fixed; unbiased measure via Metropolis","Uniform MPS ensemble: reweighting removes left-right bias","Haar-random MPS biased; new measure fixes it"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole argument stands or falls on whether every gauge-fixed matrix product state can be reached by the exponential coordinate patch used in the calculation, without missing regions or singularities, so that the local volume ratio computed at one reference point is valid everywhere.","fun_headline_variants_meta":{"raw":{"variants":["RMPS bias revealed; unbiased MPS measure via Metropolis","Unbiased MPS ensemble: fix for Haar-random bias","Sequential RMPS bias fixed; unbiased measure via Metropolis","Uniform MPS ensemble: reweighting removes left-right bias","Haar-random MPS biased; new measure fixes it"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000827,"raw_usage":{"total_tokens":3589,"prompt_tokens":893,"completion_tokens":2696,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":509,"completion_tokens_details":{"reasoning_tokens":2612}},"tokens_in":509,"tokens_out":2696,"duration_ms":18844,"temperature":1.0,"reasoning_tokens":2612,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T04:52:09.654223+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the Fubini-Study metric determinant numerically at many points of the exponential coordinates for a small chain, including points where some Schmidt coefficient vanishes, and compare the local volume element with $\\prod_i |\\Gamma_i|^{D(d-1)}$; a mismatch on any set of positive measure would falsify Eq. (12). A complementary check is to generate Metropolis samples and test for exact left-right symmetry of the entanglement entropy profile.","supporting_citations":[{"cited_title":"Wouters, N","cited_arxiv_id":null,"evidence_quote":"Cited for the claim that the exponential parameterization globally covers the gauge-fixed MPS manifold, the premise for extending the local measure ratio to all points."},{"cited_title":"Mehta, Random Matrices , Pure and Applied Mathe- matics (Academic Press, 2004)","cited_arxiv_id":null,"evidence_quote":"Provides the random-matrix Jacobian technology that converts the Haar measure on unitaries into Lebesgue measure on the exponential coordinates."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the Marchenko-Pastur law used as the RMPS baseline and as the starting approximation for the FS entanglement spectrum."}],"review_version":1}