arxiv: v3 [nucl-th] 18 Oct 2018

Transcription

1 Neutron-proton scattering with lattice chiral effective field theory at next-to-next-to-next-to-leading order Ning Li,, Serdar Elhatisari,, Evgeny Epelbaum, Dean Lee,, Bing-Nan Lu,,, 7 and Ulf-G. Meißner Facility for Rare Isotope Beams and Department of Physics and Astronomy, Michigan State University, MI, USA Helmholtz-Institut für Strahlen- und Kernphysik and Bethe Center for Theoretical Physics, Universität Bonn, D- Bonn, Germany Department of Physics, Karamanoglu Mehmetbey University, Karaman 7, Turkey Institut für Theoretische Physik II, Ruhr-Universität Bochum, D-7 Bochum, Germany Department of Physics, North Carolina State University, Raleigh, NC 79, USA Institute for Advanced Simulation, Institut für Kernphysik, and Jülich Center for Hadron Physics, Forschungszentrum Jülich, D- Jülich, Germany 7 JARA - High Performance Computing, Forschungszentrum Jülich, D- Jülich, Germany arxiv:.799v nucl-th] Oct We present a new lattice formulation of chiral effective field theory interactions with a simpler decomposition into spin channels. With these interactions the process of fitting to the empirical scattering phase shifts is simplified, and the resulting lattice phase shifts are more accurate than in previous studies. We present results for the neutron-proton system up to next-to-next-to-next-to-leading order for lattice spacings of.97,.,., and.99 fm. Our results provide a pathway to ab initio lattice calculations of nuclear structure, reactions, and thermodynamics with accurate and systematic control over the chiral nucleon-nucleon force. lini@frib.msu.edu

2 I. INTRODUCTION Chiral effective field theory (EFT) organizes the interactions of nucleons in powers of momenta and factors of the pion mass near the chiral limit where the light quarks are massless. We label terms that carry a total of n powers of nucleon momenta or factors of the pion masses as order Q n. The most important interactions at low energy are at order Q, or leading order (LO). Next-to-leading order (NLO) interactions correspond to order Q, next-to-next-to-leading order (N LO) terms ato Q, and nextto-next-to-next-to-leading order (N LO) to Q. See Ref. ] for a review of chiral EFT. Nuclear lattice simulations using chiral EFT have been used in recent years to describe the structure and scattering of atomic nuclei 7]. However, the treatment of nuclear forces at higher orders in the chiral EFT expansion is more difficult on the lattice due to the breaking of rotational invariance produced by nonzero lattice spacing, 9]. Fitting the unknown coefficients of the short-range lattice interactions to empirical phase shifts can introduce significant uncertainties. In this paper we solve these problems by introducing a new set of short-range chiral EFT interactions on the lattice with a simpler decomposition into spin channels. The angular dependence of the relative separation between the two nucleons is prescribed by spherical harmonics, and the dependence on the nucleon spins is given by the spin-orbit Clebsch-Gordan coefficients. The full details of this process are presented in this paper. We start with some definitions of the lattice operator notations used. Next we discuss the lattice Hamiltonian used in our lattice transfer matrix formalism. The short-range interactions are presented first, and we then proceed to the long-range interactions. We then compare our neutron-proton scattering results at lattice spacings of.97,.,., and.99 fm, with the empirical phase shifts. After this we compute some observable properties of the deuteron, discuss theoretical uncertainties, and present a summary and outlook. Certain interactions such as the Coulomb interaction and some isospin-breaking interactions are not directly relevant to the neutron-proton analysis that we consider here. However we include these interactions in this work for completeness and future reference. II. LATTICE OPERATOR DEFINITIONS Let us define a i,j (n) and a i,j the lattice annihilation and creation operators on lattice site n with spin i =, (up, down) and isospin j =, (proton, neutron). The operators a snl i,j (n) and asnl i,j (n) are defined via nonlocal smearing with the real parameter s NL, a snl i,j (n) = a i,j(n) + s NL a i,j (n + n ). () a snl i,j (n) = a i,j (n) + s NL n = n = Next we define the pair annihilation operators a(n)a(n )] snl S,S z,i,i z, where a i,j (n + n ). () a(n)a(n )] snl S,S z,i,i z = a snl i,j (n)m ii (S, S z)m jj (I, I z )a snl i,j (n ) () i,j,i,j with M ii (, ) = δ i, δ i, δ i, δ i,], () M ii (, ) = δ i, δ i,, () M ii (, ) = δ i, δ i, + δ i, δ i,], () We define the lattice finite difference operation l on a general lattice function f(n) as M ii (, ) = δ i, δ i,. (7) l f(n) = f(n +ˆl) f(n ˆl), ()

3 where ˆl is the spatial lattice unit vector in the l direction. It is also convenient to define the lattice finite difference operation /,l defined on points halfway between lattice sites. We use this operation only to define the Laplacian operator, /,l f(n) = f(n + ˆl) f(n ˆl). (9) / = l /,l. () Let us define the solid harmonics and their complex conjugates π R L,Lz (r) = L + rl Y L,Lz (θ, φ), () π RL,L z (r) = L + rl YL,L z (θ, φ). () Using the pair annihilation operators, lattice finite differences, and solid harmonics, we define the operator P M,sNL S,S z,l,l z,i,i z (n) = a(n) M / R L,L z ( )a(n)] snl S,S z,i,i z, () where M / and act on the second annihilation operator. More explicitly stated, this means that we act on the n in Eq. () and then set n to equal n. The even integer M gives us higher powers of the finite differences. Writing the Clebsch-Gordan coefficients as SS z, LL z JJ z, we define O M,s NL S,L,J,J z,i,i z (n) = SS z, LL z JJ z P M,sNL S,S z,l,l z,i,i z (n). () S z,l z We also define point-like density operators that depend on spin and isospin. For spin indices S =,,, and isospin indices I =,,, we define ρ(n) = i,j a i,j (n)a i,j () ρ S (n) = a i,j (n)σ S] ii a i,j () i,i,j ρ I (n) = a i,j (n)τ I] jj a i,j (7) i,j,j ρ S,I (n) = a i,j (n)σ S] ii τ I ] jj a i,j () i,i,j,j where σ S are Pauli matrices in spin space and τ I are Pauli matrices in isospin space. III. LATTICE HAMILTONIAN AND TRANSFER MATRIX FORMALISM Let a be the spatial lattice spacing and a t be the temporal lattice spacing. We work in lattice units (l.u.) where all quantities are multiplied by the powers of the spatial spacing to form a dimensionless combination. The normal-ordered transfer matrix is M =: exp Hα t ] :, (9) where the :: symbols denote normal ordering where the annihilation operators are on the right and creation operators are on the left. α t = a t /a is the ratio between the temporal lattice spacing and the spacial lattice spacing. For the temporal lattice spacing, we take a t =. fm for a =.97 fm. We rescale a t for other lattice spacings so that a /a t is fixed. We partition the lattice Hamiltonian H into a free Hamiltonian, short-range interactions, and long-range interactions,

