text:
# Constituent quark masses from a modified perturbative QCD

## I Introduction

One of the great achievements of theoretical High Energy Physics in the last thirty years has been the discovering and development of the Quantum Chromodynamics (QCD). The smallness of the coupling constant at high momenta (asymptotic freedom) has allowed to develop a perturbative framework applicable for very high energy processes. That situation strongly simplifies the study of such phenomena and the calculated quantities are in good agreement with the experimental data. However, perturbative QCD (PQCD) is far from being able to furnish even a rough description for the relevant physics at low energies. The solution of this situation is currently one of the main challenges of Particle Physics.

A relevant phenomenon related with QCD is the color confinement. Nowadays, there are strong reasons to believe that the relation between the basic quantities in QCD, like the gluon and quarks fields, and the real world characterized by a whole variety of interacting mesons and baryons, can be understood by solving the confinement problem. The basic picture which seems plausible to be derived from QCD is that the fundamental fields can not be associated to asymptotic states of quark and gluons and then the true physical states consist only of colorless composites of quarks and gluons (mesons and baryons).

The above quoted limitations of the PQCD mean in particular that the usual Fock space vacuum of the non-interacting theory, is unable to predict, even approximately, the real ground state properties of QCD Savvidy1 ; Savvidy2 ; Savvidy3 ; Reuter . This is at contrast with the case in QED where the standard perturbation theory gives a better than good concordance with the experimental data. Then, the nature of vacuum structure is one of the main problems to be clarified and naturally its solution is closely linked with the color confinement effect. Reviews of various models considered to investigate the confinement and the vacuum structure in QCD can be found in Refs. ShuryakTex ; Shuryak2 .

Another basic problem for QCD is the one related with the quark masses, that is, the so called current quarks masses which appear in the Lagrangian. As the quarks are not observed as free states, the meaning of the quark masses need to be considered with care Gell ; Gross ; Weinberg1 . The nature of these parameters is an important object of research since long time ago Weinberg2 ; Weinberg3 ; Weinberg4 . The relations among these masses and their magnitude have been determined through the methods of current algebra in combination with measurements of the light mesons masses.

In the former works 1995 ; PRD ; tesis  an attempt directed to construct a modified perturbation expansion for QCD being able to predict at least some low energy properties of this theory have been considered. First, in 1995 , the modified expansion conserving the color $`SU(3)`$ and Lorentz symmetries was proposed, aiming to solve the symmetry limitations of the earlier chromo-magnetic field models Savvidy1 ; Savvidy2 ; Savvidy3 ; Reuter  which inspired the search. A similar lack of manifest Lorentz invariance had also a modification of the Feynman rules, intending to include gluon condensation, advanced later in Hoyer1 ; Hoyer2 ; Hoyer3  in which a delta function for $`k<\mathrm{\Lambda }_{qcd}`$ was summed to the perturbative piece. The expansion proposed in 1995 , considers a change of the gluon propagator in a term associated to a condensation of zero momentum gluons. This alteration of the usual rules had the interesting property of producing a non vanishing value for the gluon condensation parameter $`g^2G^2`$ already in the first approximations Zakharov . In addition, in 1995  a non vanishing value for the effective self-energy of the gluons was obtained at the loop level. Finally, a perturbative evaluation of the effective potential as a function of the condensate parameter indicated that the condensation is spontaneously generated from a zero condensate state.

The next work PRD ; tesis , had as a main objective to justify the applicability of the Feynman expansion introduced in 1995 . This search was in need because after modifying the propagator, it was unknown whether or not the initial wavefunction generating the expansion was a physical state of the free theory. In PRD ; tesis , by making use of the operational formulation of QCD developed by Kugo and Ojima Kugo , it was possible to find a physical state of the non-interacting theory being able to generate through the Wick expansion, the sort of propagators considered in 1995 . The discussion allowed also a more precise characterization of the class of changes admitted in the diagrammatic expansion by the physical state condition on the initial state. Specifically, it was found that the $`C`$ parameter, describing the gluon condensation, must be real and positive.

The present work is an expanded version of a letter sent to publication. Its main objective is to evaluate first $`G^2`$ in terms of the condensate parameter $`C`$up to order $`g^2`$, in order to determine the quark masses after fixing the parameter $`g^2C`$ to reproduce the accepted value for $`g^2G^2`$. The essential result is the evaluation of the constituent masses for the light quarks. In order to perform it, we have made use of the recent values reported for the current masses Report . The idea inspiring the work was the possibility that the gluon condensation described by modified expansion could predict the appearance of large masses for the light quarks ($`u`$, $`d`$ and $`s`$) of the order of 1/3 of the nucleon mass (constituent quark mass), thus furnishing an explanation for the successful constituent quark models. Surprisingly, after fixing the condensate parameter $`g^2C`$ to reproduce the most accepted value of the $`g^2G^2`$ (through the here performed evaluation of this quantity in terms of $`g^2C`$), it followed that the light quarks $`u`$, $`d`$ and $`s`$ got exactly the needed constituent like values for the masses in one of the propagation modes of the modified mass shell. The free propagation modes, which are at discrepancy with the absence of free quarks, are here only arising within an approximation evaluating the main effects of the condensate. After the inclusion of all the terms completing the one loop approximation in the gluon self-energy an interesting picture arises. It follows that the previously free propagation modes of the $`u`$ and $`d`$ quarks turn to be confined. On another hand, the $`s`$, $`c`$, $`b`$ and $`t`$ quarks waves become damped. That is, the squared rest mass get complex values. It should be stressed that the peaked structure of the propagator for the $`u`$ and $`d`$ at $`\sqrt{p^2}=m=0.30915GeV`$ signals the value of $`m`$ as the one expected to appear in the kinetic energy terms $`\stackrel{}{p}^2/(2m)`$ defining the free part of the bound state equations for $`u`$ and $`d`$ constituted nucleons. Henceforth, a possible mechanism contributing to justify the drastic difference between the stability of nucleons and the higher resonances is suggested. The absence of propagating modes for the $`u`$ and $`d`$, might be assuring the absolute stability shown by the nucleons bound states. In other manner, the presence of damped modes in the rest of the quarks would contribute to the decaying channels observed for all the hadron resonances associated to these short lived states.

Next, after considering the ground state within each of the baryon resonance groups as classified in Report  and its set of three constituent quarks, and adding the calculated quark masses for each of them (which depends on the flavor through the assumed values of the current masses), a spectrum of baryon masses was estimated. The results again reproduce reasonably well the experimentally determined resonance energies Report . This outcome suggests the less relevant role of the interaction energy among quarks in determining the masses of these hadrons.

