Full-Spectrum Quantum Simulation for the Nuclear Shell Model

B. Maheshwari; P. Stevenson; P. Van Isacker

REVIEW 2 major objections 4 minor 23 references

A single quantum optimization recovers the full low-lying shell-model spectrum for few nucleons and restores total angular momentum.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · grok-4.5

2026-07-10 10:52 UTC pith:4VCMHILK

load-bearing objection Solid noiseless demo of single-run ADAPT-SSVQE recovering full few-nucleon spectra and J^{2} for j=9/2, but the systems are classically trivial and hardware path is still open. the 2 major comments →

arxiv 2607.08235 v1 pith:4VCMHILK submitted 2026-07-09 nucl-th physics.comp-phquant-ph

Full-Spectrum Quantum Simulation for the Nuclear Shell Model

B. Maheshwari , P. Stevenson , P. Van Isacker This is my paper

classification nucl-th physics.comp-phquant-ph PACS 21.60.Cs03.67.Ac

keywords nuclear shell modelquantum simulationADAPT-VQESSVQEJordan-Wigner transformationangular momentum restorationmany-body Hilbert space

verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The nuclear shell model gives a complete quantum description of nuclei, but its Hilbert space grows so large that classical computers cannot store or diagonalize the Hamiltonian for realistic valence spaces. This paper introduces a single-run quantum algorithm that finds the ground state and many excited states of the shell-model Hamiltonian at once. The Hamiltonian is mapped to qubits with the Jordan-Wigner transformation so that fermionic statistics are preserved; an adaptive ansatz built from magnetic-projection-conserving excitations then evolves a set of orthogonal reference states under one shared circuit. In noiseless simulation the method reproduces exact energies for two and three identical nucleons in a j=9/2 orbital to spectroscopic accuracy and yields states that are eigenstates of total angular momentum. A reader cares because complete low-lying spectra, not just ground states, are required for nuclear structure and applications such as decay-heat estimates, and the approach is designed to scale with qubit number rather than classical storage.

Core claim

Using ADAPT-SSVQE with an M_J-conserving double-excitation operator pool on a Jordan-Wigner mapped shell-model Hamiltonian, a single optimization simultaneously produces five (respectively ten) mutually orthogonal eigenstates for two (respectively three) identical nucleons in the 0g9/2 orbital inside a 10-qubit space, matching exact diagonalization to spectroscopic accuracy and intrinsically restoring total angular momentum J^2 within the chosen M_J subspace.

What carries the argument

ADAPT-SSVQE: a weighted subspace-search variational eigensolver whose ansatz is grown iteratively by selecting, from a pool of M_J-conserving single- and double-excitation operators, the operator of largest weighted energy gradient, then re-optimizing all parameters so that one unitary maps a set of orthogonal Hartree-Fock references onto the full targeted spectrum.

Load-bearing premise

The central claim rests on the premise that noiseless simulation of ten-qubit systems of only two or three nucleons, stopped when the gradient or energy change falls below a fixed numerical threshold, is enough evidence that the same procedure will keep spectroscopic accuracy and angular-momentum restoration for larger valence spaces or on real noisy hardware.

What would settle it

Apply the identical ADAPT-SSVQE procedure, same operator pool and convergence thresholds, to a four-nucleon or multi-j system whose exact spectrum is still classically known; check whether every targeted energy still matches the exact value to spectroscopic accuracy and whether measured J^2 values remain integer or half-integer once the circuit contains thousands of CNOT gates under a realistic noise model.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

Multiple shell-model eigenstates inside a fixed M_J subspace can be obtained from one variational optimization without separate runs or post-processing.
Total-angular-momentum quantum numbers emerge from the converged wave functions even though the operator pool never enforces J^2.
Circuit depth grows only with the number of selected operators needed to exhaust the energy gradient, giving a compact representation of nucleonic correlations.
The same single-run strategy is in principle applicable to any many-body problem that requires a full low-lying spectrum rather than a ground state alone.
The limiting resource becomes qubit count and two-qubit gate depth instead of classical Hilbert-space storage.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

