nonlinearities prevail in fluid dynamics when the reynolds number , @xmath4 , is large @xcite . for geophysical flows ,
the reynolds number is often larger than @xmath5 and for some astrophysical flows values of @xmath6 is not unreasonable .
the number of degrees of freedom ( _ dof _ ) in the flow increases as @xmath7 for @xmath8 in the kolmogorov framework @xcite ( hereafter k41 ) .
such a huge number of _ dof _ makes direct numerical simulations ( dns ) of turbulence at high @xmath4 infeasible on any existing or projected computer for decades to come . because of this intractability
, simulations of turbulence are always carried out in regions of parameter space far from the observed values , either with : ( a ) an unphysical lack of scale separation between the energy - containing , inertial , and dissipative ranges while parameterizing the missing physics , or ( b ) a study of the processes at much smaller length scales , often with periodic boundaries ( unphysical at large scales but used under the hypothesis of homogeneity of turbulent flows ) .
clearly , modeling of unresolved small scales is necessary .
two main approaches have been developed over the years to model the effects of the unresolvable small scales in turbulence on the scales resolved in the simulations .
the first approach is large eddy simulations ( les , see @xcite ) .
les is widely used in engineering , in atmospheric sciences , and to a lesser extent in astrophysics . however , in the les approach , the reynolds number is not known .
instead , one attempts modeling the behavior of the flow in the limit of very large @xmath4 .
another approach models the effects of turbulence at higher reynolds numbers than are possible with a dns on a given grid , by using a variety of techniques that can be viewed as filtering of the small scales    a rather novel approach to modeling of turbulent flows employs _ regularization modeling _ as a sgs model @xcite . unlike closures which employ eddy - viscosity concepts ( modifying the dissipative processes ) , the approach of regularization modeling modifies the spectral distribution of energy .
existence and uniqueness of smooth solutions can be rigorously proven , unlike many les models ( e.g. , eddy - viscosity ) , as well as the fact that the subgrid model recovers the navier - stokes equations in the limit of the filter width going to zero .
their robust analytical properties ensure computability of solutions . proved the existence ( but not smoothness , or uniqueness ) of solutions to the navier - stokes equations in @xmath9 , n=2,3 using the leray model @xcite .    geurts and
holm @xcite began using the leray model with ( a three - point invertible approximation of ) an inverse - helmholtz - operator filter of width @xmath2 .
later it was dubbed leray@xmath0 and an upper bound for the dimension of the global attractor was established @xcite .
the global existence and uniqueness of strong solutions for the leray model is a classical result @xcite .
leray@xmath0 has been compared to dns simulations on a grid of @xmath10 in a doubly - periodic compressible channel flow domain @xcite .
its performance was found to be superior to a dynamic mixed ( similarity plus eddy - viscosity ) model ( with an even greater reduction in computational cost ) .
however , it possessed a systematic error of a slight over - prediction of the large scales accompanied by a slight under - prediction of the small scales .
it possessed both forward- and back - scatter , but exhibited too little dissipation .
the leonard tensor - diffusivity model @xcite ( sometimes known as the clark model @xcite ) is the first term of the reconstruction series for the turbulent sub - filter stress for all symmetric filters that possess a finite nonzero second moment .
this leading order approximation of the subgrid stress is thus generic @xcite . in
_ a priori _ testing , it reconstructs a significant fraction ( @xmath11 , but not all ) of the subgrid stress , provides for local backscatter along the stretching directions while remaining globally dissipative , and possesses a better reconstruction of the subgrid stress than the scale - similarity model .
used as a les , in _ a posteriori _ testing the leonard tensor - diffusivity model required additional dissipation ( a dynamic smagorinsky term ) to achieve reasonable gains in computation speed for 3d periodic flows and for channel flows @xcite .
the leonard tensor - diffusivity model does not conserve energy in the non - viscous limit .
cao , holm , and titi @xcite developed a related ( conservative ) subgrid model which they dubbed clark@xmath0 .
the clark@xmath0 model applies an additional inverse - helmholtz filter operation to the reynolds stress tensor of the clark model .
the global well - posedness of the clark@xmath0 model and the existence and uniqueness of its solutions were demonstrated , and upper bounds for the hausdorff ( @xmath12 ) and fractal ( @xmath13 ) dimensions of the global attractor were found @xcite .    the third regularization model we will consider is the incompressible lagrangian - averaged navier - stokes ( lans@xmath0 , @xmath14model , also known as the viscous camassa - holm equation @xcite ) .
it can be derived , for instance , by applying temporal averaging to hamilton s principle , where taylor s frozen - in turbulence hypothesis ( the only approximation in the derivation ) is applied as the closure for the eulerian fluctuation velocity in the reynold decomposition , at linear order in the generalized lagrangian mean description @xcite . in this derivation , the momentum - conservation structure of the equations is retained . for scales smaller than the filter width
, lans@xmath0 reduces the steepness in gradients of the lagrangian mean velocity and thereby limits how thin the vortex tubes may become as they are transported , while the effect on larger length scales is negligible @xcite . lans@xmath0
may also be derived by smoothing the transport velocity of a material loop in kelvin s circulation theorem @xcite .
consequently , there is no attenuation of resolved circulation , which is important for many engineering and geophysical flows where accurate prediction of circulation is highly desirable .
an alternative interpretation of the @xmath14model is that it neglects fluctuations in the smoothed velocity field , while preserving them in the source term , the vorticity @xcite .
lans@xmath0 has previously been compared to direct numerical simulations ( dns ) of the navier - stokes equations at modest taylor reynolds numbers ( @xmath15 @xcite , @xmath16 @xcite , and @xmath17 @xcite ) .
lans@xmath0 was compared to a dynamic eddy - viscosity les in 3d isotropic turbulence under two different forcing functions ( for @xmath18 and @xmath19 ) and for decaying turbulence with initial conditions peaked at a low wavenumber ( with @xmath20 ) as well as at a moderate wavenumber ( with @xmath21 ) @xcite .
lans@xmath0 was preferable in these comparisons because it demonstrated correct alignment between eigenvectors of the subgrid stress tensor and the eigenvectors of the resolved stress tensor and vorticity vector .
the lans@xmath0 and the leray@xmath0 dynamic mixed model ( similarity plus eddy - viscosity ) in a turbulent mixing shear layer ( with @xmath22 ) @xcite .
lans@xmath0 was found to be the most accurate of these three les candidates at proper subgrid resolution , but the effects of numerical contamination can be strong enough to lose most of this potential . while lans@xmath0 has the greatest grid - independent accuracy of the three models , it also requires the greatest resolution . from the les perspective
, this could pose some limitations on the practical use and application of lans@xmath0 for high @xmath4 cases .
indeed , recent high - resolution simulations of lans@xmath0 showed that energy artificially accumulates in the sub - filter - scales , giving as a result only a modest computational gain at very high reynolds number @xcite .
we propose to pursue these previous studies of leray@xmath0 and lans@xmath0 further at high reynolds number , one goal is to as the three regularizations are related via truncation of sub - filter stresses , for lans@xmath0 , the predicted -scale spectra is @xmath23 @xcite .
this scaling has been observed to be subdominant to an energy spectrum @xmath3 which corresponds to `` enslaved rigid body '' or `` polymerized '' portions of the fluid @xcite .
the sub - filter scaling observed in the third - order structure function corresponded to the predicted @xmath24 scaling of the energy spectrum .
however , regions were observed in the flow where no stretching was acting in the sub - filter scales .
these regions , which give no contribution to the energy cascade , and hence do not affect the third order structure functions , are responsible for the @xmath3 scaling in the lans@xmath0 energy spectrum .
for clark@xmath0 , the correct time scale for vortex stretching is difficult to determine and its spectrum is found to range between @xmath23 and @xmath25 @xcite .
leray@xmath0 has the same difficulty and the spectrum can range between @xmath26 and @xmath27 @xcite .
the determination of these scaling laws is needed to quantify the computational gain if each model is to be used as a sgs model . as a result , to determine empirically the sub - filter scale spectra .
we present the three models and describe how they are related , derive the krmn - howarth equation for clark@xmath0 and leray@xmath0 ( from which exact scaling laws for third order quantities follow ) , and review theoretical predictions of inertial range scaling in section [ sec : models ] . in section [ sec : sfsstudy ] .
we first compute a fully resolved dns of the navier - stokes equations at a resolution of @xmath28 ( @xmath29 , @xmath30 , and @xmath31 ) .
we then perform model runs with the exact same conditions at a resolution of @xmath32 .
we take @xmath2 to be @xmath33th the box size , which was found in an earlier study to be large enough to exhibit both navier - stokes and sub - filter - scale lans@xmath0 dynamics @xcite .
this large filter case is important because it gives insight into the behavior of the models at scales much smaller than the filter width without requiring higher resolution than is feasible .
we compare the three regularizations as subgrid models guided by a previous study of lans@xmath0 @xcite , we take @xmath2 to be @xmath34th the box size .
this choice was found to produce an optimal @xmath14les ( in the sense of being optimal for the class of lans@xmath0 models , with respect to the value of @xmath2 ) .
finally , we review bounds on the size of the attractors and use these bounds to comment on the computational savings of the three regularizations viewed as sgs models .
the incompressible navier - stokes equations are given in cartesian coordinates by @xmath35 filtering these equations with a convolution filter , @xmath36 in which @xmath37 ( resp .
@xmath38 ) denotes the filtered ( resp .
unfiltered ) field , yields @xmath39 in which by convention we denote @xmath40 and the reynolds turbulence stress tensor , @xmath41 , represents the closure problem .
it has been shown @xcite that for all symmetric filters possessing a finite nonzero second moment , the first term of the reconstruction series for the turbulent sub - filter stress is @xmath42 where @xmath43 is the fourier transform of the convolution kernel ( @xmath44 is the convolution kernel where @xmath45({\mathbf}{r})=\int g({\mathbf}{r}-{\mathbf}{r { ' } } ) z({\mathbf}{r { ' } } ) d^3{\mathbf}{r{'}}$ ] ) .
this approximation of the subgrid stress is then generic and is known as the leonard tensor - diffusivity model @xcite ( or , often , the clark model @xcite ) .    related to this model , the clark@xmath0 model is @xcite @xmath46 and its subgrid stress for @xmath47 is given by @xmath48 here the filter is the inverse of a helmholtz operator , @xmath49 .
the clark@xmath0 model conserves energy in the @xmath50 norm instead of the @xmath51 norm , @xmath52 in which the clark@xmath0 energy @xmath53 is expressed as @xmath54 and the clark@xmath0 energy dissipation rate is given by @xmath55 where @xmath56 and @xmath57 . for @xmath58
, we note that @xmath59 which implies that the turbulent sub - filter stress tensor for the tensor - diffusivity model given by eq .
( [ eq : carati ] ) is @xmath60 which is proportional to the clark@xmath0 stress tensor to second order in @xmath2 .
hence , the _ a priori _ tests of @xcite should apply to clark@xmath0 , at least at this order .      in 1938 ,
krmn and howarth @xcite introduced the invariant theory of isotropic hydrodynamic turbulence , and derived from the navier - stokes equations the exact law relating the time derivative of the two - point velocity correlation to the divergence of the third - order correlation function .
the corresponding krmn - howarth theorem for lans@xmath0 in the fluid case was derived in @xcite .
the relevance of the krmn - howarth theorem for the study of turbulence can not be underestimated . as a corollary ,
rigorous scaling laws in the inertial range can be deduced . in this section ,
we derive these results for the clark@xmath0 case .    for the sake of simplicity
, we consider the case @xmath61 , since the dissipative terms may be added at any point in the derivation .
we denote @xmath62 and begin our investigation of the correlation dynamics by computing the ingredients of the partial derivative @xmath63 .
the clark@xmath0 motion equation ( [ eq : clark ] ) may be rewritten as @xmath64 combining eqs .
( [ eq : subgrid ] ) and ( [ eq : sgs_clark ] ) , we arrive at the fluctuation - velocity equation , @xmath65 where @xmath66 ( @xmath67 ) . multiplying eq .
( [ eq : clarkv2 ] ) by @xmath68 and eq .
( [ eq : clarku2 ] ) by @xmath69 , then adding the result yields @xmath70 where we have used statistical homogeneity @xmath71 we symmetrize eq .
( [ fluct - eqn ] ) in the indices @xmath72 by adding the corresponding equation for @xmath73 .
we then use homogeneity again as @xmath74 and define the tensors @xmath75 we can drop @xmath76 because the terms with the pressures @xmath77 and @xmath78 vanish everywhere , as follows from the arguments of isotropy @xcite . finally , we obtain @xmath79 this is the krmn - howarth equation for clark@xmath0 ( compare to eq .
( 3.8 ) in @xcite for lans@xmath0 ) .    by dimensional analysis , the energy dissipation rate in clark@xmath0 is @xmath80 and eq .
( [ eq : kh1 ] ) implies @xmath81 for large scales ( @xmath82 ) , we recover the navier - stokes scaling known as the four - fifths law , @xmath83 @xcite . here
, @xmath84\cdot{\mathbf}{l}/l$ ] is the longitudinal increment of @xmath85 . strictly speaking
, the four - fifths law expresses that the third - order longitudinal structure function of @xmath85 , @xmath86 , is given in the inertial range in terms of the mean energy dissipation per unit mass @xmath87 by @xmath88 or , equivalently , that the flux of energy across scales in the inertial range is constant .
we also recover the kolmogorov 1941 @xcite ( hereafter , k41 ) energy spectrum , @xmath89 or , equivalently , @xmath90 for sub - filter scales ( @xmath91 ) , we have @xmath92 and the first and third right - hand terms in eq .
( [ eq : full_clark ] ) are equivalent . in this case , we are left with two different possible scalings depending on the prefactors in eq .
( [ eq : full_clark ] ) .
if the first ( or third ) right - hand term is dominant , our scaling law becomes @xmath93 for our sub - filter scale energy spectrum we would then have @xmath94 , or , equivalently , @xmath95 this result is the same as for the @xmath14model @xcite . if , however , the second right - hand term in eq .
( [ eq : full_clark ] ) is dominant , then the k41 results are recovered , with @xmath96 substituted for @xmath85 . in that case , one finds the alternative clark@xmath0 sub - filter scale spectral energy scaling , @xmath97      we review here the derivation by dimensional analysis of the spectrum which follows the scaling ideas originally due to kraichnan @xcite and which is developed more fully in ref .
@xcite . in examining the nonlinear terms in eq .
( [ eq : clark ] ) , it is not entirely clear which of three possible scales for the average velocity for an eddy of size @xmath1 , @xmath98 @xmath99 or @xmath100 should result .
therefore , three corresponding `` turnover times '' , @xmath101 , for such an eddy may be proposed @xmath102 the term `` turnover time '' is used advisedly here , since only the velocity @xmath103 is composed of the fluid transport velocity .
we define the ( omnidirectional ) spectral energy density , @xmath104 , from the relation @xmath105 since , @xmath106 , we have @xmath107 then , the total energy dissipation rate , @xmath108 is related to the spectral energy density by @xmath109 which yields , finally , the predicted energy spectra for clark@xmath0 , @xmath110 , @xmath111 for scales much larger than @xmath2 ( @xmath112 ) the kolmogorov scaling for navier - stokes is recovered , @xmath113 whereas for scales much smaller ( @xmath114 ) , the spectrum becomes @xmath115 these arguments constrain the clark@xmath0 sub - filter scale spectrum to lie between @xmath1 and @xmath116 .
the leray model in cartesian coordinates is @xmath117 where the flow is advected by a smoothed velocity , @xmath96 . by comparison with eq .
( [ eq : subgrid ] ) we see that the leray model approximates the subgrid stress as @xmath118 , or , with @xmath119 , @xmath120 as has been noted previously @xcite , the subgrid stress of clark@xmath0 in eq .
( [ eq : sgs_clark ] ) is a truncation of the subgrid stress of leray@xmath0 , in eq .
( [ eq : sgs_leray ] ) . for leray@xmath0 ,
the @xmath121 norm is the quadratic invariant that is identified with energy , @xmath122 where @xmath123 and @xmath124 as was pointed out in @xcite , the incompressibility of the velocity field @xmath85 only implies a divergenceless filtered velocity @xmath96 under certain boundary conditions for leray@xmath0 .
when @xmath125 , the energy @xmath126 is no longer conserved ( helicity and kelvin s theorem are not conserved for leray@xmath0 ) . in our numerical study , we employ periodic boundary conditions , for which @xmath127 implies @xmath128 and leray@xmath0 conserves energy in the usual sense of @xmath51 .      in this section
we derive the krmn - howarth equation for the leray@xmath0 case . following section [ sec : kh - clark ]
, we begin our investigation of the correlation dynamics by computing the ingredients of the partial derivative @xmath129 .
( [ eq : leray ] ) may be rewritten as : @xmath130 multiplying eq .
( [ eq : lerayv2 ] ) by @xmath131 yields @xmath132 we can make this equation symmetric in the indices @xmath72 , by adding the equation for @xmath133 .
we define the tensors @xmath134 again , we may drop @xmath135 because the terms with the pressure @xmath136 vanish everywhere and thereby obtain @xmath137 this is the krmn - howarth equation for leray@xmath0 .    the energy dissipation rate for leray@xmath0
is denoted by @xmath138 and it satisfies @xmath139 . by dimensional analysis , eq .
( [ eq : kh1leray ] ) implies @xmath140 for large scales ( @xmath82 ) , we recover the navier - stokes scaling , eqs .
( [ eq : fourfifths ] ) and ( [ eq : k41 ] ) . for sub - filter scales ( @xmath91 )
our scaling law becomes @xmath141 for our small scale energy spectrum we would then have @xmath142 ( where we employed @xmath92 ) , or , equivalently , cf .
( [ eq : clarka_spectrum2 ] ) , @xmath143      we review here the derivation by dimensional analysis of the spectrum for leray@xmath0 as we did for clark@xmath0 in section [ sec : pheno_clark ] .
this analysis is developed more fully in ref .
we argue again that there are three possible scales for the average velocity for an eddy of size @xmath1 , eqs .
( [ eq : uk0 ] ) , ( [ eq : uk1 ] ) , and ( [ eq : uk2 ] ) , with the turn - over time , @xmath144 given by eq .
( [ eq : turnover ] ) . since , @xmath145 , we have @xmath146 then , the total energy dissipation rate , @xmath138 is related to the spectral energy density by @xmath147 which yields , finally , the predicted energy spectra for leray@xmath0 , @xmath148 , @xmath149 for scales much larger than @xmath2 ( @xmath112 ) the k41 spectrum is recovered , eq .
( [ eq : k41 ] ) , and for scales much smaller ( @xmath114 ) the spectrum is @xmath150 these arguments constrain the leray@xmath0 sub - filter scale spectrum to lie between @xmath151 and @xmath152 .
lans@xmath0 is given by @xmath153 for lans@xmath0 , the usual choice of filter is again @xmath154 with this filter , the subgrid stress tensor is given by @xmath155 as has been previously noted @xcite , the subgrid stress of leray@xmath0 , eq .
( [ eq : sgs_leray ] ) , is a truncation of the subgrid stress of lans@xmath0 eq .
( [ eq : sgs_lans ] ) . like clark@xmath0
, energy is conserved in the @xmath156 norm instead of the @xmath51 norm .
additionally , lans@xmath0 is the only model of the three examined here that conserves the helicity ( and kelvin s circulation theorem ) .    for lans@xmath0 in the fluid case
the krmn - howarth theorem was derived in @xcite .
we summarize here the dimensional analysis argument for the lans@xmath0 inertial range scaling that follows from this theorem , beginning from equation ( 3.8 ) in ref .
@xcite . in the statistically isotropic and homogeneous case , without external forces and with @xmath157 ,
taking the dot product of eq .
( [ eq : lans ] ) with @xmath158 yields the equation    @xmath159    the trace of this equation is the fourier transform of the detailed energy balance for lans@xmath0 ; @xmath160 is the second - order correlation tensor while @xmath161 and @xmath162 are the third - order correlation tensors for lans@xmath0 and @xmath163 is the sub - filter scale stress tensor . for @xmath164
this reduces to the well - known relation derived by krmn and howarth .
the energy dissipation rate for lans@xmath0 , @xmath165 , satisfies @xmath166 . by dimensional analysis in eq .
( [ eq : holm38 ] ) we arrive at @xmath167 for large scales ( @xmath82 ) , we recover the navier - stokes scaling eqs .
( [ eq : fourfifths ] ) and ( [ eq : k41 ] ) . for sub - filter scales ( @xmath91 )
our scaling law becomes eq .
( [ eq : lcube ] ) and our sub - filter scale spectra is given by @xmath168 in this case , by the phenomenological arguments , we know that eddies of size @xmath1 are advected by the smoothed velocity , eq .
( [ eq : uk2 ] ) .
this scaling is confirmed in ref .
but it coexists with a @xmath169 energy spectrum corresponding to `` enslaved rigid bodies '' or `` polymerized '' portions of fluid which do not contribute to the turbulent energy cascade .
to this end , we compute numerical solutions to eqs . ( [ eq : ns ] ) , ( [ eq : clark ] ) , ( [ eq : leray ] ) , and ( [ eq : lans ] ) in a three - dimensional ( 3d ) cube with periodic boundary conditions using a parallel pseudospectral code @xcite .
we employ a taylor - green forcing @xcite , @xmath170 \label{eq : tg}\ ] ] ( with @xmath171 ) , and employ dynamic control @xcite to maintain a nearly constant energy with time .
the taylor - green forcing , eq .
( [ eq : tg ] ) , is not a solution of the euler s equations , and as a result small scales are generated rapidly .
the resulting flow models the fluid between counter - rotating cylinders @xcite and it has been widely used to study turbulence , including studies in the context of the generation of magnetic fields through dynamo instability @xcite . we define the taylor microscale as @xmath172 and the mean velocity fluctuation as @xmath173 .
the taylor microscale reynolds number is defined by @xmath174 and the reynolds number based on a unit length is @xmath175 .    to compare the three regularizations ( clark@xmath0 , leray@xmath0 , and lans@xmath0 ) we compute a fully resolved dns of the navier - stokes equations at a resolution of @xmath28 ( @xmath29 , @xmath176 ) and model runs with the exact same conditions at a resolution of @xmath32 .
the details of the flow dynamics of the dns have already been given @xcite .
in particular , the reynolds number based on the integral scale @xmath177 ( where @xmath178 is the total energy ) is @xmath179 , where @xmath180 is the _ r.m.s . _
velocity and the reynolds number based on the taylor scale is @xmath31 .
the dns was run for nine turnover times @xmath181 ( in the following results , time @xmath182 is in units of the turnover time ) .
we employ a filter width of @xmath183 for which lans@xmath0 exhibits both navier - stokes and lans@xmath0 inertial ranges in the third - order structure function @xcite . from these
we hope to obtain the behavior of the models for scales much smaller than @xmath2     for leray@xmath0 and in the dns , @xmath184 for lans@xmath0 and clark@xmath0 ) .
dns ( @xmath176 ) is shown as solid black lines , lans@xmath0 as dotted red , clark@xmath0 as dashed green , and leray@xmath0 as blue dash - dotted .
the cyan long - dashed line represents a @xmath185 dns ( @xmath186 , @xmath187 ) . here
each run is calculated only until it reaches a statistical steady state .
leray@xmath0 reduces the dissipation , @xmath188 , and increases the time scale to reach a statistical turbulent steady - state .
both effects are greater as @xmath2 is increased . by comparison with the @xmath189 run ,
we see that these two effects are consistent with a reduced effective reynolds number .
a smaller reduction in flux ( but not an increase in time to steady - state ) is also observed for lans@xmath0 and is likely related to its rigid bodies.,width=338 ]    in fig .
[ fig : enstrophy ] we present the time evolution of the enstrophy ( @xmath190 for leray@xmath0 and dns , @xmath184 for lans@xmath0 and clark@xmath0 ) , which is proportional to the dissipation ( @xmath191 or @xmath192 depending on the case ) .
also shown is a well - resolved @xmath185 dns of a less turbulent flow ( @xmath193 , @xmath194 , @xmath187 , ( cyan online ) long - dashed line ) . here
each run is calculated only until it reaches a statistically steady state .
we see that the dissipation is reduced and the time - scale is increased for leray@xmath0 .
we see by comparison with the @xmath194 dns that a reduced effective reynolds number in the leray@xmath0 run . for lans@xmath0 ,
the dissipation is decreased although the time to reach steady state is not increased .
this is probably related to the enslavement of its rigid body regions which would have no internal dissipation .
of the three models , clark@xmath0 the total dissipation for a large range of @xmath2 .
indeed , as it is the order @xmath195 approximation of navier - stokes , this dissipation behavior for clark@xmath0 may continue to hold until @xmath2 becomes quite large .
figure [ fig : eavg ] shows a comparison of the energy spectrum at the turbulent steady state , for all the runs in fig .
[ fig : enstrophy ] .
the isotropic energy spectra are calculated as follows , @xmath196 where @xmath197 ( the @xmath198 norm is employed for clark@xmath0 and lans@xmath0 ) .
the length scale @xmath2 is indicated by a vertical dashed line and the plotted energy spectra are compensated by @xmath199 ( i.e. , leading to a flat k41 @xmath152 spectrum ) . from fig .
[ fig : eavg ] it is clear that the clark@xmath0 spectral behavior is close to the predicted @xmath1 spectrum , rather than the @xmath200 from eq .
( [ eq : clarka_spectrum2 ] ) , or the other possible spectrum from eq .
( [ eq : clarka_spectrumn ] ) .
likewise , the @xcite .
leray@xmath0 , on the other hand , possesses a very steep sub - filter scale spectrum as well as enhanced large - scale energy as has been previously observed @xcite .    the spectrum of leray@xmath0 in fig .
[ fig : eavg ] gives a good approximation to the @xmath189 dns in the range @xmath201 $ ] ( i.e. , to @xmath193 rather than to @xmath202 which was employed ) .
this result , leray@xmath0 s increased characteristic time scales and its reduced dissipation imply that the leray@xmath0 model is operating at a much lower effective reynolds number .
this is also clear from the rapid drop in the spectrum at small scales , shown in fig .
[ fig : eavg ] .
indeed , we can build an effective reynolds number in the large scales as @xmath203 .
since @xmath204 is controlled in this simulation by the forcing scale , the drop in the dissipation rate implies a reduced nonlinearity in leray@xmath0 .
this is also consistent with a direct comparison of the nonlinear terms in leray@xmath0 with , for instance , lans@xmath0 .
the nonlinear terms in lans@xmath0 , eq .
( [ eq : lans ] ) , may be written as @xmath205 ( where the suffix @xmath206 denotes a transposition ) , while the nonlinear term in leray@xmath0 , eq .
( [ eq : leray ] ) , is only @xmath207 .
both nonlinear terms in lans@xmath0 are of order @xmath208 ; so the absence of one of the nonlinear terms in leray@xmath0 be understood as a reduction in the nonlinearity .
validation of the krmn - howarth equation scalings , eqs .
( [ eq : lcube ] ) and ( [ eq : lcubeleray ] ) , enables us to measure scaling laws in the inertial range and , thus , compare the intermittency properties of the models .
the third - order correlations involved in the theorems , namely @xmath209 for clark@xmath0 and lans@xmath0 , and @xmath210 for leray@xmath0 , and @xmath211 for navier - stokes , are plotted versus @xmath212 in fig .
[ fig : theorem ] . in fig .
[ fig : theorem ] we can see validation of the krmn - howarth scaling for scales smaller than @xmath2 for both lans@xmath0 and clark@xmath0 . in particular , we note the observed scaling for clark@xmath0 verifies the @xmath213 scaling and not the ( theoretically possible ) @xmath214 ( @xmath215 ) scaling . the predicted scaling is not observed in leray due to its reduced effective reynolds number . with these scalings in hand
, we may proceed to observe the scaling of the longitudinal structure functions , @xmath216 where we again replace the @xmath217 norm , @xmath218 , for the @xmath219 norm , @xmath220 , in the case of clark@xmath0 and lans@xmath0 .
we utilize the extended self - similarity ( ess ) hypothesis @xcite which proposes the scaling @xmath221 and normalize the results by @xmath222 to better visualize the deviation from linearity ( which serves as a measure of intermittency ) .
the results are displayed in fig .
[ fig : exponents ] .    in fig .
[ fig : exponents ] , we may observe the intermittency properties of the models at sub - filter scales .
we note a reduced intermittency for both leray@xmath0 and the @xmath194 dns .
this is consistent with the smoother , more laminar fields ( due to the reduction of the effective @xmath4 ) possessed by both .
interestingly , though lans@xmath0 and clark@xmath0 both possess the same cascade scaling ( eq . ( [ eq : lcube ] ) ,
as confirmed in fig .
[ fig : theorem ] ) , the clark@xmath0 model is markedly more intermittent than lans@xmath0 . if artificially truncated local interactions ( in spectral space ) is taken as a cause of enhanced intermittency @xcite , then the increased intermittency observed in clark@xmath0 is the expected result of truncation of the higher - order terms in the sub - filter stress tensor . moreover , if the lans@xmath0 s @xmath3 spectrum is indeed associated with rigid bodies , these would serve to decrease the intermittency ( no internal degrees of freedom being available in a rigid body ) which is consistent with the results shown here . due to this effect , lans@xmath0 of the three regularization models the high - order intermittency of
@xmath223 corresponding to an optimal @xmath14les @xcite .
note that the value of @xmath2 has been optimized for neither clark@xmath0 nor leray@xmath0 ; as a consequence , these models might perform better in other parameter regimes than the results indicate in this study .     for leray@xmath0 and in the dns , @xmath184 for lans@xmath0 and clark@xmath0 ) .
dns ( @xmath176 ) is shown as solid black lines , lans@xmath0 as dotted red , clark@xmath0 as dashed green , and leray@xmath0 as blue dash - dotted .
an under - resolved ( @xmath32 ) navier - stokes run is shown as a pink dash - triple - dotted line.,width=338 ]    fig .
[ fig : enstrophy2 ] along with that of an under - resolved navier - stokes solution at a resolution of @xmath32 ( @xmath202 , ( pink online ) dash - triple - dotted line ) .
we see that both lans@xmath0 and clark@xmath0 reproduce the proper amount of dissipation and are within @xmath224 of the time required by the dns to reach a statistical turbulent steady state .
as has been observed before , leray@xmath0 is under - dissipative @xcite .
we also note that it takes longer than the    compensated spectra averaged over several eddy turn - over times are shown for the sgs case ( i.e. , @xmath225 ) in fig .
[ fig : eavg2 ] .
note that as the subgrid models are averaged over a different time interval , no meaningful comparison to the dns is possible for @xmath226 . even without an optimal choice for the value of @xmath2 ,
clark@xmath0 best reproduces the dns spectrum for scales larger than @xmath2 .
we compute root - mean - square spectral errors where @xmath227 is the wavenumber for the forcing scale , @xmath228 is the dns spectrum ( in the @xmath121 norm ) , and @xmath229 is the subgrid model spectrum ( in the appropriate norm ) .
we see that only clark@xmath0 improves the estimate of the power spectrum at this resolution for leray@xmath0 performs the poorest of the three regularization models , but it is also not optimized .
as previously argued , its effective reynolds number is too low to accurately model the dns flow .
either a decrease in the viscosity @xmath230 , or a decrease in the filter size @xmath2 ( and , hence , an increase in the nonlinearity ) , or both would likely improve the accuracy of leray@xmath0 as an sgs model .
due to its frozen - in ( or enslaved ) rigid - body regions , the lans@xmath0 model can not reproduce the dns spectrum unless @xmath2 is only a few times larger than the dissipation scale @xcite .
[ cols="^,^,^ " , ]     another measure of the success of a subgrid model is the reproduction of structures in the flow . in figure
[ fig : vapor ] we have 3d volume rendering of the enstrophy density @xmath231 ( @xmath232 for lans@xmath0 and clark@xmath0 ) for the dns , the three sgs - model simulations ( @xmath233 ) , the @xmath32 under - resolved navier - stokes solution , all at a reynolds number of @xmath234 , and the @xmath194 dns . due to the late times depicted ( longer than a lyapunov time )
there can be no point - by - point comparison between the simulations .
instead , we note that there are four horizontal bands where the forcing causes a maximum shear .
this large - scale feature of the flow is missing only from leray@xmath0 and the @xmath194 run .
the three other runs reproduce this feature well ( note that the apparently thicker tubes present in clark@xmath0 are vortex tube mergers ) .
the results lead again to the conclusion that the under - resolved navier - stokes , the clark@xmath0 , and the lans@xmath0 models are better subgrid models than leray@xmath0 due to its reduced effective @xmath4 .    for the sgs models ,
the predicted @xmath235 from the krmn - howarth theorem for navier - stokes is well - reproduced by all models at super - filter scales . in fig .
[ fig : exponents2 ] the sgs model intermittency results .
all models reproduce the intermittency up to the tenth - order moment within the error bars ( although there is a small decrease in intermittency for leray@xmath0 ) .
thus , we conclude that with adequately chosen values of @xmath2 ( and of @xmath230 for leray@xmath0 ) , all three models can reproduce the intermittency of the dns ( to within the error bars ) .
the rationale behind using is that it leads to adequate solutions at a reduced computational cost , because it computes fewer _ dof _ ; indeed , for an , the ratio of navier - stokes s _ dof _ to the model s _ dof _ , a prediction for memory savings and hence computation time savings for numerical simulation , is a crucial factor .
consequently , analytical bounds on the sizes of the attractors for the three regularization subgrid models may be useful indicators of their computational savings .
the _ dof _ for lans@xmath0 is derived in @xcite and confirmed in @xcite , @xmath236 where @xmath237 is the integral scale ( or domain size ) .
we may compare this to the _ dof _ for navier stokes , @xmath238 which immediately yields @xmath239 it was found , however , that to reproduce the energy spectrum of an equivalent dns , the filter - width @xmath2 must be no larger than a few times the dissipation scale , @xmath240 @xcite .
this is the result of the `` polymerization '' of the flow in lans@xmath0 , and the associated @xmath241 scaling at sub - filter scales . with this added caveat
, it follows that the reduction in _ dof _ is independent of @xmath4 ( and a net factor of about @xmath242 ) .
our study here illustrates that the high - order structure functions may be reproduced for much larger values of @xmath2 .
therefore , in applications where the spectrum is not of great concern , much greater reduction in numerical resolution would be feasible .
for clark@xmath0 there is an upper bound on the hausdorff , @xmath12 , and fractal , @xmath13 , dimensions of the attractor , @xmath243 where @xmath244 is the kolmogorov dissipation length scale corresponding to the clark@xmath0 model @xcite . from its observed @xmath1 spectrum , we may estimate @xmath244 or , equivalently , @xmath245 .
for dissipation the large wavenumbers dominate and , therefore , combining the clark@xmath0 energy balance eq .
( [ eq : lansa_balance ] ) with its sub - filter scale energy spectrum eq .
( [ eq : clarka_spectrum ] ) allows us to implicitly specify its dissipation wavenumber , @xmath246 , by @xmath247 then we have , @xmath248 it follows that @xmath249 this is similar to the prediction for lans@xmath0 , but as energy spectra are more easily reproduced for larger values of @xmath2 than with lans@xmath0 ( but not the intermittency properties ) , it may be the case that @xmath2 is not tied to the kolmogorov dissipation scale @xmath240 .
if so , then the computational saving might increase as @xmath250 which is promising for use of clark@xmath0as an les model .
this conclusion is bolstered to the extent that the results in section [ sec : super ] for @xmath233 ( @xmath251 ) are acceptable .
if , there is still a greater reduction in _ dof _ ( a factor of 20 ) for clark@xmath0 than for lans@xmath0 .    for leray@xmath0 , we have the following upper bounds on the hausdorff dimension , @xmath12 , and fractal dimension , @xmath13 , of the global attractor , @xmath252 where @xmath253 is the dissipation length scale for leray@xmath0 @xcite .
again , we estimate the dissipation wavenumber for leray@xmath0 @xmath254 . then , from eqs .
( [ eq : leray_balance ] ) and ( [ eq : leray_spectrum ] ) , that is , assuming the @xmath151 spectrum resulting from the krmn - howarth equation , we find @xmath255 consequently we have , @xmath256 it follows that @xmath257 our results suggest that for an effective les the viscosity @xmath258 must be chosen to be _ smaller _ than @xmath230 .
this leads to an upper bound on the computational savings for leray@xmath0 , @xmath259 if we further assume that @xmath2 is directly proportional to the dissipative scale @xmath240 , we arrive at @xmath260 which is not exceedingly promising for use as a les .
all such estimates are , however , purely conjectural until the proper choices of @xmath2 and @xmath258 are determined .
we derived the krmn - howarth equations for the leray@xmath0 and clark@xmath0 models .
these two models may be viewed as successive truncations of the sub - filter scale stress of the lagrangian - averaged navier - stokes @xmath14model ( lans@xmath0 ) . in the case of clark@xmath0 two different inertial range scalings
follow from the dimensional analysis of this equation .
the case of leray@xmath0 is simpler as a single scaling is predicted .
this is the case for navier - stokes and lans@xmath0 as well .
solutions for a @xmath28 dns under periodic boundary conditions ( @xmath261 ) using a @xmath32 resolution under the same exact conditions for lans@xmath0 , leray@xmath0 , clark@xmath0 , and an under - resolved @xmath32 solution of the navier - stokes equations .
we employed two different filter widths @xmath2 .
the first choice @xmath262 was used to and the second choice @xmath263 was employed to test the sgs potential of the models . in comparing these two choices , we found for leray@xmath0 that an increase in @xmath2 substantially decreases the nonlinearity ( and hence decreases the effective reynolds number @xmath4 ) .
for this reason , we were unable to confirm either the inertial range scaling from its krmn - howarth equation or its sub - filter scale energy spectrum .
for clark@xmath0 we were able to determine the dominant krmn - howarth inertial range scaling to be @xmath264 associated @xmath1 energy spectrum ,    the performance of the three regularizations as sgs models ( for a resolution of @xmath32 and @xmath225 ) was comparable to that of the under - resolved navier - stokes solution in reproducing the dns energy spectrum .
only clark@xmath0 showed a clear improvement in approximating the spectrum . from 3d volume rendering of enstrophy density we found that clark@xmath0 and lans@xmath0 were comparable to the under - resolved solution .
even at @xmath265 , leray@xmath0 s 3d spatial structures are consistent with a significantly reduced @xmath4 flow ( e.g. , comparable to a @xmath266 dns ) .
we note that the value of @xmath2 was chosen optimally for lans@xmath0 at the resolution of @xmath32 , and that for clark@xmath0 ( and especially for leray@xmath0 ) smaller resolutions may have comparable results for this value .
such a comparison is beyond the scope of the present work .
although lans@xmath0 and clark@xmath0 exhibit the same inertial range scaling arising from similarities in their krmn - howarth equations , clark@xmath0 is decidedly more intermittent than at the same time , lans@xmath0 is only slightly more intermittent than navier - stokes .
these results are consistent with the artificial truncation of local nonlinear interactions ( in spectral space ) in the sgs stress tensor of each model .
this effect is reduced for lans@xmath0 by the `` rigid - body regions '' enslaved in its larger scale flow which possess no internal degrees of freedom .
the reduced intermittency observed for leray@xmath0 is related to its smoother , more laminar fields as a result of its reduced effective @xmath4 .
finally , we analyzed the reduction in the number of _ dof _ in the models , as compared to navier - stokes ( and , hence , their les potential based on their computational savings ) .
we noted that as lans@xmath0 reproduces the intermittency properties of a dns quite well even for larger values of @xmath2 , some further reduction in numerical saving might be achieved provided the contamination due to its @xmath169 rigid - body energy spectrum were not important in a given application .
as clark@xmath0 possesses a similar reduction in _ dof _ to lans@xmath0 , its les potential is tied to the optimal value of @xmath2 for les .
our study indicates that clark@xmath0 may be applicable ( especially with regards to the energy spectrum ) for larger values of @xmath2 than lans@xmath0 .
in fact , if its optimal value is not a function of @xmath4 , the computational resolution savings increases as @xmath250 for clark@xmath0 . for the case of leray@xmath0
, the prediction is complicated by the effective reduction in @xmath4 as @xmath2 increases .
prediction of optimized values of @xmath2 and of effective dissipation @xmath258 are required to assess its les potential .
future work should include such a study for both leray@xmath0 and clark@xmath0 .
all three regularizations were shown to be successful , in that their control of the flow gradient reduces the degrees of freedom and saves computation while preserving a properly defined reynolds number ( albeit for leray@xmath0 that definition is not yet ) .
clark@xmath0 accurately reproduces the total dissipation , the time scale to obtain a turbulent statistical steady - state and the large - scale energy spectrum of a dns .
these results seem to result from clark@xmath0 being an order @xmath195 approximation of navier - stokes .
we have shown that leray@xmath0 reduces the effective reynolds number of the flow .
the last of the three models , lans@xmath0 restores kelvin s circulation theorem ( advected by a smoothed velocity ) and the conservation of a form of helicity . using spectra as a measure of the success of a subgrid model
, lans@xmath0 is less than optimal , due to its however , other measures of the success of a subgrid model are possible : for example , in regard to intermittency , lans@xmath0 may be considered a superior model . for clark@xmath0 , intermittency may be a function of filter width while for lans@xmath0 , intermittency does not vary much with @xmath2 .
we remark that the computational gain thus far achieved by any of these regularizations is insufficient for applications at very high reynolds numbers , and the three subgrid stress tensors discussed here may need to be supplemented with an enhanced effective viscosity .
this is a common practice when implementing the clark model ( see e.g. , ref .
@xcite for a study of this model with an extra smagorinsky term ) . in this light
, the present study may be useful as an analysis of the properties of the sgs tensors of the regularizations , and to pick best candidates , before the addition of enhanced dissipation .
studies similar to that in ref .
@xcite will also need to be done for the cases of clark@xmath0 and leray@xmath0 to quantify their computational savings before the addition of such dissipative terms .
computer time provided at ncar is gratefully acknowledged , as well as nsf - cmg grant 0327888 .
the work of ddh was partially supported by the us department of energy , office of science , advanced scientific computing research , and by the royal society of london wolfson award for meritorious research .
a.  leonard . .
in _ turbulent diffusion in environmental pollution ; proceedings of the second symposium , charlottesville , va .
, april 8 - 14 , 1973 .
volume a. ( a75 - 30951 13 - 47 ) new york , academic press , inc . , 1974 ,
p. 237 - 248 .
nasa - supported research .
_ , pages 237248 , 1974 .
j.  pietarila graham , d.  holm , p.  mininni , and a.  pouquet .
highly turbulent solutions of the lagrangian - averaged navier - stokes alpha model and their large - eddy - simulation potential .
, 76:056310+ , 2007 .


