# Model independent constraints from vacuum and in-medium QCD Sum Rules

<sup>1</sup><sup>1</sup>institutetext: Physik-Department, Theoretische Physik, Technische Universität München, D-85747 Garching, Germany (Received: date / Revised version: date)  

## Abstract

We discuss QCD sum rule constraints based on moments of vector meson spectral distributions in the vacuum and in a nuclear medium. Sum rules for the two lowest moments of these spectral distributions do not suffer from uncertainties related to QCD condensates of dimension higher than four. We exemplify these relations for the case of the $`\omega `$ meson and discuss the issue of in-medium mass shifts from this viewpoint.

offprints: weise@physik.tu-muenchen.de  

QCD sum rules have repeatedly been used in recent times to arrive at estimates for possible in-medium mass shifts of vector mesons 2 ; 7 . The validity of such estimates has been questioned, however, for several reasons. First, for broad structures such as the $`\rho `$ meson whose large vacuum decay width is further magnified by in-medium reactions, the QCD sum rule analysis does not provide a reliable framework to extract anything like a “mass shift” 4 ; 5 . Secondly, notorious uncertainties exist at the level of factorization assumptions commonly used to approximate four-quark condensates in terms of $`\overline{q}q^2`$, the square of the standard chiral condensate. The first objection is far less serious for the $`\omega `$ meson which may well have a chance to survive as a reasonably narrow quasi-particle in nuclear matter 4 ; 4b . The second objection, however, is difficult to overcome: factorization of four-quark condensates may indeed be questionable.

In the present note we focus on the two lowest moments ($`𝑑ss^nR(s)`$ with $`n=0,1`$) of vector meson spectral distributions, in vacuum as well as in nuclear matter, and point out that these are subject to sum rules which do not suffer from the uncertainties introduced by four-quark condensates. These sum rules are shown to provide useful, model independent constraints which we exemplify for the case of the $`\omega `$ meson spectral distribution and its change in the nuclear medium. The sum rule for the second moment, $`𝑑ss^2R(s)`$, does involve the four-quark condensate. In fact it can be used in principle to determine this particular condensate and test the factorization assumption. The detailed analysis of this question will be defered to a longer paper. In this short note we confine ourselves to conclusions that can be drawn without reference to four-quark condensates.

The starting point is the current-current correlation function

$$\mathrm{\Pi }_{\mu \nu }(q)=id^4xe^{iqx}𝒯j_\mu (x)j_\nu (0)$$
(1)   
where $`𝒯`$ denotes the time-ordered product and the expectation value is taken either in the vacuum or in the ground state of nuclear matter at rest. In vacuum the polarization tensor (1) can be reduced to a single scalar correlation function, $`\mathrm{\Pi }(q^2)={\scriptscriptstyle \frac{1}{3}}g^{\mu \nu }\mathrm{\Pi }_{\mu \nu }(q)`$. In nuclear matter there are two (longitudinal and transverse) correlation functions which coincide for a meson at rest with respect to the medium (i.e. with $`q^\mu =(\omega ,\stackrel{}{q}=0)`$).

The reduced correlation function is written as a (twice subtracted) dispersion relation,

$$\mathrm{\Pi }(q^2)=\mathrm{\Pi }(0)+\mathrm{\Pi }^{}(0)q^2+\frac{q^4}{\pi }𝑑s\frac{\mathrm{Im}\mathrm{\Pi }(s)}{s^2(sq^2iϵ)}.$$
(2)   
where $`\mathrm{\Pi }(0)`$ vanishes in vacuum but contributes in nuclear matter. At large spacelike $`Q^2=q^2>0`$ the QCD operator product (Wilson) expansion gives