If further pool truncation or hardware-efficient ansätze can cut CNOT depth, the method could become usable on near-term devices for nuclei beyond exact classical diagonalization.
The small residual J^2 deviations already visible in the three-nucleon results imply that tighter thresholds or an explicit J^2 term may be required before spin assignments are trusted for odd-mass systems.
Adding proton-neutron operators and isospin conservation would test whether single-run symmetry restoration survives in realistic even-even and odd-odd nuclei.
The two-nucleon wave-function components matching pure Clebsch-Gordan coefficients show that the algorithm reconstructs geometric angular-momentum coupling even when a realistic residual interaction is present.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Desk Editor's Note

Solid noiseless demo of single-run ADAPT-SSVQE recovering full few-nucleon spectra and J^{2} for j=9/2, but the systems are classically trivial and hardware path is still open.

read the letter

This paper cleanly shows that ADAPT-SSVQE with an M_J-conserving double-excitation pool can extract the entire low-lying spectrum for two and three identical nucleons in 0g9/2 (five and ten states) in one optimization run, matching exact JUN45 diagonalization to spectroscopic accuracy and restoring the expected J values. That is the concrete new result: simultaneous multi-state extraction plus post-convergence J^{2} and wave-function diagnostics that line up with CG/CFP structure.

What it does well is the execution. Jordan-Wigner mapping is standard and correctly preserves anti-commutation. The adaptive growth of the ansatz, weighted SSVQE loss, and L-BFGS-B re-optimization are applied without internal contradiction. Energies go to residual zero for the two-nucleon case; relative errors plunge for three nucleons after the expected early oscillations. Orthogonality is automatic from the unitary. The two-nucleon state vectors are pure geometric pairing; the three-nucleon ones show the expected seniority mixing. Citations cover the prior ADAPT-VQE shell-model work, the 58Ni ansatz paper, and the authors’ own pairing conference note. Circularity is low because everything is checked against independent classical exact diagonalization of the same Hamiltonian.

Soft spots are real but proportionate to the claim as written. The Hilbert spaces are tiny (classically trivial). Only noiseless state-vector simulation is shown; CNOT depth already hits 480–2768, and the authors themselves flag residual J^{2} deviations for the three-nucleon case as threshold artifacts. Free parameters (weights, gradient/energy cutoffs, HF reference set) exist but are conventional. No code or hardware run is released. Novelty is moderate—combination of known tools—but the full write-up with nuclear diagnostics is useful.

This is for people already working on quantum algorithms for nuclear structure or wanting a reproducible few-body benchmark. Math and data look solid within the stated scope. I would send it to peer review; the contribution is accept-shaped once residual symmetry breaking and scaling are discussed more explicitly. Engage if that is your subfield; otherwise skim the figures and move on.

Referee Report

2 major / 4 minor

Summary. The paper introduces ADAPT-SSVQE for simultaneous extraction of multiple nuclear shell-model eigenstates in one optimization run. The effective Hamiltonian (one- and two-body terms) is mapped to qubits by the Jordan-Wigner transformation. An adaptive ansatz is built from an M_J-conserving double-excitation operator pool; a weighted SSVQE loss is minimized so that a single unitary evolves a set of orthogonal Hartree-Fock references into the lowest states of a chosen M_J subspace. Benchmarks are performed in noiseless state-vector simulation for two and three identical nucleons in the 0g9/2 orbital (10 qubits, JUN45 interaction). Five (ten) mutually orthogonal states are recovered with spectroscopic accuracy relative to exact diagonalization; measured ⟨Ĵ²⟩ values restore the expected integer/half-integer spectrum (exactly for two nucleons, approximately for three).