At this point it is important to make clear that earlier works Munczek ; Burden  introduced a pure delta function at zero momentum in the gluon propagator, searching for a model of meson resonances. In these treatments the constants introduced, multiplying the delta function, were different to the one fixed in our previous work PRD ; tesis , leading also to very different results. In particular, the singularity structure of the quark propagator have no poles on the real $`p^2`$ axis Burden  at difference with our results. This difference is directly connected with the definite sign for the coefficient of the $`\delta `$-function which in our case was precisely determined in PRD ; tesis . Since our analysis in that paper was based in the direct construction of the state incorporating the condensate, it follows that the difference with the results of Munczek ; Burden  is rooted in a corresponding difference in the initial state determining the Wick expansion. It could be also possible that both procedures lead to the same exact physical predictions, but the first approximations in each approach could turn to be more appropriate for the discussion of specific physical issues. This question will be addressed elsewhere. It should be also further noticed that gluon condensate states reproducing the modified expansion of Munczek ; Burden  have been also discussed in celenza ; pavel . However, all these references are related with states showing confined quarks and massive gluons, at variance with our proposed state determining massive quarks (in the simplest approximation) and tachyonic gluons.

It is also following that another dispersion relation is obtained which predicts a vanishing value of the mass as $`m_Q0`$. This mode could be connected with the family of low mass mesons (e.g. $`\mathrm{\Pi }`$ and $`K`$ mesons). Specifically, the possibility is open that the bound states of quark excitations in these light mass states of quarks could describe such low lying mesons within the here presented framework. A circumstance pointing in this direction is that the spontaneous chirality breaking associated with this mass generation should be expected to be associated with some massless Goldstone modes. Therefore, it is possible that the bound states of such modes predicted by the Bethe-Salpeter equations turn to describe the light mesons $`\mathrm{\Pi }`$ and $`K`$. This question however, will be deferred to a further study.

The calculations of $`G^2`$ and the quark and gluon masses will be performed up order $`g^2`$. It seems worth anticipating that the approximation scheme will resemble a sort of quasi-classical limit in which the condensate seems to play a role of a macroscopic mean quantum field. Such an interpretation is in concordance with the initial motivation of our analysis as directed to construct a covariant version of the chromo-magnetic models 1995 .

Another conclusion which should be noticed is that the evaluated gluon mass determined from the allowed values of the condensate parameter $`g^2C`$ is tachyonic. The possibility for a tachyonic gluon mass has been recently discussed in the literature tachyon1 ; tachyon2 ; Hoyer3 . This result indicates that the free propagator of the expansion is not reproducing itself after the first approximation. This property, however, does not mean a limitation of the approach. This is because, at variance with Munczek ; Burden , we are not assuming that the exact propagator has the same delta function structure. This form only arose for the free propagator as determined by the Wick expansion around the considered initial state.

The presentation is organized as follows: Section II is devoted to make precise the elements characterizing the modified perturbation expansion. In the Section III the calculation of the $`G^2`$ up to $`g^2`$ order for gluodynamics is presented and the result used to estimate the condensation parameter $`g^2C`$. Further in Section IV the effects of the condensate on the quark masses are evaluated and the obtained values used to estimate the spectrum of the ground states in each baryon group as classified in Report . In Section V the calculation of the effect of the condensate on the one loop self-energy for the gluons is evaluated. Finally, the results are summarized and some open questions and a conjecture are underlined in the last Section VI.

## II The Modified Feynman Rules

In this introductory section, the main elements and conventions of the modified expansion introduced in 1995 ; PRD ; tesis  will be reviewed. The generating functional of the Green functions in the modified theory has the form

$`Z[J,\overline{\eta },\eta ,\overline{\xi },\xi ]`$ $`=`$ $`{\displaystyle \frac{1}{N}}{\displaystyle D(A^r,\overline{c}^r,c^r,\overline{\mathrm{\Psi }}^r,\mathrm{\Psi }^r)\mathrm{exp}\left\{id^4x^{sources}\right\}},`$ (1)
$`^{sources}`$ $`=`$ $`+J^{\mu ,a}A_\mu ^{ra}+\overline{c}^{ra}\eta ^a+\overline{\eta }^ac^{ra}+\overline{\mathrm{\Psi }}^{ri}\xi ^i+\overline{\xi }^i\mathrm{\Psi }^{ri},`$ (2)
$`N`$ $`=`$ $`{\displaystyle D(A^r,\overline{c}^r,c^r,\overline{\mathrm{\Psi }}^r,\mathrm{\Psi }^r)\mathrm{exp}\left\{id^4x\right\}},`$ (3)
where sources for all the renormalized fields $`(A_\mu ^{ra},\overline{c}^{ra},c^{ra},\overline{\mathrm{\Psi }}^{ri},\mathrm{\Psi }^{ri})`$ have been introduced in the usual manner, and the effective Lagrangian for the bare fields was selected in the Kugo-Ojima quantization procedure for gauge theory as

$``$ $`=`$ $`_G+_{Gh}+_Q`$ (4)
$`_G`$ $`=`$ $`{\displaystyle \frac{1}{4}}F_{\mu \nu }^a\left(x\right)F^{\mu \nu ,a}\left(x\right){\displaystyle \frac{1}{2\alpha }}\left(^\mu A_\mu ^a\left(x\right)\right)^2,`$ (5)
$`_{Gh}`$ $`=`$ $`i^\mu \overline{c}^a\left(x\right)D_\mu ^{ab}\left(x\right)c^b\left(x\right),`$ (6)
$`_Q`$ $`=`$ $`\overline{\mathrm{\Psi }}^i\left(x\right)\left(i\gamma ^\mu D_\mu ^{ij}m_Q^{}\delta ^{ij}\right)\mathrm{\Psi }^j\left(x\right).`$ (7)
The sum over the six quark flavors will be omitted everywhere in order to simplify the exposition. We think that no confusion should arise for it. The gluon field intensity has the usual form

$$F_{\mu \nu }^a\left(x\right)=_\mu A_\nu ^a\left(x\right)_\nu A_\mu ^a\left(x\right)+g^{}f^{abc}A_\mu ^b\left(x\right)A_\nu ^c\left(x\right),$$

where $`D_\mu ^{ab},D_\mu ^{ij}`$ are the covariant derivatives in the adjoint and fundamental representations of the group $`SU(3)`$ respectively

$`D_\mu ^{ab}\left(x\right)`$ $`=`$ $`_\mu \delta ^{ab}g^{}f^{abc}A_\mu ^c\left(x\right),`$
$`D_\mu ^{ij}\left(x\right)`$ $`=`$ $`_\mu \delta ^{ij}ig^{}T^{ij,a}A_\mu ^a.`$
Here, the bare fields and bare coupling constant and masses ($`g^{},m_Q^{}`$), were expressed in terms of their renormalized counterparts following the same conventions of Ref. Muta . Below, the renormalized masses and coupling will be denoted as $`m_Q`$ and $`g`$ respectively.

It should be recalled that as argued in PRD ; tesis  the physical predictions for the value $`\alpha =1`$ of the gauge parameter should have physical meaning whenever the adiabatic connection of the interaction does not lead the evolving state out the physical space. Therefore here we will fix this value of the parameter.

Then, after the standard procedure of extracting out of the functional integral the exponential of the terms of higher than second order in the fields (vertices) the generating functional takes the form to be employed in the evaluations below Faddeev :

