Abstract Quantum computers have proven to be effective in simulating many quantum systems. Simulating nuclear processes and state preparation poses significant challenges, even for traditional supercomputers. This study demonstrates the feasibility of a complete simulation of a nuclear transition, including the preparation of both ground and first excited states. To tackle the complexity of strong interactions between two and three nucleons, the states are modeled on the tritium nucleus. Both the initial and final states are represented using quantum circuits with variational quantum algorithms and inductive biases. Describing the spin‐isospin states requires four qubits, and a parameterized quantum circuit that exploits a total of 16 parameters is initialized. The estimated energy has a relative error of 2% for the ground state and 10% for the first excited state of the system. The simulation estimates the transition probability between the two states as a function of the dipole polarization angle. This work marks a first step toward leveraging digital quantum computers to simulate nuclear physics.
Abstract Quantum computers have proven to be effective in simulating many quantum systems. Simulating nuclear processes and state preparation poses significant challenges, even for traditional supercomputers. This study demonstrates the feasibility of a complete simulation of a nuclear transition, including the preparation of both ground and first excited states. To tackle the complexity of strong interactions between two and three nucleons, the states are modeled on the tritium nucleus. Both the initial and final states are represented using quantum circuits with variational quantum algorithms and inductive biases. Describing the spin‐isospin states requires four qubits, and a parameterized quantum circuit that exploits a total of 16 parameters is initialized. The estimated energy has a relative error of 2% for the ground state and 10% for the first excited state of the system. The simulation estimates the transition probability between the two states as a function of the dipole polarization angle. This work marks a first step toward leveraging digital quantum computers to simulate nuclear physics.
Abstract Quantum computers have shown promise in simulating complex quantum systems, including nuclear processes that challenge even supercomputers. We summarize the key ingredients to demonstrate the feasibility of a full nuclear transition simulation, covering also ground and excited state preparation. The tritium nucleus has been used to model strong interactions between nucleons, with quantum circuits used to represent the initial and final states. Variational quantum algorithms aid in preparing such states, requiring four qubits to describe spin-isospin states. Our results show low relative errors in the energy estimation for the obtained eigenstates. Additionally, the transition probability between these states has been estimated as a function of dipole polarization angle. We also study the transition for the neutron capture d(n, γ)t and compare the result with the tritium de-excitation.
A central challenge in nuclear physics is understanding quantum many-body systems governed by the strong nuclear force. The inherent complexity of these systems, combined with the limitations of classical computational methods, underscores the need for new approaches to study nuclear structure and dynamics. Here, we demonstrate that a spin-based digital quantum simulator using nuclear magnetic resonance, where nuclear spins simulate interacting fermions, offers a powerful tool to address this challenge. As a first step, we experimentally simulate the Agassi model, which encapsulates the interplay between collective and single-particle behaviors in finite nuclei. By representing nucleons as both bosons (nucleon pairs) and fermions (individual unpaired nucleons), the Agassi model captures highly non-linear interactions and is particularly suited for studying nuclear phase transitions, such as those between spherical and deformed shapes. We experimentally measure the correlation function as an order parameter during the evolution of the many-body system, successfully detecting a quantum phase transition. Specifically, we observe a sharp transition between the symmetric phase and the broken symmetry phase. This work underscores the potential of quantum simulation as a transformative tool in nuclear physics, particularly for exploring complex quantum many-body systems with applications in nuclear structure and reaction dynamics.
While quantum algorithms have been studied as an efficient tool for the stationary state energy determination in the case of molecular quantum systems, no similar study for analogical problems in computational nuclear physics (computation of energy levels of nuclei from empirical nucleon-nucleon or quark-quark potentials) have been realized yet. Although the difference between the above mentioned studies might seem negligible, it will be examined. First steps towards a particular simulation (on classical computer) of the Iterative Phase Estimation Algorithm for deuterium and tritium nuclei energy level computation will be carried out with the aim to prove algorithm feasibility (and extensibility to heavier nuclei) for its possible practical realization on a real quantum computer.
Abstract The nuclear binding energy problem is one of the fundamental challenges of the nuclear many-body problem, as calculating the corresponding binding energies by various approximate methods remains computationally challenging for a wide range of nuclei. In this work, a binding energy model derived from quantum computational approaches to nuclear physics is proposed and implemented in a quantum computer simulator to address the nuclear binding energy problem. The goal is to develop a framework capable of calculating the binding energies of multiple nuclei within a unified scheme. To this end, appropriate potential expressions are selected and tuned to accommodate varying nuclear properties. The defined Hamiltonian is parameterized and transformed into a form suitable for implementation on a quantum computer. Depending on the number of nucleons, the program is run on 2 to 12 qubits to obtain the minimum eigenvalue of the Hamiltonian corresponding to the ground state energy (binding energy of the nucleus) according to the variational principle of quantum mechanics. This framework provides, the binding energies of the first 29 nuclei with a total deviation of 16.86% ( 2 H, 3 H, 3 He, 4 He, 6 He, 8 He, 6 Li, 7 Li, 8 Li, 9 Li, 11 Li, 7 Be, 8 Be, 9 Be, 10 Be, 11 Be, 12 Be, 8 B, 9 B, 10 B, 11 B, 12 B, 9 C, 10 C, 11 C, 12 C, 11 N, 12 N, 12 O) obtained by quantum computation following the steps of Variational Quantum Eigensolver (VQE). These results support the feasibility of the proposed approach for nuclear binding energy calculations. In this sense, this work contributes to quantum computing and nuclear physics by providing a model able to calculate the binding energies for the first 29 nuclei using quantum computer simulator and offers a promising route for future large-scale applications.
Background: The many-body nature of nuclear physics problems poses significant computational challenges. These challenges become even more pronounced when studying the resonance states of nuclear systems, which are governed by the non-Hermitian Hamiltonian. Quantum computing, particularly for quantum many-body systems, offers a promising alternative, especially within the constraints of current noisy intermediate-scale quantum (NISQ) devices.Purpose: This work aims to simulate nuclear resonances using quantum algorithms by developing a variational framework compatible with non-Hermitian Hamiltonians and implementing it fully on a quantum simulator.Methods: We employ the complex scaling technique to extract resonance positions classically and adapt it for quantum simulations using a two-step algorithm. First, we transform the non-Hermitian Hamiltonian into a Hermitian form by using the energy variance as a cost function within a variational framework. Second, we perform $\ensuremath{\theta}$-trajectory calculations to determine optimal resonance positions in the complex energy plane. To address resource constraints on NISQ devices, we utilize Gray code (GC) encoding to reduce qubit requirements.Results: We first validate our approach using a schematic potential model that mimics a nuclear potential, successfully reproducing known resonance energies with high fidelity. We then extend the method to a more realistic $\ensuremath{\alpha}\text{\ensuremath{-}}\ensuremath{\alpha}$ nuclear potential and compute the D- and G-wave resonance energies with a basis size of N=16, using only four qubits. The quantum simulation results closely match the classical values, demonstrating the feasibility of our approach.Conclusions: This study demonstrates, for the first time, that the complete $\ensuremath{\theta}$-trajectory method can be implemented on a quantum computer without relying on any classical input beyond the Hamiltonian. The results establish a scalable and efficient quantum framework for simulating resonance phenomena in nuclear systems. This work represents a significant step toward quantum simulations of open quantum systems and lays the foundation for future investigations into resonance structures in nuclear, atomic, and molecular physics.
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.