\section{Introduction}
The electromagnetic form factors are the fundamental quantities of
theoretical and experimental interest to investigate the internal
structure of nucleon. The knowledge of internal structure of
nucleon in terms of quarks and gluons degrees of freedom of
Quantum Chromodynamics (QCD) provide a basis for understanding
more complex, strongly interacting matter.  While QCD is accepted
as the fundamental theory of strong interactions, it cannot be
solved accurately in the nonperturbative regime. A coherent
understanding of the hadron structure in this energy regime is
necessary to describe the strong interactions as they are
sensitive to the pion cloud and provide a test for the QCD
inspired effective field theories based on the chiral symmetry.
Recently, a wide variety of accurately measured data have been
accumulated for the static properties of baryons, for example,
masses, electromagnetic moments, charge radii, and low energy
dynamical properties such as scattering lengths and decay rates
etc.. The charge radii and quadrupole moments are important
observables in hadronic physics as they lie in the nonperturbative
regime of QCD and give valuable information on the internal
structure of hadrons.

The Naive Quark Model (NQM) is unable to explain the magnitude and
sign of deformation measured for the case of nucleon and
$\Delta^+$.  A promising approach offered to solve QCD in this
energy regime are the constituent quark models which can be
constructed so as to include the relevant properties of QCD coming
from the consequences of the spontaneous breaking of chiral
symmetry ($\chi$SB). One of the important nonperturbative
approaches in this energy regime is chiral constituent quark model
($\chi$CQM) \cite{wein}. The $\chi$CQM coupled with the ``quark
sea'' generation through the chiral fluctuation of a constituent
quark into a Goldstone bosons (GBs) \cite{cheng,johan,song},
successfully explains the ``proton spin crisis'' \cite{hd},
hyperon $\beta$ decay parameters \cite{nsweak}, the octet and
decuplet baryons magnetic moments \cite{hdmagnetic}. The extension
to the SU(4) symmetry successfully predicts the contribution of
{\it intrinsic} charm (IC) content in the low lying and charmed
baryon magnetic moments and their radiative decay widths
\cite{hdcharm}. In this context, it become desirable to extend the
model to other low energy properties like charge radii, quadrupole
moment, and higher order moments of the multipole expansion.

The purpose of present communication is to calculate the {\it
intrinsic} quadrupole moment of the octet and decuplet baryons
within the framework of $\chi$CQM using the general
parameterization method. In order to understand the important role
played by the pion cloud and SU(3) symmetry breaking in measuring
the quadrupole moment we would carry out the calculations with and
without symmetry breaking. The results have also been compared
with the NQM predictions and the latest available data. Further,
it would also be interesting to understand in detail the role of
GP model parameters in the determination in of baryon quadrupole
moment.

\section{Intrinsic quadrupole moment of the nucleon}
The mean square charge radii $(r^2_B)$ and quadrupole moments
(${Q}_B$) are the lowest order moments of the charge density
operator $\rho$ in a low-momentum expansion. For example, for any
baryon $|B\rangle$ with charge $e_B$, the terms up to order of
$q^2$ for the charge density are \be \langle {\mathrm B}| \rho(q)
|B \rangle = e_B - \frac{q^2}{6} r^2_B - \frac{q^2}{6} {Q}_B +
..... \,,\ee where $q$ is the photon momentum transferred to the
baryon. The first two terms arise due to the spherically symmetric
monopole part of charge density, while the third term arises due
to the quadrupole part of the charge density. They characterize
the total charge, spatial extension of charge radii, and shape of
the system, respectively.

The shape of a spatially extended particle is determined by its
quadrupole moment \cite{buch01,buch04,buch07}, corresponding to
the charge quadrupole form factor $G_{C2}(q^2)$ at zero momentum
transfer. The {\it intrinsic} quadrupole moment of a nucleus with
respect to the body frame of axis is defined as \be Q_0=\int
d^3r~\rho({\bf r}) (3 z^2 - r^2)\,. \ee If the charge density is
concentrated along the $z$-direction (symmetry axis of the
particle), the term proportional to $3z^2$ dominates, $Q_0$ is
positive, and the particle is prolate shaped. If the charge
density is concentrated in the equatorial plane perpendicular to
$z$ axis, the term proportional to $r^2$ prevails, $Q_0$ is
negative and the particle is oblate in shape. The angular momentum
selection rule however, does not allow the spin $\frac{1}{2}^+$
baryons to have any spectroscopic quadrupole moment. Therefore,
the $\gamma^* N_{\frac{1}{2}^+} \to \Delta^+_{\frac{3}{2}^+}$
transition is studied to understand the quadrupole amplitudes in
the nucleon as well as $\Delta^+$.

The spin and parity conservation in the $\Delta^+ p$ transition
require three contributing amplitudes, magnetic dipole $M1$, the
electric quadrupole moment $E2$, and the charge quadrupole moment
$C2$ photon absorption amplitudes. The information on the {\it
intrinsic} quadrupole moments can be obtained from the
measurements of $E2$ and $C2$ amplitudes \cite{tia03,ber03}. If
the charge distribution of the initial and final three-quark
states were spherically symmetric, the $E2$ and $C2$ amplitudes of
the multipole expansion would be zero \cite{becchi}. However, the
recent experiments at JLAB, SELEX Collaboration reveal that
although these quadrupole amplitudes are small compared to the
dominant magnetic dipole transition $M1$, they are clearly non
zero \cite{pdg}.  More recently, the quadrupole transition moment
($Q_{\Delta^+ N}$) measured by LEGS and Mainz Collaborations is
$-0.108 \pm 0.009 \pm 0.034 ~{\rm fm}^2$ \cite{blanquad} and
$-0.0846 \pm 0.0033 ~{\rm fm}^2$ \cite{tia03}, respectively. These
measurements lead to the conclusion that the nucleon and the
$\Delta^+$ are intrinsically deformed.

\section{General Parameterization method}

In order to predict the sign as well as magnitude of deformation
in the octet and decuplet baryons, we have used the general
parameterization (GP) method \cite{morp89}. The charge quadrupole
operator composed of a two- and three-body operator terms in
spin-flavor space is given as \bea {Q} & = & {\mathrm B} \sum_{i
\ne j}^3 e_i \left( 3 \sigma_{i \, z} \sigma_{ j \,z} -
\bf{\sigma}_i \cdot \bf{\sigma}_j \right) + {\mathrm C} \!\!
\sum_{i \ne j \ne k }^3 e_i \left( 3 \sigma_{j \, z} \sigma_{ k \,
z} - \bf{\sigma}_j \cdot \bf{\sigma}_k \right) \,, \label{quad}
\eea where the coefficients ${\mathrm B}$ and ${\mathrm C} $ are
the constants to be determined from the experimental observations
on charge radii and quadrupole moments. The quadrupole moments
${Q}$ for the octet and decuplet baryons can be calculated from
Eq. (\ref{quad}) by evaluating matrix elements of the operator
corresponding to the three-quark spin-flavor wave functions  ($
Q_{B} = \left \langle {B} \vert {Q}{\vert B} \right \rangle)$. It
is straightforward to verify that \bea \sum_{i \neq j} e_i (
\xbf{\sigma_i} \cdot \xbf{\sigma_j}) &=& 2 {\xbf J} \cdot \sum_{i}
e_i \xbf{\sigma_i} - 3 \sum_{i} e_i \,, \label{eisisj} \eea \bea
\sum_{i \neq j \neq k} e_i(\xbf{\sigma_j} \cdot \xbf{\sigma_k})
&=& \pm 3 \sum_i e_i - \sum_{i \neq j} e_i( \xbf{\sigma_i} \cdot
\xbf{\sigma_j}) \,.\label{eisjsk} \eea In Eq. (\ref{eisjsk}), $+$
sign holds for $J=\frac{3}{2}$ and $-$ sign for $J = \frac{1}{2}$
states. Using the expectation value of operator $2 J \cdot
\sum_{i} e_i \xbf{\sigma_i}$ between the baryon wavefunctions $|B
\rangle$ in the initial and final states \cite{yaoubook}, the
operators in Eqs. (\ref{eisisj}) and (\ref{eisjsk}) become \bea
\ba{c|c|c} {\rm Operator} & \sum_{{i \neq j}} e_i(\xbf{\sigma_i}
\cdot \xbf{\sigma_j}) & \sum_{i\neq j \neq k} e_i(\xbf{\sigma_j}
\cdot \xbf{\sigma_k}) \\ \hline J = \frac{1}{2} & 3 \sum_{i} e_i
\xbf{\sigma_{iz}} - 3 \sum_{i} e_i & - 3 \sum_{i} e_i
\xbf{\sigma_{iz}}\\ J = \frac{3}{2} & 5 \sum_{i} e_i
\xbf{\sigma_{iz}} - 3 \sum_{i} e_i  & 6 \sum_i e_i - 5 \sum_i e_i
\xbf{\sigma_{iz}}  \\  \ea \label{tensor}\eea