Significance. If the numerical claims hold, the work supplies a compact, symmetry-aware route to full low-lying shell-model spectra on quantum devices—an important capability because nuclear applications (decay heat, spectroscopy) require excited states, not only ground states. Strengths that should be credited are the direct comparison with exact diagonalization, the recovery of Clebsch-Gordan/CFP structure in the final state vectors, the unitary preservation of orthogonality without post-processing, and the explicit construction and measurement of Ĵ². The demonstration remains limited to noiseless few-nucleon single-j systems; circuit depths already reach hundreds to thousands of CNOTs, so hardware relevance and multi-j scalability are prospective. Within its stated scope the evidence is solid and the method is a useful, carefully scoped advance.

major comments (2)

[Abstract and §2] Abstract vs. §2: the abstract advertises a “symmetry-preserving single and double-excitation operator pool,” yet the algorithm description and all numerical work employ only M_J-conserving double excitations T_cd^ab. The discrepancy must be resolved. If singles were deliberately omitted, the authors should state why a double-only pool remains expressive enough to reach spectroscopic accuracy and approximate Ĵ² restoration for the three-nucleon system.
[§3, Table 2] §3 and Table 2 (three-nucleon case): residual ⟨Ĵ²⟩ deviations of order 1–2 units (e.g. 22.59 vs. exact 24.75 for the lowest J=9/2 candidate) are ascribed to the loose stopping criterion G_iter or ΔE < 10^{-2}. Because the central claim includes “intrinsically restores total angular momentum (Ĵ²) symmetry,” the manuscript should either tighten the threshold and re-report both energies and ⟨Ĵ²⟩, or qualify the restoration claim as approximate under the chosen convergence settings.

minor comments (4)

[§2] Repeated typographical error “ADAPT-SSQVE” / “SSQVE” appears in §2; correct to SSVQE throughout.
[Figs. 2–3] Fig. 2 and Fig. 3 captions and panels are clear, but the CNOT-count scaling (48 CNOTs per selected multi-Pauli operator) is stated only in the text; a brief note in the figure legends would help readers.
[§3] The precise list of ten initial Hartree-Fock references for the three-nucleon M_J=1/2 subspace is given, yet no statement is made about how completeness of that set was verified; a short remark would improve reproducibility.
[References] Reference [16] is listed as “Acta Physica Polonica B19, 1-A2 (2026)”; confirm the year and volume once the conference proceedings are finalized.

Circularity Check

0 steps flagged

No significant circularity: energies and J^{2} are independently benchmarked against exact diagonalization of the same JUN45 Hamiltonian; wave-function components match geometric CG/CFP structure not fitted from the optimizer.

full rationale

The paper maps a standard shell-model Hamiltonian (JUN45 interaction) via Jordan-Wigner, then applies the known ADAPT-VQE + SSVQE combination with an M_J-conserving double-excitation pool. The weighted loss (Eq. 4) and gradient selection (Eq. 5) are standard variational constructions; they do not encode the target spectrum by definition. All reported energies are compared to classical exact diagonalization of the identical Hamiltonian, residual errors are driven to zero (two-nucleon) or to spectroscopic accuracy (three-nucleon), and measured ⟨J^{2}⟩ values recover the expected integer/half-integer spectrum. State-vector components are further checked against Clebsch-Gordan coefficients and coefficients of fractional parentage that are purely geometric and independent of the interaction. The sole self-citation is a prior conference note on a four-qubit pairing Hamiltonian; it is not load-bearing for the present 10-qubit shell-model results. Within the paper’s stated scope (noiseless state-vector simulation of few-nucleon systems) the derivation is self-contained and externally falsifiable. Score 1 reflects only the minor, non-load-bearing self-citation.

Axiom & Free-Parameter Ledger

3 free parameters · 5 axioms · 0 invented entities

The central claim rests on standard second-quantized shell-model physics, the Jordan-Wigner encoding, the variational SSVQE loss, and an adaptive operator-selection heuristic. Free parameters are algorithmic thresholds and weights, not physics fits to the reported spectra. No new particles or forces are postulated; the ‘invented’ content is the specific ADAPT-SSVQE workflow and M_J-restricted pool for this nuclear problem.

