arxiv: v1 [hep-lat] 4 Nov 2016

Transcription

1 Computing the nucleon Dirac radius directly at Q = 0 arxiv:6.0383v [hep-lat] 4 Nov 06 ab, Michael Engelhardt d, Jeremy Green g, Stefan Krieg ab, Stefan Meinel eh, John Negele c, Andrew Pochinsky c and Sergey Syritsyn f h a Bergische Universität Wuppertal, Wuppertal, Germany b JSC, Jülich Supercomputing Centre, Forschungszentrum Jülich, Jülich, Germany c Center for Theoretical Physics, Massachusetts Institute of Technology, Cambridge, USA d Department of Physics, New Mexico State University, Las Cruces, NM 88003, USA e Department of Physics, University of Arizona, Tucson, AZ 857, USA f Department of Physics and Astronomy, Stony Brook University, Stony Brook, NY 794, USA g PRISMA Cluster of Excellence and Institut für Kernphysik, Johannes Gutenberg-Universität Mainz, Germany h RIKEN BNL Research Center, Brookhaven National Laboratory, Upton, NY 973, USA n.hasan@fz-juelich.de, engel@nmsu.edu, jeremy.green@desy.de, s.krieg@fz-juelich.de, smeinel@ .arizona.edu, negele@mit.edu, avp@mit.edu, ssyritsyn@quark.phy.bnl.gov We describe a lattice approach for directly computing momentum derivatives of nucleon matrix elements using the Rome method, which we apply to obtain the isovector magnetic moment and Dirac radius. We present preliminary results calculated at the physical pion mass using a HEX-smeared Wilson-clover action. For removing the effects of excited-state contamination, the calculations were done at three source-sink separations and the summation method was used. 34th annual International Symposium on Lattice Field Theory 4-30 July 06 University of Southampton, UK Speaker. Current affiliation: NIC, DESY, Zeuthen, Germany c Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License (CC BY-NC-ND 4.0).

2 Computing the nucleon Dirac radius directly at Q = 0. Introduction The proton radius puzzle refers to the 7σ discrepancy in the experimental determinations of the charge radius of the proton between the value extracted from spectroscopy of muonic hydrogen [], r p E = (39) fm, and the CODATA value determined from scattering and spectroscopy with electrons [], r p E = (5) fm. The way to determine the proton radius rp E from scattering data is to measure the slope of the electric form factor G p E (Q ) as a function of the squaredmomentum transfer Q as Q approaches zero. There is controversy surrounding finding the radius from fitting the scattering data: some claim to find a small radius from scattering data [3, 4, 5] and others argue that only a large radius is compatible with the data [6, 7, 8]. On the lattice, the same approach is usually followed. However, the smallest nonzero Q reached on the lattice with standard methods and a large volume of (5.6 fm) 3 is about 0.05 GeV, whereas the smallest Q reached in scattering experiments is below GeV. Therefore, obtaining a reliable radius from fitting to lattice form factor data may be challenging. This motivates the need for a direct calculation of the radius without fitting to form factors. For the case of a pion, it was shown in [7] that the Rome method for momentum derivatives could be used to calculate the pion charge radius with finite-volume effects that are exponentially suppressed, with asymptotic behaviour m π L e mπl. In the following, we present a method for calculating the nucleon Dirac radius directly at Q = 0. Our approach is based on calculating the momentum derivatives of the two- and threepoint functions using the Rome method [9], as explained in section. A ratio of the two- and three-point functions is then constructed and by calculating the momentum derivatives of this ratio at Q = 0 we are able to extract the anomalous magnetic moment κ = F (0) and Dirac radius r = 6 df Q F (sections 3 and 4). Our results for both quantities are shown in section 5. dq =0. Momentum derivatives of the two-point and three-point functions For calculating the momentum derivatives of the correlation functions, we need to consider quark propagators with smeared- and point-sources and sinks. This is shown in the following: Without smearing: Let N denote a proton or neutron field. The two-point function can be written as: C (p,x 0,z 0 ) αβ = e ip(x z) N α (x) N β (z) x = e ip(x z) ε abc ε de f f αγδε f βζηθ G a f γθ (x,z)gbe x = ε abc ε de f f αγδε f βζηθ G a f γθ (x,z)gbe x δη (x,z)gcd δη (x,z)gcd εζ εζ (x,z) Gac γη(x,z)g b f δθ (x,z)gcd εζ (x,z) (x,z;p) Gac γη(x,z)g b f δθ (x,z)gcd εζ (x,z;p), where f αβγδ is the spin tensor determining the quantum numbers of the nucleon operator N, G(x,z) is the quark propagator and G(x,z;p) = e ip(x z) G(x,z). As shown in [9], the first and second (.)

