Modern parallel computing algorithm has been applied to the solution of the few-body problem. The approach is based on Feynman’s continual integrals method implemented in C++ programming language using NVIDIA CUDA technology. A wide range of 3-body and 4-body bound systems has been considered including nuclei described as consisting of protons and neutrons (e.g., 3,4He) and nuclei described as consisting of clusters and nucleons (e.g., 6He). The correctness of the results was checked by the comparison with the exactly solvable 4-body oscillatory system and experimental data.
Modern parallel computing algorithm has been applied to the solution of the few-body problem. The approach is based on Feynman’s continual integrals method implemented in C++ programming language using NVIDIA CUDA technology. A wide range of 3-body and 4-body bound systems has been considered including nuclei described as consisting of protons and neutrons (e.g., 3,4He) and nuclei described as consisting of clusters and nucleons (e.g., 6He). The correctness of the results was checked by the comparison with the exactly solvable 4-body oscillatory system and experimental data.
Liquids and solids, atoms and molecules, nuclei—all these clearly are interacting many-body systems. Even a nucleon may be regarded as a manyparticle system, not just because it is now known to consist of three quarks interacting via gluons, but because of the possibility in quantum field theory of virtual excitation of many particles from the vacuum. in the table page 56 we list some of the many-particle systems we encounter in the physical world at length scales that range from a few centimeters to a few fermis. Many of these systems exhibit phenomena—superconductivity in solids and fission in nuclei, for example—whose understanding does not follow immediately from knowledge of the constituents of the system and the interactions among the constituents, but requires new concepts and ideas. Many-body physics is the branch of theoretical physics that studies the new phenomena or “emergent properties” that arise from interactions among “elementary” constituents of a many-particle system and provides means and methods for carrying out precise calculations of such characteristic properties of these systems as may be compared with experimental results to verify hypotheses about the nature of the constituents and their interactions.
We discuss computational aspects of the spherical coupled-cluster method specific to the nuclear many-body problem. Using chiral nucleon-nucleon interaction at next-to-next-to-next-to leading order (N3LO) with cutoff Λ = 500MeV, we present coupled-cluster results for the ground state of 40Ca. Scaling and performance studies are presented together with challenges we meet with when extending the coupled-cluster effort to nuclei mass hundred and beyond. The low-energy nuclear many-body problem is a challenging undertaking, how-ever, in the last decade there has been significant progress in first principle cal-culations of nuclei.1), 2), 3), 4), 5), 6), 7), 8), 9) With the advance of chiral effective field theory10), 11), 12), 13), 14), 15) which allows for a systematic and consistent derivation of
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.
We develop a computationally and numerically efficient approach to determine binding energies and corresponding wave functions of a quantum-mechanical three-body problem in low dimensions. Our approach exploits the tensor structure intrinsic to the multidimensional stationary Schr\odinger equation, which we express as a discretized eigenvalue problem. In order to obtain numerical solutions of the three-body system in one spatial dimension, we represent the Hamiltonian operator as a series of dense matrix-matrix products and propose an efficient preconditioned Jacobi-Davidson QR iteration for the resulting algebraic eigenvalue problem. This implementation allows a significantly faster computation of three-body bound states with a higher accuracy than in previous works. To investigate the feasibility of solving higher dimensional problems, we also consider the two-dimensional case where we make use of an efficient shared-memory parallel implementation. Our novel approach is of high relevance for investigating the universal behavior of few-body problems governed only by the particle masses, overall symmetries, and the dimensionality of space. Moreover, our results have straightforward applications in the physics of ultracold atomic gases that are nowadays widely utilized as quantum sensors.
Abstract Simulating quantum systems in a finite volume is a powerful theoretical tool to extract information about them. Real-world properties of the system are encoded in how its discrete energy levels change with the size of the volume. This approach is relevant not only for nuclear physics, where lattice methods for few- and many-nucleon states complement phenomenological shell-model descriptions and ab initio calculations of atomic nuclei based on harmonic oscillator expansions, but also for other fields such as simulations of cold atomic systems. This contribution presents recent progress concerning finite-volume simulations of few-body systems. In particular, it discusses details regarding the efficient numerical implementation of separable interactions and it presents eigenvector continuation as a method for performing robust and efficient volume extrapolations.
Atomic nucleus consists of nucleons (protons and neutrons) interacting each other with the so-called nucleon-nucleon (NN) forces. The structure of atomic nucleus can be studied via many-body theories. Some of them are the well-known random-phase approximation (RPA) and/or the Bardeen-Cooper-Schrieffer (BCS) theory. These theories produce very well the properties of infinite systems such as superconductors, where the quantal and thermal fluctuations are zero or very small. However, our recent studies have shown that these fluctuations are significant in finite systems such as atomic nuclei. As the result, the physical properties of nuclear systems are changed due to the effects of quantal and thermal fluctuations. Therefore, the conventional many body theories need to be modified to include the effects of these fluctuations when they are applied to atomic nucleus. Recently, we have proposed a self-consistent quasiparticle RPA (SCQRPA) taking into account the effects of quantal and thermal fluctuations. The SCQRPA has been tested and the results obtained are in reasonable agreements with the exact solutions of the model cases as well as those obtained within the quantum Monte-Carlo method for realistic nuclei.