-
The electromagnetic form factors of nucleons are important physical quantities, which are helpful to understand both the nucleon inner structure and the mechanism of strong interaction. For many years, the elastic lepton-nucleon scattering measurements have enabled us to obtain detailed information on the electromagnetic form factors of the proton and neutron over a wide range of kinematics. Numerous measurements were carried out at the Stanford Mark III accelerator [1], Cambridge Electron Accelerator [2], Stanford Linear Accelerator Center (SLAC) [3], Bonn [4], DESY [5], Mainz [6], NIKHEF [7], MIT-Bates [8], and Jefferson Lab [9]. The elastic scattering measurements have stimulated considerable activity over the past two decades in the determination of the flavor separated form factors of the up, down, and strange quarks in the nucleon. A number of measurements of the strange form factors have been successfully performed, starting with SAMPLE at Bates [10] and A4 at Mainz [11], followed by the high precision G0 [12] and HAPPEX [13] experiments at Jefferson Lab.
Numerous theoretical works have also been perfomed on the nucleon form factors. Owing to the non-perturbative property of the strong interaction, it is very difficult to study hadron properties using the QCD fundamental theory directly. Many theoretical calculations are based on phenomenological methods, such as the constituent quark model, cloudy bag model, perturbative chiral quark model, quark-diquark model, nonlocal quark meson coupling model, Nambu-Jona-Lasinio model, Dyson-Schwinger equations, vector meson dominance model, chiral quark soliton model, and AdS/QCD approach [14–30]. The lattice QCD simulation is the most rigorous method to study hadron properties. However, the simulated results have to be extrapolated to the continuum and infinite volume limits to compare them with experimental data. In addition, if the simulations are performed using heavy quark masses, the data need to be extrapolated to the physical quark masses. With the improvement in computing speed, currently, lattice simulation of nucleon form factors can be carried out at physical quark masses including the disconnected contributions [31, 32].
Besides lattice QCD, the chiral effective field theory (EFT) is another systematic approach that has been widely applied in hadron physics [33, 34]. Compared with the various phenomenological models, one advantage of the EFT is that it can provide some model-independent results. Although the chiral EFT has been a fairly successful approach, for the nucleon form factors it is only valid at relatively small
Q2 values, for example, less than 0.1 GeV2 [35]. The range can be extended up to 0.4 GeV2 by explicitly including vector meson degrees of freedom into the theory [36].In recent years, we proposed a nonlocal chiral effective theory, which makes it possible to study hadron properties at a relatively large momentum transfer [37]. The basic idea is that in the nonlocal Lagrangian, a baryon and a meson are located in different coordinates described by a correlator. To guarantee the local gauge invariance, the path integral of the gauge field is introduced, which leads to the additional interactions. The additional Feynman diagrams are crucial to obtain the corrected normalized nucleon charge. On the one hand, the correlator makes the loop integral ultraviolet convergent. On the other hand, the physical quantities at a large momentum transfer can be described very well. Perturbative calculations in chiral EFT expand the observables as series in the pseudoscalar meson mass
O(mϕ/Λχ) or small external momentumO(q/Λχ) , whereΛχ∼1 GeV is the scale associated with the chiral EFT. With the introducing of the regulator, the power counting of our nonlocal EFT is lost. The higher order terms are included as a type of resummation. This resummation of the chiral expansion induced through the introduction of a finite range cutoff in the momentum integrals of meson-loop diagrams has been shown to be a very good and effective resummation method. The nonlocal chiral EFT has been applied to study the nucleon electromagnetic form factors [38], strange form factors [39], unpolarized and polarized parton distribution functions (PDFs) [40–42], generalized parton distributions (GPDs) [43], etc. The approach has also been generalized to the fundamental interaction of QED, which provided an interesting explanation to the leptong−2 anomaly [44, 45].In the previous form factor calculations, the bubble and tadpole diagrams were not included [38, 39, 46]. In fact, in the earlier extrapolation of nucleon form factors with finite-range regularization, these diagrams were not included either [47, 48]. The bubble and tadpole diagrams have no effect on the nucleon wave function renormalization because their contributions to the nucleon charge are summed to be zero. In this work, we will apply the nonlocal EFT to study the nucleon form factors, including these bubble and tadpole diagrams, to verify whether the results are changed or not. This paper is organized in the following way. In Sec. II, we briefly introduce the nonlocal chiral effective Lagrangian. The matrix elements of the bubble and tadpole diagrams for the electromagnetic and strange form factors are presented in Sec. III. Numerical results are discussed in Sec. IV. Finally, Sec. V is a brief summary.
-
We start from the local chiral effective Lagrangian. The lowest Lagrangian that describes the pseudoscalar mesons, octet and decuplet baryons, and their interaction is expressed as
L(1)=Tr[ˉB(i⧸D−MB)B]+DTr[ˉBγμγ5{uμ,B}]+FTr[ˉBγμγ5[uμ,B]]+¯Tijkμ(iγμναDα+iMTσμν)Tijkν+C[ϵijk¯TilmμΘμν(uν)ljBmk+h.c.]+H¯Tijkμγμναγ5(uα)klTijlν+f24Tr[∂μU(∂μU)†],
(1) where D, F,
C , andH are the baryon-meson coupling constants.MB andMT are the octet and decuplet baryon masses.f≈93 MeV is the pseudoscalar decay constant. TheΘμν in the octet-decuplet transition operator is given byΘμν=gμν−Zγμγν,
(2) where Z is the off shell parameter. The tensor
γμνα is defined asγμνα=−i{σμν,γα} . The octet baryons and mesons are arranged in3×3 matrices, while the decuplet baryons are represented by the symmetric tensor with three indices. The covariant derivatives of the octet and decuplet baryon fields are given byDμB=∂μB+[Γμ,B],DμTijkν=∂μTijkν+(Γμ,Tν)ijk.
(3) The mesons couple to the baryon fields through the vector and axial vector combinations defined as
Γμ=12(u∂μu†+u†∂μu), uμ=i2(u∂μu†−u†∂μu),
(4) where
u is defined in terms of the pseudoscalar meson fieldϕ u2=exp(i√2ϕf).
(5) The octet, decuplet, and octet-decuplet transition magnetic moment operators are needed in the one loop calculation of nucleon form factors. The magnetic Lagrangian is written as
Lmag=e4MB(c1Tr[¯Bσμν{F†μν,B}]+c2Tr[¯Bσμν[F†μν,B]]+c3Tr[F†μν]Tr[¯BσμνB])−eFT24MT¯TabcμσαβFαβQceTμ,eba+iec44MBFμν(ϵijk¯BjmQilγμγ5Tν,klm+ϵijk¯Tμ,klmQliγνγ5Bmj),
(6) where
c1 ,c2 , andc3 describe the magnetic moments of octet baryons at tree levels, whereasFT2 andc4 are related to the magnetic moments of the decuplet baryons and octet-decuplet transition. With theSU(6) symmetry, one can have the following relationships [43]:c3=c2−c1,FT2=13(c1+3c2),c4=1√3c1.
(7) To calculate the contributions of the bubble diagrams, the following high order Lagrangian that provides the next-to-leading order baryon-baryon-meson-meson interaction has to be included [49]:
L(2)=i2σμν(b9Tr[¯Buμ]Tr[uνB]+b10Tr[¯B{[uμ,uν],B}]+b11Tr[¯B[[uμ,uν],B]]),
(8) where the coefficients
b9 ,b10 , andb11 are determined in Ref. [49] asb9=1.36 GeV,b10=1.24 GeV, andb11=0.46 GeV.The gauge invariant nonlocal Lagrangian can be obtained by the displacement of the coordinates of the meson and photon fields with the proper inclusion of the gauge link. For example, the gauge invariant local strong interaction
ppππ is expressed asLlocalππ=i4f2¯p(x)γμp(x)[π+(x)(∂μ−ieAμ(x))π−(x)−π−(x)(∂μ+ieAμ(x))π+(x)].
(9) The corresponding nonlocal Lagrangian is expressed as
Lnlππ=i4f2¯p(x)γμp(x)∫da∫dbF(a)F(b)[exp[ie∫x+axdzν∫dcAν(z−c)F(c)]π+(x+a)×(∂μ−ie∫dcF(c)Aμ(x−c))exp[−ie∫x+bxdzν∫dcAν(z−c)F(c)]π−(x+b)−(∂μ+ie∫dcF(c)Aμ(x−c))exp[ie∫x+axdzν∫dcAν(z−c)F(c)]π+(x+a)×exp[−ie∫x+bxdzν∫dcAν(z−c)F(c)]π−(x+b)], (10) where
F(a) is the correlation function. To guarantee the gauge invariance, the gauge links are introduced in the above Lagrangian. Similarly, the nonlocal electromagnetic interaction of the nucleon can also be obtained. For example, the local interaction between a proton and a photon is written asLlocalEM=−eˉp(x)γμp(x)Aμ(x)+(c1+3c2)e12mNˉp(x)σμνp(x)Fμν(x).
(11) The corresponding nonlocal Lagrangian is expressed as
LnlEM=−e∫daˉp(x)γμp(x)Aμ(x−a)F1(a)+(c1+3c2)e12mN∫daˉp(x)σμνp(x)Fμν(x−a)F2(a),
(12) where
F1(a) andF2(a) are the correlation functions for the nonlocal electric and magnetic interactions. As in our previous work, we assume that the charge and magnetic form factors at the tree level have the same momentum dependence as the nucleon-pion vertex, i.e.,GtreeM(p)=μpGtreeE(p)=μp˜F(p) , where˜F(p) is the Fourier transformation of the correlation functionF(a) . Therefore, the corresponding functions of˜F1(q) and˜F2(q) are expressed as˜F1(q)=˜F(q)4m2N+μpQ24m2N+Q2,˜F2(q)=˜F(q)(μp−1)4m2N4m2N+Q2.
(13) The nonlocal Lagrangian is invariant under the following gauge transformation
π+(y)→eiα(y)π+(y),p(x)→eiα(x)p(x),Aμ(x)→Aμ(x)−1e∂μα′(x),
(14) where
α(x)=∫daα′(x−a)F(a) .From the nonlocal Lagrangian, the interaction with the external photon can be obtained. There are two types of electromagnetic interactions in a nonlocal case. Besides that from the minimal substitution, which is the same as in the local interaction, there is an additional one generated from the expansion of the gauge link. For example, from Eq. (10), we can obtain the additional interaction as
Ladd=−e4f2¯p(x)γμp(x)∫da∫dbF(a)F(b)×[∫x+axdzν∫dcAν(z−c)F(c)π+(x+a)∂μπ−(x+b)−π+(x+a)∂μ(∫x+bxdzν∫dcAν(z−c)F(c)π−(x+b))−∂μ(∫x+axdzν∫dcAν(z−c)F(c)π+(x+a))π−(x+b)+∂μπ+(x+a)∫x+bxdzν∫dcAν(z−c)F(c)π−(x+b)].
(15) The additional interaction is important to acquire the renormalized proton (neutron) charge 1 (0). With the nonlocal Lagrangian, one can calculate the electromagnetic and strange form factors of the nucleon.
-
The Dirac and Pauli form factors of the nucleon are defined as
<N(p′)|Jμ|N(p)>=ˉu(p′){γμFN1(Q2)+iσμνqν2mNFN2(Q2)}u(p),
(16) where
q=p′−p andQ2=−q2 .FN1(Q2) andFN2(Q2) are the Dirac and Pauli form factors, respectively. The combination of the above form factors can generate the electric and magnetic form factors asGNE(Q2)=FN1(Q2)−Q24m2NFN2(Q2)GNM(Q2)=FN1(Q2)+FN2(Q2).
(17) In Eq. (16), the contributions of the u, d, and s quarks are all included in the current. If only the strange current is included, one can obtain the strange form factors. According to the Lagrangian, the one-loop Feynman diagrams that contribute to the nucleon electromagnetic form factors are plotted in Fig. 1, where only bubble and tadpole diagrams are shown. In Fig. 1, diagrams a and b are called bubble diagrams, while c, d, and e are called tadpole diagrams. Here, we only present the formulas for the bubble and tadpole diagrams. The formulas for the other diagrams can be found in Refs. [38, 46].
Figure 1. One-loop Feynman diagrams for the nucleon electromagnetic form factors. The solid and dashed lines are for a nucleon and meson, respectively. The crossed circle represents the interaction with the photon field from the minimal substitution, the black square represents the magnetic interaction in Eq. (6), the gray square denotes the interaction in Eq. (8), and the filled circle denotes the additional gauge link interaction with the photon field.
The contributions of the leading bubble diagram, Fig. 1(a), to the matrix elements in Eq. (16) are written as
Γμ(p)a=−14f2Iμa,π+π−−12f2Iμa,K+K−,Γμ(n)a=14f2Iμa,π+π−−14f2Iμa,K+K−,Γμ(s)a=14f2Iμa,K+K−+12f2Iμa,K0ˉK0,
(18) where the integral
Iμa,ϕ+ϕ− is expressed asIμa,ϕ+ϕ−=˜F(q)ˉu(p′)∫d4k(2π)4˜F(k)1Dϕ(k)×(2k̸+q̸)2kμ˜F(k+q)1Dϕ(k+q)u(p).
(19) Besides the matrix elements for the proton and neutron, the matrix element of the strange quark
Γμ(s)a is also written. In the above integral,Dϕ(k) is the propagator of meson ϕ. The correlators˜F(k) and˜F(k+q) make the loop integral convergent and˜F(q) gives a proper meson form factor at the tree level. The contributions of the bubble diagram at next-to-leading order, Fig. 1(b), to the matrix elements are written asΓμ(p)b=−2(b10+b11)f2Iμb,π+π−−4b11+b9f2Iμb,K+K−,Γμ(n)b=2(b10+b11)f2Iμb,π+π−−2(b10−b11)f2Iμb,K+K−,Γμ(s)b=−2(b10−b11)f2Iμb,K+K−+4b11+b9f2Iμb,K0ˉK0,
(20) where the integral
Iμb,ϕ+ϕ− is expressed asIμb,ϕ+ϕ−=˜F(q)ˉu(p′)∫d4k(2π)4σρνqρkν˜F(k)1Dϕ(k)×2kμ˜F(k+q)1Dϕ(k+q)u(p).
(21) For the tadpole diagram with electric coupling to the photon field, Fig. 1(c), the matrix elements are written as
Γμ(p)c=−4m2N+Q2(1+c1+c2)2(4m2N+Q2)f2Iμc,π+π−−4m2N+Q2(1+c2)(4m2N+Q2)f2Iμc,K+K−,Γμ(n)c=4m2N+Q2(1+c1+c2)2(4m2N+Q2)f2Iμc,π+π−−4m2N+Q2(1−c1+c2)2(4m2N+Q2)f2Iμc,K+K−,Γμ(s)c=4m2N+Q2(1−c1+c2)2(4m2N+Q2)f2Iμc,K+K−+4m2N+Q2(1+c2)(4m2N+Q2)f2Iμc,K0ˉK0,
(22) where the integral
Iμc,ϕ+ϕ− is expressed asIμc,ϕ+ϕ−=˜F(q)ˉu(p′)∫d4k(2π)4˜F2(k)1Dϕ(k)γμu(p).
(23) For the tadpole diagram with the magnetic coupling to the photon field, Fig. 1(d), the matrix elements are written as
Γμ(p)d=−2m2N(c1+c2)(4m2N+Q2)f2Iμd,π+π−−4m2Nc2(4m2N+Q2)f2Iμd,K+K−,
Γμ(n)d=2m2N(c1+c2)(4m2N+Q2)f2Iμd,π+π−+2m2N(c1−c2)(4m2N+Q2)f2Iμd,K+K−,Γμ(s)d=−2m2N(c1−c2)(4m2N+Q2)f2Iμd,K+K−+4m2Nc2(4m2N+Q2)f2Iμd,K0ˉK0,
(24) where the integral
Iμd,ϕ+ϕ− is expressed asIμd,ϕ+ϕ−=˜F(q)ˉu(p′)∫d4k(2π)4˜F2(k)1Dϕ(k)σμνqν2mNu(p).
(25) In the nonlocal chiral EFT, there is an additional interaction with the external photon, which is generated from the expansion of the gauge link. The contributions of this additional diagram, Fig. 1(e), to the matrix elements are written as
Γμ(p)e=−12f2Iμe,π+π−−1f2Iμe,K+K−,Γμ(n)e=12f2Iμe,π+π−−12f2Iμe,K+K−,Γμ(s)e=12f2Iμe,K+K−+1f2Iμe,K0ˉK0,
(26) where the integral
Iμe,ϕ+ϕ− is expressed asIμe,ϕ+ϕ−=˜F(q)ˉu(p′)∫d4k(2π)4˜F(k)1Dϕ(k)×2k̸(2k+q)μ2k⋅q+q2[˜F(k+q)−˜F(k)]u(p).
(27) After simplifying the γ matrix algebra, we can obtain the separate expressions for the Dirac and Pauli form factors. Together with the contributions from the rainbow and Kroll-Rudderman diagrams, which have been obtained in our previous work [38, 46], the nucleon form factors can be obtained. In the next section, the numerical results will be discussed.
-
In the numerical calculations, the parameters are selected as
D=0.75 ,F=23D , andH=−3D . The coupling constantC is set to be1.0 , which is the same as that in Refs. [38, 50]. The parametersc1 andc2 are determined by the experimental magnetic moments of the proton and neutron,μp=2.79 andμn=−1.91 . The covariant regulator for the baryon-meson interaction is chosen to be a dipole form\tilde{F}(k)=\left(\frac{\Lambda^{2}-m_{\phi}^{2}}{\Lambda^{2}-k^{2}}\right)^{2},
(28) where
m_\phi is the meson mass. For the photon field, the regulator is the same, except that the meson mass is replaced by zero. In our calculation, Λ is the only free parameter. It is found that when Λ is approximately 1.0 GeV, the obtained nucleon form factors are reasonable compared with the experimental data. The determined parametersc_1 andc_2 are listed in Table 1. For both cases, with or without bubble and tadpole diagrams, we can reproduce the experimental nucleon magnetic moments with proper choices ofc_i .c_1 andc_2 will be slightly smaller when the bubble and tadpole diagrams are included.Λ /GeV c_{1} c_{2} c_{1} c_{2} 0.9 2.05 0.75 1.77 0.67 1.0 1.98 0.74 1.61 0.63 1.1 1.92 0.72 1.46 0.59 Table 1. Obtained parameters
c_{1} andc_{2} . The first and second twoc_i (i=1,2) are for the cases without and with bubble and tadpole diagrams, respectively.With the determined parameters, we will show the results of the electromagnetic and strange form factors of the proton and neutron. The proton electric form factor versus momentum transfer is plotted in Fig. 2. The blue and red bands are the results for Λ
1.0 \pm 0.1 GeV with and without bubble and tadpole diagrams, respectively. The experimental data are also shown in the figure. From the figure, one can observe that the experimental data can be reasonably described in both cases. The obtained charge form factorG_E^p(Q^{2}) with bubble and tadpole contributions is slightly smaller than that without bubble and tadpole contributions. With the additional diagram generated from the interaction of gauge link, the proton charge is 1 atQ^2=0 . Compared with the results with dimensional and infrared regularizations, the nonlocal approach can describe the form factor data at a relatively large momentum transfer. The calculated magnetic form factor of the protonG_{M}^{p}(Q^{2}) is shown in Fig. 3. In both cases, the theoretical bands are also consistent with the experimental data up to 1 GeV^2 . The band with bubble and tadpole contributions is slightly smaller than that without those contributions, which is closer to the experimental data as the electric form factor case.Figure 2. (color online) Electric form factor of the proton
G_{E}^{p}(Q^{2}) versus momentum transferQ^2 with\Lambda=1.0 \pm 0.1 GeV. The blue and red bands are for the results with and without bubble and tadpole diagrams, respectively.Figure 3. (color online) Same as Fig. 2 but for the magnetic form factor.
As in the proton case, the electromagnetic form factors of the neutron,
G_{E}^{n}(Q^{2}) andG_{M}^{n}(Q^{2}) , are plotted in Fig. 4 and Fig. 5. ForG_E^n(Q^2) , both bands first increase and then smoothly decrease with the increasing momentum transfer. The shapes are both comparable with the experimental data. The band with bubble and tadpole contributions is higher than that without those contributions. Obviously, the bubble and tadpole contributions make the theoretical result more reasonable. With the additional diagrams, the neutron charge is 0 atQ^2=0 . We should mention that the neutron electric form factor in the case without bubble and tadpole diagrams is larger than our previous result [38]. This is mainly because the parameters are a little different. In particular, here the cutoff Λ is chosen to be approximately 1.0 GeV, while it is 0.85 GeV in Ref. [38]. ForG_M^n(Q^2) , there is a slight difference between these two bands, although both of them can describe the experimental data well.Figure 4. (color online) Electric form factor of the neutron
G_{E}^{n}(Q^{2}) versus momentum transferQ^2 with\Lambda=1.0 \pm 0.1 GeV. The blue and red bands are for the results with and without bubble and tadpole diagrams, respectively.Figure 5. (color online) Same as Fig. 4 but for the magnetic form factor.
From the above numerical results, one can see that for the electromagnetic form factors, even though the determined low energy constants
c_1 andc_2 are different, the calculated results with and without bubble and tadpole contributions are both comparable with the experimental data. With the determined parameters, we then present the results for the strange form factors.The strange electric form factor of the nucleon,
G_{E}^{s}(Q^{2}) , for Λ1.0 \pm 0.1 GeV is shown in Fig. 6. As for the nucleon electromagnetic case, the blue and red bands are for the results with and without bubble and tadpole contributions, respectively. The lattice data from Ref. [31] are also shown. In both cases, the strange electric form factor first increases from zero and then changes smoothly with the increasingQ^2 . The numerical difference of these two cases is obvious. The result with the bubble and tadpole contributions is much larger than that without these contributions. For example, atQ^2=0.22 GeV^2 ,G_E^s (0.22 GeV^2 )= 0.0085_{-0.0037}^{+0.0046} , which is approximately four times larger than the strange electric form factor without bubble and tadpole contributions. This is because the contributions from the bubble and tadpole diagrams are much larger than those from the diagrams with octet and decuplet intermediate states. We should mention that after we fixed the parametersc_1 andc_1 by the magnetic moments of the proton and neutron, the strange form factors are calculated directly without any parameter. Although our values are larger than those of the lattice simulation at finiteQ^2 , the numbers are consistent with the global analysis whereG_E^{s,\text{exp}}(0.22\; \text{GeV}^2) =0.035\pm 0.030 \pm 0.019 [51].Figure 6. (color online) Strange electric form factor of the nucleon
G_{E}^{s}(Q^{2}) versus momentum transferQ^2 with\Lambda=1.0 \pm 0.1 GeV. The blue and red bands are for the results with and without bubble and tadpole diagrams, respectively. The lattice data are from Ref. [31].The strange magnetic form factor is plotted in Fig. 7 for the two cases together with the lattice data. For the magnetic case, the contributions from bubble and tadpole diagrams are comparable with those from the rainbow and Kroll-Rudderman diagrams. As a result, the two bands are both negative and comparable. The strange magnetic moment with bubble and tadpole contribution is
\mu_s = -0.066_{-0.034}^{+0.028} , which is also covered by the experimental analysis\mu_s^{\text{exp}} = -0.14 \pm 0.11 \pm 0.11 [51].Figure 7. (color online) Same as Fig. 6 but for the magnetic form factor.
-
We applied the nonlocal EFT to calculate the electromagnetic and strange form factors of the nucleon. The bubble and tadpole diagrams were included in the calculation. The next-to-leading order baryon-meson interaction
{\cal L}^{(2)} was necessary, even though it was not needed in the previous calculation with only rainbow and Kroll-Rudderman diagrams. In the numerical calculation, all the parameters were predetermined except the cutoff parameter Λ in the regulator and the low energy constantsc_1 andc_2 .c_1 andc_2 were fixed by the nucleon magnetic moments, and Λ was selected to be approximately 1.0 GeV. The numerical results showed that the electromagnetic form factors with and without bubble and tadpole diagrams were close to each other with the proper choices ofc_1 andc_2 . In both cases, the electromagnetic form factors were in good agreement with the experimental data. As there was no free parameter when we calculated the strange form factors, the obtained strange form factors could have a significant difference in the two cases. With the bubble and tadpole diagrams, the magnitudes of the strange electric and magnetic form factors were both larger. The obtained magnetic moment is\mu_s =-0.066_{-0.034}^{+0.028} . The strange electric form factor is positive, andG_E^s (0.22 GeV^2 )= 0.0085_{-0.0037}^{+0.0046} whenQ^2=0.22 GeV^2 . Our results revealed that the magnitudes of the strange form factors are larger than the lattice data.
Contributions to the nucleon form factors from bubble and tadpole diagrams
- Received Date: 2022-10-17
- Available Online: 2023-02-15
Abstract: The nonlocal chiral effective theory is applied to investigate the electromagnetic and strange form factors of nucleons. The bubble and tadpole diagrams are included in the calculation. With the contributions from bubble and tadpole diagrams, the obtained electromagnetic form factors are close to the results without these contributions as long as the low energy constants