The expression for the quadrupole moment of the octet and decuplet
baryons in Eq. (\ref{quad}) can be expressed as \bea {{Q}}_{1/2}
&=& 3{\mathrm B}\sum_{i \neq j} e_i \bf{\sigma_{iz}}
\bf{\sigma_{jz}} + 3 {\mathrm C}\sum_{i \ne j \ne k } e_i
\bf{\sigma_{jz}} \bf{\sigma_{kz}} +(- 3{\mathrm B} +3 {\mathrm C}
) \sum_{i} e_i \bf{\sigma_{iz}} +3 {\mathrm B}\sum_{i} e_i \,,
\label{q1/2} \\ {{ Q}}_{3/2} &=& 3 {\mathrm B} \sum_{i \neq j} e_i
\bf{\sigma_{iz}} \bf{\sigma_{jz}} + 3{\mathrm C} \sum_{i \ne j \ne
k } e_i \bf{\sigma_{jz}} \bf{\sigma_{kz}} + (-5 {\mathrm B} + 5
{\mathrm C})\sum_{i} e_i \bf{\sigma_{iz}} + (3 {\mathrm B} -6
{\mathrm C}) \sum_{i} e_i \,,\label{q3/2} \eea where $i = (u,d,s)$
for any of the three quarks. Therefore, the calculation of
quadrupole moment reduces to the calculation of the flavor
structure $(\sum_i e_{i})$, spin structure $(\sum_i e_i
\sigma_{iz})$ and the tensor terms $(\sum_i e_i
\sigma_{iz}\sigma_{jz}$ and $\sum_i e_i \sigma_{jz} \sigma_{kz})$
for a given baryon.

The spin and flavor structure of a given baryon can be calculated
using the ${\rm SU}(6)$ spin-flavor symmetry of the wave
functions. The expectation value of the appropriate operators are
expressed as $\widehat e_{i} \equiv \langle B| \sum_i e_{i}|B
\rangle \,,$ and $ \widehat {e_{i} \sigma_{iz}} \equiv \langle B|
\sum_i {e_i} {\sigma_{iz}} |B \rangle $, where $|B \rangle$ is the
baryon wave function and $e_{i}$ and $\sigma_{i}$ are the charge
and spin operators defined as \bea \sum_i e_i &=& \sum_{q=u,d,s}
n^B_{q}q + \sum_{{\bar q} = {\bar u}, {\bar d}, {\bar s}}
n^B_{{\bar q}} {\bar q} = n^B_{u}u + n^B_{d}d + n^B_{s}s +
n^B_{\bar u}{\bar u} + n^B_{ \bar d}{\bar d} + n^B_{\bar s}{\bar
s} \label{numei} \,, \\ \sum_i {e_i} {\sigma_{iz}} &=&
\sum_{q=u,d,s} (n^B_{q_{+}}q_{+} + n^B_{q_{-}}q_{-}) =
n^B_{u_{+}}u_{+} + n^B_{u_{-}}u_{-} + n^B_{d_{+}}d_{+} +
n^B_{d_{-}}d_{-} + n^B_{s_{+}}s_{+} + n^B_{s_{-}}s_{-}
\label{numeisi} \,. \eea Here $n_q^B$$(n_{\bar q}^B)$ is the
number of quarks with charge $q$(${\bar q}$), and
$n^B_{q_{+}}$($n^B_{q_{-}}$) is the number of $q$ quarks with spin
$+$($-$). The tensor terms can be simplified and further reduced
to the calculations of spin and flavor structure as presented in
Eq. (\ref{tensor}).