4 H = H free + V short N For the free Hamiltonian we use an O(a )-improved action of the form ], H free = 9 m n + m a (n)a(n) a (n )a(n) m n,i n n i a (n )a(n) m n,i n n i n,i + V long N. () n n i a (n )a(n). () IV. SHORT-RANGE INTERACTIONS A. Order Q At order Q we have two short-range interaction operators, namely, the S-wave spin singlet which we call V, S I z=,,,,,,,i z,,,,,i z (n) () and the S-wave spin triplet V, S,,,J (n) z,, O,,,J z,, (n). () We note that since we work with interactions that act with a specified parity and specified total intrinsic spin S, when we act on two-nucleon states with total momentum equal to zero, the total isospin I is completely constrained by the requirement of overall antisymmetry of the two nucleons. However we still specify the correct total isospin I explicitly in order to remove lattice artifacts that might otherwise appear in cases when the total momentum is not zero. Wigner s SU() symmetry ] is an approximate symmetry of the low-energy nucleon-nucleon interactions where the nucleonic spin and isospin degrees of freedom can be rotated as four components of an SU() multiplet. As in previous work 7], we treat the SU() part of the short-range interactions separately. This choice allows us to control the strength of the local part of the SU() interaction, which has been shown to be important for the binding of nucleons in nuclei ]. So at leading order we also include an SU()-invariant short-range operator V with the form where f sl is defined as V = C : n,n,n a snl i,j (n )a snl i,j (n )f sl (n n)f sl (n n ) a snl i,j (n )a snl i,j (n ) :, () i,j i,j f sl (n) = for n =, = s L for n =, = otherwise. () We repeat again that, in terms of counting powers of momentum, this SU() interaction is equivalent to the SU()-invariant interaction we get by adding together V, S (n) and V, S (n). However, the separate treatment of this interaction allows us to control the strength of the local part of the SU() interaction in systems with more than a few nucleons. For the purposes of fitting operator coefficients, we keep the coefficient C fixed and tune the coefficients of V, S (n) and V, S (n) as needed to reproduce the scattering phase shifts and mixing angles. In this work we take the smearing parameter s NL for the SU() interaction to be the same as that used in the other short-range interactions. For a =.97 fm we take C =.7 l.u., s L =.7 l.u, and s NL =. l.u.. For a =. fm we use C =. l.u., s L =.9 l.u, and s NL =. l.u.. For a =. fm, we use C =. l.u., s L =.7 l.u, and s NL =. l.u.. For a =.99 fm, we use C =. l.u., s L =. l.u, and s NL =. l.u.. In future work, however, we may consider different smearing parameters for the two cases in order to accelerate the convergence of the effective field theory expansion in many-body systems.

5 B. Order Q At order Q we have the lowest radial excitations of the S-wave spin singlet which we call V, S I z=,, and the S-wave spin triplet V, S O,s NL,,,,,I z,,,,,i z (n) + O,s NL,,,J z,, (n) O,,,J z,, (n) + At order Q there is the P interaction V, P I z=,,,,,,,i z,s NL,,,,,I z (),,,J (n) z,, O,s NL,,,J z,, (n). (7),,,J (n) z,, O,,,J z,, () the P interaction V, P the P interaction V, P and the P interaction V, P I z=,, I z=,, I z=,, J z=,...,,,,,i z,,,,,i z (9),,,J z,,i z,,,j z,,i z (),,,J z,,i z,,,j z,,i z (n). () At order Q we also have the S D mixing term V,SD,,,J z,, (n) O,,,J z,, (n) +,,,J (n) z,, O,,,J z,, (n). () C. Order Q At order Q we have the next-to-lowest radial excitations of the S-wave spin singlet V, S, I z=,, O,s NL,,,,,I z,s NL,,,,,I z () and V, S, I z=,, O,s NL,,,,,I z,,,,,i z (n) + I z=,,,,,,,i z,s NL,,,,,I z () and the next-to-lowest radial excitations of the S-wave spin triplet V, S, O,s NL,,,J (n) z,, O,s NL,,,J z,, (n) ()

6 and V, S, O,s NL,,,J z,, (n) O,,,J z,, (n) +,,,J (n) z,, O,s NL,,,J z,, (n). () If we apply the on-shell equivalence condition that the magnitude of the outgoing relative momentum equals the magnitude of the incoming relative momentum, then V, S, and V, S, are equivalent and also V, S, and V, S, are equivalent. In this work we make the choice of setting the coefficients of V, S, and V, S, to. At order Q we have the first radial excitations of the P interaction V, P O,s NL,,,J z,, (n) O,,,J z,, (n) + the P interaction V, P O,s NL,,,,,I z,,,,,i z (n) + I z=,, the P interaction V, P O,s NL,,,J z,,i z,,,j z,,i z (n) + I z=,, and the P interaction V, P O,s NL,,,J z,,i z,,,j z,,i z (n) + I z=,, J z=,... I z=,, I z=,,,,,j (n) z,, O,s NL,,,J z,, (7),,,,,I z,s NL,,,,,I z () I z=,, J z=,...,,,j z,,i z,s NL,,,J z,,i z (9),,,J z,,i z,s NL,,,J z,,i z (n). () At order Q we also have the first radial excitations of the S D mixing term V,SD, O,s NL,,,J (n) z,, O,,,J (n) + z,,,,,j (n) z,, O,s NL,,,J z,, () and V,SD,,,,J z,, (n) O,s NL,,,J z,, (n) + O,s NL,,,J (n) z,, O,,,J z,, (n). () If we apply the on-shell equivalence condition, then V,SD, and V,SD, are equivalent. In this work we make the choice of setting the coefficient of V,SD, to. At order Q we have the D interaction V, D,,,J z,,i z,,,j z,,i z () the D interaction V, D the D interaction V, D and the D interaction V, D I z=,, J z=,... J z=,... J z=,... At order Q we also have the P F mixing term V,P F,,,J (n) z,, O,,,J (n) + z,, I z=,, J z=,...,,,j (n) z,, O,,,J z,, (),,,J (n) z,, O,,,J z,, (),,,J (n) z,, O,,,J z,, (n). () I z=,, J z=,...,,,j (n) z,, O,,,J z,, (n). (7)

