{"id":"063d6ad9-4432-4977-9e18-0f35c9ee1bcc","arxiv_id":"2506.06663","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A closed-form expression for the third cumulant (skewness) of von Neumann entropy over the Bures-Hall random states is obtained and validated numerically.","lead":"Researchers derived an exact formula for the third statistical moment, the skewness, of entanglement entropy of random bipartite quantum states drawn from the Bures-Hall ensemble. This lets them approximate the full distribution of entanglement more accurately than the usual Gaussian approximation.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The new summation identities (A.13)-(A.27) are under-verified and one adjacent identity appears corrupted; Proposition 1 rests on them.","rationale":"The paper's main result is an exact closed-form expression for the third cumulant, obtained through a long symbolic computation. The authors themselves identify the new summation identities as the key step, and the reader's weakest-assumption analysis flags the correctness of these identities as the load-bearing premise. Our stress-test finds no decisive internal inconsistency, but the apparent corruption in the neighboring identity (A.3) and the mismatch between the stated validity condition a>m and the parameters used in the numerical validation increase the correctness risk. A direct numerical check of the identities would settle the matter quickly and cheaply. Until such a check is performed, the conditional acceptance recommended by the reader remains the appropriate verdict; no change is needed.","tokens_in":23591,"tokens_out":8839,"duration_ms":88375,"concrete_test":"Verify numerically, to high precision, every new identity (A.13)-(A.27) for small m and for the parameter values used in Figure 2, e.g., m=3 and m=4 with n=m (α=-1/2), n=m+1 (α=1/2), n=2m, and n=3m, computing both sides and checking equality. In particular, test a values below m, such as a=m-1/2 and a=m-1, as well as half-integer a. If any identity fails, the cancellation behind (15) is invalid; if all hold, the central concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim (15) is established only if the twelve new summation identities (A.13)-(A.27) are all correct: Section 3.3 states that after processing the new anomalies with these identities, 'all anomalies cancel completely', leading to (111) and then (15). This makes the correctness of these identities the load-bearing premise. Only (A.19) is derived; the rest are asserted. Confidence in them is undercut by an apparent transcription error in the adjacent identity (A.3): its right-hand side contains a parameter a absent from the left-hand side and the term ψ0(c+1)ψ0(c+1) in place of the expected ψ0(c+1)ψ0(b+1). If any of the longer identities has a similar error, the cancellations leading to (111) and hence (15) fail. The stated condition 'a>m for (A.13)-(A.21)' also does not cover parameters used in the paper's own numerical checks: for n=m one has a=α+m=m−1/2 and a=2α+m=m−1, both below m. No analytic continuation or separate treatment is given, so the derivation as written does not cover these cases.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper derives an exact closed-form expression for the third cumulant (skewness) of the von Neumann entropy of bipartite pure states drawn from the Bures-Hall ensemble. The main result, Proposition 1, Eq. (15), expresses kappa_3 as a combination of polygamma functions with rational coefficients given in Eqs. (16)-(17). The derivation proceeds through a moment-cumulant conversion from the unconstrained Bures-Hall ensemble, Pfaffian correlation functions, Mellin transforms of Meijer G-functions, and a systematic cancellation of finite polygamma sums, with twelve new summation identities stated in Appendix A. The paper also provides the corresponding Edgeworth approximation (21), numerical comparisons for m = 3 to 12, and an asymptotic analysis showing kappa_3 / kappa_2^{3/2} = Theta(1/n), which supports the Gaussianity conjecture for the entropy distribution in the limit (23).","tokens_in":23815,"tokens_out":10487,"duration_ms":93634,"significance":"If Proposition 1 is correct, the result is a substantive exact contribution to random matrix theory and quantum information: it provides the first exact skewness formula for Bures-Hall entanglement entropy, improves the normal approximation to the entropy distribution, and gives concrete asymptotic evidence for the conjectured Gaussian limit. The derivation is parameter-free and builds on established kernel representations and previously derived first and second cumulants, which are strong points. The numerical agreement for m = 3 to 12 is encouraging. However, the central claim rests on twelve new summation identities, of which only one is derived in detail and none are independently machine-checked; because a single error in any of these identities would invalidate the cancellation leading to Eq. (111) and hence Eq. (15), this is the key point requiring verification.","major_comments":[{"comment":"The identity (A.3) as printed is corrupted: its right-hand side contains a parameter a that does not appear on the left-hand side, and the term psi_0(c+1) psi_0(c+1) should almost certainly be psi_0(c+1) psi_0(b+1) by symmetry. Since (A.3) is used in Section 3.3, e.g., in transforming the sums in Eqs. (104)-(107), this typo must be corrected and the corrected identity verified.","section":"Appendix A, Eq. (A.3)"},{"comment":"Only identity (A.19) is derived in detail; the other new identities (A.13)-(A.18), (A.20)-(A.21), (A.24), and (A.26)-(A.27) are asserted to follow by analogous re-summation or differentiation. Because the cancellation of all anomalies leading to (111), and hence the main result (15), depends on all of these identities, the authors should provide either complete derivations for all of them or an independent machine-verifiable proof, such as symbolic verification with m kept symbolic.","section":"Section 3.3 and Appendix A"},{"comment":"The stated condition 'a > m for (A.13)-(A.21)' is incompatible with the paper's own numerical checks: for n = m, one has alpha = -1/2 and the parameter values a = alpha + m = m - 1/2 and a = 2alpha + m = m - 1 are both less than m. No analytic continuation or special-case treatment is provided for these values, so the derivation as written does not cover the n = m cases displayed in Figure 2. Please supply the missing argument, even if the final rational formula extends by analyticity.","section":"Appendix A, parameter conditions before (A.1)"},{"comment":"The definitions of I_B^(3) and I_B^(4) are printed with identical integrands, both f^2(x) f(y) K_00(x,y) K_11(x,y) dxdy. The subsequent computation using B_{3,4} in Eqs. (89)-(92), with different derivative orders in beta_1 and beta_2, treats these two integrals as different, namely with f^2(x) f(y) versus f(x) f^2(y) (or an analogous kernel-argument swap). Please correct the definition of I_B^(4), because I_B enters the cumulant expression (50).","section":"Eqs. (57)-(58)"}],"minor_comments":[{"comment":"The quadruple-sum expression for B_{3,4}(beta_1,beta_2) would be easier to check if the summation ranges and Gamma-function arguments were typeset with the same level of clarity as the surrounding equations; currently the long inline fractions make verification difficult.","section":"Section 3.2, Eq. (90)"},{"comment":"The condition for (A.4) is stated as a >= m, but the identity is later used in derivations where a may equal m; please confirm whether the stated condition is necessary or whether the identity also holds at a = m by continuity.","section":"Appendix A, Eq. (A.4)"},{"comment":"The vertical axis is described as a log-linear plot, but the caption and text do not specify whether the ordinate is log |kappa_3| or log kappa_3; since kappa_3 appears to be negative in the plotted regime, please clarify the scale and the sign convention.","section":"Figure 2"},{"comment":"Reference [15] is cited as an arXiv preprint (2502.05371); since it is used for the 'summation-free framework' and the anomaly terminology, please update the reference if a journal version has appeared by the time of publication.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The core derivation is plausible and the numerical checks are reassuring, but the paper's own reliability hinges on a list of new summation identities that are largely unverified and one of which, (A.3), is visibly corrupted. The duplicated definitions in (57)-(58) and the mismatch between the stated a > m condition and the n = m numerics reinforce the need for a careful revision pass over Appendix A and Section 3.3. I would encourage the editor to ask for either full derivations or a reproducible symbolic-verification file for the new identities before acceptance."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The short version: the result is new and worth knowing. Proposition 1 gives an exact closed form for the third cumulant of von Neumann entropy over the Bures-Hall ensemble, with rational coefficients a1 and a2. The mean and variance were known; the skewness was not. The paper also shows the standardized third cumulant decays as Θ(1/n), strengthening the case for the Gaussianity conjecture, and it uses the skewness to improve the normal approximation via an Edgeworth correction. The numerics for m=3..12 at n=m, 2m, and 3m match the formula. That is real evidence.\n\nWhat I like: the derivation is parameter-free. No fitting, no invented entities. The moment-cumulant conversion takes the known first two cumulants as inputs, not as adjustable constants. The paper is honest about the difficulty, flagging the new anomalies and the twelve summation identities as the key step. The citation pattern is fine—prior works are the natural chain of the cumulant program, and the self-citations are to the exact variance and mean this paper builds on.\n\nSoft spots, in proportion. The twelve new summation identities (A.13)–(A.27) are the load-bearing wall. Only (A.19) is derived; the rest are asserted to follow from the same re-summation technique. That is a lot to take on faith for a closed form. The adjacent identity (A.3) has an obvious transcription error: an undefined parameter a on the right-hand side and ψ0(c+1)ψ0(c+1) where ψ0(c+1)ψ0(b+1) is expected. That kind of slip makes me want the other identities checked, not because I think they are wrong, but because the appendix is where typos will hide. Also, the stated condition a>m for (A.13)–(A.21) does not cover the paper's own numerical checks with n=m, where a=α+m and a=2α+m are below m. No analytic continuation is mentioned. This is a genuine gap in the written derivation, even though the numerical agreement suggests the identities survive by continuity.\n\nNet: I think the main claim is probably correct, and the numerics are reassuring, but the paper is not referee-ready as is. It needs a pass where the appendix identities are either fully derived or checked symbolically and the parameter conditions are stated correctly. If the identities check out, this is a solid contribution for the random-matrix/entanglement subfield.\n\nWho it is for: people working on exact cumulants of entanglement entropy over structured random-state ensembles, and anyone using Edgeworth approximations to entropy distributions. It deserves a serious referee despite the needed revision. If I were the editor, I would send it out with a referee who can verify the appendix with a computer algebra system.","headline":"A genuinely new third-cumulant result for Bures-Hall entropy, with a credible numerical check, but the load-bearing appendix of summation identities needs real verification before the closed form is trusted.","tokens_in":24310,"tokens_out":2818,"would_cite":true,"duration_ms":30183,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["15B52","60B20","33B15","81P40"],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper derives an exact closed-form expression for the third cumulant (skewness) of von Neumann entanglement entropy over the Bures-Hall ensemble of random quantum states.","keywords":["von Neumann entropy","Bures-Hall ensemble","skewness","third cumulant","random matrix theory","polygamma functions","summation identities","entanglement entropy"],"falsifier":"Evaluate the twelve identities (A.13)--(A.27) by direct numerical summation for several small integer values of $m$ and parameters $a, b, c$; any mismatch disproves the cancellation. Independently, simulate Bures-Hall random states for a few pairs $(m,n)$ at high precision and compare the sample third cumulant against Proposition 1.","tokens_in":23398,"feed_emoji":"🎲","tokens_out":6159,"duration_ms":52485,"temperature":0.7,"pith_summary":"The paper establishes an exact closed-form expression for the third cumulant of von Neumann entanglement entropy over the Bures-Hall ensemble of random quantum states. Earlier work had closed forms for the mean and variance; this paper supplies the skewness, showing that standardized entropy has skewness of order $\\Theta(1/n)$ as both subsystem dimensions grow with fixed ratio. The exact third cumulant feeds a skewness-corrected Edgeworth approximation to the entropy distribution, which matches simulations more closely than the plain Gaussian. This provides further quantitative evidence toward the conjecture that the standardized entropy converges to a standard normal.","feed_headline":"Third cumulant of Bures-Hall entropy is now exact","feed_subtitle":"A closed-form third cumulant sharpens the entropy distribution and supports the Gaussianity conjecture.","key_machinery":"The argument runs on the Bures-Hall measure's Pfaffian point-process structure, inherited from the Cauchy-Laguerre biorthogonal ensemble. Cumulants of the linear statistic $S$ are expressed through the three-point correlation functions of this Pfaffian process, and the integrals are evaluated via Meijer G-function kernel representations and Mellin transforms. The load-bearing combinatorial step is a set of twelve new summation identities (A.13)--(A.27) that relate finite sums of products of polygamma functions; these identities transform all unsimplifiable single sums ('anomalies') into a small set of basic sums that cancel exactly, leaving the compact closed form. The re-summation technique, which recasts a finite sum by iterating a difference recurrence, is the systematic tool used to derive the identities.","core_discovery":"The central discovery is a summation-free closed form for the third cumulant, given in Proposition 1: $\\kappa_3 = \\psi_2\\left(mn - \\frac{m^2}{2} + 1\\right) + a_1 \\psi_2\\left(\\frac{n+1}{2}\\right) + a_2 \\psi_1\\left(\\frac{n+1}{2}\\right)$, with rational coefficients $a_1, a_2$ displayed in equations (16)--(17). The formula is proven by computing the three-point correlation integrals of the Pfaffian point process associated with the Bures-Hall measure, converting the moment to a large collection of nested finite sums, and then showing that all 'anomalous' single sums cancel exactly once twelve new summation identities are applied. The result implies skewness $\\gamma_1 = \\kappa_3 / \\kappa_2^{3/2}$ and, in the limit $m, n \\to \\infty$ with $m/n = c \\in (0,1]$, yields $\\kappa_3^{(X)} = \\Theta(1/n)$, consistent with the Gaussianity conjecture.","pith_inferences":["Inference: The same summation-and-cancellation scheme likely extends to the fourth cumulant (kurtosis), though the number of anomalies will grow further; the techniques here may accelerate that computation.","Inference: The twelve identities (A.13)--(A.27) may be of independent combinatorial interest, as closed-form evaluations of families of harmonic-type sums with polygamma products.","Inference: The Edgeworth approximation with the exact skewness could be used as a statistical test for whether a measured bipartite state is typical of the Bures-Hall ensemble.","Inference: The $\\Theta(1/n)$ skewness decay, combined with known variance decay, suggests the rate of convergence to Gaussianity is driven by the same dimension parameter, which may persist for higher cumulants."],"forward_implications":["Skewness of Bures-Hall von Neumann entropy is exactly computable for any subsystem dimensions $m \\le n$, enabling precise Edgeworth approximation of the entropy distribution.","In the thermodynamic limit with $m/n = c$, the standardized third cumulant vanishes as $\\Theta(1/n)$, reinforcing the Gaussianity conjecture for Bures-Hall entanglement entropy.","The closed form provides a quantitative benchmark against which numerical sampling of random Bures-Hall states can be tested.","The twelve summation identities extend the library of polygamma finite-sum relations available for cumulant calculations in random matrix ensembles."],"supporting_citations":[{"why":"Supplies the exact variance formula for von Neumann entropy over the Bures-Hall measure and the simplification framework that this paper extends to the third cumulant.","marker":"[34]"},{"why":"Supplies the exact mean formula and the moment-conversion machinery used in the cumulant derivation.","marker":"[35]"},{"why":"Establishes the Pfaffian correlation structure and the link between the Bures measure and the Cauchy two-matrix model.","marker":"[8]"},{"why":"Provides the three-point density and kernel representations used to set up the cumulant integrals.","marker":"[17]"},{"why":"Provides the Cauchy-Laguerre kernel representations and the lemma used in evaluating quadruple sums.","marker":"[2]"},{"why":"Supplies the re-summation technique used to derive the new summation identities.","marker":"[13]"},{"why":"Gives the Hilbert-Schmidt analog whose Edgeworth-correction approach is adapted here.","marker":"[36]"}],"fun_headline_variants":["Exact skewness formula for Bures-Hall states","Bures-Hall entropy skewness closed form","Third cumulant of entropy: exact Bures-Hall","Sharpening entropy distribution via exact skewness","Bures-Hall entropy third cumulant exact"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The derivation stands on the correctness of all twelve new summation identities in Appendix A; if any one of them is wrong, the cancellation of anomalies fails and the closed-form third cumulant is not established.","fun_headline_variants_meta":{"raw":{"variants":["Exact skewness formula for Bures-Hall states","Bures-Hall entropy skewness closed form","Third cumulant of entropy: exact Bures-Hall","Sharpening entropy distribution via exact skewness","Bures-Hall entropy third cumulant exact"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000166,"raw_usage":{"total_tokens":1220,"prompt_tokens":880,"completion_tokens":340,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":496,"completion_tokens_details":{"reasoning_tokens":270}},"tokens_in":496,"tokens_out":340,"duration_ms":3739,"temperature":1.0,"reasoning_tokens":270,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T05:51:52.614644+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Evaluate the twelve identities (A.13)--(A.27) by direct numerical summation for several small integer values of $m$ and parameters $a, b, c$; any mismatch disproves the cancellation. Independently, simulate Bures-Hall random states for a few pairs $(m,n)$ at high precision and compare the sample third cumulant against Proposition 1.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the exact variance formula for von Neumann entropy over the Bures-Hall measure and the simplification framework that this paper extends to the third cumulant."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the exact mean formula and the moment-conversion machinery used in the cumulant derivation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the Pfaffian correlation structure and the link between the Bures measure and the Cauchy two-matrix model."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the three-point density and kernel representations used to set up the cumulant integrals."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the Cauchy-Laguerre kernel representations and the lemma used in evaluating quadruple sums."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the re-summation technique used to derive the new summation identities."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the Hilbert-Schmidt analog whose Edgeworth-correction approach is adapted here."}],"review_version":1}