$`Z[J,\overline{\eta },\eta ,\overline{\xi },\xi ]=N^1\mathrm{exp}\{i[{\displaystyle \frac{S_{abc}^G}{3!i^3}}{\displaystyle \frac{\delta ^3}{\delta J_a\delta J_b\delta J_c}}+{\displaystyle \frac{S_{abcd}^G}{4!i^4}}{\displaystyle \frac{\delta ^4}{\delta J_a\delta J_b\delta J_c\delta J_d}}`$
$`+{\displaystyle \frac{S_{ras}^{Gh}}{2!i^3}}{\displaystyle \frac{\delta ^3}{\delta \overline{\eta }_r\delta J_a\delta \left(\eta _s\right)}}+{\displaystyle \frac{S_{iaj}^Q}{i^3}}{\displaystyle \frac{\delta }{\delta \overline{\xi }_i\delta J_a\delta \left(\xi _j\right)}}]\left\}Z_0\right[J,\overline{\eta },\eta ,\overline{\xi },\xi ],`$ (8)
with the free functional given by

$$Z_0[J,\overline{\eta },\eta ,\overline{\xi },\xi ]=\mathrm{exp}\left\{i\frac{J_aG_G^{ab}J_b}{2}+i\overline{\eta }_rG_{Gh}^{rs}\eta _s+i\overline{\xi }_iG_Q^{ij}\xi _j\right\},$$

and $`N`$ normalizes $`Z`$ to be one at vanishing sources. Use will be made of the DeWitt compact notation Daemi , in which a Latin letter let say $`a,b,\mathrm{}`$, for a field symbolizes the space-time coordinates of it as well as all its internal quantum numbers. The same index in a source or a tensor symbolizes the same set of variables of the kind of fields associated with this specific index. For example $`A_a=A_{\mu _a}^a(x_a),c_a=c^a(x_a)\text{and}\mathrm{\Psi }^i=\mathrm{\Psi }^i(x_i)`$. Note that a particular convention has been also employed in which the same letter is employed for the global index (for all the coordinated and internal indices) and the one indicating the internal coordinates in the explicit form of the quantities. Such a procedure was useful for the calculations and we also think that will not create confusion. As usual, repeated indices represent the corresponding space time integrals and contracted, Lorentz, spinor or color components.

The following definitions also have been used in (8):

$`S_{ijk\mathrm{}}^\alpha `$ $``$ $`\left({\displaystyle \frac{\delta }{\delta \mathrm{\Phi }_i}}{\displaystyle \frac{\delta }{\delta \mathrm{\Phi }_j}}{\displaystyle \frac{\delta }{\delta \mathrm{\Phi }_k}}\mathrm{}S^\alpha \left(\mathrm{\Phi }\right)\right)_{\mathrm{\Phi }=0}\text{ for }\alpha =G,Gh,Q\text{ and }\mathrm{\Phi }=A,\overline{c},c,\overline{\mathrm{\Psi }},\mathrm{\Psi }`$
$`G_\alpha ^{ij}`$ $``$ $`S_{\alpha ,ij}^1\text{ for }\alpha =G,Gh,Q,`$
where $`G`$, $`Gh`$ and $`Q`$ mean the gluon, ghost and quark parts of the action respectively.

There is only one main element determining the difference of the usual Perturbative QCD expansion and the modified one considered in 1995 ; PRD ; tesis . It is related with the form of the gluon free propagator. As proposed in 1995  there is an additional term to the gluon free propagator which is absent in the standard expansion. In tesis ; PRD , such a term was shown to be a consequence of a Wick expansion based in a state constructed by acting on the usual vacuum with an exponential of pairs of zero momentum gluon and ghost creation operators. The ghost propagator PRD ; tesis  could remain unmodified if the parameter left after fixing the form of the wave-function is taken as a real and positive one PRD ; tesis . Such a selection will be employed here. The quark propagator is not affected in any way by the gluon condensation as introduced in 1995 ; PRD ; tesis . However, as it will be commented in the summary, the results of the present work led to the idea about a possible physical relevance of introducing quark condensates along the similar lines as it was done for the gluons. We consider the exploration of this possibility to be exposed as one of the most interesting extensions of the work.

A technical point can be noticed for precision. The operator quantization employed in PRD ; tesis  takes the ghost fields as satisfying the conjugation properties proper of the approach of Kugo and Ojima Kugo . However, at the level of the Feynman diagram expansion the difference with the standard procedure is only a change of variables.

In accordance with above remarks and the results of 1995 ; PRD ; tesis  the diagram technique defining the modified expansion has the following basic propagators for the renormalized fields

$`G_G^{ab}`$ $`=`$ $`\delta ^{ab}g^{\mu _a\mu _b}\left[{\displaystyle \frac{1}{k^2+i\epsilon }}iC\delta \left(k\right)\right],`$
$`G_{Gh}^{rs}`$ $`=`$ $`\delta ^{rs}{\displaystyle \frac{\left(i\right)}{k^2+i\epsilon }},`$
$`G_Q^{ij}`$ $`=`$ $`\delta ^{ij}{\displaystyle \frac{m+p^\mu \gamma _\mu }{\left(m^2p^2i\epsilon \right)}}.`$
and the same normal and counterterm vertices of the usual theory Muta . The dimensional regularization and minimal substraction (MS) is the renormalization procedure employed here.

Being defined the Feynman rules and the conventions to be used, below they will be applied to the calculation of various quantities of interest.

## III The mean value of $`G^2`$

The calculation of the mean value of the gluon field intensity squared is considered in this section. It extends the former evaluation done in 1995  within a simpler approximation. The result will be employed afterwards for fixing the value of the condensation parameter $`g^2C`$ to be the necessary one for furnishing the currently estimated value for $`g^2G^2`$. The calculation will be performed up to order $`g^2`$ including all the dependence on the condensate.

It is a helpful circumstance here, that the one loop renormalization procedure of the modified expansion can be performed identically as in the standard PQCD. This is a direct consequence of the fact that the $`\delta `$-functions kill the only existing integrals at the one loop level and therefore all the corresponding counterterms in the dimensional regularization are identical. This fact will be helpful in regularizing and renormalizing the composite operator $`G^2`$.

For the evaluation of $`G^2`$ the following definition of this quantity will be employed

$$G^2=\frac{0\left|S_g\right|0}{𝑑x^D}=\frac{D(A^r,\overline{c}^r,c^r,\overline{\mathrm{\Psi }}^r,\mathrm{\Psi }^r)S_g\left[A^r\right]\mathrm{exp}\left\{id^4x^{sources}\right\}}{𝑑x^DD(A^r,\overline{c}^r,c^r,\overline{\mathrm{\Psi }}^r,\mathrm{\Psi }^r)\mathrm{exp}\left\{id^4x\right\}},$$
(9)
where $`S_g\left[A^r\right]`$ represents the gluon part of the action in the absence of the gauge breaking term depending on the gauge parameter $`\alpha `$, evaluated at the renormalized gluon fields $`A^r`$, that is

$$S_g\left[A^r\right]=\frac{1}{4}d^4xF_{\mu \nu }^{ra}\left(x\right)F^{r\mu \nu ,a}\left(x\right).$$