7 7 D. Isospin-breaking short-range interactions We also include additional isospin-breaking S contact interactions for proton-proton scattering (I z = ) and neutronneutron scattering (I z = ). These are not relevant for neutron-proton scattering, but we nevertheless discuss the interactions for completeness. We define the two isospin-breaking interactions V Iz=, S,,,,,I (n) z= O,,,,,I z= (n) () and V Iz=, S,,,,,I (n) z= O,,,,,I z= (n). (9) In terms of counting momenta, these are order Q. However, they are suppressed by the small size of the isospin-breaking coefficient. Following our previous analyses, we count this correction as order Q. We do not consider higher-order isospin breaking terms in this work, but they will be included in future studies. V. LONG-RANGE INTERACTIONS The one-pion exchange interaction V OPE has the form where f S S is defined as V OPE = g A F π f S S(n n) = L q A. One-pion exchange n,n,s,s,i : ρ S,I(n )f S S(n n)ρ S,I (n) :, () Q(q S )Q(q S ) exp iq (n n) b π (q + Mπ)] q + Mπ, () where L is the length of the cubic periodic box and each lattice momentum component q S is an integer multiplied by π/l. The function Q(q S ) is defined as Q(q S ) = sin(q S) sin(q S) + sin(q S), () which equals q S up to a correction of order q 7 S. The parameter b π is included to remove short-distance lattice artifacts in the onepion exchange interaction. In the present calculation, we set b π =. in lattice units. We have used the combination q +M π in the exponential as suggested in recent work ] as a momentum-space regulator that does not affect the long-distance behavior. At leading order we take the pion mass to be the mass of the neutral pion, M π = M π = M π,i=. B. Two-pion exchange The cutoff momentum arising from the lattice regularization is Λ latt = π/a, with a the spatial lattice spacing. For coarse lattice spacings such as a =.97 and a =. fm, the corresponding lattice cutoffs are and 77 MeV, respectively. For momenta lying below these cutoff scales, the two-pion-exchange potential (TPEP) can be expanded in powers of q /(π ), resulting in operators that are exactly the same as our short-range contact terms. We conclude that TPEP at coarse lattice spacings can be replaced by retuning the low-energy constants (LECs) for these contact terms ]. For these two coarse lattice spacings, the TPEP does not have observable effects but only changes the LECs. For the two smaller lattice spacings, a =. fm and a =.99 fm, however, higher momenta can be reached and the structure of the two-pion-exchange potential can be resolved. In these cases we include the TPEP explicitly. According to the power counting of chiral EFT, the TPEP first appears at order O(Q ) or NLO, the subleading TPEP appears at order O(Q ) or N LO, and so on ]. Similarly to what we do to the one-pion exchange potential, we also regularize the TPEP by a Gaussian form factor in momentum space,

8 F (q) = exp b π (q + Mπ) ] = exp q + Mπ ] Λ, () where Λ = / b π ]. In the present calculation, we set b π =. lattice units, which equates to Λ = MeV for a =. fm and Λ = MeV for a =.99 fm. In this work, the relativistic corrections stemming from the /m N corrections to the one-pion exchange potential and /m N correction to the TPEP at order O(Q ) are not taken into account. The TPEP up to order O(Q ) or N LO is completely local and can be written in the form V TPEP = V Q TPEP + V Q TPEP + V Q TPEP = : ρ(n)v C (n n )ρ(n ) : + : ρ I (n)w C (n n )ρ I (n ) : n,n I n,,n + : ρ S (n)v σ (n n )ρ S (n ) : + : ρ S,I (n)w σ (n n )ρ S,I (n ) : () S n,n S,I n,n + : ρ S (n)(v T ) S,S (n n )ρ S (n ) : + : ρ S,I(n)(W T ) S,S (n n )ρ S,I(n ) : n,n n,n S,S S,S,I where V C/σ, (V T ) S,S, W C/σ, and (W T ) S,S are scalar functions in the coordinate space, and V C/σ (n n ) = L exp iq (n n )] V C/σ (q)f (q), () q W C/σ (n n ) = L exp iq (n n )] W C/σ (q)f (q), () q (V T ) S,S (n n ) = L exp iq (n n )] (V T ) S,S (q)f (q)q(q S )Q(q S ), (7) q (W T ) S,S (n n ) = L exp iq (n n )] (W T ) S,S (q)f (q)q(q S )Q(q S ). () q The definitions for the functions V C/σ (q), (W T ) C/σ (q), (V T ) S,S (q), and (W T ) S,S (q) are given in Refs., ]. C. Coulomb and long-range strong isospin breaking The Coulomb interaction will not be relevant for neutron-proton scattering, but we nevertheless discuss it here for completeness. The Coulomb interaction can be written as V Coulomb = α EM : ρ(n ) + ρ I= (n )] d(n n) ρ(n) + ρ I=(n)] :, (9) n,n where d(n n) is the shortest length of n n as measured on the periodic lattice, and we define the value of d at the origin to be /. Our notation ρ I= refers to the I = isospin component of ρ I. The long-range isospin-breaking correction, due to differences in the charged and neutral pion mass in one-pion exchange, has the form where fs IB SI is defined as V IB OPE = g A f π f IB S SI(n n) = L q n,n,s,s,i : ρ S,I(n )f IB S SI(n n)ρ S,I (n) :, () Q(q S )Q(q S ) exp iq (n n) b π (q + Mπ,I )] q + Mπ,I f S S(n n), () where M π, = M π, = M + π = M π and M π, = M π. As in previous analyses we count this correction as order Q.

