We report a quantum simulation of the deuteron binding energy on quantum processors accessed via cloud servers. We use a Hamiltonian from pionless effective field theory at leading order. We design a low-depth version of the unitary coupled-cluster ansatz, use the variational quantum eigensolver algorithm, and compute the binding energy to within a few percent. Our work is the first step towards scalable nuclear structure computations on a quantum processor via the cloud, and it sheds light on how to map scientific computing applications onto nascent quantum devices.
We report a quantum simulation of the deuteron binding energy on quantum processors accessed via cloud servers. We use a Hamiltonian from pionless effective field theory at leading order. We design a low-depth version of the unitary coupled-cluster ansatz, use the variational quantum eigensolver algorithm, and compute the binding energy to within a few percent. Our work is the first step towards scalable nuclear structure computations on a quantum processor via the cloud, and it sheds light on how to map scientific computing applications onto nascent quantum devices.
Accurate calculations of nuclear binding energies in theory is a long-standing problem in nuclear physics, as it provides critical insights the effective interactions and structures of atomic nuclei. Currently, there are severe limitations in theory due to the exponential growth of the Hilbert space with increasing nucleon number, i.e., curse of the dimension. There are some effort to mitigate this problem by using the state-of-art methods, such as the machine learning and quantum computing. Recent advances in quantum computing have opened new avenues for addressing these challenges, enabling the direct simulation of nuclear many-body systems on quantum hardware. In this work, a detailed study on the calculation of nuclear binding energies using the full quantum eigensolver (FQE) algorithm, implemented on a real quantum processor, i.e. Quafu, was presented. Two nuclear systems (2H and 6He) were investigated. The deuteron, as the simplest bound nucleus with just one proton and one neutron, serves as an ideal benchmark for testing quantum algorithms, FQE. 6He, a light, two-valence neutron nucleus, was used for estimating the ability of the Quafu quantum cloud platform and its “Baihua” quantum chip. Our results demonstrate that for the deuteron, the FQE algorithm on the “Baihua” quantum processor achieves a higher level of accuracy compared to previously reported results obtained using the variational quantum eigensolver (VQE) on the IBM_QX5 and IBM_19Q platforms, as discussed in Dumitrescu’s article. For 6He, the FQE-calculated binding energy shows accuracy comparable to recent VQE-based calculations reported in Yoshida’s article. These findings highlight the feasibility of performing precise nuclear structure calculations on near-term quantum devices and underscore the potential of the FQE approach for future applications in nuclear physics. Overall, this study extends the scope of quantum computing in nuclear physics by demonstrating that the FQE algorithm can serve as a practical and accurate tool for computing binding energies on real quantum hardware. The results provide a benchmark for future developments in quantum simulations of nuclei and suggest promising directions for scaling these methods to larger nuclear systems, ultimately contributing to a deeper understanding of nuclear structure and dynamics.
Calculations of nuclear properties and nuclear models are performed using quantum computing algorithms on simulated and real quantum computers. The models are a realistic calculation of deuteron binding based on effective field theory, and a simplified two-level version of the nuclear shell model. A method of reducing the number of qubits needed for practical calculation is presented, the reduction being with respect to that used when the standard Jordan-Wigner encoding is used. Its efficacy is shown in the case of the deuteron binding and shell model. A version of the variational quantum eigensolver in which all eigenstates in a spectrum are targeted on an equal basis is shown. The method involves finding the minima of the variance of the Hamiltonian, and its ability to find the full spectrum of small version of the simplified shell model is presented.
We solve the nuclear two-body and three-body bound states via quantum simulations of pionless effective field theory on a lattice in position space. While the employed lattice remains small, the usage of local Hamiltonians including two- and three-body forces ensures that the number of Pauli terms scales linearly with increasing numbers of lattice sites. We use an adaptive ansatz grown from unitary coupled cluster theory to parametrize the ground states of the deuteron and 3He, compute their corresponding energies, and analyze the scaling of the required computational resources. Our quantum simulations reproduce exact benchmarks for 2H and 3He within 100 keV, requiring at most 30 layers in the ansatz and thus resulting in modest circuit depths. Additionally, we find the number of shots required to reach a given precision scales linearly in the lattice size and more mildly in the system size. Based on the agreement with exact benchmarks and mild scaling, we conclude that this can be an efficient, scalable approach for quantum computations of nuclear ground states, particularly to prepare initial states for quantum phase estimation or other filtering algorithms.
We report results for simulating an effective field theory to compute the binding energy of the deuteron nucleus using a hybrid algorithm on a trapped-ion quantum computer. Two increasingly complex unitary coupled-cluster ansatze have been used to compute the binding energy to within a few percent for successively more complex Hamiltonians. By increasing the complexity of the Hamiltonian, allowing more terms in the effective field theory expansion, and calculating their expectation values, we present a benchmark for quantum computers based on their ability to scalably calculate the effective field theory with increasing accuracy. Our result of ${E}_{4}=\ensuremath{-}2.220\ifmmode\pm\else\textpm\fi{}0.179$ MeV may be compared with the exact deuteron ground-state energy $\ensuremath{-}2.224$ MeV. We also demonstrate an error mitigation technique using Richardson extrapolation on ion traps. The error mitigation circuit represents a record for deepest quantum circuit on a trapped-ion quantum computer.
We present the quantum simulation of the deuteron to calculate its binding energy using the variational quantum eigensolver, which is based on a hybrid quantum-classical approach. Apart from a commonly used Hamiltonian derived from pionless effective field theory, we consider the interaction which can be easily studied with conventional classical methods but the corresponding operator leading to nonzero off-tridiagonal matrix elements. To map the many-body basis states on the qubit states, three encodings are explored, namely, one-hot, Bravyi-Kitaev, and the Gray code. We perform a generalized transformation for many-body operators in Gray code encoding, and simulate the Hamiltonians with nonzero off-tridiagonal matrix elements. Furthermore, the analyses of the relative efficiency of all encodings and corresponding transformations are carried out using the noise model of a real IBM quantum device. We demonstrate that the Gray code is more efficient for a large basis, irrespective of the form of potential and the presence of hardware noise.
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.