Therefore, $`G^2`$ will be considered as defined by the mean value of the gluon action divided by the space-time volume in dimensional regularization.

After substituting in (9) the perturbative expression for the mean value in terms of the Generating Functional (8) it follows

$$0\left|S_g\right|0=\left\{\left[\frac{S_{ab}^g}{2i^2}\frac{\delta ^2}{\delta J_a\delta J_b}+\frac{S_{abc}^g}{3!i^3}\frac{\delta ^3}{\delta J_a\delta J_b\delta J_c}+\frac{S_{abcd}^g}{4!i^4}\frac{\delta ^4}{\delta J_a\delta J_b\delta J_c\delta J_d}\right]Z[J,\overline{\eta },\eta ,\overline{\xi },\xi ]\right\}_{J,\overline{\eta },\eta ,\overline{\xi },\xi =0}.$$
(10)
In the momentum representation the kernel $`S_{ab}^g`$ has the form

$$S_{ab}^g=\left(2\pi \right)^4\delta \left(k_a+k_b\right)\delta ^{ab}\left[\left(g_{\mu _a\mu _b}k_a^2k_{a,\mu _a}k_{a,\mu _a}\right)\right],$$

and the $`S_{abc}^g`$ and $`S_{abcd}^g`$ ones are the usual three and four legs gluon vertices Muta .

Let us consider below the evaluation of $`G^2`$ up to order $`g^2`$. As any Feynman diagram can be expressed as a polynomial in the constant $`C`$, after associating a new diagrammatic line for this part of the free propagator, all the Feynman diagrams can be divided in two groups. The first, let say $`G_1`$ including the diagrams not showing the new propagator line and another group, $`G_2`$ having one or more, let us call them, “condensate” lines. Such lines have associated to them a multiplicative Dirac delta function at zero momentum.

After writing out all the expressions of type $`G_2`$ (up to second order in the renormalized coupling $`g^2`$) contributing to (10), it follows that the only non-vanishing diagram is the one appearing in the Figure 1 in which both gluon lines are substituted by “condensate” lines. The remaining terms vanish thanks to the zero momentum evaluation produced by the Dirac Delta functions introduced by the “condensate” lines and also due to the vanishing of tadpoles in dimensional regularization. Also, it should be underlined that the retention of the $`iϵ`$ factors defining the Feynman propagator, up to the evaluation of the integrals was also important in obtaining the zero result for some of the diagrams in $`G_2`$.

The contribution related with Figure 1 had been calculated before and the result is 1995 ; PRD

$$T_{G_2}=\frac{1}{4}\frac{288g^2C^2}{\left(2\pi \right)^8}.$$

Next, the diagrams in $`G_1`$ are the ones contributing to $`G^2`$ up to order $`g^2`$ in standard perturbative QCD. They are illustrated in Figure 2. The terms $`T_6`$ and $`T_7`$ represent the only non-vanishing and non-tadpole like second order contributions. It can be noticed that the appearance of two tadpoles in the term $`T_7`$ leads to the exact vanishing of this contribution in dimensional regularization.

For the diagram $`T_6`$ corresponds the expression

$$E_{T_6}=\frac{1}{3!i^3}\frac{dp^D}{(2\pi )^Di}\frac{g^{\mu \nu }\delta ^{ab}}{p^2+i\epsilon }\mathrm{\Pi }_{\mu \nu }^{(g)ab}(p),$$

where $`\mathrm{\Pi }^{(g)}`$ is the gluonic contribution to the polarization tensor in second order, which has been explicitly calculated. After substituting the known result for $`\mathrm{\Pi }`$ and performing the Wick rotation, the previous expression turns out to be of the form

$$E_{T_6}\frac{dq^D}{(2\pi )^D}(q^2)^{\frac{D}{2}2}$$
(11)
Then considering the formula narison

$`{\displaystyle \frac{dq^D}{(2\pi )^D}\frac{1}{(q^2)^\alpha }\frac{1}{(qk)^{2\beta }}}`$ $`=`$ $`{\displaystyle \frac{i(1)^{\alpha \beta }(k^2)^{\alpha \beta +n/2}}{(16\pi ^2)^{n/4}}}{\displaystyle \frac{\mathrm{\Gamma }(\alpha +\beta n/2)}{\mathrm{\Gamma }(\alpha )\mathrm{\Gamma }(\beta )}}B({\displaystyle \frac{n}{2}}\beta ,{\displaystyle \frac{n}{2}}\alpha ),`$
$`B(n,m)`$ $`=`$ $`{\displaystyle \frac{\mathrm{\Gamma }(m)\mathrm{\Gamma }(n)}{\mathrm{\Gamma }(m+n)}},`$
for the case $`\alpha +\beta =2\frac{D}{2}`$ and fixing $`k_\mu =0`$, it follows the vanishing of diagram $`T_6`$ within the dimensional regularization for $`D>2`$.

It rest now to discuss the contributions of the tadpole like diagrams. It is worth noticing at this point that they are quite similar in their form to the one corresponding to the square of the electromagnetic field strength in QED. In fact they differ from it only in some numerical color factors. An important quantitative difference between the expressions for both theories is related with magnitude of the coupling constant which is very much higher in QCD. Let us consider the $`T_1T_5`$ terms in Figure 2 embodying the gluon, ghost and fermion contributions with the inclusion of the one loop counterterms rendering finite the polarization operator at one loop Muta . This expression turns to be of the form

$$E_{T_1T_5}\frac{dp^D}{(2\pi )^D}\frac{(p_\mu p_\nu p^2g_{\mu \nu })(g^{\mu \nu }\frac{p^\mu p^\nu }{p^2+i\epsilon })}{p^2+i\epsilon }\delta ^{aa}\left(1\mathrm{\Pi }(p^2)\right),$$
(12)
where $`\mathrm{\Pi }`$ as usual is defined through the polarization tensor

$$\mathrm{\Pi }_{\mu \nu }^{ab}(p)=\delta ^{ab}(p_\mu p_\nu p^2g_{\mu \nu })\mathrm{\Pi }(p^2),$$

and it can be expressed as the superposition of the pure gluodynamical contribution and the quark one as Muta

$`\mathrm{\Pi }(p^2)`$ $`=`$ $`\mathrm{\Pi }_G(p^2)+\mathrm{\Pi }_Q(p^2)`$ (13)
$`=`$ $`{\displaystyle \frac{g^2\mu ^{4D}}{4\pi ^{D/2}}}(p^2)^{2D/2}\mathrm{\Gamma }(2{\displaystyle \frac{D}{2}})B({\displaystyle \frac{D}{2}},{\displaystyle \frac{D}{2}})\times `$
$`\times \left({\displaystyle \frac{C_G(3D2)}{\frac{D}{2}1}}+{\displaystyle \frac{8T_R}{B(\frac{D}{2},\frac{D}{2})}}{\displaystyle \underset{Q}{}}{\displaystyle _0^1}𝑑xx(1x)\left(x(1x){\displaystyle \frac{m_Q^2i\epsilon }{p^2}}\right)^{(2D/2)}\right)`$
$`+(Z_31),`$
where $`C_G=3`$ and $`T_R=1/2`$ for $`SU(3)`$ and $`\mu `$ is the scale mass in the dimensional regularization.