9 9 VI. GALILEAN INVARIANCE RESTORATION (GIR) Galilean invariance is the statement that the laws of Newtonian physics for a non-relativistic system are independent of the velocity of the center of mass. In a lattice regularized system, however, the effect of the cutoff is different in moving frames, and this leads to the breaking of Galilean invariance 9]. There is also some breaking of Galilean invariance caused by the nonlocal smearing parameter s NL that we use in the construction of our interactions. This arises from the residual dependence of the interactions on the velocity of the center of mass. Fortunately, in many cases of interest these two Galilean invariance breaking effects have the tendency to partially cancel. In order to restore Galilean invariance in the two-nucleon system, we include the two-nucleon nearest-neighbor hopping operator, where and V GIR = C GIR V GIR = C GIR n,i,j,i,j V GIR = V GIR + V GIR, () n,i,j,i,j a i,j (n)a i,j (n)a i,j (n)a i,j(n) () n a i,j (n + n )a i,j (n + n )a i,j (n)a i,j(n). () Let us write P tot as a two-body bound-state wave function with total momentum P tot. We note that is independent of P tot, and so we have P tot V GIR P tot () P tot V GIR P tot = V GIR, () where is the two-body bound-state wave function with zero total momentum. Furthermore, Therefore P tot V GIR P tot = V GIR n e iptot n = P tot V GIR P tot cos(p tot,) cos(p tot,) ] cos(p tot,) VGIR. (7) = cos(p tot,) cos(p tot,) cos(p tot,)] V GIR. () In this manner we can restore Galilean invariance up to order Q by tuning the coefficient of V GIR according to the dispersion relation of the deuteron and the S ground state at finite volume. The deuteron dispersion relation is far more useful for this purpose, however, since the S ground state is a continuum state that shows negligible Galilean invariance breaking in its dispersion relation. This is true for all continuum states, and this is why the amount of Galilean invariance breaking seen in the higher partial waves are also negligible. We will consider Galilean invariance breaking effects beyond order Q in future work. As an example of how to determine the Galilean invariance restoration operator coefficient, we show the results for lattice spacing a =.97 fm in Fig. (). The left and right panels are the deuteron dispersion relations before and after including the Galilean invariance restoration operator respectively. We see that the amount of correction is relatively small. The fitted coefficient for V GIR is found to be C GIR =.. The amount of Galilean invariance breaking is somewhat smaller than this for the smaller lattice spacings. VII. SCATTERING ON THE LATTICE In order to calculate the scattering phase shifts and mixing angles, we first construct radial wave functions through the spherical harmonics with quantum numbers (L, L z ), ], r L,Lz = r Y L,Lz (ˆr )δ r =r r, (9)

10 Lattice Continuum W/O GIR Lattice Continuum W/ GIR Energy (MeV) Energy (MeV) - P tot (MeV) - P tot (MeV) FIG.. (Color online) Dispersion relation of the deuteron. Left: Without any Galilean invariance restoration (GIR). Right: With GIR as provided by the operator V GIR with coefficient C GIR =.. where r runs over all lattice grid points having the same radial lattice distance. We group together data into a large number of radial bins so that in each bin, r δr/ < r < r + δr/, with a very small width parameter δr. Using this definition of the radial wave function, the Hamiltonian matrix over a three-dimensional lattice can be reduced to a one-dimensional radial Hamiltonian, H r,r H r,r. We follow the method described in Ref. ], which uses an auxiliary radial potential. We extract the phase shifts as well as the mixing angles from the radial wave functions in the region where the NN force and auxiliary potentials are vanishing. In this range, the wave function has the form, A L h L (kr) B Lh + L (kr), (7) where h L (kr) and h+ L (kr) are the spherical Bessel functions, k = µe, µ is the reduced mass, and E is the energy. The scattering coefficients A L and B L satisfy the relations, B L = S L A L, (7) where S L = exp (iδ L ) is the S-matrix and δ L is the phase shift. The phase shift is determined by setting δ L = ( ) i log BL. (7) A L In the case of the coupled channels with J >, both of the coupled partial waves, L = J and L = J +, satisfy Eq. (7), and the S-matrix couples the two channels together. Throughout this work we adopt the so-called Stapp parametrization of the phase shifts and mixing angles for the coupled channels ], S = cos(ɛ) exp ( iδ J J i sin ɛ exp ( iδ J J + iδj J+ ) i sin(ɛ) exp ( ) iδ J cos(ɛ) exp ( ) iδj+ J J + iδj J+ ) ]. (7) VIII. RESULTS FOR THE NEUTRON-PROTON PHASE SHIFTS Different lattice spacings introduce different lattice artifacts. We make calculations using four different lattice spacings, a =.97,.,. and.99 fm to study the lattice spacing effects. We choose these values because the corresponding lattice momentum cutoffs, Λ latt = π/a, remain below the estimated breakdown scale of chiral effective field theory and the orderby-order convergence has been demonstrated to be favorable in few-body and many-body calculations. As noted in discussion above, we do not include the TPEP for the two coarse lattice spacings, a =.97 fm and. fm. For the two smaller lattice spacings, a =. fm and.99 fm, we present results both with and without the TPEP in order to discern the effect of the TPEP. In previous lattice studies we had to contend with interactions that had an effect in all channels. With these new lattice interactions this problem is now completely solved. We need only to consider the interactions that participate in a given channel, and our labeling of the operators makes clear which channels these are. We determine the LECs by reproducing the neutronproton scattering phase shifts and mixing angles of the Nijmegen partial wave analysis (NPWA) ]. Since the NPWA provides