free parameters (3)

SSVQE state weights w_i = w_i = k - i
Chosen by hand as w_i = k − i (descending integers) to prioritize lower states; not derived from a uniqueness theorem.
ADAPT gradient / energy convergence thresholds = 10^{-2}
Iteration stops when max weighted gradient falls below 10^{-2} or consecutive energy change ΔE < 10^{-2}; these cutoffs are simulation choices that affect reported residual J² errors.
Initial Hartree-Fock reference set
Specific orthogonal m-scheme Slater determinants are selected by hand as the k starting states for each M_J subspace; different choices could change early gradients and circuit growth.

axioms (5)

domain assumption Effective shell-model Hamiltonian truncated to one- and two-body terms in the m-scheme is an adequate description of the low-lying spectrum for the chosen valence space.
Eq. (1) and JUN45 adoption; standard nuclear theory premise, not re-derived here.
standard math Jordan-Wigner mapping with Z-strings preserves fermionic anticommutation on qubits.
Eqs. (2)–(3); classical result used without re-proof.
domain assumption Minimizing the weighted sum of expectation values over k orthogonal inputs evolved by a shared unitary yields the lowest k eigenstates of H (SSVQE principle).
Eq. (4) and citation to Nakanishi et al.; assumed valid for the nuclear spectrum.
ad hoc to paper An operator pool restricted to M_J-conserving double excitations is expressive enough to reach spectroscopic accuracy and restore full SU(2) J² within the target manifold.
Stated in §2; empirically supported for two nucleons, only approximately for three nucleons.
ad hoc to paper Noiseless state-vector simulation of the hybrid quantum-classical loop is a valid proxy for the algorithm’s correctness (hardware noise deferred).
All results in §3 are exact simulation; conclusion flags hardware noise as future work.

pith-pipeline@v1.1.0-grok45 · 15729 in / 3598 out tokens · 41671 ms · 2026-07-10T10:52:24.401334+00:00 · methodology

0 comments

read the original abstract

The nuclear shell model is a general way of expressing the many-body nuclear Hamiltonian and deciphering the underlying nuclear structure. In today's era of modern and high-power computation, the primary limitation of the nuclear shell model is the enormous dimensionality of its Hilbert space, which far exceeds available storage capacity and prevents the diagonalization of the full Hamiltonian matrix in that space. Quantum computing offers a scalable solution to bypass this curse of dimensionality. In this work, we introduce a single-run quantum simulation capable of obtaining multiple shell-model eigenstates simultaneously. The nuclear Hamiltonian is transformed from a bit to a qubit basis using the Jordan-Wigner transformation, explicitly preserving fermionic anti-commutation. We employ a Subspace Search Variational Quantum Eigensolver (SSVQE) along with an Adaptive Derivative-Assembled Pseudo-Trotter (ADAPT) ansatz to construct the quantum circuit required to solve the shell-model problem. The ADAPT-SSVQE algorithm uses a symmetry-preserving single and double-excitation operator pool and optimizes a weighted energy sum to obtain the simultaneous convergence of all eigenstates within a targeted MJ subspace, eliminating the need for post-processing efforts to extract excited spectra. We benchmark this approach by solving the problem for two and three identical nucleons in a j = 9/2 orbital, successfully extracting five and ten mutually orthogonal states, respectively, within a 10-qubit active space. The algorithm achieves spectroscopic accuracy, in simulation, relative to exact diagonalization and intrinsically restores total angular momentum (\hat{J}^2) symmetry.

Figures

Figures reproduced from arXiv: 2607.08235 by B. Maheshwari, P. Stevenson, P. Van Isacker.