For the case of octet baryons, Eq. (\ref{q1/2}) can be solved for
the spin-flavor symmetric SU(6) octet baryon wave function by
using the operators  defined in Eqs. (\ref{numei}) and
(\ref{numeisi}). For ready reference, the quadrupole moment of
proton, neutron and $\Sigma^+$ are expressed as \bea {{Q}}_p &=&
3{\mathrm B} \left(- 2u_+ - d_+ + 2 u + d \right) + {\mathrm C}
\left(-4u + d + 4u_+ - d_+ \right) \,, \label{qp}
\\ {{Q}}_n &=& 3 {\mathrm B} \left( -u_+ - 2d_+ + u + 2d \right) +
{\mathrm C} \left( u - 4 d - u_+ + 4d_+ \right) \,, \label{qn} \\
{{Q}}_{\Sigma^+} &=& 3 {\mathrm B} \left( -2 u_+ - s_+ + 2u + s
\right) + {\mathrm C} \left( -4u + s + 4u_+ - s_+ \right) \,.
\label{qs} \eea Similarly, the quadrupole moment of the decuplet
baryon $\Delta^+$ can be calculated from Eq. ({\ref{q3/2}}) by
solving the matrix elements corresponding to the decuplet baryon
wave function and by using the operators defined in Eqs.
({\ref{numei}) and ({\ref{numeisi}). We have \bea {{
Q}}_{{\Delta^+}} &=& {\mathrm B}\left( 2u_+ + d_+ + 3(2 u +
d)\right) + {\mathrm C}\left( -3(2u + d) + 5(2 u_+ + d_+ )\right)
\,. \label{qd} \eea

In the case of naive quark model (NQM), the calculations lead to
$u= \frac{2}{3}$, $d = -\frac{1}{3}$, $s = -\frac{1}{3}$, $u_+ =
\frac{2}{3}$, $d_+ =- \frac{1}{3}$, and $s_+ = -\frac{1}{3}$. The
quadrupole moments of $p$, $n$, $\Sigma^+$, and $\Delta^+$ are now
given as \be {{Q}}_p =0\,, ~~~ {{Q}}_n =0 \,, ~~~{{Q}}_{\Sigma^+}
=0 \,, ~~~ {{Q}}_{{\Delta^+}} = 4{\mathrm B} + 2{\mathrm C}\,.\ee
These results suggest that the octet baryons are spherically
symmetric and there is no deformation. However, the experimental
results suggest that the nucleon is intrinsically deformed.
Therefore, there is a need to go beyond this model to include the
effects of ``quark sea''.

\section{Chiral constituent quark model}
The basic process in the $\chi$CQM is the emission of a GB by a
constituent quark which further splits into a $q \bar q$ pair as $
q_{\pm} \rightarrow {\rm GB}^{0} + q^{'}_{\mp} \rightarrow (q \bar
q^{'}) +q_{\mp}^{'}\,, \label{basic}$ where $q \bar q^{'} +q^{'}$
constitute the ``quark sea'' \cite{cheng,johan,song}. The
effective Lagrangian describing interaction between quarks and a
nonet of GBs is ${\cal L} = g_8 \bar q \Phi' q \,,$ where $g_8$ is
the coupling constant for the GBs. The GB field $\Phi'$ is
expressed in terms of the quark contents as \bea {\Phi'} = \left(
\ba{ccc} \phi_{uu} u \bar u+ \phi_{ud} d \bar d +\phi_{us} s \bar
s& \varphi_{ud} u \bar d & \varphi_{us} u \bar s
\\ \varphi_{du} d \bar u & \phi_{du}u \bar u+ \phi_{dd} d \bar d
+\phi_{ds} s \bar s & \varphi_{ds} d \bar s
\\ \varphi_{su} s \bar u & \phi_{sd} s \bar d & \phi_{su} u \bar u
+ \phi_{sd} d \bar d +\phi_{ss} s \bar s \\ \ea \right)\,, \eea
where \bea \phi_{uu} &=& \phi_{dd}= \frac{1}{2} +\frac{\beta}{6} +
\frac{\zeta}{3}\,,~~~~~ \phi_{ss} =\frac{2\beta}{3} +
\frac{\zeta}{3}\,,~~~~~ \phi_{us} = \phi_{ds}=
\phi_{su}=\phi_{sd}= -\frac{\beta}{3} + \frac{\zeta}{3}\,,
\nonumber \\ \phi_{du} &=& \phi_{ud}= -\frac{1}{2}
+\frac{\beta}{6} + \frac{\zeta}{3}\,,~~~~~ \varphi_{ud} =
\varphi_{du} = 1\,, ~~~~\varphi_{us} = \varphi_{ds}=\varphi_{su} =
\varphi_{sd}= \alpha \,. \label{phi1} \eea SU(3) symmetry breaking
is introduced by considering $M_s > M_{u,d}$ as well as by
considering the masses of GBs to be nondegenerate $(M_{K, \eta} >
M_{\pi}$ and $M_{\eta^{'}} > M_{K, \eta})$
{\cite{cheng,johan,song,hd}. The parameter $a(=|g_8|^2$) denotes
the transition probability of chiral fluctuation of the splitting
$u(d) \rightarrow d(u) + \pi^{+(-)}$, whereas $\alpha^2 a$,
$\beta^2 a$ and $\zeta^2 a$ respectively, denote the probabilities
of transitions of~ $u(d) \rightarrow s  + K^{-(o)}$, $u(d,s)
\rightarrow u(d,s) + \eta$, and $u(d,s) \rightarrow u(d,s) +
\eta^{'}$.

Since the quadrupole moment operators for the spin $\frac{1}{2}$
and spin $\frac{3}{2}$ baryons involve the knowledge of spin and
flavor structure of baryons, it is important to mention here that
redistribution of flavor and spin take place among the ``quark
sea''. The modified flavor and spin structure of the baryon in
$\chi$CQM due to the chiral symmetry breaking can be calculated by
substituting for every constituent quark \bea q \to P_q q +
|\psi(q)|^2  ~~~~~~~~ {\rm and} ~~~~~~~~q_{\pm} \to P_q q_{\pm}+
|\psi(q_{\pm})|^2 \,, \label{spin} \eea where $P_q = 1 - \sum P_q$
is the transition probability of no emission of GB from any of the
$q$ quark, $|\psi(q)|^2$ is the transition probability of the $q$
quark, and $|\psi(q_{\pm})|^2$ is the probability of transforming
a $q_{\pm}$ quark \cite{hd}.

\section{Results and Discussion}

In order to calculate the quadrupole moments of the octet and
decuplet baryons in the $\chi$CQM, we substitute Eq. (\ref{spin})
in Eqs. (\ref{q1/2}) and (\ref{q3/2}) . For example, on
substituting Eq. (\ref{spin}) in Eqs. (\ref{qp}), (\ref{qn}) and
(\ref{qs}), the quadrupole moments of $p$, $n$, and $\Sigma^+$ in
$\chi$CQM can now be expressed as \bea {{Q}}_p &=& {\mathrm B} a
\left(6 + \beta^2 + 2 \zeta^2 \right) - {\mathrm C} a \left(4 + 2
\alpha^2 + \beta^2 + 2 \zeta^2 \right) \,, \label{qp1}
\\ {{Q}}_n &=& {\mathrm B} a (3 - 3 \alpha^2) + {\mathrm C}
\frac{a}{3}(3 + 9 \alpha^2 + 2 \beta^2 + 4 \zeta^2) \,,
\label{qn1} \\ {Q}_{\Sigma^+}  &=& {\mathrm B}a (6 + \alpha^2 + 2
\zeta^2) -{\mathrm C} \frac{a}{3} (5 \alpha^2 + 4 \beta^2 + 12 + 6
\zeta^2) \,.\label{qs1} \eea From the above equations, we can
directly estimate the effects of SU(3) symmetry breaking and pion
cloud in the quadrupole moments of the octet baryons. It is clear
from the equations that, for the non zero value of the GP
parameters ($\mathrm B$ and $\mathrm C$), there is a significant
contribution of the $\chi$CQM parameters ($a, \alpha, \beta,
\zeta$) to the quadrupole moments. Similarly, the quadrupole
moment of $\Delta^+$ and $\Xi^{*0}$ in $\chi$CQM, after
substituting the contribution coming from the ``quark sea'', can
be expressed as \bea {{ Q}}_{{\Delta^+}} &=& {\mathrm B}\left(4 -
\frac{1}{3} a \left(6 + \beta^2 + 2 \zeta^2 \right) \right) +
{\mathrm C} \left(2 - \frac{5}{3} a \left(6 + \beta^2 + 2\zeta^2
\right) \right) \,, \label{qd1} \\ {Q}_{\Xi^{*0}} &=& {\mathrm B}
\frac{a}{3} (-3 + \alpha^2 + 2 \beta^2) + {\mathrm C} \frac{5a}{3}
(-3 + \alpha^2 + 2 \beta^2)\,. \eea

Before giving the numerical results, we would like to discuss the
input parameters involved in the calculations of baryon quadrupole
moments in $\chi$CQM. The calculations involve the two set of
parameters, the ones corresponding to the GP method (${\mathrm B}$
and ${\mathrm C}$, where ${\mathrm B}> {\mathrm C}$), and other
corresponding to the $\chi$CQM symmetry breaking ($a$, $\alpha$,
$\beta$, and $\zeta$). In order to fix the values of GP
parameters, we have performed a fit to the available experimental
values of the baryon charge radii leading to $ {\mathrm B} =
-0.0525 \,, ~~ {\mathrm C} = -0.0158 \,,$ as the best fit. For the
$\chi$CQM parameters, we have used the same set of parameters as
discussed in our earlier publication \cite{nsweak}. The values
used are $ a = 0.12\,, ~~~ \alpha = 0.45\,, ~~~ \beta = 0.45\,,
~~~ \zeta = -0.15$.

Using the above discussed set of parameters, we have calculated
the quadrupole moments of octet and decuplet baryons in $\chi$CQM
and results have been presented in Tables \ref{quad1/2} and
\ref{quad3/2}. In order to understand the role of chiral symmetry
breaking and SU(3) symmetry breaking, we have also presented the
results of NQM in the tables.  Most of the models in literature
are unable to estimate the {\it intrinsic} quadrupole moment of
the baryons. In the case of octet baryons there are indirect
evidences of small deformation in the nucleon. This deformation
can easily be observed in our results. For the case of decuplet
baryons also, $\chi$CQM is able to give a fairly good prediction
of the sign as well as magnitude of the quadrupole moments.

In order to understand the effect of three-quark contributing term
and to make our calculations more responsive, we have also
presented the results by neglecting this contribution which can be
obtained by fixing the coefficient $\mathrm C =0$. From Table
\ref{quad1/2}, it is clear that the results are affected to a very
small extent by the inclusion of three-quark term. The three-quark
terms do not seem to play any significant role in the case of
octet baryons. However, in the case of decuplet baryons, the
results in Table \ref{quad3/2} reveal that the inclusion of
three-quark term increases the quadrupole moments thus making them
significant in this case. It is interesting to observe that the
effect of three-quark contribution is even more in the case of
neutral baryons. For example, in the case of charged baryons, the
predictions are increased about 10\% whereas for the case of
neutral particles the variation is more then 50\% making the
effect of three-quark term significant.

Since the electric quadrupole moment of the octet baryons as well
as that of octet-decuplet transitions are amenable to measurement,
any experimental information would have important implications for
the basic tenets of $\chi$CQM and the effects of SU(3) symmetry
breaking.

\begin{table}
\begin{center}
\begin{tabular}{|c|c|c|}
\hline
Baryon&\multicolumn{2}{c|}{$\chi$CQM} \\ \hline
&  ${\mathrm C}$ = 0 & ${\mathrm C}$ = $-$0.016 \\ \hline ${
Q}_{p}$ & $-$0.041& $-$0.032  \\ ${ Q}_{n}$ &$-$0.016 & $-$0.019
\\ ${ Q}_{\Sigma^{+}}$ &$-$0.041 & $-$0.032 \\ ${ Q}_{\Sigma^{-}}$
&0.010 & 0.009 \\ ${ Q}_{\Sigma^{0}}$ &$-$0.016& $-$0.012 \\ ${
Q}_{\Xi^{0}}$ & $-$0.016 & $-$ 0.019\\ ${ Q}_{\Xi^-}$ & 0.010&
0.009 \\
\hline
\end{tabular}
\end{center} \caption{Quadrupole moments of spin $\frac{1}{2}^+$
baryons in $\chi$CQM. The deformation is zero for all the cases in
NQM. }\label{quad1/2}
\end{table}

\begin{table}
\begin{center}
\begin{tabular}{|c|c|c|c|}
\hline Baryon &  NQM& \multicolumn{2}{c|}{$\chi$CQM}  \\ \hline &
& ${\mathrm C}$ = 0 &${\mathrm C}$ = $-$0.016 \\ \hline ${
Q}_{\Delta^{++}}$  & $-$0.484 & $-$0.398 &$-$0.428 \\
${Q}_{\Delta^{+}}$ & $-$0.242 &$-$0.196 &$-$0.208 \\ ${
Q}_{\Delta^{0}}$ & 0.0 & 0.005 & 0.013 \\ ${ Q}_{\Delta^{-}}$
&0.242 & 0.207 & 0.208 \\ ${Q}_{\Sigma^{*+}}$ &$-$0.242 &$-$0.196
&$-$0.208 \\ ${ Q}_{\Sigma^{*-}}$ &0.242 & 0.207 & 0.234\\ ${
Q}_{\Sigma^{*0}}$ &0.0 & 0.005 &0.013 \\
 ${Q}_{\Xi^{*0}}$ & 0.0& 0.005 & 0.013 \\
${ Q}_{\Xi^{*-}}$ &0.242 & 0.207 & 0.234\\ ${ Q}_{\Omega^-}$ &
0.315 & 0.296 &0.296 \\
  \hline
\end{tabular}
\end{center} \caption{Quadrupole moments of spin $\frac{3}{2}^+$ baryons in $\chi$CQM.}
\label{quad3/2}
\end{table}


\section{Introduction} \label{S:intro}

The definite integral
\begin{equation}
N_{0,4}(a;m) = \int_{0}^{\infty} \frac{dx}{(x^{4}+2ax^{2}+1)^{m+1}},
\end{equation}
\noindent
with $m \in \mathbb{N}$ and $a>-1$ has the value
\begin{equation}
N_{0,4}(a;m)  =   \frac{\pi}{2^{m+3/2} (a+1)^{m+1/2} }
P_{m}(a)  \label{identity-1}
\end{equation}
\noindent
where
\begin{equation}
P_{m}(a) = 2^{-2m} \sum_{k=0}^{m} 2^{k} \binom{2m-2k}{m-k} \binom{m+k}{k}
(a+1)^{k}.
\label{polyP-def}
\end{equation}

The reader will find in \cite{amram} a survey of the many proofs of
(\ref{identity-1}) that have appeared in the literature. After this survey
was written, other proofs have appeared; see
\cite{apagodu10,berg-vignat, gonzalez3,koutschan1}.  The goal of this note
is to present some additional proofs: the first one is based in
the {\em Schwinger  parametrization}, a second one using {\em Schur
functions} and the last one involves a representation involving  integrals
of {\em Bessel functions}.

\section{Schwinger parametrization} \label{S:Fey}

The method of Schwinger  parameters,  described in classical quantum
theory books such as \cite{peskin}, is utilized here to give a proof of
(\ref{identity-1}). A preliminary discussion is presented first.

Starting from \begin{equation}
\frac{1}{A_{1}A_{2}}=
\int_{0}^{1}\frac{dx}{\left(xA_{1}+\left(1-x\right)A_{2}\right)^{2}}
=\int_{0}^{1} \int_{0}^{1}\frac{\delta\left(1-x_{1}-x_{2}\right)}
{\left(x_{1}A_{1}+x_{2}A_{2}\right)^{2}}dx_{1}dx_{2}
\label{eq:A1A2}\end{equation}
\noindent
where $\delta$ is the Dirac measure, it
can be deduced by induction that
\begin{equation}
\frac{1}{\prod_{i=1}^{n}A_{i}}=
\left(n-1\right)!\int_{[0,1]^{n}}
\frac{\delta\left(1-\sum_{i=1}^{n}x_{i}\right)}
{\left(\sum_{i=1}^{n}x_{i}A_{i}\right)^{n}}\prod_{i=1}^{n}dx_{i}.
\label{eq:A1A2An}\end{equation}
The induction step uses the $n$-times differentiated version of \eqref{eq:A1A2},
namely
\begin{equation}
\frac{1}{A_{1}A_{2}^{\nu_{2}}}=
\int_{0}^{1}\int_{0}^{1}\frac{\delta\left(1-x_{1}-x_{2}\right)
\nu_{2}x_{2}^{\nu_{2}-1}}{\left(x_{1}A_{1}+x_{2}A_{2}\right)^{\nu_{2}+1}}
dx_{1}dx_{2}.
\end{equation}
Repeated differentiations of \eqref{eq:A1A2An} yields the more general
result\begin{equation}
\frac{1}{\prod_{i=1}^{n}A_{i}^{\nu_{i}}}=
\frac{\Gamma\left(\nu\right)}{\prod_{i=1}^{n}\Gamma\left(\nu_{i}\right)}
\int_{[0,1]^{n}}
\frac{\delta\left(1-\sum_{i=1}^{n}x_{i}\right)
\prod_{i=1}^{n}x_{i}^{\nu_{i}-1} \, dx_{i}}
{\left(\sum_{i=1}^{n}x_{i}A_{i}\right)^{\nu}}
\label{eq:general formula}
\end{equation}
with $\nu= \nu_{1} + \cdots + \nu_{n}$ and $A_{i} \in \mathbb{C}$.

An alternative proof has been provided by Shapiro
in his lecture notes \cite{shapiro}. This direct proof of \eqref{eq:A1A2An}
starts with
\begin{equation}
\frac{1}{A}=\int_{0}^{\infty}e^{-Ax}dx
\end{equation}
for $\mathop{\rm Re}\nolimits{A}>0$.  It follows that

\begin{equation}
\prod_{i=1}^{n}\frac{1}{A_{i}}=
\int_{\mathbb{R}_{+}^{n}}e^{-\sum_{i=1}^{n}A_{i}x_{i}}\prod_{i=1}^{n}dx_{i}
\end{equation}
Now let $x=x_{1}+ \cdots x_{n}$ and $y_{i}=\frac{x_{i}}{x}$ to obtain
\begin{equation}
\prod_{i=1}^{n} dx_{i}
=x^{n-1}dx\prod_{i=1}^{n}\delta\left(1-\sum_{i=1}^{n}y_{i}\right)
\, dy_{i}
\end{equation}
\noindent
so that
\begin{equation}
\prod_{i=1}^{n}\frac{1}{A_{i}}=
\int_{\mathbb{R}_{+}^{n}}\prod_{i=1}^{n}
\delta\left(1-\sum_{i=1}^{n}y_{i}\right)x^{n-1}dx\,
e^{-x\sum_{i=1}^{n}A_{i}y_{i}} \, dy_{i}.
\end{equation}
After a change of variable, the latter integral is evaluated as
\begin{equation}
\int_{0}^{\infty} x^{n-1} e^{-x \sum_{i=1}^{n} A_{i}y_{i}} dx =
\Gamma(n) \left( \sum_{i=1}^{n} A_{i}y_{i} \right)^{-n}
\end{equation}
\noindent
resulting in
\begin{equation}
\prod_{i=1}^{n}\frac{1}{A_{i}}
=\Gamma\left(n\right) \int_{\mathbb{R}_{+}^{n}}\prod_{i=1}^{n}
\frac{\delta\left(1-\sum_{i=1}^{n}y_{i}\right)}
{\left(\sum_{i=1}^{n}A_{i}y_{i}\right)^{n}} \, dy_{i}.
\end{equation}

Shapiro \cite{shapiro} states:
\begin{quotation}
An interesting anecdote of physics history is that Schwinger remained
bitter that a virtually identical mathematical trick became commonly
known as Feynman parameters. Why two brilliant physicists, each of
whom had an appropriately won a Nobel prize, should fight over what
is essentially a trivial mathematical trick, is an interesting question
in the sociology of physicists.
\end{quotation}\
An additional remark is that this representation is also known in
probability theory, and was apparently introduced by Mauldon \cite{mauldon}.
The $n$-dimensional Dirichlet distribution
with parameters $\left\{ \nu_{i}\right\} _{1\le i\le n}$
of a random vector $X$ reads\[
f_{X}\left(x_{1},\dots,x_{n}\right)=\frac{\Gamma\left(\nu\right)}{\prod_{i=1}^{n}\Gamma\left(\nu_{i}\right)}\prod_{i=1}^{n}x_{i}^{\nu_{i}-1}\delta\left(1-\sum_{i=1}^{n}x_{i}\right)\]
with $\nu = \nu_{1} + \ldots + \nu_{n}$. Mauldon
proved that a Dirichlet distributed random vector $X$ satisfies the following
equality
\begin{equation}
E_{X}\left(\sum_{i=1}^{n}\lambda_{i}X_{i}\right)^{-\nu}
=\sum_{i=1}^{n}\lambda_{i}^{-\nu_{i}},
\end{equation}
which is exactly \eqref{eq:general formula}.

\medskip

This procedure is now applied to the integral $N_{0,4}(a;m)$ written as
\begin{equation}
N_{0,4}(a;m) = \int_{0}^{\infty}
\frac{dx}{(1+ e^{i \theta} x^{2})^{m+1} \, \,
(1+ e^{-i \theta} x^{2})^{m+1}} \, dx
\end{equation}
\noindent
with $a = \cos \theta$ and $\theta \in (-\pi,\pi)$. The corresponding
parameters
are $n=2$ and $\nu_{1}=\nu_{2} = m+1$.

Recall the expression for the beta function
\begin{equation}
B(a,b)  = \int_{0}^{1} x^{a-1}(1-x)^{b-1} \, dx
\end{equation}
\noindent
and its relation to the gamma function
\begin{equation}
B(a,b) = \frac{\Gamma(a+b)}{\Gamma(a) \, \Gamma(b)}
\end{equation}
\noindent
to write
\begin{equation}
N_{0,4}(a;m) = \frac{1}{B(m+1,m+1)}
\int_{0}^{\infty} \, dx \int_{0}^{1}
\frac{x_{1}^{m} (1-x_{1})^{m} \, dx_{1}}
{ \left[ 1 + x^{2} \times  (x_{1}e^{i \theta} +
(1-x_{1})e^{-i \theta} )\right]^{2m+2}}.
\nonumber
\end{equation}
\noindent
Formula $3.194.3$ in \cite{gr} gives
\begin{equation}
\int_{0}^{\infty} \frac{dx}{(1+ux^{2})^{\alpha}} =
\frac{1}{2 \sqrt{u}} B \left( \tfrac{1}{2}, \alpha - \tfrac{1}{2} \right)
\end{equation}
\noindent
valid for  $u \in \mathbb{C}$ with
$|\text{Arg}(u) | < \pi$. This yields
\begin{equation}
N_{0,4}(a;m) = \frac{B(\tfrac{1}{2},2m+\tfrac{3}{2})}{2 B(m+1,m+1)}
\int_{0}^{1} \frac{x_{1}^{m} (1-x_{1})^{m} \, dx_{1}}
{ \sqrt{x_{1}e^{i \theta} + (1-x_{1})e^{-i \theta}}}.
\label{formula-10}
\end{equation}

The next lemma provides an evaluation of the integral in (\ref{formula-10}).

\begin{Lem}
Define
\begin{equation}
I_{m}(z) := \int_{0}^{1} \frac{x^{m}(1-x)^{m} \, dx}{\sqrt{1 + (z^{2}-1)x}}.
\end{equation}
\noindent
Then
\begin{equation}
I_{m}(z) = \frac{2^{2m+2}}{(m+1) \binom{2m+1}{m} \binom{4m+2}{2m+1}
(z+1)^{2m+1}} \sum_{k=0}^{m} \binom{2m-2k}{m-k} \binom{m+k}{k} (z+1)^{2k}
z^{m-k}.
\nonumber
\end{equation}
\end{Lem}
\begin{proof}
Expand the integrand in powers of $x$ and integrate. The result is simplified
using the value for the beta function
\begin{equation}
B(a,b) = \frac{a+b}{ab} \binom{a+b}{a}^{-1}
\end{equation}
\noindent
for $a, \, b \in \mathbb{N}$. It follows that
\begin{eqnarray}
I_{m}(z) & = & \sum_{k=0}^{\infty} \binom{-\tfrac{1}{2}}{k} (z^{2}-1)^{k}
\int_{0}^{1} x^{m+k} (1-x)^{m} \, dx \nonumber \\
& = & \sum_{k=0}^{\infty} \binom{-\tfrac{1}{2}}{k} B(m+k+1,m+1) (z^{2}-1)^{k}
\nonumber \\
& = & \frac{1}{m+1} \sum_{k=0}^{\infty} \binom{2k}{k} \binom{2m+k+1}{m+1}^{-1}
\left( \frac{1-z^{2}}{4} \right)^{k}.
\nonumber
\end{eqnarray}
\noindent
The result now follows from the identity
\begin{multline}
 \sum_{k=0}^{\infty} \binom{2k}{k} \binom{2m+k+1}{m+1}^{-1}
\left( \frac{1-z^{2}}{4} \right)^{k} =  \label{formula-00} \\
\frac{2^{2m+2}}{\binom{2m+1}{m} \binom{4m+2}{2m+1}
(z+1)^{2m+1}} \sum_{k=0}^{m} \binom{2m-2k}{m-k} \binom{m+k}{k} (z+1)^{2k}
z^{m-k}.
\end{multline}
\indent
To establish (\ref{formula-00}), the Wilf-Zeilberger automated method
\cite{aequalsb} is employed. It provides the recurrence
\begin{multline}
(z^{2}-1)^{2} (4m+9)(4m+7)a_{m+2} - (z^{4}-6z^{2}+1)(2m+6)(2m+3)a_{m+1} -
\\
4z^{2}(m+2)(m+3)a_{m} = 0
\nonumber
\end{multline}
\noindent
that both sides of (\ref{formula-00}) satisfy. Checking the initial values
at $m=0$ and $m=1$ is elementary and the proof is complete.
\end{proof}

\smallskip

The evaluation of (\ref{identity-1}) is restated in a slightly different
form, as the content of the next statement.

\begin{Thm}
Let $m \in \mathbb{N}, \, z = e^{i \theta}$ and $a  = \cos \theta$. Then
\begin{equation}
N_{0,4}(a;m) = \frac{B(\tfrac{1}{2}, 2m+ \tfrac{3}{2})}{2 B(m+1,m+1)}
\sqrt{z}I_{m}(z)
\end{equation}
\noindent
holds.
\end{Thm}

\section{A connection with Schur functions} \label{S:schur}

The integral
\begin{equation}
G_{n}(\mathbf{q}) = \frac{2}{\pi} \int_{0}^{\infty} \prod_{k=1}^{n}
\frac{1}{x^{2}+q_{k}^{2}} \,dx,
\end{equation}
\noindent
is clearly a symmetric function of the parameters $q_{1}, \ldots, q_{n}$. It has
been expressed in \cite{aems-1} as
\begin{equation}
G_{n}(\mathbf{q}) = \frac{s_{\lambda(n-1)}(\mathbf{q})}
                {e_{n}(\mathbf{q}) \, s_{\lambda(n)}(\mathbf{q})},
\label{schur-1}
\end{equation}
\noindent
where $\lambda(n)$ is the partition $\lambda(n) = (n-1,n-2,\ldots,1)$,
$e_{n}(\mathbf{q}) = q_{1}q_{2}\cdots q_{n}$ and the Schur function
corresponding to a partition $\mu$ defined by
\begin{equation}
s_{\mu}(\mathbf{q}) = \frac{a_{\mu+\lambda(n)}(\mathbf{q})}
                           {a_{\lambda(n)}(\mathbf{q})}
\end{equation}
\noindent
with $a_{\mu}(\mathbf{q}) := \text{det}(q_{i}^{\mu_{j}})_{1\leq i,j
\leq n}$.

\smallskip

In the current problem, let $a = \tfrac{1}{2}(w^{2}+w^{-2})$ so that
\begin{equation}
N_{0,4}(a;m) =  \int_{0}^{\infty} \frac{dx}{(x^{2}+w^{2})^{m+1} \,
(x^{2}+w^{-2})^{m+1}}.
\end{equation}
\noindent
Now take the partition $\lambda$ as above
and the $(2m+2)-$tuple $\mathbf{q}$ given by
\begin{equation}
\lambda(n) = (n-1, \, n-2, \ldots, 1), \quad \mathbf{q}
 = (w, \ldots, w,w^{-1},\ldots,
w^{-1})
\end{equation}
\noindent
where $\mathbf{q}$ has $m+1$ copies of $w$ and
$w^{-1}$. The expression (\ref{schur-1}), that appears as
Theorem 5.1 of \cite{aems-1},
provides the next result.

\begin{Thm}
\label{thm-schur}
The quartic integral is given by
\begin{equation}
N_{0,4}(a;m) = \frac{\pi}{2} \frac{s_{\lambda(2m+1)}(\mathbf{q})}
{s_{\lambda(2m+2)}(\mathbf{q})}.
\end{equation}
\end{Thm}

As a corollary, this implies a specialization for the Schur functions.

\begin{Cor}
If $q = (w, \ldots,w,w^{-1},\ldots,w^{-1})$ is of length $2m+2$, then
\begin{equation}
\frac{
s_{\lambda(2m+1)}(w,\ldots,w,w^{-1},\ldots,w^{-1}) }
{s_{\lambda(2m+3)}(w,\ldots,w,w^{-1},\ldots,w^{-1}) }
=
2^{-2m} \sum_{k=0}^{m} \binom{2k}{k} \binom{2m-k}{m}
\frac{1}{(w+w^{-1})^{2k+1}}.
\nonumber
\end{equation}
\end{Cor}

Theorem \ref{thm-schur} has a natural generalization to integrals of the form
\begin{equation}
\int_{0}^{\infty} \frac{dx}{Q_{2n}(x)^{m}},
\nonumber
\end{equation}
\noindent
where
\begin{equation}
Q_{2n}(x) = x^{2n} + a_{1}x^{2n-2} + a_{2}x^{2n-4} + \cdots + a_{2}x^{4} +
a_{1}x^{2} + 1
\end{equation}
\noindent
is a {\em palindromic polynomial}. These polynomials factor as
\begin{equation}
Q_{2n}(x) = (x^{2}+w_{1}^{2})(x_{2}+w_{1}^{-2}) \cdots
(x^{2}+w_{n}^{2})(x_{2}+w_{n}^{-2}).
\end{equation}

Consider the partition $\lambda$ and the $n$-concatenation
${\mathbf{q}} = (q_{1}, \ldots, q_{n})$, of the $2m$-tuples $q_{j}$, by
\begin{equation}
\lambda(s) = (s-1,s-2,\ldots,1), \text{   and   }
q_{j} = (w_{j},\ldots,w_{j};w_{j}^{-1},\ldots,w_{j}^{-1})
\end{equation}
\noindent
where $q_{j}$ has $m$ copies of each $w_{j}$ and $w_{j}^{-1}$. A direct
application of (\ref{schur-1})  gives the next result.

\begin{Prop}
Preserving the notation from above,
\begin{equation}
\int_{0}^{\infty} \frac{dx}{Q_{2n}(x)^{m}} =
\frac{\pi}{2} \frac{s_{\lambda(2nm-1)}(q)}{s_{\lambda(2nm)(q)}}.
\end{equation}
\end{Prop}

\section{A connection with Bessel functions} \label{S:Bessel}

The quartic integral (\ref{identity-1})
is evaluated in this section using the modified
Bessel function $K_{\alpha}(x)$. The classical Bessel function
$J_{\alpha}(x)$ is usually defined via the differential equation
\begin{equation}
\frac{d^{2}y}{dx^{2}} + \frac{1}{x} \frac{dy}{dx} -
\left( 1 + \frac{\alpha^{2}}{x^{2}} \right) y = 0.
\label{diff-bes}
\end{equation}
\noindent
The modified function $I_{\alpha}(x) = e^{-\alpha \pi i/2} J_{\alpha}
(xe^{\pi i/2})$ is also real valued. In the case $\alpha \not \in
\mathbb{Z}$, the functions $I_{\alpha}(x)$ and $I_{-\alpha}(x)$ form a basis
for the solutions of (\ref{diff-bes}). To deal with the situation of
integer parameters, it is convenient to define
\begin{equation}
K_{\alpha}(x) = \frac{\pi}{2 \, \sin \pi \alpha}
\left[ I_{-\alpha}(x) - I_{\alpha}(x) \right].
\end{equation}
\noindent
This is the modified Bessel function of second kind, also called
MacDonald function.

The proof of (\ref{identity-1}) proceeds along the following lines: the first
step is to obtain an expression for $N_{0,4}(a,m)$ as an integral involving
$K_{1/4}(x)$. An entry in \cite{prudnikov1}  gives a hypergeometric form of
this integral. The final step is to prove a hypergeometric identity that
transforms this form to a result in \cite{bomohyper}. Details about Bessel
functions can be found in \cite{andrews3} and chapter 10 of \cite{nist}.

\smallskip

The identity
\begin{equation}
\int_{0}^{\infty} e^{-c u} u^{m} \, du  =
\frac{\Gamma(m+1)}{c^{m+1}},
\end{equation}
\noindent
valid for $\mathop{\rm Re}\nolimits{c}>0$, is used with $c = x^{4}+2ax^{2}+1$ to produce
\begin{equation}
\frac{1}{(x^{4} + 2ax^{2} + 1)^{m+1}} =
\frac{1}{\Gamma(m+1)} \int_{0}^{\infty} e^{-u} u^{m}
e^{-u(x^{4}+2ax^{2})} \, du.
\end{equation}
\noindent
Integration yields
\begin{equation}
N_{0,4}(a;m) = \frac{1}{\Gamma(m+1)} \int_{0}^{\infty} e^{-u}
u^{m} \int_{0}^{\infty} e^{-u(x^{4}+2ax^{2})} \, du \, dx.
\end{equation}
\noindent
The restriction $\mathop{\rm Re}\nolimits{c}>0$ is satisfied by taking $-1 <  a < 1$.

Entry $3.469.1$ in \cite{gr} gives
\begin{equation}
\int_{0}^{\infty} e^{-\mu x^{4} - 2 \nu x^{2}} \, dx =
\frac{1}{4} \sqrt{\frac{2 \nu}{\mu}} \text{exp} \left( \frac{\nu^{2}}
{2 \mu} \right) K_{\small{1/4}} \left( \frac{\nu^{2}}{2 \mu} \right).
\end{equation}

The choice
$\mu = u$ and $\nu = a u$ yields:

\begin{Lem}
\label{lemma-bes1}
The quartic integral $N_{0,4}(a;m)$ is given by
\begin{equation}
N_{0,4}(a;m) = \frac{2^{m - \tfrac{1}{2}}}{\Gamma(m+1) \, a^{2m+ \tfrac{3}{2}}}
\int_{0}^{\infty} t^{m} e^{-bt} K_{\small{1/4}}(t) \, dt,
\end{equation}
\noindent
where $b = 2/a^{2}-1$.
\end{Lem}

This expression can be evaluated using \cite[vol.2, 2.16.6.2]{prudnikov1}
\begin{equation}
\int_{0}^{\infty}x^{\alpha-1}e^{-px}K_{\nu}\left(cx\right)dx
=\frac{\left(2c\right)^{\nu}\sqrt{\pi}}{\left(p+c\right)^{\alpha+\nu}}
\frac{\Gamma(\alpha-\nu) \, \Gamma(\alpha+\nu)}{\Gamma(\alpha + \tfrac{1}{2})}
\,_{2}F_{1}\left(\begin{array}{c}
\alpha+\nu,\nu+\frac{1}{2}\\
\alpha+\frac{1}{2}\end{array}\Big{\vert}\frac{p-c}{p+c}\right)
\nonumber
\end{equation}
valid for $\mathop{\rm Re}\nolimits{(c+p)}>0,\,\,\mathop{\rm Re}\nolimits{\alpha} >
\vert \mathop{\rm Re}\nolimits{\nu} \vert.$  This yields the expression
\begin{equation}
N_{0,4}(a;m) =
\frac{\sqrt{\pi} a}{2 \sqrt{2} m!}
\frac{\Gamma\left(m+\tfrac{3}{4}\right) \Gamma(m + \tfrac{5}{4})}
{\Gamma\left(m+ \tfrac{3}{2} \right)}
\,_{2}F_{1}\left(\begin{array}{c}
m+\frac{5}{4},\frac{3}{4}\\
m+\frac{3}{2}\end{array} \Big{\vert} \, 1-a^{2}\right).
\end{equation}

Using \cite[15.3.24]{abramowitz1} this can be written as

\begin{equation}
N_{0,4}\left(a;m\right)
=\frac{\sqrt{\pi}}{2}\frac{\Gamma\left(2m+\frac{3}{2}\right)}
{\Gamma\left(2m+2\right)}\,_{2}F_{1}\left(\begin{array}{c}
2m+\frac{3}{2},\frac{1}{2}\\
m+\frac{3}{2}\end{array}\Big{\vert}\frac{1-a}{2}\right).
\label{original}
\end{equation}

\smallskip

The proof of (\ref{identity-1}) is now based on the representation
\begin{equation}
N_{0,4}\left(a;m\right)
= \frac{2^{m-\tfrac{1}{2}}}{(a+1)^{m+1/2}} B\left( 2m + \tfrac{3}{2},
\tfrac{1}{2} \right)
\,_{2}F_{1}\left(\begin{array}{c}
-m,m+1\\
m+\frac{3}{2}\end{array}\Big{\vert}\frac{1-a}{2}\right)
\end{equation}
\noindent
described in \cite{bomohyper}. In particular, the polynomial $P_{m}(a)$
is identified there as the Jacobi polynomial $P_{m}^{(m+1/2,-m-1/2)}(a)$.

Matching  both representations for $N_{0,4}(a;m)$ shows that (\ref{identity-1})
follows from the next result.

\begin{Prop}
\label{hyper-iden1}
The identity
\begin{equation}
_{2}F_{1}\left(\begin{array}{c}
-m,m+1\\
m+\frac{3}{2}\end{array}\Big{\vert} z \right) =
(1-2z) (1-z)^{m+\tfrac{1}{2}} \, \,
_{2}F_{1}\left(\begin{array}{c}
m + \tfrac{5}{4}, \tfrac{3}{4}\\
m+\frac{3}{2}\end{array}\Big{\vert} \, 4z(1-z) \right)
\label{hyp-0}
\end{equation}
\noindent
holds.
\end{Prop}
\begin{proof}
The WZ-method shows that both sides satisfy the recurrence
\begin{multline}
(4m+7)(4m+9)za_{m+2} + (2m+3)(2m+5)(4z^{2}-4z-1)a_{m+1} \label{recur-11} \\
-(2m+3)(2m+5)(z-1)a_{m} = 0.
\end{multline}

To complete the proof, it suffices to check two initial values. The required
identities for $m=0$ and $m=1$ are written in terms of
$t = 4z(1-z)$ in the form
\begin{equation}
\label{init-1}
_{2}F_{1}\left(\begin{array}{c}
\tfrac{5}{4}, \tfrac{3}{4}\\
\frac{3}{2}\end{array}\Big{\vert} \, t \right)
= \frac{\sqrt{2}}{(1-t)^{1/2} \, \sqrt{1 + \sqrt{1-t}}}
\end{equation}
\noindent
and
\begin{equation}
\label{init-2}
_{2}F_{1}\left(\begin{array}{c}
\tfrac{9}{4}, \tfrac{3}{4}\\
\frac{5}{2}\end{array}\Big{\vert} \, t \right)
= \frac{2 \sqrt{2} (3 + 2 \sqrt{1-t})}{5(1+\sqrt{1-t})^{3/2} \, \sqrt{1-t}}.
\end{equation}

The proofs of these identities are elementary using the fact that the
hypergeometric differential equation
\begin{equation}
t(1-t) y'' - \left[ c - (a+b+1)t \right]y' - aby = 0.
\end{equation}
\noindent
has a unique solution analytic at $t=0$ with initial value $y(0)=1$. Indeed,
both sides of (\ref{init-1}) satisfy
\begin{equation}
t(1-t)y'' + \left( \tfrac{3}{2} - 3t \right)y' - \tfrac{15}{16}y = 0
\end{equation}
\noindent
with value $y(0)=1$. This establishes (\ref{init-1}). The same holds for
(\ref{init-2}) using
\begin{equation}
t(1-t)y'' + \left( \tfrac{5}{2} - 4t \right)y' - \tfrac{27}{16}y = 0
\end{equation}
\end{proof}

\begin{Note}
The identity in Proposition \ref{hyper-iden1} can also be established using
C. Koutschan package \texttt{HolonomicFunctions} \cite{Koutschan10b}. Denote
the left hand side of (\ref{hyp-0}) by $F1$ and the right hand side by $F2$.
The command
\begin{center}
\texttt{Annihilator}$[F1,S[m]]$
\end{center}
\noindent
 finds the operator
\begin{equation}
(7+4m)(9+4m)zS_{m}^{2} +
(3 +2m)(5+2m)(-1-4z+4z^{2})S_{m} -
(3+2m)(5+2m)(-1+z)
\nonumber
\end{equation}
\noindent
with $S_{m}$ being the shift in the discrete parameter $m$; that is,
$S_{m}g(m;x) = g(m+1;x)$. The
package claims that this
operator annihilates the left hand side of (\ref{hyp-0}). This is
precisely the recurrence (\ref{recur-11}) obtained before. The command
\begin{center}
\texttt{Annihilator}$[F1,S[m]]$ ==
\texttt{Annihilator}$[F2,S[m]]$
\end{center}
\noindent
returns \texttt{True}, showing that the right
hand side $F2$ satisfies the same recurrence. The initial conditions can also
be determined automatically.
\end{Note}

\section{A second proof of the main identity via  Bessel functions} \label{S:Bessel1}

The previous section contains the identity (\ref{lemma-bes1}) giving the
quartic integral $N_{0,4}(a,m)$ as an integral involving Bessel functions:
\begin{equation}
N_{0,4}(a;m) = \frac{2^{m - \tfrac{1}{2}}}{\Gamma(m+1) \, a^{2m+ \tfrac{3}{2}}}
\int_{0}^{\infty} t^{m} e^{-bt} K_{\small{1/4}}(t) \, dt.
\label{form-00}
\end{equation}
\noindent
The table of integrals \cite{gr} contains, as Entry $6.611.3$, the
special case $m=0$:
\begin{equation}
\int_{0}^{\infty} e^{-bt} K_{\small{1/4}}(t) \, dt =
\frac{\pi}{\sqrt{2} \, \sqrt{b^{2}-1}}
\left[ (b + \sqrt{b^{2}-1})^{1/4} - (b+\sqrt{b^{2}-1})^{-1/4} \right]
\end{equation}
\noindent
that can be written as
\begin{equation}
\int_{0}^{\infty} e^{-bt} K_{\small{1/4}}(t) \, dt =
\frac{\pi}{2} \frac{a^{3/2}}{\sqrt{1+a}},
\label{bes-11}
\end{equation}
\noindent
using $b = a^{2}/2-1$. Differentiating
$m$ times with respect to $b$, it follows that
\begin{equation}
N_{0,4}(a,m) = \frac{(-1)^{m} 2^{m-1/2}}{m! \, a^{2m+3/2}}
\frac{\partial^{m}}{\partial b^{m}} \left( \frac{\pi}{2}
\frac{a^{3/2}}{\sqrt{1+a}}\right).
\label{bes-12}
\end{equation}

The proof of (\ref{identity-1}) is thus reduced to finding an analytic
expression for the derivatives in (\ref{bes-12}).

\begin{Prop}
There exists a polynomial $Q_{m}(a)$ such that
\begin{equation}
\frac{\partial^{m}}{\partial b^{m}}
\left(\frac{a^{3/2}}{\sqrt{1+a}} \right)
= \frac{(-1)^{m} m!}{2^{2m}} \frac{a^{2m+3/2}}
{(1+a)^{m+1/2}} Q_{m}(a).
\end{equation}
\end{Prop}
\begin{proof}
The case $m=0$ holds with $Q_{0}(m) = 1$.  To complete
the inductive step, it is shown that if $Q_{m}(a)$ is a polynomial then the
function $Q_{m+1}(a)$ defined by the relation
\begin{equation}
\frac{\partial}{\partial b}
\left(
\frac{(-1)^{m} m!}{2^{2m}} \frac{a^{2m+3/2}}
{(1+a)^{m+1/2}} Q_{m}(a) \right) =
\frac{(-1)^{m+1} (m+1)!}{2^{2m+2}} \frac{a^{2m+7/2}}
{(1+a)^{m+3/2}} Q_{m+1}(a),
\label{def-tn}
\end{equation}
\noindent
is also a polynomial.

The chain rule and
$\frac{\partial a}{\partial b}
= - \tfrac{1}{4}a^{3}$ show that (\ref{def-tn}) is equivalent to
\begin{multline}
2(m+1)Q_{m+1}(a) =
\left[2(m+1)(a+1) + (2m+1) \right]Q_{m}(a) + \label{rec-Tm} \\
\left[ 2(a+1)^{2} - 2(a+1) \right]Q_{m}'(a).
\end{multline}
\noindent
Thus $Q_{m+1}(a)$ is a polynomial in $a$.
\end{proof}

The proof of (\ref{identity-1}) is now reduced to checking that the polynomial
$Q_{m}$ in the previous lemma is given by $P_{m}(a)$ defined
in (\ref{polyP-def}).

\begin{Thm}
The polynomial $Q_{m}(a)$ is the same as $P_{m}(a)$.
\end{Thm}
\begin{proof}
It suffices to check that
$P_{m}(a)$ satisfies the same recurrence as $Q_{m}(a)$. To this end, compare
the coefficients of $(1+a)^{k}$ on both sides of (\ref{rec-Tm}). The result
is equivalent to
\begin{eqnarray}
2^{k} (m+1) \binom{2m-2k+2}{m-k+1} \binom{m+k+1}{m+1} & = &
2^{k+1} (m+1)\binom{2m-2k+2}{m-k+1} \binom{m+k-1}{m}  \nonumber \\
& + &  2^{k+1} (2m+1)\binom{2m-2k}{m-k} \binom{m+k}{m}  \nonumber \\
& + &  2^{k+1} (k-1)\binom{2m-2k+2}{m-k+1} \binom{m+k-1}{m}  \nonumber \\
& - &  2^{k+2} k\binom{2m-2k}{m-k} \binom{m+k}{m}. \nonumber
\end{eqnarray}

\noindent
This however can routinely be verified. Simply
divide through by $\binom{2m-2k}{m-k} \binom{m+k}{m}$
and the statement reduces to a simple polynomial identity.
\end{proof}

\medskip

A small variation of this  proof of (\ref{identity-1}) is obtained by
differentiating
\begin{equation}
\int_{0}^{\infty} e^{-bt} K_{\tfrac{1}{4}}(t) \, dt =
\frac{\pi}{2^{1/4}(b+1)^{1/2} \, \sqrt{\sqrt{2} + \sqrt{b+1}}}.
\end{equation}
directly with respect to the parameter $b$. The first few
examples suggest the next result.

\begin{Lem}
There are polynomials $S_{m}, \, T_{m}$ such that
\begin{multline}
\frac{d^{m}}{db^{m}}
\left( \frac{1}{(b+1)^{1/2} \, \sqrt{\sqrt{2} + \sqrt{b+1}}} \right) =
\label{formula-77} \\
(-1)^{m} \frac{S_{m}(b) + \sqrt{b+1} T_{m}(b)}
{2^{2m} (b+1)^{m+1/2} \, (\sqrt{2} + \sqrt{b+1})^{m+1/2}}.
\end{multline}
\end{Lem}
\begin{proof}
The proof is by induction on $m$, and is obtained upon differentiating the
stated expression for the $m$-th derivative. The base case $m=0$ is obvious.
Assume (\ref{formula-77}) holds for $m$. Then differentiating the right
hand side of (\ref{formula-77}) generates the recurrence
\begin{eqnarray}
S_{m+1}(b) & = & 2 \sqrt{2}(2m+1) S_{m}(b) \label{recu-ST} \\
 & & -(b+1) \left[ 4 \sqrt{2} S_{m}'(b) -
(1+6m) T_{m}(b) + 4(b+1)T_{m}'(b) \right] \nonumber \\
T_{m+1}(b) & = & 3(2m+1)S_{m}(b) + 4 \sqrt{2}m T_{m}(b) \nonumber \\
 & & - 4(b+1) \left( S_{m}'(b) + \sqrt{2} T_{m}'(b) \right), \nonumber
\end{eqnarray}
\noindent
where $b = 2/a^{2}-1$. The initial condition $S_{0}(b)=1$ and $T_{0}(b)=0$
yield the result.
\end{proof}

\begin{Lem}
Let $m' = \lfloor{ m/2 \rfloor}$. Define
\begin{equation}
U_{m}(b) = \frac{2^{-2m}}{(1+b)^{m'}}
\sum_{k=0}^{m} 2^{k} \binom{2m-2k}{m-k} \binom{m+k}{m}
\sum_{j=0}^{\lfloor{k/2 \rfloor}} \binom{k}{2j} 2^{j} (1+b)^{m'-j}
\nonumber
\end{equation}
\noindent
and
\begin{equation}
V_{m}(b) = \frac{2^{-2m}}{(1+b)^{m'}}
\sum_{k=0}^{m} 2^{k} \binom{2m-2k}{m-k} \binom{m+k}{m}
\sum_{j=0}^{\lfloor{(k-1)/2 \rfloor}} \binom{k}{2j+1} 2^{j} (1+b)^{m'-j}.
\nonumber
\end{equation}
\noindent
Then
\begin{equation}
P_{m}(a) = U_{m}(b) + a V_{m}(b).
\end{equation}
\end{Lem}
\begin{proof}
The polynomial $P_{m}(a)$
is decomposed into its even and odd part using
\begin{equation}
(1+a)^{k} = \frac{1}{2} \left[ (1+a)^{k} + (1-a)^{k} \right] +
            \frac{1}{2} \left[ (1+a)^{k} - (1-a)^{k} \right].
\end{equation}
\noindent
The result follows.
\end{proof}

\smallskip

The proof of (\ref{identity-1}) now reduces to expressing the polynomials
$S_{m}(b)$ and $T_{m}(b)$ in terms of $U_{m}(b)$ and $V_{m}(b)$. This is
given by
\begin{equation}
T_{m}(b) = \frac{m! 2^{3m/2}}{\sqrt{2} a^{m-1}}
\times \begin{cases}
U_{m}(b) & \quad \text{ if } m \text{ is odd} \\
aV_{m}(b) & \quad \text{ if } m \text{ is even}
\end{cases}
\end{equation}
\noindent
and
\begin{equation}
S_{m}(b) = \frac{m! 2^{3m/2}}{\sqrt{2} a^{m}}
\times \begin{cases}
aV_{m}(b) & \quad \text{ if } m \text{ is odd} \\
U_{m}(b) & \quad \text{ if } m \text{ is even}.
\end{cases}
\end{equation}
\noindent
The details are elementary and are left to the reader.

\begin{Note}
The identity (\ref{bes-12}) yields
\begin{equation}
N_{0,4}(a,m) = \frac{\pi \, a^{m-1}}{m! 2^{5m/2+3/2} \, (a+1)^{m+1/2}}
\left[ \sqrt{2} T_{m}(b) + S_{m}(b) \right].
\label{form-99}
\end{equation}
\end{Note}

\smallskip
\medskip

\section{Conclusions} \label{S:conc}

The value of the definite integral
\begin{equation}
N_{0,4}(a;m) = \int_{0}^{\infty} \frac{dx}{(x^{4}+2ax^{2}+1)^{m+1}},
\nonumber
\end{equation}
\noindent
is given by
\begin{equation}
N_{0,4}(a;m)  =   \frac{\pi}{2^{m+3/2} (a+1)^{m+1/2} }
P_{m}(a)
\nonumber
\end{equation}
\noindent
where
\begin{equation}
P_{m}(a) = 2^{-2m} \sum_{k=0}^{m} 2^{k} \binom{2m-2k}{m-k} \binom{m+k}{k}
(a+1)^{k}.
\nonumber
\end{equation}

\medskip

Several proofs of this identity appeared in the literature. The current
work includes proofs relating this identity to Bessel functions, Schur
polynomials and a method of Schwinger for the evaluation of definite
integrals.

\medskip

\noindent
{\bf Acknowledgements}. The authors wish to thank Armin Straub for
comments on C. Koutschan package. The work of the second author was
partially supported by NSF-DMS 0070567. C. Vignat thanks V. Moll for his
invitation to visit Tulane in July 2010. \\