11 only the statistical errors, and not the systematic errors, we use the procedure described in Ref. ] to account for the systematic uncertainties. Specifically, we use ( ) X = max NPWA X, δ NijmI X δx NPWA, δ NijmII X δx NPWA, δx Reid9 δx NPWA, (7) where δx i are the phase shifts (mixing angles ) in channel X based on different NPWA potentials, while NPWA X are the statistical errors of the phase shifts (mixing angles) of the NPWA. For the coupled channel, S D, we define the χ as χ = i (δ Latt i δ NPWA i ) i + (ELatt b E Exp E Exp b ), (7) with the deuteron binding energy E Exp b channels, we define χ as =.7 MeV and corresponding error E Exp = 9 MeV. For the other χ = i (δ Latt δ NPWA i i ). (7) In our fits we choose energy ranges that are appropriate for the chiral order and lattice spacing used. Specifically, for the coarser lattice spacings, a =.97 and. fm, we take the energy range E lab MeV for the LO, NLO/N LO, and N LO fits. In those cases we use five points, E lab =,,,, and MeV, to compute the corresponding χ. For the fits with the smaller lattice spacings, a =. and.99 fm, we take the energy range E lab MeV for the LO, NLO and N LO fits, and E lab MeV for the N LO fits. Thus we determine the χ for the N LO fits using six points, E lab =,,,, and MeV. The LECs determined by the N LO fits are listed in Table I for each of the lattice spacings. In Fig. () we show the phase shifts and mixing angles versus the relative momenta calculated using the coarsest lattice spacings, a =.97 fm. We plot the results for relative momenta up to p rel = MeV. The error bars we quote in this plot and in the following plots indicate uncertainties from the fitting procedure only. A more comprehensive analysis that includes systematic errors due to the truncated chiral EFT expansion is presented later in our discussion. From the results, it is clear that with the new lattice operators the N LO calculations reproduce the NPWA phase shifts and mixing angles for most of the S, P, and D waves with a good accuracy for relative momenta less than MeV. One can also see clearly that the agreement improves with chiral order. Unfortunately, the mixing angle ɛ bends up for the relative momenta p rel at around MeV, which indicates that higher order corrections, e.g., NLO terms, or smaller lattice spacings would be needed to get the proper behavior for ɛ at higher momenta. In Fig. () we show the neutron-proton scattering phase shifts and mixing angles versus the relative momenta calculated using lattice spacing a =. fm. The results are very similar to those at a =.97 fm, but the mixing angle ɛ stays accurate up to higher momenta compared with that at a =.97 fm. The smaller errors for the channels, P and ɛ, indicate the results at a =. fm have smaller lattice artifacts than those at a =.97 fm, as one might expect. For calculations involving the two smaller lattice spacings, a =. and.99 fm, we use the full NN interactions up to chiral order O(Q ) or N LO. The results are presented in Figs. () and (), respectively. We plot the results for relative momenta up to p rel = MeV. Compared to the results using the larger lattice spacings, one can see clear improvement. Again good convergence is observed with increasing chiral order. With the full NN interactions up to order O(Q ), the calculation using a =.99 fm can describe the S, P, and D waves with a good accuracy over the whole momentum range, < p rel < MeV. To study the importance of the long-range part of the TPEP in the calculations, we also redo the same fits without the TPEP for a =. and a =.99 fm. Our results are shown in Fig. () and (7). For the calculations using a =. and a =.99 fm, the phase shifts and mixing angles without the TPEP are very similar to those with the TPEP, though the LECs are quite different. This indicates that the TPEP can be emulated by a retuning of the LECs. At the rather low scattering energies we probe, we do not see a clear improvement due to TPEP from the phase shifts and mixing angles. However, we do expect that this will change at higher scattering energies. IX. PROPERTIES OF THE DEUTERON In this section, we calculate the properties of the deuteron using the radial deuteron wave function obtained in the calculations with a =.99 fm and the full NN interactions up to order Q or N LO. At distance r beyond the range of the interaction, the radial wave function for the deuteron in the S channel behaves as, u(r) = A S e γr, (77)

12 TABLE I. Low-energy constants determined by N LO fits using lattice spacings, a =.97 fm,. fm,. fm and.99 fm. For calculations using a =. fm and.99 fm, full NN interactions are used. All LECs are given in lattice units. LECs a =.97 fm a =. fm a =. fm a =.99 fm C, S. ±..79 ±.. ±.. ±. C, S. ±.. ±.. ±.9. ±. C, S.7 ±.. ±..7 ±.7. ±. C, S.7 ±.. ±.. ±.79. ±. C,SD. ±.. ±.7.9 ±.7.9 ±. C, P.9 ±..7 ±..9 ±.. ±. C, P. ±..9 ±.7.9 ±..9 ±. C, P. ±..7 ±.7. ±.. ±.7 C, P. ±.. ±.. ±.. ±. C, S.7 ±.. ±.. ±.. ±. C, S.79 ±.7. ±..9 ±.9. ±. C,SD. ±.. ±.. ±.9.9 ±. C, P. ±..7 ±.. ±.. ±. C, P. ±.. ±.. ±.. ±.7 C, P. ±.. ±..9 ±.. ±. C, P. ±..9 ±.. ±..7 ±. C,P F.7 ±.. ±..7 ±.. ±. C, D. ±.. ±.. ±..9 ±. C, D.7 ±..9 ±.9. ±..7 ±. C, D. ±.7. ±.9.7 ±.. ±. C, D. ±.. ±.7.9 ±..7 ±. where A S is the S-wave asymptotic normalization coefficient. Here γ = m E d with E d denoting the deuteron binding energy. In the D channel, the radial wave function behaves as w(r) = ηa S + γr + ] (γr) e γr. (7) In Fig. (), we show the radial wave functions of the deuteron calculated using a =.99 fm with the full NN interaction up to chiral order O(Q ). The left panel shows the S-wave radial wave function while the right panel is for the D-wave radial wave function. In calculating the asymptotic normalization factors, we take the range < r < fm. From the plots, one can see clearly that when the neutron and proton are well separated the S and D waves behave as the asymptotic forms in Eq. (77) and (7) respectively. The numerical values for A S and η are shown in Table II. Using the radial wave functions, we can compute the root-mean-square radius of the deuteron, r d = δrr u (r) + w (r) ]] /, (79) where δr is the small separation between the radial bins we are using for the radial deuteron wave function, and the summation is over all the radial bins. In the same manner, we can also compute the quadrupole moment of the deuteron, Q d = ] δrr w(r) u(r) w(r), () In addition to the deuteron properties, we can also compute the S-wave effective range parameters at very low energies. The effective range expansion has the form p cot(δ) = a + rp + O(p ), ()

13 δ( S ) degrees] 7 δ( S ) degrees] NPWA LO NLO/NLO NLO ε degrees] - ε degrees] p rel MeV] p rel MeV] - p rel MeV] - p rel MeV] δ( P ) degrees] p rel MeV] δ( P ) degrees] - p rel MeV] δ( P ) degrees] p rel MeV] δ( P ) degrees] - p rel MeV] δ( D ) degrees] - p rel MeV] δ( D ) degrees] p rel MeV] δ( D ) degrees] p rel MeV] δ( D ) degrees] - - p rel MeV] FIG.. (Color online) Neutron-proton scattering phase shifts and mixing angles versus relative momenta. The lattice spacing a =.97 fm is used. The TPEP is not included explicitly as discussed in the text. TABLE II. Deuteron properties and S-wave parameters calculated with the full NN interaction up to chiral order O(Q ) using a =.99 fm. Error bars we list here indicate uncertainties from the fitting procedure only. LO NLO N LO N LO Empirical E d (MeV). ±..7 ±..7 ±..7 ±..7(9)] A s(fm / ). ±.7.77 ±..777 ±..7 ±..(9)] η. ±.. ±..7 ±.. ±..() ] Q d (fm ). ±.. ±.. ±..97 ±..9() 7] r d (fm).9 ±..9 ±..9 ±..9 ±..97() ] a S. ±.. ±.. ±..7 ±..() 9] r S. ±..79 ±..7 ±..7 ±..79()9] a S. ±..9 ±.. ±..7 ±..7()9] r S. ±..7 ±..9 ±..7 ±..7() 9] P D(%).9... where p is the relative momenta between the neutron and the proton, while a and r are the scattering length and effective range respectively. Using these formula, we can extract the scattering length and effective range for S and S. In Table II, we present the properties of the deuteron and S-wave parameters obtained using the a =.99 fm and the full NN interactions up to order Q. In order to provide some insight into the nature of the lattice wave functions, we also list the computed D-wave probabilities of the deuteron, P D. We note, however, that P D is strongly dependent on short-distance physics and the scale at which it is regulated ]. From the results in Table II, it is clear that the deuteron properties can be reproduced accurately at lattice spacing a =.99 fm.

