Sea quark contributions to nucleon electromagnetic form factors with the nonlocal chiral effective Lagrangian

The sea quark contributions to the nucleon electromagnetic form factors of the up, down and strange quarks are studied with the nonlocal chiral effective Lagrangian. Both octet and decuplet intermediate states are included in the one loop calculations. Compared with the strange quark form factors, although their signs are the same, the absolute value of the light quark form factors are much larger. For both the electric and magnetic form factors, the contribution of the d quark is larger than of the u quark. The current lattice simulations of the light sea quark form factors are in between our results for the u and dquarks.


Introduction
Nucleon structure is one of the hottest research topics in hadron physics, and has attracted considerable experimental and theoretical effort. With the upgrade of the experimental facilities, more information about the structure of the nucleon will become available. By increasing the energy, experimentalists can obtain the parton distribution functions (PDFs) from deep inelastic scattering, as well as the form factors at relatively large momentum transfer from elastic scattering. In addition to the valence quark, both the analytical and numerical information about the sea quark contribution to nucleon properties is also very important, as it is crucial for determining the leading non-analytic behavior of physical quantities.
The strange quark contribution to the nucleon form factor is of special interest because it is purely from sea quarks. Experiments by different collaborations have been carried out to measure this quantity precisely [1][2][3][4]. There have also been many theoretical discussions of the strange quark form factors [5][6][7]. In 2002, we studied the strange quark form factors with the perturbative chiral quark model (PCQM), and it was found that the strange quark electric form factor is positive while the magnetic form factor is negative [8]. At that time, the theoretical predictions were quite different as there were no precise experimental measurements available.
It is impossible to study the nucleon structure using quantum chromodynamics directly because of its nonperturbative nature. Besides the phenomenological quark models, two methods are systematically used in hadron physics, the lattice simulations and the effective field theory or the chiral perturbation theory. Traditionally, the chiral perturbation theory with dimensional regularization (DR) is valid only at low momentum transfer GeV 2 [9]. It can be applied up to 0.4 GeV 2 if vector mesons are included [10]. In lattice simulations, many quantities simulated on the lattice are for large quark (pion) mass, and it is necessary to extrapolate the lattice results to the physical pion mass. Instead of DR, we applied the effective field theory with finite range regularization (FRR). The vector meson mass, magnetic moments, magnetic form factors, strange quark form factors, charge radii, first moments of GPDs, nucleon spin etc., could be successfully extrapolated to the physical mass [11][12][13][14][15][16][17][18][19][20]. The strange quark form factors obtained with FRR are also consistent with our previous results with PCQM [12,15]. Recently, the nucleon form factors, as well as the sea contribution of the light and strange quarks, were simulated on the lattice for the physical pion mass [21][22][23][24]. Therefore, it would be interesting to compare these result with those calculated in the framework of the effective field theory.
In recent years, we proposed a nonlocal chiral effective Lagrangian which makes it possible to study the hadron properties at relatively large [25][26][27][28][29][30]. The nonlocal interaction generates both the regulator which makes the loop integral convergent, and the dependence of the form factors at tree level. The obtained electromagnetic form factors and strange quark form factors of the nucleon are very close to the experimental data [28,29]. In this paper, we apply the nonlocal Lagrangian to investigate the light sea quark contribution to the nucleon form factors. With the quark flow method as in Ref. [30], we obtain the sea and valence quark contributions separately. This method is equivalent to the quenched chiral perturbation theory. In Sec. 2, we introduce the nonlocal chiral Lagrangian. The sea quark contributions to the nucleon form factors are derived in Sec. 3. Numerical results are shown in Sec. 4, and finally, a short summary is given in Sec. 5.