Considering first the pure gluonic contributions associated $`(1\mathrm{\Pi }_G(p^2))`$ in (12), it follows that they reduce to a sum of dimensional regularized integrals of the same form as (11). Thus, the second order gluodynamical contribution of the tadpole like contributions vanish. Let us analyze now the quark contributions associated to $`\mathrm{\Pi }_Q`$. It can be noticed that for massless quarks, that is $`m_Q=0`$, the $`\mathrm{\Pi }_Q`$ has the same analytical dependence of $`p^2`$ as the gluonic part $`\mathrm{\Pi }_G`$ and all the quark contributions vanish after the regularization is employed. However, for massive quarks the result for the second order dimensionally regularized contribution along the same lines employed for the previous terms gives rise to a non-vanishing contribution of the following form after the substraction of the divergent pole terms at $`D=4`$. The expression for the regularized term before substraction is

$$g^2\mu ^{4D}E_{T_6}=\left(\underset{Q}{}(m_Q^2)^{D2}\right)\frac{64\sqrt{\pi }g^4\mu ^{82D}T_R(D1)\mathrm{\Gamma }(2D)\mathrm{\Gamma }(4D)}{2^{72D}(4\pi )^D\mathrm{\Gamma }(\frac{9}{2}D)}.$$

This expression has simple as well as double poles which after substracted (in the $`D=4`$ limit following the rules of the MS scheme for composite operators collins ) give a finite correction to $`g^2G^2`$ for the gluon condensate parameter. It can be evaluated and the numerical expression becomes

$$\left(g^2\mu ^{4D}E_{T_6}\right)^{fin}=0.00064T_Rg^4\left(\underset{Q}{}m_Q^4\left[55.3237.4\mathrm{ln}\left(\frac{\mu ^2}{m_Q^2}\right)+6.0\mathrm{ln}^2\left(\frac{\mu ^2}{m_Q^2}\right)\right]\right).$$

However, it should be noticed that this finite expression has been obtained for the second order correction of the polarization operator. The summation of the infinite ladder of self-energy insertions leading to the one loop propagator can be suspected to produce a noticeable modification of this result. After also considering the high indefiniteness in the quark current masses, and the technical difficulty we had in evaluating the dimensionally regularized expression for the above mentioned self-energy insertions, we decided to employ for the further consideration here the purely gluodynamical contribution. In support of this way of proceeding there is also another circumstance. It is needed to notice that a refined version of the definition of the composite operator $`G^2(x)`$ could be suspected to exist in such a way that the terms associated to a tadpole diagram of the exact propagator vanish. This property would assure the exact vanishing of the $`F^2`$ in QED a property that to our knowledge, is expected. Because, at one loop, both QED and QCD tadpole expressions are basically the same, the above written quark contribution would not really appear. The existence of the mentioned definition, opens the possibility for the second order gluodynamical contribution to be exact one for QCD in this approximation. However, we shift the discussion of this question to further analysis to be considered elsewhere and simply adopt the gluodynamical result.

Therefore, the total second contribution to the mean value of $`G^2`$ in gluodynamics is defined by the only non-vanishing term

$$G^2=\frac{288g^2C^2}{\left(2\pi \right)^8},$$

coming from the condensate contributions arising from the diagram $`T_6`$ which is defined by the four gluon interactions. Finally, let us determine in a phenomenological way a value for the constant $`g^2C`$ by making use of a currently estimated value for $`g^2G^2`$:

$$g^2G^20.5\left(GeV/c^2\right)^4.$$

Therefore, the parameter $`g^2C`$ takes the value:

$$g^2C=64.9394\left(GeV/c^2\right)^2.$$
(14)
This fixing of the condensation parameter allows to investigate the physical predictions of the modified expansion in the next section.

## IV Constituent Quarks masses

Let us consider now the effect on the current quark masses produced by the condensate. The full inverse propagator for the quarks in terms of the corresponding mass operator has the general expression

$$G_Q^1\left(p\right)=\left(m_Qp^\mu \gamma _\mu \mathrm{\Sigma }\left(p\right)\right),$$
(15)
where the color index identity matrix has been omitted and $`\mathrm{\Sigma }\left(p\right)`$ is the mass operator. Its expression is related with the diagram in Figure 3.

The form for $`\mathrm{\Sigma }\left(p\right)`$ up to order $`g^2`$, within the modified Feynman rules can be written as

$$\mathrm{\Sigma }\left(p\right)=g^2C_F\frac{d^4k}{\left(2\pi \right)^4i}\frac{\gamma _\mu \left(m_Q+\left(pk\right)^\alpha \gamma _\alpha \right)\gamma _\nu G_G^{\mu \nu }\left(k\right)}{\left(m_Q^2\left(pk\right)^2\right)},$$
(16)
in which the color factor $`C_F`$ for the $`SU\left(N\right)`$ group is $`C_F=\frac{N^21}{2N}`$ which for the case of interest here $`SU\left(3\right)`$, reduces to $`C_F=\frac{4}{3}`$ Muta .

Taking into account the modified gluon propagator and the standard relations for the $`\gamma `$ matrices (which here are adopted in the same way as in Ref. Muta ) the one loop mass operator expression is simplified to be

$$\mathrm{\Sigma }\left(p\right)=g^2C_F\frac{d^4k}{\left(2\pi \right)^4i}\frac{2\left(2m_Q\left(pk\right)^\alpha \gamma _\alpha \right)}{\left(m_Q^2\left(pk\right)^2\right)}\left(\frac{1}{k^2}iC\delta \left(k\right)\right).$$

Let us disregard for the moment the term not involving the condensate parameter $`C`$ in order to first study the sole influence on the quark masses of the condensation effect. Then, the expression for $`\mathrm{\Sigma }`$ reduces to

$$\mathrm{\Sigma }\left(p\right)=\frac{g^2C_FC}{\left(2\pi \right)^4}\frac{2\left(2m_Qp^\alpha \gamma _\alpha \right)}{\left(m_Q^2p^2\right)}.$$

Therefore, the inverse propagator for the quark takes the simple form

$$G_Q^1\left(p\right)=m_Q\left(1+2\frac{M^2}{\left(m_Q^2p^2\right)}\right)p^\mu \gamma _\mu \left(1+\frac{M^2}{\left(m_Q^2p^2\right)}\right),$$
(17)
where the constant $`M`$ is defined by

$$M^2=\frac{2g^2C_FC}{\left(2\pi \right)^4}=0.1111\left(GeV/c^2\right)^2,$$

in which the numerical value has been obtained by using (14).

The zeros of the determinant associated to the inverse propagator (17) then allows to determine the effects of the condensate (reflected by the parameter $`M`$), on the effective mass of the quarks. Let below $`s,s=1\mathrm{}6`$ an index indicating each kind of quark characterized by its particular current mass $`m_{Q_s}`$.