14 δ( S ) degrees] 7 δ( S ) degrees] NPWA LO NLO/NLO NLO ε degrees] - ε degrees] p rel MeV] p rel MeV] - p rel MeV] - p rel MeV] δ( P ) degrees] p rel MeV] δ( P ) degrees] - p rel MeV] δ( P ) degrees] p rel MeV] δ( P ) degrees] - p rel MeV] δ( D ) degrees] - p rel MeV] δ( D ) degrees] p rel MeV] δ( D ) degrees] p rel MeV] δ( D ) degrees] - - p rel MeV] FIG.. (Color online) Neutron-proton scattering phase shifts and mixing angles versus relative momenta. The lattice spacing a =. fm is used. The TPEP is not included explicitly as discussed in the text. There are, however, still some small systematic discrepancies that suggest additional corrections are needed. While these could be due to corrections beyond N LO in the lattice Hamiltonian, they could also be due to missing corrections to the observables themselves such as the r operator. Such corrections are needed to cancel ambiguities on how the operators are defined on a discrete lattice. For example, the nucleons could be regarded as exactly localized as delta functions at the lattice sites or they could be viewed as having some other distribution with a width comparable to the lattice spacing. While numerically small, these corrections to the operator observables are required for a full accounting of all lattice and regularization artifacts. See, for example, Ref. ]. This is an interesting but extensive subject that requires further investigation in future studies. X. THEORETICAL UNCERTAINTIES It is necessary also to address the convergence of the effective field theory expansion on the lattice and their associated systematic errors. These important topics have generated much recent interest, ]. We follow the prescription in Refs., ] where the theoretical uncertainty for some observable X(p) at order N m LO and momentum p is given by ( X NmLO (p) = max Q m+ X LO (p), Q m X LO (p) X NLO (p) ),, Q X Nm LO (p) X NmLO (p). () Here Q is the estimated expansion parameter controlling the rate of convergence, Q = max (p/λ b, M π /Λ b ), () and Λ b the breakdown momentum scale. On the lattice, cubic symmetry replaces the rotational symmetry of the continuum, and the (L + )-dimensional irreducible representation of SO() decomposes into irreducible representations of the rotational octahedral group O. For example, L = corresponds to the A of O, and L = corresponds to the T of O. However L = splits into the E and T representations of O, and similar splittings occur in all of the larger L representations. As a result

15 δ( S ) degrees] 7 δ( S ) degrees] NPWA LO NLO NLO NLO ε degrees] - - ε degrees] - - p rel MeV] p rel MeV] - p rel MeV] - p rel MeV] δ( P ) degrees] p rel MeV] δ( P ) degrees] - p rel MeV] δ( P ) degrees] p rel MeV] δ( P ) degrees] - p rel MeV] δ( D ) degrees] - - p rel MeV] δ( D ) degrees] p rel MeV] δ( D ) degrees] p rel MeV] δ( D ) degrees] - - p rel MeV] FIG.. (Color online) Neutron-proton scattering phase shifts and mixing angles versus relative momenta. The lattice spacing a =. fm and the full NN interactions are used. the breaking of rotational symmetry for L is numerically larger than that for L <. This leads to a lower momentum breakdown scale for D waves and above compared to the S and P waves. To account for this in our calculations, we take Λ b to be the lattice momentum cutoff Λ latt = π/a for the lower partial waves, and we take Λ b = (/)Λ latt for ɛ, D waves, and higher partial waves. We will study the dependence of the lattice breakdown scale on L and J in more detail in future work. The theoretical error bands for the neutron-proton scattering phase shifts and mixing angles versus the relative momenta for a =.97,.,. and.99 fm are shown in Figs. (9-), in which we see a systematic decrease in the uncertainties for the S- and D- wave phase shifts with decreasing lattice spacing. The unexpectedly small NLO uncertainties for the P phase shifts at a coarse lattice spacing are caused by the rather good but accidental accuracy of the P phase shifts at LO. We also show the estimated theoretical uncertainties for the neutron-proton scattering phase shifts and mixing angles for a =. and.99 fm without the long-range TPEP in Figs. () and (). With only a few exceptions, the error bands for each order generally overlap with each other and cover the empirical phase shifts. This is a promising sign of convergence of the chiral effective field theory expansion on the lattice. XI. SUMMARY AND OUTLOOK We have proposed a new lattice formulation of the chiral NN force which is easily decomposed into partial waves. The new lattice operators work as projection operators, which only survive in particular channels. This advantage greatly simplifies the fitting procedure. Instead of fitting the phase shifts and mixing angles for all the channels simultaneously, only one uncoupled channel or two coupled channels are needed to be computed for each calculation. To study the dependence on the lattice spacing, we have computed neutron-proton phase shifts and mixing angles using four different lattice spacings, a =.97 fm,. fm,. fm and.99 fm. For two coarser lattice spacings, a =.97 fm or. fm, we did not explicitly include the TPEP, whereas for those using a =. fm or.99 fm, we did. Our numerical results

16 δ( S ) degrees] 7 p rel MeV] δ( S ) degrees] NPWA LO NLO NLO NLO p rel MeV] ε degrees] p rel MeV] ε degrees] p rel MeV] δ( P ) degrees] p rel MeV] δ( P ) degrees] - p rel MeV] δ( P ) degrees] p rel MeV] δ( P ) degrees] - p rel MeV] δ( D ) degrees] - - p rel MeV] δ( D ) degrees] p rel MeV] δ( D ) degrees] p rel MeV] δ( D ) degrees] - - p rel MeV] FIG.. (Color online) Neutron-proton scattering phase shifts and mixing angles versus relative momenta. The lattice spacing a =.99 fm and the full NN interactions are used. indicate a good convergence with respect to chiral order. One also observes an obvious improvement when the lattice spacing is decreased Comparing the results obtained with and without the TPEP at lattice spacings of a =. and.99 fm, we did not find significant differences. This may, however, be a consequence of the rather low scattering energies we probe, and we expect that differences would appear at higher scattering energies. We have also studied the properties of the deuteron wave function and the S-wave effective range parameters obtained with the full NN interaction at lattice spacing a =.99 fm. The numerical values are very close to the empirical values, which indicates that the current version of NN interactions is quite accurate, and a very significant improvement over previous lattice studies. Some small discrepancies remain, but these may well be fixed in studies that reach a higher order in the chiral effective field theory expansion. In summary, the new lattice interactions are far more efficient and accurate in reproducing physical data than previous lattice interactions. We have begun studying the properties of light and medium-mass nuclei using these interactions, and the results are promising. These interactions were specifically designed to facilitate very efficient Monte Carlo simulations of few- and many-body systems using auxiliary fields. The results of these studies using these new interactions will be reported in several future publications. XII. ACKNOWLEDGEMENT We acknowledge partial financial support from the Deutsche Forschungsgemeinschaft (SFB/TR, Symmetries and the Emergence of Structure in QCD ), the BMBF (Grant No. P), the U.S. Department of Energy (DE-SC), and the Scientific and Technological Research Council of Turkey (TUBITAK project no. F). Further support was provided by the Chinese Academy of Sciences (CAS) President s International Fellowship Initiative (PIFI) (Grant No. DM) and by VolkswagenStiftung (Grant No. 9). The computational resources were provided by the Jülich Supercomputing Centre at Forschungszentrum Jülich, Oak Ridge Leadership Computing Facility, RWTH Aachen, North Carolina State University, and

