Electric quadrupole form factors of singly heavy baryons with spin 3/2

We study the electromagnetic form factors of the lowest-lying singly heavy baryons in a pion mean-field approach, which is also known as the SU(3) chiral quark-soliton model. In the limit of the heavy-quark mass, the dynamics inside a singly heavy baryon is governed by the $N_c-1$ valence quarks, while the heavy quark remains as a mere static color source. In this framework, a singly heavy baryon is described by combining the colored soliton with the singly heavy quark. In the infinitely heavy-quark mass limit, we can compute the electric quadrupole form factors of the baryon sextet with spin 3/2 with the rotational $1/N_c$ and linear corrections of the explicit flavor SU(3) symmetry breaking taken into account. We find that the sea-quark contributions or the Dirac-sea level contributions dominate over the valence-quark contributions in lower $Q^2$ region. We examined the effects of explicit flavor SU(3) symmetry breaking in detail. The numerical results are also compared with the recent data from the lattice calculation with the unphysical value of the pion mass considered, which was used in the lattice calculation.


I. INTRODUCTION
Conventional lowest-lying singly heavy baryons consist of a heavy quark and two light valence quarks. In the limit of the infinitely heavy-quark mass (m Q → ∞), the physics of singly heavy baryons becomes simple: The spin of the heavy quark J Q is conserved in this limit and hence it leads also to the conservation of the spin of the light-quark degrees of freedom, i.e. J L = J − J Q . This is known as the heavy-quark spin symmetry [1,2]. In the m Q → ∞ limit, we do not distinguish a charm quark from a bottom quark, which gives up heavy-quark flavor symmetry. On the other hand, chiral symmetry and its spontaneous breakdown still play an important part in describing the singly heavy baryons because of the presence of the light quarks inside a singly heavy baryon [3]. Then the singly heavy baryons consisting of two light valence quarks can be represented in terms of irreducible representations of flavor SU(3) symmetry: 3 ⊗ 3 =3 ⊕ 6, thus we have the two representations for the lowest-lying singly heavy baryons, i.e. the baryon antitriplet and sextet. The baryon antitriplet has the total spin J = 1/2 that comes from J Q = 1/2 and J L = 0, whereas the baryon sextet can have either J = 1/2 or J = 3/2 with J L = 1 and J Q = 1/2.
In a pion mean-field approach, which is also known as the SU(3) chiral quark-soliton model (χQSM), a singly heavy baryon can be viewed as the N c − 1 valence quarks bound by the pion mean fields that are created from the presence of the N c − 1 valence quarks [4,5]. In fact, this idea is taken from Witten's seminal paper on baryons in the large N c limit [6]. This pion mean-field approach has successfully reproduced the mass spectra of the lowest-lying singly heavy baryons [5] and even explained the nontrivial isospin mass splittings of them [7]. Interestingly, the corrections from the heavy quark mass are indeed negligible in the description of the isospin mass splittings of the singly heavy baryons as shown in Ref. [7], although they provide hyperfine interactions to remove the spin degeneracy of the baryon sextet.
Recently, the electromagnetic (EM) form factors of singly heavy baryons have been studied for the first time within lattice QCD [8,9]. While there are no experimental data on the EM form factors of the singly heavy baryons to date, the results from the lattice calculation provide a clue to the internal structure of singly heavy baryons. Thus, it is also of great importance to investigate the EM form factors of the singly heavy baryons. In Ref. [12,13], we have studied the electric monopole and magnetic dipole form factors of the singly heavy baryons in detail, based on the χQSM. Since we consider the limit of the infinitely heavy-quark mass, there is no physical difference between the heavy baryons with spin 1/2 and those with 3/2 except for the value of the spin. On the other hand, the baryon sextet with spin 3/2 has yet another structure that arises from its higher spin, which is revealed by the electric quadrupole (E2) form factor. The E2 form factor of a baryon exhibits how the baryon is deformed. It is also known that the pion clouds play a significant role in understanding this deformation [14]. This will be also discussed in the present work. We will also examine the effects of flavor SU(3) symmetry breaking on the E2 form factors of the baryon sextet with spin 3/2. The numerical results for Ω * 0 c will be compared with that from the lattice calculation. The present work is organized as follows: In Section II, we briefly recapitulate the general formalism for the electric quadrupole form factors within the framework of the chiral quark-soliton model. In Section III, we present the numerical results and discuss them in detail. The final Section is devoted to the summary and conclusion.

