Chiral extrapolation of nucleon axial charge gA in effective field theory

The extrapolation of nucleon axial charge gA is investigated within the framework of heavy baryon chiral effective field theory. The intermediate octet and decuplet baryons are included in the one loop calculation. Finite range regularization is applied to improve the convergence in the quark-mass expansion. The lattice data from three different groups are used for the extrapolation. At physical pion mass, the extrapolated gA are all smaller than the experimental value.


Introduction
The nucleon axial charge, g A , is a fundamental property of the nucleon, which reveals how the up and down quark intrinsic spins contribute to the spin of the proton and neutron, governing β decay and providing a quantitative measure of spontaneous chiral symmetry breaking in low energy hadronic physics. The axial charge g A is of great importance to any further calculation of hadron structure.
The axial charge g A is defined as the axial vector form factor at zero four-momentum transfer, g A = G A (0). The axial vector form factor is given by the nucleon matrix element of the axial vector current, A a µ = ψγ µ γ 5 (τ a /2) ψ, with u, d quark doublet ψ, N (p , s )| A 3 µ |N (p, s) = iū(p , s )[γ µ γ 5 G A (q 2 ) + q µ 2M N γ 5 G P (q 2 )] τ 3 2 u(p, s), where G P is the induced pseudoscalar form factor, τ a is an isospin Pauli matrix, and q µ = p µ − p µ is the momentum transfer. At zero momentum transfer, the axial charge g A is the spin difference between u and d quarks in the proton, i.e. g A = ∆u − ∆d.
Experimentally, g A has been obtained very precisely through neutron β decay, with the Particle Data Group value g A = 1.27±0.003 [1]. Theoretically, there are many calculations in different methods, such as the cloudy bag model [2], the perturbative chiral quark model [3], the relativistic constituent quark model [4], Schwinger-Dyson formalism [5], chiral perturbation theory [6], etc. There are also many lattice simulations of axial charge [7][8][9][10][11][12]. Due to the limitations of computing ability, all the simulations of g A are at large quark mass. The obtained g A at large quark mass are smaller than the experimental data. Therefore, it is interesting to see how the axial charge g A changes at low pion mass.
In this paper, we will extrapolate nucleon axial charge g A in the framework of heavy baryon chiral perturbation theory with finite range regularization (FRR). FRR has been applied in the extrapolation of nucleon mass, magnetic form factors, strange form factors, charge radii, first moments, etc [13][14][15][16][17][18][19][20][21][22][23][24][25][26][27]. It is proved that FRR can provide good convergent behaviour of pion mass expansion. Therefore, it is expected to have a good description of the pion mass dependence of axial charge g A over a wide range of pion mass.