17 7 δ( S ) degrees] 7 δ( S ) degrees] NPWA LO NLO/NLO NLO ε degrees] - ε degrees] p rel MeV] p rel MeV] - p rel MeV] - p rel MeV] δ( P ) degrees] p rel MeV] δ( P ) degrees] - p rel MeV] δ( P ) degrees] p rel MeV] δ( P ) degrees] - p rel MeV] δ( D ) degrees] - p rel MeV] δ( D ) degrees] p rel MeV] δ( D ) degrees] p rel MeV] δ( D ) degrees] - - p rel MeV] FIG.. (Color online) Neutron-proton scattering phase shifts and mixing angles versus relative momenta. The lattice spacing a =. fm is used. The TPEP is not included in this case for comparison. Michigan State University. ] E. Epelbaum, H.-W. Hammer, and U.-G. Meißner, Rev. Mod. Phys., 77 (9), arxiv:. nucl-th]. ] E. Epelbaum, H. Krebs, D. Lee, and U.-G. Meißner, Phys. Rev. Lett., 9 (),.7. ] E. Epelbaum, H. Krebs, T. Lähde, D. Lee, and U.-G. Meißner, Phys. Rev. Lett. 9, (),.. ] E. Epelbaum, H. Krebs, T. A. Lähde, D. Lee, U.-G. Meißner, and G. Rupak, Phys. Rev. Lett., (),.77. ] S. Elhatisari, D. Lee, G. Rupak, E. Epelbaum, H. Krebs, T. A. Lähde, T. Luu, and U.-G. Meißner, Nature, (),.. ] S. Elhatisari et al., Phys. Rev. Lett. 7, (),.9. 7] S. Elhatisari, E. Epelbaum, H. Krebs, T. A. Lähde, D. Lee, N. Li, B.-n. Lu, U.-G. Meißner, and G. Rupak, Phys. Rev. Lett. 9, (7), ] J. M. Alarcón, D. Du, N. Klein, T. A. Lähde, D. Lee, N. Li, B.-N. Lu, T. Luu, and U.-G. Meißner, Eur. Phys. J. A, (7), ] N. Klein, S. Elhatisari, T. A. Lähde, D. Lee, and U.-G. Meißner, Eur. Phys. J. A, (),.. ] D. Lee, Prog. Part. Nucl. Phys., 7 (9),.. ] E. Wigner, Phys. Rev., (97). ] P. Reinert, H. Krebs, and E. Epelbaum, Eur. Phys. J. A, (), 7.. ] B. Borasoy, E. Epelbaum, H. Krebs, D. Lee, and U.-G. Meißner, Eur. Phys. J. A, (), ] S. Weinberg, Phys. Lett. B, (99). ] N. Kaiser, Phys. Rev. C, 7 (), nucl-th/7. ] E. Epelbaum, H. Krebs, and U.-G. Meißner, Eur. Phys. J. A, (),.. 7] D. R. Entem, N. Kaiser, R. Machleidt, and Y. Nosyk, Phys. Rev. C9, (),.. ] E. Epelbaum, W. Glöckle, and U.-G. Meißner, Nucl. Phys. A77, (), nucl-th/. 9] D. Lee and R. Thomson, Phys. Rev. C7, (7), nucl-th/7.

18 δ( S ) degrees] 7 δ( S ) degrees] NPWA LO NLO/NLO NLO ε degrees] - ε degrees] p rel MeV] p rel MeV] - p rel MeV] - p rel MeV] δ( P ) degrees] p rel MeV] δ( P ) degrees] - p rel MeV] δ( P ) degrees] p rel MeV] δ( P ) degrees] - p rel MeV] δ( D ) degrees] - p rel MeV] δ( D ) degrees] p rel MeV] δ( D ) degrees] p rel MeV] δ( D ) degrees] - - p rel MeV] FIG. 7. (Color online) Neutron-proton scattering phase shifts and mixing angles versus relative momenta. The lattice spacing a =.99 fm is used. The TPEP is not included in this case for comparison. ] B.-N. Lu, T. A. Lähde, D. Lee, and U.-G. Meißner, Phys. Lett. B7, 9 (),.. ] S. Elhatisari, D. Lee, U.-G. Meißner, and G. Rupak, Eur. Phys. J. A, 7 (),.. ] H. P. Stapp, T. J. Ypsilantis, and N. Metropolis, Phys. Rev., (97). ] V. G. J. Stoks, R. A. M. Klomp, M. C. M. Rentmeester, and J. J. de Swart, Phys. Rev. C, 79 (99). ] C. Van Der Leun and C. Alderliesten, Nucl. Phys. A, (9). ] T. E. O. Ericson and M. Rosa-Clot, Nucl. Phys. A, 97 (9). ] N. L. Rodning and L. D. Knutson, Phys. Rev. C, 9 (99). 7] D. M. Bishop and L. M. Cheung, Phys. Rev. A, (979). ] A. Huber, T. Udem, B. Gross, J. Reichert, M. Kourogi, K. Pachucki, M. Weitz, and T. W. Hansch, Phys. Rev. Lett., (99). 9] O. Dumbrajs, R. Koch, H. Pilkuhn, G. c. Oades, H. Behrens, J. j. De Swart, and P. Kroll, Nucl. Phys. B, 77 (9). ] J. L. Friar, Phys. Rev. C, (979). ] N. Klein, D. Lee, and U.-G. Meißner (), 7.. ] E. Epelbaum, H. Krebs, and U.-G. Meißner, Phys. Rev. Lett., (),.. ] R. J. Furnstahl, D. R. Phillips, and S. Wesolowski, J. Phys. G, (), 7.7. ] R. J. Furnstahl, N. Klco, D. R. Phillips, and S. Wesolowski, Phys. Rev. C9, (),..

