because of its nonperturbative nature , studying low energy hadron physics in qcd is a long standing difficult problem . for low lying pseudoscalar mesons ,
a widely used approach is the theory of the effective chiral lagrangian based on the consideration of the global symmetry of the system and the momentum expansion without dealing with the nonperturbative dynamics of qcd @xcite@xcite . in the chiral lagrangian approach ,
the coefficients in the lagrangian are all unknown phenomenological parameters which should be determined by experimental inputs .
the number of unknown parameters increases rapidly with the increase of the precision in the momentum expansion .
therefore studying the relation between the chiral langrangian and the fundamental principles of qcd will not only be theoretically interesting for a deeper understanding of the chiral lagrangian , but will also be helpful for reducing the number of unknown parameters and increasing the predictive power of the chiral lagrangian .    in a previous paper , ref.@xcite , certain techniques were developed , with which the chiral lagrangian was formally derived from the first principles of qcd without taking approximations .
the chiral lagrangian coefficients are contributed both by the anomaly part ( from the quark functional measure ) and the normal part ( from the qcd lagrangian ) . in ref.@xcite
, all the chiral lagrangian coefficients contributed from the normal part of the theory are expressed in terms of certain green s functions in qcd , which can be regarded as exact qcd definitions of the chiral lagrangian coefficients .
after expanding the effective action in powers of the rotated sources ( momentum expansion ) , the effective action , up to @xmath2 , contributed from the normal part is of the form @xcite
@xmath3 ^ 2   -{\cal k}_2(d^{\mu}a_{\omega}^{\nu}-d^{\nu}a_{\omega}^{\mu } ) ( d_{\mu}a_{\omega,\nu}-d_{\nu}a_{\omega,\mu } ) \nonumber\\ & & \hspace{0.5 cm } + { \cal k}_3[a_{\omega}^2]^2 + { \cal k}_4a_{\omega}^{\mu}a_{\omega}^{\nu}a_{\omega,\mu}a_{\omega,\nu } + { \cal k}_5a_\omega^2{\rm tr}_f[a_{\omega}^2 ] + { \cal k}_6a_{\omega}^{\mu}a_{\omega}^{\nu } { \rm tr}_f[a_{\omega,\mu}a_{\omega,\nu}]+{\cal k}_7s_{\omega}^2 \nonumber\\ & & \hspace{0.5 cm } + { \cal k}_8s_{\omega}{\rm tr}_f[s_{\omega } ] + { \cal k}_9p_{\omega}^2 + { \cal k}_{10}p_{\omega}{\rm tr}_f[p_{\omega } ] + { \cal k}_{11}s_{\omega}a_{\omega}^2+{\cal k}_{12}s_{\omega}{\rm tr}_f [ a_{\omega}^2 ] -{\cal k}_{13}v_{\omega}^{\mu\nu}v_{\omega,\mu\nu } \nonumber\\ & & \hspace{0.5 cm } + i{\cal k}_{14}v_{\omega}^{\mu\nu}a_{\omega,\mu}a_{\omega,\nu } + { \cal k}_{15}p_{\omega}d_{\mu}a^{\mu}_{\omega}\bigg ] + o(p^6 ) , \label{p4}\end{aligned}\ ] ] where @xmath4 is related to the nonlinearly realized meson field @xmath5 by @xmath6 ; @xmath7 , and @xmath8 are , respectively , the external scalar , pseudoscalar , vector , and axial - vector sources rotated by @xmath4 ; and the @xmath9s are terms with different lorentz structures in the relevant qcd green s functions .
the obtained expressions for the chiral lagrangian coefficients are    @xmath10 : @xmath11[\bar{\psi}^b(x)\gamma_\mu\gamma_5\psi^a(x)]\rangle -\frac{1}{n_f}\langle [ \bar{\psi}^a(0)\gamma^\mu\gamma_5 \psi^a(0 ) ] [ \bar{\psi}^b(x)\gamma_\mu\gamma_5\psi^b(x)]\rangle\nonumber\\ & & \hspace{0.8cm}-\langle \bar{\psi}^a(0)\gamma^\mu\gamma_5\psi^b(0)\rangle \langle\bar{\psi}^b(x)\gamma_\mu\gamma_5\psi^a(x)\rangle + \frac{1}{n_f}\langle \bar{\psi}^a(0)\gamma^\mu\gamma_5\psi^a(0)\rangle \langle\bar{\psi}^b(x)\gamma_\mu\gamma_5\psi^b(x)\rangle\bigg ] , \nonumber\\ & & f_0 ^ 2b_0=-\frac{1}{n_f}\langle \bar{\psi}\psi\rangle , \label{f_0b_0}\end{aligned}\ ] ]      in the literature , the anomaly part contributions are usually calculated by means of the heat kernel regularization technique @xcite . however
, this technique is difficult to implement in the calculation of the normal part contributions which contain complicated functions of the momentum , say the quark self - energy @xmath1 , reflecting nonperturbative qcd dynamics ( which are even unspecified in the analytical part of the calculation ) . in order to treat the anomaly part and the normal part contributions on equal footing , a certain new regularization technique feasible for the calculations of both parts
should be developed . in this paper
, we use the generalized schwinger proper time regularization technique developed in ref.@xcite to regularize the system , which keeps the local chiral symmetry at every step in the calculation , and can be applied to the calculations of both the contributions from the anomaly part and from the normal part .
thus , in this paper , the contributions from the anomaly part and the normal part are calculated by means of the same technique . as the first conclusion of this unified treatment ,
we show that _ the anomaly contributions to the chiral lagrangian coefficients given in ref.@xcite , which are independent of qcd interactions , will actually be cancelled by certain terms in the normal part contributions , and the final expressions for the coefficients concern only the remaining terms from the normal part contributions related to qcd interactions_. it should be so since the coefficients indicate meson interactions which should be residual interactions between quarks and gluons , and thus should depend on qcd interactions .
these contributions have not been carefully calculated in the literature .
it has been shown in ref.@xcite that in the approximations of large-@xmath0 limit , leading order in dynamical perturbation , and improved ladder approximation , the formula for @xmath14 in eqs.([f_0b_0 ] ) reduces to the well - known pagels - stokar formula @xcite in which all dynamical effects from qcd are represented by the quark self - energy @xmath1 in the formula . in this paper , we take the same approximations to calculate the chiral lagrangian coefficients ( the relevant qcd green s functions ) as an illustration of the main feature of how qcd predicts the chiral lagrangian coefficients .
similar to the case of the pagels - stokar formula , the relevant qcd green s functions can all be expressed as functions of the quark self - energy @xmath1 . by solving the schwinger - dyson equation ,
we obtain @xmath1 , and thus the approximate qcd predicted values of the coefficients .
we shall see that the obtained coefficients @xmath15 and quark condensate are consistent with the experimental values .
the calculation is checked by the absence of divergences in the large-@xmath0 limit as it should be since the divergent meson - loop contributions are of next - to - the - leading order in the @xmath16 expansion .
although the present approximation is rather crude , _ it reveals the main feature of qcd predictions for the chiral lagrangian coefficients_.    this paper is organized as follows : in sec .
ii , we calculate the anomaly part contributions to the @xmath2 coefficients using the schwinger proper time regularization technique , and the results coincide with those in ref.@xcite in the chiral limit .
then , in sec .
iii , we apply the same technique to the normal part , and show generally that , in the large-@xmath0 limit , the anomaly part contributions to the chiral lagrangian coefficients are exactly cancelled by the contributions from a piece in the normal part independent of the quark self - energy , and the contributions from the remaining piece in the normal part depending on the quark self - energy play the real role in the chiral lagrangian coefficients .
specific approximations in the calculation of the normal part contributions and the formulae for the complete chiral lagrangian coefficients in terms of the quark self - energy are given in sec .
iv . in sec .
v , we present the numerical calculations of the quark self - energy and the obtained values of the chiral lagrangian coefficients .
section vi is a concluding remark .
in order to see the relation between the anomaly part and the normal part contributions to the chiral lagrangian coefficients , we present here the calculation of the anomaly part contributions by means of the schwinger proper time regularization .
we shall see that the obtained results exactly coincide with those obtained from the heat kernel technique @xcite .
our present approach is different from that in ref.@xcite in the sense that the constant constituent quark mass @xmath17 is not put in by hand as is done in ref.@xcite but is naturally included in the normal part solution through the dynamical quark mass reflecting chiral symmetry breaking
. therefore our result of the anomaly part contribution is to compare with that in ref.@xcite in the chiral limit .    in the schwinger proper time regularization
, the anomaly part does not contribute to the coefficients of the @xmath10 terms in the case corresponding to the result with @xmath18 in ref.@xcite to regularize the divergent integrals as was done in ref.@xcite before putting in the constituent quark mass @xmath17 , the @xmath10 coefficient @xmath14 will be proportional to @xmath19 ( cf .
ref.@xcite ) . as has been pointed out in ref.@xcite this term is exactly cancelled by a corresponding term in the normal part contribution [ cf .
eq.(74 ) in ref.@xcite ] . ] .
therefore we are only going to calculate the anomaly contribution to the coefficients of the @xmath2 terms .
the anomaly term in the path - integral formalism is @xmath20\nonumber\\ & = & in_c[{\rm tr}\ln(i\partial\!\!\!\ ! /\;+j_{\omega } ) + \omega-{\rm independent~term } ] .
\label{sanom}\end{aligned}\ ] ] the @xmath4-independent term is independent of the @xmath5 field , so that it is irrelevant to the chiral lagrangian coefficients .
we shall only evaluate the @xmath4-dependent term in eq.([sanom ] ) . to have a unified parametrization , we can parametrize the anomaly contributed effective action similar to that in eq.([p4 ] ) , i.e. , @xmath21 ^ 2   -{\cal k}_2^{\rm ( anom)}(d^{\mu}a_{\omega}^{\nu}-d^{\nu}a_{\omega}^{\mu } ) ( d_{\mu}a_{\omega,\nu}-d_{\nu}a_{\omega,\mu } ) + { \cal k}_3^{\rm ( anom)}[a_{\omega}^2]^2 \nonumber\\ & & + { \cal k}_4^{\rm ( anom)}a_{\omega}^{\mu}a_{\omega}^{\nu } a_{\omega,\mu}a_{\omega,\nu } + { \cal k}_5^{\rm ( anom)}a_\omega^2{\rm tr}_f[a_{\omega}^2 ] + { \cal k}_6^{\rm ( anom)}a_{\omega}^{\mu}a_{\omega}^{\nu } { \rm tr}_f[a_{\omega,\mu}a_{\omega,\nu}]+{\cal k}_7^{\rm ( anom)}s_{\omega}^2 \nonumber\\ & & + { \cal k}_8^{\rm ( anom)}s_{\omega}{\rm tr}_f[s_{\omega } ] + { \cal k}_9^{\rm ( anom)}p_{\omega}^2 + { \cal k}_{10}^{\rm ( anom)}p_{\omega}{\rm tr}_f[p_{\omega } ] + { \cal k}_{11}^{\rm ( anom)}s_{\omega}a_{\omega}^2 \nonumber\\ & & + { \cal k}_{12}^{\rm ( anom)}s_{\omega}{\rm tr}_f[a_{\omega}^2 ] -{\cal k}_{13}^{\rm ( anom)}v_{\omega}^{\mu\nu}v_{\omega,\mu\nu } + i{\cal k}_{14}^{\rm ( anom)}v_{\omega}^{\mu\nu}a_{\omega,\mu}a_{\omega,\nu } + { \cal k}_{15}^{\rm ( anom)}p_{\omega}d_{\mu}a^{\mu}_{\omega}\bigg ]   \nonumber\\ & & + o(p^6)+u-\mbox{independent~source~terms}. \label{p4anom}\end{aligned}\ ] ] the @xmath4-dependent term in eq.([sanom ] ) suffers from ultraviolet divergence , and we take the schwinger proper time regularization with an ultraviolet cutoff parameter @xmath22 to regularize it . to apply this regularization
, we first work in the euclidean space - time , and analytically continue the results to the minkowskian space - time after the evaluation . the main procedure of evaluating the general functional determinant including the quark self - energy @xmath23 is described in appendix a. in the case of @xmath24
, there is no @xmath23-dependence in eq.([sanom ] ) .
however , for regularizing the infrared divergence , we should replace the @xmath23 in eqs.([lnsigma ] ) and ( [ momexp ] ) by an infrared cutoff parameter @xmath25 . the momentum integration in eq.([momexp ] ) can be explicitly carried out with a lengthy but elementary calculation . after expanding in powers of the external sources
, we can identify the expressions for @xmath26 by comparing with the form of eq.([p4anom ] ) , and we obtain @xmath27    comparing with the standard form of momentum expansion to identify the @xmath2 chiral lagrangian coefficients , we obtain the anomaly contribution to these coefficients @xmath28 these exactly coincide with the results with @xmath18 in ref.@xcite .
note that the final expressions of the coefficients @xmath29 are independent of the infrared cutoff parameter @xmath25 and the ultraviolet cutoff @xmath22 although these cutoff parameters appear in @xmath30 , while @xmath31 and @xmath32 depend on the cutoff parameters .
this implies that @xmath31 and @xmath32 are not measurable quantities . with @xmath33 ,
the values of the coefficients are ( in units of @xmath34 ) @xmath35 these are to be compared with the experimental values ( in units of @xmath34 ) @xcite @xmath36 the numbers in eqs.([anomalyl ] ) are close to the experimental results of eqs.([exptl ] ) except @xmath37 and @xmath38 are of wrong signs .
this gives people an impression that the coefficients @xmath39 might mainly be contributed by the anomaly part , and the normal part might only contribute small corrections @xcite .
however , we note that the results in eqs.([anomalyresults ] ) are independent of qcd interactions , i.e. , these terms remain unchanged when we switch off the qcd gauge coupling constant @xmath40 .
this is somewhat confusing since these coefficients indicate meson interactions which should be residual interactions between quarks and gluons .
we shall see in the next section that _ these terms will actually be completely cancelled by the terms independent of the quark self - energy in the normal part contribution , so that they do not really appear in the final form of the coefficients .
what appear in the coefficients are the remaining terms in the normal part contribution which depend on the quark self - energy and hence on the qcd interactions as it should be . _
another feature of the terms in eqs.([anomalyresults ] ) indicating that they should be exactly cancelled and should not appear in the final formulae for the coefficients is the divergence of @xmath31 and @xmath32 when taking @xmath41 .
we know from ref.@xcite that the ultraviolet divergences in the @xmath2 chiral lagrangian coefficients come merely from the meson - loop corrections with the @xmath10 interactions . in the @xmath16 expansion , the meson - loop corrections belonging to @xmath42 will not take place in the large-@xmath0 limit .
therefore , in the large-@xmath0 limit , the final expressions for the @xmath2 coefficients should be finite when @xmath41 .
now the ultraviolet divergences in @xmath31 and @xmath32 in eqs.([anomalyresults ] ) have nothing to do with the meson - loop corrections , so that they should be exactly cancelled by other terms and should not appear in the final expressions for the @xmath2 coefficients .
in this section we use the same regularization technique as in sec .
ii to calculate the normal part contributions to the chiral lagrangian coefficients .
we start from the effective action @xmath43 given in ref.@xcite , @xmath44 } \nonumber\\ & & \hspace{0.4cm}=\int{\cal d}\xi{\cal d}\phi_{\omega}\exp\bigg\ { i\gamma_0[j_{\omega},\phi_{\omega},\pi_{\omega c}]+i\gamma_i[\phi_{\omega } ] + in_c\int d^4x{\rm tr}_{lf}\{\xi(x ) [ -i\sin\frac{\vartheta(x)}{n_f } + \gamma_5\cos\frac{\vartheta(x)}{n_f } ] \phi_{\omega}^t(x , x)\}\bigg\ } \label{seff}\end{aligned}\ ] ] which satisfies a useful relation @xcite @xmath45 the symbols are defined in ref.@xcite .    in the large-@xmath0 limit , the integrations in eq.([seff ] ) can be carried out by the saddle point approximation with the saddle point equations @xmath46^{(b\zeta)(a\eta)}(y , x ) , \label{pieq}\\ & & \pi_{\omega c}^{\sigma\rho}(x , y)=-\tilde{\xi}^{\sigma\rho}(x)\delta^4(x - y ) -\sum^{\infty}_{n=1}{\int}d^{4}x_1\cdots{d^4}x_{n } d^{4}x_{1}'\cdots{d^4}x_{n}'\frac{(-i)^{n+1}(n_c g_s^2)^n}{n!}\nonumber\\ & & \hspace{0.2cm}\times\overline{g}^{\sigma\sigma_1\cdots\sigma_n}_{\rho\rho_1 \cdots\rho_n}(x , y , x_1,x'_1,\cdots , x_n , x'_n ) \phi_{\omega c}^{\sigma_1\rho_1}(x_1 , x'_1)\cdots   \phi_{\omega c}^{\sigma_n\rho_n}(x_n , x'_n ) , \label{fineqnc}\\ & & { \rm tr}_l\bigg[\bigg(-i\sin\frac{\vartheta(x)}{n_f } + \gamma_5\cos\frac{\vartheta(x)}{n_f}\bigg)\phi^t_{\omega c}(x , x)\bigg]=0 , \label{constraint}\end{aligned}\ ] ] where @xmath47 is a short notation for the following quantity @xmath48\phi_{\omega c}^t(y , y)\ } \bigg|_{\xi_c\mbox { fixed}}.\label{xitilde}\end{aligned}\ ] ] then the obtained @xmath43 in this approximation is @xmath49 \nonumber\\ & = & -in_c{\rm tr}\ln[i\partial\!\!\!/+j_{\omega}-\pi_{\omega c } ] + n_c\int d^{4}xd^{4}x ' \phi^{\sigma\rho}_{\omega c}(x , x')\pi^{\sigma\rho}_{\omega c}(x , x ' ) + n_c\sum^{\infty}_{n=2}{\int}d^{4}x_1\cdots d^4x_{n}'\nonumber\\ & & \times \frac{(-i)^{n}(n_c g_s^2)^{n-1}}{n!}\bar{g}^{\sigma_1\cdots\sigma_n}_{\rho_1 \cdots\rho_n}(x_1,x'_1,\cdots , x_n , x'_n ) \phi^{\sigma_1\rho_1}_{\omega c}(x_1 , x'_1)\cdots   \phi^{\sigma_n\rho_n}_{\omega c}(x_n , x'_n)\nonumber\\ & & + in_c \int d^4x~{\rm tr}_{lf}\{\xi_c(x)[-i\sin\frac{\vartheta_c(x)}{n_f } + \gamma_5\cos\frac{\vartheta_c(x)}{n_f}]\phi_{\omega , c}^t(x , x)\ } , \label{snorm}\end{aligned}\ ] ] in which the @xmath42 term @xmath50 is neglected .
note that the last term in eq.([snorm ] ) actually vanishes due to eq.([constraint ] ) .
we keep it here for showing the relation between the effective action @xmath43 and its stationary conditions eqs.([pieq])@xmath51 ( [ constraint ] ) .    in the large-@xmath0 limit , @xmath52 on the right hand side of eq.([seffdiff ] ) . the left hand side of eq.([seffdiff ] )
can be carried out from eq.([snorm ] ) using eqs.([pieq])@xmath51([constraint ] ) .
then the explicit form of eq.([seffdiff ] ) in this approximation is @xmath53^{\rho\sigma}(x , x ) = \phi^{\sigma\rho}_{\omega c}(x , x ) .
\label{eq}\end{aligned}\ ] ] we see that @xmath54 and @xmath55 play the roles of the quark self - energy and the quark propagator , respectively , in the case with @xmath56
.    now we decompose @xmath43 into a part _ independent of @xmath54 _ and a part _ depending on @xmath54_. the part independent of @xmath54 can be extracted from @xmath43 by setting @xmath57 , i.e. , @xmath58 + n_c\bigg[\sum^{\infty}_{n=2}{\int}d^{4}x_1\cdots d^4x_{n } ' \frac{(-i)^{n}(n_c g_s^2)^{n-1}}{n!}\bar{g}^{\sigma_1\cdots\sigma_n}_{\rho_1 \cdots\rho_n}(x_1,x'_1,\cdots , x_n , x'_n)\nonumber\\ & & \times\phi^{\sigma_1\rho_1}_{\omega c}(x_1 , x'_1)\cdots   \phi^{\sigma_n\rho_n}_{\omega c}(x_n ,
x'_n)\bigg]_{\pi_{\omega c}=0}. \label{spi=0}\end{aligned}\ ] ] here we have ignored the last term in eq.([snorm ] ) which actually vanishes due to eq.([constraint ] ) .
we show in appendix b that the last term in eq.([spi=0 ] ) is actually @xmath4-independent
. therefore , eq.([spi=0 ] ) can be written as @xmath59 .
\label{spi=0'}\end{aligned}\ ] ] comparing the @xmath60-dependent terms in eqs.([sanom ] ) and ( [ spi=0 ] ) , we see that they are of the same form but with an opposite sign .
thus _ their contributions to the chiral lagrangian coefficients exactly cancel each other to all orders in the momentum expansion_. the cancellation in the case of the @xmath10 coefficient @xmath61 has been described in footnote 1 in sec .
ii . for the @xmath2 coefficients
, we have @xmath62 thus , in the large-@xmath0 limit , the anomaly part contributed chiral lagrangian coefficients in eqs.([anomalyresults ] ) do not really appear in the final results of the chiral lagrangian coefficients although their values eqs.([anomalyl ] ) are close to the experimental values . _ the chiral lagrangian coefficients are actually contributed from the @xmath63 part of @xmath43 _ , @xmath64 _ which leads to the @xmath65 part of @xmath66 _ , @xmath67 _ this is our first new conclusion in this study_.    the final chiral lagrangian coefficients are then @xmath68 and @xmath69\nonumber,\\ & & l_9 = \frac{1}{8}(4{\cal k}_{13}^{({\rm norm},\pi_{\omega c}\ne 0 ) } -{\cal k}_{14}^{({\rm norm},\pi_{\omega c}\ne 0)})\nonumber,\\ & & l_{10 } = \frac{1}{2}({\cal k}_2^{({\rm norm},\pi_{\omega c}\ne 0 ) } -{\cal k}_{13}^{({\rm norm},\pi_{\omega c}\ne 0)})\nonumber,\\ & & h_1 = -\frac{1}{4}({\cal k}_2^{({\rm norm},\pi_{\omega c}\ne 0 ) } + { \cal k}_{13}^{({\rm norm},\pi_{\omega c}\ne 0)})\nonumber,\\ & & h_2 = \frac{1}{8}[-{\cal k}_1^{({\rm norm},\pi_{\omega c}\ne 0 ) } + \frac{1}{b_0 ^ 2}{\cal k}_7^{({\rm norm},\pi_{\omega c}\ne 0 ) } + \frac{1}{b_0 ^ 2}{\cal k}_9^{({\rm norm},\pi_{\omega c}\ne 0 ) } -\frac{1}{b_0}{\cal k}_{15}^{({\rm norm},\pi_{\omega c}\ne 0 ) } ] .
\label{p4full}\end{aligned}\ ] ] since @xmath70 , these @xmath2 chiral lagrangian coefficients will vanish if we switch off the qcd coupling constant @xmath71 as it should be .
we see that to calculate the chiral lagrangian coefficients from @xmath72 , we should mainly deal with @xmath43 given in eq.([snorm ] ) which has never been carefully calculated in the literature . ignoring the vanishing last term in eq.([snorm ] )
, there are still rather complicated terms in it . for example , the third term includes various ranks of gluon green s functions which concern very complicated calculations of qcd dynamics . as the first time of doing this kind of calculation , we shall take futher approximations to simplify the evaluation of @xmath43 .
we know that the pion decay constant @xmath73 has been studied from qcd in ref.@xcite by taking the approximation of keeping only the leading order in dynamical perturbation , i.e. , taking into account only the qcd interaction in the schwinger - dyson equation leading to the nonperturbative solution of chiral symmetry breaking , and neglecting other qcd corrections in positive powers of @xmath71 ( perturbative ) .
this approximation leads to the widely used pagels - stokar formula which is reasonable though not perfect . in the large-@xmath0 limit
, @xmath73 is just the @xmath10 chiral lagrangian coefficient @xmath74 given in eqs.([f_0b_0 ] ) .
so , as in ref.@xcite , we take the approximation of keeping only the leading order in dynamical perturbation to calculate the chiral lagrangian coefficients from @xmath43 . in this spirit
, we neglect the complicated third term in eq.([snorm ] ) which contains only positive powers of @xmath71 .
furthermore , we see from eq.([fineqnc ] ) that the second term in eq.([snorm ] ) is of the same order as the third term , so that we neglect the second term in eq.([snorm ] ) as well . with this approximation ,
@xmath43 is simplified as @xmath75 .
\label{snormsimple}\end{aligned}\ ] ] now the concerned qcd dynamics resides in the @xmath54 term which is related to the quark self - energy [ cf .
eq.([eq ] ) ] .
we expect such a simple approximation may also lead to reasonable results of the @xmath2 chiral lagrangian coefficients since reasonable values of the @xmath2 chiral lagrangian coefficients have been obtained in a model by holdom @xcite considering only the quark self - energy contribution with certain phenomenological ansatz .
we shall see in sec .
v that our obtained @xmath2 chiral lagrangian coefficients are indeed reasonable .
although this approximation is crude , it provides a simple illustration of the main feature of how qcd predicts the chiral lagrangian coefficients . further improved study beyond this simple approximation taking into account the second and third terms in eq.([snorm ] ) is of course needed .
that will be presented in a later paper .
now we need to calculate @xmath76 and carry out the explicit expression for @xmath43 in eq.([snormsimple ] ) .
we have noticed that @xmath77 is related to the quark self - energy .
if we find out the relation between @xmath76 and the conventional quark self - energy @xmath78 , then we can obtain @xmath78 by solving the well - known schwinger - dyson equation .
as in the literature , we shall write down the schwinger - dyson equation in the landau gauge which is stable against the gauge parameter . in the same approximation of taking the leading order in dynamical perturbation theory and with the improved ladder approximation ,
the schwinger - dyson equation in the euclidean space - time reads @xcite @xmath79}{(p - q)^2 } \frac{\sigma(q^2)}{q^2+\sigma^2(q^2 ) } = 0 .
\label{eq0}\end{aligned}\ ] ] this equation can be solved numerically , and the details will be presented in sec .
v. naively , we may expect that @xmath80 .
but this is not correct . under a local chiral transformation @xmath81 ( hidden symmetry transformation @xcite ) , @xmath77 transforms as @xmath82 while @xmath83 does not transform like this
. the correct relation can be found by replacing the ordinary derivative @xmath84 by the covariant derivative @xmath85 since the external source @xmath86 transforms as @xmath87,\end{aligned}\ ] ] the covariant derivative @xmath88 transforms as @xmath89 thus the correct identification is @xmath90^{\sigma\rho}\delta^4(x - y ) .
\end{aligned}\ ] ] then the effective action ( [ snormsimple ] ) can be written as @xmath91 .
\label{sgl}\end{aligned}\ ] ]    next we evaluate the effective action ( [ sgl ] ) using the schwinger proper time regularization as before [ cf . appendix a ] to obtain the expressions for the chiral lagrangian coefficients .
this is not trivial since usually this regularization scheme is used in the case that @xmath23 is a constant , and thus the momentum integration can be explicitly carried out to check the local gauge invariance of the result .
now we leave @xmath92 as an unspecified function in eq.([sgl ] ) , so that the momentum integration can not be carried out explicitly .
organizing terms to guarantee local chiral invariance is rather tedious and the details are given in ref.@xcite .
our obtained expressions in the minkowskian space - time are @xmath93 , \label{f0}\\ & & { \cal k}^{({\rm norm})}_{1}=2\int d\tilde{p}\bigg[-2a_px_p^3 + 2a_p\frac{x_p^2}{\lambda^2}-a_p\frac{x_p}{\lambda^4 } + \frac{p^2}{2}\sigma'^2_p\frac{x_p}{\lambda^2}-\frac{p^2}{2}\sigma'^2_px_p^2 , \bigg ] , \nonumber\\ & & { \cal k}^{({\rm norm})}_{2}=\int d\tilde{p}\bigg[-2b_px_p^3   + 2b_p\frac{x_p^2}{\lambda^2}-b_p\frac{x_p}{\lambda^4 } + \frac{p^2}{2}\sigma^{\prime 2}_p\frac{x_p}{\lambda^2 } , -\frac{p^2}{2}\sigma^{\prime 2}_px_p^2\bigg ] , \nonumber\\ & & { \cal k}^{({\rm norm})}_{3}=2\int d\tilde{p}\bigg [ ( \frac{4\sigma^4_p}{3}-\frac{2p^2\sigma^2_p}{3}+\frac{p^4}{18 } ) (   6x_p^4-\frac{6x_p^3}{\lambda^2}+\frac{3x_p^2}{\lambda^4 } -\frac{x_p}{\lambda^6 } ) , \nonumber\\ & & \hspace{0.4cm}+(-4\sigma^2_p+\frac{p^2}{2})(-2x_p^3 + \frac{2x_p^2}{\lambda^2}-\frac{x_p}{\lambda^4 } ) -\frac{x_p}{\lambda^2}+x_p^2\bigg ] , \nonumber\\ & & { \cal k}^{({\rm norm})}_{4}=\int d\tilde{p}\bigg [ ( \frac{-4\sigma^4_p}{3}+\frac{2p^2\sigma^2_p}{3}+\frac{p^4}{18 } ) ( 6x_p^4-\frac{6x_p^3}{\lambda^2}+\frac{3x_p^2}{\lambda^4 } -\frac{x_p}{\lambda^6 } ) + 4\sigma^2_p(-2x_p^3+\frac{2x_p^2}{\lambda^2 } \nonumber\\ & & \hspace{0.4cm}-\frac{x_p}{\lambda^4})+\frac{x_p}{\lambda^2}-x_p^2\bigg ] , \nonumber\\ & & { \cal k}^{({\rm norm})}_5={\cal k}^{({\rm norm})}_6=0 , \nonumber\\ & & { \cal k}^{({\rm norm})}_7=2\int d\tilde{p}\bigg[(3\sigma^2_p+2p^2\sigma_p\sigma'_p)x_p^2 + [ -2\sigma^2_p - p^2(1 + 2\sigma_p\sigma'_p)]\frac{x_p}{\lambda^2}\bigg ] , \nonumber\\ & & { \cal k}^{({\rm norm})}_8=0 , \nonumber\\ & & { \cal k}^{({\rm norm})}_9=2\int d\tilde{p}\bigg[(\sigma^2_p+2p^2\sigma_p\sigma'_p)x_p^2 -p^2(1 + 2\sigma_p\sigma'_p)\frac{x_p}{\lambda^2}\bigg ] , \nonumber\\ & & { \cal k}^{({\rm norm})}_{10}=0 , \nonumber\\ & & { \cal k}^{({\rm norm})}_{11}=4\int d\tilde{p } \bigg[(-4\sigma^3_p+p^2\sigma_p)x_p^3 + ( 4\sigma^3_p - p^2\sigma_p)\frac{x_p^2}{\lambda^2 } -(2\sigma^3_p-\frac{1}{2}p^2\sigma_p)\frac{x_p}{\lambda^4 } + 3\sigma_p\frac{x_p}{\lambda^2 } \nonumber\\ & & \hspace{0.4cm}-3\sigma_p x_p^2\bigg ] , \nonumber\\ & & { \cal k}^{({\rm norm})}_{12}=0 , \nonumber\\ & & { \cal k}^{({\rm norm})}_{13}=\int d\tilde{p}\bigg [ ( \frac{1}{3}p^2\sigma'_p\sigma''_p+\frac{1}{3}\sigma_p\sigma''_p)x_p + ( c_p - d_p)\frac{x_p}{\lambda^2 } -(c_p - d_p)x_p^2 - 2e_px_p^3 \nonumber\\ & & \hspace{0.4cm}+2e_p\frac{x_p^2}{\lambda^2 } -e_p\frac{x_p}{\lambda^4}\bigg ] , \nonumber\\ & & { \cal k}^{({\rm norm})}_{14}=-4\int d\tilde{p}\bigg [ -2f_px_p^3 + 2f_p\frac{x_p^2}{\lambda^2 } -f_p\frac{x_p}{\lambda^4 } + \frac{p^2}{2}\sigma_p^{\prime 2}\frac{x_p}{\lambda^2 } -\frac{p^2}{2}\sigma^{\prime 2}_px_p^2\bigg ] , \nonumber\\ & & { \cal k}^{({\rm norm})}_{15}=-4\int d\tilde{p}\bigg [ -(\sigma_p+\frac{1}{2}p^2\sigma'_p)\frac{x_p}{\lambda^2 } + ( \sigma_p+\frac{1}{2}p^2\sigma'_p)x_p^2\bigg ] , \label{kresult}\end{aligned}\ ] ] in which the short notations ( in the minkowskian space - time ) are @xmath94    for the coefficient @xmath14 , eq.([f0 ] ) is just the well - known pagels - stokar formula @xcite when taking the regularization cutoff parameter @xmath95 .
it is easy to check that these @xmath96 ( @xmath97 ) do contain the @xmath54-independent ( @xmath98-independent ) piece which exactly cancel the anomaly contributions in eqs.([anomk ] ) mentioned in sec .
this can be done by taking a constant @xmath98 to carry out the momentum integrations , and picking up the @xmath98-independent terms which are just the @xmath54-independent terms mentioned in sec .
subtracting these @xmath54-independent terms from the obtained @xmath96 in eqs.([kresult ] ) , we get the desired @xmath99 in eq.([kpi ] ) , which is needed in the final expressions for the chiral lagrangian coefficients in eqs.([p4full ] ) .    we can also check that the regularization cutoff @xmath22 does not appear in @xmath99 , so that the obtained chiral lagrangian coefficients @xmath100 are all finite as it should be since there is no divergence in the large-@xmath0 limit [ the divergent meson - loop corrections are of @xmath42 ] .
@xmath125    with boundary conditions : @xmath126 where @xmath127 is a momentum cutoff regularizing the integral .
we shall eventually take @xmath128 .
we know that the asymptotic behavior of @xmath1 reflecting chiral symmetry breaking is @xmath129 where @xmath130 .
we have found the numerical solution of eqs.([fineq1 ] ) , ( [ f01 ] ) and ( [ fp01 ] ) satisfying this asymptotic behavior . the obtained solution with @xmath128 ( a large enough number which can be regarded as infinity ) in the three cases
are plotted in fig .
2 . again they are different mainly in the low momentum region .    with the obtained @xmath1
, we can calculate the @xmath2 chiral lagrangian coefficients from eqs.([p4full ] ) , ( [ kpi ] ) , and ( [ kresult ] ) .
the obtained values of @xmath15 are listed in table i together with the experimental values @xcite for comparison .
note that there is no running of @xmath15 in this simple approach since the meson - loop effects causing the running of @xmath15 @xcite are of the order of @xmath16 , and are neglected in this approach .
thus the predicted numbers of @xmath15 can be directly compared with the experimental values .
we see from table i that :    \(i ) these coefficients are not so sensitive to the forms of @xmath104 ;    \(ii ) all the obtained @xmath15 are of the right orders of magnitude and the right signs ;    \(iii ) @xmath131 and @xmath132 are consistent with the experiments at the @xmath133 level ;    \(iv ) @xmath134 and @xmath38 are consistent with the experiments at the @xmath135 level ; and    \(v ) only @xmath136 deviates from the experimental value by @xmath137 .    considering the large theoretical uncertainty in this simple approach , _ the obtained @xmath15
are consistent with the experiments_. we see that the nonperturbative quark self - energy plays an important role in qcd contributions to the chiral lagrangian coefficients .
this supports the phenomenological model of holdom @xcite .
in addition to @xmath15 , we can also calculate the quark condensate @xmath138 from the @xmath10 coefficient @xmath139 in eq.([f0b0 ] ) . in the simple approach in this paper ,
the relation between @xmath138 and @xmath139 is @xcite @xmath140 we know that , in this simple approach , @xmath109 mev is finite . but @xmath139 in eq.([f0b0 ] )
is divergent , @xmath141 so that it needs to be renormalized .
we take a simple renormalization scheme by taking the counter term as @xmath142 , in which @xmath115 is the renormalization scale  ) scheme @xcite . ] .
thus the renormalized quantity @xmath143 then the renormalized ( [ condensate ] ) is @xmath144 we take the renormalization scale to be @xmath145 gev to define the quark condensate .
the obtained values of @xmath146 for the three forms of @xmath104 are @xmath147 these are to be compared with the experimentally determined value @xmath148 from the qcd sum rule at the scale of the typical hadronic mass @xcite . considering the large theoretical uncertainty in this calculation , the predicted quark condensate is also consistent with the experiment .
the above results show that _ the present simple approach does reveal the main feature of the qcd predictions for the chiral lagrangian coefficients _ although the approximations in this approach are rather crude . of course , further improvements of the approximations beyond this simple approach are needed .
this kind of study is in progress .    finally , we would like to mention that , in our calculation , we have taken the ultraviolet cutoff parameters @xmath149 , i.e. , we have taken account of the qcd contributions in the whole momentum range .
note that this has nothing to do with the validity range of the chiral lagrangian determined by the range in which the expansion in the meson momentum makes sense , i.e. , up to @xmath150 . to see the role of the qcd contributions from the high momentum region , say above 1 gev
, we have made a check by doing the calculations with the same @xmath1 but taking @xmath151 gev instead of @xmath41 .
the results are listed in table ii . comparing the nonvanishing results in table ii with the corresponding @xmath149 results in table i
, we see that this change of @xmath22 does not cause much difference in @xmath152 and @xmath74 , while it causes @xmath153 and @xmath132 to reduce by at least a factor of 2 .
therefore , we see that @xmath152 and @xmath74 are mainly contributed by the qcd dynamics in the low momentum region , while high momentum region contributions to @xmath154 , and @xmath132 are not negligible .
in this paper , we have calculated the coefficients in the gasser - leutwyler lagrangian from the underlying theory of qcd in a simple approach with the approximations of taking the large-@xmath0 limit , the leading order in dynamical perturbation theory , and the improved ladder approximation based on the qcd formulae given in ref.@xcite to illustrate the main feature of how qcd predicts the chiral lagrangian coefficients . in the calculation
, we use the same regularization technique , the generalized schwinger proper time regularization , in the calculations of the contributions from both the anomaly part and the normal part , so that the relation between the contributions from the two parts can be clearly seen .
we first take the large-@xmath0 limit to evaluate the effective action in qcd .
our first conclusion in this study is that , _ in the large-@xmath0 limit , to all orders in momentum expansion , the anomaly part contributions to the chiral lagrangian coefficients [ cf .
eqs.([anomalyresults ] ) ] given in the literature @xcite from the effective action @xmath155 [ cf .
eq.([sanom ] ) ] are exactly cancelled by the contributions from the piece of the effective action @xmath156 [ cf .
eq.([spi=0 ] ) ] in the normal part contributions _ , so that _ the chiral lagrangian coefficients are eventually contributed by the remaining piece of the normal part effective action @xmath157 [ cf .
eq.([spi])]_. the final qcd expressions for the @xmath2 coefficients are given in eqs.([p4full ] ) .
to simplify @xmath157 , we further make the approximation of taking the leading order in dynamical perturbation theory .
then @xmath157 is reduced to the simple form in eq.([sgl ] ) , and all the chiral lagrangian coefficients are approximately expressed in terms of the quark self - energy @xmath1 shown in eqs.([f0b0])@xmath51([kresult ] ) . to solve the schwinger - dyson equation for @xmath1
, we further take the improved ladder approximation .
lacking of the knowledge about the running coupling constant @xmath104 in the nonperturbative region , we take certain models for it from the literature @xcite [ cf .
eqs.([omega1 ] ) , ( [ omega2 ] ) and ( [ omega3 ] ) ] , and we further take the usual approximation @xmath158 to simplify the calculation .
the quark self - energy reflecting chiral symmetry breaking is obtained by solving the simplified schwinger - dyson equation numerically .
the obtained results of the @xmath2 coefficients are listed in table i. compared with the experimental values of @xmath159 @xcite , _ the agreement of @xmath160 , and @xmath132 is of the level of @xmath133 , and that of @xmath134 and @xmath38 is of the level of @xmath135_. only @xmath136 deviates from the experimental value by @xmath137 . considering the large theoretical uncertainty in this simple approach , _ all the obtained coefficients @xmath159 are consistent with the experiments_. we have also calculated the renormalized quark condensate @xmath146 from the obtained @xmath10 coefficient [ cf .
eq.([condensaterenorm ] ) ] _ which is also consistent with the experiment_.    although the approximations in this simple approach are rather crude , the above results show that _ this simple approach does reveal the main feature of qcd predictions for the chiral lagrangian coefficients_. for studying physics not requiring high precision , this simple approach may already be useful .
of course further improvements of the approximations beyond this simple approach ( reflecting more about qcd dynamics ) are needed .
this kind of study is in progress and will be presented in another paper .
the approach can also be applied to electroweak theories to study how the coefficients in the electroweak chiral lagrangian are predicted by various kinds of underlying gauge theories of the electroweak symmetry breaking mechanism .
this kind of study is also in progress , and will be presented in separate papers .
this work is supported by the national natural science foundation of china , the foundation of fundamental research grant of tsinghua university , and a special grant from the ministry of education of china .
in this appendix , we take the schwinger proper time regulation to regularize the one - loop functional determinant in which the quark self - energy @xmath161 reflecting chiral symmetry breaking takes place .
the matrix element in eq.([lnsigma ] ) can be evaluated in the momentum representation @xmath164}|x\rangle\nonumber\\ & & = \int\frac{d^4p}{(2\pi)^4}\exp\bigg\{-\tau\bigg[\overline{e}(x ) -\nabla_x^2 - 2ip\cdot\nabla_x+p^2 + \sigma^2(-\overline{\nabla}^2 - 2ip\cdot\overline{\nabla}_x+p^2 ) \nonumber\\ & & \hspace{0.5 cm } + \hat{i}_\omega\sigma(-\overline{\nabla}_x^2 - 2ip\cdot\overline{\nabla}_x+p^2 ) + \sigma(-\overline{\nabla}_x^2 - 2ip\cdot\overline{\nabla}_x+p^2)\tilde{i}_\omega -d\!\!\ !
/\;\sigma(-\overline{\nabla}_x^2 - 2ip\cdot\overline{\nabla}_x+p^2 ) \bigg]\bigg\}. % \nonumber\\ % & & \hspace{0.5 cm } % -i[p\!\!\ ! /\;,\sigma(-\overline{\nabla}_x^2 - 2ip\cdot\overline{\nabla}_x+p^2 ) ] % \bigg]\bigg\}. \label{momexp}\end{aligned}\ ] ] then after lengthy but elementary calculations and expanding in powers of the external sources , we can identify the expressions for @xmath165 by comparing with the form of eqs.([p4 ] ) , and the obtained results in the minkowskian space - time are just those given in eqs.([f0 ] ) , ( [ f0b0 ] ) and ( [ kresult ] ) in the text . the details are given in ref.@xcite .    for the evaluation of the effective action @xmath24 in eq.([sanom ] ) in the minkowskian space - time , we note that there is no @xmath166 term in eq.([lnsigma ] ) , but we still have to replace @xmath166 by an infrared cutoff parameter @xmath25 in eq.([lnsigma ] ) to regularize the infrared divergence . then the momentum integration can be explicitly carried out , and we obtain the results in eqs.([anomk ] ) in the text .
here we show that the last term in eq.([spi=0 ] ) , @xmath167 with the saddle point equation [ cf .
eq.([pieq ] ) ] @xmath168^{\rho\sigma}(y , x ) \label{sta1}\end{aligned}\ ] ] is @xmath4-independent .
the @xmath4-rotated quantities in eqs.([3rdterm ] ) and ( [ sta1 ] ) are defined by @xcite @xmath169 ~[j(x)+i\partial\!\!\!\ ! /\;][\omega(x)p_r+\omega^{\dagger}(x)p_l],\nonumber\\ & & \phi_{\omega}^t(x , y)=[\omega^{\dagger}(x)p_r+\omega(x)p_l]~\phi^t(x , y ) [ \omega^{\dagger}(y)p_r+\omega(y)p_l ] , \label{rotation}\end{aligned}\ ] ] and the @xmath170-point green s function @xmath171 is define by @xcite @xmath172\cdots [ \overline{\psi}^{a_n}_{{\alpha}_n}(x_n ) ( \frac{\lambda_{i_n}}{2})_{\alpha_n\beta_n } \gamma^{\mu_n}{\psi}^{a_n}_{\beta_n}(x_n)]\nonumber   \\ & & = { \int}d^{4}x'_1\cdots{d^{4}x'_n } g_s^{n-2}\overline{g}^{\sigma_1\cdots\sigma_n}_{\rho_1 \cdots\rho_n}(x_1,x'_1,\cdots , x_n , x'_n ) \overline{\psi}^{\sigma_1}_{\alpha_1}(x_1){\psi}^{\rho_1 } _ { \alpha_1}(x'_1)\cdots\overline{\psi}^{\sigma_n }   _ { \alpha_n}(x_n){\psi}^{\rho_n}_{\alpha_n}(x'_n ) .
\label{bargdef}\end{aligned}\ ] ]    first we see from eqs.([rotation ] ) that eq.([sta1 ] ) can be written as @xmath173~\phi^t_c(y , x ) [ \omega^{\dagger}(x)p_r+\omega(x)p_l]\bigg ] ^{\rho\sigma}\bigg|_{\pi_c=0}\nonumber\\ & & = -i\bigg[[\omega^{\dagger}(y)p_r+\omega(y)p_l]~ [ ( i\partial\!\!\ !
/\;+j)^{-1}](y , x ) [ \omega^{\dagger}(x)p_r+\omega(x)p_l]\bigg]^{\rho\sigma}\nonumber\\ & & = -i[\omega^{\dagger}(y)p_r+\omega(y)p_l]^{\rho\rho ' } [ ( i\partial\!\!\ ! /\;+j)^{-1}]^{\rho'\sigma'}(y , x ) [ \omega^{\dagger}(x)p_r+\omega(x)p_l]^{\sigma'\sigma}\nonumber\\ & & = -i[\gamma_0v_{\omega}\gamma_0]^{\dagger\rho\rho'}(y ) [ ( i\partial\!\!\ ! /\;+j)^{-1}]^{\rho'\sigma'}(y , x ) v_{\omega}^{\sigma'\sigma}(x ) , \label{phio}\end{aligned}\ ] ] in which @xmath174 satisfies @xmath175 and @xmath176(y , x ) \label{phic}\end{aligned}\ ] ] is @xmath4-independent .
next , we look at this transformed green s function @xmath181 . from the definition ( [ bargdef ] ) and the property ( [ vcondition ] ) we have @xmath182\cdots\nonumber\\ & & \cdots[(\overline{\psi}\gamma_0v^{\dagger}_\omega\gamma_0)^{a_n}_{{\alpha}_n}(x_n ) ( \frac{\lambda_{i_n}}{2})_{\alpha_n\beta_n } \gamma^{\mu_n}(v_\omega{\psi})^{a_n}_{\beta_n}(x_n)]\nonumber   \\   & = & g_{\mu_1\cdots\mu_n}^{i_1\cdots{i_n } }   ( x_1,\cdots , x_n)[\overline{\psi}^{a_1}_{{\alpha}_1}(x_1 )   ( \frac{\lambda_{i_1}}{2})_{\alpha_1\beta_1}\gamma^{\mu_1 }   { \psi}^{a_1}_{{\beta}_1}(x_1)]\cdots [ \overline{\psi}^{a_n}_{{\alpha}_n}(x_n ) ( \frac{\lambda_{i_n}}{2})_{\alpha_n\beta_n } \gamma^{\mu_n}{\psi}^{a_n}_{\beta_n}(x_n)]\nonumber   \\ & = & { \int}d^{4}x'_1\cdots{d^{4}x'_n } g_s^{n-2}\bar{g}^{\sigma_1\cdots\sigma_n}_{\rho_1 \cdots\rho_n}(x_1,x'_1,\cdots , x_n , x'_n ) \overline{\psi}^{\sigma_1}_{\alpha_1}(x_1){\psi}^{\rho_1 } _ { \alpha_1}(x'_1)\cdots\overline{\psi}^{\sigma_n }   _ { \alpha_n}(x_n){\psi}^{\rho_n}_{\alpha_n}(x'_n),\nonumber\end{aligned}\ ] ] i.e. , @xmath183 thus the transformed green s function @xmath184 in eq.([sgv ] ) can be replaced by @xmath185 , and eq.([sgv ] ) becomes @xmath186 which is independent of @xmath4 .
1[biblio ] s. weinberg , physica , * 96a * , 327 ( 1979 ) . j. gasser and h. leutwyler , ann
( n.y . ) * 158 * , 142 ( 1984 ) ; nucl . phys .
* b 250 * , 465 ( 1985 ) .
q. wang , y .- p .
kuang , x .-
wang , and m. xiao , phys .
d * 61 * , 054011 ( 2000 ) , hep - ph9903201 .
d. espriu , e. de rafael , and j. taron , nucl . phys .
* b 345 * , 22 ( 1990 ) ; * b 355 * , 278(e ) ( 1991 ) .
q. lu , h. yang , and q. wang , tsinghua university preprint tuhep - th-02133 , hep - ph/0202078 , to appear in commun .
h. pagels and s. stokar , phys .
d * 20 * , 2947 ( 1979 ) .
simic , phys .
d * 34 * , 1903 ( 1986 ) .
b. holdom , phys . rev .
d * 45 * , 2534 ( 1992 ) ; phys .
b * 294 * , 293 ( 1992 ) ; b. holdom , j. terning , and k. verbeek , _
ibid_. , * 245 * , 612 ( 1990 ) ; _
ibid_. , * 273 * , 549 ( e ) ( 1991 ) ; b. holdom and q. wang , phys .
d * 55 * , 3101 ( 1997 ) .
aoki , m. bando , t. kugo , m.g .
mitchard , and h. nakatani , prog .
. phys . * 84 * , 683(1990 ) .
p. jain and h.j .
munczek , phys .
d * 44 * , 1873(1991 ) .
b . dai , c .- s .
huang , and d .- s .
liu , phys .
d * 43 * , 1717 ( 1991 ) ; y .- b .
dai , y .- b .
ding , c .- s .
huang , c .-
wang , and y .- l zhu , phys . rev .
d * 47 * , 1256 ( 1993 ) . s .-
h. tye , e. tomboulis , and e. c. poggio , phys .
d * 11 * , 2839 ( 1975 ) ; v. a. miransky , sov .
. phys . * 38 * , 280 ( 1983 ) ; k. higashijima , phys .
d * 29 * , 1228 ( 1984 ) .
for instance , h .- j .
he and y .- p .
kuang , z. phys .
c * 47 * , 565 ( 1990 ) .
shifman , a.i .
vainshtein , and v.i .
zakharov , nucl .
* b 147 * , 385 ( 1979 ) .
roberts , r.t .
cahill and j. praschifka , ann .
phys . * 188 * , 20 ( 1988 ) .
e. witten , nucl . phys . * b 223 * , 422 ( 1983 ) .      .
the obtained values of the @xmath2 coefficients @xmath159 for model a [ eq.([omega1 ] ) ] , model b [ eq.([omega2 ] ) ] and model c [ eq.([omega3 ] ) ] with @xmath149 together with the experimental values [ eq.([exptl ] ) ] for comparison .
@xmath106 is in mev , and the coefficients are in units of @xmath34 .
the obtained @xmath1 from the schwinger - dyson equation ( [ eq0 ] ) with model a [ eq.([omega1 ] ) ] , model b [ eq.([omega2 ] ) ] and model c [ eq .
( [ omega3 ] ) ] for the running coupling constant @xmath40 .
the solid , dashed , and dotted lines are for models a , b , and c , respectively .