Nucleon axial charge
The lowest-order chiral Lagrangian including the octet and decuplet baryons is expressed as where S µ v is the covariant spin operator defined as Here, v ν is the nucleon four-velocity. In the rest frame, we have v ν = (1, 0, 0, 0)). D, F and C are the standard SU (3)-flavour coupling constants. According to the Lagrangian, the one-loop Feynman diagrams, which contribute to axial charge g A of the proton, are plotted in Fig. 1. The axial charge is the spin difference between u and d quark. The contributions of the u-and d-quark sector to the proton spin, from Fig. 1a, are expressed as where the coefficients, C are expressed as The tree level contributions to the proton spin from u and d quark of intermediate octet baryons are used in the above formulas. For example, for the intermediate proton and neutron, their spins are expressed as s u and s d are the single quark spin of u and d quark.
With the SU (2) symmetry, s u = s d = s q and s q can be written as where c 0 , c 2 and c 4 are the low energy constants. The contribution of u-, d-quark sector to the proton spin, described by diagram (b) of Fig. 1, are expressed as where the coefficients C ∆π and C Σ * K are (14) Similar to the case of the octet intermediate state, the tree level quark contributions to the spin of decuplet baryons are also used. For example Diagrams (c) and (d) of Fig. 1 provide contributions from intermediate states involving an octet-decuplet transition. The u-, d-quark-sector contribution to the proton spin from these diagrams are expressed as where The integrals in the above equations, I αβ 2j , I αβγ 5j and I αβ 3j are defined in Ref. [15]. The total u-, d-quark sector contributions to the spin of the proton are written as ∆u = Z 4 3 (c 0 + c 2 m 2 π + c 4 m 4 π ) + ∆u a + ∆u b + ∆u c+d , where Z is the wave function renormalization constant, expressed as The above coefficients are expressed as The K-and η-meson masses have relationships with the pion mass as This enable a direct relationship between the nucleon axial charge and the pion mass. By fitting the lattice data with different pion mass, we can get the low energy constants c 0 , c 2 and c 4 . In our calculation, the one-gluon-exchange is also included. Although it lies outside the framework of chiral effective field theory, the effect of one-gluon-exchange (OGE) is particularly important for spin dependent quantities. Hogaason and Myhrer [28] showed that the incorporation of the exchange current correction arising from the effective one-gluon-exchange (OGE) force shifts the tree-level non-singlet charge, g A , from 5 3 s q to 5 3 s q −G, where G is about 0.05. The OGE correction shifts the tree-level singlet charge g 0 from s q to s q − 3G. In other words, the spin of each constituent quark gain a OGE correction −G at tree level.
In Fig. 2, the pion mass dependence of g A with Λ = 0.8 GeV is shown for lattice data of Ref. [7]. The dotted, dashed and solid lines are for tree level, loop and total contribution, respectively. At large pion mass, the axial charge g A changes little. At small pion mass, g A decreases with the decreasing pion mass. Compared with the pion mass dependence of proton magnetic form factors [15,23], at low pion mass, the curvature is small and opposite. This is because the leading diagram in the case of magnetic form factor has no contribution for g A . At physical pion mass, the extrapolated g A is 1.10, which is smaller than the experimental value 1.27. To provide an estimate of the uncertainty in these results, we vary the regulator parameter, Λ, governing the size of meson cloud contributions to proton structure. Considering Λ = 0.8 ± 0.2 GeV, the obtained low energy constants c 0 , c 2 , c 4 as well as the quark spin at physical pion mass are listed in Table 1. By varying Λ, we can provide an error bar for g A . For example, the highest and lowest g A at physical pion mass are 1.14 (0.805 − (−0.333)) and 1.07 (0.772 − (−0.302)). From the table, one can see that the loop/tree contribution increases/decreases with the increasing Λ. The highest and lowest value of g A versus pion mass as well as the central value of g A are shown in Fig. 3. It is clear that the extrapolated g A with error bar is still smaller than the experimental value.
There are also other lattice groups simulating the axial charge g A . Figure 4 and Fig. 5 are results for the lattice data from Refs. [8][9]. The same as in Fig. 3, the lines in the middle are for Λ = 0.8 GeV. The upper and lower lines are obtained by varying Λ from 0.6 to 1 GeV. The extrapolated g A from Ref. [8] at physical pion mass is 1.12 +0.03 −0.03 . The lattice data from Ref. [9] varied a lot with the change of the pion mass though the extrapolated g A at physical pion mass is a little larger than the other two lattice groups. At large pion mass, Fig. 4 and Fig. 5 show that g A changes quickly with the increasing pion mass for the data of ETMC and data from Ref. [9]. This is because, different from the data of LHPC, there is no con-straint from these lattice data at large pion mass. Overall, one can see the results from different lattice groups are comparable and all the extrapolated g A at physical pion mass are smaller than the experimental values, within the error bars. The obtained results with central Λ = 0.8 GeV for these three lattice groups are listed in Table 2.

Summary
In summary, we extrapolated the axial charge g A in chiral effective field theory with finite range regularisation. The dipole regulator is used as our previous extrapolation for nucleon mass, form factors, first moments, etc. The lattice data are from three lattice groups where the volume corrections are given explicitly. Different from the proton magnetic form factor, the axial charge g A decreases with decreasing pion mass when m π is small. The lattice data over a wide pion mass range can be well described with the FRR chiral effective field theory. At physical pion mass, the extrapolated g A are comparable to each other and all of them are smaller than the experimental value. To estimate the error bars for the extrapolation, we vary Λ in the regulator from 0.6 to 1 GeV. The upper limit of the extrapolated g A at physical pion mass is still smaller than the experimental value.