The question of the conformal solution

May 4, 2009

There is currently a lot of activity to understand the behavior of the two-point functions in a Yang-Mills theory when quarks are not considered. The relevance of these results relies on the possibility to get working tools to manage low-energy phenomenology of QCD. E.g. we know quite well that the Nambu-Jona-Lasinio model is very successful to describe the behavior of hadronic matter and this model can be easily derived from full QCD if one knows the gluon propagator and the behavior of the ghost field. This idea dates back to a paper by Terry Goldman and Richard Haymaker (see here) on 1981. This can be summed up by saying that: given the gluon propagator one can get back a Nambu-Jona-Lasinio model. One could be able to derive it directly from QCD and this would be a great achievement explaining a lot of work done in several decades by a lot of smart researchers.

In these days it seems that we have to cope with a couple of scenarios and long lasting debate is still alive about them. One of this, the conformal solution, has been supported by a lot of researchers for a long period of time, making sometime very difficult for other scenarios to be commonly accepted and being published on archival journals. The conformal solution is easy to describe with the following sentences:

  • The gluon propagator goes to zero at lower momenta.
  • The ghost propagator goes to infinity at lower momenta faster than the free case.

Supporters of this solution also claimed that the running coupling reaches a fixed point at zero momenta. This fixed point was strongly lowered as lattice computations easily showed that things do not stay this way.

Indeed, lattice computations showed without doubt, increasing the volume, as the gluon propagator reaches a finite value at zero momenta and never converges toward zero. The ghost propagator is seen to behaves exactly as that of a free particle. But notwithstanding such an evidence, much effort is still devoted to understand why the conformal solution is not seen in such computations. A lot of hypotheses were put forward to explain why Yang-Mills propagators do not behave in such a well acquired way.

Axel Maas had a bright idea on this way to understand. He turned his attention to the D=2 model on the lattice (see here). Nobody did this before as this model is known to be trivial as it has no dynamics at all. This was proved by ‘t Hooft long ago. Maas showed that the 2D model on the lattice gives back the conformal solution. This means that a Yang-Mills model without dynamics has two-point functions behaving like those of the conformal solution. Most people saw this as a strong support to the conformal solution but things should be seen the other way round. Removing dynamics makes the conformal solution appear and lattice computations in 3 and 4d, that use the same code, are telling us the right behavior of the two-point functions. This is a serious setback for this approach and my personal view should be that any effort to further support this kind of solution should be  abandoned as is a loss of resources and precious time. But this could be a very difficult choice for a lot of people that for a long time worked on this and contributed to create this scenario.

In my opinion, the most severe drawback of this solution is that it goes against the original idea by Goldman and Haymaker that the gluon propagator should give back the Nambu-Jona-Lasinio model. Indeed, the conformal solution does not even give back a mass gap and it is practically useless to answer to a lot of open questions in the low-energy phenomenology. In some way it remembers to me the position  of the bootstrap model against the gauge paradigm that in the end proved to be the more fruitful approach. The last paper by Cucchieri and Mendes (see here) proves without doubt that the right form is a sum of Yukawa propagators and this gives back immediately a Nambu-Jona-Lasinio model.

As I cannot see reasons for people supporting the conformal solution to surrender, I think a lot of time will pass yet before truth will be acquired. Meantime, we stay here looking at the fight between such  strong contenders.


A formula I was looking for

September 23, 2008

As usual I put in this blog some useful formulas to work out computations in quantum field theory. My aim in these days is to compute the width of the \sigma resonance. This is a major aim in QCD as the nature of this particle is hotly debated. Some authors think that it is a tetraquark or molecular state while others as Narison, Ochs, Minkowski and Mennessier point out the gluonic nature of this resonance. We have expressed our view in some posts (see here and here) and our results heavily show that this resonance is a glueball in agreement with the spectrum we have found for a pure Yang-Mills theory.

Our next step is to understand the role of this resonance in QCD. Indeed, we have shown in our recent paper (see here) that, once the gluon propagator is known, it is possible to derive a Nambu-Jona-Lasinio model from QCD with all parameters properly fixed. We have obtained the following:

S_{NJL} \approx \int d^4x\left[\sum_q \bar q(x)(i\gamma\cdot\partial-m_q)q(x)\right.

-\frac{1}{2}g^2\sum_{q,q'}\bar q(x)\frac{\lambda^a}{2}\gamma^\mu q(x)\times

\left.\int d^4yG(x-y)\bar q'(y)\frac{\lambda^a}{2}\gamma_\mu q'(y)\right]

being G(x-y) the gluon propagator with all indexes in color and space-time already saturated. This in turn means that we can use the following formula (see my paper here and here):

e^{\frac{i}{2}\int d^4xd^4yj(x)G(x-y)j(y)}\approx {\cal N}\int [d\sigma]e^{-i\int d^4x\left[\sigma\left(\frac{1}{2}\partial^2\sigma+\frac{Ng^2}{4}\sigma^3\right)-j\sigma\right]}

being again G(x-y) the gluon propagator for SU(N) and {\cal N} a normalization factor. This formula does hold only for infrared limit, that is when the theory is strongly coupled. We plan to extract physical results from this formula and define in this way the role of \sigma resonance.


QCD and gluon propagator

September 8, 2008

As I have pointed out in a preceding post, I have got a paper of mine accepted for publication (see here) where I derive a Nambu-Jona-Lasinio model from QCD. I am able to fix the following parameters

G_{NJL}\approx 3.76 \frac{g^2}{\sigma}

and

\Lambda_{NJL}=\frac{\pi}{2K(i)}\sqrt{\sigma}

being as always \sigma=(440MeV)^2 the string tension and g the gauge coupling. Indeed, this model is very similar to the one presented in a paper by Kurt Langfeld, Christiane Kettner, and  Hugo Reinhardt (see here and here). They solved the model through the numerical solution of the Dyson-Schwinger equations. One can verify that in this case a mean field approximation does not work. In our case, with the above parameters, one gets

\frac{G_{NJL}\Lambda^2_{NJL}}{16\pi^2}\approx 1.117

when \alpha_s\approx 2.6 and we see that we have built a consistent model. Indeed, the authors of the aforementioned paper get the following table

These results are really striking. For \alpha_s=5.3 our model obtain a substantial agreement with experiments. In the near future we hope to find some analytical treatment to work out these values. At this stage we can be really satisfied of this achievement.

Indeed, this comparison is somewhat rough. We have to obtain the Dyson-Schwinger equation for our model and carry out the above computations by ourselves. This is a program to work out in the future. But already at this stage, the consistency with values of \alpha_s we obtain is really satisfactory.