the standard model of particle physics is known to have many problems with observations , such as the matter abundance of our universe , and many candidates of new physics have been invented so far . among them ,
the supersymmetry @xcite is widely noticed , due to the possibility to resolve many phenomenological problems encountered in the standard model .
one of the most important experimental observables to probe the @xmath2 violation in the supersymmetric models is the electric dipole moment ( edm ) .
the edm is sensitive to the violation of parity and time reversal symmetry ( or equivalently the @xmath2 symmetry ) , and is known to have very small standard model prediction @xcite .
this fact gives a direct evidence of the @xmath2 violation beyond the standard model if any finite value is observed experimentally .
due to the stringent experimental limits on the edm measured in many systems @xcite , many @xmath2 phases of the supersymmetric models with @xcite and without @xmath0-parity @xcite were constrained so far .    in this work
we focus on the @xmath2 violation of the @xmath0-parity violating ( rpv ) sector @xcite .
the conservation of the @xmath0-parity is often assumed to keep the lightest supersymmetric particle stable , but this is only an _ ad hoc _ assumption , and phenomenologically , the existence of the rpv interactions is still allowed .
the ( trilinear ) rpv interactions contribute to the fermion edm from the two - loop level @xcite , and the barr - zee type process ( see fig . [ fig : barr - zee ] ) is currently known to give the dominant effect @xcite .
barr - zee type diagrams within @xmath0-parity violation .
, width=302 ]    recently , it was pointed that the two - loop level _ rainbowlike _ contribution to the fermion edm may have comparable size with the leading barr - zee type process in the @xmath0-parity conserving sector @xcite .
this contribution has a fermion - sfermion inner loop , connected to the fermion external line with a higgsino and a gaugino , and the flavor structure is exactly the same as for the barr - zee type diagram , with the same coupling constants . in the rpv sector ,
similar rainbowlike diagrams can be drawn by replacing the higgsino with the neutrino ( see fig . [
fig : rpv_rainbow ] ) , and it is therefore of necessity to evaluate and analyze this potentially sizable contribution .
example of rainbow diagram contributing to the fermion edm within @xmath0-parity violation .
the neutralino is given by @xmath3 and the gluino by @xmath4 .
, width=264 ]    this paper is organized as follows .
we first introduce the supersymmetric and rpv interactions relevant in this discussion in the next section .
we then calculate in section [ sec : rainbow ] the rainbowlike contribution to the fermion edm and to the chromo - edm within @xmath0-parity violation and give the explicit formulae . in section [ sec : analysis ] , we analyze the effect of the rpv supersymmetric rainbow diagrams by comparing with the edm and the chromo - edm generated by the rpv barr - zee type diagrams .
the last section is devoted to the summary .
we have also collected the detailed derivation of important formulae in the appendix .
let us first give the lagrangian of the particles relevant in this discussion .
the rpv interactions relevant in this discussion are given by the following superpotential : @xmath5 with @xmath6 indicating the generation , @xmath7 the @xmath8 indices .
the @xmath9 indices have been omitted .
@xmath10 and @xmath11 denote the lepton doublet and singlet left - chiral superfields .
@xmath12 , @xmath13 and @xmath14 denote respectively the quark doublet , the up - type quark singlet and the down - type quark singlet left - chiral superfields .
the bilinear term has been omitted in our discussion .
also baryon number violating rpv interactions ( @xmath15 ) will be omitted from now , to avoid rapid proton decay .
this rpv superpotential gives the following lepton number violating yukawa interactions : @xmath16 + ( { \rm h.c.})\nonumber\\ & & -\lambda'_{ijk } \left [ \tilde \nu_i \bar d_k p_l d_j + \tilde d_{lj } \bar d_k p_l \nu_i + \tilde d_{rk}^\dagger \bar \nu_i^c p_l d_j \right . \nonumber\\ & & \hspace{3.5em } \left .
-\tilde e_{li } \bar d_k p_l u_j - \tilde u_{lj } \bar d_k p_l e_i - \tilde d_{rk}^\dagger \bar e_i^c p_l u_j \right ] \nonumber\\ & & \hspace{4em }   + ( { \rm h.c . } )   \ .\end{aligned}\ ] ] the soft breaking terms of the rpv sector will not be considered in this discussion
.    we should also give the supersymmetric interactions of the @xmath0-parity conserving sector .
the neutralino mass matrix is given as follows : @xmath17 with @xmath18 where the corresponding neutralino field vector @xmath3 has as the first two components the up- and down - type higgsinos , and the last two components refer to the @xmath19 and the @xmath8 gauginos . in this work , we assume that the supersymmetry breaking scale is beyond tev , which is strongly suggested by the recent result of the lhc experiment @xcite , and we neglect the mixing between gauginos and higgsinos . in our discussion ,
the higgsino is thus irrelevant .
the components of the mass matrix have @xmath2 phases in the general case , and they can manifest themselves as observable effect when the mass insertion occurs in the process .
for the gluino , we remove its @xmath2 phase as usual .
the sfermion mass matrix is given by @xmath20 where @xmath21 for up - type squarks , and @xmath22 for down - type squarks and charged sleptons .
note that in this discussion , only the down - type squarks and the charged sleptons are relevant . to obtain the sfermion mass eigenbasis , we introduce the following unitary matrix with the sfermion mixing angle @xmath23 and the @xmath2 phase @xmath24 : @xmath25 with @xmath26 and @xmath27 for @xmath28 . in this approximation ,
the two mass eigenvalues are given by @xmath29 , and @xmath30 .
the @xmath2 phase is given by @xmath31 .
we should also give the gaugino - sfermion - fermion interactions .
they are given by @xmath32 where @xmath33 is the gauge index of the fermion @xmath34 ( but will be omitted from now ) .
the convention for the sign of the gauge coupling is @xmath35 where @xmath36 is the generator of the gauge group .
the index @xmath37 gives the gauge group of the gaugino @xmath38 .
we define @xmath39 as the @xmath19 gaugino , @xmath40 as the neutral @xmath8 gaugino and @xmath41 as the gluino .
the fermion - sfermion - gaugino coupling constants are given as follows : for the gluino couplings , we have @xmath42 and @xmath43 . for the couplings of @xmath39 , we have @xmath44 , @xmath45 , @xmath46 and @xmath47 .
we have finally for the @xmath40 couplings @xmath48 , and @xmath49 .
let us now calculate the rpv rainbow diagrams . we can classify the rainbow diagrams into two types .
the first type is the insertion of the effective one - loop level neutrino - gaugino - gauge boson vertex [ fig .
[ fig : rainbow_type ] ( 1 ) ] .
the second type is the insertion of the effective one - loop level neutrino - gaugino transition [ fig .
[ fig : rainbow_type ] ( 2 ) ] .
possible insertions of one - loop effective vertex for the rainbow diagram in rpv supersymmetry .
the grey blob represents the one - loop effective vertex .
the solid lines attached to the grey blob indicate the neutrino and gaugino .
, width=302 ]    the derivation of the rpv rainbow diagrams is very similar to that of the rainbow diagrams in the minimal supersymmetric standard model @xcite .
the detail of the computation of the type 1 and 2 insertions is presented in appendices [ sec : type1rainbow ] and [ sec : type2rainbow ] , respectively .    to obtain the edm contribution ,
the above rainbow diagrams are calculated to the first order in the external momentum carried by the gauge boson ( this manipulation is suitable since the edm is the first order coefficient of the multipolar expansion ) . in our calculation , we have neglected the fermion mass originating from the mass insertion of the inner loop .
the formulae of the edm of the fermion @xmath50 ( @xmath51 is the flavor of the external fermion ) given by the type 1 [ fig . [ fig : rainbow_type ] ( 1 ) ] and the type 2 [ fig .
[ fig : rainbow_type ] ( 2 ) ] diagrams are respectively    @xmath52 \ , , \label{eq : edmnugg } \\ d^2_{f_k } & = & -\sum_{n=1,2 }   { \rm i m } ( \hat \lambda_{ijj } \tilde \lambda^*_{ikk } e^{i(\theta_n -\delta_{f_j } ) } ) \frac{q_f n_c } { 64 \pi^3 } |m_{\lambda_n } | \frac{e ( g^{(n)}_{\tilde f_l } + g^{(n)}_{\tilde f_r } ) } { 4\pi } \sin \theta_{f_j } \cos \theta_{f_j }   \nonumber\\ & & \hspace{10em } \times \sum_{\tilde f= \tilde f_l , \tilde f_r } g^{(n)}_{\tilde f } m_{\tilde f_{k}}^2 \left [ f''(|m_{\lambda_n}|^2 , 0 , m_{\tilde f_{k}}^2 , m_{\tilde f_{1j}}^2 \ , ) -f''(|m_{\lambda_n}|^2 , 0 , m_{\tilde f_{k}}^2 , m_{\tilde f_{2j}}^2 \ , ) \right ] \
, , \label{eq : edmnug}\end{aligned}\ ] ]    where @xmath53 for inner quark - squark loop and @xmath54 for inner lepton - slepton loop . here
@xmath55 and @xmath56 are the flavor indices of the intermediate neutrino and the fermion - sfermion of the inner loop , respectively .
the rpv couplings are given by @xmath57 for the charged lepton - slepton inner loop , @xmath58 in the case of the quark - squark inner loop , @xmath59 for the external lepton , and @xmath60 for the external quark .
the constant @xmath61 is @xmath62 for the left - handed sfermion @xmath63 and @xmath64 for the right - handed sfermion @xmath65 .
the electric charge of the inner loop fermion ( sfermion ) and the external fermion in the unit of @xmath66 is denoted , respectively , by @xmath67 and @xmath68 .
the functions @xmath69 and @xmath70 are defined by @xmath71 -\frac{f}{ac } \frac{\pi^2}{6 } \
, , \label{eq : f ' } \\ % % % % % % % % % % % % % % f '' ( a,0,c , f ) % % % % % % % % % % % % % % f''(a,0,c , f ) & \equiv & \frac{(a - f)^2 } { a(a - c)^3 } \left [ { \rm li}_2 ( 1-a / f ) -{\rm li}_2 ( 1-c / f ) \right ]   - \frac{f^2}{ac^3 } \left [ { \rm li}_2 ( 1-c / f ) - \frac{\pi^2}{6 } \right ] \nonumber\\ & & + \frac{1}{a(c - f ) } \left [ \frac{(a - f)^2}{(a - c)^2 } -\frac{f^2}{c^2 } \right ] \ln \frac{c}{f } - \frac{c \ln c -   c+2f } { 2c^2 ( a - c ) }   + \frac{f+a \ln f}{(a - c)^2 } \left [ \frac{1}{c } -\frac{1}{a - c } \ln \frac{a}{c } \right ] \ , , \label{eq : f''}\end{aligned}\ ] ] with @xmath72 denoting the dilogarithm function .
the original functions @xmath73 and @xmath74 were defined in ref .
similarly , the chromo - edm is given as follows : @xmath75 \ , .
\label{eq : cedm}\ ] ]    note that the gluino mass parameter @xmath76 can be negative , and the overall sign of @xmath77 may change .
we see that the formulae ( [ eq : edmnugg ] ) , ( [ eq : edmnug ] ) , and ( [ eq : cedm ] ) are very similar to the supersymmetric rainbow contribution with @xmath0-parity conservation @xcite ( this is just a replacement of the matter - higgs yukawa couplings by the rpv couplings , and the higgsino mass @xmath78 by the negligible mass of the neutrino ) .
the rpv rainbow diagrams vanish if the mass eigenvalues of the inner loop sfermion @xmath79 and @xmath80 are degenerate .
this property can also be seen for the supersymmetric rainbow process with @xmath0-parity conservation .
we now analyze the contribution of the rpv rainbow diagrams to the edm and to the chromo - edm of the fermions . in this discussion , we only consider processes with the third generation fermion - sfermion inner loop , since they give the dominant effect due to @xmath81 .
we choose the following values for the sparticle masses in the function of the supersymmetry breaking scale @xmath82 : @xmath83 , @xmath84 , @xmath85 , @xmath86 , and @xmath87 .
the down squark masses for both gauge eigenstates are given by @xmath88 , and @xmath89 for the selectron .
the @xmath19 and @xmath8 gaugino mass is given by @xmath90 , and the sneutrino mass by @xmath91 . in our discussion , we fix the @xmath92 parameter to @xmath93 tev . for @xmath1 , we assign @xmath94 .
the large value of @xmath1 is important , since the rainbow diagram contribution will be enhanced for large @xmath1 , due to the factor @xmath95 .
the rainbow diagram will be compared with the leading rpv barr - zee type diagram @xcite .
the fermion edm and the quark chromo - edm generated from the barr - zee type diagram is given respectively as : @xmath96    the first case to consider is the rainbow diagram with external down - quark and quark - squark inner loop .
in this case , the @xmath2 violating combinations of parameters im(@xmath97 ) , im(@xmath98 ) , and im(@xmath99 ) ( @xmath100 ) contribute to the neutron edm via the edm and the chromo - edm .
the dependence of the quark edm and the chromo - edm on the neutron edm is given by the following relation , obtained from the analysis using qcd sum rules @xcite : @xmath101 where we have assumed the peccei - quinn symmetry @xcite .
the rainbow diagram effect ( for both edm and chromo - edm ) with external down - quark and bottom quark - bottom squark inner loop , together with that of rpv barr - zee diagram , is plotted in fig .
[ fig : rpv_q - qrbnedm ] .
contribution of the rpv rainbow diagram with external down quark and @xmath102 inner loop ( absolute value ) to the neutron edm plotted in the function of the supersymmetry breaking scale [ we have set @xmath103 .
the barr - zee type contribution is also plotted for comparison .
, width=302 ]    we see that the chromo - edm generated by the rpv barr - zee diagram dominates over all other contributions .
the dominance of the chromo - edm over the quark edm can be explained by the large value of the strong coupling @xmath104 over the electromagnetic coupling @xmath105 .
we remark also that the rainbow and barr - zee type contributions have very close size near @xmath106 tev for the quark edm , even though the two effects have a difference by about one order of magnitude for the chromo - edm .
this is due to the suppression of the barr - zee type diagram of the quark edm , due to the double factor of the fractional charge of the down - type quarks ( @xmath107 ) .
we can see also the fast decrease of the rainbow diagrams against the barr - zee type effect in high energy region of @xmath82 .
this can be explained by simple dimensional analysis of the scaling in @xmath82 . from eqs .
( [ eq : barr - zee ] ) and ( [ eq : chromo - barr - zee ] ) , we see that the rpv barr - zee type diagram is scaled as @xmath108 .
the rpv rainbow diagram is however scaled as @xmath109 , due to @xmath110 , where we have fixed @xmath92 .
if the @xmath92-parameter have a supersymmetric dynamical origin @xcite and the @xmath92-parameter is scaled as @xmath82 , the rainbow diagram contribution will be scaled as @xmath108 .
we must note that the rainbow diagram and the barr - zee type diagram can interfere constructively or destructively , depending on the phases @xmath111 , @xmath112 and the signs of the gluino mass parameter @xmath76 and of @xmath113 .
the rpv barr - zee type contribution to the quark chromo - edm however dominates over the rainbow chromo - edm with external down quark and @xmath114 inner loop by one order of magnitude , so it is difficult to enhance or cancel the total contribution , and the result will not be much changed from the analysis taking into account the barr - zee type diagrams alone . in fig .
[ fig : rpv_q - qrbnedm ] , we can see that the rpv contribution to the neutron edm ( plotted with @xmath115 ) largely exceed the currently known experimental data of the neutron edm ( @xmath116 ) @xcite .
it is thus possible to constrain this combination of rpv couplings .
if @xmath117 and the interference is constructive , we obtain the following limit @xmath118 for @xmath106 tev .     contribution of the rpv rainbow diagram with external electron and @xmath119 inner loop ( absolute value ) to the electron edm plotted in the function of the supersymmetry breaking scale [ we have set @xmath120 .
the barr - zee type contribution is also plotted for comparison .
, width=302 ]    the next case to consider is the rainbow diagram with external electron and charged lepton - slepton inner loop .
in this case , the @xmath2 violating combinations of parameters im(@xmath121 ) and im(@xmath122 ) contribute to the electron edm .
the rainbow diagram effect with external electron and @xmath123-@xmath124 inner loop , together with that of the rpv barr - zee type diagram , is plotted in fig .
[ fig : rpv_l - lrbeedm ] . here
we find that the rainbow and the barr - zee type contributions have similar size .
this means that important enhancement or cancellation may occur , depending on the phase of @xmath125 and the sign of @xmath126 . in this case
we do not have the suppression of the barr - zee type process by fractional charge .
the comparable size of these two contributions is natural , since both involve no colored sparticles which were assumed to be heavier than uncolored sparticles .
here we find again that the rainbow diagram decreases faster than the barr - zee type diagram in heavier supersymmetry breaking region .
the reason is the same as for the down - quark edm and chromo - edm with quark - squark inner loop . as we can see in fig .
[ fig : rpv_l - lrbeedm ] , this rpv contribution to the electron edm ( for @xmath127 ) exceeds largely the limit given by the experimental limit of the electron edm ( @xmath128 ) @xcite .
if the interference is constructive , the limit of the rpv couplings @xmath129 is improved by about factor 1.6 from the analysis of the barr - zee type diagram alone ( for @xmath106 tev ) : @xmath130 if the interference is destructive , the limit to the rpv couplings will of course be loosened .
contribution of the rpv rainbow diagram with external electron and @xmath102 inner loop ( absolute value ) to the electron edm plotted in the function of the supersymmetry breaking scale [ we have set @xmath131 .
the barr - zee type contribution is also plotted for comparison .
, width=302 ]    the third case to consider is the rainbow diagram with external electron and quark - squark inner loop . in this case , the @xmath2 violating combinations of parameters im(@xmath132 ) and im(@xmath133 ) ( @xmath134 ) contribute to the electron edm .
the rainbow diagram effect with external electron and bottom quark - squark inner loop is plotted in fig .
[ fig : rpv_l - qrbeedm ] , with also the barr - zee type contribution for comparison . here
again we find that the rainbow and barr - zee type contributions have similar size .
this is due to the suppression of the barr - zee type diagram with a factor of the fractional charge of the bottom quark . as for
the electron edm with lepton - slepton loop , the interference can enhance or suppress the total contribution . in this case
, we can constrain the rpv couplings im(@xmath133 ) with the experimental data of the electron edm . if the interference is constructive and maximal
, we obtain the following limit to the rpv couplings for @xmath135 tev : @xmath136 this limit is tighter than the constraint given by the experimental data of the @xmath137hg atom @xcite via the @xmath138 , @xmath2-odd electron - nucleon interaction @xcite @xmath139 this result is important since , the barr - zee type contribution alone could not give the tightest constraint .
if the interference is destructive , the limit to the rpv couplings is loosened , and may even disappear .
contribution of the rpv rainbow diagram with external down quark and @xmath140 inner loop ( absolute value ) to the neutron edm plotted in the function of the supersymmetry breaking scale [ we have set @xmath141 . the barr - zee type contribution is also plotted for comparison . , width=302 ]    the final case to consider is the rainbow diagram with external down - quark and charged lepton - slepton inner loop . in this case , im(@xmath142 ) and im(@xmath143 ) ( @xmath134 ) contribute to the neutron edm . by again using the relation ( [ eq : qcdsumrule ] ) ,
we obtain the plot of fig .
[ fig : rpv_q - lrbnedm ] .
we see that although being comparable , the rpv rainbow diagram is smaller than the barr - zee type contribution .
this is because the effect of the type 2 rainbow diagram @xmath144 is suppressed by the fractional charge of the quark [ see eq . ( [ eq : edmnug ] ) ] .
as for the other cases , if the interference is constructive , the limit on the rpv couplings is tightened as @xmath145 for @xmath106 tev . on the contrary , the limit will be loosened for the destructive interference .
in this paper , we have calculated a new sizable rpv contribution to the fermion edm , the rainbow diagram .
the calculation was similar to that of the rainbow diagram of the minimal supersymmetric standard model , and we have given the formulae for it .
it was found that the rpv rainbow diagrams have sizable effects compared with the currently known rpv barr - zee type diagrams , for the electron edm with @xmath102 or @xmath119 inner loops , and for the neutron edm with @xmath119 loop .
we have also given new limits on the related combinations of rpv couplings when the interference is maximally constructive .
the rpv rainbow diagrams have also many interesting properties , such as the large enhancement due to large @xmath1 , the cancellation for the equal mass eigenvalues of the inner loop sfermions , the @xmath109 scaling , etc .
the interference between the barr - zee type and rainbow diagrams is dependent on the phases @xmath146 , @xmath147 and the signs of the gluino mass parameter and of @xmath148 .    in this work
, we did not considered the rpv rainbow diagrams with chargino propagation .
as the matter - chargino interactions have comparable coupling with the matter - neutralino ones , it is highly probable that the rainbow diagrams with chargino have sizable effects .
we should study this topic in the next work .
the author thanks h. kamano for useful discussion and comments .
in this section , we calculate the type 1 rainbow diagram with the insertion of the one - loop effective neutrino - gaugino - gauge boson vertex [ see fig . [
fig : rainbow_type ] ( 1 ) ] .
the inner one - loop diagrams are depicted in fig .
[ fig : nugg ] .     one - loop contribution to the effective neutrino - gaugino - gauge boson vertex .
, width=680 ]    this contribution can generate both edm and chromo - edm .
the amplitude of the effective neutrino - gaugino - gauge boson vertex is given by the sum of the diagrams shown in fig .
[ fig : nugg ] : @xmath149 the off - shell amplitudes @xmath150 , @xmath151 , @xmath152 , and @xmath153 are given as : @xmath154\nonumber\\ & & \hspace{1.5em}+ 2g^{(n)}_{\tilde f_r } \bar \lambda_n ( k\hspace{-.5em}/\ , + q\hspace{-.45em}/_2 ) p_l \nu_i \cdot k^\mu \sin \theta_{f_j } \cos \theta_{f_j } e^{i\frac{\theta_n}{2 } - i\delta_{f_j } } \left [ d ( m_{\tilde f_{2j}}^2 , m_{\tilde f_{2j}}^2 , m _ { f_j}^2 ) - d ( m_{\tilde f_{1j}}^2 , m_{\tilde f_{1j}}^2 , m
_ { f_j}^2 ) \right ] \biggl\ } \ , \label{eq : m(a)}\end{aligned}\ ] ] @xmath155\nonumber\\ & & \hspace{1.5em } -2g^{(n)}_{\tilde f_l } \bar \lambda_n ( k\hspace{-.5em}/\ , + q\hspace{-.5em}/_2 ) p_l \nu_i \cdot k^\mu \sin \theta_{f_j } \cos \theta_{f_j } e^{i\frac{\theta_n}{2 } - i\delta_{f_j } } \left [ d ( m_{\tilde f_{2j}}^2 , m_{\tilde f_{2j}}^2 , m _ { f_j}^2 ) - d ( m_{\tilde f_{1j}}^2 , m_{\tilde f_{1j}}^2 , m _ { f_j}^2 ) \right ] \biggl\ } \ , \label{eq : m(a')}\end{aligned}\ ] ] @xmath156 \nonumber\\ & & \hspace{1.5em } - g^{(n)}_{\tilde f_l } \bar \lambda_n ( k\hspace{-.5em}/\ , \gamma^\mu ( k\hspace{-.5em}/\ , - q\hspace{-.5em}/_1 ) + m_{f_j}^2 \gamma^\mu ) p_l \nu_i \cos \theta_{f_j } \sin \theta_{f_j }   e^{i\frac{\theta_n}{2 } - i\delta_{f_j } } \left [ d ( m _ { f_j}^2 , m _ { f_j}^2 , m_{\tilde f_{2j}}^2 ) - d ( m _ { f_j}^2 , m _ { f_j}^2 , m_{\tilde f_{1j}}^2   ) \right ] \biggr\ } \ , , \nonumber\\ \label{eq : m(b ) } \\ % % % % % % % % % % % % % % % %   m_(b ' ) % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % i{\cal m}_{\rm ( b ' ) } & = & \hat \lambda_{ijj } \sqrt{2 } g_f t_n \epsilon_\mu^ * ( q_1 ) \int \frac{d^4 k}{(2\pi ) ^4 } \nonumber\\ & & \times \biggl\ { -g^{(n)}_{\tilde f_l } m_{f_j } \bar \lambda_n ( 2k^\mu - \gamma^\mu q\hspace{-.5em}/_1 ) p_l \nu_i   \ , e^{-i\frac{\theta_n}{2 } } \left [ \cos^2 \theta_{f_j } d ( m _ { f_j}^2 , m _ { f_j}^2 , m_{\tilde f_{1j}}^2 ) + \sin^2 \theta_{f_j }
d ( m _ { f_j}^2 , m _ { f_j}^2 , m_{\tilde f_{2j}}^2   ) \right ] \nonumber\\ & & \hspace{1.5em}+   g^{(n)}_{\tilde f_r } \bar \lambda_n ( k\hspace{-.5em}/\ , \gamma^\mu ( k\hspace{-.5em}/\ , - q\hspace{-.5em}/_1 ) + m_{f_j}^2 \gamma^\mu ) p_l \nu_i \cos \theta_{f_j } \sin \theta_{f_j } e^{i\frac{\theta_n}{2 } - i\delta_{f_j } }    \left [ d ( m _ { f_j}^2 , m _ { f_j}^2 , m_{\tilde f_{2j}}^2 ) - d ( m _ { f_j}^2 , m _ { f_j}^2 , m_{\tilde f_{1j}}^2   ) \right ] \biggr\ } \ , , \nonumber\\ \label{eq : m(b')}\end{aligned}\ ] ] where @xmath157 is the electromagnetic gauge coupling of the fermion @xmath34 or sfermion @xmath158 of the inner loop ( the strong gauge coupling @xmath159 for the chromo - edm ) , @xmath160 the gauge coupling to the fermion - sfermion pair [ see eq . ( [ eq : gaugino - int ] ) for the lagrangian ] . here
@xmath161 is the rpv coupling , with @xmath162 when charged lepton - slepton runs in the loop , and @xmath58 for the case of down - type quark - squark loop .
the factor @xmath163 is given by @xmath164 for the contribution with external photon and @xmath165 for external gluon .
the momenta carried by the exiting gauge boson and gaugino are denoted by @xmath166 and @xmath167 , respectively .
the gaugino and the neutrino spinors are given by @xmath38 and @xmath168 , respectively .
note that the gaugino @xmath38 is majorana .
the polarization vector of the external gauge boson is given by @xmath169 .
here we have used the following notation to simplify the above equations : @xmath170 \left [ ( k - q_1)^2 - y \right ] \left [ ( k+q_2)^2 - z \right ] } \ .\ ] ] terms with @xmath171 and @xmath172 are the contribution from the mass insertion in the fermion propagator , and terms with @xmath173 are given by the mass insertion in the sfermion propagator .    by integrating eqs .
( [ eq : m(a ) ] ) , ( [ eq : m(a ) ] ) , ( [ eq : m(b ) ] ) , and ( [ eq : m(b ) ] ) , we obtain the following one - loop level effective neutrino - gaugino - gauge boson vertex ( expanded up to the first order in the external momentum carried by the gauge boson ) : @xmath174 in writing the above off - shell amplitude , we have omitted several other dirac bilinears with different lorentz structure , such as @xmath175 , @xmath176 .
these terms do not contribute to the edm operator @xmath177 after insertion into the second loop ( they cancel with their complex conjugates ) , so we do not consider them further .
we have also neglected the quark and lepton masses .
similarly , the off - shell effective antineutrino - gaugino - gauge boson amplitude is given by @xmath178 this amplitude is given by the complex conjugated rpv coupling @xmath179 .
note that we also have to consider the amplitudes generated by transposed fermion interactions .
they are given by @xmath180 here the momentum @xmath167 is of course that carried by the exiting neutrino or antineutrino .
we should note that the one - loop level effective vertices ( [ eq : nulamg ] ) , ( [ eq : antinulamg ] ) , ( [ eq : lamnug ] ) , and ( [ eq : lamantinug ] ) are all first order in @xmath166 .
we now insert the above effective vertices ( [ eq : nulamg ] ) , ( [ eq : antinulamg ] ) , ( [ eq : lamnug ] ) , and ( [ eq : lamantinug ] ) to the second loop , as shown in fig .
[ fig : rainbow_type ] ( 1 ) .
we obtain then the following contribution for the type 1 insertion ( expanded in the first order of @xmath166 ) : @xmath181 \left [ k^2 - m_{\tilde f_k}^2 \right ] } \int_0 ^ 1 dz \,\frac{z } { z k^2 - m_{\tilde f_{2j}}^2 } - ( m_{\tilde f_{2j}}^2 \leftrightarrow m_{\tilde f_{1j}}^2 ) \nonumber\\ & \approx & \frac{i}{64 \pi^3 } { \rm i m } ( \hat \lambda_{ijj } \tilde \lambda^*_{ikk } e^{i(\theta_n -\delta_{f_j } ) } )    \nonumber\\ & & \hspace{2em } \times \epsilon_\mu^ * ( q_1 ) \bar f_k \sigma^{\mu \nu } ( q_1)_\nu \gamma_5 f_k \sum_{\tilde f= \tilde f_l , \tilde f_r } s\ , g^{(n)}_{\tilde f }   \left [ f'(|m_{\lambda_n}|^2 , 0 , m_{\tilde f_{k}}^2 , m_{\tilde f_{1j}}^2 \ , ) -f'(|m_{\lambda_n}|^2 , 0 , m_{\tilde f_{k}}^2 , m_{\tilde f_{2j}}^2 \ , ) \right ] \ , , \end{aligned}\ ] ] where we have omitted terms which do not contribute to the edm operator @xmath177 .
the rpv coupling @xmath182 is given by @xmath59 for external lepton and @xmath183 for external quark .
the constant @xmath61 is @xmath62 for left - handed sfermion @xmath184 and @xmath64 for right - handed sfermion @xmath185 .
the functions @xmath69 and @xmath70 are defined by eqs .
( [ eq : f ] ) and ( [ eq : f  ] ) , respectively . from the above equation ,
we obtain the formulae for the edm ( [ eq : edmnugg ] ) and for the chromo - edm ( [ eq : cedm ] ) .
in this section , we calculate the second type rainbow diagram with the insertion of the one - loop effective neutrino - gaugino transition ( see fig . [
fig : rainbow_type ] ( 2 ) ) .
the contributing one - loop diagrams are shown in fig . [
fig : nug ] .
the one - loop level effective neutrino - gaugino transition within @xmath0-parity violation is given by @xmath186 \left [ ( k+q_2)^2 - m_{f_j}^2 \right ] } + \frac{g^{(n)}_{\tilde f_l } \sin^2 \theta_{f_j } + g^{(n)}_{\tilde f_r } \cos^2 \theta_{f_j } } { \left [ k^2 - m_{\tilde f_{2i}}^2 \right ] \left [ ( k+q_2)^2 - m_{f_j}^2 \right ] } \right\ } \nonumber\\ & & - \sqrt{2 } \hat \lambda_{ijj } n_c ( g^{(n)}_{\tilde f_l } + g^{(n)}_{\tilde f_r } ) \sin \theta_{f_j } \cos \theta_{f_j } e^{i\frac{\theta_n}{2 } - i\delta_{f_j } } \int \frac{d^4 k}{(2\pi ) ^4 } \left\ { \frac{\bar \lambda_n ( k\hspace{-.5em}/\ , + q\hspace{-.5em}/_2 ) p_l \nu_i } { \left [ k^2 - m_{\tilde f_{2i}}^2 \right ] \left [ ( k+q_2)^2 - m_{f_j}^2 \right ] } - ( m_{\tilde f_{2i}}^2 \leftrightarrow m_{\tilde f_{1i}}^2 ) \right\ } , \nonumber\\\end{aligned}\ ] ] where @xmath167 is the momentum carried by the incident fermion .
the other notations are the same as for appendix [ sec : type1rainbow ] .
it is important to note that the first line of the above equation is divergent .
this divergence occurs because the first term is a one - loop correction to the tree - level rpv higgs - slepton transition ( note the mass insertion @xmath187 ) . as we have assumed that the bilinear rpv interaction can be rotated away , the divergent term is renormalized to zero .
the second line is however finite and can not be renormalized away .
this contribution does not give the same effective interaction as the renormalizable first term ( note the lorentz structure @xmath188 of the effective operator ) .
it is only generated at the loop level through the trilinear rpv interactions .
after renormalization , we obtain @xmath189 where we have also written the off - shell amplitudes given by the transposed and hermitian conjugated interactions . here
again we have neglected the quark and lepton masses . by inserting the above effective vertices to the second loop as shown in fig .
[ fig : rainbow_type ] ( 2 ) , we obtain the following contribution for the type 2 insertion ( expanded in the first order of the momentum @xmath190 carried by the external photon ) : @xmath191 \left [ ( k+q - p)^2 -m_{\tilde f_k}^2 \right ] \left [ ( k - p)^2 -m_{\tilde f_k}^2 \right ] } \int_0 ^ 1 dx \ , ( 1-x )   \ln \left ( \frac{xk^2 -m_{\tilde f_{1i}}^2 } { xk^2 -m_{\tilde f_{2i}}^2 } \right ) \nonumber\\ % % % % % % % % % % % % % % & \approx & \frac{2i}{(4\pi ) ^2 }
{ \rm i m } ( \hat \lambda_{ijj } \tilde \lambda^*_{ikk } e^{i(\theta_n -\delta_{f_j } ) } ) n_c q_f e ( g^{(n)}_{\tilde f_l } + g^{(n)}_{\tilde f_r } ) \sin \theta_{f_j } \cos \theta_{f_j } |m_{\lambda_n}| \epsilon^*_\mu ( q ) \nonumber\\ & & \hspace{4em}\times   \sum_{\tilde f = \tilde f_l , \tilde f_r } \int \frac{d^4 k}{(2\pi)^4 } \frac{-g^{(n)}_{\tilde f } m_{\tilde f_k}^2 \bar f_k \gamma_5 f_k \cdot ( 2p^\mu -q^\mu ) } { \left [ k^2 - |m_{\lambda_n}|^2 \right ] \left [ k^2 -m_{\tilde f_k}^2 \right]^3 } \int_0 ^ 1 dx \ , ( 1-x )   \ln \left ( \frac{xk^2 - m_{\tilde f_{1i}}^2 } { xk^2 - m_{\tilde f_{2i}}^2 } \right ) \nonumber\\ % % % % % % % % % % % % % % % & \approx & \frac{i } { ( 4\pi ) ^4 } { \rm i m } ( \hat \lambda_{ijj } \tilde \lambda^*_{ikk } e^{i(\theta_n -\delta_{f_j } ) } ) n_c q_f e ( g^{(n)}_{\tilde f_l } + g^{(n)}_{\tilde f_r } ) \sin \theta_{f_j } \cos \theta_{f_j } |m_{\lambda_n}| \epsilon^*_\mu ( q ) \bar f_k \sigma^{\mu \nu } q_\nu \gamma_5 f_k \nonumber\\ & & \hspace{8em } \times \sum_{\tilde f= \tilde f_l,\tilde f_r }   g^{(n)}_{\tilde f } m_{\tilde f_k}^2   \left [   f''(|m_{\lambda_n}|^2 , 0 , m_{\tilde f_k}^2 , m_{\tilde f_{1i}}^2 ) -f''(|m_{\lambda_n}|^2 , 0 , m_{\tilde f_k}^2 , m_{\tilde f_{2i}}^2 ) \right ] \ , \end{aligned}\ ] ] where the function @xmath70 is defined in eq .
( [ eq : f  ] ) , and the external momentum brought by the fermion @xmath50 is given by @xmath192 .
this leads to the formula ( [ eq : edmnug ] ) .
h. e. haber and g. l. kane , phys .
rept . * 117 * , 75 ( 1985 ) ; j. f. gunion and h. e. haber , nucl .
b * 272 * , 1 ( 1986 ) ; s. p. martin , in _
perspectives on supersymmetry ii _ , edited by g. l. kane ( world scientific , singapore , 2010 ) , p. 1
[ arxiv : hep - ph/9709356 ] .
j. ellis and m. k. gaillard , nucl .
b * 150 * , 141 ( 1979 ) ; e. p. shabalin , sov .
. phys . * 31 * , 864 ( 1980 ) ; i. b. khriplovich and a. r. zhitnitsky , phys .
b * 109 * , 490 ( 1982 ) .
i. b. khriplovich , sov . j. nucl
* 44 * , 659 ( 1986 ) ; phys .
b * 173 * , 193 ( 1986 ) ; m. e. pospelov and i. b. khriplovich , sov . j. nucl
. phys . * 53 * , 638 ( 1991 ) [ yad .
fiz . * 53 * , 1030 ( 1991 ) ] ; x .- g . he and b. mckellar , phys .
d * 46 * , 2131 ( 1992 ) ; x .- g . he , b. h. j. mckellar , s. pakvasa , phys .
b * 283 * , 348 ( 1992 ) ; a. czarnecki and b. krause , phys .
lett . * 78 * , 4339 ( 1997 ) ; t. mannel and n. uraltsev , phys .
d * 85 * , 096002 ( 2012 ) .
j. r. ellis , s. ferrera and d. v. nanopoulos , phys .
b * 114 * , 231 ( 1982 ) ; w. buchm@xmath193ller and d. wyler , phys .
b * 121 * , 321 ( 1983 ) ; j. polchinski and m. b. wise , phys .
b * 125 * , 393 ( 1983 ) ; f. del aguila , m. b. gavela , j. a. grifols and a. mendez , phys .
b * 126 * , 71 ( 1983 ) ; d. v. nanopoulos and m. srednicki , phys . lett .
b * 128 * , 61 ( 1983 ) ; m. dugan , b. grinstein and l. j. hall , nucl .
b * 255 * , 413 ( 1985 ) ; p. nath , phys .
* 66 * , 2565 ( 1991 ) ; y. kizukuri and n. oshimo , phys .
d * 45 * , 1806 ( 1992 ) ; phys .
d * 46 * , 3025 ( 1992 ) ; w. fischler , s. paban and s. d. thomas , phys .
b * 289 * , 373 ( 1992 ) ; t. inui , y. mimura , n. sakai and t. sasaki , nucl .
b * 449 * , 49 ( 1995 ) ; t. ibrahim and p. nath , phys .
d * 57 * , 478 ( 1998 ) ; phys .
b * 418 * , 98 ( 1998 ) ; phys .
d * 58 * , 111301 ( 1998 ) ; s. pokorski , j. rosiek and c. a. savoy , nucl .
b * 570 * , 81 ( 2000 ) ; s. y. ayazi and y. farzan , phys .
d * 74 * , 055008 ( 2006 ) ; jhep * 0706 * , 013 ( 2007 ) ;    t. h. west , phys .
d * 50 * , 7025 ( 1994 ) ; t. kadoyoshi and n. oshimo , phys . rev .
d * 55 * , 1481 ( 1997 ) ; d. chang , w .- y .
keung and a. pilaftsis , phys .
lett . * 82 * , 900 ( 1999 ) ; a. pilaftsis , phys .
b * 471 * , 174 ( 1999 ) ; d. chang , w .- f .
chang and w .- y .
keung , phys .
b * 478 * , 239 ( 2000 ) ; d. chang , w .- f . chang and w
keung , phys .
d * 66 * , 116008 ( 2002 ) ; t .- f .
feng , x .- q .
li , j. maalampi and x. zhang , phys .
d * 71 * , 056005 ( 2005 ) ; n. arkani - hamed , s. dimopoulos , g.f . giudice and a. romanino , nucl .
b * 709 * , 3 ( 2005 ) ; d. chang , w .- f .
chang and w .- y .
keung , phys .
d * 71 * , 076006 ( 2005 ) ; g. f. giudice and a. romanino , phys . lett .
b * 634 * , 307 ( 2006 ) ; t .- f .
feng , x .- q .
li , l. lin , j. maalampi and h .- s .
song , phys .
d * 73 * , 116001 ( 2006 ) ; y. li , s. profumo and m. j. ramsey - musolf , phys .
d * 78 * , 075009 ( 2008 ) .
j. dai , h. dykstra , r. g. leigh , s. paban and d. dicus , phys .
b * 237 * , 216 ( 1990 ) ; t. falk , k. a. olive , m. pospelov and r. roiban , nucl .
b * 560 * , 3 ( 1999 ) ; m. brhlik , g. j. good and g. l. kane , phys . rev .
d * 59 * , 115004 ( 1999 ) ; s. abel , s. khalil and o. lebedev , nucl .
b * 606 * , 151 ( 2001 ) ; a. pilaftsis , nucl .
b * 644 * , 263 ( 2002 ) ; o. lebedev and m. pospelov , phys .
lett . * 89 * , 101801 ( 2002 ) ; d. demir , o. lebedev , k. a. olive , m. pospelov and a. ritz , nucl . phys .
b * 680 * , 339 ( 2004 ) ; o. lebedev , k. a. olive , m. pospelov and a. ritz , phys . rev .
d * 70 * , 016003 ( 2004 ) ; k. a. olive , m. pospelov , a. ritz and y. santoso , phys .
d * 72 * , 075001 ( 2005 ) ; g. degrassi , e. franco , s. marchetti and l. silvestrini , jhep * 0511 * , 044 ( 2005 ) ; s. abel and o. lebedev , jhep * 0601 * , 133 ( 2006 ) ; j. ellis , j. s. lee and a. pilaftsis , jhep * 0810 * , 049 ( 2008 ) ; jhep * 1010 * , 049 ( 2010 ) ; jhep * 1102 * , 045 ( 2011 ) ; y. li , s. profumo and m. j. ramsey - musolf , jhep * 1008 * , 062 ( 2010 ) ; t. fukuyama , int . j. mod
a * 27 * , 1230015 ( 2012 ) .
j. hisano and y. shimizu , phys .
d * 70 * , 093001 ( 2004 ) ; m. endo , m. kakizaki and m. yamaguchi , phys .
b * 583 * , 186 ( 2004 ) ; g .- c .
cho , n. haba and m. honda , mod .
a * 20 * , 2969 ( 2005 ) ; j. hisano , m. nagai and p. paradisi , phys . lett .
b * 642 * , 510 ( 2006 ) ; phys .
d * 78 * , 075019 ( 2008 ) ; phys .
d * 80 * , 095014 ( 2009 ) ; w. altmannshofer , a. j. buras and p. paradisi , phys .
b * 688 * , 202 ( 2010 ) .
r. m. godbole , s. pakvasa , s. d. rindani and x. tata , phys .
d * 61 * , 113003 ( 2000 ) ; s. a. abel , a. dedes and h. k. dreiner , jhep * 0005 * , 13 ( 2000 ) ; d. chang , w .- f .
chang , m. frank and w .- y .
keung , phys .
d * 62 * , 095002 ( 2000 ) .
k. choi , e. j. chun and k. hwang , phys .
d * 63 * , 013002 ( 2000 ) ; y.y .
keum and otto c.w .
kong , phys . rev .
d * 63 * , 113012 ( 2001 ) ; c .- c .
chiou , o. c. w. kong and r. d. vaidya , phys .
d * 76 * , 013003 ( 2007 ) .
g. bhattacharyya , arxiv : hep - ph/9709395 ; h. k. dreiner , _ in perspectives on supersymmetry ii _ , edited by g. l. kane ( world scientific , singapore , 1997 ) , p. 565 ; r. barbier _
et al . _ ,
. rept . * 420 * , 1 ( 2005 ) ; m. chemtob , prog . part .
phys . * 54 * , 71 ( 2005 ) .
atlas collaboration ( georges aad _ et al .
_ ) , phys . rev .
* 108 * , 261804 ( 2012 ) ; arxiv:1208.4688 [ hep - ex ] ; arxiv:1210.1314 [ hep - ex ] ; cms collaboration ( vardan khachatryan _
_ ) , phys . rev .
lett . * 107 * , 221804 ( 2011 ) ; arxiv:1210.2052 [ hep - ex ] ; arxiv:1210.8115 [ hep - ex ] .