Recent Progress in Nuclear Lattice Simulations
- Facility for Rare Isotope Beams and Department of Physics & Astronomy, Michigan State University, East Lansing, MI, United States
We review several recent results on lattice simulations by the Nuclear Lattice Effective Field Theory Collaboration. In the first part we discuss the implementation of nuclear forces on the lattice using chiral effective field theory. The new development we highlight is the use of non-local lattice operators to achieve a simpler spin channel decomposition, in contrast with previous studies that considered only local interactions. In the second part, we present evidence that nuclear physics is close to a quantum phase transition. This development is also linked to the study of the differences between local and non-local interactions. In the final part we further explore the link between the nuclear forces and nuclear structure. We consider the simplest possible nuclear interaction which can accurately reproduce the ground state energies of neutron matter, light nuclei, and medium-mass nuclei. We discuss what these recent developments say about the emergence of nuclear structure from nuclear forces and the road ahead for nuclear lattice simulations.
In this article we review several recent advances by the Nuclear Lattice Effective Field Theory Collaboration that combine of chiral effective field theory with lattice methods, an approach that we call nuclear lattice effective field theory. See Epelbaum et al.  and Machleidt and Entem  for reviews of chiral effective field theory (EFT) and Lee  and Lähde and Meißner  for reviews of nuclear lattice EFT. We begin the article by presenting some technical developments regarding the nuclear lattice interactions. The new technology that makes this possible is the use of non-local smearing for the operator interactions. These are important for producing interactions that can give accurate reproductions of the empirical nucleon-nucleon phase shifts.
In the second part of the article we consider the features of the nuclear interactions which are responsible for the bulk binding properties of atomic nuclei. From ab initio lattice simulations we present numerical data showing that nature sits near a quantum phase transition. It turns out that the differences between local and non-local interactions will again be important. In the last part of the review we continue with our first principles investigations of nuclear forces and nuclear structure. This time we discuss a simple nuclear interaction which gives the correct ground state energies of neutron matter, light nuclei, and medium-mass nuclei with no more than a few percent error. We then conclude by discussing the connection between nuclear structure and nuclear forces and some possible future directions for nuclear lattice simulations.
2. Improved Lattice Interactions
Simulations using lattice chiral EFT have been applied to predict the scattering and structure of nuclei [5–7]. The application of nuclear forces at higher chiral orders are problematic on the lattice due to the breaking of rotational invariance produced by the lattice regularization [8, 9]. In this review we focus only on the short-range interactions since this is where the new developments have been made. We use the terms leading order (LO), next-to-leading order (NLO), next-to-next-to-leading order (N2LO), and next-to-next-to-next-to-leading order (N3LO) for the successive orders in chiral EFT.
In Li et al. , we introduce a set of short-range interactions in lattice chiral EFT that decompose more naturally into spin channels. In the following we explain how to construct two-nucleon operators on the lattice with intrinsic spin S, Sz, orbital angular momentum L, Lz, and isospin I, Iz. Since orbital angular momentum is not exactly conserved on the lattice, we enforce the required orbital angular momentum projection by hand using smeared distributions of annihilation and creation operators whose angular dependence are given by spherical harmonics. The use of these smeared annihilation and creation operators has the effect of controlling the range of the interaction as well as reducing lattice artifacts. Since the annihilation and creation operators are no longer at the same point in space, the resulting interaction depends on the particle velocities and is therefore a non-local interaction. This is in contrast with local interactions, where the annhiliation and creation operators are at the same point in space. Previous lattice studies had only considered local interactions .
We will work with lattice units where quantities are multiplied by the corresponding power of the spatial lattice spacing a to make the object dimensionless. Let 〈n′n〉 denote the nearest neighbors n′ of lattice site n. Our starting point is ai,j(n), the nucleon annihilation operator for lattice site n with spin i and isospin j. We then add lattice operators on neighboring lattice sites with coefficient, sNL. This defines the smeared annihilation operator
The overall normalization of the smeared annihilation operators is not important since any normalization factors can be absorbed into the definition of the interaction coefficients that we build with these smeared operators.
This process can easily be extended to lattice sites beyond the nearest neighbors of n. After this we construct combinations with spin S, Sz, and isospin I, Iz,
Let us define the lattice derivative ∇l along the l spatial direction as
We can also define the lattice derivative ∇1/2,l with half steps in the forward and backward directions,
This yields a well-defined function on the lattice sites when we take double derivatives,
The lattice Laplacian operator is defined as
We select orbital angular momenta with the help of solid spherical harmonics,
We work with polynomials of the lattice derivatives that act on one of the annihilation operators, with coefficients prescribed by the solid spherical harmonics,
This ensures the lattice operators have the correct rotational properties in the continuum limit. We use Clebsch-Gordan coefficients to put together the required spin and orbital angular momentum combinations,
In Li et al. , we consider the neutron-proton system up to N3LO. The chosen lattice spacings are 1.97, 1.64, 1.32, and 0.99 fm. In Figure S1 we show the phase shifts and mixing angles for neutron-proton scattering vs. relative momenta for lattice spacing a = 1.32 fm. In Figure S2 we show the phase shifts and mixing angles for neutron-proton scattering for lattice spacing a = 0.99 fm. The estimated uncertainties at NLO, N2LO, and N3LO are labeled with the blue, green, and red bands, respectively. The Nijmegen partial wave analysis results are shown with black solid lines, and the lattice results at N3LO are shown with diamonds. These results show marked improvement over previous studies of the lattice chiral EFT interactions . We see that in most cases the error bands at LO, NLO, N2LO, and N3LO are overlapping with each other and decreasing in width with each successive order. This indicates the consistency of the effective field theory expansion as well as its convergence rate. There are however some exceptions such as the high-momentum region of the 3P2 and 3D2 phase shifts where the convergence is not yet optimal. We have work in progress now that indicates these problems are resolved with a better lattice treatment of the tensor force in the one-pion exchange potential and will be discussed in a forthcoming publication. As one would expect, we also see a general improvement of the lattice phase shifts as we decrease the lattice spacing. We see both better agreement with the empirical data as well as smaller error bands.
3. Nuclear Physics Near a Quantum Phase Transition
In Elhatisari et al. , we investigate two different LO interactions in chiral EFT, which we call interactions A and B. Both are defined with lattice spacing a = 1.97 fm and have nearly the same nucleon-nucleon scattering phase shifts and binding energies for the three- and four-nucleon systems. However, for heavier nuclei the binding energies are very different. Interactions A and B each have the same one-pion exchange interaction as well as the same Coulomb potential between protons. The only distinction between the interactions A and B is the form of the short-range interactions, and we limit our discussion here to the short-range interactions.
We write σS with S = 1, 2, 3 for the spin Pauli matrices, and τI with I = 1, 2, 3 for the isospin Pauli matrices. We will write a(n) for the column vector of annihilation operators ai,j(n), and we will write a†(n) for the row vector of creation operators . For our chosen non-local smearing parameter sNL, we construct the smeared operators as in Equation (1),
The point-like density operators are defined as
The smeared non-local densities are defined as
while the smeared local densities for local smearing parameter sL are
The non-local short-range interactions have the form
while the local short-range interactions are given by
The :: symbol indicates normal ordering where the annihilation operators are pushed to the right and the creation operators are pulled to the left.
For interaction A we take the short-range interaction to have the purely non-local form described in Equation (24). For interaction B we take the short-range interaction to have a combination of the non-local and local forms described in Equations (24), (25), respectively. In Table 1, we show results for the ground state energies of 8Be, 12C, 16O, and 20Ne. We present results for interactions A and B at LO, including Coulomb interactions, and the comparison with experimental data. While the ground state energies for B are close to the experimental data, the binding energies for A are not.
Table 1. Results for the ground energies of 8Be, 12C, 16O, and 20Ne using interactions A and B at LO with Coulomb interactions and the comparison with experimental data.
To identify the reason for the disagreement, we compute the ground state energies of 8Be, 12C, 16O, and 20Ne for interaction A with the Coulomb interactions turned off. When we turn off Coulomb and divide the ground state energy for each nucleus with that of 4He, we find the ratios
These consecutive integer numbers are an indication of a Bose condensate of alpha particles filling the periodic box. In Elhatisari et al.  this interpretation is confirmed by showing that the alpha-alpha interactions are indeed very weak for interaction A.
These results shows that nature sits near a quantum phase transition. Along the N = Z line where the number of neutrons equals the number of protons, there is a first-order quantum transition separating a Bose gas of alpha particles from a nuclear liquid. The strength of the alpha-alpha interactions determines whether one has a Bose gas of alpha particles or a nuclear liquid. In turn, the alpha-alpha interactions are impacted by the strength and range of the local part of the nucleon-nucleon interactions. The local part of interaction B is stronger than the local part of interaction A, and hence the different behaviors that we see.
The phase diagram of nuclear matter at zero temperature is in Figure S3. The parameter λ indicates the strength of the local part of the nucleon-nucleon interactions. This can be constructed explicitly by taking a linear combination of the interactions used for interactions A and B. Nuclear matter is in a Bose gas phase of alpha particles when λ is below the critical value, and nuclear matter is in a nuclear liquid phase if λ crosses above the critical value. When λ increases even further, finite A-body nuclei become stable as their energy fall below the threshold value EαA/4 associated with multi-alpha continuum states.
4. Essential Elements for Nuclear Binding
In Lu et al. , we investigate a basic question connecting nuclear forces and nuclear structure. What are the essential elements for nuclear binding? For this analysis we work with a simple leading order effective field theory without pions. This simple pionless EFT theory is SU(4)-invariant, where the SU(4) is Wigner's approximate symmetry where the four nucleon degrees of freedom can be rotated into each other . The lattice Hamiltonian has the form
where Hfree is the free nucleon Hamiltonian and
and ρNL(n) was defined in Equation (16). The local part of the interaction is adjusted using the parameter sL, while the non-local part of the interaction is controlled by the parameter sNL, which appears implicitly in the definition of ρNL(n). The parameter C2 controls the strength of the two-nucleon interaction, and C3 controls the three-nucleon interaction. For these calculations the lattice spacing is a = 1.32 fm.
The two-nucleon interaction coefficient C2 and interaction range, controlled jointly by sNL and sL, are set by fitting the scattering length a0 and effective range r0 averaged over the two S-wave channels, and . The three-nucleon coupling strength C3 is set according to the binding energy of 3H. At the empirical binding energy B(3H) = 8.48 MeV, the 4He binding energy is 28.9 MeV, and this is near the empirical value B(4He) = 28.3 MeV. This fitting process is carried out for several pairs of values for sNL and sL. We calculate several nuclear ground states for each pair using auxiliary-field lattice Monte Carlo simulations. We find that sNL = 0.5 and sL = 0.061 gives the best agreement overall. We note that the rather large value of sL as compared with sNL is an indication that the local interaction plays an important role in nuclear binding. The corresponding couplings are MeV−2 and MeV−5. Overall about 20% of the binding is coming from the three-body interaction. This is consistent with the expected hierarchy of forces in effective field theory, with three-body forces being less important than two-body forces.
We present binding energies and charge radii in Table 2 for selected nuclei together with the experimental values and the computed Coulomb energy. Although the 3H energy is a constraint in the fitting process, the other values are predictions. We have included the charge radius of the proton for calculations of the nuclear charge radii, but have not included smaller effects arising from the charge distribution of the neutron, relativistic corrections, and spin-orbit terms. In Figure S4 we show the binding energies for 86 bound nuclei with up to A = 48 nucleons in comparison with empirical data. The Monte Carlo error bars are not visible as they are smaller than the symbol size. The errors associated with Euclidean time extrapolation and volume extrapolation are less than 1% relative error, and these errors are not shown. In Figure S4 one can see that the agreement with empirical results is fairly good. The remaining discrepancies are a sign of missing effects such as interactions which are spin dependent.
In Figure S5 we show the charge densities of 16O and 40Ca. These densities are calculated using the pinhole algorithm  where a barrier is placed in the middle of the Euclidean time evolution, and the amplitude vanishes unless each nucleon passes through a pinhole. As the number of pinholes are set to equal the number of nucleons, the sampling over pinholes yields a classical distribution of the nucleon positions. We show the comparison with electron scattering data, as well as the lattice results with the Coulomb interaction included via first order perturbation theory. The Coulomb force suppresses the densities near the nucleus center, bringing the results in better agreement with the experimental data. The results are surprisingly accurate in view of the simple nuclear interaction.
The ground state energy of pure neutron matter is shown in Figure S6 as a function of the neutron density. We show the results of several other calculations for comparison. For the lattice results the number of neutrons ranges from 14 to 66. The calculations are done with three box sizes L = 5, 6, 7. Our lattice results, shown as the filled red polygons, are generally in good agreement with the other calculations for densities higher than 0.05 fm−3. For lower densities the disagreement is larger due to the neutron-neutron scattering length being incorrect. The open red polygons correspond to an improved calculation where a contact interaction is included to give the correct neutron-neutron scattering length. There is also a correction included to restore invariance with respect to Galilean boosts . In spite of the simplicity of the interaction, the results are quite good.
We have reviewed several recent results by the Nuclear Lattice EFT Collaboration. We presented the improved description of scattering phase shifts in chiral effective field theory up to N3LO. We then showed evidence that nuclear physics is close to a quantum phase transition. After this we described the minimal nuclear interaction that can accurately reproduce the ground state energies of neutron matter, light nuclei, and medium-mass nuclei. A common theme flowing through all aspects of our review was the notable difference between local and non-local interactions. We now put some of these findings into context.
Numerous calculations show the reliability of chiral EFT for the properties of light nuclei [18–22]. However, the binding energies as well as the charge radii of medium mass nuclei are sometimes not well reproduced [20, 23–29]. One notable case is that the charge radius of 16O, which is often too small [23, 25–28]. Without further input, chiral EFT calculations do not yet provide accurate predictions at higher nuclear densities.
One practical approach is to put constraints on the nuclear force by fitting the properties of medium mass nuclei and nuclear matter saturation . This strategy has been used in several calculations [30–33]. The approach we have pursued in nuclear lattice studies is to focus on the microscopic origins of the problem. In Elhatisari et al. , we find that nuclear matter resides near a quantum phase transition that lies at the boundary between a Bose condensate of alpha particles and a nuclear liquid. There we show that local interactions are especially important for nuclear binding. In Lu et al. , we construct a simple nuclear interaction that can produce, with no more than a few percent error, the ground state energies of neutron matter, light nuclei, and medium-mass nuclei. From that analysis we see the importance of the range and locality of the SU(4)-invariant forces in determining the bulk properties of nuclei.
The dominance of SU(4)-invariant interactions can be explained in terms of coherent enhancement. Upon summation over nucleonic spin configurations, much of the effect of spin-dependent forces cancel. In a similar manner, isospin-dependent forces will cancel in symmetric nuclear matter due to the protons and neutrons being equal in number. One important exception though is the Coulomb interaction. SU(4) symmetry and large scattering length universality have a long history in nuclear physics. It is well-known that the Tjon line connecting 3H and 4He binding energies is a manifestation of universality in nuclear systems [9, 34–36].
SU(4)-symmetric short-range interactions are now being used with local and non-local smearing and one-pion exchange. These improved calculations of light and medium-mass nuclei will use chiral forces up to N3LO. One of the central questions that we seek to address is why the straightforward application of chiral EFT does not give reliable and accurate predictions at higher nuclear densities. While more investigations are needed, it seems that part of the answer to this question is related to the emergence of a physical length scale relevant to many-body nuclear systems, the size of the alpha particle. The amount of fine tuning needed to reproduce the alpha-alpha and alpha-nucleon scattering phase shifts can be understood as the competition between Pauli repulsion and attractive nucleon-nucleon interactions. This important and delicate balance is amplified by the fact that the range of the nucleon-nucleon interaction is comparable to the size of the alpha particle. Hence small differences in the range or locality of the interaction can have a large impact on the interactions of the alpha particle. This is turn has consequences for nuclear systems with increasing numbers of nucleons.
Data Availability Statement
The datasets analyzed in this article are not publicly available. Requests to access the datasets should be directed to firstname.lastname@example.org.
The author confirms being the sole contributor of this work and has approved it for publication.
Conflict of Interest
The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
I thank my wonderful collaborators Serdar Elhatisari, Evgeny Epelbaum, Hermann Krebs, Timo A. Lähde, Ning Li, Bing-Nan Lu, Thomas Luu, Ulf-G. Meißner, Gautam Rupak, for their insight, work, and determination that made on the collaborative work reported here possible. Partial support provided by the U.S. Department of Energy (DE-SC0018638 and DE-AC52-06NA25396). The computational resources were provided by the Julich Supercomputing Centre at Forschungszentrum Jülich, Oak Ridge Leadership Computing Facility, RWTH Aachen, and Michigan State University.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphy.2020.00174/full#supplementary-material
Figure S1. Results for the neutron-proton scattering phase shifts and mixing angles vs. the relative momenta for lattice spacing a = 1.32 fm. The blue, green and red bands show the estimated uncertainties at NLO, N2LO and N3LO respectively. The black solid line and diamonds denote phase shift or mixing angle from the Nijmegen partial wave analysis and lattice calculation at N3LO, respectively. Figure from Li et al. .
Figure S2. Results for the neutron-proton scattering phase shifts and mixing angles vs. the relative momenta for lattice spacing a = 0.99 fm. The blue, green and red bands show the estimated uncertainties at NLO, N2LO and N3LO respectively. The black solid line and diamonds denote phase shift or mixing angle from the Nijmegen partial wave analysis and lattice calculation at N3LO, respectively. Figure from Li et al. .
Figure S3. Phase diagram of nuclear matter at zero temperature. λ controls the strength of the local part of the interactions. Figure from Elhatisari et al. .
Figure S4. Calculated binding energies from 3H to 48Ca. The solid symbols are lattice results and the open symbols are experimental values. The experimental values are from Wang et al. . Figure adapted from Lu et al. .
Figure S5. Computed 16O and 40Ca charge densities compared with experimental results. The circles denote the results without Coulomb interaction. The squares are the results with the Coulomb interaction included perturbatively. Experimental results are from De Vries et al. . Figure adapted from Lu et al. .
Figure S6. Ground state energy of pure neutron matter as a function of neutron density for box sizes L = 5 (upright triangles), L = 6 (squares), L = 7 (right-pointing triangles). The filled red polygons are results for the leading-order interaction. The open red polygons show an improved calculation with a short-range interaction tuned to the physical neutron-neutron scattering length as well as a correction to restore Galilean invariance. Also shown are results calculated with N3LO chiral interactions (EM 500 MeV, EGM 450/500 MeV and EGM 450/700 MeV) , as well as results from variational (APR)  and auxiliary-field diffusion MC calculations (GCR) . Figure adapted from Lu et al. .
8. Alarcón JM., Du D, Klein N, Lähde TA, Lee D, Li N, et al. Neutron-proton scattering at next-to-next-to-leading order in Nuclear Lattice Effective Field Theory. Eur Phys J A. (2017) 53:83. doi: 10.1140/epja/i2017-12273-x
10. Li N, Elhatisari S, Epelbaum E, Lee D, Lu BN, Meißner UG. Neutron-proton scattering with lattice chiral effective field theory at next-to-next-to-next-to-leading order. Phys Rev C. (2018) 98:044002. doi: 10.1103/PhysRevC.98.044002
15. De Vries H, De Jager CW, De Vries C. Nuclear charge-density-distribution parameters from elastic electron scattering. Atom Data Nuclear Data Tabl. (1987) 36:495–536. doi: 10.1016/0092-640X(87)90013-1
16. Elhatisari S, Epelbaum E, Krebs H, Lähde TA, Lee D, Li N, et al. Ab initio calculations of the isotopic dependence of nuclear clustering. Phys Rev Lett. (2017) 119:222505. doi: 10.1103/PhysRevLett.119.222505
18. Hupin G, Quaglioni S, Navrátil P. Unified description of 6Li structure and deuterium-4He dynamics with chiral two- and three-nucleon forces. Phys Rev Lett. (2015) 114:212502. doi: 10.1103/PhysRevLett.114.212502
20. Lonardoni D, Carlson J, Gandolfi S, Lynn JE, Schmidt KE, Schwenk A, et al. Properties of nuclei up to A = 16 using local chiral interactions. Phys Rev Lett. (2018) 120:122502. doi: 10.1103/PhysRevLett.120.122502
22. Epelbaum E, Golak J, Hebeler K, Hüther T, Kamada H, Krebs H, Maris P, et al. Few- and many-nucleon systems with semilocal coordinate-space regularized chiral two- and three-body forces. Phys Rev C. (2019) 99:024313. doi: 10.1103/PhysRevC.99.024313
25. Hagen G, Hjorth-Jensen M, Jansen GR, Machleidt R, Papenbrock T. Continuum effects and three-nucleon forces in neutron-rich oxygen isotopes. Phys Rev Lett. (2012) 108:242501. doi: 10.1103/PhysRevLett.108.242501
29. Ekström A, Jansen JR, Wendt KA, Hagen G, Papenbrock T, Carlsson BD, et al. Accurate nuclear radii and binding energies from a chiral interaction. Phys Rev C. (2015) 91:051301. doi: 10.1103/PhysRevC.91.051301
37. Tews I, Krüger T, Hebeler K, Schwenk A. Neutron matter at next-to-next-to-next-to-leading order in chiral effective field theory. Phys Rev Lett. (2013) 110:032504. doi: 10.1103/PhysRevLett.110.032504
Keywords: nuclear structure, nuclear forces, lattice effective field theory, nuclear lattice simulations, quantum phase transition
Citation: Lee D (2020) Recent Progress in Nuclear Lattice Simulations. Front. Phys. 8:174. doi: 10.3389/fphy.2020.00174
Received: 26 December 2019; Accepted: 24 April 2020;
Published: 21 May 2020.
Edited by:Saori Pastore, Los Alamos National Laboratory (DOE), United States
Reviewed by:Robert B. Wiringa, Argonne National Laboratory (DOE), United States
Roelof Bijker, National Autonomous University of Mexico, Mexico
Copyright © 2020 Lee. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Dean Lee, email@example.com