In order to find the dependence of the effective quark mass $`m_{q_{s,i}}^2`$ as a function of the current mass parameters $`m_{Q_s}`$ an analysis was done for the zeros of the determinant of the inverse green function (17). It follows that there is only one solution having a squared mass being positive for arbitrary values of the current quark masses. The existence of this real and positive solution for the squared mass becomes possible thanks to the real and positive character of $`C`$, as was shown previously PRD ; tesis . There is another interesting branch of the zeros of the determinant (dispersion relations) showing vanishing values for the effective mass in the light current mass limit. For the purely real solution the value of the effective quark mass $`m_{q_s}=f\left(m_{Q_s}\right)`$ as a function of $`m_{Q_s}`$ is shown in the Figure 4. The graph is plotted for the region $`m_{Q_s}<2GeV/c^2`$ which contains the current mass values of the $`u,d,s`$ and $`c`$ quarks.

As it can be appreciated in the picture, the calculated effective quark masses for light flavors ($`u,d`$ and $`s`$) are clearly predicting the values of the quark masses being in use in the constituent quark models of the baryons. That is, the light quarks get a weight of near one third of the nucleon mass. In the table 1 the mass values obtained for each of the quark flavors are shown Report . From the global properties of table 1 it can be observed that the main effect of the gluon condensate seems to be the dressing the light quarks with a cloud gluons having a total mass of one third of the nucleon mass. Then, the results point in the direction of the mainly glue nature of the constituent quark masses of the $`u,d`$ and $`s`$ quarks within many baryon resonance. These results then support the idea that a modified perturbative expansion like the one being considered, in which the effect of the gluon condensate have been incorporated, would be able to predict with reasonably good approximation properties of the low energy strong interactions. Similar values for the constituent quarks masses have been obtained by different methods in Steele1 ; Steele2 .

We will now investigate the influence on the mass spectrum produced by including the standard one loop contribution in the dispersion relations. After substituting the known narison  quark one-loop self-energy expression, the one loop inverse quark propagator can be written as

$`G_Q^1\left(p\right)`$ $`=`$ $`\left(m_Qp^\mu \gamma _\mu \mathrm{\Sigma }\left(p\right)\right)=A(p^2)m_QB(p^2)p^\mu \gamma _\mu `$
$`=`$ $`m_Q[1+{\displaystyle \frac{4C_Fg^2C}{(2\pi )^4}}{\displaystyle \frac{1}{m_Q^2p^2}}{\displaystyle \frac{C_Fg^2(\mu )}{(4\pi )^2}}\times `$
$`(6+4\gamma \mathrm{ln}(4\pi )+4\mathrm{ln}\left({\displaystyle \frac{m_Q^2p^2}{\mu ^2}}\right)4{\displaystyle \frac{m_Q^2}{p^2}}\mathrm{ln}\left({\displaystyle \frac{m_Q^2p^2}{m_Q^2}}\right))]`$
$`p^\mu \gamma _\mu [1+{\displaystyle \frac{2C_Fg^2C}{(2\pi )^4}}{\displaystyle \frac{1}{m_Q^2p^2}}+{\displaystyle \frac{C_Fg^2(\mu )}{(4\pi )^2}}\times `$
$`(1\gamma +\mathrm{ln}(4\pi )+{\displaystyle \frac{m_Q^2}{p^2}}\mathrm{ln}\left({\displaystyle \frac{m_Q^2p^2}{\mu ^2}}\right)+({\displaystyle \frac{m_Q^2}{p^2}})^2\mathrm{ln}\left({\displaystyle \frac{m_Q^2p^2}{m_Q^2}}\right))],`$
where the dimensionless coupling constant $`g^2`$ has been substituted by its running expression as a function of the scale parameter $`\mu `$ as

$$g(\mu )^2=\frac{4\pi 0.12}{1+\frac{\mathrm{7\; 0.24}}{4\pi }\mathrm{ln}\left(\frac{\mu }{91.18}\right)}.$$

That is, the value of $`\alpha ^2(\mu )=g^2(\mu )/(4\pi )`$ is fixed to be $`0.12`$ at the scale of the $`Z`$ mass $`\mu =m_Z=91.18GeV`$. In order to automatically fix the scale $`\mu `$ to the one associated to the magnitude of the obtained solution for the mass, the value of $`\mu `$ was taken to satisfy $`\mu ^2=p^2`$. That is, to reproduce the same value of the solution for the squared mass.

Before proceeding, it should be noticed that the inverse propagator, after the introduction of the one loop terms, gets a branch cut at the real axes for the considered complex variable $`p^2`$ for Re$`(p^2)>m_Q^2`$. However, the $`iϵ`$ prescription selects one of the sides of the cut as giving a well defined value of the inverse propagator at real $`p^2`$ values.

After multiplying (IV) by $`A(p^2)m_QB(p^2)p^\mu \gamma _\mu `$ the modified mass shells can be most easily found by solving

$$D(q^2)=A^2(p^2)m_Q^2B^2(p^2)p^2=0.$$

Concretely, we estimated the corrected values of the masses in Table 1 by searching for the peaks of the inverse for the absolute value of $`D(q^2)`$ as a 2D-function of the real and imaginary parts of $`p^2`$. The new values are reported in Table 2. As it can be observed the masses were not drastically changed. They all show a relatively small reduction in their magnitude, then indicating the consistency of the predicted constituent masses for the light quarks at the one loop level.

An interesting outcome following from the inclusion of the standard one loop terms, should be stressed. It consists, in that for the smallest $`m_Q`$ masses associated to the lower and higher bounds for the $`u`$ and $`d`$ quarks, $`\frac{1}{D(q^2)}`$ although showing a definite peak at $`0.30915GeV`$ does not really diverge near this value. Therefore, the propagator is not showing a definite pole associated to a propagating quark mode at this $`p^2`$ value. However, at the level of the Bethe-Salpeter equation for bound states of quarks, it should be expected that $`p^2`$ at this peak should behave as the non-relativistic value for the square of the mass entering the kinetic energy expression $`\frac{p^2}{2m_Q}`$. Henceforth, in spite of the presence of free propagating modes in the simplest approximation, the inclusion of the quantum one loop correction makes the $`u`$ and $`d`$ modes non-propagating (confined) in the same way as it occurs in the approach of Munczek and Nemirovsky Munczek ; Burden .

For greater values of $`m_Q`$ in the range presently estimated for the $`s`$ quark current mass ($`0.0600.170GeV`$) in the neighborhood of the non-diverging peak of $`\frac{1}{D(q^2)}`$, this quantity starts showing divergency points at complex values of $`p^2`$ which represent damped but propagating quark modes. The just described structure of the spectrum suggests a possible explanation for the drastic difference of stability between the $`u`$ and $`d`$ composed hadrons and the resonances including $`s`$, $`c`$, $`b`$ and $`t`$ quarks. The point is that the presence of poles at complex $`p^2`$ for the heavier quarks could lead to the existence of disintegration channels for the higher mass resonances in which the $`s`$, $`c`$, $`b`$ and $`t`$ quarks can escape in a damped propagating mode to be after hadronized by some separate mechanism. On another hand, as the $`u`$ and $`d`$ quarks do not show propagating states, that could remain in a perennial confinement within an undisturbed nucleon. A further study of the consistency of this explanation is expected to be done within the planned extension of this work to investigate bound state equations for hadrons.