the low energy processes with a hadron participation are permanently in the center of theoretical and experimental activity . new measurements and more precise experimental data force us to look for new theoretical approaches .
there are two known and widely used methods of effective lagrangian s approach in hadron physics : the effective lagrangian deriving from the qcd principles immediately @xcite-@xcite , and the use of various dynamical symmetries to build the basic lagrangian in a phenomenological way @xcite-@xcite .
the linear sigma - model @xmath1 is the most popular and examined part of the second way .
it is the main approach to the effective analysis of nucleon - nucleon @xcite-@xcite and quark - meson interaction @xcite together with the vector dominance model @xcite-@xcite .
an exact theory of the meson interactions should be nonperturbative .
that is why a lot of hadron characteristics ( effective couplings , formfactors etc . )
are calculated , in particular , in the qcd sum rules method operating with some nonperturbative vacuum parameters . and
that is precisely why we need in effective meson interaction theory giving an information on the low - energy effects and hadron structure .    an effective model of @xmath2
chiral interactions can be deduced from the fundamental qcd @xcite , using the bosonization procedure .
in such a way , the gauge structure of meson interactions follows from the quark level qft @xcite .
thus , the initial quarks with current masses transform to the constituent ones with effective masses @xmath3 , including some contribution of the internal gluon substructure .
therewith , quark loops simulate an uprising of nonperturbative model parameter .
so , the proportionality of pion mass to the `` quark condensate '' is restored in the framework of such effective models  @xcite .
further , to incorporate electromagnetic and strong ( vectorlike ) interactions into the effective theory an appropriate vector fields are introduced as the gauge fields .
the vector meson treatment as the gauge fields realizing some dynamical symmetry , is useful to diminish some theoretical uncertainty of phenomenological description of the hadron interactions . in refs .
@xcite-@xcite the gauge models were successfully used for the consideration of some low - energy aspects of baryon - meson interaction .
it means that the fundamental quantum field principles can be applied to the describing of interaction at the different hierarchical levels . in the case under consideration
, it provides a transition from the quark - level sigma model ( @xmath4 ) to the nucleon - level sigma model ( @xmath5 ) .    to construct an effective meson interaction model
, we will use the gauge scheme based on @xmath0 group .
this group is the simplest one to consider the light unflavored vector mesons , @xmath6 and @xmath7 , together with the photon , as the gauge fields .
moreover , it has the necessary symmetry to describe the electromagnetic and strong meson - meson and quark - meson interactions , which are insensitive to the chiral structure .
so , to analyze the electromagnetic and strong effects only , it is sufficient to localize the diagonal sum of the global chiral @xmath8 .
this sum corresponds to @xmath9 subgroup of the total gauge group @xmath0 .
the extra @xmath10 groups are introduced for the gauge realization of the vector meson dominance ( see the next section ) .
note that in the model we deal with the constituent quarks not with the fundamental current ones .
obviously , the electromagnetic interaction violates the isospin symmetry .
after the spontaneous breaking of the symmetry , higgs fields can be associated with the scalar mesons @xmath11 and @xmath12 ( see the next section ) .    as it was noted above
, this model reproduces the relation @xmath13 , if the quark condensate is taken into account in the equation for the vacuum shift .
further , quarks arise in the model as the gauge group representations , and these degrees of freedom account for the internal meson structure .
in fact , quark - meson models ( see also @xcite ) are necessary hybrid approaches to reproduce both nonperturbative and structure effects of the low - energy meson interactions .    in this paper ,
the gauge scheme interactions of @xmath7 , @xmath6 and @xmath14 with quarks and mesons are analyzed .
in particular , this scheme describes some important features of radiative decays of vector mesons .
namely , the electromagnetic quanta mixing with @xmath7- and @xmath6- fields leads to the gauge variant of vector dominance , i.e. all these interactions are considered in the same way .
so , we deal with the @xmath15 which incorporates gauge interactions of vector mesons with the constituent quarks
. then , the model accounts for hadronic internal degrees of freedom , which are described in the quantum field approach .
really , two - level ( quark - meson ) structure of the model realizes some analog of the `` bootstrap '' idea .
the model contains a number of free parameters in the gauge sector , which can be fixed from the well measured two - particle decays of vector mesons .
as it will be shown later , the gauge model predictions for the radiative vector meson decays are strict ( without any free parameters ) .    as a rule
, the radiative decays @xmath16 and @xmath17 were described in a phenomenological way . within the framework of the gauge model ,
these processes arise at the tree level .
there is a good agreement between the theory and the experimental data on the vector meson radiative decays .
so , these decays can be successfully studied in the gauge field approach .
internal meson structure is described by the quark - meson sector of the model .
these aspects of the model have been tested in @xmath18 , @xmath19 and @xmath20 decays and results of calculations are in a good agreement with the experimental data .
analysis of the decay @xmath21 @xcite , which arises at the loop level due to the quark - meson interaction only , is more cumbersome .
calculations done show that the three - pion decays of @xmath22 and @xmath23 mesons are well descripted in this model .
the obtained results are in preparation now .
the structure of the paper is as follows . in section 2 the gauge field model
is described .
numerical results for tree decays @xmath24 and @xmath25 are presented in section 3 .
we consider loop decays @xmath26 , @xmath27 and @xmath20 in section 4 .
the gauge model of the low - energy meson interactions is based on the conceptual assumption concerning the transfer of gauge principles from the fundamental quark level to the effective hadron one .
just the same transfer was deduced in ref .
@xcite for the case of quark - meson hierarchy levels .
the applicability of the gauge approach to the baryon - meson interaction was demonstrated in refs .
@xcite-@xcite . in this work ,
we show the validity of the gauge scheme for the case of meson - meson and quark - meson effective interactions .
we describe @xmath7- and @xmath6- mesons together with @xmath14- field within the framework of the gauge approach .
thus , the strong and electromagnetic interactions are united in the gauge scheme .
the simplest variant of corresponding dynamical symmetry is based on @xmath28 group .
as it must be , the global symmetry of the chiral model takes place here , and , in a general case , the total chiral group of symmetry should be localized .
however , for the application of the gauge approach to the @xmath23 and @xmath22 strong and electromagnetic decays , it is reasonably to localize the diagonal sum of @xmath29 subgroups together with the additional @xmath10 groups .
the traditional @xmath30- model is a part of the model involving nonlinear terms of self - action .
the @xmath30-model symmetry can be realized in various representations of the global groups : @xmath31 , where @xmath9 means the diagonal sum of @xmath29 subgroups .
for all these cases we have six independent parameters of transformations .
transformation of the quark doublet is : @xmath32 where @xmath33 are the parameters of the @xmath34 group and @xmath35 are the parameters of the @xmath36 group .
thus , here we use @xmath37 representation for @xmath38 and @xmath30 fields .
then , the corresponding transformation properties fo @xmath39 are : @xmath40    having two sets of independent parameters @xmath41 and @xmath35 , to describe strong and electromagnetic interactions it is sufficient to localize the @xmath41 set together with the parameters of @xmath10 groups .
( to consider weak interactions the localization of the @xmath35 parameters should be added . ) as it is seen , the triplet of pions is the adjoined representation of the gauge group . at the same time , the quark doublet @xmath42 is the fundamental representation of this group .
thus , we consider the triplet @xmath43 , the singlet @xmath30-meson and @xmath44-quarks as fields of matter .    the physical scalar fields , emerging from the initial higgs multiplets ,
can be associated with the scalar mesons @xmath11 and @xmath12 .    from these considerations
the model lagrangian is : @xmath45 here @xmath46 - is the first generation quark doublet ; @xmath47 - two scalar fields doublets with hypercharges @xmath48 , @xmath49 and @xmath50 .
the gauge derivatives and field strengthes are : @xmath51 in principle , the gauge invariant term @xmath52 can be added to the lagrangian .
however , this term can be diagonalized from the very beginning by the orthogonal transformation of the initial @xmath53 and @xmath54 fields .
so , we get the same lagrangian with the redefined @xmath53 and @xmath54 fields .    to describe the realistic processes of decays and scattering we introduce also the interactions of leptons with @xmath55 and @xmath10 fields in the gauge form : @xmath56 where interaction of @xmath54 field with leptons is driven by the phenomenological coefficient @xmath57 .
note that the origin of the last term in the above formula is specific .
namely , this term is caused by the nondiagonal @xmath52 part of the initial lagrangian .
thus , the interaction of leptons with the @xmath54 vector field occurs allowing to describe the decays @xmath58 ( see the section 4 ) .
the introducing of the vector fields to the theory in a gauge way provides universality of coupling with vector fields , i.e. it strictly bounds the number of free parameters .
so , this approach raises the predictability of the model .
the analogous universality of vector and pseudoscalar meson interactions was analyzed in @xcite in a phenomenological way .
physical states are formed by the primary fields mixing when mass quadratic forms of scalar and vector fields are diagonalized . at the tree level
the mass forms arise as a result of vacuum shifts : @xmath59    note , the loop contributions to the mass matrix can play an essential role in the mixing of mesons .
so , in a general case this matrix has the form : @xmath60 where @xmath61 is formed by the shifts and @xmath62 by the self - energy insertions . in our case
these insertions are significant for the description of the @xmath63 mixing .    as the result of the diagonalization of vector and scalar mass forms , we get the model spectrum of the vector and scalar multiplets .
the gauge status of vector mesons @xmath64 and @xmath65 is confirmed by the calculation of their model mass spectrum and decay properties .
these model parameters are in agreement with the experimental data for dominant decay channels of vector mesons  @xcite .
moreover , due to the presence of free parameters in the scalar sector , it is possible to describe the mass spectrum and decay channels of the scalar mesons @xmath30 ( or @xmath66 ) , @xmath11 and @xmath12 , which have the status of higgs scalars in the model ( see also eqs.([e:14 ] ) and comment ) .
an emergence of primary hadron component in the photon is an important consequence of the mixing in the gauge sector .
the structure of the vector boson physical states can be illustrated in a tree approximation ( when @xmath67 ) by the following expressions : @xmath68 where the mixing angle @xmath69 is determined under the diagonalization of the vector fields quadratic form .
due to the mixing ( [ e:11 ] ) the processes with initial `` unphysical '' photon , @xmath70 , give contributions to the processes with the intermediate hadron states : @xmath71    some parameters of the mixing can be fixed from the experimental data on the vector meson decays .
it will allows to describe some decays properties of vector mesons considered in this paper .
the tree form of the mixing ( [ e:11 ] ) is caused by the diagonalization of the real mass matrix ( this matrix is generated by the vacuum shifts only , without an account of any self - energy insertions ) with the help of the real orthogonal matrix .
together with the tree real gauge couplings , it leads to the absence of the relative phase in the @xmath72 and @xmath73 channels at the tree level ( see ( [ e:12a ] ) ) .
experiments , however , indicate that there is a nonzero relative phase shift between amplitudes of @xmath74 and @xmath75 processes ( it is interpreted as @xmath76 mixing effect in the pion formfactor ) .    in a phenomenological approach
, the complexity is introduced into the superposition of pure isospin @xmath77 and @xmath78 states .
so , the needed relative phase is a free parameter which is determined from the fit to the experiment  @xcite .
note , such important feature of vector mesons physics as the @xmath79 mixing , can be reproduced in the model at the loop level . as a rule ,
physical vector states of @xmath23 and @xmath22 mesons are described as linear combinations of pure isospin states with complex coefficients .
it results to the relative @xmath80 mixing phase . in our model the complexity
can not be directly introduced into the superposition of operator fields ( [ e:11 ] ) , because it describes the physical neutral fields . to reproduce the relative phase
, we would treat the complexity in a following way .
namely , instead of tree @xmath81 and @xmath82 we should have absolute values of the renormalized parameters of mixing in the superposition ( [ e:11 ] ) . then
, the corresponding phase should be introduced into the renormalized couplings as the vertex factors ( see ( [ e:12a ] ) and ( [ e:12b ] ) ) .
this treatment is equivalent to the known approach when physical meson states are the complex superpositions of pure isospin vector states  @xcite .
the necessary phase value can be provided in the framework of this model due to presence of free parameters , in particular , in the scalar sector ( @xmath83 couplings ) .
certainly , these free parameters should be fixed simultaneously from the analysis of the @xmath84- decays .
for the processes , which are under consideration in this paper , the relative @xmath85 phase is inessential , so we do not use it here .
an analysis of the vector mesons mass matrix with an account of loop contributions and the study of the @xmath86 mixing in the model will be the subject of a forthcoming paper .
now , we give the main part of the physical lagrangian which will be used for calculations : @xmath87 in eqs.([e:12a ] ) : @xmath88 where @xmath89 and @xmath90 are the complex ( renormalized ) parameters of the @xmath91 mixing with account of the self - energy insertion to the mass matrix .
there are also another parts of the interaction lagrangian , which have the structures @xmath92 and @xmath93 . here
@xmath94 and @xmath95 are the sets of scalar and vector physical fields .
the used identification of the gauge and higgs fields components provides the invariance of the physical lagrangian under the @xmath96 transformations .
some relations arise as a direct consequences of the model : @xmath97    the value of @xmath98 can be found from the second relation @xmath99 .
values of @xmath100 and @xmath101 ( or @xmath102 ) are fixed from the experimental values of @xmath103 and @xmath104 .
it should be noted , the last decay takes place due to the mixing ( [ e:11 ] ) only .
thus , we extract the values of main model parameters which we use to describe another decay properties of mesons : @xmath105 these values were applied for calculation of the vector meson radiative decay widths to verify the gauge vector dominance approach . in our strategy of calculations
the strong couplings are extracted from the above mentioned processes as the effective final values .
so , we do not need in loop corrections to these couplings . at the same time
, electromagnetic vertices should be renormalized by the strong interactions ( for details , see the decay @xmath20 which is considered in the fourth section ) .
in addition to ( [ e:12a ] ) we give also a part of the lagrangian , describing an interaction of the scalar mesons with @xmath38-mesons : @xmath106 this part of the lagrangian contains free coupling constants and makes it possible to describe the dominant decay channels of scalar mesons - @xmath107 and @xmath108  @xcite .
the decay channel @xmath109 is not observed and from the model coupling @xmath110 it follows that @xmath111 .
so , in the scalar sector the residual global @xmath9- symmetry approximately takes place after the shift .
exact symmetry , @xmath112 , is not possible because it forbids the decay channel @xmath113 , which is observable experimentally .
evidently , the inequality @xmath114 , violating the global isotopic @xmath9 symmetry , is related with an account of the electromagnetic interaction in the model .
there are no two - photon decays of scalar mesons at the tree level , i.e. these channels are suppressed in the frame of the model .
thus , the status of scalars as the higgs mesons does not contradict to their observed mass spectra and decay properties .
because of the model two - level structure corresponding to the meson - meson and quark - meson interactions , both of them are tested independently in two - particle and three - particle channels .
it is seen from ( [ e:12a ] ) that the meson - meson lagrangian describes both tree and loop level processes , while the quark - meson interactions occur in the model at the loop level only .
radiative decays of neutral vector mesons of the type @xmath115 are an object of steady attention during some decades ( see @xcite and also @xcite-@xcite )
. experimental investigation of these reactions and accompanied theoretical speculations contribute to the understanding of hadron intermediate states and low energy dynamic of meson interactions .    in the channels where the charged pair @xmath116 is formed , dominant contribution comes from tree diagrams corresponding to the vector dominance approximation @xcite .
radiative decays @xmath117 and @xmath118 are described by meson - meson sector of the gauge model .    at the tree level
, the former decay is represented by the diagrams in fig.1 .
the total amplitude for the process is : @xmath119.\ ] ] here @xmath120 and @xmath121 are polarization vectors for @xmath122 - meson and photon , @xmath123 are 4-momenta for all particles in the process . in ( [ e:1 ] )
we omit all terms which are equal to zero on the mass shell in the transversal gauge . for comparison with the experimental spectrum of photons ( see the work of dolinsky , @xcite ) the differential width is presented in the form : @xmath124 where : @xmath125.\end{aligned}\ ] ]    our numerical results , which follow from ( [ e:2 ] ) and ( [ e:3 ] ) , agree with the results of  @xcite given by the vector dominance approach . in fig.2 the theoretical spectrum of photons in comparison with the experimental data ( from  @xcite ) are represented .
the curve ( 2 ) in fig.2 describes the spectrum @xmath126 , normalized by the total width . here
@xmath127 , and it does not depend on the model coupling constants .
the curve ( 1 ) represents the model fit of the experimental data with the help of a single free multiplicative parameter .
as it is seen , this fitted curve ( 1 ) improves the prediction at the low energy range .
this effect can be caused by the renormalization of @xmath128 coupling .    to improve the theoretical spectrum near @xmath129 it was suggested to consider the loop corrections ( see @xcite and @xcite ) .
however , this consideration is reasonable only if the more detailed and reliable experimental data on the photon spectrum are available .
integration of ( [ e:2 ] ) from @xmath130 up to @xmath131 gives the value of partial @xmath6- meson branching @xmath132 which slightly exceeds the experimental value @xmath133 from  @xcite .
an account of the loop contributions ( with the phenomenological couplings ) leads to the result : @xmath134  @xcite .
it should be noted that there is some discrepancy between the experimental branching and the one following from the integration of spectrum ( see fig.2 ) .
namely , an excess of @xmath135 over @xmath136 can be caused by the deviation of the theoretical spectrum from the experimental one at the energy range @xmath137 ( to the left of the first point of the experimental data in fig.2 ) .
decay characteristics of the process @xmath138 at the tree level are computed analogously with the following replacement in ( [ e:1 ] ) - ( [ e:3 ] ) : @xmath139 in @xmath140 , and @xmath141 in @xmath142 and @xmath143 .
the partial width for the decay is damped by the small mixing parameter , @xmath144 , so we have @xmath145 and @xmath146 for @xmath147 and @xmath148 , respectively .
these estimations do not contradict to the experimental restriction @xmath149  @xcite and agree with the theoretical results of  @xcite .
thus , the same set of the fixed parameters ( the gauge constant and two mixing tree angles ) describes two different decays at the tree level .
certainly , loop corrections can be important for the case due to the smallness of the tree contribution , see , for example  @xcite . in analogy with the results of these papers
, loop corrections can increase @xmath150 up to @xmath151 , which does not contradict again to the upper limit of the experiment .
as it was noted above , consideration of the loop corrections to the decay @xmath152 in @xcite-@xcite was intended to describe the photon spectrum fine structure .
however , it increases the discrepancy between the model and experimental values of the total width .
quantitatively , this effect depends on the underlying model .
moreover , in effective theories an account of loop diagrams ( for some process arising at the tree level ) has some subtleties connected with the compensation of divergencies and renormalizability . for the processes which occur at the loop level only , these problems are absent  all divergencies are summed to zero when all external lines are on the mass shell  @xcite .
an examples of such loop processes will be given in the next section .
the usage of parameters from ( [ e:13 ] ) together with the constants @xmath153 and @xmath154 extracted from the widths of @xmath155 and @xmath156 decays , leads to the correct estimation of the cross section in the peaks vicinity , @xmath157 .    as it was noted above , the tree approximation ( without an account of the @xmath76 mixing at the loop level ) does not allow us to describe the resonance curve in this region in details .
a relative @xmath85 phase consideration is a very essential for the resonance curve description .
in the most papers these decays are defined by the phenomenological vertices which are introduced at the tree level @xcite . in the model considered , radiative decays @xmath161 and three - particle decays @xmath162 occur via quark loops with the gauge vertices .
one - loop diagrams for the decays @xmath163 are shown in fig.3 , where @xmath164 denote quark fields .
total amplitude for the process @xmath165 has the following form @xmath166 here @xmath167 is the color factor , vertex constant @xmath168 ( @xmath169 is the constant of @xmath170-interaction , see ( [ e:12a ] ) ) and @xmath171 is the three - point passarino - veltman function @xcite . for the constituent quark mass
we suppose @xmath172 . substituting the function @xmath173 into ( [ e:4 ] )
, we have : @xmath174 where @xmath175 and the decay constant @xmath176 .
in ( [ e:5 ] ) it was used the goldberger - treiman relation @xmath177 ( see , for example , @xcite ) .
it defines the coupling @xmath169 through the constituent quark mass .
then the constituent quark mass value , which is taken as a free parameter , can be found from the widths fit .
the decay @xmath178 is described by the diagrams in fig.3 with the corresponding replacement of the coupling constant .
the expression for the width is @xmath179 follows from the isotopic structure of @xmath180 , @xmath181 and @xmath182 vertices ( [ e:12a])-([e:12b ] ) ( here we omit small terms proportional to @xmath183 and @xmath89 ) .
both the widths are in a good agreement with the experimental data @xcite when the effective quark mass value is @xmath184 : @xmath185 thus , the fitting of two widths by one model parameter - mass of the constituent quark @xmath186 , gives @xmath187 .
this quark mass value can be understood in a connection with the schematic representation of the quark - gluon content of hadrons .
specifically , for the nucleon and @xmath6- meson masses we can write approximately @xmath188    here we use an effective quark mass , @xmath186 , and suppose the same gluon component of meson and baryon structure , @xmath189
. from this system of equations it follows @xmath190 and @xmath191 .
the latest value nearly coincides with the quark mass assessment resulting from the decay channels analysis .
note that a sufficiently low quark mass value @xmath192 together with the @xmath193 are followed from the composite - meson model with the four - quark interaction  @xcite .
so , the model discussed introduces the constituent quark which is differed from the traditional effective quark component of a hadron with @xmath194 ( see , however ,  @xcite ) .
moreover , the value of @xmath189 makes it possible to interpret the meson @xmath195 as the ( mostly ) glueball state  @xcite .
it does not contradict to some known ideas on the gluon nature of @xmath196 meson ( see  @xcite and references therein for the review of scalar meson properties ) .
if it is the case , the low value of the constituent quark mass can be treated as the consequence of a model separation of the quark and gluon degrees of freedom in a hadron .
hence , the value of @xmath186 , which is essentially lower than @xmath197 , provides the agreement of the model results with the experiment for the processes under consideration .    of course
, if we explanate the @xmath30-meson as the effective glueball state , we need in an detailed analysis of its properties and manifestations in the model .
this studying will be presented in a forthcoming paper .
note also that the @xmath196 meson was carefully analyzed , supposing the meson is a @xmath198 state , in  @xcite , for instance .    in an analogy with the process @xmath199 ,
the decays under consideration can contain an anomaly contribution caused by the @xmath200 vertex .
this contribution is small here due to the smallness of the ratio @xmath201 , that is approximation @xmath202 takes place with a good accuracy .
it is known that the same approximation is not admissible for the decay @xmath203 , where the anomaly contribution dominates .
an analysis of this decay in our model is coincide with the one in  @xcite , for example .
namely , due to the dominant anomaly contribution , the decay amplitude is again @xmath204 i.e. it depends on the @xmath205 and does not contain the quark mass .
the constant @xmath205 arises in the amplitude from pion - quark coupling which is @xmath206 , as it dictates by the goldberger - treiman relation .
so , our value of the quark mass corresponds to the ( lower ) pion - quark coupling , @xmath169 , providing , however , the reasonable values for @xmath207 widths .
note also that in the consideration of loop three - particle decays @xmath208 we will have in this model a non - negligible ratio @xmath209 as a consequence of this relatively low quark mass ( cf . with  @xcite ) .
nevertheless , an accurate calculations give for @xmath210 the values which are in agreement with the data ( these results will be discussed in the forthcoming work too ) .
one of the most important features of the model is the implantation of electromagnetic ( em ) interaction into the gauge scheme of strong interactions .
the mixing of the gauge fields leads to the effect of mutual renormalization of the interaction couplings .
now we discuss the decay @xmath20 which is caused mainly by the em interaction . at the tree level this process is described by the vertex following from ( [ e:10 ] ) and the mixing ( [ e:11 ] ) .
corresponding decay coupling is @xmath211 , so it gives @xmath212 , while the experimental decay rate is @xmath213 @xcite .
here we show that this discrepancy can be eliminated due to renormalization of the decay coupling @xmath153 by the strong interactions . in the one - loop approximation
there are two types of diagrams which give the main contribution to the renormalization ( see fig.4 ) . in this figure
, solid lines correspond to @xmath214 and @xmath215 ( or all set of scalar particles , @xmath216 ) , and wavy lines correspond to @xmath14 and @xmath122 .
all extra diagrams are damped by the small mixing parameters .    as a result
, we get full amplitude of the process : @xmath217,\ ] ] where @xmath218 and @xmath219 is the sum of the loop contributions . finally , the decay rate has the form : @xmath220 where @xmath221 and @xmath143 is the renormalization parameter .
the function @xmath219 in the @xmath222 scheme is : @xmath223 here we use normalized quantities : @xmath224 in the ms scheme an additional term arises in ( [ e:18 ] ) : @xmath225 where @xmath226 . with the help of the expressions ( [ e:16])-([e:20 ] )
we have found that @xmath227 when the renormalization parameter @xmath228 in the @xmath222 scheme and @xmath229 in the ms scheme .
so , strong interactions lead to the large renormalization of the em coupling and the physical result @xmath230 is obtained within the framework of @xmath222 scheme when the renormalization scale @xmath143 is in the range of typical meson masses .
the value of @xmath231 depends on the masses of quarks and scalar mesons very weakly .
the quark mass changing from @xmath232 to @xmath233 leads to the variation of @xmath143 lower than @xmath234 .
variation of @xmath143 due to the increasing of scalar mass from @xmath235 up to @xmath236 is approximately @xmath237 .
an analogous effect takes place for the decay @xmath238 , but in this case we have use free parameter ( @xmath57 from ( [ e:10 ] ) ) .
the presence of this parameter makes it possible to get the accordance @xmath239 for the same values of @xmath143 .
the linear sigma - model is the most popular approach in the description of low - energy hadron interaction .
we have considered the gauge generalization of this model and include quark degrees of freedom explicitly . in the gauge scheme
the vector meson dominance occurs in the tree - level processes .
the quark - meson sector describes some loop - level processes , when quark structure of mesons plays a noticeable role .
the model is applied to some radiative decays of mesons which are intensively discussed in the literature .
decays @xmath24 and @xmath25 have been considered within the framework of the vector dominance gauge variant .
the results are in good accordance with the experimental data .
the radiative decays @xmath158 and @xmath240 take place at the loop level only , when the quark - meson interaction is important . for these cases
we get results in agreement with the experiment too .
decay channel @xmath20 can be described by the model when we take into account the electromagnetic decay constant renormalization by strong interactions .    from the analysis we fulfilled
, it follows that the radiative decay rates of vector mesons can be calculated with a good accuracy within the framework of the gauge model .
quantum field approach to the vector dominance in meson interaction can be also supplemented by the quark - meson interaction which expands the linear sigma - model .
99 y.  nambu and g.  jona - lasinio , _ dynamical model of elementary particles based on an analogy with superconductivity .
* 122 * ( 1961 ) 345 .
y.  nambu and g.  jona - lasinio , _ dynamical model of elementary particles based on an analogy with superconductivity .
ii , phys .
rev . _ * 124 * , ( 1961 ) 246 .
m.  k.  volkov , _ low - energy meson physics in the quark model of superconductivity type , sov .
and nuclei . _
* 17 * ( 1986 ) 186 ; _ effective chiral lagrangians and the nambu - jona - lasinio model , sov
. j. part .
and nuclei . _ * 24 * ( 1993 ) 81 .
m.  k.  volkov , a.  e.  radzhabov , _ the nambu - jona - lasinio model and its development , sov . j.
phys - usp .
_ * 49(6 ) * ( 2006 ) 551 .
g.  v.  efimov , m.  a.  ivanov and v.  e.  lyubovitskiy , _
predictions for semileptonic decay rates of charmed baryons in the qurk confinement model , z.phys .
c - particles and fields _ * 52 * ( 1991 ) 149 [ jinr preprint e2 - 90 - 24 , dubna , 1990 ] . n.  i.  karchev and a.  a.  slavnov , _ effective chiral lagrangians from @xmath241 quantum chromodynamics , teor . mat . fiz .
( ussr ) _ * 65*(1985 ) 192 .
t.  goldman and r.  w.  haymaker , _ dynamically broken chiral symmetry with bag confinement phys
* d24 * ( 1981 ) 724 .
d.  espriu , e. de rafael and j.  taron , _ the qcd effective action at long distances , nucl .
* b345 * ( 1990 ) 22 .
d.  ebert and h.  reinhardt , _ effective chiral hadron lagrangian with anomalies and skyrme terms from qurk flavour dynamics , nucl .
phys . _ * b271 * ( 1986 ) 188 .
hsieh , _ effective meson field theory from qcd _ [ hep - ph/0306143 ]
. s.  s.  olean , _ effective lagrangian for baryons and baryon - meson interaction _ [ hep - ph/0310005 ]
. s.  l.  adler , _ axial vector vertex in spinor electrodynamics , phys .
* 177 * ( 1969 ) 2426 .
j.  s.  bell and r.  jackiw , _ a pcac puzzle : @xmath242 gamma gamma in the sigma model , nuovo cim . _ * 60 * ( 1969 ) 47 .
j.  d.  walecka , _ a theory of highly condensed matter , ann .
phys . _ * 83 * ( 1974 ) 491 .
b.  d.  serot , _ a relativistic nuclear field theory with pi and rho mesons , phys .
b _ * 86 * ( 1979 ) 146 .
k.  sailor , i.  lovas , _ gauge - invariant qhd lagrangian , acta phys . hung .
_ * 64 * ( 1988 ) 123 .
r.  raczka , _ a spontaneously broken quintet baryon model for strong interactions , czech .
b _ * 37 * ( 1987 ) 421 .
g.  m.  vereshkov , v.  i.  kuksa , _
u(1)su(3 ) gauge model of barion - meson interactions , yad . fiz .
_ * 54 * ( 1991 ) 1700 .
j.  c.  polkinghorne , _ renormalization of axial vector coupling , il nuovo cim .
_ * 8 * ( 1958 ) 179 .
j.  schwinger , _ field theory commutators , phys .
lett . _ * 3 * ( 1959 ) 296 . m.  gell - mann , m.  levy , _ the axial vector current in beta decay , il nuovo cim .
_ * 16 * ( 1960 ) 705 .
r.  delbourgo and m.  d.  scadron , _ dynamical generation of the gauged @xmath9 linear sigma model , mod .
lett . _ * a10 * ( 1995 ) 251 [ hep - ph/9910242 ] .
a.  bramon , riazuddin and m.  d.  scadron , _
double counting ambiguities in the linear sigma model , j. phys . _
* c24 * ( 1998 ) 1 [ hep - ph/9709274 ] .
r.  delbourgo , d.  liu and m.  d.  scadron , _ meson pvv interactions are determined by quark loops , int .
* a14 * ( 1999 ) 4331 ( 1999 ) [ hep - ph/9905501 ] .
r.  p.  feynman , _ photon - hadron interactions , w.  a.  benjamin , inc .
_ ( 1972 ) .
a.  bramon , a.  grau and g.  pancheri , _ intermediate vector meson contributions to @xmath243 decays , phys .
_ * b283 * ( 1992 ) 416 ; _ chiral perturbation theory and radiative @xmath243 decays , phys.lett._ * b289 * ( 1992 ) 97 ; _ radiative vector - meson decays in @xmath244 broken effective chiral lagrangians , phys.lett._ * b344 * ( 1995 ) 240 .
d.  black , m.  harada and j.  schechter , _ vector meson dominance model for radiative decays involving light scalar mesons , phys .
* 88 * ( 2002 ) 181603 [ hep - ph/0202069 ] .
baboukhadia , v.  elias , m.d .
scadron , _ linear sigma model linkage with nonperturbative qcd , j. phys . _ * g23 * ( 1997 ) 1065 [ hep - ph/9708431 ] . m.d .
scadron , _ quark level linear @xmath30 model ( l@xmath30 m ) via loop graphs _ [ hep - ph/0007184 ] .
j.  l.  lucio m. and m.  napsuciale , m.  d.  scadron , v.  m.  villanueva , _ the @xmath245 transition and @xmath246 decay , phys .
* d61 * ( 2000 ) 034013 [ hep - ph/9902349 ] .
p.  singer , _ decay mode @xmath247 , phys .
* 128 * ( 1962 ) 2789 ; _ radiative @xmath6-meson decay , phys . rev . _ * 130 * ( 1963 ) 2441 .
m.  gell - mann , d.  sharp and w.  g.  wagner , _ decay rates of neutral mesons , phys .
* 8 * ( 1962 ) 261 .
e.  a.  kuraev , z.  k.  silagadze , _ once more about the @xmath248 contact term , phys .
* 58 * ( 1995 ) 1589 [ hep - ph/9502406 ] .
yao et al . , _ review of particle physics , particle data group j. phys . _ * g33 * ( 2006 ) 1 . f.  klingl , n.  kaiser , w.  weise , _ effective lagrangian approach to vector mesons , their structure and decays@xmath249 , z.phys._ * a356 * ( 1996 ) 193 [ hep - ph/9607431 ] .
r.  urech , _ @xmath85 mixing in chiral perturbation theory , phys.lett._ * b355 * ( 1995 ) 308 [ hep - ph/9504238 ] .
oconnell , b.c .
pearce , a.w .
thomas , a.g .
williams , _ @xmath85 mixing , vector meson dominance and the pion form - factor , prog.nucl.part.phys._ * 39 * ( 1997 ) 201 [ hep - ph/9501251 ] .
yan , j .- h .  jiang , x .- j .
wang , _ on the mixed propagator approach to @xmath76 mixing , commun.theor.phys._ * 38 * ( 2002 ) 195 [ hep - ph/0111281 ] .
m.  benayoun et.al .
_ new results in @xmath122 meson physics , eur.phys.j._ * c2 * ( 1998 ) 269 [ hep - ph/9707509 ] . m.  benayoun et.al .
_ the pion form factor within the hidden local symmetry model , eur.phys.j._ * c29 * ( 2003 ) 397 [ nucl - th/0301037 ] .
a.  bramon , r.  escribano , _ a consistent description of @xmath250 decays including @xmath251- meson effects , jhep _ * 0402 * ( 2004 ) 048 [ hep - ph/0305043 ] .
a.  gokalp , s.  solmaz and o.  yilmaz , _ scalar sigma meson effects in radiative @xmath122-meson decays , phys.rev._ * d67 * ( 2003 ) 073007 [ hep - ph/0302129 ] .
y.  oh and h.  kim , _ higher meson resonances in @xmath252 and @xmath253 , phys
* d68 * ( 2003 ) 094003 [ hep - ph/0307286 ] .
s.  i.  dolinsky et al . , _ summary of experiments with the neutral detector at the @xmath254 storage ring vepp-2 m , phys . rep .
_ * 202 * ( 1991 ) 99 .
a.  gokalp , a.  kucukarslan and o.  yilmaz , _ vmd , chiral loops , @xmath30-meson , and @xmath255 mixing in @xmath253 decay , phys .
* d67 * ( 2003 ) 073008 [ hep - ph/0302240 ] .
a.  gokalp , a.  kucukarslan , s.  solmaz and o.  yilmaz , _ sigma - meson and @xmath255 mixing effects in @xmath256 decay , acta .
polon . _ * b34 * ( 2003 ) 4095 [ hep - ph/0306044 ] .
v.  kuksa , _ compensation of divergences in rare processes , yad . fiz . _
* 66 * ( 2003 ) 1785 .
d.  bardin and g.  passarino , _ standard model in the making , oxford university press _ ( 1999 ) .
d.  ebert , m.k .
volkov , _ composite - meson model with vector dominance based on @xmath257 invariant four - quark interactions , z.phys.c-particles and fields _ * 16 * ( 1983 ) 205
. k.  hagiwara et al . , _ review of particle physics , particle data group , phys .
* d66 * ( 2002 ) 1 . m.d .
scadron , f.  kleefeld , g.  rupp , _ constituent and current quark masses at low chiral energies , eurphys .
* 80 * ( 2007 ) 51001 [ hep - ph/0710.2273 ]
. h.  jin , x.  zhang , _ scalar glueball decay into pions in effective theory , phys .
_ * d66 * ( 2002 ) 057505 [ hep - ph/0208120 ] .
anisovich , _ systematics of @xmath258- states , scalar mesons and glueball , aip .
* 619 * ( 2002 ) 197 [ hep - ph/0110326 ] .
scadron , _ scalar sigma meson via chiral and crossing dynamics , eur .
phys . j. _ * c6 * ( 1999 ) 141 [ hep - ph/9710317 ] .
e.  van beveren , f.  kleefeld , g.  rupp , m.d .
scadron , _
remarks on the @xmath259 scalar meson as the dynamically generated chiral partner of the pion , mod .
phys.lett._ * a17 * ( 2002 ) 1673 [ hep - ph/0204139 ] .