arxiv:11.8v1 [nucl-th] 8 Nov 0 Lattice effective field theory for nuclei from A = to A = 8, a Evgeny Epelbaum, b Hermann Krebs, b Dean Lee, c Ulf-G. Meißner, ade and Gautam Rupak f a Institute for Advanced Simulation, Institut für Kernphysik, and Jülich Center for Hadron Physics, Forschungszentrum Jülich, D Jülich, Germany b Institut für Theoretische Physik II, Ruhr-Universität Bochum, D 80 Bochum, Germany c Department of Physics, North Carolina State University, Raleigh, NC 6, USA d Helmholtz-Institut für Strahlen- und Kernphysik and Bethe Center for Theoretical Physics, Universität Bonn, D 11 Bonn, Germany e JARA High Performance Computing, Forschungszentrum Jülich, D Jülich, Germany f Department of Physics and Astronomy, Mississippi State University, Mississippi State, MS 6, USA E-mail: t.laehde@fz-juelich.de, evgeny.epelbaum@rub.de, hermann.krebs@rub.de, dean_lee@ncsu.edu, u.meissner@fz-juelich.de, grupak@u.washington.edu We present an overview of the extension of Nuclear Lattice Effective Field Theory simulations to the regime of medium-mass nuclei. We focus on the determination of the ground-state energies of the alpha nuclei O, 0 Ne, Mg and 8 Si by means of Euclidean time projection. PACS: 1.10.Dr, 1.0.-x, 0.De 1st International Symposium on Lattice Field Theory LATTICE 0 July August, 0 Mainz, Germany Speaker. c Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. http://pos.sissa.it/
1. Introduction Nuclear Lattice Effective Field Theory (NLEFT) is a first-principles approach, in which Chiral EFT for nucleons is combined with numerical Auxiliary-Field Quantum Monte Carlo (AFQMC) lattice simulations. NLEFT differs from other ab initio methods [1 6] in that it is an unconstrained Monte Carlo calculation, which does not rely on truncated basis expansions or many-body perturbation theory, nor on prior information about the structure of the nuclear wave function.. Nuclear Lattice EFT at NNLO As in Chiral EFT, our calculations are organized in powers of a generic soft scale Q associated with factors of momenta and the pion mass []. We denote O(Q 0 ) as leading order (LO), O(Q ) as next-to-leading order (NLO), and O(Q ) as next-to-next-to-leading order (NNLO) contributions. The present calculations are performed up to NNLO. We define H LO as the LO lattice Hamiltonian, and H SU() as the equivalent Hamiltonian with the pion-nucleon coupling g A = 0 and contact interactions that respect Wigner s SU() symmetry. In our NLEFT calculations (see Ref. [8] for a review), H LO is treated non-perturbatively. The NLO contribution to the two-nucleon force (NF), the electromagnetic and strong isospin-breaking contributions (EMIB), and the three-nucleon force (NF) which first enters at NNLO, are all treated as perturbations. It should be noted that our LO calculations use smeared short-range interactions that capture much of the corrections usually treated at NLO []. The NF at NNLO over-binds nuclei with A > due to a clustering instability which involves four nucleons on the same lattice site. The long-term objective of NLEFT is to remedy this problem by decreasing the lattice spacing and including the NLO corrections in Chiral EFT. In the mean time, the over-binding problem has been rectified by means of a N contact interaction, tuned to the empirical binding energy of either He or 8 Be [10]. While this provides a good description of the alpha nuclei up to A = 1 including the Hoyle state [10 1], the over-binding is found to increase more rapidly for A. Therefore, in Ref. [] a non-local N interaction which accounts for all possible configurations of four nucleons on adjacent lattice sites was introduced, and adjusted to the empirical binding energy of Mg. A detailed study of the spectrum of O will be reported separately [].. Euclidean time projection The NLEFT calculations reported here (see also Ref. []) are performed with a (spatial) lattice spacing of a=1. fm in a periodic cube of length L=1 fm. Our trial wave function Ψ init A is a Slater-determinant state composed of delocalized standing waves, with A nucleons and the desired spin and isospin. First, we project Ψ init A for a time t using the Euclidean-time evolution operator of the SU() Hamiltonian, giving the trial state Ψ A (t ) exp( H SU() t ) Ψ init A. Second, we use the full Hamiltonian H LO to construct the Euclidean-time projection amplitude Z A (t) Ψ A (t ) exp( H LO t) Ψ A (t ), E A (t)= [lnz A (t)]/ t, (.1) and the transient energy E A (t). If we denote by Ψ A,0 the lowest (normalizable) eigenstate of H LO which has a non-vanishing overlap with the trial state Ψ A (t ), we obtain the corresponding energy E A,0 as the t limit of E A (t).
0... E (LO) [MeV] 0 0 E (NLO) [MeV]. E (EMIB) [MeV] 1. E (NF) [MeV].. 1. 1. 0 1 11.. E / C 0 [l.u.] 11 10. 10. 8. 8. 6. - -. -8-8. - -... E / C Q 1 1 11 10 E / C S - -...... E / C S I -0 - - -6-8 -0 - - -8 E / C (QS) - -. -. -.6 -.8 - -. -. -.6 -.8 (QS) E / C S - -. -. -.6 -.8 - -. -. -.6 -.8 S E ( x π ) [MeV] 6 0 8 6 0 8 E / D (OPE) [l.u.].8.6. 1. E / D (TPEq) [l.u.] -0.08-0.0-0.1-0.11-0.1-0. -0. -0.1 E / D (TPEm) [l.u.] -0.06-0.08-0.1-0.1-0. -0. -0. E / D (TPExx) [l.u.] -0. -0. -0. -0.6-0.8 6 8-0. E / D (contact) [l.u.].... E / D (contact) [l.u.] 1. 1. 1. 1. 1.1 1 0. 0.8 0. 0.6 0. Figure 1: NLEFT results for O. The LO energy is E LO =.() MeV, and the result at NNLO including N interactions is E NNLO+N = 1.() MeV. The empirical binding energy is 1.6 MeV. The NLO and NNLO contributions are evaluated in perturbation theory. We compute operator expectation values using Z O A(t) Ψ A (t ) exp( H LO t/)o exp( H LO t/) Ψ A (t ), (.) for any operator O. Given the ratio XA O(t)=ZO A (t)/z A (t), the expectation value of O for the desired state Ψ A,0 is obtained as XA,0 O Ψ A,0 O Ψ A,0 =lim t XA O(t). Sign oscillations make it difficult to reach sufficiently large values of the projection time t. It is helpful to note that the closer the trial state Ψ A (t ) is to Ψ A,0, the less the necessary projection time t. Ψ A (t ) can be optimized by adjusting both the SU() projection time t and the strength of the coupling C SU() of H SU(). The accuracy of the extrapolation t can be further improved by simultaneously incorporating data from trial states that differ in C SU().
0 0. 1. 1 E 0 (LO) [MeV] 0 0 0 E 0 (NLO) [MeV] 1. E 0 (EMIB) [MeV]. E 0 (NF) [MeV] 0 1-0. 1. E 0 / C 0 [l.u.]. 1. 1 11. 11 10. 10. E 0 / C I...... E 0 / C Q 1 0 1 1 1 E 0 / C I E 0 / C S E 0 / C S I - -8-0 - - -6-8 0 E 0 / C (QS) -. -. -.6 -.8 - -. -. -.6 -.8. E 0 / C I (QS) E 0 / C S -. -. -.6 -.8 - -. -. -.6 -.8. E 0 / C I S E 0 ( x π ) [MeV] 60 8 6 0 8 6 E 0 / D (OPE) [l.u.].8.6...8.6. E 0 / D (TPEq) [l.u.] -0.1-0. -0. -0.1-0.1-0.1-0. -0.1 E 0 / D (TPEm) [l.u.] -0. -0. -0. -0.6-0.8 E 0 / D (TPExx) [l.u.] -0. -0. -0. -0. E 0 / D (contact) [l.u.]. 6. 6.. E 0 / D (contact) [l.u.] 1. 1. 1 0.8 Figure : NLEFT results for 0 Ne. The LO energy is E LO = 1.() MeV, and the result at NNLO including N interactions is E NNLO+N =.() MeV. The empirical binding energy is 0.6 MeV. The large-time behavior of Z A (t) and ZA O(t) is controlled by the low-energy spectrum of H LO. Let E label the eigenstates of H LO with energy E, and let ρ A (E) denote the density of states for a system of A nucleons. We then express Z A (t) and ZA O (t) in terms of their spectral representations, Z A (t)= de ρ A (E) E Ψ A (t ) exp( Et), (.) ZA O (t)= de de ρ A (E)ρ A (E ) Ψ A (t ) E E O E E Ψ A (t ) exp( (E+ E )t/), (.) from which we construct the spectral representations of E A (t) and XA O (t). We can approximate these to arbitrary accuracy over any finite range of t by taking ρ A (E) to be a sum of energy delta functions, ρ A (E) i max i=0 δ(e E A,i ), where we take i max = for the He ground state, and i max = for A 8. Using data obtained for different values of C SU(), we perform a correlated fit of E A (t)
0 0 1. 6.8 0.6 E (LO) [MeV] 0 0-0 E (NLO) [MeV] 1 0. E (EMIB) [MeV].. E (NF) [MeV] -8 0 0.8-0.6-0 1. - E / C 0 [l.u.]. 1. 1.. E / C I..... E / C Q... 1. 1 0. 0 1. 1 E / C I - - -8 E / C S.... E / C S I 0 6 8-0 - - E / C (QS)...6.8 - -. -. E / C I (QS)..... E / C S...6.8 - -. -. E / C I S..... E ( x π ) [MeV] 0 68 66 6 6 60 8 E / D (OPE) [l.u.].8.6...8.6. E / D (TPEq) [l.u.] -0.1-0. -0. -0. -0. -0. -0. -0.6-0.6 E / D (TPEm) [l.u.] 1 6 8 E / D (TPExx) [l.u.] -0. -0.6-0.8-0. -0. -0. -0.6-0.8 E / D (contact) [l.u.] 8. 8. 6. E / D (contact) [l.u.]....1 1. 1. 1. Figure : NLEFT results for Mg. The LO energy is E LO = () MeV, and the result at NNLO including N interactions is E NNLO+N = () MeV. The empirical binding energy is.6 MeV. and XA O (t) for all operators O that contribute to the NLO and NNLO energy corrections. We find that the use of 6 trial states allows for a much more precise determination of E A,0 and XA,0 O than hitherto possible. In particular, we may triangulate XA,0 O using trial states that correspond to functions XA O(t) which converge both from above and below, thereby bracketing XO A,0.. Results The NLEFT results for O are given in Fig. 1, for 0 Ne in Fig., for Mg in Fig., and for 8 Si in Fig.. The curves show a correlated fit for all trial states, using the same spectral density ρ A (E). The upper row in each figure shows the LO energy, the total isospin-symmetric NF correction (NLO), the electromagnetic and isospin-breaking corrections (EMIB) and the total
0-80.. - -6 E 8 (LO) [MeV] -00-0 E 8 (NLO) [MeV]. E 8 (EMIB) [MeV].8.6. E 8 (NF) [MeV] -8-0 -0. - -60 - E 8 / C 0 [l.u.]. 1. 1. E 8 / C I.... E 8 / C Q. 6. 6... E 8 / C I -.. - -. -8-8. - -. -0-0. E 8 / C S.... - E 8 / C S I -0 - - -6-8 -80-8 -8 E 8 / C (QS) -. -. -.6 -.8-8 -8. -8. E 8 / C I (QS).. - -.. - -. E 8 / C S -. -. -.6 -.8-8 -8. -8. E 8 / C I S.. - -.. - -. E 8 ( x π ) [MeV] 8 6 1 0 6 68 E 8 / D (OPE) [l.u.]..8..6.....1 E 8 / D (TPEq) [l.u.] -0.6-0. -0. -0.8-0.8-0. -0. 0 1 1 E 8 / D (TPEm) [l.u.] -0.1-0. -0. -0. -0.6-0. -0.8-0. -0. E 8 / D (TPExx) [l.u.] -0. -0.1-0. -0. -0. -0. -0.6-0. -0.8-0. -0.6-0.61 E 8 / D (contact) [l.u.]. 8.8 8.6 8. 8. 8.8.6.. E 8 / D (contact) [l.u.].8.6. 1. Figure : NLEFT results for 8 Si. The LO energy is E LO = 0() MeV, and the result at NNLO including N interactions is E NNLO+N = () MeV. The empirical binding energy is 6. MeV. NF correction. The remaining panels show the matrix elements XA O (t) that form part of the NLO and NF terms. The operators E A / C i give the contributions of the NLO contact interactions, and E A ( x π ) denotes the energy shift due the O(a )-improved pion-nucleon coupling. The operators E A / D i give the individual contributions to the total NF correction. To summarize, we have reported on the extension of NLEFT to the regime of medium-mass nuclei. While the NNLO results are good up to A = 1, an increasing over-binding (associated with the momentum-cutoff scale and neglected higher-order contributions) manifests itself for A. While the long-term objectives of NLEFT are to decrease the lattice spacing and include higher orders in the EFT expansion, we also find that the missing physics can be approximated by an effective N interaction. The current exploratory results represent an important step towards more comprehensive NLEFT simulations of medium-mass nuclei in the future. 6
Acknowledgments We are grateful for the help in automated data collection by Thomas Luu. Partial financial support from the Deutsche Forschungsgemeinschaft (Sino-German CRC 110), the Helmholtz Association (Contract No. VH-VI-1), BMBF (Grant No. 06BN006), and the U.S. Department of Energy (DE-FG0-0ER0) is acknowledged. This work was further supported by the EU HadronPhysics project, and funds provided by the ERC Project No. NUCLEAREFT. The computational resources were provided by the Jülich Supercomputing Centre at the Forschungszentrum Jülich and by RWTH Aachen. References [1] G. Hagen, M. Hjorth-Jensen, G. R. Jansen, R. Machleidt, and T. Papenbrock, Evolution of shell structure in neutron-rich calcium isotopes, Phys. Rev. Lett. 10, 00 (01). [] E. D. Jurgenson et al., P-shell nuclei using Similarity Renormalization Group evolved three-nucleon interactions, Phys. Rev. C 8, 01 (0). [] R. Roth, J. Langhammer, A. Calci, S. Binder, and P. Navratil, Similarity-Transformed Chiral NN+N Interactions for the Ab Initio Description of 1-C and -O, Phys. Rev. Lett. 10, 001 (011). [] H. Hergert et al., In-Medium Similarity Renormalization Group with Chiral Two- Plus Three-Nucleon Interactions, Phys. Rev. C 8, 00 (0). [] A. Lovato et al., Charge form factor and sum rules of electromagnetic response functions in Carbon, Phys. Rev. Lett. 111, 001 (0). [6] V. Somà, C. Barbieri, and T. Duguet, Ab initio Gorkov-Green s function calculations of open-shell nuclei, Phys. Rev. C 8, 010(R) (0). [] E. Epelbaum, H.-W. Hammer, and Ulf-G. Meißner, Modern Theory of Nuclear Forces, Rev. Mod. Phys. 81, 1 (00). [8] D. Lee, Lattice simulations for few- and many-body systems, Prog. Part. Nucl. Phys. 6, 11 (00). [] B. Borasoy, E. Epelbaum, H. Krebs, D. Lee, and Ulf-G. Meißner, Lattice Simulations for Light Nuclei: Chiral Effective Field Theory at Leading Order, Eur. Phys. J. A 1, 10 (00). [10] E. Epelbaum, H. Krebs, D. Lee, and Ulf-G. Meißner, Lattice effective field theory calculations for A =,,6,1 nuclei, Phys. Rev. Lett. 10, 01 (010); ibid., Lattice calculations for A =,, 6, 1 nuclei using chiral effective field theory, Eur. Phys. J. A, (010). [11] E. Epelbaum, H. Krebs, D. Lee, and Ulf.-G. Meißner, Ab initio calculation of the Hoyle state, Phys. Rev. Lett. 106, 101 (011). [1] E. Epelbaum, H. Krebs, T. A. Lähde, D. Lee, and Ulf-G. Meißner, Structure and rotations of the Hoyle state, Phys. Rev. Lett. 10, 01 (01); ibid., Viability of carbon-based life as a function of the light quark mass, Phys. Rev. Lett. 110, 110 (0); ibid., Dependence of the triple-alpha process on the fundamental constants of nature, Eur. Phys. J. A, 8 (0). [] T. A. Lähde, E. Epelbaum, H. Krebs, D. Lee, Ulf-G. Meißner, and G. Rupak, Lattice Effective Field Theory for Medium-Mass Nuclei [arxiv:11.0 [nucl-th]]. [] E. Epelbaum, H. Krebs, T. A. Lähde, D. Lee, Ulf-G. Meißner, and G. Rupak, in preparation.