The nuclear shell model is one of the prime many-body methods to study the structure of atomic nuclei, but it is hampered by an exponential scaling on the basis size as the number of particles increases. We present a shell-model quantum circuit design strategy to find nuclear ground states by exploiting an adaptive variational quantum eigensolver algorithm. Our circuit implementation is in excellent agreement with classical shell-model simulations for a dozen of light and medium-mass nuclei, including neon and calcium isotopes. We quantify the circuit depth, width and number of gates to encode realistic shell-model wavefunctions. Our strategy also addresses explicitly energy measurements and the required number of circuits to perform them. Our simulated circuits approach the benchmark results exponentially with a polynomial scaling in quantum resources for each nucleus. This work paves the way for quantum computing shell-model studies across the nuclear chart and our quantum resource quantification may be used in configuration-interaction calculations of other fermionic systems.
The nuclear shell model is one of the prime many-body methods to study the structure of atomic nuclei, but it is hampered by an exponential scaling on the basis size as the number of particles increases. We present a shell-model quantum circuit design strategy to find nuclear ground states by exploiting an adaptive variational quantum eigensolver algorithm. Our circuit implementation is in excellent agreement with classical shell-model simulations for a dozen of light and medium-mass nuclei, including neon and calcium isotopes. We quantify the circuit depth, width and number of gates to encode realistic shell-model wavefunctions. Our strategy also addresses explicitly energy measurements and the required number of circuits to perform them. Our simulated circuits approach the benchmark results exponentially with a polynomial scaling in quantum resources for each nucleus. This work paves the way for quantum computing shell-model studies across the nuclear chart and our quantum resource quantification may be used in configuration-interaction calculations of other fermionic systems.
The nuclear shell model accurately describes the structure and dynamics of atomic nuclei. However, the exponential scaling of the basis size with the number of degrees of freedom hampers a direct numerical solution for heavy nuclei. In this work, we present a quantum annealing protocol to obtain nuclear ground states. We propose a tailored driver Hamiltonian that preserves a large gap and validate our approach in a dozen nuclei with basis sizes up to 10^5 10 5 using classical simulations of the annealing evolution. We explore the relation between the spectral gap and the total time of the annealing protocol, assessing its accuracy by comparing the fidelity and energy relative error to classical benchmarks. While the nuclear Hamiltonian is non-local and thus challenging to implement in current setups, the estimated computational cost of our annealing protocol on quantum circuits is polynomial in the single-particle basis size, paving the way to study heavier nuclei.
Abstract In this work, we introduce a new qubit mapping strategy for the Variational Quantum Eigensolver (VQE) applied to nuclear shell model calculations, where each Slater determinant (SD) is mapped to a qubit, rather than assigning qubits to individual single-particle states. While this approach may increase the total number of qubits required in some cases, it enables the construction of simpler quantum circuits that are more compatible with current noisy intermediate-scale quantum (NISQ) devices. We apply this method to seven nuclei: Four lithium isotopes 6−9Li from the \textit{p}-shell, 18F from the \textit{sd}-shell, and two heavier nuclei (210Po, and 210Pb). We run circuits representing their ground states on a noisy simulator (IBM's \textit{FakeFez} backend) and quantum hardware (ibm_pittsburgh). For heavier nuclei, we demonstrate the feasibility of simulating 210Po and 210Pb as 22- and 29-qubit systems, respectively. Additionally, we employ Zero-Noise Extrapolation (ZNE) via two-qubit gate folding to mitigate errors in both simulated and hardware-executed results. Post-mitigation, the best results show less than 4 \% deviation from shell model predictions across all nuclei studied. This SD-based qubit mapping proves particularly effective for lighter nuclei and two-nucleon systems, offering a promising route for near-term quantum simulations in nuclear physics. PACS numbers: 21.60.Cs, 21.30.Fe, 21.10.Dr, 27.20.+n
In this work, we compare three qubit-mapping strategies to study the structure of the nuclear ground state within the shell model description employing the Variational Quantum Eigensolver (VQE) approach. Although the initial point for different mappings is a Hamiltonian matrix in many-body particle basis or Slater determinant (SD) basis, the structure of the trial wavefunction and resource counts are different for each mapping. These three mappings are tested for a mid p-shell nucleus 10B and compared the quantum resources required to find the ground state for each mapping. Further, we extend the qubit-efficient mapping to study the ground state of one more mid p-shell nucleus 12C. We run circuits up to 26-qubits representing their ground states on a noisy simulator (IBM's FakeFez backend) and quantum hardware (ibm_fez). The best post-error mitigated results from the hardware for 10B ground state is obtained following SD to qubit mapping with a percent error of 0.21 \%. The percent errors for the same state following cSD and pnSD mapping are 3.37 and 8.88 \%, respectively. On the other hand, following the cSD mapping, the post-error mitigated ground state energy of 12C is 6.82 \% away from the exact result. We further evaluate the fidelity of the VQE wavefunctions obtained from hardware with respect to the shell model wavefunctions for the cSD mapping. This cSD mapping can be useful for scaling the VQE algorithm for complex nuclei across different mass regions in terms of qubit efficiency.
Abstract In this work, we introduce a new qubit mapping strategy for the Variational Quantum Eigensolver (VQE) applied to nuclear shell model calculations, where each Slater Determinant (SD) is mapped to a qubit, rather than assigning qubits to individual single-particle states. While this approach may increase the total number of qubits required in some cases, it enables the construction of simpler quantum circuits that are more compatible with current noisy intermediate-scale quantum (NISQ) devices. We apply this method to seven nuclei: Four lithium isotopes 6−9 Li from the p -shell, 18 F from the sd -shell, and two heavier nuclei ( 210 Po, and 210 Pb). We run circuits representing their ground states on a noisy simulator (IBM’s FakeFez backend) and quantum hardware ( ibm_pittsburgh ). For heavier nuclei, we demonstrate the feasibility of simulating 210 Po and 210 Pb as 22- and 29-qubit systems, respectively. Additionally, we employ Zero-Noise Extrapolation (ZNE) via two-qubit gate folding to mitigate errors in both simulated and hardware-executed results. Post-mitigation, the best results show less than 4 % deviation from shell model binding energy predictions across all nuclei studied. This SD-based qubit mapping proves particularly effective for lighter nuclei and two-nucleon systems, offering a promising route for near-term quantum simulations in nuclear physics.
To maximize the value of fault-tolerant quantum computers, it is essential to develop concrete applications beyond well-established domains such as chemistry and condensed-matter physics. Here we construct and compile quantum algorithms to simulate the structure of atomic nuclei -- a topic that has received relatively little attention from the quantum computing community despite its similarities to the electronic structure problem in chemistry -- via effective shell-model Hamiltonians and no-core-shell-model Hamiltonians with three-body interactions derived from chiral effective field theory. Furthermore, we provide quantum resource estimates, in terms of Toffoli gate and qubit counts, for these algorithms, which, to our knowledge, are the first such estimates for fault-tolerant quantum simulation of atomic nuclei. Notably, the estimates for 32Mg and 219At shell-model Hamiltonians are comparable to recent estimates of Femoco simulations, a standard benchmark in chemistry. For no-core-shell-model Hamiltonians suitable for light nuclei (up to 40Ca or so), we find that resource requirements are significantly higher, suggesting that more bespoke strategies are required to make such simulations practicable. Throughout this work, we draw upon the similarities between nuclear and electronic structure problems, while also highlighting challenges that are specific to the former. We hope this work will spur long-term collaborations between the nuclear and quantum computing community with the ultimate goal of realizing useful nuclear simulations on quantum computers.
Quantum computing is one of the most promising technologies of the near future, and the simulation of quantum many-body systems is a natural application. In this work, we present classical simulations of the ground states of light atomic nuclei within the p shell, from He 6 to B 10 , calculated within the nuclear shell model. We compare the performance of two leading variational quantum eigensolver algorithms: the unitary coupled cluster (UCC) and the adaptive derivative-assembled pseudo-Trotter (ADAPT) methods, introducing a new metric to quantify the use of quantum resources in each simulation. We find that Slater determinants are the most useful reference states for both approaches. Our analysis suggests that ADAPT is more efficient for nuclei close to magic numbers, while UCC tends to require fewer resources toward the mid shell. This work lays the groundwork for robust benchmarking of quantum algorithms in nuclear structure studies.