3 Computing the nucleon Dirac radius directly at Q = 0 momentum derivatives of a quark propagator at zero momentum are given by: where p j G(x,y;p) p=0 = i G(x,z)Γ j V G(z,y), (.) z ( p j ) G(x,y;p) p=0 = G(x,z)Γ j V G(z,z )Γ j V G(z,y) z,z z G(x,z)Γ j T G(z,y), (.3) Γ j V /T G(z,y) U + γ j j (z ĵ) G(z ĵ,y) U j (z) γ j G(z + ĵ,y). (.4) For connected diagrams, the three-point function, with current O Γ = qγq and zero sink momentum p = 0, can be written as: C 3 (p,x 0,y 0,z 0 ) αβ = e ip(y z) N α (x)o Γ (y) N β (z) = G S (y)γg(y,z;p), (.5) x,y y where G S (y) is the sequential backward propagator, which is independent of p. Only the forward propagator G(y, z; p) needs to be expanded using. and.3. Hence, no additional backward propagators are needed. The Rome method can be understood as doing a calculation with twisted boundary conditions and then taking the derivative with respect to the twist angle at zero twist [7]. With smearing: In the two-point function, we have the smeared-source smeared-sink propagator G(x,y;p) = e ip(x y) x,y K(x,x )G(x,y )K(y,y) = e ip(x x ) K(x,x ) e ip(x y ) G(x,y ) e ip(y y) K(y,y), x,y }{{}}{{}}{{} K(x,x ;p) G(x,y ;p) K(y,y;p) where K is the smearing kernel. The momentum derivatives can then be calculated using the product rule along with. and.3. Denoting the momentum derivative with for shorter notation, (KGK) = K GK + K(GK), (.6) (KGK) = K GK + K (GK) + K(GK). (.7) For the smeared-source point-sink propagator, which is needed for the three-point function, we get: (GK) = G[ iγ V GK + K ], (.8) (GK) = G[ iγ V (GK) Γ T GK + K ]. (.9) Gaussian Wuppertal smearing is given by K(x,y;p) = x,x,... K 0 (x,x ;p)k 0 (x,x ;p)...k(x...,y;p), with }{{} N W ( K 0 (x,y;p) = e ip(x y) 3 [ ] ) δ x,y + α + 6α U j (x)δ x+ĵ,y +U j (x ĵ)δ x ĵ,y j= = + 6α ( δ x,y + 3 [ α e ip j U j (x)δ x+ĵ,y + e ip j U j x ĵ,y] ) (x ĵ)δ. (.0) j= The mth derivative of K 0 at zero momentum is equal to [ ] ( K (m) ) mk0 0 (x,y) p j (x,y;p) = α i m U j (x)δ x+ĵ,y + ( i) m U j + 6α (x ĵ)δ x ĵ,y. (.) p=0

4 Computing the nucleon Dirac radius directly at Q = 0 The first and second derivatives of smearing with N W iterations, K = K N W 0, can be computed iteratively using (K0 N) = K 0 KN 0 + K 0 (K0 N ), (K0 N) = K 0 KN 0 + K 0 (KN 0 ) + K 0 (K0 N ). 3. Ground-state contribution to correlation functions We will be tracing the correlators with polarization matrices that contain the projector ( + γ 0 )/, so that we can effectively write the overlap matrix elements as 0 N α (0) N(p,s) = Z(p)u α (p,s), [5, 6]. Here and in the following we use Minkowski-space gamma matrices. The Dirac and Pauli form factors, F (N,q) (Q ) and F (N,q) (Q ), parametrize matrix elements of the vector current N(p,s ) qγ µ q N(p,s) = ū(p,s )F [γ µ,p,p]u(p,s), (3.) with the short-hand notation F [γ µ,p,p] = F (N,q) γ µ + F (N,q) iσ µν (p p) ν m, where Q = (p p). Having T = x 0 z 0 and τ = y 0 z 0, the ground-state contributions to the two- and three-point functions (for p = 0) are: C (p,t ) αβ = Z (p) E(p) e E(p)T (m + /p) αβ, (3.) C 3 (p,τ,t ) αβ = Z(0) m Z(p) E(p) e m(t τ) e E(p)τ m [ F (γ µ,0,p)(m + /p) ] αβ. (3.3) 4. Momentum derivatives of the ratio Because we don t know how Z(p) depends on the momentum, we need to compute derivatives of the ratio of three-point and two-point functions. We set p = 0 and p = ke j, where e j is the unit vector in j-direction. We compute the following ratio: R N (k,τ,t ) αβ = R αβ (k,τ,t ) = R N (k,τ,t ) αβ R A (k,τ,t ), with (4.) C 3 (k,τ,t ) αβ C (p = 0,T )C (k,t ), R A(k,τ,T ) = C (k,t τ)c (p = 0,τ) C (p = 0,T τ)c (k,τ). (4.) For computing the first and second momentum derivatives of this ratio we need: R N(k) αβ = C (k)c 3(k) αβ + C (k)c 3 (k) αβ C (0)C (k) 3, (4.3) R N(k) αβ = (3[C (k)] C (k)c (k))c 3(k) αβ + 4C (k)( C (k)c 3 (k) αβ +C (k)c 3 (k) αβ ) 4 C (0)C (k) 5. (4.4) where, for more readability we suppress τ,t parameters, use the notation C (k) = Tr ( +γ 0 C (k) ) and denote the derivatives with a prime e.g. C (k) dc (k) dk. We know that C (0) = 0 in the infinitestatistics limit. Hence, we can eliminate this from the ratios. Similarly, we can calculate R A (k) and R A (k) which can be used together with 4.3 and 4.4 to calculate the first and second derivatives of the ratio R αβ. These derivatives are computed on the lattice directly at k = 0 as discussed earlier in 3

5 Computing the nucleon Dirac radius directly at Q = 0 section. The ground-state contributions are equal to: R αβ (k) = [F (k)(m + Eγ0 kγ j )] αβ 8 E(E + m) R αβ (k) = [ F (k)(m + Eγ 0 kγ j ) ] αβ + [ F (k)(e γ 0 γ j ) ] αβ 8 E(E + m), (4.5) [ F (k)(m + Eγ 0 kγ j ) ] (E + m)e αβ 6, [E(E + m)] 3/ and R αβ (k) can be calculated in a similar way. We use the continuum dispersion relation E(k) = m + k, which implies Q = m m + k m, and find that at k = 0, the second derivative is needed to obtain the slope of F : df dk (4.6) = dq df Q k=0 dk k=0 dq = 0, d F k=0 =0 dk = df Q dq. (4.7) =0 Furthermore, we have at k = 0, E(0) = m, E (0) = 0, E (0) = /m and, F (0) = F (0)γ µ, F (0) = F (0) iσ µ j m, F (0) = d dq F (0)γ µ iσ µ0 F (0) m. (4.8) Using Γ pol = ( + γ 3 γ 5 ) +γ0, we find nonzero results for the following values of the index µ (labeling the components of the vector current), Tr[R(µ = 0)Γ pol ] = F, Tr[ R(µ = )Γ pol ] = i m (F + F ), (4.9) Tr[ R(µ = )Γ pol ] = i m (F + F ), Tr[,,3R(µ = 0)Γ pol ] = 4m (F + F ) 3 F r, (4.0) with i = and r p i = 6 df Q F. From equations 4.9 and 4.0 we find the following relations dq =0 for the anomalous magnetic moment κ and Dirac radius r : κ = m Im(Tr[R (µ = ) Γ pol ]) Tr[R(µ = 0) Γ pol ], (4.) r = mim[r (µ = ) Γ pol ] + 3Tr[R(µ = 0) Γ pol ] m Tr[R (µ = 0) Γ pol ] 4m, (4.) Tr[R(µ = 0) Γ pol ] where we average over equivalent vector components and directions: Tr[R (µ = )Γ pol ] = (Tr[ R(µ = ) Γ pol ] Tr[ R(µ = ) Γ pol ]). Tr[R (µ = 0)Γ pol ] = 3 (Tr[ R(µ = 0) Γ pol ] + Tr[ R(µ = 0) Γ pol ] + Tr[ 3 R(µ = 0) Γ pol ]). (4.3) 5. Results and conclusions We perform lattice QCD calculations using a tree-level Symanzik-improved gauge action [8, ] and + flavors of tree-level improved Wilson-clover quarks, which couple to the gauge links via two levels of HEX smearing, at the physical pion mass m π = 35 MeV, lattice spacing a = fm, and a large volume L 3 s L t = 64 4 satisfying m π L = 4. We are measuring the isovector combination u d of the three-point functions, where the disconnected contributions cancel out. Furthermore, we perform measurements using three source-sink separations T /a {0, 3, 6} 4

6 Computing the nucleon Dirac radius directly at Q = 0 Quantity Traditional method Derivative method κ v 3.74(4) 3.7(35) [r v] [fm] 0.55(9) 0.45(8) Table : A comparison between the derivative method and the traditional one ranging from 0.9 fm to.5 fm, and we are using the summation method [3, 4] for removing contributions from excited states. We apply our analysis on 44 gauge configurations, using all-mode-averaging [0] with 64 approximate samples and one high-precision bias correction per configuration. The plateau plots in Figure, show F v(0) (left-hand side, using 4.) and [rv ] /a ( right-hand side, using 4.) from the ratio method, as well as the summation method. Figure shows a comparison between the results we get for the anomalous magnetic moment κ = F v (0) and the isovector Dirac radius [r v] using the momentum derivative approach and what we get using the traditional approach of measuring Pauli and Dirac form factors for each value of Q (on the same ensemble) and then applying the z-expansion fit [, ]. Preliminary results are given in Table. We confirm that our approach produces results consistent with those obtained using the traditional method. However, we found that this approach yields large statistical uncertainties, especially for the Dirac radius, which requires two momentum derivatives applied to a single quark line. Therefore, putting one momentum derivative on each of two different quark lines, as was done in [7], might be less noisy F(0) v (r v ) /a T = 0a -0 T = 0a.5 T = 3a T = 3a T = 6a T = 6a summation -40 summation (t T /)/a 0 (t T /)/a Figure : Anomalous isovector magnetic moment (left) and isovector Dirac radius (right). For both κ and [r v ] /a, results from ratio method are shown using source-sink separations T /a {0,3,6}, as well as the summation method. Acknowledgments We thank the Budapest-Marseille-Wuppertal collaboration for making their configurations [8, ] available to us. This research used resources at Forschungszentrum Jülich and on CRAY XC40 (Hazel Hen) at HLRS. SM is supported by NSF grant PHY-50996, SM and SS are supported by RBRC, JN is supported by the Office of Nuclear Physics of the U.S. Department of Energy (DOE) under Contract DE-SC00090, ME is supported by DOE grant DE-FG0-96ER40965 and AP is supported in part by DOE under grant DE-FC0-06ER4444. Calculations for this project were done using the Qlua software suite [0]. 5

7 Computing the nucleon Dirac radius directly at Q = 0 Isovector F PRELIMINARY derivative method z-expansion fit standard method Kelly Q (GeV ) Isovector F PRELIMINARY z-expansion fit standard method derivative method Kelly Q (GeV ) Figure : Isovector Dirac (left) and Pauli (right) form factors. The blue points show results from the standard method and the red bands show a z-expansion fit to those points. The green band (left) and point (right) show the slope and value of the respective form factor at Q = 0, computed using the momentum derivative method. The black curves result from a phenomenological fit to experimental data by Kelly [9]. References [] A. Antognini, F. Nez, K. Schuhmann, F. D. Amaro, F. Biraben, et al., Science 339, 47 (03). [] P. J. Mohr, B. N. Taylor, D. B. Newell, Rev. Mod. Phys.84, 57(0) [arxiv:03.545] [3] I. T. Lorenz, Ulf-G. Meißner, Phys. Lett. B737, 57 (04) [arxiv:406.96] [4] K. Griffioen, C. Carlson, S. Maddox, Phys. Rev. C93, (06) [arxiv: ] [5] D. W. Higinbotham, A. A. Kabir, V. Lin, et al., Phys. Rev. C93, (06) [arxiv:50.093] [6] G. Lee, J. R. Arrington, R. J. Hill, Phys. Rev. D9, 0303 (05) [arxiv: ] [7] M. O. Distler, T. Walcher, J. C. Bernauer (05) [arxiv: ] [8] J. C. Bernauer, M. O. Distler, (06) [arxiv: ] [9] G. M. de Divitiis, R. Petronzio, N. Tantalo, Phys. Lett. B 78, 589 (0) [arxiv:08.594]. [0] T. Blum, T. Izubuchi, and E. Shintani, Phys. Rev. D88, (0) [arxiv: ]. [] R. J. Hill, G. Paz, Phys.Rev. D8, 3005 (00) [arxiv: ]. [] Z. Epstein, G. Paz, J. Roy, Phys.Rev. D90, (04) [arxiv: ]. [3] S. Capitani, B. Knippschild, et al., PoS LATTICE00, 47 (00) [arxiv:0.358] [4] J. Bulava, M.A. Donnellan, R. Sommer, PoS LATTICE00, 303 (00), [arxiv:0.4393] [5] K. C. Bowler, R. D. Kenway, et al., Phys. Rev. D57, 6948 (998) [arxiv: hep-lat/970908] [6] S. Capitani, M. Della Morte, et al., Phys. Rev. D9, 0545 (05), [arxiv: ] [7] B. C. Tiburzi, Phys. Rev. D90, (04), [arxiv: ] [8] S. Dürr, Z. Fodor, C. Hoelbling, S. Katz, S. Krieg, et al., JHEP, 08, 48 (00) [arxiv: 0.7] [9] J. J. Kelly, Simple parametrization of nucleon form factors, Phys. Rev. C 70, 0680 (004). [0] A. Pochinsky, Qlua. [] S. Dürr, Z. Fodor, C. Hoelbling, S. Katz, S. Krieg, et al., Phys. Lett. B70, 65 [arxiv: 0.403]. 6