REVIEW 3 major objections 5 minor 48 references
Explicit Formulas for Estimating Trace of Reduced Density Matrix Powers via Single-Circuit Measurement Probabilities
T0 review · 3 major / 5 minor · reviewed 2026-08-06 · deepseek-v4-flash
Pith's one-line read A single controlled-SWAP circuit can estimate every power trace of a reduced density matrix from one set of measurement probabilities.
desk verdict Useful single-circuit power-trace formulas with a real proof gap in the general case and misstated Tsallis/q-concurrence applications; worth referee time. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the controlled-SWAP test on $n$ copies with $n-1$ control qubits, whose effect is encoded by Kraus operators $K_z=2^{-(n-1)}\prod_{k=1}^{n-1}[I+(-1)^{z_k}S_k]$, where $S_k$ swaps copies $k$ and $k+1$. The argument runs on two mechanisms: trace lemmas that evaluate products of SWAP operators against identical copies, and a parity filter that, when probabilities are summed over outcomes with $c^z_{1,\ldots,k-1}$ even, selects the cyclic term $\operatorname{tr}(\rho^k)$ while cancelling the remaining 'waste' terms $S_{\mathrm{waste}}$ collected in the expansion. This cancellation is what lets one set of $2^{n-1}$ outcome probabilities yield $n-1$ separate trace estimates.
What would settle it
Pick a specific $n\ge 5$ and a state with a known reduced spectrum, such as a random mixed state on two qubits; expand $K_z^\dagger K_z$ exactly through all $2^{2n-2}$ terms, compute $2\sum_{c^z_{1,\ldots,k-1}\ \mathrm{even}} p(z)-1$ for every $k$, and compare with $\operatorname{tr}(\rho_A^k)$ obtained by exact diagonalization. Any mismatch for some $k$ or some state disproves Theorem 3.
Extended reading notes
Core claim
Theorem 3 states that for any $N$-partite state $\rho$ and any $k=2,\dots,n$, the power trace of the reduced state on the first subsystem obeys $\operatorname{tr}(\rho_{A_1}^k)=2\sum_{c^z_{1,2,\ldots,k-1}\ \mathrm{even}} p(z)-1$, where $p(z)$ is the probability of reading the bitstring $z$ on the $n-1$ control qubits and the sum runs over outcomes in which the first $k-1$ control bits have an even number of $1$s. The paper derives this from the Kraus-operator expansion of the circuit and uses it to claim simultaneous estimation of all these traces from a single set of measurement statistics. A technical point of the derivation is that neighbouring SWAP operators do not commute, so $K_z^\dagger K_z\neq K_z$; the paper shows that earlier work missed the resulting cross terms and supplies lemmas that account for them.
Load-bearing premise
The proof of Theorem 3 assumes that all leftover 'waste' terms in the Kraus expansion cancel when the measurement probabilities are summed over outcomes whose first $k-1$ control bits contain an even number of $1$s, for every $n$ and every state; the paper verifies this cancellation only for $n=3$ and $n=4$ and for larger $n$ asserts it 'by a similar discussion'.
Editorial extensions
If this is right
- One fixed circuit with $n$ copies replaces up to $n-1$ purpose-built circuits for estimating $\operatorname{tr}(\rho_A^2),\dots,\operatorname{tr}(\rho_A^n)$.
- The same measurement data feed directly into entanglement quantifiers such as concurrence, q-concurrence, Tsallis-q entanglement, and informationally complete entanglement measures.
- Nonlinear functionals like $\operatorname{tr}(e^{\beta\rho})$ and von Neumann entropy can be approximated by truncating expansions in the estimated power traces.
- The hybrid algorithm needs only $r$ copies when the rank $r$ of $\rho_A$ is known, and the sample bound $O(\epsilon^{-2}\log(n/\delta))$ holds simultaneously for all $k$.
Reading between the lines
- If the waste-term cancellation holds generally, the same probability data could also estimate products of traces such as $\operatorname{tr}(\rho^2)\operatorname{tr}(\rho^3)$, which appear as distinct parity classes in the expansion; the paper does not use this.
- A direct numerical test on random mixed states with $n\ge 5$ copies would be a stronger check of Theorem 3 than the symmetric GHZ and W simulations reported here.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes estimating tr(rho_A^k) for k=2..n using a single controlled-SWAP circuit with n copies of the state. The main result, Theorem 3, claims that tr(rho_A^k) = 2 * sum_{z: c^z_{1..k-1} even} p(z) - 1, and Corollary 1 says one set of measurement data suffices for all powers. The paper derives the n=3 and n=4 cases explicitly, develops a purely quantum algorithm and a hybrid Newton-Girard iteration algorithm, gives a Hoeffding-based sample complexity O(1/epsilon^2 log(n/delta)), applies the formulas to nonlinear functionals and entanglement measures, and reports Qiskit simulations for GHZ and W states.
Significance. If Theorem 3 holds, the single-circuit simultaneous estimation of all power traces is a practical advance over multi-circuit approaches, and the identification that [25]'s probability representation neglects the non-commutativity of neighboring SWAP operators is a useful correction. The n=3,4 derivations are explicit and internally consistent, the formulas are derived from first principles with no fitted parameters, and the Qiskit simulations match the theoretical probabilities for GHZ and W states. However, the general proof of Theorem 3 rests on an unproved cancellation of the remainder S_waste, so the central claim is currently conditional on a nontrivial combinatorial identity that should be established rigorously.
major comments (3)
- [Appendix B, Theorem 3 proof (Eqs. S4-S6)] The proof of Theorem 3 is incomplete at its load-bearing step. After expanding p(z) into cyclic terms and a remainder S_waste, the text asserts that the contribution of S_waste vanishes after summing over z with c^z_{1,...,k-1} even, because each remaining term has a parity structure 'other than' the cyclic ones. This cancellation is verified only for k=2 and k=n; for all intermediate k the proof says 'by a similar discussion,' and the full expansion of K_z^dagger K_z, which has 2^{2n-2} terms, is explicitly not given. Since S_waste contains cross terms such as (tr rho_A^2)^2, the explicit formula tr(rho_A^k) = 2 sum p(z) - 1 depends entirely on this unproved identity. Please supply a general combinatorial or representation-theoretic proof of the cancellation, or give a complete expansion for general n.
- [Section VI (Numerical simulations)] The numerical validation uses only the GHZ and W states, whose reduced density matrices are highly symmetric rank-2 states. This tests a very special slice of the parameter space and cannot substitute for a proof of the S_waste cancellation for generic rho_A or for intermediate k. Adding simulations on random pure states with generic spectra, and on random mixed states if the mixed-state claim is retained, would materially strengthen the evidence and could reveal whether the formula holds beyond the symmetric cases.
- [Section V.B, Definition 4 and following display] The probabilistic representation of q-concurrence is inconsistent with the definition. From Definition 4, C_q = 1 - tr(rho_A^q), and Theorem 3 gives tr(rho_A^q) = 2 sum_{even} p(z) - 1, so C_q = 2 - 2 sum_{even} p(z). The paper instead writes C_q = 2 sum p(z), which would imply tr(rho_A^q) = 1 - 2 sum p(z), contradicting Theorem 3. Please correct the displayed formula.
minor comments (5)
- [Remark after Theorem 3] The remark extending the result to mixed states via convex decomposition and linearity of the trace is not correct as stated: p(z) = tr(K_z^dagger K_z rho^{otimes n}) is a degree-n function of rho, not an affine function of rho. The Kraus-operator proof itself already applies to arbitrary rho through the spectral decomposition used in the lemmas, so please delete or rephrase this remark.
- [Section V.A, von Neumann entropy] The Taylor expansion of tr(rho ln rho) around the identity requires ||rho - I|| < 1, which holds only for full-rank rho. Please state this condition explicitly or use a series valid on the relevant part of the spectrum.
- [Section IV, Table I and Section VI, Table II] The table references are inconsistent: Section IV introduces a comparison table as Table I, while the error table in Section VI is also called Table I in the text but labeled Table II in the caption. The sentence introducing the comparison table is also incomplete. Please renumber and clean up the references.
- [Appendix B, Eq. (S4)] There are several typos in the proof of Theorem 3, including the definition of sigma_n (the text says permutations of 1,...,n-2 and 1,...,n-1 while sigma_n appears in the notation), and garbled expressions in the proof of Lemma 4 such as '|iim+1>' and 'Sis'. These should be corrected in a careful copyedit.
- [Section III, hybrid algorithm] The improved Newton-Girard algorithm is stated to require r quantum copies and assumes prior knowledge of r = rank(rho_A). This assumption should be stated explicitly in the algorithm description in Section III, not only in the conclusions.
Circularity Check
No circular derivation found; the general-n proof has an unproved S_waste cancellation step, which is a correctness gap rather than a circularity.
full rationale
The central formula tr(ρ^k_{A1}) = 2 Σ_{c^z_{1,...,k-1} even} p(z) − 1 is derived from the Kraus expansion of the controlled-SWAP channel and from trace identities for products of SWAP operators (Lemmas 1-4). The probabilities p(z) are not fitted to the target traces, and the numerical validation uses analytically known GHZ and W values as external benchmarks, so the result is not equivalent to its inputs by construction. The author-overlapping citations, such as [25] and [27], are used for context, for circuit inspiration, or as the target of correction; they do not carry the proof of Theorem 3. The genuine weakness is the S_waste cancellation asserted in the proof of Theorem 3. In the main text and again in Appendix B, the paper states that expanding K_z†K_z fully is "impractical and unnecessarily complex" and denotes the remaining terms by S_waste, then asserts that these terms cancel over the even-parity sums without exhibiting the general expansion or proving the cancellation for general n; only n = 3 and n = 4 are worked out, and intermediate k is dismissed with "By a similar discussion." This is an omitted proof and a correctness risk, not a circular reduction: S_waste is not defined as the cancellation conclusion, and the asserted cancellation is a separate combinatorial-mathematical claim about the expansion rather than a restatement of the theorem's input. Accordingly, the circularity score is low.
Assumptions & free parameters
assumptions (3)
- standard math Any multipartite density matrix admits an expansion as a linear combination of pure separable states (Lemma 5).
- domain assumption The controlled-SWAP channel on n copies has Kraus operators K_z = 2^{-(n-1)} ∏_{k=1}^{n-1} [I + (-1)^{z_k} S_k], so p(z) = tr(K_z^† K_z ρ^{⊗n}).
- domain assumption The Taylor expansion tr(ρ ln ρ) = tr[ρ Σ_{n=1}∞ (-1)^{n+1}/n (ρ-I)^n] requires ||ρ-I|| < 1 for some matrix norm.
Cite this review
Pith. "Pith review of Explicit Formulas for Estimating Trace of Reduced Density Matrix Powers via Single-Circuit Measurement Probabilities." pith.science (2026). https://pith.science/paper/ML5UEFA7
@misc{pith2026250717117,
author = {Pith},
title = {Pith review of: Explicit Formulas for Estimating Trace of Reduced Density Matrix Powers via Single-Circuit Measurement Probabilities},
year = {2026},
howpublished = {\url{https://pith.science/paper/ML5UEFA7}},
note = {Machine review of arXiv:2507.17117}
}
abstract
In the fields of quantum mechanics and quantum information science, the traces of reduced density matrix powers play a crucial role in the study of quantum systems and have numerous important applications. In this paper, we propose a universal framework to simultaneously estimate the traces of the $2$nd to the $n$th powers of a reduced density matrix using a single quantum circuit with $n$ copies of the quantum state. Specifically, our approach leverages the controlled SWAP test and establishes explicit formulas connecting measurement probabilities to these traces. We further develop two algorithms: a purely quantum method and a hybrid quantum-classical approach combining Newton-Girard iteration. Rigorous analysis via Hoeffding inequality demonstrates the method's efficiency, requiring only $M=O\left(\frac{1}{\epsilon^2}\log(\frac{n}{\delta})\right)$ measurements to achieve precision $\epsilon$ with confidence $1-\delta$. Additionally, we explore various applications including the estimation of nonlinear functions and the representation of entanglement measures. Numerical simulations are conducted for two maximally entangled states, the GHZ state and the W state, to validate the proposed method.
Figures
Figures from the paper (14 more)
Reference graph
Works this paper leans on
-
[25]
Z. X. Jin, X. Li-Jost, S. M. Fei, C. F. Qiao, Entanglement measures based on the complete information of reduced states. Phys. Rev. A 107(1), 012409 (2023)
work page 2023
-
[1]
R. Horodecki, P. Horodecki, M. Horodecki, K. Horodecki, Quantum entanglement. Reviews of modern physics 81(2), 865-942 (2009)
work page 2009
-
[2]
V. Vedral, Quantum entanglemen. Nature Physics10(4), 256-258 (2014)
work page 2014
-
[3]
J. Barrett, Nonsequential positive-operator-valued mea- surements on entangled mixed states do not always vio- late a Bell inequality. Phys. Rev. A65(4), 042302 (2002)
work page 2002
- [4]
-
[5]
Bose, Quantum communication through an unmodu- lated spin chain
S. Bose, Quantum communication through an unmodu- lated spin chain. Phys. Rev. Lett.91(20), 207901 (2003)
work page 2003
- [6]
-
[7]
D. J. Bernstein, T. Lange, Post-quantum cryptography. Nature 549(7671), 188-194 (2017)
work page 2017
Show all 48 references
-
[8]
D. P. DiVincenzo, Quantum computation. Science 270(5234), 255-261 (1995)
1995
-
[9]
H. J. Briegel, D. E. Browne, W. D¨ ur, R. Raussendorf, M. Van den Nest, Measurement-based quantum compu- tation. Nature Physics 5(1), 19-26 (2009)
2009
-
[10]
Horodecki, Entanglement measures
M. Horodecki, Entanglement measures. Quantum Inf. Comput. 1(1), 3-26 (2001)
2001
-
[11]
G¨ uhne, M
O. G¨ uhne, M. Reimpell, R. F. Werner, Estimating en- tanglement measures in experiments. Phys. Rev. Lett. 98(11), 110502 (2007)
2007
-
[12]
Szalay, Multipartite entanglement measures
S. Szalay, Multipartite entanglement measures. Phys. Rev. A 92(4), 042329 (2015)
2015
-
[13]
W. K. Wootters, Entanglement of formation and concur- rence. Quantum Inf. Comput. 1(1), 27-44 (2001)
2001
-
[14]
H. Fan, K. Matsumoto, H. Imai, Quantify entanglement by concurrence hierarchy. Journal of Physics A: Mathe- matical and General 36(14), 4151 (2003)
2003
-
[15]
W. K. Wootters, Entanglement of formation of an arbi- trary state of two qubits. Phys. Rev. Lett. 80(10), 2245 (1998)
1998
-
[16]
K. Chen, S. Albeverio, S. M. Fei, Entanglement of for- mation of bipartite quantum states. Phys. Rev. Lett. 95(21), 210501 (2005)
2005
-
[17]
Henderson, V
L. Henderson, V. Vedral, Information, relative entropy of entanglement, and irreversibility. Phys. Rev. Lett. 84(10), 2263 (2000)
2000
-
[18]
Moitra, R
U. Moitra, R. M. Soni, S. P. Trivedi, Entanglement en- tropy, relative entropy and duality. Journal of High En- ergy Physics 2019(8), 1-24 (2019)
2019
-
[19]
A. R. Cullen, P. Kok, Calculating concentratable entan- glement in graph states. Phys. Rev. A 106(4), 042411 (2022)
2022
-
[20]
Schatzki, G
L. Schatzki, G. Liu, M. Cerezo, E. Chitambar, Hierar- chy of multipartite correlations based on concentratable entanglement. Physical Review Research 6(2), 023019 (2024)
2024
-
[21]
J. L. Beckey, G. Pelegr ´ ı, S. Foulds, N. J. Pearson, Mul- tipartite entanglement measures via Bell-basis measure- ments. Phys. Rev. A 107(6), 062425 (2023)
2023
-
[22]
J. L. Beckey, N. Gigena, P. J. Coles, M. Cerezo, Com- putable and operationally meaningful multipartite en- tanglement measures. Phys. Rev. Lett. 127(14), 140501 (2021)
2021
-
[23]
H. Li, T. Gao, F. Yan, Parametrized multipartite entan- glement measures. Phys. Rev. A 109(1), 012213 (2024)
2024
-
[24]
Zhou, Z.X
W. Zhou, Z.X. Shen, D.P. Xuan, Z.X. Wang, S.M. Fei, Parameterized bipartite entanglement measures and en- tanglement constraints. Adv. Quantum Technol. 2400707 (2025)
2025
-
[26]
Foulds, V
S. Foulds, V. Kendon, T. Spiller, The controlled SW AP test for determining quantum entanglement. Quantum Science and Technology 6(3), 035002 (2021)
2021
-
[27]
R. Q. Zhang, Y.D. Qu, S.Q. Shen, M. Li, J. Wang, The controlled SW AP test for entanglement of mixed quan- tum states. Europhysics Letters 146(1), 18001 (2024)
2024
-
[28]
A. K. Ekert, C. M. Alves, D. K. Oi, M. Horodecki, P. Horodecki, L. C. Kwek, Direct estimations of linear and nonlinear functionals of a quantum state. Phys. Rev. Lett. 88(21), 217901 (2002). 17
2002
-
[29]
Johri, D
S. Johri, D. S. Steiger, M. Troyer, Entanglement spec- troscopy on a quantum computer. Phys. Rev. B 96(19), 195136 (2017)
2017
-
[30]
Suba¸ sı, L
Y. Suba¸ sı, L. Cincio, P. J. Coles, Entanglement spec- troscopy with a depth-two quantum circuit. J. Phys. A 52(4), 044001 (2019)
2019
-
[31]
Yirka, Y
J. Yirka, Y. Suba¸ sı, Qubit-efficient entanglement spec- troscopy using qubit resets. Quantum 5, 535 (2021)
2021
-
[32]
M. Shin, J. Lee, S. Lee, K. Jeong, Rank Is All You Need: Estimating the Trace of Powers of Density Ma- trices. arXiv:2408.00314 (2024)
2024 arXiv
-
[33]
H. Y. Huang, R. Kueng, J. Preskill, Predicting many properties of a quantum system from very few measure- ments. Nature Physics 16(10), 1050-1057 (2020)
2020
-
[34]
R. Yang, W. Zhao, T. Chen, Scalable algorithms for cal- culating power functions of random quantum states in the noisy intermediate-scale quantum era. Phys. Rev. A 109(6), 062421 (2024)
2024
-
[35]
Z. Liu, Y. Tang, H. Dai, P. Liu, S. Chen, X. Ma, De- tecting entanglement in quantum many-body systems via permutation moments. Phys. Rev. Lett.129(26), 260501 (2022)
2022
-
[36]
Neven, J
A. Neven, J. Carrasco, V. Vitale, C. Kokail, A. Elben, M. Dalmonte, P. Calabrese, P. Zoller, B. Vermersch, R. Kueng, B. Kraus, Symmetry-resolved entanglement de- tection using partial transpose moments. npj Quantum Information, 7(1), 152 (2021)
2021
-
[37]
Natalini, P
P. Natalini, P. E. Ricci, Newton sum rules of polynomials defined by a three-term recurrence relation. Computers & Mathematics with Applications 42(6-7), 767-771 (2001)
2001
-
[38]
Hoeffding, Probability inequalities for sums of bounded random variables
W. Hoeffding, Probability inequalities for sums of bounded random variables. J. Am. Stat. Assoc. 58, 13(1963)
1963
-
[39]
Rotter, Dynamics of quantum systems
I. Rotter, Dynamics of quantum systems. Phys. Rev. E 64(3), 036213 (2001)
2001
-
[40]
Lloyd, M
S. Lloyd, M. Mohseni, P. Rebentrost, Quantum princi- pal component analysis. Nature physics 10(9), 631-633 (2014)
2014
-
[41]
P. Boes, J. Eisert, R. Gallego, M. P. M¨ uller, H. Wilming, Von Neumann entropy from unitarity. Phys. Rev. Lett. 122(21), 210402 (2019)
2019
-
[42]
Y. Wang, G. Li, X. Wang, Variational quantum Gibbs state preparation with a truncated Taylor series. Phys. Rev. Appl. 16(5), 054035 (2021)
2021
-
[43]
Consiglio, J
M. Consiglio, J. Settino, A. Giordano, C. Mastroianni, F. Plastina, S. Lorenzo, S. Maniscalco, J. Goold, T. JG Apollaro, Variational Gibbs state preparation on noisy intermediate-scale quantum devices. Phys. Rev. A 110(1), 012445 (2024)
2024
-
[44]
San Kim, Generalized entanglement constraints in multi-qubit systems in terms of Tsallis entropy
J. San Kim, Generalized entanglement constraints in multi-qubit systems in terms of Tsallis entropy. Annals of Physics 373, 197-206 (2016)
2016
-
[45]
Y. Luo, T. Tian, L. H. Shao, Y. Li, General monogamy of Tsallis q-entropy entanglement in multiqubit systems. Phys. Rev. A 93(6), 062340 (2016)
2016
-
[46]
X. Yang, M. X. Luo, Y. H. Yang, S. M. Fei, Parametrized entanglement monotone. Phys. Rev. A 103(5), 052423 (2021)
2021
-
[47]
J. I. De Vicente, Separability criteria based on the bloch representation of density matrices. Quantum Inf Comput 7(7), 624-638 (2007)
2007
-
[48]
cz 1,··· ,s is even
M. S. Byrd, N. Khaneja, Characterization of the positiv- ity of the density matrix in terms of the coherence vector representation. Phys. Rev. A 68(6), 062322 (2003). APPENDIX In this section, we will present the proof of the lemmas and theorems in the main text. For the conve...
2003
Reviewed August 6, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.