II. ELECTRIC QUADRUPOLE FORM FACTORS IN THE χQSM
We start with the EM current for a singly heavy baryon, which is defined by where ψ(x) stands for the light-quark field ψ = (u, d, s) in SU(3) flavor space and Ψ denotes the heavy-quark field for the charmed or bottom quark. The charge operator Q is expressed aŝ The e Q in the second term in Eq. (1) denotes the charge corresponding to a heavy quark, which has the value 2/3 for the charm quark and −1/3 for the bottom quark. The matrix element of J µ between baryons with spin 3/2 can be parametrized in terms of four different real form factors as follows: where M B denotes the mass of a singly heavy baryon in the baryon sextet with spin 3/2. The metric tensor η αβ of Minkowski space is defined as η αβ = diag(1, −1, −1, −1). q α represents the momentum transfer q α = p α − p α and its square is written as q 2 = −Q 2 with Q 2 > 0. u α (p, s) means the Rarita-Schwinger spinor for a singly heavy baryon with spin 3/2, carrying the momentum p and the spin component s projected along the direction of the momentum. σ µν designates the antisymmetric tensor σ µν = i[γ µ , γ ν ]/2. Note that when one takes the limit of the infinitely heavy quark mass (m Q → ∞), the heavy-quark current given in the second part of Eq. (1) can be safely neglected for the EM form factors. It gives only a constant contribution to the electric form factors as already shown in Ref. [12]. It is more convenient to introduce the Sachs-type form factors or the multipole EM form factors, in particular, when the EM structure of a baryon with spin 3/2 is examined. The electric quadrupole form factor reveals how the shape of a baryon with spin 3/2 is deviated from the rotationally symmetric one. The Sachs-type form factors can be expressed in terms of F B i given in Eq. (3) have been already investigated in Ref. [12], we will focus on the electric quadrupole form factors of the baryon sextet with spin 3/2, i.e., G B E2 in the present work. At Q 2 = 0, G E2 (0) yields the electric quadrupole moment which determines how a baryon with spin 3/2 undergoes deformation from a spherical shape. If Q B has a negative value (Q B < 0), then the baryon takes a cushion shape, whereas if Q B is positive (Q B > 0), then it looks like a rugby-ball shape. We want to mention that the M 3 form factors vanish in the present work. In fact, any chiral solitonic approaches yield the null results of the M 3 form factors because of the hedgehog structure [24]. However, the experimental data on M 3 is absent to date and its value should be very tiny even if it is measured. In fact, one could compute the M 3 form factors if one goes beyond 1/N c corrections. This means that the M 3 form factors should be suppressed in the large N c limit. We will focus in the present work on the E2 form factors of the baryon sextet with spin 3/2. The SU(3) χQSM is constructed based on the following low-energy effective partition function in Euclidean space, defined by where ψ and U denote respectively the quark and pseudo-Nambu-Goldstone boson fields. Having integrated over quark fields, we can express the partition function in terms of the effective chiral action S eff , which is defined by where Tr represents the functional trace running over spacetime and all relevant internal spaces. The N c denotes the number of colors. M is the dynamical quark mass that arises from spontaneous symmetry breaking of chiral symmetry. U γ5 represents the chiral field that consists of the pseudo-Nambu-Goldstone (pNG) fields π a , a = 1, · · · 8, which is expressed as with We assume isospin symmetry, i.e., m u = m d . The average mass of the up and down quarks is defined by m = (m u + m d )/2. Then, the matrix of the current quark masses is written asm = diag(m, m, m s ) = m + δm. δm is written as where m 1 and m 8 denote the singlet and octet components of the current quark masses, defined by m 1 = (−m+m s )/3 and Since the pion field has flavor indices, one has to combine a minimal symmetric ansatz will be the hedgehog ansatz with which the flavor indices can be coupled to three-dimensional spatial axes. The pion fields are then expressed in terms of a single function P (r), which is called the profile function, as follows π a (x) = n a P (r) (12) with n a = x a /r. Then the SU(2) chiral field is written as with U SU(2) = exp(in · τ P (r)). The SU (3) which preserves the hedgehog ansatz. Integration over U in Eq. (6) quantizes the pNG fields. In the large N c limit, the meson mean-field approximation is justified [6,11]. Thus, we can carry out the integration over U in Eq. (6) around the saddle point, where δS eff /δP (r) = 0 is satisfied. This saddle-point approximation yields the equation of motion that can be solved self-consistently. The solution provides the self-consistent profile function P c (r) of the chiral soliton. A detailed method of the self-consistent procedure can be found in Ref. [17].
While the quantum fluctuations of the self-consistent pion fields can be ignored by the large N c argument, the fluctuations along the direction of both the rotational and translational zero modes cannot be ignored, since they are not at all small. Note that rotational and translational zero modes are related to rotational and translational symmetries. Thus, the zero modes can be taken into account by the following rotational and translational transformations where A(t) is an SU(3) unitary matrix. So, the functional integral over U can be approximated by those over zero modes The integration over translational zero modes will naturally give the Fourier transform of the EM densities. We refer to Ref. [18] for a detailed description of the zero-mode quantization in the present scheme.
Having carried out the zero-mode quantization, we obtain the collective Hamiltonian as where I 1 and I 2 denote the moments of inertia for the soliton. The parameters α, β, and γ for heavy baryons arise from the breaking of flavor SU(3) symmetry, which are defined by where K 1, 2 are the anomalous moments of inertia. Note that the number of light valence quarks for a singly heavy baryon is N c − 1. This means that the expression for the valence part of Σ πN contains also N c − 1 in place of N c . It can be related to the πN sigma term as follows: Σ πN = (N c − 1)N −1 c Σ πN . The detailed expressions for the moments of inertia and Σ πN are given in Ref. [20].
The presence of the symmetry-breaking part in the collective Hamiltonian, H sb , have the baryon wavefunctions mixed with those in higher SU (3) representations. In the present case, the collective wavefunctions for the baryon antitriplet (J = 0) and the sextet (J = 1) are obtained respectively as [20] with the mixing coefficients in the basis [Λ Q , Ξ Q ] for the antitriplet and Σ Q , Ξ Q , Ω Q for the sextets with both spin 1/2 and 3/2. The parameters p 15 , q 15 , and q 24 are explicitly written as The collective wavefunction for the soliton with (N c − 1) valence quarks is then obtained in terms of the SU (3) Wigner D functions where dim(ν) represents the dimension of the representation ν and Q S a charge corresponding to the soliton state S, i.e., Q S = J 3 + Y /2. F and S stand for the flavor and spin quantum numbers corresponding to the soliton for the singly heavy baryon. Finally, the complete wavefunction for a singly heavy baryon can be derived by coupling the soliton wavefunction to the heavy quark spinor where χ J Q3 denote the Pauli spinors for the heavy quark and C J J 3 J,J3 J Q J Q3 the Clebsch-Gordan coefficients. The matrix elements of the EM current (3) can be computed within the χQSM by representing them in terms of the functional integral in Euclidean space, where the baryon states |B, p and B, p | are, respectively, defined by The heavy baryon current J B can be constructed from the N c − 1 valence quarks where α 1 · · · α Nc−1 represent spin-flavor indices and i 1 · · · i Nc−1 color indices. The matrices Γ are taken to consider the quantum numbers JJ 3 T T 3 Y of the N c − 1 soliton. The creation operator J † B can be constructed in a similar way. The calculation of the baryonic correlation function given in Eq. (25) is a tedious one, so we will present here only the final expressions for the E2 form factor, As for the detailed formalism, we refer to Refs. [17,18].
The final expressions for the electric quadrupole form factors of the baryon sextet with spin 3/2 can be written as where j 2 (|q||z|) stands for the spherical Bessel function with order 2 and the corresponding density of the E2 form factors is given as The densities of E2 form factors I 1E2 and K 1E2 can be found in Appendix A. In the limit of m Q → ∞, the charge distribution of the heavy quark becomes a point-like static charge given as ρ Q (r) = e Q δ (3) (r). This leads to Q ij = d 3 rρ Q (r)(3r i r j − r 2 δ ij ) = 0. This implies that the E2 form factors of the singly heavy baryons are solely governed by the light quarks in the m Q → ∞ limit.
Having calculated the matrix elements of the collective operators in Eq. (29), we arrive at the final expressions for the E2 form factors of the baryon sextet with spin 3/2 where G , and G