**Figure 2.** Figure 2: Variation of (a) simulated energies (in MeV) for ground and excited states shown in di [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗

**Figure 3.** Figure 3: Variation of (a) energies (in MeV) for ground and excited states shown in di [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗

discussion (0)

Reference graph

Works this paper leans on

23 extracted references · 23 canonical work pages

[1]

P. J. Brussaard and P. M. W. Glaudemans, Shell Model Applications in Nuclear Spectroscopy, North-Holland Publishing Company (1977)

work page 1977
[2]

J. R. McClean, J. Romero, R. Babbush, and A. Aspuru- Guzik, New J. Phys.18, 023023 (2016)

work page 2016
[3]

Bharti,et al.Rev

K. Bharti,et al.Rev. Mod. Phys.94, 015004 (2022)

work page 2022
[4]

Lacroix, Phys

D. Lacroix, Phys. Rev. Lett.125, 230502 (2020). 6

work page 2020
[5]

E. A. Ruiz Guzman and D. Lacroix, Phys. Rev. C105, 024324 (2022)

work page 2022
[6]

E. A. Ruiz Guzman and D. Lacroix, Eur. Phys. J. A60, 112 (2024)

work page 2024
[7]

Gibbs, Z

J. Gibbs, Z. Holmes, and P. Stevenson Quantum Machine Intelligence7, 14 (2025)

work page 2025
[8]

Anandet al.Chem

A. Anandet al.Chem. Soc. Rev.51, 1659 (2022)

work page 2022
[9]

H. R. Grimsley, S. E. Economou, E. Barnes,et al.Nat Commun.10, 3007 (2019)

work page 2019
[10]

Pérez-Obiol, A

A. Pérez-Obiol, A. M. Romero, J. Menéndez,et al.Sci Rep13, 12291 (2023)

work page 2023
[11]

Bhoy and P

B. Bhoy and P. Stevenson, New Journal of Physics, 26(7), 075001 (2024)

work page 2024
[12]

Jiang and J

C. Jiang and J. Pei, Phys. Rev. C107, 044308 (2023)

work page 2023
[13]

Zhang, D

J. Zhang, D. Lacroix, and Y . Beaujeault-Taudière, Phys. Rev. C110, 064320 (2024)

work page 2024
[14]

Zhang and D

J. Zhang and D. Lacroix, Phys. Lett. B869, 139841 (2025)

work page 2025
[15]

D. Liu, W. Du, L. Lin, J. P. Vary, and C. Yang, J. Comput. Sci.85, 102480 (2025)

work page 2025
[16]

Maheshwari, P

B. Maheshwari, P. Stevenson and P. Van Isacker, Acta Physica Polonica B19, 1-A2 (2026)

work page 2026
[17]

Algora, B

A. Algora, B. Rubio, J. Luis Tan and W. Gelletly, To- tal Absorption Technique for Nuclear Structure and Ap- plications, Springer Briefs in Physics, Springer Nature (2025)

work page 2025
[18]

Nakanishi, K

Ken M. Nakanishi, K. Mitarai and K. Fujii, Phys. Rev. Research1, 033602 (2019)

work page 2019
[19]

Jordan and E

P. Jordan and E. Wigner, Zeitschrift für Physik47 (9), 631 (1928)

work page 1928
[20]

Wong, Introduction to Classical and Quantum Com- puting, Rooted Grove, Omaha (2022)

T. Wong, Introduction to Classical and Quantum Com- puting, Rooted Grove, Omaha (2022)

work page 2022
[21]

github.io/qiskit-algorithms/stubs/qiskit_ algorithms.optimizers.L_BFGS_B.html

L-BFGS-B algorithm,https://qiskit-community. github.io/qiskit-algorithms/stubs/qiskit_ algorithms.optimizers.L_BFGS_B.html

work page
[22]

Honma, T

M. Honma, T. Otsuka, T. Mizusaki, and M. Hjorth- Jensen, Phys. Rev. C80, 064323 (2009)

work page 2009
[23]

De-Shalit and I

A. De-Shalit and I. Talmi, Nuclear Shell Theory, Aca- demic Press, New York (1963). 7

work page 1963

[1] [1]

P. J. Brussaard and P. M. W. Glaudemans, Shell Model Applications in Nuclear Spectroscopy, North-Holland Publishing Company (1977)

work page 1977

[2] [2]

J. R. McClean, J. Romero, R. Babbush, and A. Aspuru- Guzik, New J. Phys.18, 023023 (2016)

work page 2016

[3] [3]

Bharti,et al.Rev

K. Bharti,et al.Rev. Mod. Phys.94, 015004 (2022)

work page 2022

[4] [4]

Lacroix, Phys

D. Lacroix, Phys. Rev. Lett.125, 230502 (2020). 6

work page 2020

[5] [5]

E. A. Ruiz Guzman and D. Lacroix, Phys. Rev. C105, 024324 (2022)

work page 2022

[6] [6]

E. A. Ruiz Guzman and D. Lacroix, Eur. Phys. J. A60, 112 (2024)

work page 2024

[7] [7]

Gibbs, Z

J. Gibbs, Z. Holmes, and P. Stevenson Quantum Machine Intelligence7, 14 (2025)

work page 2025

[8] [8]

Anandet al.Chem

A. Anandet al.Chem. Soc. Rev.51, 1659 (2022)

work page 2022

[9] [9]

H. R. Grimsley, S. E. Economou, E. Barnes,et al.Nat Commun.10, 3007 (2019)

work page 2019

[10] [10]

Pérez-Obiol, A

A. Pérez-Obiol, A. M. Romero, J. Menéndez,et al.Sci Rep13, 12291 (2023)

work page 2023

[11] [11]

Bhoy and P

B. Bhoy and P. Stevenson, New Journal of Physics, 26(7), 075001 (2024)

work page 2024

[12] [12]

Jiang and J

C. Jiang and J. Pei, Phys. Rev. C107, 044308 (2023)

work page 2023

[13] [13]

Zhang, D

J. Zhang, D. Lacroix, and Y . Beaujeault-Taudière, Phys. Rev. C110, 064320 (2024)

work page 2024

[14] [14]

Zhang and D

J. Zhang and D. Lacroix, Phys. Lett. B869, 139841 (2025)

work page 2025

[15] [15]

D. Liu, W. Du, L. Lin, J. P. Vary, and C. Yang, J. Comput. Sci.85, 102480 (2025)

work page 2025

[16] [16]

Maheshwari, P

B. Maheshwari, P. Stevenson and P. Van Isacker, Acta Physica Polonica B19, 1-A2 (2026)

work page 2026

[17] [17]

Algora, B

A. Algora, B. Rubio, J. Luis Tan and W. Gelletly, To- tal Absorption Technique for Nuclear Structure and Ap- plications, Springer Briefs in Physics, Springer Nature (2025)

work page 2025

[18] [18]

Nakanishi, K

Ken M. Nakanishi, K. Mitarai and K. Fujii, Phys. Rev. Research1, 033602 (2019)

work page 2019

[19] [19]

Jordan and E

P. Jordan and E. Wigner, Zeitschrift für Physik47 (9), 631 (1928)

work page 1928

[20] [20]

Wong, Introduction to Classical and Quantum Com- puting, Rooted Grove, Omaha (2022)

T. Wong, Introduction to Classical and Quantum Com- puting, Rooted Grove, Omaha (2022)

work page 2022

[21] [21]

github.io/qiskit-algorithms/stubs/qiskit_ algorithms.optimizers.L_BFGS_B.html

L-BFGS-B algorithm,https://qiskit-community. github.io/qiskit-algorithms/stubs/qiskit_ algorithms.optimizers.L_BFGS_B.html

work page

[22] [22]

Honma, T

M. Honma, T. Otsuka, T. Mizusaki, and M. Hjorth- Jensen, Phys. Rev. C80, 064323 (2009)

work page 2009

[23] [23]

De-Shalit and I

A. De-Shalit and I. Talmi, Nuclear Shell Theory, Aca- demic Press, New York (1963). 7

work page 1963