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.
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.
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.
Quantum computing is increasingly offering concrete solutions toward the simulation of nuclear structure, with the potential to overcome the exponential scaling that limits classical diagonalization methods in large spaces. A particularly efficient encoding scheme, based on collective like-nucleon pairing modes, reduces the qubit requirements by half and avoids the non-local operator strings of standard fermion-to-qubit mappings. While this quasiparticle framework provides accurate results for semimagic nuclei, it does not adequately describe open-shell systems where proton-neutron correlations become important. In this work, we apply Brillouin-Wigner perturbation theory to systematically improve the quasiparticle description of open-shell nuclei in the sd shell, reaching an energy relative error below 0.2% compared to the nuclear shell model. Furthermore, to make the effective Hamiltonian suitable for quantum simulation, we introduce a mean-field Hartree-Fock approximation of the non-quasiparticle resolvent, achieving ground-state energies typically within 2% of the exact shell-model result. This represents a systematic improvement over the bare quasiparticle Hamiltonian while remaining within the reach of near-term quantum devices.
Quantum computing is emerging as a promising tool in nuclear physics. However, the cost of encoding fermionic operators hampers the application of algorithms in current noisy quantum devices. In this work, we analyze an encoding scheme based on pairing nucleon modes. This approach significantly reduces the complexity of the encoding, while maintaining a high accuracy for the ground states of semimagic nuclei across the sd and pf shells and for tin isotopes. In addition, we also explore the encoding ability to describe open-shell nuclei within the above configuration spaces. When this scheme is applied to a trotterized quantum adiabatic evolution, our results demonstrate a computational advantage of up to three orders of magnitude in CNOT gate count compared to the standard Jordan-Wigner encoding. Our approach paves the way for efficient quantum simulations of nuclear structure using quantum annealing, with applications to both digital and hybrid quantum computing platforms.
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.
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