B(wf) E2
denote respectively the symmetric terms, the flavor SU(3) symmetry-breaking ones from the effective chiral action, and those from the mixed collective wavefunctions, expressed explicitly as where Q B stands for the charge of the light-quark components of the corresponding baryons. We can derive similar sum rules for the electric quadrupole moments of singly heavy baryons with spin 3/2 as follows [24] B∈sextet Q B = 0, Even though the flavor SU(3) symmetry is broken, we still can find the following sum rules

III. RESULTS AND DISCUSSION
In the χQSM, there are several parameters to fix. Since the sea-quark or Dirac-sea contributions contain divergent integrals, one has to introduce a regularization to tame the divergences. In the present work, we introduce the propertime regularizations with the cutoff mass Λ. This can be fixed by using the pion decay constant f π = 93 MeV. The average mass of the up and down current quarks m is determined by the physical pion mass m π = 140 MeV. While the mass of the strange current quark m s can be also fixed by reproducing the kaon mass, which gives m s = 150 MeV, we preferably use m s = 180 MeV, since this value of m s yields the best results for the hyperon mass splittings [17,25]. The remaining parameter is the dynamical quark mass M , which is the only free parameter of the model. However, M = 420 MeV is known to be the best value in reproducing various observables in the light baryon sector [17]. Thus, we will use this value also in the present calculation.
It was shown that in the calculation of the E2 form factors of the baryon decuplet the sea-quark contributions turn out to be rather important, we will first examine the valence-and sea-quark contributions separately. In Fig. 1, we draw the numerical results for the E2 form factors of the baryon sextet with spin 3/2. As expected, the general behaviors of the valence-and sea-quark contributions to the E2 form factors of the heavy singly baryons are rather similar to those of the baryon decuplet. As shown in Fig. 1, the valence-quark contributions decrease off mildly as Q 2 increases, whereas the sea-quark or Dirac-sea contributions fall off drastically in the smaller Q 2 region, so that they govern the Q 2 dependence of the E2 form factors. In particular, the magnitudes of the sea-quark contributions are quite larger than in the region of smaller Q 2 . Thus, they are the main contributions to the electric quadrupole moments of the baryon sextet with spin 3/2. Considering the fact that the electric quadrupole moment exhibits how the corresponding baryon is deformed, the present results provide certain physical implications. Recent investigations on the gravitational form factors of baryons within the χQSM indicate that the valence quarks are mainly located in the inner part of a baryon, while the sea quarks lie in its outer part [27,28]. Thus, the sea-quark contributions, which can be also interpreted as pion clouds, mainly describe how a singly heavy baryon with spin 3/2 is deformed. The present results are in line with what was discussed in Ref. [14].
In Fig. 2 (Q 2 ) arises from the current-quark mass term in the effective chiral action given in Eq. (7), whereas the second one comes from the wavefunction corrections (20). each correction affects E2 form factors in a different way, as shown in Fig. 3. The wavefunction corrections to the E2 form factor of Σ * ++ c are negligibly tiny and the corrections from the current-quark mass term is also small. As a result, the m s corrections turn out to be negligible, as shown in the upper left panel of Fig. 2. On the other hand, the wavefunction corrections contribute noticeably to the E2 form factors of Σ * + c , while those from the current-quark mass term are of the same order as in the case of Σ * ++ c . In the case of Σ * 0 c and Ξ * 0 c , the wavefunction corrections to G are even larger than those from the mass term. This can be understood by examining Eqs. (32) and (33).
In the left panel of Fig. 4, we compare the results for the E2 form factors of the Ω * 0 c baryon with that from the lattice calculation. We employ for this comparison the unphysical pion mass m π = 156 MeV that is used in the lattice calculation. Note that there is only one lattice data with large uncertainty. The present results are still larger than the lattice one even with the larger pion mass used. We anticipate more accurate lattice data in the near future, so that one can draw a clear conclusion. In the right panel of Fig. 4, we depict the results of G Ω * 0 c E2 as a function of the pion mass m π with Q 2 = 0.183 GeV 2 fixed. As expected, the present results fall off slowly as m π increases.  For completeness, we present the results for the electric quadrupole moments of the baryon sextet with spin 3/2. Table I lists those of the Q B in the second and third rows, which correspond to the SU(3) symmetric and breaking cases, respectively. As already shown in Fig. 2, those of the charged baryon sextet have negative values of Q B , which indicates that the positively charged singly heavy baryons with spin 3/2 take oblate shapes. On the other hand, those of the neutral ones get positive values, so they are distorted in prolate forms. It is interesting to see that the Q B of the doubly positive-charged Σ * c is approximately 8 times larger than that of the singly positive-charged one. This can be understood by examining Eqs. (30).

