{"id":"32deec0b-7372-43fd-aec2-12a483c8b8ad","arxiv_id":"2509.21846","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Exact closed-form averages of relative entropy between independent random density matrices are derived for Hilbert-Schmidt and Bures-Hall ensembles, including cross-ensemble cases.","lead":"This paper derives exact formulas for the average relative entropy of two independent random quantum states, for the Hilbert-Schmidt and Bures-Hall ensembles and for cross-ensemble pairs. The formulas are exact, not asymptotic, and they complement and extend a previous leading-order result based on the replica method.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Bures-Hall formulas rest on an entropy formula that violates the maximum-entropy bound; Eq. (64) has a sign error.","rationale":"The factorization E[tr(ρ ln σ)] = (1/m)E[tr(ln σ)] in Eq. (53) is secure for the unitarily invariant Hilbert-Schmidt and Bures-Hall ensembles: eigenvector unitaries are Haar-distributed and independent of the eigenvalues, so the Weingarten average is exact. The reader's identified weakest assumption is therefore not where the argument breaks. Instead, the Bures-Hall branch contains an internal contradiction: Eq. (7) violates the universal bound S(ρ) ≤ ln m already at m=2, n=2. Since Propositions 2 and 3 explicitly use Eq. (7), they cannot be correct as stated. The sign error in Eq. (64) is concrete and testable. This is a stronger failure than the m=1 boundary issue noted by the reader, and it changes the central formulas rather than merely adding a caveat.","tokens_in":9027,"tokens_out":44235,"duration_ms":372615,"concrete_test":"Re-derive Eq. (64) directly from Eq. (26) by differentiating ln C_BH with respect to α and check the sign of the mψ₀(m(m+2α)/2) term. Independently, run a Monte Carlo for m=2, n=2 Bures-Hall states, e.g. via Eq. (27), estimating E[−tr(ρ ln ρ)] and E[D(ρ||σ)]; if the entropy estimate is near 0.2 nats while Eq. (7) gives 0.886 nats, the Bures-Hall propositions are numerically falsified.","verdict_should_be":"REJECT","load_bearing_attack":"Propositions 2 and 3 both use Eq. (7), the claimed mean von Neumann entropy of the Bures-Hall ensemble. As printed, Eq. (7) is impossible: for m=2, n=2 it gives ψ₀(3)−ψ₀(3/2) ≈ 0.8863, but every 2×2 density matrix has entropy at most ln2 ≈ 0.6931. Hence Eq. (7) cannot be a mean entropy, and Eqs. (15), (16), (17), and (18) inherited from it are suspect. The likely root cause is in the derivative of the Bures-Hall normalization constant: differentiating Eq. (26) gives d ln C_BH/dα = −2m ln2 + mψ₀(m(m+2α)/2) + 2Σψ₀(i+2α) − Σψ₀(i+α), but Eq. (64) puts a minus sign on the mψ₀ term. This is not a boundary artifact: it affects the generic m ≥ 2 regime and invalidates the claimed exact Bures-Hall results.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes exact closed-form formulas for the average quantum relative entropy E[D(ρ||σ)] for two independent random density matrices drawn from the Hilbert–Schmidt (HS) and Bures–Hall (BH) ensembles, including the mixed HS-vs-BH case. The main mathematical step is a unitary-integral factorization E[tr(ρ ln σ)] = (1/m)E[tr(ln σ)] (Eq. 53), which reduces the problem to known single-state entropy formulas and to derivatives of the ensemble normalization constants. The results are stated as Propositions 1–3 with limiting corollaries, and are compared with numerical simulations.","tokens_in":9242,"tokens_out":11952,"duration_ms":94381,"significance":"The Hilbert–Schmidt part appears correct and useful: the factorization argument via zonal polynomials and Weingarten calculus is clean, Proposition 1 passes the m=1 consistency check, and its large-dimension limit reproduces the known asymptotic result of Kudler-Flam. The Bures–Hall part, however, is not currently reliable: two independent checks (the maximum-entropy bound and the deterministic m=1 limit) fail, and the normalization derivative in Eq. (64) has a sign error. If the BH computations are corrected, the factorization method could still be a valuable contribution, but the paper as it stands does not establish the BH formulas.","major_comments":[{"comment":"The derivative of ln C_BH printed in Eq. (64) has the wrong sign on the term mψ0(m(m+2α)/2). Starting from Eq. (26), d ln C_BH/dα = −2m ln2 + mψ0(m(m+2α)/2) + 2Σψ0(i+2α) − Σψ0(i+α), whereas Eq. (64) uses −mψ0(...). This sign error propagates into Eq. (15), and through Propositions 2 and 3 into Eqs. (16), (17), (18), and (69).","section":"Proof of Proposition 2, Eq. (64)"},{"comment":"The Bures–Hall mean entropy quoted in Eq. (7) cannot be correct: for m=2, n=2 it gives ψ0(3)−ψ0(3/2) ≈ 0.8863, which exceeds the absolute upper bound ln2 ≈ 0.6931 for any 2×2 density matrix. Since Propositions 2 and 3 are built on Eq. (7), their claimed exact formulas and all BH limiting formulas are unsubstantiated as they stand.","section":"Section 1, Eq. (7)"},{"comment":"At m=1, n1=n2=1 both ensembles give the deterministic state ρ=σ=1, so D(ρ||σ)=0. Eq. (15) evaluates to approximately 0.386 and Eq. (17) to approximately −0.614, contradicting both this limiting value and the nonnegativity (3). This internal inconsistency confirms that the Bures–Hall computation of E[tr(ln σ)] or the entropy input is wrong, independent of the maximum-entropy bound.","section":"Proposition 2 and Proposition 3, m=1 special case"}],"minor_comments":[{"comment":"The symbol \"ρ_BS\" appears where \"σ_BH\" or \"ρ_BH\" seems intended; please correct this notation.","section":"Section 1, text after Eq. (10)"},{"comment":"The captions say the plots compare exact formulas with simulations, but they do not state the values of n1, n2 (or c1, c2) used for each curve; please add this information for reproducibility.","section":"Figures 1 and 2"},{"comment":"Reference [13] is cited as arXiv:2502.05371; if it has been published, please provide the journal reference.","section":"Reference [13]"},{"comment":"The notation n_i = m + α_i is first used in the propositions without an explicit definition in the introduction; please define it before presenting Proposition 1.","section":"Around Eq. (11)"}],"recommendation":"major_revision","confidential_remarks":"The Bures–Hall defects are load-bearing but appear fixable: correct the sign in Eq. (64), replace Eq. (7) with the correct mean BH entropy, and rederive Propositions 2–3. The editor may wish to ask the author to verify the revised formulas against the m=1 deterministic limit and the m=2, n=2 maximum-entropy bound before further review."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: the Hilbert-Schmidt half of this paper is real; the Bures-Hall half is not. Proposition 1, the exact average relative entropy for two independent HS states, checks out. I spot-checked the m=1 limit and a nontrivial m=n=2 case; it behaves. The limiting formula in Corollary 1 recovers Kudler-Flam's asymptotic as it should, which is the right sanity check. The factoring step E[tr(ρ ln σ)] = (1/m)E[tr(ln σ)] via unitary invariance is clean and correct, and the derivative of the HS normalization constant is fine. That single result is a genuine closed-form upgrade over the existing asymptotic and is worth publishing on its own.\n\nThe same cannot be said for Propositions 2 and 3. Equation (7), the input BH mean entropy, is impossible as printed: for m=2, n=2 it gives ψ0(3)−ψ0(3/2) ≈ 0.886, while every 2×2 density matrix has entropy at most ln2 ≈ 0.693. No interpretation of the parameters fixes that. The BH formulas inherited from it fail the trivial m=1 check, where both states are deterministic and D=0 but Eq. (15) gives roughly 0.3865. Independent of that, the derivative calculation in Eq. (64) has a sign error: differentiating the logarithm of (26) gives +m ψ0(m(m+2α)/2), not −m ψ0(...). So the BH results are not missing a caveat—they are wrong, and the wrongness is load-bearing.\n\nThe paper's machinery and the HS result are unaffected. A reader working on random-state distinguishability or quantum hypothesis testing will want the HS formula; the BH claims should be ignored until the author tracks down the correct BH entropy input (the cited literature may not actually say what Eq. (7) says) and redoes Eqs. (15)–(18).\n\nFor peer review: send it out, but with instructions for major revision. One advertised contribution is genuinely new and correct; two are currently false. A competent referee can verify the HS part quickly and give the author the algebra errors on the BH side. I would not cite the BH formulas in their current form, but the paper deserves referee time.\n\nOverall: conditional accept in spirit—major revision, not rejection.","headline":"Hilbert-Schmidt half is solid and new; the Bures-Hall half rests on an impossible entropy formula and a sign error—worth a major-revision round, not a desk reject.","tokens_in":9726,"tokens_out":7713,"would_cite":false,"duration_ms":65697,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["60B20","81P45","94A17"],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper derives exact closed-form formulas for the average relative entropy of pairs of independent random states from the Hilbert-Schmidt and Bures-Hall ensembles.","keywords":["relative entropy","random quantum states","Hilbert-Schmidt ensemble","Bures-Hall ensemble","entanglement entropy","unitary integration","zonal polynomials","digamma function"],"falsifier":"Draw many independent pairs of random density matrices for a small explicit case, such as $m=2$, $n_1=3$, $n_2=4$, from the Hilbert-Schmidt and Bures-Hall constructions, and compare the sample average of $D(\\rho\\|\\sigma)$ with (11), (15), and (17); a disagreement beyond Monte Carlo error would refute the formulas. A sharper test isolates the factorization by checking whether the sample average of $\\operatorname{tr}(\\rho\\ln\\sigma)$ equals $(1/m)$ times the sample average of $\\operatorname{tr}(\\ln\\sigma)$ for independent pairs, and that this equality breaks when the eigenbasis of $\\sigma$ is forced to coincide with that of $\\rho$.","tokens_in":8805,"feed_emoji":"🎲","tokens_out":12459,"duration_ms":92257,"temperature":0.7,"pith_summary":"Relative entropy measures how well one quantum state can stand in for another, so its average over random states quantifies typical distinguishability. This paper establishes exact, finite-size formulas for that average when the two states are drawn independently from the Hilbert-Schmidt ensemble, the Bures-Hall ensemble, or one from each. The core step is a factorization: after averaging over the unitary relating the two eigenbases, the cross term $E[\\operatorname{tr}(\\rho\\ln\\sigma)]$ collapses to $(1/m)E[\\operatorname{tr}(\\ln\\sigma)]$, reducing a two-state problem to known single-state entropy averages plus a log-determinant average. That produces the closed forms in Propositions 1–3, valid for arbitrary subsystem dimension $m$ and parameters $n_1=m+\\alpha_1$, $n_2=m+\\alpha_2$, with large-dimension limits that complement the earlier asymptotic replica-method result. Exact formulas of this kind give finite-size predictions for quantum hypothesis testing and for studies of eigenstate thermalization.","feed_headline":"Exact average relative entropy of random states derived","feed_subtitle":"Finite-size formulas cover Hilbert-Schmidt and Bures-Hall ensembles and sharpen the asymptotic replica estimates.","key_machinery":"The load-bearing identity is the unitary-integral factorization\n$$\\int_{U(m)}\\operatorname{tr}\\bigl(\\Lambda_\\rho U\\ln\\Lambda_\\$\\sigma$ U^\\dagger\\bigr)\\,dU=\\frac{1}{m}\\operatorname{tr}(\\ln\\$\\sigma$),$$\nwhich follows from the zonal-polynomial integral $\\int_{U(m)}C_\\kappa(XUYU^\\dagger)\\,dU=C_\\kappa(X)C_\\kappa(Y)/C_\\kappa(I_m)$ and, for this first moment, can also be obtained by Weingarten calculus. It reduces the cross term in relative entropy to a single-state log-determinant average. The remaining task is to compute $E[\\operatorname{tr}(\\ln\\sigma)]$ for each ensemble by differentiating the normalization constants of the Hilbert-Schmidt and Bures-Hall densities with respect to $\\alpha$, which produces the digamma-function combinations in (11), (15), and (17).","core_discovery":"On the paper's own terms, the discovery is that the average relative entropy of two independent random density matrices is controlled by a factorization identity rather than by the joint spectrum of the pair. Writing $\\rho=V\\Lambda_\\rho V^\\dagger$ and $\\sigma=W\\Lambda_\\sigma W^\\dagger$, the unitary invariance of the Hilbert-Schmidt and Bures-Hall ensembles lets the average over $U=V^\\dagger W$ be performed first, giving $E[\\operatorname{tr}(\\rho\\ln\\sigma)]=(1/m)E[\\operatorname{tr}(\\ln\\sigma)]$. Combined with the exact mean entanglement entropy of a single random state from each ensemble, this yields explicit digamma-function expressions: Proposition 1 for Hilbert-Schmidt versus Hilbert-Schmidt, Proposition 2 for Bures-Hall versus Bures-Hall, and Proposition 3 for Bures-Hall versus Hilbert-Schmidt, with the complementary Hilbert-Schmidt versus Bures-Hall case recorded as well. The formulas are exact for arbitrary $m$, $n_1=m+\\alpha_1$, and $n_2=m+\\alpha_2$, and their large-dimension limits recover the previously known asymptotic behavior to leading order.","pith_inferences":["Because the factorization step uses only unitary invariance and independence, the same reduction $E[\\operatorname{tr}(\\rho\\ln\\sigma)]=(1/m)E[\\operatorname{tr}(\\ln\\sigma)]$ should hold for any two independent states from other unitarily invariant ensembles, so analogous exact formulas could be derived for fermionic Gaussian states once the relevant log-determinant averages are known.","The zonal-polynomial machinery used here is suited to higher powers of $\\operatorname{tr}(\\rho\\ln\\sigma)$, which suggests the method can be pushed to exact variance or higher-cumulant formulas for relative entropy, not just the mean.","If two states are drawn with correlated eigenbases or from a non-unitarily-invariant ensemble, the factorization fails; measuring that failure could serve as a quantitative probe of how eigenvector alignment changes typical distinguishability, a testable extension of the independent-state setting."],"forward_implications":["For any fixed finite subsystem dimension $m$ and bipartite parameters $n_1,n_2$, the average relative entropy can now be evaluated exactly instead of approximated, removing finite-size error in regimes where only the asymptotic replica formula was available.","In the limit $m\\to\\infty$ with $n_i/m=c_i$ fixed, the limiting formulas (13), (16), and (18) show that the average relative entropy depends only on the ratios $c_1,c_2$ and decreases monotonically as either parameter grows, the maximum occurring at $c_1=c_2=1$.","Comparing the two same-ensemble formulas shows that the Bures-Hall average relative entropy exceeds the Hilbert-Schmidt value for the same parameters, which the paper attributes to the wider spectral width of the Bures-Hall ensemble.","The cross-ensemble formula quantifies how much a Bures-Hall state typically differs from a Hilbert-Schmidt model, and its asymmetry under exchanging parameters indicates different robustness of the two ensembles as reference models."],"supporting_citations":[{"why":"This supplies the mean entanglement entropy formula for the Hilbert-Schmidt ensemble, which is the first term in (5) and in Proposition 1.","marker":"[20]"},{"why":"This provides the proof of the Hilbert-Schmidt mean entropy conjecture that underlies the single-state term used for Proposition 1.","marker":"[6]"},{"why":"This supplies spectral densities and average entropies for the Bures-Hall ensemble, from which the single-state entropy term in Proposition 2 is taken.","marker":"[22]"},{"why":"This proves the Bures-Hall average entropy conjectures used as the single-state term in Proposition 2.","marker":"[26]"},{"why":"This gives the prior asymptotic replica-method formula for equal-dimension Hilbert-Schmidt relative entropy, which the exact result complements and reproduces in the large-dimension limit.","marker":"[15]"},{"why":"This provides the zonal-polynomial unitary integral identity used to derive the factorization (39)-(43).","marker":"[7]"},{"why":"This contains an earlier derivation of the zonal-polynomial unitary integral that the factorization step relies on.","marker":"[14]"}],"fun_headline_variants":["Exact relative entropy for random quantum states","Random-state relative entropy: exact finite-size formulas","Factorization unlocks exact relative entropy","Arbitrary-dimension exact relative entropy derived","Beyond replicas: exact relative entropy formulas"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The chain of proofs depends on the two random states being independent and on each state's eigenvector basis being uniformly random according to the invariant unitary measure and statistically independent of its eigenvalues; if either ensemble lost its unitary invariance, or if the states were not independent, the identity $E[\\operatorname{tr}(\\rho\\ln\\sigma)]=(1/m)E[\\operatorname{tr}(\\ln\\sigma)]$ would fail and all three propositions would fall.","fun_headline_variants_meta":{"raw":{"variants":["Exact relative entropy for random quantum states","Random-state relative entropy: exact finite-size formulas","Factorization unlocks exact relative entropy","Arbitrary-dimension exact relative entropy derived","Beyond replicas: exact relative entropy formulas"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000365,"raw_usage":{"total_tokens":1931,"prompt_tokens":878,"completion_tokens":1053,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":494,"completion_tokens_details":{"reasoning_tokens":988}},"tokens_in":494,"tokens_out":1053,"duration_ms":7908,"temperature":1.0,"reasoning_tokens":988,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T15:47:08.303030+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Draw many independent pairs of random density matrices for a small explicit case, such as $m=2$, $n_1=3$, $n_2=4$, from the Hilbert-Schmidt and Bures-Hall constructions, and compare the sample average of $D(\\rho\\|\\sigma)$ with (11), (15), and (17); a disagreement beyond Monte Carlo error would refute the formulas. A sharper test isolates the factorization by checking whether the sample average of $\\operatorname{tr}(\\rho\\ln\\sigma)$ equals $(1/m)$ times the sample average of $\\operatorname{tr}(\\ln\\sigma)$ for independent pairs, and that this equality breaks when the eigenbasis of $\\sigma$ is forced to coincide with that of $\\rho$.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"This supplies the mean entanglement entropy formula for the Hilbert-Schmidt ensemble, which is the first term in (5) and in Proposition 1."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"This provides the proof of the Hilbert-Schmidt mean entropy conjecture that underlies the single-state term used for Proposition 1."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"This supplies spectral densities and average entropies for the Bures-Hall ensemble, from which the single-state entropy term in Proposition 2 is taken."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"This proves the Bures-Hall average entropy conjectures used as the single-state term in Proposition 2."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"This gives the prior asymptotic replica-method formula for equal-dimension Hilbert-Schmidt relative entropy, which the exact result complements and reproduces in the large-dimension limit."},{"cited_title":"Princeton University Press, Princeton (2010)","cited_arxiv_id":null,"evidence_quote":"This provides the zonal-polynomial unitary integral identity used to derive the factorization (39)-(43)."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"This contains an earlier derivation of the zonal-polynomial unitary integral that the factorization step relies on."}],"review_version":1}