Chiral effective Lagrangian
The lowest order chiral Lagrangian for baryons, pseudo-scalar mesons and their interactions can be written as [28, [31][32][33], where D, F and C are the coupling constants. The chiral covariant derivative is defined as . The pseudo-scalar meson octet couples to the baryon field via the vector and axial vector combinations as where is the real charge matrix . and B are the matrices of pseudo-scalar fields and octet baryons.
is the photon field. The covariant derivative in the is the chiral connection defined as . , are the antisymmetric matrices expressed as The octet, decuplet and octet-decuplet transition magnetic moment operators are needed in the one loop calculations of the nucleon electromagnetic form factors. The baryon octet anomalous magnetic Lagrangian is written as where, At the lowest order, the contribution of quark q to the nucleon magnetic moments can be obtained by replacing the charge matrix with the corresponding diagonal quark matrices . After expansion of the above equation, it is found that Comparing with the results of the constituent quark model where F p,s 2 = 0 and F n,s we get The transition magnetic operator is The effective decuplet anomalous magnetic moment operator can be expressed as The anomalous magnetic moments of baryons can also be expressed in terms of quark magnetic moments . For example, , , . Using the SU(3) symmetry, , and as well as can be written in terms of or . For example, , , .
π The gauge invariant nonlocal Lagrangian can be obtained using the method given in [26,28,29]. For instance, the local interaction including the meson can be written as Chinese Physics C Vol. 44, No. 5 (2020) 053101 The corresponding nonlocal Lagrangian is expressed as where is the correlation function. To guarantee gauge invariance, the gauge link is introduced in the above Lagrangian. The regulator can be generated automatically with the correlation function.
With the same procedure, the nonlocal electromagnetic interaction can also be obtained. For example, the local interaction between the proton and photon is written as The corresponding nonlocal Lagrangian is expressed as where and are the correlation functions for the nonlocal electric and magnetic interactions.
The form factors at tree level which are momentum dependent can be easily obtained with the Fourier transformation of the correlation function. As in our previous work [28,29], the correlation function is chosen such that the charge and magnetic form factors at tree level have the same momentum dependence as the nucleon-pion vertex, i.e.
, where is the Fourier transformation of the correlation function . The corresponding functions and are then expressed as where is the momentum transfer. From Eq. (13), two kinds of couplings between hadrons and photons can be obtained. One is the normal coupling, expressed as This interaction is similar to the traditional local Lagrangian except for the correlation function. The other is the additional interaction obtained by the expansion of the gauge link, expressed as

Electromagnetic form-factors
The contribution of the quark flavor f to the Dirac and Pauli form factors of the nucleon are defined as where . The electromagnetic form factors are defined as combinations of the above form factors for each flavor as In this work, we investigate the contribution of the u , d and s sea quarks to the nucleon electromagnetic form factors. According to the Lagrangian, the one loop Feynman diagrams which contribute to the nucleon electromagnetic form factors are shown in Fig. 1.
The coupling constants between baryons and mesons (coefficients) in Fig. 1 can be obtained from the Lagrangian. For each diagram in Fig. 1, there exist quenched and disconnected diagrams. In order to obtain the pure sea quark contribution, we need to get the coefficients for the disconnected diagrams. The coefficients for the quenched and disconnected loop diagrams can be obtained separately as in Ref. [30], using the quark flows shown in Fig. 2. The obtained coefficients are the same as those extracted with the graded symmetry formalism in the quenched chiral perturbation theory [34]. In Fig. 2, we plot the rainbow diagram using quark flows, as an example of the method of separating the quenched and sea quark contributions. The coefficients for the loop diagram in full QCD is . The coefficient in Fig.  2(b) for the sea quark contribution is the same as in Fig.  2(c)  The sea quark contributions in Fig. 1a for u d and s quarks are written as where the integral is expressed as where the integral is written as where is expressed as Figures 1d and 1e are the Kroll-Ruderman diagrams. The contributions of these two diagrams are written as where where γ Using Package-X [35] to simplify the matrix algebra, we can get the expressions for the Dirac and Pauli form factors. In the next section, we discuss the numerical results.

Numerical results
In the numerical calculations, the parameters are chosen as and ( ). The coupling constant is chosen as 1, the same as in Refs. [28,36]. The off-shell parameter z is [37]. The low energy constants and are determined as and , giving the experimental moments and . The covariant regulator is of the dipole form [28, 29, 31]F where is the meson mass for the baryon-meson interaction and is zero for the hadron-photon interaction. It was found that for around 0.90 , the results are very close to the experimental nucleon form factors. Λ = 1.0, 0.9 0.8 GeV In Fig. 3, we plot the sea contribution of the u quark with unit charge to the proton electric form factor. The three blue lines from top to bottom are for and . As a comparison, the central result for the strange quark is also plotted in the figure with the red line. The solid dots with error bars are the lattice simulations from Ref. [22]. Since we do not include the valence contribution of the u quark in the proton, the electric form factor of the u quark is zero for . It then increases with increasing . For larger than about 0.3 , it deceases with . From the figure, one can see that the strange quark form factor can be described very well. The u quark result is slightly smaller than the lattice simulations. We should mention that the lattice results are for the light quarks and it was assumed that u and d have the Chinese Physics C Vol. 44, No. 5 (2020) 053101 053101-5 same sea contributions. Therefore, the lattice results for the light quarks can be approximately treated as an average of the u and d contributions.
The d quark sea contribution to the proton electric form factor is plotted in Fig. 4. Similar to the u quark, the electric form factor first increases from 0 and then decreases with increasing . It can be clearly seen that the calculated sea quark contribution is larger than the lattice results. The larger sea contribution of the d quark than of the u quark is due to the fact that there is no intermediate octet contribution of the u quark, for which the only contribution is from the decuplet intermediate states. Similar results are found for the asymmetry in the proton, where is in excess with respect to [38][39][40][41]. Although there is an obvious difference between the sea quark contributions of u and d, both are much larger than the strange quark contribution. The strange quark electric form factor is about 5-10 times smaller due to the sup-pression of the K meson loop.
The sea contribution to the proton magnetic form factor of the u and d quarks with unit charge are plotted in Fig. 5 and Fig. 6. The calculated strange quark magnetic form factor is again in good agreement with the lattice results. All quark magnetic form factors increase monotonously with increasing . The absolute values of the u quark contribution are smaller than the lattice simulation results for the light quarks, while for the d quark contribution, the absolute values are larger than the lattice results. The absolute contributions of both u and d quarks are larger than of the strange quark, especially for small . For , the magnetic moments of the u and d sea quarks are and , respectively, while the strange quark magnetic moment is about .
π 0 π + From the above figures, it can be seen that the u and d sea quark contributions are quite different. For both the electric and magnetic form factors, the absolute value of the d sea quark contribution is much larger than of u, because there are two up quarks and one down quark in a proton. The u quark in the loop diagram can only form a decuplet state and there is no intermediate octet contribution to the u sea quark form factors. We should mention that the difference between the light sea quark form factors is not due to the mass difference of u and d quarks. In fact, in our calculation the masses of , and  GeV. The points with error bars are from lattice simulations [22].  GeV. The points with error bars are from lattice simulations [22].   It is straightforward to see that the sea quark form factor of u (d) in the proton is the same as that of d (u) in the neutron, if the masses of proton and neutron are the same. The mass difference between proton and neutron leads to a small charge symmetry violation, i.e. a small difference between and ( and ). The large difference between and is due to the effect of non-perturbative valence quark rather than of the mass difference between u and d. If the mass of the three u quark state is taken to be degenerate with respect to the nucleon mass, and the mass is taken to be degenerate with the mass. the sea contribution of u and d quarks in the proton would be the same. This is an artifact of the present lattice simulations. Physically, three u quarks cannot form an octet baryon, and the mass of is much larger than of . Therefore, it would be very interesting and challenging to get the flavor asymmetry from the full lattice QCD simulations.

Summary
In this work, we applied the nonlocal chiral effective Lagrangian to study the sea quark contribution of light quarks to the proton electromagnetic form factors. Since the sign of the u, d and s quark form factors are the same, this calculation helps to understand the experimental values of strange quark form factors. It is also interesting to compare our results with the lattice simulations. In our calculations, the parameter in the regulator is the same as in the lattice simulations, and was determined by fitting the nucleon form factors. The low energy constants and were determined from the experimental magnetic moments of the proton and neutron. Therefore, in calculating the sea quark form factors, there are no free parameters to be adjusted. Our results show that the electric form factors of light sea quarks with unit charge are positive, while the magnetic form factors are negative. Compared with the strange quark form factors, the absolute values of the light quark form factors are much larger. For both the electric and magnetic form factors, the contribution of the d quark is larger than of the u quark. The current lattice simulations for the light sea quark form factors are in between our results for u and d. Therefore, it would be interesting if this flavor asymmetry could be obtained from the full lattice QCD simulations.