$$12\pi ^2\mathrm{\Pi }(q^2=Q^2)=c_0Q^2\mathrm{ln}\left(\frac{Q^2}{\mu ^2}\right)+c_1+\frac{c_2}{Q^2}+\frac{c_3}{Q^4}+\mathrm{}$$
(3)   
We specify the coefficients $`c_i`$ for the isoscalar current $`j^\mu ={\scriptscriptstyle \frac{1}{6}}(\overline{u}\gamma ^\mu u+\overline{d}\gamma ^\mu d)`$, the case of the $`\omega `$ meson that we wish to use here for explicit evaluations. In vacuum we have:

$`c_0`$ $`=`$ $`{\displaystyle \frac{1}{6}}\left(1+{\displaystyle \frac{\alpha _S}{\pi }}\right),c_1={\displaystyle \frac{1}{2}}(m_u^2+m_d^2),`$ (4)   
$`c_2`$ $`=`$ $`{\displaystyle \frac{\pi ^2}{18}}{\displaystyle \frac{\alpha _S}{\pi }}𝒢^{\mu \nu }𝒢_{\mu \nu }+{\displaystyle \frac{2\pi ^2}{3}}m_uu\overline{u}+m_dd\overline{d},`$ (5)   
while $`c_3`$ involves combinations of four-quark condensates of (mass) dimension 6. The quark mass term $`c_1`$ is small and can be dropped in the actual calculations. For the gluon condensate we use $`\frac{\alpha _S}{\pi }𝒢^2=(0.36\mathrm{GeV})^4`$ 9 , and the (chiral) quark condensate is given by $`m_u\overline{u}u+m_d\overline{d}dm_q\overline{u}u+\overline{d}d=m_\pi ^2f_\pi ^2(0.11\mathrm{GeV})^4`$ through the Gell-Mann, Oakes, Renner relation.

In the nuclear medium with baryon density $`\rho `$ we have $`c_i(\rho )=c_i(\rho =0)+\delta c_i(\rho )`$ with $`c_i(0)`$ given by eqs.(4,5), and

$$\delta c_2(\rho )=\frac{\pi ^2}{3}\left[\frac{4}{27}M_N^{(0)}+2\sigma _N+A_1M_N\right]\rho $$
(6)   
to linear order in $`\rho `$. The first term in brackets is the leading density dependent correction to the gluon condensate and involves the nucleon mass in the chiral limit, $`M_N^{(0)}0.75\mathrm{GeV}`$ 6 . The second part proportional to the nucleon sigma term $`\sigma _N45\mathrm{MeV}`$ is the first order correction of the quark condensate, and the third term introduces the first moment of the quark distribution function in the nucleon:

$$A_1=2𝑑xx\left[u(x)+\overline{u}(x)+d(x)+\overline{d}(x)\right].$$
(7)   
It represents twice the fraction of momentum carried by quarks in the proton. We take $`A_11`$ as determined by deep-inelastic lepton scattering at $`Q2\mathrm{GeV}`$. Note that $`\delta c_2(\rho _0)410^3\mathrm{GeV}^4`$ at $`\rho =\rho _0=0.17\mathrm{fm}^3`$, the density of nuclear matter, and almost all of this correction comes from the term proportional to $`A_1`$.

Next we introduce the Borel transform of eq. (3):

$$12\pi ^2\mathrm{\Pi }(0)+_0^{\mathrm{}}𝑑sR(s)e^{s/^2}=c_0^2+c_1+\frac{c_2}{^2}+\frac{c_3}{2^4}+\mathrm{}$$
(8)   
with $`R(s)=\frac{12\pi }{s}\mathrm{Im}\mathrm{\Pi }(s)`$ and $`\mathrm{\Pi }(0)=\rho /4M_N`$, the vector meson analogue of the Thomson term in photon scattering.

We separate the spectrum $`R(s)`$ into a resonance part with $`ss_0`$ and a continuum $`R_c(s)`$ which must approach the perturbative QCD limit for $`s>s_0`$:

$$R_c(s)=\frac{1}{6}\left(1+\frac{\alpha _S}{\pi }\right)\mathrm{\Theta }(ss_o).$$
(9)   
The factor $`{\scriptscriptstyle \frac{1}{6}}`$ is again specific for the isoscalar channel. The Borel mass parameter $``$ must be sufficiently large so that eq.(8) converges rapidly, but otherwise it is arbitrary. We choose $`>\sqrt{s_0}`$ so that $`e^{s/^2}`$ can be expanded in powers of $`s/^2`$ for $`s<s_0`$. The remaining integral $`_{s_0}^{\mathrm{}}𝑑sR_c(s)e^{s/^2}`$ is evaluated inserting the running coupling strength $`\alpha _S(s_o)`$ in eq.(9). Then the term-by-term comparison in eq.(8) gives the following set of sum rules for the moments of the spectrum $`R(s)`$ (see also refs. 7 ; 8 )

$`{\displaystyle _0^{s_0}}𝑑sR(s)`$ $`=`$ $`s_0c_0+c_112\pi ^2\mathrm{\Pi }(0),`$ (10)   
$`{\displaystyle _0^{s_0}}𝑑ssR(s)`$ $`=`$ $`{\displaystyle \frac{s_0^2}{2}}c_0c_2,`$ (11)   
$`{\displaystyle _0^{s_0}}𝑑ss^2R(s)`$ $`=`$ $`{\displaystyle \frac{s_0^3}{3}}c_0+c_3.`$ (12)   
Note that the first two sum rules are well determined and represent useful constraints for the spectrum $`R(s)`$. Only the third sum rule (12) involves four-quark condensates which are uncertain. In this short paper we concentrate on eqs. (10,11). A detailed analysis of eq.(12) will be presented in a forthcoming longer paper. It is instructive to illustrate the sum rules (10,11) for the $`\omega `$ meson in vacuum using Vector Meson Dominance (VMD) for the resonant part of $`R(s)`$. In this model we have

$$R(s)=12\pi ^2\frac{m_\omega ^2}{g_\omega ^2}\delta (sm_\omega ^2)+\frac{1}{6}\left(1+\frac{\alpha _S}{\pi }\right)\mathrm{\Theta }(ss_o).$$
(13)   
with $`g_\omega =3g16.8`$ (using the vector coupling constant $`g=5.6`$). We can neglect the small quark mass term $`c_1`$ and find from eq. (10):

$$\frac{8\pi ^2}{g^2}\frac{m_\omega ^2}{s_0}=1+\frac{\alpha _S}{\pi },$$
(14)   
which fixes $`\sqrt{s_0}=1.16\mathrm{GeV}`$ using $`\alpha _S(s_0)0.4`$ and $`m_\omega =0.78\mathrm{GeV}`$. It is interesting to identify the spectral gap $`\mathrm{\Delta }=\sqrt{s_0}`$ with the scale for spontaneous chiral symmetry breaking, $`\mathrm{\Delta }=4\pi f_\pi `$, where $`f_\pi =92.4MeV`$ is the pion decay constant. In the VMD model, taking the zero width limit, eq.(14) holds for both the $`\omega `$ and $`\rho `$ meson, with equal mass $`m_V=m_\rho =m_\omega `$. Inserting $`s_0=16\pi ^2f_\pi ^2`$ in eq.(14) one recovers the famous KSFR relation $`m_V=\sqrt{2}gf_\pi `$ up to a small QCD correction.

The sum rule (11) for the first moment gives

$$\frac{8\pi ^2}{g^2}m_\omega ^4=\frac{s_0^2}{2}\left(1+\frac{\alpha _S}{\pi }\right)\frac{\pi ^2}{3}\left[\frac{\alpha _S}{\pi }𝒢^2+12m_u\overline{u}u+m_d\overline{d}d\right].$$
(15)   
Inserting the values for the gluon and quark condensates we find indeed perfect consistency. Given a model for the $`\omega `$ meson spectral function in the vacuum and in the nuclear medium, the sum rules (10,11) therefore provide useful constraints to test the calculated spectra.