IV. SUMMARY AND CONCLUSION
In the present work, we have investigated the electric quadrupole form factors of the lowest-lying singly heavy baryons with spin 3/2 in a pion mean-field approach, also knwon as the SU(3) chiral quark-soliton model. In the limit of an infinitely heavy quark, a heavy quark inside a singly heavy baryon can be regarded as a mere static color source. This means that the N c − 1 light valence quarks govern the quark dynamics inside a heavy baryon. The presence of the N c − 1 light valence quarks make the vacuum polarized, which produces the pion mean fields. The N c − 1 valence quarks are bound by the attraction provided by the pion mean fields self-consistently, from which a colored soliton consisting of the N c − 1 valence quarks arises. We call this soliton as a N c − 1 soliton. Then the singly heavy baryon can be constructed by coupling the N c − 1 soliton with a heavy quark. This is called the pion mean-field approach for the singly heavy baryons. Based on this pion mean-field approach, we computed the electric qudrupole form factors of the baryon sextet with spin 3/2, taking into account the rotational 1/N c and linear m s corrections.   We first examined the valence-and sea-quark contributions separately. As in the case of the baryon decuplet, the contributions from the sea quarks or the Dirac-sea level quarks govern the electric quadrupole form factors, in particular, in the smaller Q 2 region. Considering the fact that the electric quadrupole moment of a baryon provides information on how the baryon is deformed, we can draw the following physical implications: the deformation of a singly heavy baryon is also mainly governed by the sea-quark contributions or the pion cloud effects. We found a similar feature in the case of the baryon decuplet. The effects of the explicit flavor SU (3)