Universität Bonn & Forschungszentrum Jülich (Germany)
Local Project ID:
HPC Platform used:
JUQUEEN of JSC
How do neutrons and protons bind to form atomic nuclei? Why do we observe alpha-particle clustering in light and medium-mass nuclei but not in heavy ones? These questions can be tackled in the framework of nuclear lattice effective field theory. These investigations have revealed some intriguing features of nuclei related to much discussed quantum phenomena such as entanglement and quantum phase transitions.
How do protons and neutrons bind to form nuclei? This is the central question of ab initio nuclear structure theory. While the answer may seem as simple as the fact that nuclear forces are attractive, the full story is more complex and interesting. Using the approach of nuclear lattice effective field theory, we have presented numerical evidence showing that nature is near a quantum phase transition, a zero-temperature transition driven by quantum fluctuations .
Before discussing this phase transition, let us briefly discuss the underlying computational tool. Nuclear lattice effective field theory combines the modern approach to the nuclear force problem based on a chiral effective field theory with high-performance computing methods. It defines a completely new method to exactly solve the nuclear A-body problem (with A the number of nucleons, that is protons and neutrons, in a nucleus). The first ingredient of this method is a systematic and precise effective field theory description of the forces between two and three nucleons, that has been worked out in the last decades by various groups worldwide. To go beyond atomic number four, one has to devise a method to exactly solve the A-body problem. Such a method is given by nuclear lattice simulations. Space-time is discretized with spatial length Ls and temporal length Lt , and nucleons are placed on the lattice sites. The minimal length on the lattice, the so-called lattice spacing a, entails a maximum momentum, pmax = π/a. Such a lattice representation is ideally suited for parallel computing. Breakthrough calculations have been the ab initio calculation of the Hoyle state in 12C  and of low-energy alpha-alpha scattering .
Let us come back to the quantum phase transition. Initially, we consider two types of two-nucleon interactions with a different degree of non-locality (velocity-dependence). Both interactions give an identical description of neutron-proton scattering and the properties of the lightest nucleus, the deuteron. However, when one computes the groundstate energies of alpha-like nuclei (8Be, 12C, 16O, 20Ne), a surprise emerges. The groundstate for the purely non-local interaction VA in each case turns out to be a weakly-interacting Bose gas of alpha-particles, whereas the binding energies for the more local interaction VB are within a few percent of the empirical values. In order to discuss the many-body limit, we switch off the Coulomb interaction and define a one-parameter family of interactions, Vλ = (1 − λ)VA + λVB, 0 ≤ λ ≤ 1. While the properties of the two, three, and four nucleon systems vary only slightly with λ, the many-body ground state of Vλ undergoes a quantum phase transition from a Bose-condensed gas to a nuclear liquid, as sketched in the phase diagram in Fig. 1. The phase transition occurs when the alpha-alpha s-wave scattering length aαα crosses zero, and the Bose gas collapses due to the attractive interactions. At slightly larger λ, finite alpha-like nuclei also become bound, starting with the largest nuclei first. The last alpha-like nucleus to be bound is 8Be at the so-called unitarity point where |aαα | = ∞. Superimposed on the phase diagram, we have sketched the alpha-like nuclear ground state energies EA for A nucleons up to A = 20 relative to the corresponding multi-alpha threshold Eα A/4.
Another intriguing phenomenon in nuclear physics is the observation of the formation of alpha- clusters in light and medium-mass nuclei and the disappearance of it in heavier systems. The alpha-particle is formed as a substructure because of its relatively strong binding and its resistance against external perturbations. While clustering is observed in many approaches such as density functional theory  or fermion molecular dynamics , a truly quantitative measure of clustering was not available until the work of Ref. .
Measuring the distribution of three or four nucleons on the lattices sites allows one to test the amount of clustering. Indeed, in case of non-interacting clusters, this measure called ρ4 would lead to integer numbers, like e.g. four for oxygen and its neutron-rich even-even isotopes. However, it was found that with increasing atomic number, the value calculated for ρ4 was increasing steadily when going from systems with one alpha-particle (the He isotopes) to the oxygen chain (four alpha-particles). This increase can be interpreted in terms of quantum entanglement, that is the nucleons feel the presence of the other nucleons and this eventually leads to the disappearance of the individual clusters. Stated differently, these results show that the transition from cluster-like states in light systems to nuclear liquid-like states in heavier systems should not be viewed as a simple suppression of multi-nucleon short-distance correlations, but rather an increasing entanglement of the nucleons involved in the multi-nucleon correlations, as depicted schematically in Fig. 2.
This project was made possible through the HPC system JUQUEEN of the Jülich Supercomputer Centre (JSC) under project number hfz02.
I thank my collaborators Serdar Elhatisari, Dean Lee, Gautam Rupak, Evgeny Epelbaum, Hermann Krebs, Timo A. Lähde, Ning Li and Bing-Nan Lu as part of the NLEFT collaboration.
 S. Elhatisari, N. Li, A. Rokash, J. M. Alarcon, D. Du, N. Klein, B.-N. Lu, U.-G. Meißner, E. Epelbaum, H. Krebs, T. A. L¨ahde, D. Lee and G. Rupak, Phys. Rev. Lett. 117 (2016) no.13, 132501 [arXiv:1602.04539 [nucl-th]].
 E. Epelbaum, H. Krebs, D. Lee and U.-G. Meißner, Phys. Rev. Lett. 106 (2011) 192501 [arXiv:1101.2547 [nucl-th]].
 S. Elhatisari, D. Lee, G. Rupak, E. Epelbaum, H. Krebs, T. A. Lähde, T. Luu and U.-G. Meißner, Nature 528 (2015) 111 [arXiv:1506.03513 [nucl-th]].
 J.-P. Ebran, E. Khan, T. Niksic and D. Vretenar, Nature 487 (2012) 341.
 M. Chernykh, H. Feldmeier, T. Neff, P. von Neumann-Cosel and A. Richter, Phys. Rev. Lett.98 (2007) 032501.
 S. Elhatisari, E. Epelbaum, H. Krebs, T. A. L¨ahde, D. Lee, N. Li, B.-N. Lu, Ulf-G. Meißner and G. Rupak, Phys. Rev. Lett. 119 (2017) no.22, 222505 [arXiv:1702.05177 [nucl-th]].
Prof. Dr. Ulf-G. Meißner
Universität Bonn und Forschungszentrum Jülich
D-53115 Bonn (Germany)
e-mail: meissner [at] hiskp.uni-bonn.de
JSC Project ID: hfz02