REVIEW 3 major objections 5 minor 38 references
Compressed sensing quantum state tomography for qudits: A comparison of Gell-Mann and Heisenberg-Weyl observable bases
T0 review · 3 major / 5 minor · reviewed 2026-08-15 · deepseek-v4-flash
Pith's one-line read Compressed sensing quantum state tomography for qudits is more measurement-efficient with the Heisenberg-Weyl observable basis than with the generalized Gell-Mann basis, and the advantage grows with the qudit dimension.
desk verdict Solid, clearly described numerical comparison of GGM vs HWO bases for two-qudit CS-QST on rank-1 states, but the abstract overgeneralizes beyond what was simulated. 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 minimum coherence $\nu_{\min}\equiv d\max_a\|B_a\|_2$ defined in Eq. (17) from compressed sensing theory, which measures how concentrated tomographic weight is across the basis elements. For the GGM basis, the diagonal operators have spectral norms that grow with $k$, giving $\nu_{\min}=(k-1)^N$, while for the HWO basis all non-identity elements share the same spectral norm, giving $\nu_{\min}\approx 2^N$ (exactly $2^N$ when $k$ is a multiple of 8, and exactly $1$ for $k=4$). The paper uses this coherence difference to predict and explain why HWO requires fewer random measurements, alongside the singular value thresholding (SVT) algorithm used for reconstruction.
What would settle it
Repeat the same protocol for two-qudit states of ranks $r=2$ and $r=3$ at $k=9,12,15$ and for $N=3$ or $4$ qudits, keeping the same noise model, and compute $m_{95\%}^{\mathrm{GGM}}-m_{95\%}^{\mathrm{HWO}}$; if the gap stops growing as $k^2$, or reverses, the central claim holds only within the simulated regime.
Extended reading notes
Core claim
The central claim is that in CS-QST of two-qudit systems, the HWO basis becomes more efficient than the GGM basis as the qudit dimension $k$ increases: the minimum number $m_{95\%}$ of randomly selected measurement operators needed to keep mean fidelity minus one standard deviation at or above 95% is smaller for HWO, and the gap between the two bases grows approximately as $k^2$ (Eqs. 28-29, Fig. 3). The paper explains this through the minimum coherence $\nu_{\min}=d\max_a\|B_a\|_2$: for GGM it grows as $(k-1)^N$, whereas for HWO it is nearly constant in $k$, taking values close to $2^N$ and equalling $2^N$ exactly when $k$ is a multiple of 8. The same coherence argument accounts for the paper's fixed-dimension comparison at $d=16$, where HWO-based qudit tomography matches Pauli-based qubit tomography while GGM degrades with increasing $k$.
Load-bearing premise
The central claim rests on simulations limited to rank-1 pure states of two qudits with a fixed 5% depolarizing channel and Gaussian measurement noise of standard deviation $0.1/d$; the general conclusion that HWO becomes more efficient as qudit dimension increases assumes this ordering survives for mixed low-rank states, more qudits, and other noise levels.
Editorial extensions
If this is right
- For two-qudit systems, the number of measurements needed to reach 95% mean fidelity minus one standard deviation grows more slowly with $k$ for the HWO basis than for the GGM basis, and the gap grows roughly as $k^2$.
- At a fixed Hilbert space dimension $d=16$, qudit tomography using the HWO basis matches the efficiency of four-qubit tomography using the Pauli basis, while GGM efficiency degrades as $k$ increases.
- The coherence analysis predicts that GGM becomes less reliable at small $m/d^2$ because its diagonal operators have larger spectral norms, and the simulations confirm larger error bars for GGM in that regime.
- Both bases eventually achieve high fidelity, so an experiment that can more easily implement GGM measurements can still use GGM at low or moderate $k$; high-dimensional implementations should prefer HWO.
Reading between the lines
- A direct testable extension is to run the same comparison for mixed low-rank states ($r=2,3$); if the coherence argument is the true driver, the HWO advantage should persist because coherence is a property of the basis rather than of the state rank.
- The coherence analysis suggests a mitigation for GGM: sampling measurement operators with probabilities weighted by their spectral norms, or always including the diagonal operators, could reduce the gap; this is worth testing.
- The same basis-efficiency logic likely applies to other compressed-sensing recovery tasks in quantum information, such as channel or process tomography, where the measurement basis is also chosen rather than fixed.
- Because the paper fixes the SVT algorithm and its hyperparameters, the size of the gap may shift with a different reconstruction algorithm; repeating the comparison with another matrix-completion method would separate basis effects from algorithm effects.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper numerically compares compressed sensing quantum state tomography (CS-QST) for qudit systems using two operator bases: the generalized Gell-Mann (GGM) basis and the Heisenberg-Weyl observable (HWO) basis. The authors generate rank-1 Ginibre states, apply 5% depolarizing noise and Gaussian measurement noise, reconstruct via singular value thresholding, and evaluate fidelity and trace distance. For two-qudit systems (N=2) with qudit dimension k=3,6,9,12,15, they report that the HWO basis requires fewer randomly selected measurements to achieve a mean fidelity minus one standard deviation of at least 0.95, and they fit the measured m_95% to functions of k (Eqs. 28-29), claiming the GGM-HWO gap grows approximately as k^2. They also compare fixed-dimension d=16 systems and attribute the observed differences to the coherence properties of the two bases, defined via a 'minimum coherence' ν_min (Eq. 17). The abstract states without qualification that 'the HWO basis becomes more efficient as the qudit dimension increases.'
Significance. If the reported trend is robust, the paper offers practically useful guidance for choosing measurement bases in qudit CS-QST, an area relevant to high-dimensional quantum information processing. The numerical protocol is clearly described (Ginibre states, SVT, AIC model selection), and the paper provides explicit fitting formulas that can be tested against other simulations or experiments. The main theoretical support, however, rests on a coherence calculation that, as written, is inconsistent with the paper's own definition of ν_min, and the empirical claims are based on a restricted set of simulations (rank-1, N=2) without error bars on the m_95% values. These issues do not, by themselves, invalidate the numerical comparison, but they make the central claim stronger than the evidence presented.
major comments (3)
- [Eq. (17), Eq. (31), Table I, Fig. 1] The coherence measure ν_min defined in Eq. (17) as d max_a ||B_a||_2 is inconsistent with the values reported in Eq. (31) and Table I for the HWO basis. With B_a = (ω_{i1}⊗...⊗ω_{iN})/√(k^N) as in Eq. (2), which the text calls an orthonormal basis, the spectral norm of B_a is (max_i ||ω_i||_2)^N / k^{N/2}, so ν_min = k^{N/2} (max_i ||ω_i||_2)^N, not (max_i ||ω_i||_2)^N as stated in Sec. II C. For the HWO basis, using the eigenvalue expression in Eq. (30), this gives ν_min = k^{N/2} (max_n ε_n^2)^{N/2}, not (max_n ε_n^2)^N as in Eq. (31). For example, for k=3, max_n ε_n^2 = 1.866 but ||W(0,1)||_2 = 1.366, so the value 1.866 reported in Table I is a squared eigenvalue, not a spectral-norm coherence. Under the corrected definition, HWO ν_min grows as k^{N/2}(√2)^N for k a multiple of 8, rather than remaining constant. This invalidates the statement in Sec. II C that ν_min for the HWO basis is 'nearly constant regardless of k' and undercuts the coherence-based explanation in Sec. III C. The authors should correct the normalization in Eq. (31) and revise Fig. 1, Table I, and the related discussion accordingly.
- [Abstract and Secs. III B, III C 1] The abstract and Sec. IV state without qualification that 'the HWO basis becomes more efficient as the qudit dimension increases,' but the numerical evidence is obtained only for N=2 two-qudit systems with rank-1 pure states (Sec. III B: 'we consider only pure states with rank r=1'), a fixed 5% depolarizing channel, and Gaussian measurement noise with standard deviation 0.1/d. No simulations for mixed states (r>1) or for N>2 are reported, so the unqualified claim is an extrapolation from a narrow regime. The abstract and summary should either explicitly restrict the claim to the simulated conditions (e.g., 'for two-qudit pure states under the noise model considered') or the authors should add simulations for other ranks and numbers of qudits to support the general statement.
- [Sec. III C 1, Eqs. (28)-(29), Fig. 3] The quantitative claim that the GGM-HWO gap in m_95% grows roughly as k^2 is based on least-squares fits to m_95% values, but the m_95% data points are reported without error bars. Each m_95% is determined from 50 repetitions (Sec. III B), so the sample mean and standard deviation used in the condition 'mean fidelity minus one standard deviation ≥ 0.95' carry statistical uncertainty that propagates to m_95%. Since no uncertainties are reported for the fitted parameters in Eqs. (28)-(29), the fitted k^2 scaling is not statistically grounded. The authors should provide error estimates for m_95% (e.g., via bootstrapping) and report confidence intervals for the fit parameters, or otherwise demonstrate that the observed trend is robust to the sampling uncertainty.
minor comments (5)
- [Sec. III B] The text describes the target state as a rank-1 pure state but then applies depolarizing noise to obtain \tilde{ρ}, which has full rank. Since the fidelity in Eq. (25) is computed against the original ρ, the effective reconstruction target is a pure state corrupted by noise, not the depolarized state itself; this is a reasonable protocol, but it would be clearer to state explicitly that ρ, not \tilde{ρ}, is the target for the fidelity and trace distance.
- [Sec. II C] The sentence 'Section II C presents our numerical simulations' is a typo; it should refer to Section III.
- [Acknowledgments] The sentence 'The work of was supported' is incomplete; the author names are missing.
- [Sec. III C 1] The connection between the fitted m_95% curves and the theoretical scaling m ∼ O(rd log^2 d) is only qualitative; the fitted forms have p_1 >> p_0, so the leading behavior is m_95% ∼ p_1 k^2, and the paper would benefit from a brief explanation of why the logarithmic factor is subleading in the simulated range.
- [Appendix B] Equation (31) uses the notation (max_n ε_n^2)^N without explicitly stating that ε_n are the eigenvalues of the single-qudit HWO operator; relating this to the corrected definition of ν_min discussed above would improve clarity.
Circularity Check
No significant circularity: the HWO-vs-GGM comparison is produced by symmetric numerical simulation, and the m95% scaling curves are presented as descriptive fits, not predictions.
full rationale
The central derivation is self-contained. In Sec. III B-C, GGM and HWO reconstructions are generated from the same Ginibre rank-1 states with the same 5% depolarizing noise and Gaussian measurement noise, and processed by the same SVT algorithm, so the observed HWO advantage is an emergent simulation result, not an input. Eqs. (28)-(29) are explicitly least-squares/AIC fits to the simulated m95% data (Sec. III C 1 and Appendix V D), and the paper describes them as fitting curves rather than as predictions, so no fitted parameter is renamed as a prediction. The coherence mechanism in Sec. II C is imported from external references (Gross [27]; Asadian et al. [24]) and is used as an explanatory hypothesis tested by the independent numerics; there is no self-citation chain and no uniqueness theorem from the authors' prior work forcing the choice of basis. The paper itself states the scope limitation to rank-1, N=2 states (Sec. III B: 'we consider only pure states with rank r=1'), which narrows generality but is not a circular dependency. The apparent normalization mismatch between the nu_min definition in Eq. (17) and its evaluation in Eq. (31) is a consistency/correctness concern in the theoretical explanation, but it does not make the numerical comparison circular, since the simulations do not use nu_min as an input.
Assumptions & free parameters
free parameters (8)
- GGM p0 =
-0.117
- GGM p1 =
18.865
- GGM p2 =
-53.044
- GGM p3 =
42.449
- HWO p0 =
-0.123
- HWO p1 =
13.897
- HWO p2 =
-31.882
- HWO p3 =
20.894
assumptions (5)
- domain assumption The SVT algorithm (Ref [5]) reconstructs the low-rank state from the random subset P_Omega(M) under the simulations' sampling model.
- domain assumption Rank-1 pure states generated from the Ginibre ensemble, with 5% depolarizing noise and Gaussian measurement noise N(0,0.1/d), are representative targets for low-rank CS-QST.
- domain assumption The minimum coherence nu_min (Eq. 17) is a meaningful basis-dependent predictor of practical reconstruction efficiency.
- standard math The normalized basis B_a=(omega_{i1} tensor ... tensor omega_{iN})/sqrt(d) with Tr[omega_i omega_j]=k delta_{ij} is a valid orthonormal Hermitian basis for expanding density matrices (Eq. 2).
- domain assumption The AIC-selected functional form in Eqs. (28)-(29) is an adequate model for the m95% data over the simulated range.
Cite this review
Pith. "Pith review of Compressed sensing quantum state tomography for qudits: A comparison of Gell-Mann and Heisenberg-Weyl observable bases." pith.science (2026). https://pith.science/paper/6JKVEUHM
@misc{pith2026250510462,
author = {Pith},
title = {Pith review of: Compressed sensing quantum state tomography for qudits: A comparison of Gell-Mann and Heisenberg-Weyl observable bases},
year = {2026},
howpublished = {\url{https://pith.science/paper/6JKVEUHM}},
note = {Machine review of arXiv:2505.10462}
}
read the original abstract
Quantum state tomography (QST) is an essential technique for reconstructing the density matrix of an unknown quantum state from measurement data, crucial for quantum information processing. However, conventional QST requires an exponentially growing number of measurements as the system dimension increases, posing a significant challenge for high-dimensional systems. To mitigate this issue, compressed sensing quantum state tomography (CS-QST) has been proposed, significantly reducing the required number of measurements. In this study, we investigate the impact of basis selection in CS-QST for qudit systems, which are fundamental to high-dimensional quantum information processing. Specifically, we compare the efficiency of the generalized Gell-Mann (GGM) and Heisenberg-Weyl observable (HWO) bases by numerically reconstructing density matrices and evaluating reconstruction accuracy using fidelity and trace distance metrics. Our results demonstrate that, while both bases allow for successful density matrix reconstruction, the HWO basis becomes more efficient as the qudit dimension increases. Furthermore, we find the best fitting curves that estimate the number of measurement operators required to achieve a fidelity of at least 95%. These findings highlight the significance of basis selection in CS-QST and provide valuable insights for optimizing measurement strategies in high-dimensional quantum state tomography.
Figures
Reference graph
Works this paper leans on
-
[1]
Generalized Gell-Mann basis The generalized Gell-Mann (GGM) basis [23], is formed by the standard set of generators of the SU(k) Lie algebra. The elements of the GGM basis are divided into three types as follows: i):Symmetric matrices Λjj ′ s =|j⟩⟨j′|+|j ′⟩⟨j|,(0≤j <j ′≤k−1) (3) ii):Anti-symmetric matrices Λjj ′ a =−i|j⟩⟨j ′|+i|j ′⟩⟨j|,(0≤j <j ′≤k−1) (4) ...
-
[2]
Heisenberg-Weyl (HW) ma- trices
Heisenberg-Weyl observable basis The matrices known as “Heisenberg-Weyl (HW) ma- trices” or “Sylvester’s generalized Pauli matrices” are a non-Hermitian extension of the Pauli matrices to higher dimensions [31]. Let us consider the “shift operator”S and “clock operator”C, which act on the state bases as S|j⟩=|j+ 1 modk⟩andC|j⟩=e i2πj/k|j⟩. Using these ope...
-
[3]
Two-qudit states Figures 2(a-e) show the efficiency of CS-QST for an ensemble of random pure states of two qudits (N= 2) using two different bases, as a function of the numberm of randomly selected measurement operators, for qudit dimensionsk= 3,6,9,12, and 15, respectively. The re- sults provide insights into how the choice of basis affects reconstructio...
-
[4]
Fixed Hilbert space dimension Finally, we examine how the choice of quantum in- formation unit—qubit or qudit—affects the efficiency of CS-QST when the Hilbert space dimensiondis fixed. By varyingNandkwhile keepingd=k N constant, we com- pare the relative advantages of different representations. Specifically, we focus on the case of a density matrix withd...
-
[5]
H¨ affner et al., Nature (london), Nature438, 643 (2005)
H. H¨ affner et al., Nature (london), Nature438, 643 (2005)
work page 2005
- [6]
-
[7]
J. J. Renema, G. Frucci, M. J. A. de Dood, R. Gill, A. Fiore, and M. P. van Exter, Tomography and state re- construction with superconducting single-photon detec- tors, Phys. Rev. A86, 062113 (2012)
work page 2012
-
[8]
Gross, Y.-K
D. Gross, Y.-K. Liu, S. T. Flammia, S. Becker, and J. Eis- ert, Quantum state tomography via compressed sensing, Phys. Rev. Lett.105, 150401 (2010)
2010
Show all 38 references
-
[9]
J.-F. Cai, E. J. Candes, and Z. Shen, A singular value thresholding algorithm for matrix completion (2008), arXiv:0810.3286 [math.OC]
2008 arXiv
-
[10]
Shabani, R
A. Shabani, R. L. Kosut, M. Mohseni, H. Rabitz, M. A. Broome, M. P. Almeida, A. Fedrizzi, and A. G. White, Efficient measurement of quantum dynamics via com- pressive sensing, Phys. Rev. Lett.106, 100401 (2011)
2011
-
[11]
Riofr´ ıo, D
C. Riofr´ ıo, D. Gross, S. Flammia, et al., Experimental quantum compressed sensing for a seven-qubit system, Nature Communications8, 15305 (2017)
2017
-
[12]
Y. Wang, Z. Hu, B. C. Sanders, and S. Kais, Qudits and high-dimensional quantum computing, Frontiers in Physics8, 10.3389/fphy.2020.589504 (2020)
2020
-
[13]
Babazadeh, M
A. Babazadeh, M. Erhard, F. Wang, M. Malik, R. Nouroozi, M. Krenn, and A. Zeilinger, High- dimensional single-photon quantum gates: Concepts and experiments, Phys. Rev. Lett.119, 180510 (2017)
2017
-
[14]
H.-H. Lu, Z. Hu, M. S. Alshaykh, A. J. Moore, Y. Wang, 13 P. Imany, A. M. Weiner, and S. Kais, Quantum phase estimation with time-frequency qudits in a single photon, Advanced Quantum Technologies3, 1900074
-
[15]
P. Liu, R. Wang, J.-N. Zhang, Y. Zhang, X. Cai, H. Xu, Z. Li, J. Han, X. Li, G. Xue, W. Liu, L. You, Y. Jin, and H. Yu, Performing SU(d) operations and rudimentary al- gorithms in a superconducting transmon qudit ford= 3 andd= 4, Phys. Rev. X13, 021028 (2023)
2023
-
[16]
Ringbauer, M
M. Ringbauer, M. Meth, L. Postler, R. Stricker, R. Blatt, P. Schindler, and T. Monz, A universal qudit quantum processor with trapped ions, Nat. Phys.18, 1053 (2022)
2022
-
[17]
P. Hrmo, B. Wilhelm, L. Gerster, M. W. van Mourik, M. Huber, R. Blatt, P. Schindler, T. Monz, and M. Ring- bauer, Native qudit entanglement in a trapped ion quantum processor, NATURE COMMUNICATIONS14, 10.1038/s41467-023-37375-2 (2023)
2023 doi
-
[18]
M. Meth, J. Zhang, J. F. Haase, C. Edmunds, L. Postler, A. J. Jena, A. Steiner, L. Dellantonio, R. Blatt, P. Zoller, T. Monz, P. Schindler, C. Muschik, and M. Ring- bauer, Simulating two-dimensional lattice gauge theo- ries on a qudit quantum computer, NATURE PHYSICS 10.1038/s...
2025 doi
-
[19]
Chicco, G
S. Chicco, G. Allodi, A. Chiesa, E. Garlatti, C. D. Buch, P. Santini, R. De Renzi, S. Piligkos, and S. Carretta, Proof-of-concept quantum simulator based on molecular spin qudits, Journal of the American Chemical Society 146, 1053–1061 (2023)
2023
-
[20]
Lindon, A
J. Lindon, A. Tashchilina, L. W. Cooke, and L. J. LeBlanc, Complete unitary qutrit control in ultracold atoms, Phys. Rev. Appl.19, 034089 (2023)
2023
-
[21]
Gedik, I
Z. Gedik, I. A. Silva, B. Cakmak, G. Karpat, E. L. G. Vidoto, D. O. Soares-Pinto, E. R. deAzevedo, and F. F. Fanchini, Computational speed-up with a single qudit, SCIENTIFIC REPORTS5, 10.1038/srep14671 (2015)
2015 doi
-
[22]
T. Roy, Z. Li, E. Kapit, and D. Schuster, Two-qutrit quantum algorithms on a programmable superconducting processor, Phys. Rev. Appl.19, 064024 (2023)
2023
-
[23]
S. Cao, M. Bakr, G. Campanaro, S. D. Fasciati, J. Wills, D. Lall, B. Shteynas, V. Chidambaram, I. Rungger, and P. Leek, Emulating two qubits with a four-level transmon qudit for variational quantum algorithms, Quantum Sci- ence and Technology9, 035003 (2024)
2024
-
[24]
E. O. Kiktenko, A. S. Nikolaeva, and A. K. Fedorov, Realization of quantum algorithms with qudits (2023), arXiv:2311.12003 [quant-ph]
2023 arXiv
-
[25]
R. T. Thew, K. Nemoto, A. G. White, and W. J. Munro, Qudit quantum-state tomography, Physical Review A 66, 10.1103/physreva.66.012303 (2002)
2002 doi
-
[26]
Badveli, V
R. Badveli, V. Jagadish, R. Srikanth, and F. Petruccione, Compressed-sensing tomography for qudits in hilbert spaces of non-power-of-two dimensions, Phys. Rev. A 101, 062328 (2020)
2020
-
[27]
R. A. Bertlmann and P. Krammer, Bloch vectors for qu- dits, Journal of Physics A: Mathematical and Theoretical 41, 235303 (2008)
2008
-
[28]
Asadian, P
A. Asadian, P. Erker, M. Huber, and C. Kl¨ ockl, Heisenberg-weyl observables: Bloch vectors in phase space, Physical Review A94, 10.1103/phys- reva.94.010301 (2016)
2016 doi
-
[29]
Agnew, J
M. Agnew, J. Leach, M. McLaren, F. S. Roux, and R. W. Boyd, Tomography of the quantum state of photons en- tangled in high dimensions, Phys. Rev. A84, 062101 (2011)
2011
-
[30]
A. M. P˘ alici, T.-A. Isdrail˘ a, S. Ataman, and R. Ion- icioiu, Oam tomography with heisenberg–weyl observ- ables, Quantum Science and Technology5, 045004 (2020)
2020
-
[31]
Gross, Recovering low-rank matrices from few coef- ficients in any basis, IEEE Transactions on Information Theory57, 1548 (2011)
D. Gross, Recovering low-rank matrices from few coef- ficients in any basis, IEEE Transactions on Information Theory57, 1548 (2011)
2011
-
[32]
S. T. Flammia, D. Gross, Y.-K. Liu, and J. Eisert, Quan- tum tomography via compressed sensing: error bounds, sample complexity and efficient estimators, New Journal of Physics14, 095022 (2012)
2012
-
[33]
Kalev, R
A. Kalev, R. L. Kosut, and I. H. Deutsch, Quantum to- mography protocols with positivity are compressed sens- ing protocols, npj Quantum Information1, 15018 (2015)
2015
-
[34]
Gell-Mann, Symmetries of baryons and mesons, Phys
M. Gell-Mann, Symmetries of baryons and mesons, Phys. Rev.125, 1067 (1962)
1962
-
[35]
J. J. Sylvester, Papers in johns hopkins univer- sity circulars,https://archive.org/details/ collectedmathema03sylvuoft(1882), johns Hop- kins University Circulars I: 241–242 (1882); II: 46 (1883); III: 7–9 (1884). Summarized inThe Collected Mathematical Papers of James Joseph...
1909
-
[36]
J. A. Miszczak, Generating and using truly random quan- tum states in mathematica, Computer Physics Commu- nications183, 118–124 (2012)
2012
-
[37]
III C, except for one instance in which they are 8%, namely GGM basis,k= 9 two-qudit states,m/d 2 = 0.02
Typically they are less than 4% for the lowest one or two values ofm/d 2 in the simulations presented in Sec. III C, except for one instance in which they are 8%, namely GGM basis,k= 9 two-qudit states,m/d 2 = 0.02
-
[38]
Akaike, A new look at the statistical model identifica- tion, IEEE Transactions on Automatic Control19, 716 (1974)
H. Akaike, A new look at the statistical model identifica- tion, IEEE Transactions on Automatic Control19, 716 (1974)
1974
Reviewed August 15, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.