the precise shape of the quark - gluon vertex plays an essential role in the understanding of many important properties of the infrared behavior of qcd @xcite .
this vertex is in particular fundamental in the analysis of the origin of dynamical chiral symmetry breaking in the context of truncated dyson - schwinger ( ds ) equations .
indeed , the ds equation for the quark self - energy requires the knowledge of this vertex , together with the gluon two - point vertex function .
even if this kind of object depends on the gauge - fixing , one expects that the gauge - invariant physical content of the theory is encoded on its behavior .
for example , the meson spectrum can be studied by injecting such vertices and two - point functions in appropriate bethe - salpeter equations , see for example @xcite . for these reasons ,
the determination of the quark - gluon vertex has been a topic of intensive work in the past .
it has been first studied in the perturbative domain .
a general parametrization based on its symmetries and a one - loop calculation of the abelian part in feynman gauge was performed in @xcite .
some years later , the complete one loop calculation in an arbitrary linear covariant gauge and in arbitrary dimension was performed in @xcite ( where previous partial results are reviewed also ) .
more recently the complete 2-loop off - shell calculation has been performed in @xmath0 @xcite .
these perturbative studies lead to a satisfactory description of this vertex at high momenta . however , the infrared sector , which is of direct interest for the understanding of chiral symmetry breaking , is not accessible to perturbation theory . to cope with problem ,
mainly two methods have been used : lattice simulations and ds equations .    in both of these approaches ,
the most convenient gauge is the landau gauge . from the lattice simulation side ,
it is well - adapted because the gauge - fixing operation can be expressed in terms of an extremization , which can be implemented efficiently on the lattice . from the analytic side , this gauge has several advantages .
first , it is a linear gauge condition that manifestly preserves lorentz invariance .
second , among all linear covariant gauges , it is the only one that includes the symmetries of the curci - ferrari gauge @xcite .
this implies , in particular , the existence of some non - renormalization theorems , that simplify considerably the algebraic work @xcite .
third , in this gauge , the gluon propagator is transverse .
consequently , only the vertices where the gluon lines are multiplied by transverse projectors contribute to the correlation functions .
accordingly , all consideration below will be done in this gauge fixing .
there have been many analysis of the infrared behavior of the quark - gluon vertex based in the study of truncated ds equations ( for a review , see @xcite ; for more recent work on the subject , see @xcite ) . in parallel the intermediate and infrared regime of the vertex has been studied by lattice simulations .
many tensorial components have been analyzed @xcite for some particular configuration of momenta .
moreover , for the tensorial component that is already present at bare level , which is expected to be `` dominant '' , more general configurations of momenta have been considered @xcite .    in continuum analytical approaches , the gauge fixing is generally implemented through the faddeev - popov construction , which is a perfectly valid procedure in the ultraviolet regime , where standard perturbation theory works well and where the faddeev - popov construction is well justified . however , for momenta smaller or of the order of 1 gev ( which is the typical energy scale of the strong interaction ) this procedure is not fully justified .
the standard lore is that it becomes inappropriate because the coupling grows and perturbation theory eventually becomes a very bad approximation scheme .
this argument must be tempered because the real expansion parameter of qcd is of the order of @xmath1 and because lattice simulations , at least in some renormalization schemes , indicate that this parameter remains moderate ( see , for example @xcite ) .
accordingly it is reasonable to believe that some sort of perturbation theory should give results that are in qualitative agreement with experiments or lattice simulations in that regime .
there is , however , another reason why standard perturbation theory should become inappropriate in the infrared regime : the faddeev - popov procedure itself becomes questionable .
this originates in the existence of gribov copies @xcite in covariant gauge - fixings of non - abelian gauge theories that are ignored in the faddeev - popov construction .
these copies play no significant role in the ultraviolet regime but a systematic treatment of the associated ambiguity remains an open problem in the infrared .
the most developed procedure which aims at taking into account the effects of gribov copies is probably the gribov - zwanziger ( gz ) construction @xcite .
the purpose of this approach is to restrict the functional integral to the first gribov region where the faddeev - popov operator is positive definite .
gribov showed that in this region there are no copies that are an infinitesimal transformation of one another .
unfortunately it was shown later that there are _ finite _ gauge transformations that preserve the landau condition @xcite .
consequently the gz procedure does not remove all copies .
moreover , the implementation of the restriction of the functional integral to the first gribov region is not completely justified from first principles .    despite these difficulties , first studies of
the gz predicted that the gluon propagator would vanish at zero momentum , instead of diverging as in the bare theory .
a similar behavior was found also as solutions of truncated ds equations @xcite and in non - perturbative renormalization group ( nprg ) calculations @xcite ( for a general review on this formalism , see @xcite ) .
this kind of solution is called `` scaling solution '' .
the suppression of the gluon propagator was rapidly confirmed by lattice simulations . however , when lattice simulations became more precise they showed that the gluon propagator in landau gauge did not behave precisely as the `` scaling solution '' that ds had suggested , neither in quenched lattice simulations @xcite nor in the presence of dynamical quarks @xcite .
in fact , lattice simulations favored a gluon propagator which tends to a finite value at low momentum in @xmath0 . this massive behavior ( sometimes called `` decoupling '' solution )
gives a propagator which is suppressed with respect to the bare one but not as strongly as in the `` scaling solution '' of ds equations . in parallel such decoupling solutions
were also obtained from ds and nprg equations @xcite , as well as in a `` refined '' version of the gz construction @xcite .
nowadays , all approaches finally converge toward a massive - like gluon propagator in yang - mills theory .
moreover , in what concerns the ghost sector , there is also a consensus in favor of a massless behavior , with a slight enhancement at low energy but with the tree level power - law .
finally , we note that ds and nprg equations favor a coupling constant bounded in the infrared at not too large values , see for instance @xcite .
all these observations led two of us to test wether a simple perturbative analysis where a massive gluon propagator was introduced at the bare level @xcite could reproduce quantitatively the lattice data for several correlation functions .
this lagrangian is invariant under a modified becchi - rouet - stora - tuytin ( brst ) symmetry that ensures the renormalizability of the model @xcite .
however , because of the gluon mass , the associated transformation is not nilpotent . as a consequence
, the textbook construction of a physical space based on the co - homology of the brst charge does not apply .
it must be stressed that this difficulty is common to all methods that go beyond the standard perturbative fp method . in particular
, the standard definition of the physical space does not apply to the gz model but even stronger , does not apply to the procedure implemented on landau - gauge lattice simulations .
in fact the standard perturbative definition based on the cohomology of the brst operator includes in the physical space the transverse gluons .
but this is not satisfactory for two reasons .
first , on physical grounds , we know that gluons are confined and are not part of the physical spectrum .
second , if trasverse gluons were in the physical space , the spectral density associated with their two - point function should be positive and lattice simulations cleary show positivity violations @xcite .
all these considerations show that a proper definition of a physical space is an open and important problem . in this sense
, the lack of a nilpotent brst charge in the curci - ferrari should not prevent us from studying it in relation with the lattice simulations of correlation functions .
it is an open question to know if one can define a physical space for the curcci - ferrari model where all states would have a positive norm .    in the present article
we exploit this massive extension of the faddeev - popov lagrangian in order to study the quark - gluon vertex in landau gauge .
our aim is twofold .
first , we want to compare our findings with lattice simulation in order to have a further test of the ability of the curci - ferrari model to reproduce lattice simulations .
as we show below , we are able to compute the _ full _ vertex at one loop in any dimension , for all tensorial structures and for any configuration of momenta .
the comparison of this one loop result is in reasonable agreement with all available lattice data ( in many cases the agreement is very good ) . in particular , we discuss in detail the results of a certain scalar function ( @xmath2 in the terminology of @xcite ) where a significant discrepancy has been observed between lattice and ds results @xcite .
we discuss this discrepancy in detail giving a very simple explanation of its origin .
second ( and more importantly ) as explained before , this vertex is very important in various physical applications of ds and bethe - salpeter equations .
we expect that the present result will be useful in future physical applications in these equations .
we want to stress here that , as it stands , the simple one - loop calculation presented in this article does not encode a spontaneous breaking of the chiral symmetry .
however , the renormalization - group flow leads to a strong enhancement of the quark mass in the infrared , which mimicks , for all practical purposes , the behavior induced by this symmetry breaking .
this is a key ingredient for reproducing the infrared behavior of the quark - gluon vertex .
the rest of the article is organized as follows . in sec .
[ sec_model ] we describe in more details our model and shortly review the systematic comparison of one - loop computations with lattice data .
we give some details on the one - loop calculations in sec .
[ sec_one_loop ] and finally present our results in sec .
[ sec_results ] .
in this article , we pursue our systematic analysis of the correlation functions of qcd in the landau gauge , in a scheme where perturbative calculations can be performed in a controlled way for all energies . as explained in the introduction we consider the ( euclidean ) curci - ferrari lagrangian @xcite : @xmath3 where @xmath4 is the gauge field , @xmath5 and @xmath6 are ghost and antighost fields
, @xmath7 is a lagrange multiplier that enforces the landau gauge condition @xmath8 , @xmath9 is the coupling constant and @xmath10 it is important to recall that mass is introduced at the level of the gauge - fixed action . as a consequence , the bare propagators decrease at large momenta as @xmath11 and
the mass induces a soft breaking of the brst symmetry . as a consequence ,
all the ultraviolet properties of the standard faddeev - popov action are preserved to all orders in perturbation theory , in particular renormalizability and the behavior of correlation functions .
one of the main consequences of the mass term is that the low energy properties of the system are more regular than within the usual faddeev - popov approach . in particular
, it is possible to find renormalization schemes in which the coupling constant remains finite down to the deep infrared @xcite , at odds with the findings of standard perturbation theory which presents a divergence of the coupling constant at some energy of the order of 1  gev .
the absence of such a landau pole opens the way to a consistent treatment of the low energy regime of the theory within perturbation theory .
the general idea was tested in @xcite where the quenched gluon and ghost propagators were computed at one - loop order .
the comparison with lattice data is very satisfactory , with a discrepancy of at most 10% in the whole range of energy .
the analysis of the 3-point correlation functions ( three gluons and ghost - gluon ) was performed under the same line of investigation in @xcite .
the comparison with lattice data is also very good although the lattice data are more noisy in this case .    more recently ,
we have considered the influence of dynamical quarks on the theory , which are governed by the lagrangian : @xmath12 where @xmath13 are euclidean dirac matrices satisfying @xmath14 , the flavor index @xmath15 runs over the @xmath16 quark flavors . ] and the covariant derivative of the quark field reads @xmath17 in @xcite we have included the one - loop contribution of the quarks to the previously computed gluon propagator and have computed the quark propagator to the same order .
again , the gluon propagator compares very well with the lattice results .
functions , which gives an indirect effect on the ghost propagator .
again , the comparison of the ghost propagator with lattice data is very satisfactory . ]
the quark propagator can be decomposed in two independent structures in spinor indices , often called the mass ( proportional to the identity in spinor indices ) and field renormalization ( proportional to a dirac matrix ) .
the one - loop results for the quark mass is in good agreement with the lattice data . in particular
, it reproduces the strong enhancement at energies smaller than roughly 2 gev , which is reminiscent of the chiral symmetry breaking . below .
] however , we do not expect that the present simple one - loop calculation can break spontaneously the chiral symmetry . on the other hand , the field renormalization of the quark
is not reproduced correctly , even at a qualitative level .
this can be attributed to the fact that the one - loop contribution to this quantity is unusually small @xcite .
indeed , when the gluon mass is set to zero , the field renormalization has no contribution at one loop but it does have contributions at two - loops , that , in consequence , become dominant .
this model has also been used in another line of research in a recent series of articles : the extension of yang - mills theory at finite temperature was considered within the same framework with the goal of studying the deconfinement transition @xcite . it was found that a simple one - loop calculation can reproduce the known phenomenology [ second order transition for the su(2 ) group and first order transition for su(3 ) ] .
a next - to leading calculation in the su(2 ) case leads to a quantitative agreement for the critical temperature .
the determination of the quark - gluon vertex at one loop requires to compute the two diagrams depicted in fig .
[ fig_diags ] .
the result can be decomposed in twelve independent tensorial structures .
we follow here the convention of @xcite and express the quark - gluon vertex as @xmath18 with @xmath19 we give in table [ tab_struc ] the tensorial structures @xmath20 and @xmath21 , together with some informations concerning the scalar functions @xmath22 and @xmath23 .    [ cols="^,<,^,^",options="header " , ]     the calculations are performed by following the same procedure as described in @xcite .
the feynman rules are implemented in a mathematica procedure and the integrals over momentum written in terms of the passarino - veltman integrals @xcite : @xmath24[(p+q)^2+m_2 ^ 2]}\\    & c_0(p_1,p_2,m_1,m_2,m_3)=\\ & \nonumber\int \frac { d^dq}{(2\pi)^d}\frac 1{[q^2+m_1 ^ 2][(p_1+q)^2+m_2 ^ 2][(p_1+p_2+q)^2+m_3 ^ 2]}\end{aligned}\ ] ] our expressions for the different scalar functions expressed in terms of these integrals are available in the supplemental material @xcite .
in actual calculations , it is convenient to use the feynman trick to perform the integrals . in order to reduce the number of feynman parameters
, we use the identity : @xmath25 as many times as possible to reduce the number of propagators appearing in the diagrams . the price to pay is that there appears differences of masses in the denominator . for general momentum configurations , after performing the integrals over the internal momentum , we are left with a double integral over feynman parameters .
however , for integer space dimension and for some momentum configurations that are considered in next section , all integrals can be realized analytically .
we thus obtain analytic expressions for the one - loop vertex for those configurations .
we have made several checks of our one - loop expressions .
we have compared our results with those of @xcite which were obtained in the limit of vanishing gluon mass .
our check consists in comparing numerically our results for vanishing gluon mass with those of @xcite .
this was done for several momentum configurations and we find a perfect agreement .
therefore , in the regime where the external momenta are large compared with the gluon mass , our results reproduce those of standard perturbation theory ( including anomalous dimensions and all angular dependences ) . as a consequence , we do not study this regime further below .
another check comes from the behavior of the coupling constants in the chiral limit ( limit of vanishing quark mass ) . in this limit
, we easily derive the ward identity associated with the chiral symmetry : @xmath26 where @xmath27 . by using the anticommutation relations of dirac matrices ,
it is easy to check that , in the chiral limit , the structures @xmath28 , @xmath29 , @xmath30 , @xmath31 , @xmath32 and @xmath33 are not compatible with this ward identity .
we therefore conclude that the associated coupling constants must vanish in the chiral limit .
we have explicitly verified that our one - loop expressions are compatible with this result .
we regularize the theory by using the infrared - safe ( is ) scheme @xcite . as usual , we introduce renormalized quantities , which are related to the bare ones by renormalization factors : @xmath34 the different renormalization factors are defined by imposing : @xmath35 the last two conditions are consequences of non - renormalization theorems which relate the divergent parts of the different renormalization factors @xcite . here , these relations are imposed also for the finite parts .
the last renormalization condition corresponds to the taylor scheme and is obtained by defining the coupling constant through the ghost - gluon vertex function with vanishing momentum for the external ghost .
the condition ( [ nonrenmcond ] ) can not be expressed directly in terms of vertex ( or correlation ) functions . in particular , the mass parameter thus defined is not directly related to the value of the gluon propagator at vanishing momentum .
expressions for the renormalization factors can be found in @xcite .    in the following , we use the renormalization - group improvement for the quark - gluon vertex .
it is obtained by integrating the rg equation for the quark - gluon vertex function @xmath36\gamma_\mu=0 , \end{split}\ ] ] where @xmath37 by integrating eq .
( [ eq_rg_eq ] ) we relate the quark - gluon vertex at different scales through the equation : @xmath38 here @xmath39 , @xmath40 and @xmath41 are obtained by integrating the beta functions with initial conditions given at some scale @xmath42 and : @xmath43 observe that the running quark mass appears explicitely in the previous expressions .
this leads to a back - reaction of this running mass on the vertex functions . as a consequence
, the mass takes very different values at @xmath44 and @xmath45  gev . even at a perturbative level , this is a large effect and , as discussed in the next section , is at the origin of the enhancement of some scalar components of the vertex .    in
the is scheme , for the initial conditions of interest for comparing with the lattice data , the coupling constant does not present a landau pole . in order to avoid large logarithms
, we evaluate the right - hand - side of eq .
( [ eq_int_rg ] ) with @xmath46 of the order of the external momenta when they are larger than the mass . in practice
, we used @xmath47 .
it is important to stress that contrarily to what happens in the pure glue sector , there is an significative dependence on the scheme in the present calculation .
in fact , it is important to chose the renormalization point at a scale that remain fixed when all momenta go to zero as done in @xcite .     in the following ,
we compare our results with lattice data which were obtained in the quenched approximation .
consequently , the beta functions for @xmath9 and @xmath48 as well as the anomalous dimensions for @xmath49 and @xmath50 are evaluated at @xmath51 .
therefore , the quark mass does not back - react on these quantities , but it does back - react on the vertex .
the results presented in the previous section were obtained for an arbitrary momentum configuration and for all tensor structures .
there have been lattice studies of the quark gluon vertex for general kinematics @xcite for a particular tensorial structure in the quenched approximation .
our results are in general qualitative agreement with those of @xcite at the level of precision of the data .    for some particular momentum configurations ,
lattice data were obtained with larger statistics and therefore smaller error bars . in the following ,
we concentrate on these configurations which give more stringent tests of the quality of our results .
we perform a comparison with lattice data in three momentum configurations .
the first one corresponds to configurations where the gluon has zero momentum with arbitrary momentum of the quark . in this configuration
many tensorial structures were obtained in the quenched approximation @xcite and we compare below our results for all available functions .
the second configuration corresponds to equal momenta for the quark and antiquark ( and , accordingly , gluon momentum equal to minus the double of them ) .
again , many tensorial structures were obtained in the quenched approximation for this configuration @xcite and we study all of them below . on top of these two families of configurations ,
we consider a more general configuration with equal momentum in modulus of quarks and anti - quarks but for various possible angles between them ( or , equivalently , different values of the modulus of the gluon momentum ) .
this configuration includes as particular cases the previous ones but for general angles only a single tensorial structure has been simulated @xcite ( also in the quenched approximation ) . for the comparison with the simulations ,
we focus on two values of modulus of the gluon momenta and also for the completely symmetric configuration where the incoming momentum of the gluon is equal in modulus to the quark and antiquark momenta .
it must be pointed out that the employed data are more that ten years old .
more recent data with improved lattice parameter could allow to perform a more precise test of our curves but they do not exist for the moment .    in the end of this section
we finally discuss the effect of including dynamical quarks for the present vertex .
even if there are no simulations of this vertex in the unquenched case , there are some simulations of 2-point functions @xcite that allows to fix the parameters and , in consequence , we are able to make a prediction of the behaviour of the quark - gluon vertex for the conditions that are present in those simulations .
these could be compared to future unquenched simulations of the quark - gluon vertex .
let us first consider the case of zero gluon momentum . in this case
, the tensor structure simplifies to the following form @xmath52.\ ] ] moreover , all integrals appearing in the vertex functions can be performed analytically . by introducing @xmath53 @xmath54,\ ] ]
@xmath55\]]and @xmath56 we obtain in @xmath57    @xmath58\nonumber\\ & + \frac{g^2}{192 \pi ^2 m^4
n p^4 } \log\left(\frac{m}{m}\right )     \left\{-6 m^8 \left(n^2 - 1\right)+m^6 \left[m^2 \left(10 n^2 - 9\right)-\left(14 n^2 + 3\right ) p^2\right]\right.\nonumber\\ & \left.-3 m^4     \left[m^4 n^2 - 14 m^2 p^2 + 3 \left(n^2 + 2\right ) p^4\right]+3 m^2 \left(m^6 - 7 m^4 p^2 + 15 m^2     p^4-p^6\right)\right.\nonumber\\ & \left .- n^2 \left(m^2+p^2\right)^4\right\}-\frac{g^2}{192 \pi ^2 m^4 n p^4}\log \left(\frac{m^2+p^2}{m^2}\right)\left(m^2+p^2\right)^3   \left[n^2 \left(m^2+p^2\right)-3     m^2\right]\nonumber\\ & + \frac{g^2}{384 \pi ^2 m^4 n p^4 y_1}\left\ { m^2 p^2 y_1 \left[12 m^4      \left(n^2 - 1\right)+m^2 m^2 \left(6 - 8 n^2\right)+m^2 p^2\left(19 n^2 + 6\right )     \right.\right.\nonumber\\ & \left.-6
\mathbb{l}_2 \left(2 m^2 - 5 m^2+p^2\right )     \left(m^2-m^2+p^2\right)-2 n^2 \left(m^2+p^2\right)^2\right]\nonumber\\ & + \mathbb{l}_1 \left[m^4 + 2 m^2     \left(p^2-m^2\right)+\left(m^2+p^2\right)^2\right ] \left[6 m^6     \left(n^2 - 1\right)+ m^4 m^2     \left(3 - 4 n^2\right)\right.\nonumber\\ & + m^4\left.\left(8 n^2 - 3\right ) p^2\right)\left.\left .-
m^4m^2     \left(n^2 - 3\right)+6m^2 m^2 p^2+m^2\left(n^2 + 3\right ) p^4-n^2     \left(m^2+p^2\right)^3 \right]\right\},\end{aligned}\ ] ]    @xmath59\nonumber\\ & -\frac{g^2}{384 \pi ^2 m^4 n p^6}\log\left(\frac{m}{m}\right)\left\{12 m^8 \left(n^2 - 1\right)-m^6 \left(10 n^2 - 9\right ) \left(2 m^2-p^2\right)+6 m^4 \left(m^4     n^2+p^4\right)\right.\nonumber\\ & + m^2 \left[-6 m^6 - 9 m^4 p^2 + 6 m^2 \left(5 n^2 - 1\right ) p^4 + 9 p^6\right]\nonumber\\ & \left.+2 n^2     \left(m^8 + 4 m^6 p^2 - 9 m^4 p^4 + 19 m^2 p^6+p^8\right)\right\}\nonumber\\ & -\frac{g^2}{384 \pi ^2 m^4 n p^6}\log \left(\frac{p^2}{m^2}+1\right)\left(m^2+p^2\right)^2   \left[m^2 \left(3 p^2 - 6 m^2\right)+2 n^2     \left(m^2+p^2\right)^2\right]\nonumber\\ & + \frac{g^2}{768 \pi ^2 m^4 n p^8}\frac{\mathbb{l}_1}{y_1}\left\{12 m^{10 } \left(n^2 - 1\right ) p^2+m^8 p^2 \left[m^2 \left(30 - 32 n^2\right)+\left(22 n^2 - 21\right )     p^2\right]\right.\nonumber\\ & + m^6 p^2 \left[2 m^4 \left(13 n^2 - 9\right)+3 m^2 \left(1 - 2 n^2\right ) p^2+\left(10     n^2 - 9\right ) p^4\right]\nonumber\\ & -m^4 \left[150 m^8 n^2 + 6 m^6 \left(n^2 + 1\right ) p^2+m^4 \left(10
n^2 - 3\right ) p^4 + 2 m^2 \left(2 n^2 - 3\right ) p^6 + 3 p^8\right]\nonumber\\ & + m^2 \left[300 m^{10 } n^2 + 2 m^8     \left(3 - 149 n^2\right ) p^2+m^6 \left(4 n^2 + 15\right ) p^4 + 9 m^4 p^6-m^2 \left(4 n^2 + 3\right )     p^8\right.\nonumber\\ & \left.\left.-\left(2 n^2 + 3\right ) p^{10}\right]+150 m^8 n^2 y_1 ^ 2 - 2 n^2 \left(m^2+p^2\right)^2 \left(75     m^8+m^6 p^2 + 3 m^4 p^4 + 3 m^2 p^6+p^8\right)\right\}\nonumber\\ & -\frac{g^2}{128 \pi ^2 m^4 n p^2}\mathbb{l}_2 \left(m^2 + 5 m^2 n^2\right ) \left(m^2-m^2+p^2\right ) ,   \end{aligned}\ ] ]    and @xmath60 - 6 m^4     \left(n^2 - 1\right ) p^2\nonumber\\ & -\frac{3 \mathbb{l}_1}{y_1 } m^4 \left(n^2 - 1\right ) \left[p^2
\left(m^2+m^2\right)+\left(m^2-m^2\right)^2\right]\nonumber\\ & -p^2 \mathbb{l}_2   \left(m^2-m^2+p^2\right ) \left[m^2+n^2 \left(p^2 - 3 m^2\right)\right]\left.\right\}.   \end{aligned}\ ] ]    when performing lattice simulations , these functions are obtained from the simulated vertex @xmath61 by projecting in the various tensor structures in the following way @xcite : @xmath62 @xmath63+\lambda_1(p^2)\right)\ ] ] and @xmath64\ ] ] where , in these expressions , no implicit sum over repeated indices is meant .
it is interesting to note that , in order to extract the function @xmath65 , a contribution of @xmath66 must be added .
this issue will be discussed in more detail below .    to compare our results with the lattice simulation of @xcite we renormalize the model in the is scheme and implement the corresponding renormalization group as discussed in sect .
[ sec_one_loop ] .
it is important to stress that the simulations of @xcite have been done in the quenched approximation .
accordingly , we use the beta - functions of the pure yang - mills theory derived in @xcite for the gluon mass and coupling constant .
this implies that the coupling constant and the gluon mass can be fixed once the gluon and ghost 2-point functions have been fitted .
the only remaining free parameter is the valence quark mass ( and the overall normalization factor ) . to determine the best fit parameters for the ghost and gluon propagators , we consider the errors functions : @xmath67 we represent these functions in fig .
[ fig_error_aacc ] .     for the gluon propagator ( upper pannel , for 19% , 22% and 24% ) and @xmath68 for the ghost propagator ( lower pannel , for 4% , 7% and 10%).,title="fig : " ]   for the gluon propagator ( upper pannel , for 19% , 22% and 24% ) and @xmath68 for the ghost propagator ( lower pannel , for 4% , 7% and 10%).,title="fig : " ]    the best fits are obtained for @xmath69 and mass @xmath70  gev at the scale 1 gev which compare well with those found in @xcite .
as in that reference , this value has been compared to other values of the coupling by taking into account the rg running , finding a good agreement with best estimates up to a 10 % error .
the valence quark mass has been fixed in order to fit properly the quenched quark propagator .
more precisely we choose the quark mass in order for the constituent quark mass to match the lattice value .
doing so , one obtains a mass function ( the scalar part of the quark propagator , see @xcite ) that reasonably agree with the lattice data from @xcite as can be seen in fig  [ fig_mass ] .
one see that this curve mimick the enhancement of the mass in the infrared but not the spontaneous chiral symmetry breaking . indeed , the constituent quark mass vanishes when the current mass at 1 gev tends to zero .
the corresponding value of the quark mass is @xmath71  gev at the scale @xmath72 , and it gives a constituent quark mass @xmath73  gev .
these values tend to overestimate the mass function @xmath74 in the uv by approximatively 30% ( for example , at 3 gev we obtain 0.13 gev instead of 0.10 gev ) and to underestimate the propagator in the ir ( below 1 gev ) by a maximum value of approximatively 18% .
having said that , the overall behavior is approximatively correct . it must be pointed out when fixing the parameters in order to compare both to the quark mass function and to the vertex functions , a practical difficulty takes place .
the simulations have been done for both quantities with different mass parameters . for the quark mass function
many bare masses were employed . in the vertex case ,
the employed values for the bare valence quark mass were @xmath75  mev and @xmath76  mev @xcite .
however the authors claim that the vertex functions are almost insensitive to the bare quark mass and employ the data corresponding to @xmath75  mev . in order to have a single set of data we choose the bare valence quark nearest to this value that has been used for propagators which is @xmath77  mev
this introduces a systematic error but that do not seem to be large because for relatively small masses like those , the mass function do not seem to depend on @xmath78 too much ( see @xcite and fig .
[ fig_mass ] ) .
one can also analyze the dependence on the valence quark mass of the vertex by studying the error level for the function @xmath66 defined as @xmath79 ^ 2}\left[\lambda_{1,\rm lt.}(p_i^2)-\lambda_{1,\rm      th.}(p_i^2)\right]^2,\ ] ] we also consider the analogous error level for the function @xmath80 .
for this last function the total error may be misleading as discussed below because it average some momenta where the the curve is essentially within error bars with the lowest momenta where errors are much larger .
is problematic .
consequently , we do not use this function for fitting the parameters . ]
they are both shown in fig .
[ errorl1 ] .
contour levels for the quantity @xmath81 ( upper panel ) and @xmath82 ( lower panel ) in the plane of gluon mass vs. quark mass .
the contour lines correspond to @xmath83 , @xmath84 and @xmath85 for @xmath81 and @xmath86 , @xmath87 and @xmath88 for @xmath82.,title="fig : " ]   contour levels for the quantity @xmath81 ( upper panel ) and @xmath82 ( lower panel ) in the plane of gluon mass vs. quark mass .
the contour lines correspond to @xmath83 , @xmath84 and @xmath85 for @xmath81 and @xmath86 , @xmath87 and @xmath88 for @xmath82.,title="fig : " ]     we observe that , even if the optimum value of the masses do not coincide when fitting @xmath66 and @xmath80 , there is a region of parameters where both errors are reasonably small ( that is less than @xmath86 ) .
moreover , the value of the gluon mass extracted from the ghost and gluon propagators ( @xmath70  gev ) is within an error region with less than @xmath87 .
it is important to stress that once the parameters @xmath9 , @xmath48 and @xmath89 have been fixed , they are used for the calculation of the quark - gluon vertex for all momentum configurations .
accordingly , they become pure predictions without free parameters of the present calculation .    , @xmath70  gev , and @xmath71  gev .
lattice data @xcite are marqued with dots and parametrized by their bare masses : @xmath90  mev ( squares ) , @xmath91  mev ( triangles ) , @xmath77  mev ( circles ) , @xmath92  mev ( diamond ) . ]    the resulting curves for the functions @xmath66 , @xmath65 and @xmath80 can bee seen in fig .
[ funesck0 ] .
scalar functions of the quark - gluon vertex : @xmath93 ( top ) , @xmath2 ( middle ) and @xmath94 ( bottom ) for a vanishing gluon momentum , as a function of the quark momentum .
the line are the one - loop results with @xmath71  gev .
the dots correspond to the lattice data of @xcite.,title="fig : " ]   scalar functions of the quark - gluon vertex : @xmath93 ( top ) , @xmath2 ( middle ) and @xmath94 ( bottom ) for a vanishing gluon momentum , as a function of the quark momentum .
the line are the one - loop results with @xmath71  gev .
the dots correspond to the lattice data of @xcite.,title="fig : " ]   scalar functions of the quark - gluon vertex : @xmath93 ( top ) , @xmath2 ( middle ) and @xmath94 ( bottom ) for a vanishing gluon momentum , as a function of the quark momentum .
the line are the one - loop results with @xmath71  gev .
the dots correspond to the lattice data of @xcite.,title="fig : " ]    we find a good agreement for the function @xmath66 .
all the point are below a 16% of error except the first point ( @xmath95  mev ) where a larger error is observed ( 27% ) .
let us point out , however that our curve is systematically below the lattice one .
the behavior for @xmath80 is qualitatively reproduced , with a strong increase ( in absolute value ) at momenta of the order of 1  gev this last curve depends more strongly on the quark mass .
this can be understood by the observation that @xmath94 tends to zero in the chiral limit , see table  [ tab_struc ] .
the function is very well reproduced ( essentially within error bars ) except for the two lowest momenta ( @xmath96  140 and 400 mev ) where the errors are larger ( 37% and 33% respectively when comparing to the lattice central value ) .
let us , however , point out that the statistical error bars of these lattice points are large ( 22% in both points ) .
we therefore find that for a unique choice of parameters @xmath9 , @xmath48 and @xmath89 , we reproduce very well the ghost and gluon propagators and reasonably well the mass function @xmath74 and the vertex functions @xmath93 and @xmath94 , which is a highly nontrivial result .
the situation looks qualitatively different for the function @xmath65 , which tends to a constant at zero momentum in our approach while it diverges in the lattice data .
this disagreement was already observed in previous analytical treatments of this function @xcite and was interpreted as a consequence of the non - perturbative features of the theory .
we want to propose a completely different explanation for this phenomenon .
let us go back to eq .
( [ eql2 ] ) which is used to extract @xmath2 from the lattice data .
we observe that we must subtract from the function @xmath66 a quantity @xmath97.\ ] ] which is directly extracted from lattice data .
the difference is then divided by @xmath98 . in our calculation , the fact that @xmath2 tends to a constant originates from a compensation of these two terms in the limit @xmath99 , the difference being of order @xmath98 .
we can interpret this constant as a compensation of large numbers , a situation which is difficult to treat numerically . indeed , a small error in one of these two terms would lead to a divergence of @xmath65 when @xmath99 but this would be certainly an artifact .
note that such a divergence in @xmath65 would imply a nonanalytic behavior of the vertex function at small momenta [ see eq .
( [ zerogluonmomentum ] ) ] since the @xmath99 limit would depend on the direction ( in momentum space ) which is used to perform this limit . in our calculations , which are done in the framework of the action ( [ eq_cf])([eq_quarks ] ) , the only non - analytic behaviors observed so far are a consequence of ghost loops ( see , for example , @xcite ) .
such contributions are not present at one - loop order in the present vertex and most probably their contributions of higher order should be extremely small .    if the disagreement takes it origin in partial compensations in the extraction of the @xmath65 function from the formula ( [ eql2 ] ) , the comparison should be better for the quantity that is directly extracted from lattice data ( that is @xmath100 ) than that for @xmath2 itself . in fig .
[ funesck0l1l2comb ] we compare our prediction for @xmath100 to the corresponding lattice values .
we observe a very good agreement in both expressions : our curve is almost within the lattice error bars except for the lowest momentum where the error is around the 17% with respect to the central value of simulations .
we must , however point out that the curve is systematically above the central values of lattice data .
the agreement between the @xmath100 with lattice data suggests that the true origin of the disagreement between analytical and lattice data for @xmath2 lies in the difficulty of obtaining a function from differences of large numbers .
let us point out that the origin of the error is certainly not statistical .
if it were so , the signal obtained from the lattice data for @xmath65 would be extremely noisy .
most probably the origin of the error is more systematical .
now it is important to recall that the lattice extraction of these scalar function requires a careful treatment of various lattice artifacts ( see @xcite ) .
a mismatch in the treatment of lattice artifact of the longitudinal and transverse parts of the vertex could easily introduce a systematic error that , when divided by @xmath98 , gives a large spurious contribution in the estimate of @xmath65 . in any case ,
[ funesck0l1l2comb ] shows clearly that not only the component transverse to the external momentum but also the longitudinal one are correctly reproduced by the present one - loop calculation .
comparison of the quantities @xmath93 and @xmath101 as a function of the momentum [ see eq .
[ eq_def_lambda2 t ] ] .
our one - loop results for @xmath93 ( full line ) are almost degenerate with those for @xmath102 ( dashed line ) .
the higher ( blue ) points correspond to @xmath93 and the lower ( grey ) ones to @xmath103 , extracted from @xcite .
in the latter case , our estimation of the error bar is probably pessimistic . ]      in this section we present the results for the kinematic configuration given by equal momenta of quark and anti - quark ( @xmath104 ) and , accordingly , @xmath105 . in this kinematics
the quark - gluon vertex simplifies considerably @xmath106\end{aligned}\ ] ] and the transverse projected vertex ( that is @xmath107 ) < has the form @xmath108\ ] ] where we have introduced @xmath109 .    in spite of this simplification , for this configuration of momenta
, integrals can not be done analytically and some remaining feynman parameter integrals must be performed numerically .
fortunately , having fixed the parameters from propagators , the calculation must be done only once for a given choice of parameters .    in this kinematic configuration ,
the scalar functions @xmath110 and @xmath111 are extracted from the vertex as @xmath112 and @xmath113 in this case , no subtraction is needed to extract @xmath110 and @xmath111 from the lattice data .    in fig .
[ funessim ] we compare our one - loop results with the lattice data using the same initial condition of the coupling constant and masses as in the previous case , @xmath69 , @xmath70 gev and for the value of @xmath71 gev at @xmath114 gev .
we observe that the function @xmath115 is well reproduced by our results ( mostly within the lattice error bars , except for the lowest momentum point where the error is of the order of 25@xmath116 ) .
the other function @xmath111 is seen to be very sensitive to the choice of the quark mass at 1 gev .
again , this is a consequence of the fact that , in the chiral limit , this function tends to zero , see table  [ tab_struc ] .
we observe that the qualitative behavior for this function is reproduced but the agreement is not quantitative .
in particular , our one - loop results do not reach such large values in the deep infrared where for the lowest value of momenta the error reach 45% .
moreover , the curve is systematically above the lattice curve .
this is similar to what has already been observed in the ds results of @xcite .
scalar function of the quark - gluon vertex : @xmath115 ( top ) and @xmath111 ( bottom ) for equal momentum ( p ) for the quark and antiquark .
the full line are the one - loop results and the dots correspond to the lattice data of  @xcite.,title="fig : " ]   scalar function of the quark - gluon vertex : @xmath115 ( top ) and @xmath111 ( bottom ) for equal momentum ( p ) for the quark and antiquark .
the full line are the one - loop results and the dots correspond to the lattice data of  @xcite.,title="fig : " ]      in this subsection we consider the case where the quark and antiquark have momenta with the same modulus but can have arbitrary angle between them ( or , equivalently , that one can vary the value of the gluon momentum ) .
it is clear that the two previous kinematic configurations correspond to particular cases of the more general situation considered here .
however , more general configurations of momenta were studied in the lattice @xcite , but only for a particular tensorial structure .
more precisely , in terms of the tensorial decomposition ( [ eq : decomp_2 ] ) , the function that is extracted in the lattice is @xmath109 as in the previous case .
we limit to three particular kinematic configurations . in the first two ,
the modulus of the gluon momentum taking values @xmath117 gev and @xmath118 gev .
the third corresponds to the completely symmetric case where the modulus of the three momenta are equal @xmath119 .
as before , initial conditions of the renormalization group equations correspond again to @xmath69 , @xmath70 gev and @xmath71 gev at @xmath114 gev .
the results are presented in figs .
[ fun_fully_sym ] and [ fun_fully_sym_2 ] . for @xmath117
gev the maximal error is of 12@xmath116 . for @xmath118  gev ,
our results are within lattice error bars with the exception of the point of lowest momentum ( @xmath120 gev ) where the error is of 18@xmath116 . in the completely symmetric configuration
the curve is almost within error bars except for the lowest momentum where the error is 25% . at the end of this section
we discuss on possible origins of these discrepancies .
scalar function of the quark - gluon vertex @xmath109 for @xmath121 and @xmath122 gev ( top ) and @xmath123 gev ( bottom ) .
the full line are the one loop results and the dots the lattice data of @xcite.,title="fig : " ]   scalar function of the quark - gluon vertex @xmath109 for @xmath121 and @xmath122 gev ( top ) and @xmath123 gev ( bottom ) .
the full line are the one loop results and the dots the lattice data of @xcite.,title="fig : " ]     scalar function of the quark - gluon vertex @xmath109 for the completely symmetric configuration @xmath119 .
the full line are the one loop results and the dots the lattice data of @xcite interpolated to reconstruct the symmetric configuration . ]
we conclude our study by introducing the effect of dynamic quarks in our calculations .
this influences the renormalization factor for the gluon @xmath124 and for the gluon mass @xmath125 , which both have a contribution from a fermion closed loop . on the contrary ,
the first contribution of dynamic quarks to @xmath126 and to the quark - gluon vertex comes at two loops .
we conclude that the effect of dynamic quarks only influence indirectly our calculation , through the modification of the @xmath127 functions for the coupling constant and the gluon mass and for the field renormalization @xmath128 .
the unquenched propagators for the quarks , gluons and ghosts where studied in @xcite , in the framework of the curci - ferrari action .
we use here the parameters that were obtained in that article , as best fits when comparing with the lattice data in the only case where the mass quark can be fixed from the a quark mass function , that is @xmath129 .
the best fit parameters in this case is @xmath130 , @xmath131  gev and @xmath132  gev at @xmath114  gev .
the corresponding values of the couplings were compared in @xcite with standard estimates of the coupling giving an agreement with an error of the order of 17 % , which is similar to the typical overall error in the infrared of the one - loop approximation .
using these values of the parameter we can now predict the vertex in the unquenched case .
it is interesting to note that in the present procedure this is an extremely simple calculation once the beta functions have been calculated and the parameters fixed by comparing to the 2-point functions .
this is in contrast to lattice simulations that are extremely costly in the unquenched case .
this is probably the reason why there are , to our knowledge , no lattice data for this vertex .
we focus here on the zero gluon momentum but other tensorial structures and kinematic configuration could be presented as well .
we present in fig .
[ unquenchedcurves ] the three scalar functions @xmath93 , @xmath103 and @xmath94 for different number of flavors .
the presence of dynamic quarks clearly tends to increase these functions at low energy .
scalar functions of the quark - gluon vertex in the unquenched case .
@xmath93 ( top ) , @xmath133 ( middle ) and @xmath94 ( bottom ) for a vanishing gluon momentum , as a function of the quark momentum .
the curves corresponds to for @xmath51 ( full line ) and @xmath16=2 + 1 ( small dashes).,title="fig : " ]   scalar functions of the quark - gluon vertex in the unquenched case .
@xmath93 ( top ) , @xmath133 ( middle ) and @xmath94 ( bottom ) for a vanishing gluon momentum , as a function of the quark momentum .
the curves corresponds to for @xmath51 ( full line ) and @xmath16=2 + 1 ( small dashes).,title="fig : " ]   scalar functions of the quark - gluon vertex in the unquenched case .
@xmath93 ( top ) , @xmath133 ( middle ) and @xmath94 ( bottom ) for a vanishing gluon momentum , as a function of the quark momentum .
the curves corresponds to for @xmath51 ( full line ) and @xmath16=2 + 1 ( small dashes).,title="fig : " ]      the overall agreement of the present one - loop calculation is reasonably good .
most of the lattice data are well reproduced with a good accuracy .
it is important to stress that most of the calculations are pure predictions without free parameters ( that are already adjusted via 2-point functions ) and that they are all obtained by one loop calculations .
however , some discrepancies are present and we analyze them in the present subsection . putting aside the @xmath134 , the largest discrepancy concerns the @xmath135 function in the equal momenta configuration .
even if in this case the order of magnitude of the function and general behavior is well reproduced , there is not quantitative agreement .
finally , one observes some quantitative disagreement for the function @xmath136 function in the symmetric configuration at small momenta .
we can say that there is a good general agreement with lattice data but the quantitative agreement is not as good as for gluon and ghost vertices .
a similar behavior is observed in the study of propagators @xcite .
we give now two possible ( probably complementary ) explanations for this .
first , one of the reasons of the success of a perturbative analysis in the gluon and ghost sector is that the relevant couplings ( the value of the 3-gluon and ghost - gluon couplings ) are moderate except , for the 3-gluon vertex for very small values of the coupling where the presence of the gluon mass , has already suppressed fluctuations .
however , in the quark sector , even if the coupling does not presents a landau pole , it reaches significantly larger values than in the pure - glue sector .
one measure of this is the value of the coupling @xmath137 or @xmath136 .
one can define a natural running coupling function for this functions as @xmath138 .
it is important to stress that even if for @xmath139 , all the definitions of the running coupling give a universal beta function , this is not the case for momenta @xmath140 . in consequence
, the moderate values of the coupling in a certain scheme does not exclude that it may be large in some other scheme .
in particular , one observes that for small values of @xmath141 the quark - gluon coupling is larger by a factor of 2.4 with respect to the ghost - gluon coupling ( giving a larger expansion parameter by a factor @xmath142 . in a sense , this is to be expected because it is well established ( see , for example , @xcite ) that the quark - gluon coupling that is required for generating dynamical chiral symmetry breaking ( dscb ) is larger that the one that is extracted from the ghost - gluon vertex .
now , the expansion parameter of the perturbative expansion in the present model in the ghost - gluon sector ( precisely discussed in @xcite ) may reach values of the order of 0.2 . in consequence ,
the expansion parameter in the quark - gluon sector seem to put the perturbative expansion in trouble at very low momenta .
this may indicate that even in the present model where the ghost - gluon sector seem to be treatable perturbatively , some sort of partial resummation is needed in the quark sector ( or , at least , in part of it ) .
these facts also reflects on a certain tension when the parameters have to be chosen in order to fit all the data at the same time . even if the overall agreement is good with a single set of parameters , a precise renormalization scheme had to be chosen because at very low momenta an important sensitivity is observed .
this last point is related to a second possible origin of discrepancies . in the quark - gluon vertex ,
the quantities that are not sensitive to the dcsb are much better es- timated that those that are sensitive to it .
the function @xmath137 , gives better agreement than the @xmath143 or the @xmath135 ( which are sensitive to the dsbc ) .
this manifests also in the fact that those quantities seem to have a larger sensitivity on the renormalization conditions .
now , it is not clear in what measure we have properly taken into account the role of the dsbc . even if rg effect gives an important enhancement of the quark mass in the infrared that quantitatively reasonably reproduces the mass function measured in the lattice , it is unclear at the moment if a proper treatment dsbc will require to use more sophisticated techniques as the introduction of composite fields for the chiral condensate .
this may be at the origin of the lower quality of quantities that are sensitive to the dsbc .
in this article , we have presented the calculation of the quark - gluon vertex at one loop in the massive ( curci - ferrari ) extension of the landau gauge .
the results were decomposed on the twelve independent tensorial structures .
we implement the renormalization - group improvement so as to avoid large logarithms .
we find that these analytic results compare qualitatively well with the quenched lattice data . in one case ( @xmath2 for a vanishing gluon mass )
, the difference between lattice data ( that predicts a divergence at low momentum ) and our ( which are finite at small momenta ) may be attributed to a difference of large numbers when extracting this function from lattice data .    at the quantitative level , we observe that the quantities that are not forbidden by chiral symmetry are reproduced at the level of 10% except for momenta of the order of 100 - 200 mev where larger errors are found , as detailed in the results section .
the others are much more difficult to reproduce and show strong dependence on the renormalization - group scheme .
this is probably to be attributed to the rather crude treatment of chiral symmetry breaking in our calculation .
moreover , in order to fit simultaneously both the mass function coming from the quark propagator and the vertex function , a certain tension in the choose of parameters is observed and in order for this to be fulfilled it is necessary to choose a precise renormalization scheme .
it is interesting to point out that even if , at the end , the agreement of the present one - loop calculation is reasonably good , the precision is not as good as in the ghost - gluon sector .
this may be related to a larger value of the associated coupling and some sort of resummation in this sector could be useful .
this issues will be addressed in the future .
m.  s.  bhagwat , a.  holl , a.  krassnigg , c.  d.  roberts and p.  c.  tandy , phys .
c * 70 * , 035205 ( 2004 ) .
r.  alkofer , c.  s.  fischer , f.  j.  llanes - estrada and k.  schwenzer , annals phys .
* 324 * , 106 ( 2009 ) .
h.  h.  matevosyan , a.  w.  thomas and p.  c.  tandy , phys .  rev .
c * 75 * ( 2007 ) 045201 .
v.  a.  karmanov and p.  maris , few body syst .
* 46 * , 95 ( 2009 ) [ arxiv:0811.1100 [ hep - ph ] ] .
a.  krassnigg , phys .
d * 80 * , 114010 ( 2009 ) [ arxiv:0909.4016 [ hep - ph ] ] .
d.  nicmorus , g.  eichmann and r.  alkofer , phys .
d * 82 * , 114017 ( 2010 ) [ arxiv:1008.3184 [ hep - ph ] ] .
h.  l.  l.  roberts , l.  chang , i.  c.  cloet and c.  d.  roberts , few body syst .   *
51 * , 1 ( 2011 ) [ arxiv:1101.4244 [ nucl - th ] ] .
h.  sanchis - alepuz , r.  williams and r.  alkofer , phys .
d * 87 * , no .
9 , 096015 ( 2013 ) [ arxiv:1302.6048 [ hep - ph ] ] . c.  s.  fischer , s.  kubrak and r.  williams , eur .
phys .  j.  a * 50 * , 126 ( 2014 ) [ arxiv:1406.4370 [ hep - ph ] ] .
d.  dudal , h.  verschelde and s.  p.  sorella , phys .
b * 555 * , 126 ( 2003 ) .
n.  wschebor , int .  j.  mod
.  phys .
a * 23 * , 2961 ( 2008 ) .
j.  a.  gracey , phys .
b * 552 * , 101 ( 2003 ) .
r.  williams , arxiv:1404.2545 [ hep - ph ] .
a.  c.  aguilar , d.  binosi , d.  ibaez and j.  papavassiliou , phys .
d * 90 * , no . 6 , 065027 ( 2014 ) .
m.  mitter , j.  m.  pawlowski and n.  strodthoff , phys .
d * 91 * , 054035 ( 2015 ) [ arxiv:1411.7978 [ hep - ph ] ] .
j.  skullerud and a.  kizilersu , jhep * 0209 * , 013 ( 2002 ) .
j.  i.  skullerud , p.  o.  bowman , a.  kizilersu , d.  b.  leinweber and a.  g.  williams , jhep * 0304 * , 047 ( 2003 ) .
d.  zwanziger , nucl .
b * 323 * , 513 ( 1989 ) .
d.  zwanziger , nucl .
b * 399 * , 477 ( 1993 ) .
d.  zwanziger , phys .  rev .
d * 65 * , 094039 ( 2002 ) .
d.  dudal , j.  a.  gracey , s.  p.  sorella , n.  vandersickel and h.  verschelde , phys .
d * 78 * , 125012 ( 2008 ) .
j.  berges , n.  tetradis and c.  wetterich , phys .
* 363 * , 223 ( 2002 ) [ hep - ph/0005122 ] .
a.  sternbeck , l.  von smekal , d.  b.  leinweber and a.  g.  williams , pos lat * 2007 * , 340 ( 2007 ) .
a.  cucchieri and t.  mendes , phys .
lett .   * 100 * , 241601 ( 2008 ) .
a.  cucchieri and t.  mendes , phys .
d * 78 * , 094503 ( 2008 ) .
a.  sternbeck and l.  von smekal , eur .
j.  c * 68 * , 487 ( 2010 ) .
a.  cucchieri and t.  mendes , phys .
d * 81 * , 016005 ( 2010 ) .
i.  l.  bogolubsky , e.  m.  ilgenfritz , m.  muller - preussker and a.  sternbeck , phys .
b * 676 * , 69 ( 2009 ) .
d.  dudal , o.  oliveira and n.  vandersickel , phys .
d * 81 * , 074505 ( 2010 ) .
o.  bowman , u.  m.  heller , d.  b.  leinweber , m.  b.  parappilly and a.  g.  williams , phys .
d * 70 * , 034509 ( 2004 ) .
m.  q.  huber and l.  von smekal , jhep * 04 * , 149 ( 2013 ) . c.  becchi , a.  rouet and r.  stora , commun .  math .
phys .   * 42 * , 127 ( 1975 ) . c.  becchi , a.  rouet and r.  stora , annals phys .
* 98 * , 287 ( 1976 ) .
i.  v.  tyutin , arxiv:0812.0580 [ hep - th ] .
a.  cucchieri , t.  mendes and a.  r.  taurines , phys .
d * 71 * , 051902 ( 2005 ) [ hep - lat/0406020 ] .
o.  bowman , u.  m.  heller , d.  b.  leinweber , m.  b.  parappilly , a.  sternbeck , l.  von smekal , a.  g.  williams and j.  b.  zhang , phys .
d * 76 * , 094505 ( 2007 ) [ hep - lat/0703022 [ hep - lat ] ] . m.  tissier and n.  wschebor , phys .
d * 84 * ( 2011 ) 045018 .
j.  i.  skullerud , p.  o.  bowman , a.  kizilersu , d.  b.  leinweber and a.  g.  williams , nucl .
suppl .   * 141 * , 244 ( 2005 ) .
a.  sternbeck , k.  maltman , m.  muller - preussker and l.  von smekal , pos lattice * 2012 * ( 2012 ) 243 .


understanding the origin of spiral arms , bars , and rings is key for discerning the long - term evolution of disk galaxies .
the density wave ( dw ) theory , originally conceived by  @xcite and developed in a linear form by  @xcite , has been commonly brought up to explain spiral structure .
one fundamental prediction of this theory is the presence of offsets between tracers of the different stages of star formation and evolution  @xcite .
a manifestation of these offsets would be age / color gradients across arms , as a result of the interaction between a spiral pattern with an angular speed that is nearly constant with radius , and stars and gas with differential rotation .
the gas overtakes the spiral pattern inside corotation ( cr , i.e. , the radius where the pattern and orbital angular velocities are the same ) , and is left behind outside cr . as gas feels the gravitational potential near the spiral arms , it shocks  @xcite .
this produces the piling up of dust and molecular gas that is commonly found in the concave side of the arms , if the process occurs inside cr and spiral arms trail .
if the shock triggers star formation , an age ( or color ) gradient is expected across spiral arms as the new stars age and drift away from the arms .
young stars would be located immediately downstream the dust lanes , followed by increasingly older stars in the direction of rotation . in part due to the difficulties encountered to detect such gradients , the relevance of spiral dws for disk galaxies has been repeatedly challenged .
recent examples are @xcite ,  @xcite , and  @xcite .    in a series of papers , we have detected and studied azimuthal color / age gradients across spiral arms and bars , thus establishing a relation between spiral or bar dynamics and star formation  @xcite .
our studies suggest that azimuthal gradients can be found in almost every galaxy , but are difficult to detect and do not appear widespread within a single galaxy .
@xcite has noticed an inverse correlation between the signatures of a spiral shock , and the presence of chains of star complexes ( `` beads on a string '' ) in spiral arms ; he has also observed that spiral shocks seem to be associated with an irregular magnetic field  @xcite , whereas the  beads on a string " configuration correlates with the existence of a regular magnetic field . pc , irrespective of central column density  @xcite . ]    the `` infant mortality '' of star clusters  @xcite , that leads to their disruption in the first 10 - 30 myr , may prevent the formation of age gradients even if shocks trigger star formation .
stellar winds and supernovae may drive out the gas not used in star formation , causing a change in the gravitational potential .
this would lead to the disintegration of the clusters , as stars become unbound .
coherent age gradients must first survive the `` infant mortality '' process ; then , after @xmath4 myr , the `` dissolution of stellar groups '' scenario proposed by  @xcite can also take place near spiral arms .
these last process can explain the apparent `` downstream decline '' of the gradients already noticed by us  ( cf .
* ; * ? ? ?
* ) .
[ omegavar ]    recently , based on the analysis with the radial tremaine - weinberg ( rtw ) method of co and hi data of several galaxies , it has been proposed that the spiral pattern speed @xmath5 may actually increase with decreasing radius in some objects  @xcite . according to  @xcite and  @xcite
, this behavior is also seen in simulations
. occurs only in simulations with flocculent , transient arms . ]
if this is indeed the case , there would be no real long - lived pattern with solid body rotation , and hence no significant offsets between the tracers of the different stages of star formation would be expected ( see , for example , * ? ? ?
in particular , azimuthal color gradients across the arms would not be observed .    in  @xcite , we investigated the effect of having ( and neglecting ) non - circular streaming motions , that is , of using a purely circular dynamic model to derive spiral pattern speeds from observations of color gradients . summarizing , based on semi - analytical solutions and mhd simulations ,
gas orbits were obtained for models with a fixed pattern speed .
stellar population synthesis ( sps ) models of age gradients were incorporated to these gas orbits , and synthetic photometric observations were performed .
the method of  ( * ? ? ?
* hereafter gg96 ) was then applied to measure @xmath5 . in spite of having a real pattern speed , constant with radius
, higher pattern speeds were measured at smaller radii ( see figure  [ noncircular ] ) . as mentioned in  @xcite ,
the reason for this bias is that gas streamlines in a steady rotating spiral shock turn somewhat along the arm after passing through the shock .
hence , stars take longer to move away from the arm , and the observer will think there is a smaller difference between the pattern speed and the orbital frequency .
@xmath5 will seem to follow @xmath6 , just like in the previously mentioned applications of the radial rtw method , except that in our case we can be sure that the trend is caused by the streaming motions of the young stars , rather than by a radial variation of @xmath5 .
the bias is also found in real observations , as shown in  @xcite .    in an ideal scenario where @xmath5 is indeed constant for all radii and spiral arms trigger star formation in all their extension , color gradients are predicted to be located all along the arms ( except near cr ) . in reality ,
age gradients are detected only in a few small regions per galaxy , but they do seem linked to disk dynamics ( gg96 ; * ? ? ? * ; * ? ? ?
* ; * ? ? ?
so far , we have detected gradients through visual inspection . in order to locate gradients in a more objective way , and to attempt to identify some that do not meet the theoretical expectations of coherence or smoothness , in what follows
we will adopt fourier techniques to study normal ( or weakly barred ) spirals .
our data consist of deep photometric images of the 13 sa and sab spirals . ]
described in  @xcite .
the data include mosaics in the optical @xmath7 , @xmath8 , and @xmath9 bands  ( * ? ? ?
* ; * ? ? ?
* see table  [ tbl - filters ] ) , the near - infrared ( nir ) @xmath0 band , and the nir @xmath10 , @xmath11 , @xmath12 or @xmath13 filters .
we have now added six new objects , whose observation log is shown in table  [ tbl - obslog ] .
ngc  1703 and ngc  6907 , classified as sb by rc3 , were included in the analysis because ( apparently ) they have no prominent bar structure , as appreciated in the images .
the data were reduced and calibrated with the same standard procedures used in  @xcite , and  @xcite .
the optical ( @xmath7 , @xmath8 , and @xmath9 ) calibration was done in the thuan - gunn system  @xcite .
the zero point of this photometric system is chosen such that the standard star bd+17@xmath144708 has @xmath15 mag .
the nir @xmath0 data were calibrated with images from the two micron all sky survey  ( 2mass ; * ? ? ? * ) .
@xmath7 _ & 5000 & 830 + _ @xmath8 _ & 6800 & 1330 + _ @xmath9 _ & 7800 & 1420 + _ @xmath0 _ & 1.25 & 0.29 +    cccccl    ngc  1703 & @xmath7 & 5100 & ctio 0.9 m & 3/94 , 3/95 & type : sb(r)b + & @xmath8 & 4200 &  &  & pa ( deg ) : 104 + & @xmath9 & 5400 &  &  & @xmath16 ( deg ) : 27.0 @xmath17 12.1 + & @xmath0 & 2717 & ctio 1.5 m & 2/94 & @xmath18 ( km s@xmath19 ) : 55 @xmath17 20 + & @xmath11 & 1281 &  &  & dist ( mpc ) : 20.4 @xmath17 1.7 + & @xmath10 & 2758 &  &  & +    ngc  3001 & @xmath7 & 4500 & ctio 0.9 m & 3/94 , 3/95 & type : sab(rs)bc + & @xmath8 & 5100 &  &  & pa ( deg ) : 6 + & @xmath9 & 4500 &  &  & @xmath16 ( deg ) : 47.5 @xmath17 4.9 + & @xmath0 & 2718 & ctio 1.5 m & 2/94 & @xmath18 ( km s@xmath19 ) : 248 @xmath17 20 + & @xmath10 & 2700 &  &  & dist ( mpc ) : 35.6 @xmath17 3.0 +    ngc  3338 & @xmath7 & 3600 & lick 1 m & 4/94 , 11/94 & type : sa(s)c + & @xmath8 & 4500 &  &  & pa ( deg ) : 100 + & @xmath9 & 4200 &  & 2/94 , 4/94 , 11/94 & @xmath16 ( deg ) : 51.9 @xmath17 2.1 + & @xmath0 & 1200 & lick 1 m & 12/94 , 2/95 & @xmath18 ( km s@xmath19 ) : 194 @xmath17 9 + & @xmath0 & 1002 & kitt peak 1.3 m & 3/94 , 11/94 & dist ( mpc ) : 23.7 @xmath17 2.1 + & @xmath20 & 926 &  &  & +    ngc  4603 & @xmath7 & 5400 & ctio 0.9 m & 3/94 , 3/95 & type : sa(s)c + & @xmath8 & 4500 &  &  & pa ( deg ) : 27 + & @xmath9 & 4200 &  &  & @xmath16 ( deg ) : 43.6 @xmath17 5.7 + & @xmath0 & 1393 & ctio 1.5 m & 2/94 & @xmath18 ( km s@xmath19 ) : 242 @xmath17 32 + & @xmath10 & 1399 &  &  & dist ( mpc ) : 29.4 @xmath17 2.6 +    ngc  6907 & @xmath7 & 3600 & ctio 0.9 m & 9/94 & type : sb(s)bc + & @xmath8 & 3900 &  &  & pa ( deg ) : 46 + & @xmath9 & 3600 &  &  & @xmath16 ( deg ) : 35.6 @xmath17 3.7 + & @xmath0 & 570 & ctio 1.5 m & 9/94 , 9/95 & @xmath18 ( km s@xmath19 ) : 244 @xmath17 41 + & @xmath20 & 355 &  &  & dist ( mpc ) : 47.7 @xmath17 4.1 +    ngc  7083 & @xmath7 & 3600 & ctio 0.9 m & 9/94 & type : sa(s)bc + & @xmath8 & 3900 &  &  & pa ( deg ) : 5 + & @xmath9 & 3900 &  &  & @xmath16 ( deg ) : 52.9 @xmath17 3.0 + & @xmath0 & 840 & ctio 1.5 m & 9/94 & @xmath18 ( km s@xmath19 ) : 226 @xmath17 10 + & @xmath20 & 392 &  &  & dist ( mpc ) : 45.4 @xmath17 3.9 +
the dw theory predicts the existence of age gradients across spiral arms .
multi - armed or flocculent galaxies , however , are not described by this theory .
the short arms of the latter are probably unbound star complexes sheared by galactic differential rotation , whereas the `` knee - like '' . ] shaped arms of m101-type galaxies might be formed by the gravitational instability of galactic disks ; if so , they are transient features  @xcite .
the azimuthal profiles ( intensity versus @xmath21 ) of such arms should be symmetric about the gravitational potential well , where the gas is yet to be transformed into stars ; this configuration is indeed sometimes observed in certain galaxies ( for example , in ngc  4535 and m61 , a narrow dark lane is seen right at the center of a stellar arm ) .    to analyze the validity of the dw theory
, one should study galaxies whose structure is likely to be explained by it . mainly , from the point of view of their spiral arms , there are three different types of galaxies  ( e.g. , * ? ? ?
these are ( 1 ) symmetric grand - design spirals ( e.g. , m81 ) , with mass arms due to dws ; ( 2 ) multi - armed or `` knee - like '' spirals ( e.g. , m101 ) ; and ( 3 ) flocculent spiral galaxies ( e.g. , ngc  2841 ) . however , even in the same galaxy , different arms may be formed by different mechanisms .
thus , in the classical multi - armed galaxy m101 the inner arms are rather regular , and in our own galaxy the quite symmetrical inner arms transform into a polygonal multi - armed structure in the outer regions  @xcite .
albeit flocculent and multi - armed structures have a low mass - to - light ratio ( @xmath22 ) and are seen mainly in the optical ( see , for example , the classical example of ngc  309 in block et al .
1994 ) , in order to determine whether the arms in the nir are indeed mass dws it is important to disentangle the contribution of young stars and clusters at these longer wavelengths ; such contribution can reach up to 20%-30% ( e.g. , rix & rieke 1993 ; gg96 ; rhoads 1998 ; james & seigar 1999 ; patsis et al .
2001 ; grosbl et al .
2006 ) . with the aim of determining the actual structural type of the spiral arms , resolved stellar mass maps are required .
mass maps get around the problem that the @xmath22 , at any wavelength , is not constant across the disk and near spiral arms .    in order
to build resolved maps of stellar mass of our sample galaxies , we use the method of  @xcite .
the method relies on a monte carlo library of 50,000 stellar population spectra , constructed from the sps models of  @xcite , and charlot  & bruzual ( 2007 , private communication ) .
each spectrum is computed by randomly drawing the model parameters ( star formation history , metallicity , and dust attenuation ) from adequate physical distributions  ( see also * ? ? ?
we use the `` magphys '' code of  @xcite to extract colors from the spectral library .
each element in a theoretical two - dimensional color - color plot can be produced by several combinations of model parameters .
the median @xmath22 of all these combinations is then found , and attributed to the data with the same position in the observed two - dimensional color - color plot .
resolution elements in two - dimensional color space ( in our case , @xmath2 versus @xmath3 ) have a size of @xmath23 mag@xmath24 , in order to account for the observational error .
we compute the median @xmath22 in the @xmath0-band ( see figure  [ colorsml ] ) , and compare the resulting distributions with the colors obtained from the observations , zibetti / software/ adaptsmooth.html the code enhances the s / n with a minimum loss of effective resolution , and keeps the photometric fluxes unaltered .
the models are reddened before comparing with the observations , to take into account foreground galactic extinction  @xcite . ] on a pixel - by - pixel basis .
this comparison produces a resolved map of @xmath1 .
finally , a relative  mass map " is obtained , equal to the product of the @xmath1 map by the @xmath0-band image .
absolute mass maps can be worked out , taking into account the distance to the galaxies , but are not necessary for our purposes .
the results are shown in figures  [ j_vs_mass ] and  [ j_vs_mass_cont ] , where we compare the @xmath7 and @xmath0 frames with the mass maps at the same relative scale .
we then proceed to select the objects in our sample that retain a spiral structure that may be explained by the dw theory . from this
_ visual _ inspection we reject the following objects prior to the fourier analysis :    1 .
_ ngc  3162_. there is no two - armed spiral structure appreciable in the mass map .
ngc  3938_. the mass map has a very low quality , due to a low signal - to - noise ratio ( s / n ) of the @xmath7-band mosaic .
ngc  5371_. ring structures rather than spirals are discerned in the mass map .
ngc  7083_. the spiral structure looks flocculent in the mass map .
_ ngc  7126_. no clear spirals can be appreciated in the mass map .
_ ngc  7753_. spiral structure is not observed in the mass map .
the remaining 13 objects were analyzed as described below .
[ fourier ]      [ q_fourier ]    the first step in our fourier analysis consists in generating @xmath25 mosaics , following the method of gg96 .
@xmath26 is defined as follows :    @xmath27    where @xmath28 is the color excess term .
this photometric index is reddening - free for screen absorption , and also for a mixture of dust and stars , as long as @xmath29 . at the same time
, @xmath25 is sensitive to supergiant stars , and therefore traces star formation .
sps models  @xcite predict that , following a continuous burst of star formation that lasts for @xmath30 yr , the @xmath25 index increases its value during @xmath31 yr and then starts to decline  ( cf .
, gg96 ; * ? ? ?
* ) . by mass ) and
old stars ( with an age of @xmath32 yr ) .
both populations have solar metallicity and a salpeter initial mass function ( imf ) , with @xmath33 0.1 @xmath34 , and @xmath35 10 @xmath34 .
the choice of @xmath36 is justified by the observed inverse correlation between the locations where color gradients are detected , and those with important @xmath37 emission  ( gg96 ; * ? ? ?
* ; * ? ? ?
] @xmath25 mosaics of our sample were transformed to polar coordinates . vs.  ln @xmath38 ; logarithmic arms appear as straight lines in such a map , with slope @xmath39 , where @xmath9 is the arm pitch angle . ]
spiral arms were then `` straightened '' , by adding to the @xmath21 coordinate a phase as a function of @xmath40 , until the arms appear horizontal . unlike our previous studies , where the `` straightening slope '' was different depending on the specific arm and region under analysis , here we adopt the same `` straightening slope '' for both ( or all the ) arms in the same object . for each galaxy , this slope is obtained from the median of the pitch angles measured in different wavelengths  ( see section  [ pitch_angles ] and table  [ tbl - pitchs ] ) . , @xmath8 , @xmath9 , and @xmath0 from the lowest to the highest value , and averaging together the two middle values . ]
we divide the `` straightened '' images in vertical sections , equivalent to concentric annuli in a deprojected image of the galaxy , in order to perform a fourier analysis . here
we aim to compare the observed @xmath25 index profiles with the sps models adopting fourier techniques .
sps models give @xmath25 as a function of age , @xmath41 , while observations provide @xmath25 versus azimuthal angle ( or azimuthal distance after assuming a distance to the galaxy ) .
the age gradient signal can be thought of as an asymmetric `` pulse '' ( see figure  [ qmodel ] ) , whose width is a function of radius ; the width approaches zero at cr and has larger values at smaller radii .
the width of this pulse is analogous to a period , @xmath42 , of a sinusoidal signal .
hence , the pulse can also be characterized by its angular frequency @xmath43 , and analyzed with fourier techniques .
we compute the fourier transform of the @xmath25 index profiles in the annuli ,    @xmath44    where    @xmath45    in general , @xmath46 is complex , say , of the form @xmath47 .
the fourier spectrum is obtained as @xmath48 .
we search the spectrum for frequencies with amplitudes between 0.02 and 0.06 mag .
this range of amplitudes corresponds to the ones achieved by the @xmath25 index of our models , ) .
model `` a '' is obtained with the sps models only .
model `` b '' includes both the dissolution of stellar groups  @xcite after 50 myr , and stellar orbit diffusion . according to  @xcite
, the diffusion of stellar orbits can enhance the dissolution of young stellar groups by increasing their internal velocity dispersion .
] including the photometric error of the observed @xmath25 profiles , @xmath49 .
if more than one frequency is found , we then select the one with the highest frequency value .
profile , whereas higher frequencies correspond to small - scale structure . ]    the expected radial behavior of @xmath50 after applying the procedure to the mhd simulation of  @xcite is shown in figure  [ wplot_prediction ] , with a dotted line .
we remind the reader that the simulation has a spiral pattern with a constant angular speed that does not depend on radius .
the plot also shows ( solid line ) the behavior of @xmath50 when newborn stars move in perfectly circular orbits . in both cases ,
since the width of the pulse goes to zero at the cr radius , its reciprocal , @xmath50 , diverges there ; @xmath51 .    figures  [ wplusphases_1 ] through  [ wplusphases_7 ] , left panels , show the results of this analysis when applied to our sample of objects .
data @xmath52 values are graphed versus @xmath53 , where @xmath54 is the mean galactocentric radius , and @xmath55 is the spiral maximum radial extent . the latter is determined visually in the nir band indicated in table  [ tbl - armend ] .
llr    ngc  578 &  95.7 @xmath17 2.9 ( @xmath12 ) & 10.5 @xmath17 0.9 + ngc  918 &  75.4 @xmath17 5.8 ( @xmath13 ) & 7.9 @xmath17 0.7 + ngc  1417 &  60.0 @xmath17 5.0 ( @xmath0 ) & 16.6 @xmath17 1.4 + ngc  1421 &  92.8 @xmath17 2.9 ( @xmath12 ) & 13.2 @xmath17 1.1 + ngc  1703 &  27.5 @xmath17 1.4 ( @xmath10 ) & 2.7 @xmath17 0.3 + ngc  3001 &  44.0 @xmath17 2.8 ( @xmath10 ) & 7.6 @xmath17 0.8 + ngc  3338 & 145.0 @xmath17 5.0 ( @xmath20 ) & 16.7 @xmath17 1.6 + ngc  4254 & 157.5 @xmath17 7.5 ( @xmath12 ) & 12.6 @xmath17 0.8 + ngc  4603 &  45.4 @xmath17 1.4 ( @xmath10 ) & 6.5 @xmath17 0.6 + ngc  4939 & 145.0 @xmath17 5.0 ( @xmath12 ) & 32.7 @xmath17 2.8 + ngc  6907 &  55.1 @xmath17 2.9 ( @xmath20 ) & 12.7 @xmath17 1.3 + ngc  6951 & 105.0 @xmath17 5.0 ( @xmath12 ) & 12.7 @xmath17 1.1 + ngc  7125 &  95.7 @xmath17 2.9 ( @xmath0 ) & 20.7 @xmath17 1.8 +      [ phases ]    if dws indeed produce shocks and induce star formation , we expect increasingly older stellar populations ( with intensity peaks at longer wavelengths ) to be sequentially located downstream ( in the direction of rotation ) from the shock position within cr , and to be sequentially placed upstream ( opposite rotation ) outside cr .
no gradients should be observed at cr .
we use also fourier techniques to analyze our sample in a systematic manner , at all galactic radii , in search of this effect .
we adopt the method of  @xcite to determine the phases ( basically , the conjugate of the azimuthal positions ) of stellar populations of different ages , as traced in the @xmath7 , @xmath8 , @xmath9 , and @xmath0 bands .
the method is based on computing the fourier transform of the form    @xmath56    where @xmath57 is the intensity of radiation , and with phase    @xmath58}{\mathrm{im}[\hat{f}(m ) ] } \right\}.\ ] ]    given that we have selected objects with a strong two - armed pattern , @xmath59 .
this means that 2@xmath60 radians include two periods of the analyzed signal , i.e. , there is a @xmath60 radian symmetry .
although this may not be exactly true for real arms in spirals ( p.
grosbl 2012 , private communication ) , it is a very good approximation for our purposes .
we analyze the same concentric annuli in the `` straightened '' images as before .    for the data analysis
we assume all spirals in our sample trail , and hence that the spiral opens in the direction opposite rotation .
furthermore , we adopt the convention that the angular coordinate increases with rotation , regardless of whether the arms have an `` s '' or a `` z '' on - the - sky view . under these premises , for a wave of the form @xmath61 , where @xmath62 is the phase , @xmath62 should decrease with wavelength inside cr and increase outside cr .
there should be no phase difference with wavelength at cr ( see figure  [ phaseplot_prediction ] ) .
the results of this exercise are shown in figures  [ wplusphases_1 ] through  [ wplusphases_7 ] , right panels .
for seven objects ( ngc  918 , ngc  1417 , ngc  1421 , ngc  3001 , ngc  3338 , ngc  4254 , and ngc  7125 , see table  [ tbl - crll ] ) , we find similarities between the theoretical expectations and the observations , regarding both the radial distribution of @xmath50 and the radial run of intensity phases with wavelength , at least for some range of radii .
however , all the objects also show discrepancies . for example , in ngc  918 , ngc  1417 , ngc  1421 , ngc  3338 , ngc  4254 , and ngc  7125 , the phases of all four passbands overlap between @xmath63 0.2 and 0.4 @xmath64 , but then behave in accordance to theory up to what is likely the cr radius , as judged by the position where @xmath50 diverges and where , again , phases are the same for all wavelengths . in ngc  918 and ngc  4254 , phase differences agree with theory within cr , but then do not change direction
. this could be due to a weaker spiral potential and , hence , weaker shocks at radii beyond cr .
conversely , in ngc  3001 the phases show the expected behavior only outside the cr position , if it is indeed located where the @xmath52 plot shows a clear divergence . the intensity phases in ngc  1703 ( which may be an sb type galaxy ) behave as predicted by theory , but @xmath52 does not diverge at the cr position expected from the phases .    from this analysis ,
there is not even one object for which both diagnostics agree perfectly with theory at all radii .
these findings confirm more objectively a previous result from visual inspection : age / color textbook gradients that run coherently along entire arms and change direction at cr are very hard to detect in practice .
possible causes may be a low s / n of the data , or masking of the gradients by the intense and clumpy ( continuum and line ) radiation from massive star clusters . in the case of the intensity phases method , dusty environments can attenuate the shorter wavelengths , hindering the detection of an ordered age sequence .. ] it is also possible that the gradients are perturbed by the `` infant mortality '' of star clusters  @xcite , or by the destruction of an ordered age front by supernovae shocks and stellar winds @xcite . on the other hand , even in the presence of dws , it is possible that no gradients are produced . for example , there may be inadequate physical conditions for the creation of the gradients , related to the strengths of the shock or of the magnetic field @xcite .
there is , of course , the possibility that a real constant pattern speed for all radii does not exist , i.e. , that the spiral structure is not explained by the dw theory .
our subsample , however , was selected explicitly to minimize this eventuality .
furthermore , we had previously established the link of the gradients in the sample , taken as an ensemble , to disk dynamics . as confirmed by the comparison with an mhd simulation with a constant spiral pattern speed , an apparent radial dependency of @xmath5
is caused by neglecting non - circular streaming motions in the analysis of the gradients .
conversely , there could be color gradients that are not linked to the dynamics of the disk .
one speculative scenario is provided by the igimf ( the integrated galaxial igimf ) theory of  @xcite and  @xcite .
this theory is able to explain the h@xmath16 cutoff in galaxy disks without a corresponding cutoff in star formation ( as deduced from the non - ionizing far - ultraviolet flux given out by b - type stars ) .
the igimf theory posits that maximum star cluster mass depends on gas surface density ( or , equivalently , on star formation rate ) , and that maximum stellar mass depends on cluster mass .
if this is indeed the case , one could obtain a color gradient from coeval star clusters across a region with a gas surface density gradient .    at any rate , a statistical analysis , tied to spiral dynamics , of color gradients in several objects
is necessary to identify the gradients that are linked to dws .
this is even more important since by general rule gradients are painstakingly identified only in selected regions and not along entire arms .
it is , however , important to stress that our new results imply that long - lived modes seem to have produced gradients that run for a sizeable part of a spiral arm in at least @xmath65 of the analyzed objects .
this is , in itself , remarkable .
for the rest of the objects ( ngc  578 , ngc  1703 , ngc  4603 , ngc  4939 , ngc  6907 , and ngc  6951 ) , the fourier methods failed to detect signs of gradients besides the ones that had already been identified by eye .
it is interesting to point out , though , that ngc  4603 shows a multi - armed structure in the optical and nir bands , although it was kept in the sample because of its two - armed structure in the mass map . for ngc  6951 ( type sab in the rc3 ) and ngc  6907 ( type sb in the rc3 )
, the cr radius may lie around @xmath66 .
this is possibly an indication of color gradients within the bar region , as already found in  @xcite .
we find no relation between hubble type and success rate in our fourier tests .
lccr    ngc  918 & 0.63 & @xmath67 & @xmath68 + ngc  1417 & 0.67 & @xmath69 & @xmath70 + ngc  1421 & 0.49 & @xmath71 & @xmath72 + ngc  3001 & 0.55 & @xmath73 & @xmath74 + ngc  3338 & 0.38 & @xmath75 & @xmath76 + ngc  4254 & 0.46 & @xmath77 & @xmath78 + ngc  7125 & 0.67 & @xmath79 & @xmath80 +      [ pitch_angles ]    as a complementary test of dw presence , we examine the pitch angles of the spiral arms in different wavebands ( @xmath7 , @xmath8 , @xmath9 , and @xmath0 ) , and in the `` resolved mass maps '' previously obtained with the method of  @xcite . as already mentioned in the introduction , according to  @xcite , age gradients are observed across arm segments that fulfill the conditions for a strong shock to occur .
these segments seem to have an irregular magnetic field , while the regions with a regular magnetic field harbor instead weak shocks , and chains of star complexes with a nearly regular spacing .
strong shocks are associated with a larger pitch angle of the spiral segment .
another prediction involving pitch angles in the dw theory paradigm comes from the local dispersion relation . at a fixed mass surface density ,
the pitch angle increases proportionally to the square of the velocity dispersion  ( see , e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ) , that is , with the age of the stellar population and , hence , with image wavelength in the optical and nir .
pitch angles were measured with fourier techniques  @xcite , in the radial range where the spiral perturbation seems more prominent ( see table  [ tbl - pitchs ] ) .
this method assumes a logarithmic geometry for the spirals , although this might not be the case for all objects at all radii  @xcite .
figure  [ histo_j_gri ] shows histograms of @xmath81 - @xmath82 , where @xmath83 is the pitch angle and @xmath84 can be @xmath7 , @xmath8 , @xmath9 , or `` mass '' , depending on the image where @xmath83 has been measured . were adopted for analysis . however , @xmath81 values were discarded if their sign does not agree with that of @xmath82 .
( signs indicate whether the galaxy has an `` s '' or a `` z '' on - the - sky view . )
the pitch angles measured for ngc  578 with the fourier method were also rejected , since their nominal sign is at odds with the observed on - the - sky view . ]
negative values of the differences indicate that the pitch angle is larger in the nir than in the optical . in the histograms of @xmath85 , @xmath86 , and @xmath87 , there are indeed more objects with negative differences .
the median values of each histogram are , respectively , -0.9@xmath14 , -1.1@xmath14 , and 0.0@xmath14 .
this result is in agreement with  @xcite , who find tighter arms in bluer colors in images of four spirals , suggesting the presence of dws .
interestingly , the median value of the @xmath88 histogram is 2.2@xmath14 , so that the mass arms are statistically more open than the arms in nir light .
an inspection of the mass maps indicates that the spiral arms loci in the nir images do not always coincide with the location of the arms in the mass maps  ( see figures  [ j_vs_mass ] and  [ j_vs_mass_cont ] ) , a fact that can be explained by the contribution of young stars to the nir light .
finally , we compare the median pitch angle distribution in those galaxies that we have classified as having radially extended azimuthal color gradients , against the distribution in the galaxies where we have detected localized gradients only .
figure  [ histo_longlived ] displays the pitch angle histograms for each data set .
the statistical significance is low , due to the small number of objects , but the objects with signs of widespread color gradients may have more open spiral arms ( larger pitch angles ) .    crrrrrrc    ngc  578 & ( 20.3 - 63.8 ) & -34.70 ( @xmath89=2 ) & -35.72 ( @xmath89=2 ) & -37.92 ( @xmath89=2 ) & -41.68 ( @xmath89=2 ) & -49.45 ( @xmath89=2 ) & 40.97 + ngc  918 & ( 18.6 - 57.2 ) & 20.60 ( @xmath89=2 ) & 21.77 ( @xmath89=2 ) & 22.63 ( @xmath89=2 ) & 23.54 ( @xmath89=2 ) & 29.42 ( @xmath89=2 ) & 22.20 + ngc  1417 & ( 20.0 - 60.0 ) & 36.82 ( @xmath89=2 ) & 33.76 ( @xmath89=2 ) & 36.82 ( @xmath89=2 ) & 41.71 ( @xmath89=2 ) & 39.13 ( @xmath89=2 ) & 36.82 + ngc  1421 & ( 34.4 - 91.6 ) & 36.17 ( @xmath89=2 ) & 37.17 ( @xmath89=2 ) & 38.21 ( @xmath89=2 ) & 39.31 ( @xmath89=2 ) & -81.66 ( @xmath89=2 ) & 37.69 + ngc  1703 & ( 12.4 - 27.5 ) & 16.55 ( @xmath89=2 ) & 16.29 ( @xmath89=2 ) & 15.79 ( @xmath89=2 ) & -36.16 ( @xmath89=1 ) & -55.62 ( @xmath89=2 ) & 16.04 + ngc  3001 & ( 19.2 - 44.0 ) & 31.55 ( @xmath89=2 ) & 31.55 ( @xmath89=2 ) & 32.45 ( @xmath89=2 ) & 32.45 ( @xmath89=2 ) & 34.40 ( @xmath89=2 ) & 32.00 + ngc  3338 & ( 45.0 - 145.0 ) & -13.43 ( @xmath89=2 ) & -13.60 ( @xmath89=2 ) & -13.77 ( @xmath89=2 ) & -13.77 ( @xmath89=2 ) & 83.87 ( @xmath89=2 ) & -13.68 + ngc  4254 & ( 30.0 - 157.5 ) & -29.30 ( @xmath89=3 ) & 18.13 ( @xmath89=1 ) & -29.30 ( @xmath89=3 ) & -29.30 ( @xmath89=3 ) & 22.62 ( @xmath89=1 ) & -29.30 + ngc  4603 & ( 13.8 - 45.4 ) & 25.61 ( @xmath89=2 ) & 26.87 ( @xmath89=2 ) & 27.54 ( @xmath89=2 ) & 28.25 ( @xmath89=2 ) & 30.59 ( @xmath89=2 ) & 27.20 + ngc  4939 & ( 50.0 - 145.0 ) & 10.96 ( @xmath89=2 ) & 11.07 ( @xmath89=2 ) & 11.07 ( @xmath89=2 ) & 10.74 ( @xmath89=2 ) & 29.89 ( @xmath89=3 ) & 11.02 + ngc  6907 & ( 34.8 - 55.1 ) & 24.89 ( @xmath89=2 ) & 24.34 ( @xmath89=2 ) & 27.34 ( @xmath89=2 ) & 27.34 ( @xmath89=2 ) & 45.15 ( @xmath89=2 ) & 26.12 + ngc  6951 & ( 45.0 - 105.0 ) & -31.48 ( @xmath89=2 ) & -33.26 ( @xmath89=2 ) & -32.35 ( @xmath89=2 ) & -34.23 ( @xmath89=2 ) & -36.30 ( @xmath89=2 ) & -32.80 + ngc  7125 & ( 20.3 - 95.7 ) & -23.88 ( @xmath89=2 ) & -26.76 ( @xmath89=2 ) & -35.99 ( @xmath89=2 ) & -40.84 ( @xmath89=2 ) & -52.36 ( @xmath89=2 ) & -31.38 +
[ conclusions ]    from a fourier analysis of 13 objects , with a spiral arm structure that can be likely explained by dw theory , we find that @xmath90 of them show evidence of color gradients that run for at least a range of radii . for the remaining objects ,
the evidence of extended gradients is insufficient .
this notwithstanding , our previous studies  @xcite suggest that azimuthal age / color gradients ( or candidates ) can be found for most galaxies , at least in some regions of the spiral arms .
a relation of such very localized gradients with the dynamics of the disk has been previously established by the comparison with an mhd simulation with a constant pattern speed at all radii  @xcite .
our results are consistent with other observational studies of age patterns ( or offsets ) across spiral arms  @xcite . by means of an h@xmath16 to far - uv flux ratio method ,  @xcite found age gradients across the spiral arms of the grand - design spirals m74 and m100 ( ngc  628 and ngc  4321 , respectively ) .
likewise , the `` breaks '' in the radial metallicity distribution found near cr in a sample of @xmath63 20 galaxies by  @xcite and  @xcite imply that spiral arms must be long - lived structures .
a dominant pattern speed must exist with a unique cr radius , or otherwise any discontinuities in the radial metallicity profiles would be smoothed out  @xcite .
on the other hand , simulations have shown mutually exclusive results .
long - lived quasi - steady features have been obtained for some models  ( e.g. , * ? ? ?
* ; * ? ? ?
other works , however  ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ) , have found a lack of a persistent pattern , and instead recurrent short - lived transient spirals only .    any model that intends to explain spiral structure has to avoid the winding dilemma , that ultimately led to the idea of stationary dws with a fixed pattern speed for all radii .
for this same reason , the observational results of  @xcite and  @xcite for some objects , where @xmath91 decreases with radius , imply the existence of transient spiral modes that emerge and disappear at different stages in the disk evolution .
such transient arms can not produce stellar age / color gradients , though .    a theory that reconciles the existence of age / color gradients with a radially varying pattern speed is still missing .
however , as argued in section  [ omegavar ] , the observation of a pattern speed that varies with radius may be an artifact produced by the non - circular streaming motions of the stars newly born in shocked material .
this effect can explain the coexistence of azimuthal age / color gradients across spiral arms , and of a spiral pattern speed @xmath5 that seems to vary with radius in the same way as the orbital frequency @xmath6 .
we acknowledge dr .
yuri efremov , the referee , for his comments and suggestions that have greatly improved this paper .
special thanks go to ivanio puerari , for helpful discussions about fourier analysis of spiral galaxies .
e.m.g . acknowledges a conacyt grant for a postdoctoral fellowship at inaoe , and former postdoctoral financial support from unam ( dgapa ) , mxico .
is also grateful for financial support from both conacyt and dgapa - unam .
arzoumanian , d. , andr , p. , didelon , p. , et al .
2011 , , 529 , l6 athanassoula , e. , romero - gmez , m. , bosma , a. , & masdemont , j.  j.  2010 , , 407 , 1433 block , d.  l. , bertin , g. , stockton , a. , et al .
1994 , , 288 , 365 bruzual , g. , & charlot , s.  2003 , , 344 , 1000 clarke , c. , & gittins , d.  2006 , , 371 , 530 charlot , s. , & fall , s.  m.  2000 , , 539 , 718 considere , s. , & athanassoula , e.  1982 , , 111 , 28 considere , s. , & athanassoula , e.  1988 , , 76 , 365 da cunha , e. , charlot , s. , & elbaz , d.  2008 , , 388 , 1595 davis , b.  l. , berrier , j.  c. , shields , d.  w. , et al .
2012 , , 199 , 33 de vaucouleurs , g. , de vaucouleurs , a. , corwin , h.  g. , jr . , et al .  1991 , third reference catalogue of bright galaxies
( berlin : springer ) , ( rc3 ) dobbs , c.  l. , & bonnell , i.  a.  2008 , , 385 , 1893 dobbs , c.  l. , & price , d.  j.  2008 , , 383 , 497 donner , k.  j. , & thomasson , m.  1994 , , 290 , 785 efremov , y.  n.  1985 , soviet astronomy letters , 11 , 69 efremov , y.  n.  2010 , , 405 , 1531 efremov , y.  n.  2011 , astronomy reports , 55 , 108 egusa , f. , kohno , k. , sofue , y. , nakanishi , h. , & komugi , s.  2009 , , 697 , 1870 ferreras , i. , cropper , m. , kawata , d. , page , m. , & hoversten , e.  a.  2012 , , 424 , 1636 foyle , k. , rix , h .- w . ,
dobbs , c.  l. , leroy , a.  k. , & walter , f.  2011 , , 735 , 101 fujii , m.  s. , baba , j. , saitoh , t.  r. , et al .
2011 , , 730 , 109 gerola , h. , & seiden , p.  e.  1978 , , 223 , 129 gittins , d.  m. , & clarke , c.  j.  2004 , , 349 , 909 goodwin , s.  p. , & bastian , n.  2006 , , 373 , 752 gonzlez , r. a. , & graham , j. r. 1996 , , 460 , 651 ( gg96 ) gonzlez - lpezlira , r.  a. , pflamm - altenburg , j. , & kroupa , p.  2012
, , 761 , 124 grand , r.  j.  j. , kawata , d. , & cropper , m.  2012a , , 421 , 1529 grand , r.  j.  j. , kawata , d. , & cropper , m.  2012b , , 426 , 167 grosbl , p.  j. , & patsis , p.  a.  1998 , , 336 , 840 grosbl , p. , & dottori , h.  2009 , , 499 , l21 grosbl , p. , dottori , h. , & gredel , r.  2006 , , 453 , l25 hozumi , s.  2003 , in galaxies and chaos , ed . g. contopoulos & n. voglis ( lecture notes in physics , vol .
626 ; berlin : springer ) , 380 james , p.  a. , & seigar , m.  s.  1999 , , 350 , 791 lada , c.  j. , & lada , e.  a.  2003 , , 41 , 57 lin , c.  c. , & shu , f.  h.  1964 , , 140 , 646 lindblad , b.  1963 , stockholms observatoriums annaler , 22 , 5 martnez - garca , e.  e. , gonzlez - lpezlira , r.  a. , & bruzual - a , g.  2009a , , 694 , 512 martnez - garca , e.  e. , gonzlez - lpezlira , r.  a. , & gmez , g.  c.  2009b , , 707 , 1650 martnez - garca , e.  e. , & gonzlez - lpezlira , r.  a.  2011 , , 734 , 122 martnez - garca , e.  e.  2012 , , 744 , 92 meidt , s.  e. , rand , r.  j. , & merrifield , m.  r.  2009 , , 702 , 277 merrifield , m.  r. , rand , r.  j. , & meidt , s.  e.  2006 , , 366 , l17 mould , j. r.  et al .
2000 , , 529 , 786 mueller , m.  w. , & arnett , w.  d.  1976 , , 210 , 670 patsis , p.  a. , hraudeau , p. , & grosbl , p.  2001
, , 370 , 875 paturel , g. , fang , y. , petit , c. , garnier , r. , & rousseau , j.  2000 , , 146 , 19 paturel , g. , theureau , g. , bottinelli , l. , gouguenheim , l. , coudreau - durand , n. , hallet , n. , petit , c. 2003 , , 412 , 57 pflamm - altenburg , j. , & kroupa , p.  2008
, , 455 , 641 pflamm - altenburg , j. , weidner , c. , & kroupa , p.  2009
, , 395 , 394 pilbratt , g.  l. , riedinger , j.  r. , passvogel , t. , et al .
2010 , , 518 , l1 puerari , i. , & dottori , h.  a.  1992 , , 93 , 469 puerari , i. , & dottori , h.  1997 , , 476 , l73 rix , h.w . , & rieke , m. j. 1993 , , 418 , 123 rhoads , j.  e.  1998 , , 115 , 472 roberts , w.  w.  1969 , , 158 , 123 roca - fbrega , s. , valenzuela , o. , figueras , f. , et al .
2013 , , submitted ringermacher , h.  i. , & mead , l.  r.  2009 , , 397 , 164 snchez - gil , m.  c. , jones , d.  h. , prez , e. , et al .
2011 , , 415 , 753 saraiva schroeder , m.  f. , pastoriza , m.  g. , kepler , s.  o. , & puerari , i.  1994 , , 108 , 41 savchenko , s.  s.  2012 , astrophysical bulletin , 67 , 310 scarano , s. , jr .
, lpine , j.  r.  d. , & marcon - uchida , m.  m.  2011 , , 412 , 1741 scarano , s. , & lpine , j.  r.  d.  2013 , , 428 , 625 schlegel , d.  j. , finkbeiner , d.  p. , & davis , m.  1998 , , 500 , 525 sellwood , j.  a.  2011 , , 410 , 1637 speights , j.  c. , & westpfahl , d.  j.  2012 , , 752 , 52 skrutskie , m.  f. , cutri , r.  m. , stiening , r. , et al .
2006 , , 131 , 1163 tamburro , d. , rix , h .- w . ,
walter , f. , et al .
2008 , , 136 , 2872 thuan , t. x. , gunn , j. e. 1976 , , 88 , 543 wada , k. , baba , j. , & saitoh , t.  r.  2011 , , 735 , 1 wade , r. a. , hoessel , j. g. , elias , j. h. , huchra , j. p. 1979
, , 91 , 35 wielen , r.  1977 , , 60 , 263 zhang , x.  1998 , , 499 , 93 zibetti , s. , charlot , s. , & rix , h .- w .
2009 , , 400 , 1181 zibetti , s.  2009 , arxiv:0911.4956