Correlated Purification for Restoring N-Representability in Quantum Simulation
Pith reviewed 2026-05-17 21:53 UTC · model grok-4.3
The pith
Correlated purification restores N-representability to noisy two-electron RDMs via bi-objective semidefinite programming.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The correlated purification framework based on semidefinite programming restores accuracy to noisy, unphysical two-electron RDMs by performing a bi-objective optimization that minimizes both the many-electron energy and the nuclear norm of the change in the measured 2-RDM, yielding substantial reductions in energy and 2-RDM error and achieving chemical accuracy in applications to large hydrogen chains.
What carries the argument
Bi-objective semidefinite-programming optimization that minimizes the many-electron energy together with the nuclear norm of the correction to the measured 2-RDM, thereby promoting low-rank physically valid adjustments.
If this is right
- Substantial reductions occur in both energy and 2-RDM error for fermionic shadow tomography of large hydrogen chains.
- Chemical accuracy is reached across the full dissociation curve of the studied systems.
- The same procedure applies to excited and non-stationary states simply by decreasing the energy weight relative to the error norm.
- The approach supplies a robust post-processing strategy for any many-body quantum simulation that produces noisy 2-RDM estimates.
Where Pith is reading between the lines
- The method may be combined with hardware-specific error-mitigation techniques to improve performance on actual quantum processors.
- Testing the same purification pipeline on molecules other than hydrogen chains would reveal how broadly the nuclear-norm regularizer succeeds.
- Varying the relative weight between energy and nuclear-norm terms could be optimized per observable to minimize error on properties beyond the ground-state energy.
Load-bearing premise
The bi-objective optimization with tunable energy weight and nuclear-norm regularization can recover physically valid 2-RDMs from realistic noise levels without introducing artifacts that invalidate the energy or other observables.
What would settle it
Applying the procedure to shadow-tomography data from a hydrogen chain and finding that the output energy lies outside chemical accuracy or that the purified 2-RDM violates a basic N-representability condition such as positive semidefiniteness would falsify the central claim.
Figures
read the original abstract
Classical shadow tomography offers a scalable route to estimating properties of quantum states, but the resulting reduced density matrices (RDMs) often violate constraints that ensure they represent $N$-electron states -- known as $N$-representability conditions -- because of statistical and hardware noise. We present a correlated purification framework based on semidefinite programming to restore accuracy to these noisy, unphysical two-electron RDMs. The method performs a bi-objective optimization that minimizes both the many-electron energy and the nuclear norm of the change in the measured 2-RDM. The nuclear norm, often employed in matrix completion, promotes low-rank, physically meaningful corrections to the 2-RDM, while the energy term acts as a regularization term that can improve the purity of the ground state. While the method is particularly effective for the ground state, it can also be applied to excited and non-stationary states by decreasing the weight of the energy relative to the error norm. In an application to fermionic shadow tomography of large hydrogen chains, correlated purification yields substantial reductions in both energy and 2-RDM error, achieving chemical accuracy across dissociation curves. This framework provides a robust strategy for tomography in many-body quantum simulations.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript proposes a correlated purification framework using semidefinite programming to restore N-representability to noisy 2-RDMs from fermionic shadow tomography. It solves a bi-objective optimization minimizing the many-electron energy computed from the corrected 2-RDM plus a nuclear-norm penalty on the correction matrix Delta. The method is applied to hydrogen-chain simulations, claiming substantial reductions in energy and 2-RDM errors that achieve chemical accuracy across dissociation curves, with the energy term acting as regularization that is particularly effective for ground states.
Significance. If the central claims hold without bias, the framework could provide a practical post-processing tool for improving accuracy in shadow tomography for quantum many-body simulations, leveraging standard SDP and nuclear-norm techniques to enforce physical constraints on noisy data. The approach extends matrix-completion ideas to fermionic systems and offers a tunable method for ground versus excited states.
major comments (1)
- [Abstract and optimization formulation] Abstract and optimization formulation: the bi-objective SDP minimizes the energy (computed from the corrected 2-RDM using the Hamiltonian) plus lambda * nuclear_norm(Delta). This explicitly incorporates the Hamiltonian as a prior during purification rather than enforcing N-representability constraints alone. For the central claim of chemical accuracy on hydrogen chains, the reported energy-error reductions could arise from variational pulling toward lower energies instead of faithful recovery from noise; the nuclear-norm term alone does not guarantee unbiased observables when the energy weight is nonzero.
minor comments (2)
- [Results on hydrogen chains] Results on hydrogen chains: the performance gains are stated without error bars, details on the number of shadow tomography shots, or explicit baseline comparisons to standard N-representability projections that omit the energy term.
- [Method and parameter selection] Parameter selection: the relative weight between the energy term and nuclear-norm error is described as tunable but no systematic procedure or sensitivity analysis for its choice is provided.
Simulated Author's Rebuttal
We thank the referee for their careful reading and constructive feedback on our manuscript. We address the major comment regarding the optimization formulation below.
read point-by-point responses
-
Referee: [Abstract and optimization formulation] Abstract and optimization formulation: the bi-objective SDP minimizes the energy (computed from the corrected 2-RDM using the Hamiltonian) plus lambda * nuclear_norm(Delta). This explicitly incorporates the Hamiltonian as a prior during purification rather than enforcing N-representability constraints alone. For the central claim of chemical accuracy on hydrogen chains, the reported energy-error reductions could arise from variational pulling toward lower energies instead of faithful recovery from noise; the nuclear-norm term alone does not guarantee unbiased observables when the energy weight is nonzero.
Authors: We agree that the objective explicitly includes the many-electron energy, which incorporates information from the Hamiltonian. However, the semidefinite program enforces N-representability conditions as hard constraints on the corrected 2-RDM, so that the output is guaranteed to correspond to a valid N-electron state. The energy term functions as a regularization that selects, among the feasible N-representable matrices, the one closest to the ground state, while the nuclear-norm penalty on Delta ensures the correction remains minimal and low-rank. This combination is intentional for ground-state applications, as stated in the manuscript, and the weight on the energy term can be reduced or set to zero for excited or non-stationary states. To address the concern of bias, we note that the reported improvements include direct reductions in 2-RDM error metrics (independent of the Hamiltonian), which are not part of the objective and therefore provide evidence of improved fidelity beyond variational energy lowering. We will add a clarifying paragraph in the Methods section and an ablation with zero energy weight to make this distinction explicit. revision: yes
Circularity Check
No significant circularity; method is a standard SDP construction with independent numerical validation
full rationale
The paper presents a bi-objective SDP that enforces N-representability constraints while using nuclear-norm regularization on the correction and an energy term as tunable regularization. The reported error reductions on hydrogen-chain shadow tomography are numerical outcomes of applying this procedure to noisy input 2-RDMs; they are not obtained by renaming or re-fitting quantities already defined in the input data or by a self-citation chain that supplies the central result. The framework remains self-contained against external benchmarks because the N-representability conditions and nuclear-norm objective are drawn from established convex optimization literature rather than from quantities defined inside the present work.
Axiom & Free-Parameter Ledger
free parameters (1)
- relative weight of energy term versus nuclear-norm error
axioms (1)
- domain assumption Nuclear norm promotes low-rank, physically meaningful corrections to the 2-RDM
Reference graph
Works this paper leans on
-
[1]
B. P. Lanyon, J. D. Whitfield, G. G. Gillett, M. E. Goggin, M. P. Almeida, I. Kassal, J. D. Biamonte, M. Mohseni, B. J. Powell, M. Barbieri,et al., Towards quantum chemistry on a quantum computer, Nat. Chem. 2, 106 (2010)
work page 2010
-
[2]
I. M. Georgescu, S. Ashhab, and F. Nori, Quantum sim- ulation, Rev. Mod. Phys.86, 153 (2014)
work page 2014
-
[3]
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
- [4]
- [5]
-
[6]
A. J. Coleman and V. I. Yukalov,Reduced Density Ma- trices: Coulson ’s Challenge(Springer, Berlin, 2000)
work page 2000
-
[7]
D. A. Mazziotti,Reduced-Density-Matrix Mechan- ics: With Application to Many-Electron Atoms and Molecules, Advances in Chemical Physics (Wiley, 2007)
work page 2007
-
[8]
S. Aaronson, Shadow tomography of quantum states, in Proceedings of the 50th annual ACM SIGACT sympo- sium on theory of computing(2018) pp. 325–338
work page 2018
- [9]
-
[10]
A. Zhao, N. C. Rubin, and A. Miyake, Fermionic par- tial tomography via classical shadows, Physical Review Letters127, 110504 (2021)
work page 2021
-
[11]
G. H. Low, Classical shadows of fermions with parti- cle number symmetry, arXiv preprint arXiv:2208.08964 10.48550/arXiv.2208.08964 (2022)
-
[12]
A. J. Coleman, Structure of fermion density matrices, Rev. Mod. Phys.35, 668 (1963)
work page 1963
-
[13]
R. M. Erdahl, Representability, Int. J. Quantum Chem. 13, 697 (1978)
work page 1978
-
[14]
J. E. Harriman, Geometry of density matrices. ii. reduced density matrices andNrepresentability, Phys. Rev. A 17, 1257 (1978)
work page 1978
-
[15]
D. A. Mazziotti, Structure of fermionic density matrices: CompleteN-representability conditions, Phys. Rev. Lett. 108, 263002 (2012)
work page 2012
- [16]
-
[17]
McWeeny, Some recent advances in density matrix theory, Reviews of Modern Physics32, 335 (1960)
R. McWeeny, Some recent advances in density matrix theory, Reviews of Modern Physics32, 335 (1960)
work page 1960
-
[18]
Goedecker, Linear scaling electronic structure meth- ods, Rev
S. Goedecker, Linear scaling electronic structure meth- ods, Rev. Mod. Phys.71, 1085 (1999)
work page 1999
-
[19]
D. A. Mazziotti, Purification of correlated reduced den- sity matrices, Phys. Rev. E65, 026704 (2002)
work page 2002
-
[20]
E. Pescoller, M. Eder, and I. Bˇ rezinov´ a, Projective pu- rification of correlated reduced density matrices, Phys. Rev. Res.7, 013211 (2025)
work page 2025
-
[21]
E. Candes and B. Recht, Exact matrix completion via convex optimization, Communications of the ACM55, 111 (2012)
work page 2012
-
[22]
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
-
[23]
S. E. Smart, J.-N. Boyn, and D. A. Mazziotti, Resolving correlated states of benzyne with an error-mitigated con- 7 tracted quantum eigensolver, Phys. Rev. A105, 022405 (2022)
work page 2022
-
[24]
N. C. Rubin, R. Babbush, and J. McClean, Application of fermionic marginal constraints to hybrid quantum al- gorithms, New J. Phys.20, 053020 (2018)
work page 2018
-
[25]
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
- [26]
-
[27]
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
-
[28]
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
-
[29]
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
-
[30]
The Electronic Ground State Energy Problem: a New Reduced Density Matrix Approach
E. Canc` es, G. Stoltz, and M. Lewin, The electronic ground-state energy problem: A new reduced density matrix approach, J. Chem. Phys.125, 064101 (2006), quant-ph/0602042
work page internal anchor Pith review Pith/arXiv arXiv 2006
-
[31]
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
-
[32]
N. Shenvi and A. F. Izmaylov, Active-SpaceN- representability Constraints for Variational Two-Particle Reduced Density Matrix Calculations, Phys. Rev. Lett. 105, 213003 (2010)
work page 2010
-
[33]
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
-
[34]
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
-
[35]
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
-
[36]
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
- [37]
-
[38]
R. R. Li, M. D. Liebenthal, and A. E. DePrince, Chal- lenges for variational reduced-density-matrix theory with three-particleN-representability conditions, J. Chem. Phys.155, 174110 (2021)
work page 2021
- [39]
- [40]
-
[41]
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
-
[42]
M. Piris, Global Natural Orbital Functional: Towards the Complete Description of the Electron Correlation, Phys. Rev. Lett.127, 233001 (2021), 2112.02119
-
[43]
Y. Wang, I. Avdic, and D. A. Mazziotti, Shadow ansatz for the many-fermion wave function in scalable molecular simulations on quantum computing devices, Phys. Rev. A112, 022432 (2025)
work page 2025
- [44]
-
[45]
K. Head-Marsden, J. Flick, C. J. Ciccarino, and P. Narang, Quantum Information and Algorithms for Correlated Quantum Matter, Chem. Rev.121, 3061 (2021)
work page 2021
-
[46]
A. Peruzzo, J. McClean, P. Shadbolt, M.-H. Yung, X.-Q. Zhou, P. J. Love, A. Aspuru-Guzik, and J. L. O’Brien, A variational eigenvalue solver on a photonic quantum processor, Nature Commun.5, 4213 (2014)
work page 2014
-
[47]
Y. Cao, A. Papageorgiou, I. Petras, J. Traub, and S. Kais, Quantum algorithm and circuit design solving the Pois- son equation, New Journal of Physics15, 013021 (2013)
work page 2013
-
[48]
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
-
[49]
C. Garrod and J. K. Percus, Reduction of theN-Particle Variational Problem, J. Math, Phys.5, 1756 (1964)
work page 1964
-
[50]
L. Vandenberghe and S. Boyd, Semidefinite program- ming, SIAM Rev.38, 49 (1996)
work page 1996
-
[51]
J. J. Foley IV 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
-
[52]
I. Avdic and D. A. Mazziotti, Fewer Measurements from Shadow Tomography withN-Representability Condi- tions, Phys. Rev. Lett.132, 220802 (2024)
work page 2024
-
[53]
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
- [54]
-
[55]
Y. Wang and D. A. Mazziotti, Electronic excited states from a variance-based contracted quantum eigensolver, Phys. Rev. A108, 022814 (2023)
work page 2023
- [56]
-
[57]
Q. Sun, T. C. Berkelbach, N. S. Blunt, G. H. Booth, 8 S. Guo, Z. Li, J. Liu, J. D. McClain, E. R. Sayfutyarova, S. Sharma, S. Wouters, and G. K. Chan, PySCF: the Python-based simulations of chemistry framework, Wi- ley Interdisciplinary Reviews: Computational Molecular Science8, 10.1002/wcms.1340 (2018)
-
[58]
The ffsim developers, ffsim: Faster simulations of fermionic quantum circuits,https://github.com/ qiskit-community/ffsim
-
[59]
Qiskit contributors, Qiskit: An open-source framework for quantum computing (2025)
work page 2025
-
[60]
D. A. Mazziotti, Large-Scale Semidefinite Programming for Many-Electron Quantum Mechanics, Phys. Rev. Lett. 106, 083001 (2011). [61]Maple(Maplesoft, Waterloo, 2025). [62]Quantum Chemistry Toolbox in Maple(Maplesoft, Wa- terloo, 2025)
work page 2011
-
[61]
W. J. Hehre, R. F. Stewart, and J. A. Pople, Self- consistent molecular-orbital methods. i. use of gaussian expansions of slater-type atomic orbitals, J. Chem. Phys. 51, 2657 (1969)
work page 1969
-
[62]
Both numerical data and a sample computational script are openly available athttps://github.com/damazz/ CorrelatedPurification
-
[63]
P. Jordan and E. Wigner, ¨Uber das Paulische ¨Aquivalenzverbot, Z. Physik47, 631 (1928)
work page 1928
- [64]
-
[65]
D. A. Mazziotti, Variational method for solving the con- tracted Schr¨ odinger equation through a projection of the N-particle power method onto the two-particle space, J. Chem. Phys.116, 1239 (2002)
work page 2002
-
[66]
Propagating two-particle reduced density matrices without wavefunctions
F. Lackner, I. Bˇ rezinov´ a, T. Sato, K. L. Ishikawa, and J. Burgd¨ orfer, Propagating two-particle reduced den- sity matrices without wave functions, Phys. Rev. A91, 023412 (2015), 1411.0495
work page internal anchor Pith review Pith/arXiv arXiv 2015
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.