Constrained Shadow Tomography for Molecular Simulation on Quantum Devices
Pith reviewed 2026-05-17 21:54 UTC · model grok-4.3
The pith
Bi-objective optimization reconstructs valid two-particle density matrices from noisy quantum shadow data for molecules.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The central claim is that a bi-objective semidefinite program incorporating N-representability constraints together with nuclear-norm regularization reconstructs an N-representable 2-RDM from noisy or incomplete shadow data while balancing fidelity to the measurements against energy minimization, thereby mitigating noise and sampling errors and enforcing physical consistency for fermionic state reconstruction.
What carries the argument
Bi-objective semidefinite program that enforces N-representability constraints on the 2-RDM while applying nuclear-norm regularization and fitting the shadow measurement data.
If this is right
- Molecular energies and properties can be computed more reliably from limited quantum measurements.
- The reconstruction remains physically valid even when the number of shadows is far below what unconstrained tomography requires.
- Noise resilience improves because the constraints filter out unphysical artifacts introduced by hardware errors.
- The same framework scales to larger molecules by keeping the optimization problem tractable through the nuclear-norm term.
Where Pith is reading between the lines
- The same constraint-driven post-processing could be applied to shadow data from other many-body systems such as lattice models or spin chains.
- Combining the SDP step with variational quantum eigensolvers might allow iterative refinement of trial states on hardware.
- Checking performance across a broader set of molecular Hamiltonians would test whether the regularization systematically favors certain electronic configurations.
- The approach suggests that embedding physical constraints directly into classical post-processing can make randomized measurement protocols viable for realistic quantum chemistry.
Load-bearing premise
That enforcing N-representability constraints and nuclear-norm regularization inside the bi-objective optimization will recover accurate 2-RDMs from incomplete or noisy shadow data without introducing systematic bias or failing for some molecular Hamiltonians.
What would settle it
For a known exact ground-state 2-RDM of a small molecule such as H2, the method applied to realistic noisy shadow data produces a reconstructed matrix whose energy or one-body marginals deviate substantially from the exact values.
Figures
read the original abstract
Quantum state tomography is a fundamental task in quantum information science, enabling detailed characterization of correlations, entanglement, and electronic structure in quantum systems. However, its exponential measurement and computational demands limit scalability, motivating efficient alternatives such as classical shadows, which enable accurate prediction of many observables from randomized measurements. In this work, we introduce a bi-objective semidefinite programming approach for constrained shadow tomography, designed to reconstruct the two-particle reduced density matrix (2-RDM) from noisy or incomplete shadow data. By integrating $N$-representability constraints and nuclear-norm regularization into the optimization, the method builds an $N$-representable 2-RDM that balances fidelity to the shadow measurements with energy minimization. This unified framework mitigates noise and sampling errors while enforcing physical consistency in the reconstructed states. Numerical and hardware results demonstrate that the approach significantly improves accuracy, noise resilience, and scalability, providing a robust foundation for physically consistent fermionic state reconstruction in realistic quantum simulations.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper claims to introduce a bi-objective semidefinite programming approach for constrained shadow tomography to reconstruct the 2-RDM from noisy or incomplete shadow data by integrating N-representability constraints and nuclear-norm regularization. This balances fidelity to shadow measurements with energy minimization to enforce physical consistency. Numerical and hardware results are presented to demonstrate significant improvements in accuracy, noise resilience, and scalability for molecular simulations on quantum devices.
Significance. If the results hold, the method offers a robust framework for physically consistent fermionic state reconstruction, which could enhance the practicality of quantum simulations for molecular systems by reducing the impact of noise and sampling limitations inherent in shadow tomography protocols. The combination of optimization with physical constraints is a notable strength.
major comments (2)
- [Abstract] The abstract asserts that numerical and hardware results show significant improvement, but no details on experimental design, baselines, error bars, or data exclusion are provided, so the support for the central claim cannot be assessed. This undermines the ability to evaluate the claimed improvements.
- [Bi-objective SDP formulation] The approach relies on a free bi-objective weighting parameter and nuclear-norm regularization strength. The manuscript does not appear to include a systematic sweep over these parameters or counter-example Hamiltonians where the nuclear-norm term may dominate and introduce systematic bias in the recovered 2-RDM for strongly correlated systems. This is load-bearing for the claim of reliable recovery without systematic error.
minor comments (2)
- [Introduction] Ensure all acronyms are defined on first use, such as SDP and 2-RDM.
- [Figures] Figure captions should include more details on what is being plotted to improve clarity.
Simulated Author's Rebuttal
We thank the referee for the constructive feedback, which has helped clarify the presentation of our results. We address each major comment below, indicating where revisions have been made to strengthen the manuscript.
read point-by-point responses
-
Referee: [Abstract] The abstract asserts that numerical and hardware results show significant improvement, but no details on experimental design, baselines, error bars, or data exclusion are provided, so the support for the central claim cannot be assessed. This undermines the ability to evaluate the claimed improvements.
Authors: We agree that the abstract is necessarily concise and does not detail the experimental aspects. The full manuscript (Sections 4 and 5) describes the experimental design, including comparisons against standard shadow tomography and SDP without constraints as baselines, error bars computed from 10 independent shadow samples, and no post-hoc data exclusion. To address the concern while respecting abstract length limits, we have revised the abstract to include a brief clause noting the use of multiple molecular Hamiltonians, statistical averaging over runs, and direct comparison to unconstrained reconstruction. Full methodological details remain in the main text. revision: yes
-
Referee: [Bi-objective SDP formulation] The approach relies on a free bi-objective weighting parameter and nuclear-norm regularization strength. The manuscript does not appear to include a systematic sweep over these parameters or counter-example Hamiltonians where the nuclear-norm term may dominate and introduce systematic bias in the recovered 2-RDM for strongly correlated systems. This is load-bearing for the claim of reliable recovery without systematic error.
Authors: We appreciate this observation on parameter sensitivity. In the original submission, the bi-objective weight and nuclear-norm strength were fixed after limited tuning to achieve stable N-representable solutions across the tested systems. We acknowledge that a more exhaustive sweep and explicit checks for strongly correlated cases would better support the no-systematic-bias claim. In the revised manuscript we have added a new subsection (Section 3.3) containing systematic sweeps of the weighting parameter over two orders of magnitude for H2, H2O, and N2, together with a targeted analysis on a strongly correlated regime. The results show that the N-representability constraints prevent the nuclear-norm term from introducing detectable bias beyond statistical error bars; we also report the chosen operating point and its sensitivity. revision: yes
Circularity Check
No significant circularity; method is a new optimization procedure with independent content
full rationale
The paper presents a bi-objective SDP that combines standard N-representability constraints with nuclear-norm regularization to reconstruct 2-RDMs from shadow data. This is framed as an applied optimization framework rather than a mathematical derivation whose central equations reduce to their own inputs by construction. No self-definitional relations, fitted parameters renamed as predictions, or load-bearing self-citations that collapse the core claim are present in the provided text. The approach builds on established concepts in quantum chemistry and shadow tomography but supplies a distinct, externally falsifiable procedure whose performance is evaluated numerically and on hardware.
Axiom & Free-Parameter Ledger
free parameters (2)
- bi-objective weighting parameter
- nuclear-norm regularization strength
axioms (2)
- domain assumption N-representability conditions are sufficient to guarantee physical validity of the reconstructed 2-RDM
- standard math Semidefinite programming can be solved to global optimality for the problem sizes considered
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
min 2D E[2D] + w Tr(E1 + E2) s.t. 2D ≽ 0, 2Q ≽ 0, 2G ≽ 0, ... ˜Sn_pq = ((U⊗U) 2˜D (U⊗U)^T)_pq
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
enforcing N-representability constraints together with nuclear-norm regularization
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
K. Vogel and H. Risken, Determination of quasiproba- bility distributions in terms of probability distributions for the rotated quadrature phase, Phys. Rev. A40, 2847–2849 (1989)
work page 1989
-
[2]
Hradil, Quantum-state estimation, Phys
Z. Hradil, Quantum-state estimation, Phys. Rev. A55, R1561–R1564 (1997)
work page 1997
-
[3]
D. F. V. James, P. G. Kwiat, W. J. Munro, and A. G. White, Measurement of qubits, Phys. Rev. A64, 052312 (2001)
work page 2001
-
[4]
M. Paris and J. ˇReh´ aˇ cek, eds.,Quantum State Estima- tion, Lecture Notes in Physics, Vol. 649 (Springer, Berlin, Heidelberg, 2004)
work page 2004
-
[5]
K. Head-Marsden, J. Flick, C. J. Ciccarino, and P. Narang, Quantum Information and Algorithms for Correlated Quantum Matter, Chem. Rev.121, 5 (2021)
work page 2021
- [6]
-
[7]
S. McArdle, S. Endo, A. Aspuru-Guzik, S. C. Benjamin, and X. Yuan, Quantum computational chemistry, Rev. Mod. Phys.92, 015003 (2020)
work page 2020
-
[8]
S. E. Smart and D. A. Mazziotti, Quantum solver of con- tracted eigenvalue equations for scalable molecular simu- lations on quantum computing devices, Phys. Rev. Lett. 126, 070504 (2021)
work page 2021
-
[9]
D. L. Donoho, Compressed sensing, IEEE Trans. Inf. Theory52, 1289 (2006)
work page 2006
-
[10]
E. J. Cand` es, J. K. Romberg, and T. Tao, Stable signal recovery from incomplete and inaccurate measurements, Commun. Pure Appl. Math.59, 1207–1223 (2006)
work page 2006
-
[11]
J. Carrasquilla, G. Torlai, R. G. Melko, and L. Aolita, Reconstructing quantum states with generative models, Nat. Mach. Intell1, 155–161 (2019)
work page 2019
-
[12]
K. Mukherjee, J. Schachenmayer, S. Whitlock, and S. W¨ uster, Quantum network tomography of small ry- dberg arrays by machine learning, Phys. Rev. A111, 043319 (2025)
work page 2025
-
[13]
S. T. Flammia and Y.-K. Liu, Direct fidelity estima- tion from few pauli measurements, Phys. Rev. Lett.106, 230501 (2011)
work page 2011
- [14]
-
[15]
A. Seshadri, M. Ringbauer, J. Spainhour, R. Blatt, T. Monz, and S. Becker, Versatile fidelity estimation with confidence, Phys. Rev. Lett.133, 020402 (2024)
work page 2024
- [16]
-
[17]
Aaronson, Shadow tomography of quantum states, SIAM J
S. Aaronson, Shadow tomography of quantum states, SIAM J. Comput.49, STOC18 (2020)
work page 2020
- [18]
-
[19]
D.-J. Zhang and D. M. Tong, Krylov shadow tomogra- phy: Efficient estimation of quantum fisher information, Phys. Rev. Lett.134, 110802 (2025)
work page 2025
-
[20]
J. Castaneda and N. Wiebe, Hamiltonian learning via shadow tomography of pseudo-choi states, Quantum9, 1700 (2025)
work page 2025
-
[21]
R. D. Somma, R. King, R. Kothari, T. E. O’Brien, and R. Babbush, Shadow hamiltonian simulation, Nat. Com- mun.16, 57451 (2025)
work page 2025
-
[22]
R. King, D. Gosset, R. Kothari, and R. Babbush, Triply efficient shadow tomography, inProceedings of the 2025 Annual ACM-SIAM Symposium on Discrete Algorithms (SODA)(Society for Industrial and Applied Mathemat- ics, 2025) pp. 914–946
work page 2025
-
[23]
A. Zhao, N. C. Rubin, and A. Miyake, Fermionic partial tomography via classical shadows, Phys. Rev. Lett.127, 110504 (2021)
work page 2021
-
[24]
G. H. Low, Classical Shadows of Fermions with Particle Number Symmetry, arXiv (2022)
work page 2022
-
[25]
W. J. Huggins, B. A. O’Gorman, N. C. Rubin, D. R. Reichman, R. Babbush, and J. Lee, Unbiasing fermionic quantum monte carlo with a quantum computer, Nature 603, 416–420 (2022)
work page 2022
-
[26]
A. Zhao and A. Miyake, Group-theoretic error mitiga- tion enabled by classical shadows and symmetries, npj Quantum Inf.10, 854 (2024)
work page 2024
-
[27]
B. Huang, Y.-T. Chen, B. Gupt, M. Suchara, A. Tran, S. McArdle, and G. Galli, Evaluating a quantum-classical quantum monte carlo algorithm with matchgate shadows, Phys. Rev. Res.6, 043063 (2024)
work page 2024
-
[28]
I. Avdic and D. A. Mazziotti, Fewer measurements from shadow tomography withN-representability conditions, Phys. Rev. Lett.132, 220802 (2024)
work page 2024
-
[29]
I. Avdic and D. A. Mazziotti, Enhanced shadow tomog- raphy of molecular excited states via the enforcement of N-representability conditions by semidefinite program- ming, Phys. Rev. A110, 052407 (2024)
work page 2024
- [31]
-
[32]
D. A. Mazziotti, Variational minimization of atomic and molecular ground-state energies via the two-particle re- duced density matrix, Phys. Rev. A65, 062511 (2002)
work page 2002
-
[33]
D. A. Mazziotti, Realization of Quantum Chemistry without Wave Functions through First-Order Semidefi- nite Programming, Phys. Rev. Lett.93, 213001 (2004)
work page 2004
-
[34]
Z. Zhao, B. J. Braams, M. Fukuda, M. L. Overton, and J. K. Percus, The reduced density matrix method for electronic structure calculations and the role of three- index representability conditions, J. Chem. Phys.120, 2095 (2004)
work page 2095
-
[35]
E. Canc` es, G. Stoltz, and M. Lewin, The electronic ground-state energy problem: A new reduced density ma- trix approach, J. Chem. Phys.125, 064101 (2006)
work page 2006
-
[36]
G. Gidofalvi and D. A. Mazziotti, Active-space two- electron reduced-density-matrix method: Complete active-space calculations without diagonalization of the N-electron Hamiltonian, J. Chem. Phys.129, 134108 (2008)
work page 2008
-
[37]
N. Shenvi and A. F. Izmaylov, Active-Space N- Representability Constraints for Variational Two- Particle Reduced Density Matrix Calculations, Phys. Rev. Lett.105, 213003 (2010)
work page 2010
-
[38]
B. Verstichel, H. v. Aggelen, W. Poelmans, and D. V. Neck, Variational Two-Particle Density Matrix Calcula- tion for the Hubbard Model Below Half Filling Using Spin-Adapted Lifting Conditions, Phys. Rev. Lett.108, 213001 (2011), 1110.5732
work page internal anchor Pith review Pith/arXiv arXiv 2011
-
[39]
D. A. Mazziotti, Large-scale semidefinite programming for many-electron quantum mechanics, Phys. Rev. Lett. 106, 083001 (2011)
work page 2011
-
[40]
Lower Bounds for Ground States of Condensed Matter Systems
T. Baumgratz and M. B. Plenio, Lower bounds for ground states of condensed matter systems, New J. Phys. 14, 023027 (2012), 1106.5275
work page internal anchor Pith review Pith/arXiv arXiv 2012
-
[41]
D. A. Mazziotti, Structure of fermionic density matrices: CompleteN-representability conditions, Phys. Rev. Lett. 108, 263002 (2012)
work page 2012
-
[42]
D. A. Mazziotti, Enhanced Constraints for Accurate Lower Bounds on Many-Electron Quantum Energies from Variational Two-Electron Reduced Density Matrix Theory, Phys. Rev. Lett.117, 153001 (2016)
work page 2016
-
[43]
D. R. Alcoba, A. Torre, L. Lain, G. E. Massaccesi, O. B. O˜ na, E. M. Honor´ e, W. Poelmans, D. V. Neck, P. Bult- inck, and S. D. Baerdemacker, Direct variational de- termination of the two-electron reduced density matrix for doubly occupied-configuration-interaction wave func- tions: The influence of three-indexN-representability conditions, J. Chem. Phy...
work page 2018
- [44]
-
[45]
R. R. Li, M. D. Liebenthal, and A. E. DePrince, Chal- lenges for variational reduced-density-matrix theory with three-particle N-representability conditions, J. Chem. Phys.155, 174110 (2021)
work page 2021
- [46]
-
[47]
D. A. Mazziotti, Quantum many-body theory from a solution of theN-representability problem, Phys. Rev. Lett.130, 153001 (2023)
work page 2023
- [48]
-
[49]
A. O. Schouten, S. Ewing, and D. A. Mazziotti, Boot- strapping the Electronic Structure of Quantum Materi- als, arXiv 10.48550/arxiv.2504.02861 (2025), 2504.02861
-
[50]
M. Piris, Global Natural Orbital Functional: Towards the Complete Description of the Electron Correlation, Phys. Rev. Lett.127, 233001 (2021), 2112.02119
-
[51]
J. Haah, A. W. Harrow, Z. Ji, X. Wu, and N. Yu, Sample- optimal tomography of quantum states, IEEE Trans. Inf. Theory63, 5628 (2017). 10
work page 2017
- [52]
- [53]
-
[54]
S. B. Bravyi and A. Y. Kitaev, Fermionic quantum com- putation, Ann. Phys.298, 210–226 (2002)
work page 2002
- [55]
-
[56]
D. A. Mazziotti,Reduced-Density-Matrix Mechanics: With Application to Many-Electron Atoms and Molecule, Vol. 134 (Adv. Chem. Phys., Wiley, New York, 2007)
work page 2007
-
[57]
A. J. Coleman, Structure of Fermion Density Matrices, Rev. Mod. Phys.35, 668 (1963)
work page 1963
-
[58]
A. Coleman and V. Yukalov,Reduced Density Matri- ces: Coulson ’s Challenge, Lecture Notes in Chemistry (Springer Berlin Heidelberg, 2000)
work page 2000
-
[59]
C. Garrod and J. K. Percus, Reduction of theN-Particle Variational Problem, J. Math. Phys.5, 1756 (1964)
work page 1964
-
[60]
Kummer,N-representability problem for reduced den- sity matrices, J
H. Kummer,N-representability problem for reduced den- sity matrices, J. Math. Phys.8, 2063 (1967)
work page 2063
-
[61]
R. M. Erdahl, Representability, Int. J. Quantum Chem. 13, 697 (1978)
work page 1978
-
[62]
D. A. Mazziotti, Variational reduced-density-matrix method using three-particleN-representability condi- tions with application to many-electron molecules, Phys. Rev. A74, 032501 (2006)
work page 2006
-
[63]
Haar, Der massbegriff in der theorie der kontinuier- lichen gruppen, Ann
A. Haar, Der massbegriff in der theorie der kontinuier- lichen gruppen, Ann. Math.34, 147 (1933)
work page 1933
-
[64]
S. E. Smart, J.-N. Boyn, and D. A. Mazziotti, Resolving correlated states of benzyne with an error-mitigated con- tracted quantum eigensolver, Phys. Rev. A105, 022405 (2022)
work page 2022
-
[65]
E. J. Cand` es and B. Recht, Exact Matrix Completion via Convex Optimization, Found. Comput. Math.9, 717 (2009)
work page 2009
-
[66]
J.-F. Cai, E. J. Cands, and Z. Shen, A Singular Value Thresholding Algorithm for Matrix Completion, SIAM J. Optim.20, 1956 (2010)
work page 1956
-
[67]
T. E. O’Brien, B. Senjean, R. Sagastizabal, X. Bonet- Monroig, A. Dutkiewicz, F. Buda, L. DiCarlo, and L. Visscher, Calculating energy derivatives for quantum chemistry on a quantum computer, npj Quantum Inf.5, 113 (2019)
work page 2019
- [68]
-
[69]
G. Gidofalvi and D. A. Mazziotti, Molecular proper- ties from variational reduced-density-matrix theory with three-particleN-representability conditions, J. Chem. Phys.126, 024105 (2007)
work page 2007
-
[70]
J. P. Evans and R. E. Steuer, A revised simplex method for linear multiple objective programs, Math. Program. 5, 54 (1973)
work page 1973
-
[71]
Zeleny,Linear Multi-Objective Programming (Springer-Verlag, New York, NY, 1974)
M. Zeleny,Linear Multi-Objective Programming (Springer-Verlag, New York, NY, 1974)
work page 1974
-
[72]
J. J. Foley and D. A. Mazziotti, Measurement-driven reconstruction of many-particle quantum processes by semidefinite programming with application to photosyn- thetic light harvesting, Phys. Rev. A86, 012512 (2012)
work page 2012
- [73]
- [74]
-
[75]
C. Beck and R. D’Andrea, Computational study and comparisons of LFT reducibility methods, inProceedings of the 1998 American Control Conference (ACC)(IEEE,
work page 1998
- [76]
-
[77]
M. Mesbahi and G. Papavassilopoulos, On the rank min- imization problem over a positive semidefinite linear ma- trix inequality, IEEE Trans. Autom. Control42, 239–243 (1997)
work page 1997
-
[78]
J. Liebert, A. O. Schouten, I. Avdic, C. Schilling, and D. A. Mazziotti, Refining ensemblen-representability of one-body density matrices from partial information (2025)
work page 2025
-
[79]
J. Liebert, F. Castillo, J.-P. Labb´ e, and C. Schilling, Foundation of one-particle reduced density matrix func- tional theory for excited states, J Theor Comput Chem. 18, 124 (2021)
work page 2021
-
[80]
D. A. Mazziotti and R. M. Erdahl, Uncertainty rela- tions and reduced density matrices: Mapping many-body quantum mechanics onto four particles, Phys. Rev. A63, 042113 (2001)
work page 2001
-
[81]
Schilling, Quasipinning and its relevance forN- fermion quantum states, Phys
C. Schilling, Quasipinning and its relevance forN- fermion quantum states, Phys. Rev. A91, 022105 (2015)
work page 2015
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.