We now continue on from VMD to a more realistic approach. In refs. 4 ; 4b  we have used an effective Lagrangian based on chiral $`SU(3)SU(3)`$ symmetry with inclusion of vector mesons as well as anomalous couplings from the Wess-Zumino action in order to calculate the $`\omega `$ meson spectrum both in the vacuum and in nuclear matter. The resulting vacuum spectrum reproduces the observed $`e^+e^{}`$hadrons$`(I=0)`$ data very well 4  (see Fig. 1a). The predicted in-medium mass spectrum (for $`\omega `$ excitations with $`\stackrel{}{q}=0`$) shows a pronounced downward shift of the $`\omega `$-meson peak and a substantial, but not overwhelming increase of its width from reactions such as $`\omega N\pi N,\pi \pi N`$ etc. (see Fig. 1b). At large $`s>s_0`$, both spectra should approach the QCD limit (9). The consistency test of these calculated spectral distributions with the sum rules (10) and (11) goes as follows:

* the vacuum case:   
the two sides of eq. (10),

$$_0^{s_0}𝑑sR(s)=\frac{s_0}{6}\left(1+\frac{\alpha _S(s_0)}{\pi }\right),$$
(16)   
now match at $`\sqrt{s_0}=1.25\mathrm{GeV}`$, with $`_0^{s_0}𝑑sR(s)=0.29\mathrm{GeV}^2`$. The sum rule for the first moment gives $`_0^{s_0}𝑑ssR(s)=0.19\mathrm{GeV}^4`$, to be compared with $`\frac{1}{2}s_0^2(1+\alpha _s/\pi )c_2=0.22\mathrm{GeV}^4`$, so there is consistency at the 10% level.
* the in-medium case:   
now we have to match the moments of the density dependent spectral distributions,

$$_0^{s_0}𝑑sR(s,\rho )=\frac{s_0}{6}\left(1+\frac{\alpha _S(s_0)}{\pi }\right)+\frac{3\pi ^2}{M_N}\rho ,$$
(17)   
together with

$$_0^{s_0}𝑑ssR(s,\rho )=\frac{s_0^2}{12}\left(1+\frac{\alpha _S(s_0)}{\pi }\right)c_2(0)\delta c_2(\rho ).$$
(18)   
Using our calculated spectrum 4b  shown in Fig. 1b, we find $`\sqrt{s_0}=1.08\mathrm{GeV}`$ at $`\rho =\rho _0=0.17\mathrm{fm}^3`$, with $`_0^{s_0}𝑑sR(s,\rho _0)=0.26\mathrm{GeV}^2`$. Then $`_0^{s_0}𝑑ssR(s,\rho _0)=0.11\mathrm{GeV}^4`$ is to be compared with the right hand side of eq. (18) which gives $`0.12\mathrm{GeV}^4`$, so there is again excellent consistency.

Note again that these tests do not involve uncertain four-quark condensates. Furthermore, if the in-medium spectrum shows a reasonably narrow quasi particle excitation, the quantity $`\overline{m}^2=_0^{s_0}𝑑ssR(s)/_0^{s_0}𝑑sR(s)`$ can indeed be interpreted as the square of an in-medium “mass” of this excitation. For our $`\omega `$ meson case we find $`\overline{m}=0.65\mathrm{GeV}`$ at $`\rho =\rho _0`$, a substantial downward mass shift as discussed in refs. 4 ; 4b . (For the broad $`\rho `$ meson spectrum, on the other hand , the interpretation of $`\overline{m}`$ as an in-medium mass is not meaningful as demonstrated in ref. 4 ).

Amusingly, the spectral gap $`\mathrm{\Delta }=\sqrt{s_0}`$ decreases by about 15 percent when replacing the vacuum by nuclear matter. This is in line with the proposition that this gap reflects the order parameter for spontaneous chiral symmetry breaking and scales like the pion decay constant $`f_\pi `$ (or, equivalently, like the square root of the chiral condensate $`\overline{q}q`$).

In summary, we have shown that the combination of sum rules (10) and (11) for the lowest moments of the spectral distributions does serve as a model-independent consistency test for calculated spectral functions.