It can be concluded that the constituent quark masses evaluated from the simplest approximation for the dispersion relation considered above are stable under the inclusion of the standard one loop corrections to the self-energy.

Next, in the Table 3, the masses of the ground state resonance selected each of them from one of the groups of baryons as classified in Report , are shown in comparison with the values obtained for them by a theoretical estimate consisting in adding the here calculated masses for each of their known quark constituents. The quark masses of the Table 1 were employed for this evaluation. Also, in the Table 4, the masses for some vector mesons reported also in Report  are shown in conjunction with the result of the sum of the calculated masses for their constituent quarks (Table 1).

As it can be appreciated the results obtained for the vector mesons and baryon resonance are only estimates. However, the obtained values are in a reasonable well correspondence with the reported experimental values. It can be also noticed that the validity of the present calculations implies that the binding energy contribution to the baryon rest masses should be small. This property then could justify a non relativistic description of such resonances.

As it was remarked before, another form of the dependence of the effective quark mass with $`m_{Q_s}`$ is also following from the dispersion relations. It, by contrast shows a vanishing effective mass when $`m_{Q_s}0`$. The existence of this branch opens the possibility for the existence an appropriate description of the light mesons ($`\mathrm{\Pi }`$ and $`K`$) within the proposed scheme. They could be associated to bound states of light quarks propagating in such modes. The analysis of this question deserves a separate treatment to be examined elsewhere. The graphic for the mentioned branch is represented in Figure 5, for the region of $`m_{Q_s}<0.1GeV/c^2`$, a point at which the change in the analytic dependence signals where the solution gets an imaginary part when the current mass value $`m_Q`$ grows.

## V  Tachyonic gluon Self-Energy

Let us consider now the evaluation of the gluon self-energy. As it is known the exact gluon propagator can be written in the form

$$G_G^{ab}\left(p\right)=\frac{\delta ^{ab}}{p^2}\left(\frac{g_{\mu _a\mu _b}\frac{p_{\mu _a}p_{\mu _b}}{p^2}}{1+\mathrm{\Pi }\left(p^2\right)}+\alpha \frac{p_{\mu _a}p_{\mu _b}}{p^2}\right),$$
(19)
where the function $`\mathrm{\Pi }`$ is defined at the one loop level, by the expression

$`\mathrm{\Pi }^{ab}\left(p^2\right)`$ $`=`$ $`\mathrm{\Pi }_G^{ab}+\mathrm{\Pi }_T^{ab}+\mathrm{\Pi }_{Gh}^{ab}+\mathrm{\Pi }_Q^{ab}`$
$`=`$ $`\delta ^{ab}\left(p_{\mu _a}p_{\mu _b}g_{\mu _a\mu _b}p^2\right)\mathrm{\Pi }\left(p^2\right),`$
which is diagrammatically represented in Figure 6.

Again, expanding in powers of the condensate parameter $`C`$, the total contribution can be decomposed in the standard one (associated to the ghost and the quark loop integrals) and the $`C`$ dependent ones $`\mathrm{\Pi }_G^{ab}`$ and $`\mathrm{\Pi }_T^{ab}`$ which can be readily calculated to be

$`\mathrm{\Pi }_G^{ab}`$ $`=`$ $`{\displaystyle \frac{g^2CN}{\left(2\pi \right)^4}}\delta ^{ab}\left(5g_{\mu _a\mu _b}+2{\displaystyle \frac{p_{\mu _a}p_{\mu _b}}{p^2}}\right),`$
$`\mathrm{\Pi }_T^{ab}`$ $`=`$ $`{\displaystyle \frac{3g^2CN}{\left(2\pi \right)^4}}\delta ^{ab}g_{\mu _a\mu _b}.`$
Summing up these particular two terms leads up to

$$\mathrm{\Pi }_C^{ab}\left(p^2\right)=\delta ^{ab}\left(p_{\mu _a}p_{\mu _b}g_{\mu _a\mu _b}p^2\right)\frac{2g^2CN}{\left(2\pi \right)^4p^2}.$$

Considering again our case the $`SU\left(3\right)`$ group, that is $`N=3`$, and disregarding in a first consideration the standard one loop contribution, the function $`\mathrm{\Pi }`$ reduces to

$$\mathrm{\Pi }\left(p^2\right)=\frac{6g^2C}{\left(2\pi \right)^4}\frac{1}{p^2},$$
(20)
which in turns leads out to the following condition for the poles of (19)

$$p^2m_G^2=0,$$

where

$$m_G^2=\frac{6g^2C}{\left(2\pi \right)^4}=0.25\left(GeV/c^2\right)^2.$$

Therefore, as the parameter $`C`$ has been defined before as real and positive PRD ; tesis , it follows that the transverse gluon mass correction becomes tachyonic. The ability of a tachyonic mass in producing improvements in models for the inter quark potential has been recently argued in the literature tachyon1 ; tachyon2 . The effect is linked with the introduction of a linearly rising term as a first correction to the Coulomb potential in the massless case tachyon2 . Therefore, the tachyonic result arising here for the gluon mass appears to be of interest for a derivation, within the present framework, of former successful phenomenological bound state models for mesons sommer1 ; sommer2 ; El-Hady . We are grateful to J.P. Vary and J. Qiu for the suggestion of this possibility. It should be also underlined that the tachyonic character of the gluon self-energy arising here, is fully compatible with the more general result of Fukuda for the normal Green functions in QCD, when gluon condensation is present fukuda . This circumstance then suggests, that the modified expansion under discussion could be a sort of perturbative representation of the Fukuda non-perturbative treatment of the gluon condensation in QCD. In a different manner, the approach of Celenza and Shakin, and other authors can be expected to be related with a different initial state also incorporating gluon condensation. It could be the point that, at the end, both pictures can appear to be equivalent at the level of exact results. This could be the case because the physical degrees of freedom are not the gluons and quarks ones, but their colorless composite combinations which physics is determined by the bound state dispersion relations.

Finally in this section it deserves to be stressed that if a modification of the ghost propagator is present, then a mass terms for the longitudinal propagator could also appear and the polarization operator (gluon self-energy) becomes non-transverse. This result means a violation of a Ward identity representing the gauge invariance by the presence of ghost (and fermion) condensates. Thus, it follows that the selection of the parameter of the initial state in PRD ; tesis  was the appropriate one for the satisfaction of a gauge invariance condition at this stage. The physical meaning of the alternative procedures of constructing the initial state needs however for clarification.

## VI Summary

The following results can be enumerated:

1) A calculation of the mean value of $`G^2`$ up to order $`g^2`$ in gluodynamics, was performed. The result coincided with a previous one obtained in a simpler approximation 1995 ; PRD  and which signals that quartic gluon self-interaction is relevant in fixing a non-vanishing gluon condensation parameter already in the first order of the perturbation theory Zakharov .