19 9. lattice Asymptotic Normalization. lattice Asymptotic Normalization.... µ(r) fm -/ ]. ω(r) fm -/ ] r fm] r fm] FIG.. (Color online) Radial wave functions of the deuteron and its asymptotic behavior at large r. Wave functions are calculated using a =.99 fm with the full NN interaction up to order O(Q ). Left: S wave. Right: D wave. δ( S ) degrees] δ( S ) degrees] Error(NLO) Error(NLO) NPWA Lattice (NLO) ε degrees] - - ε degrees] - - p rel MeV] p rel MeV] - p rel MeV] - p rel MeV] δ( P ) degrees] p rel MeV] δ( P ) degrees] p rel MeV] δ( P ) degrees] p rel MeV] δ( P ) degrees] p rel MeV] δ( D ) degrees] - p rel MeV] δ( D ) degrees] p rel MeV] δ( D ) degrees] p rel MeV] δ( D ) degrees] - - p rel MeV] FIG. 9. (Color online) Theoretical error bands for the neutron-proton scattering phase shifts and mixing angles versus the relative momenta for a =.97 fm. Blue and red bands signify the estimated uncertainties at NLO and N LO respectively. The black solid line and diamonds denote the phase shift or mixing angle from the Nijmegen partial-wave analysis (NPWA) and lattice calculation at N LO, respectively.

20 δ( S ) degrees] δ( S ) degrees] Error(NLO) Error(NLO) NPWA Lattice (NLO) ε degrees] - - ε degrees] - - p rel MeV] p rel MeV] - p rel MeV] - p rel MeV] δ( P ) degrees] p rel MeV] δ( P ) degrees] p rel MeV] δ( P ) degrees] p rel MeV] δ( P ) degrees] p rel MeV] δ( D ) degrees] - p rel MeV] δ( D ) degrees] p rel MeV] δ( D ) degrees] p rel MeV] δ( D ) degrees] - - p rel MeV] FIG.. (Color online) Theoretical error bands for the neutron-proton scattering phase shifts and mixing angles versus the relative momenta for a =. fm. Blue and red bands signify the estimated uncertainties at NLO and N LO, respectively. The black solid line and diamonds denote phase shift or mixing angle from the Nijmegen partial-wave analysis (NPWA) and lattice calculation at N LO, respectively.

21 δ( S ) degrees] δ( S ) degrees] Error(NLO) Error(NLO) Error(NLO) NPWA Lattice (NLO) ε degrees] - - ε degrees] - - p rel MeV] p rel MeV] - p rel MeV] - p rel MeV] δ( P ) degrees] p rel MeV] δ( P ) degrees] p rel MeV] δ( P ) degrees] p rel MeV] δ( P ) degrees] p rel MeV] δ( D ) degrees] - p rel MeV] δ( D ) degrees] p rel MeV] δ( D ) degrees] p rel MeV] δ( D ) degrees] - - p rel MeV] FIG.. (Color online) Theoretical error bands for neutron-proton scattering phase shifts and mixing angles versus the relative momenta for a =. fm with the full NN interaction. Blue, green, and red bands signify the estimated uncertainties at NLO, N LO, and N LO, respectively. The black solid line and diamonds denote the phase shift or mixing angle from the Nijmegen partial-wave analysis (NPWA) and lattice calculation at N LO, respectively.

22 δ( S ) degrees] δ( S ) degrees] Error(NLO) Error(NLO) Error(NLO) NPWA Lattice (NLO) ε degrees] - - ε degrees] - - p rel MeV] p rel MeV] - p rel MeV] - p rel MeV] δ( P ) degrees] p rel MeV] δ( P ) degrees] p rel MeV] δ( P ) degrees] p rel MeV] δ( P ) degrees] p rel MeV] δ( D ) degrees] - p rel MeV] δ( D ) degrees] p rel MeV] δ( D ) degrees] p rel MeV] δ( D ) degrees] - - p rel MeV] FIG.. (Color online) Theoretical error bands for the neutron-proton scattering phase shifts and mixing angles versus the relative momenta for a =.99 fm with the full NN interaction. Blue, green, and red bands signify the estimated uncertainties at NLO, N LO and N LO respectively. The black solid line and diamonds denote the phase shift or mixing angle from the Nijmegen partial-wave analysis (NPWA) and lattice calculation at N LO, respectively.

23 δ( S ) degrees] δ( S ) degrees] Error(NLO) Error(NLO) NPWA Lattice (NLO) ε degrees] - - ε degrees] - - p rel MeV] p rel MeV] - p rel MeV] - p rel MeV] δ( P ) degrees] p rel MeV] δ( P ) degrees] p rel MeV] δ( P ) degrees] p rel MeV] δ( P ) degrees] p rel MeV] δ( D ) degrees] - p rel MeV] δ( D ) degrees] p rel MeV] δ( D ) degrees] p rel MeV] δ( D ) degrees] - - p rel MeV] FIG.. (Color online) Theoretical error bands for the neutron-proton scattering phase shifts and mixing angles versus the relative momenta for a =. fm without the TPEP. Blue and red bands signify the estimated uncertainties at NLO and N LO respectively. The black solid line and diamonds denote phase shift or mixing angle from the Nijmegen partial-wave analysis (NPWA) and lattice calculation at N LO, respectively.

24 δ( S ) degrees] δ( S ) degrees] Error(NLO) Error(NLO) NPWA Lattice (NLO) ε degrees] - - ε degrees] - - p rel MeV] p rel MeV] - p rel MeV] - p rel MeV] δ( P ) degrees] p rel MeV] δ( P ) degrees] p rel MeV] δ( P ) degrees] p rel MeV] δ( P ) degrees] p rel MeV] δ( D ) degrees] - p rel MeV] δ( D ) degrees] p rel MeV] δ( D ) degrees] p rel MeV] δ( D ) degrees] - - p rel MeV] FIG.. (Color online) Theoretical error bands for the neutron-proton scattering phase shifts and mixing angles versus the relative momenta for a =.99 fm without the TPEP. Blue and red bands signify the estimated uncertainties at NLO and N LO respectively. The black solid line and diamonds denote phase shift or mixing angle from the Nijmegen partial-wave analysis (NPWA) and lattice calculation at N LO, respectively.