3) The one loop corrections for the quark masses under the influence of the condensate were evaluated. When only the terms reflecting the effect of the condensate are considered, the results for the constituent masses showed a branch which is real as a function of all the values of the quark current mass. Another solution which is real up to a critical value of the Lagrangian masses, and having a vanishing proper mass in the limit of small current masses is also obtained. The former branch predicts values of one third of the nucleon mass for the masses of the light quarks. The effective quark masses for all the quark flavors were evaluated as determined by their known current mass in the mentioned approximation.

4) After the inclusion of the rest of the terms completing the full one loop contribution to the quark self-energy, an interesting outcome arises. The $`u`$ and $`d`$ quarks, although showing high peak values of the propagator components near the former values of the squared proper masses, become confined. That is, their propagators do not show pole structures. On the other hand the $`s`$, $`c`$, $`b`$ and $`t`$ quarks acquire poles but now not at real values of $`p^2`$. Therefore, these free modes are damped. These properties suggest a possibility for the explanation of the disparity of stability properties in the nucleons and higher resonances. It would work as follows, the absence of propagation solutions for the $`u`$ and $`d`$ quarks could assure the unlimited life time of nucleons. Moreover, the damped oscillations present for the rest of the quarks can create decay modes for the bound states composed of them.

5) By mean of the obtained mass values, the ground state energies selected within each of the groups hadron resonance reported in Report , were estimated through the simple addition of the masses of their known quark constituents. It could be remarked that the obtained values are reasonable ones in comparison with the data, after taking into account the high errors in the experimentally determined Lagrangian masses. This fact supports the hypothesis about the main glue nature of the constituent mass values of the light quarks. The small contribution of the binding potential energy of such resonance to their total mass is also indicated and thus suggests the applicability of non relativistic approximations in their study.

6) The contribution of the condensate to the gluon self-energy of order $`g^2`$ was calculated. A tachyonic result for the gluon mass arises for the parameter values employed in the construction of the initial state defined in PRD ; tesis . This conclusion is in concordance with the analysis of Fukuda for gluon condensation fukuda . As it was remarked before, there are recent arguments claiming the possibility of a tachyonic gluon mass tachyon1 ; tachyon2 ; Hoyer3  and its role in improving the quark antiquark potential in bound state models of mesons tachyon2 .

It seems useful to comment here some basic issues related with the modified expansion being considered:

a) The Feynman expansion depending on the condensate parameter and taken for the gauge parameter $`\alpha =1`$ corresponds with the Wick expansion around a physical state of the interaction free QCD. In this somewhat limited way the arbitrariness of the gauge parameter in the result for the gluon mass in 1995  was solved in PRD ; tesis . However, in order to furnish a fully consistent picture including renormalization, a physically equivalent version of the expansion for arbitrary values of $`\alpha `$ should developed. This question will be addressed in future works.

b) In PRD ; tesis  the condensation parameter $`C`$ was defined as a real and positive one. This result determined the tachyonic character for the gluon mass found here. Further, the sign of $`C`$ also fixes the real values of the order of the one third of the nucleon mass for the constituent quark masses also calculated in the present work. The opposite sign produces confined quark modes.

c) An aspect which is important to underline is that the selection of the parameters defining the initial state in PRD ; tesis  was designed also in order to impose the absence of a modification for the free ghost propagator. However, as it followed from the present analysis, this property in turns, is related with the fact that a condensate modification of the ghost propagator can produce a longitudinal contribution to the self-energy. But, such a term should not exist for the transversality condition of the polarization operator (a Ward identity) to be obeyed. Thus, its appearance could break manifestly the gauge invariance. Therefore, the present work is also given foundation to the non-modified ghost propagator choice considered in 1995 ; PRD ; tesis .

At the present stage we identify some important questions to be addressed in the further extension of the work:

1) To explore the possibility of showing the exact gauge parameter independence of the physical results within the modified expansion.

2) To address the general proof of the renormalizability in the expansion.

3) To search for a derivation of existing successful bounded state models for the heavy quark mesons sommer1 ; sommer2 ; El-Hady , as realized by the a ladder approximation for the Bethe-Salpeter equation within the proposed modified expansion. The presence of a tachyonic gluon propagator in the approach (which is argued to have the effect of introducing a linearly rising component in the inter-quark potential tachyon2 ) and the obtained constituent values for the light quark masses already support such a possibility.

4) Improve the study of the effective action as a functional of the condensate parameter done in 1995  by also including the gauge field in order to search for a variant of the leading logarithm model useful for the investigation of field configurations associated to inter-quark strings, nucleon-nucleon potentials, etc.

5) It can be underlined that the considered expansion seems that could furnish a general way for the perturbative description of pair condensates in Many Body theory. An important example coming to the mind is the BCS state of the usual superconductivity. Pair states being more closely connected with QCD were discussed in Mishra ; Diakonov . In addition, in Diakonov  it was advanced before the idea of using diquark condensates as generating an alternative for the Higgs mechanism.

Finally, we would like to conjecture about a possibility suggested by the results of this work. It is related with the question about whether the mass spectrum of the whole three generations of fundamental fermions could be predicted by a slight generalization of the modified perturbation expansion under consideration. In this sense, the given arguments led us to the idea that (after the introduction of quark condensates along the same lines as it was done for gluon ones in 1995 ; PRD ; tesis ) the received perturbation expansion can have the chance to predict both, the Lagrangian mass and the constituent quark mass spectra of the three families of fundamental fermions. The fermion condensates as described in the proposed perturbative way, would have the role of producing the Lagrangian quark masses, through the chiral symmetry breaking. The gluonic condensates, in one hand, and as illustrated here, could be responsible of generating states of large constituent mass for the low mass quarks states. In another hand, it seems feasible that the higher order radiative corrections (including color interactions with the condensate as well as chiral symmetry corrections) could also determine the mass spectra for leptons and neutrinos. The smaller scale for the masses of these particles could be produced by the lack of lower order either color or electric-weak interaction terms in their self-energy. Therefore, the possibility that the Lagrangian mass spectrum of the three generations of the fundamental fermions could be predicted by a modified perturbation expansion of the sort being proposed is suggested. Work directed to investigate this possibility will be considered elsewhere. The modified expansion under study could be a useful technique of implementing this idea.

###### Acknowledgements.

The authors wish to deeply acknowledge the helpful advice and comments of Profess. J. Pestieau, J. Alfaro, M. Lowe, C. Friedli, M. Hirsch and A. González. One of the authors (A. C. M.) also would like to express his gratitude to the Abdus Salam ICTP, and in particular to its Associate-Membership Program for the whole general support. The support of the International Institute for Theoretical and Applied (IITAP) (UNESCO and Iowa State University) and the Christopher Reynolds Foundation during the completion of this work at IITAP is also greatly acknowledged. The suggestion of the possibility of applying the procedure for describing bound state model for mesons is acknowledged from Dr. J. P. Vary and Dr. J. Qiu. Finally, we are grateful to K. Tornés for the help in the preparation of the manuscript.