in the past two decades , cosmology has acquired a new standing due to a wide range of observations of high precision data @xcite .
these observations seem to reveal that the two dominant components of the universe are a non - baryonic form of matter , the so - called dark matter @xcite , and an exotic form of energy , dark energy @xcite .
the former leads to structure formation while the latter is responsible for the present accelerated expansion of the universe .
the search for an explanation for the existence and properties of these dark components of the universe has prompted a strong debate about their origin and nature .
the behaviour of the galaxy rotation curves and the dynamical mass in clusters of galaxies suggest the existence of an exotic form of matter , dark matter , at galactic and extra galactic scales , respectively , where the presence of this sort of matter is revealed through its gravitational effects .
notwithstanding , recent results demonstrated that dark matter may interact with ordinary matter , with the production of patterns of annihilations or decays of dark matter particles in the fluxes of cosmic rays @xcite .
evidence for the existence of dark matter lies in the velocity dispersion of clusters that , combined with its morphological features , leads one to conclude that the overall mass should be far greater than the visible mass .
in fact , the dark matter component of these clusters is well described by dark matter models , where the interaction of dark matter with baryonic matter determines the density profile of galaxy clusters @xcite .
an alternative approach to the above mentioned problems is to consider modified theories of gravity . amongst modified gravity theories , those that introduce higher order curvature terms in the einstein
hilbert action  and , in particular , the @xmath0 theories of gravity @xcite  offer alternative explanations for these problems , along with other cosmological , astrophysical , and high - energy physics motivations @xcite .
recently , it has been shown that , in the context of power - law theories , @xmath1 , the rotation curves can be explained as a curvature effect @xcite ; a similar result @xcite can be obtained in another extension of gravity that relies not only on a non - trivial curvature term , but also on a non - minimal coupling between matter and geometry @xcite .
clusters are the largest astronomical structures whose masses can be measured in a reliable form , since they constitute the largest configuration of objects that passed through gravitational relaxation and entered into virial equilibrium @xcite .
the measuring methods , x - rays and gravitational lensing for visible and dark matter respectively , reveal that the masses of clusters are approximately seven to ten times larger than the total combined mass of stars and hot gas in the intracluster medium ( icm ) .
this missing mass is predominant in the inner regions of clusters , but it is also extended beyond the core radius , as defined by the gas density distribution .    in the context of modified newtonian dynamics @xcite ,
another popular ( albeit incomplete @xcite ) alternative to general relativity ( gr ) , it is difficult to explain cluster observations without the assumption of `` real '' dark matter , because the mass of the dark component is otherwise significantly reduced but not completely removed @xcite .    the problem of clusters has recently been discussed in @xmath0 theories , by considering a generalized version of the virial theorem .
previously discussed in the context of a possible unification and interaction of dark matter and dark energy @xcite . assuming a steady state , the virial theorem can be used to deduce the mean density of astrophysical objects , and therefore the corresponding total mass , such as galaxies and clusters , by observing the velocities of test particles moving around them .
the overall mass of a sample of clusters can also be estimated through corrections in the gravitational potential that emerge in the weak field limit of these theories @xcite .
along with these works , there were also conducted the first cluster abundance constraints on a modified gravity model , specifically the modified action @xmath0 model , which predict deviations in the abundance of massive dark matter halos @xcite .    in the context of the non - minimal coupling between the ricci scalar curvature and matter @xcite , it has been shown in several works that both dark energy and dark matter can be effectively described within this framework , along with other cosmological and astrophysical implications : these include the perturbed hydrostatic equilibrium of the sun @xcite , the acceleration of the universe @xcite , the galaxy rotation curves @xcite , the reheating scenario after inflation @xcite , the mimicking of a cosmological constant @xcite , and the change in the gravitational potential @xcite .    in order to extend this non - minimally coupled model one tries to obtain the large mass difference of clusters through a dark matter mimicking , an extra matter component that emerges due to the presence of the non - minimal coupling , and a subsequent modification of gravity .
one first explores here the scenario of the abell cluster a586 , which is a nearby massive strong - lensing cluster , with a high x - ray luminosity , related with the mass of the intracluster gas .
its morphology is well suited to the assumption of spherical symmetry , making it an ideal test bed for modifications of gravity at extra - galactic scales .
one then extends this analysis to a larger set of galaxy clusters with the same features .
this work is organized as follows .
first , one shortly reviews the formalism leading to the total cluster mass and the field equations that result in the presence of the non - minimal coupling and stablishes the strategy to address the problem of mimicking the dark matter component of a cluster .
one then studies the particular case of the abell cluster a586 , discusses the coupling function that leads to the dark matter mimicking , and extends the model to a larger sample of galaxy clusters .
finally , the numerical results are presented and discussed .
galaxy clusters are usually considered as closed gravitational systems with spherical symmetry and in hydrostatic equilibrium when virial equilibrium is attained  despite the fact that recent observations reveal that clusters have more evolved structures with strong interactions and dynamical activity , particularly in their inner regions @xcite .    under these assumptions ,
the structure equation can be derived from the collisionless boltzmann equation for an isotropic system @xcite ,    ( _ g _ r^2)=-_g ddr  , [ boltzmann_equation ]    where @xmath2 is the newtonian gravitational potential of the cluster , @xmath3 the mass - weighted velocity dispersion in the radial direction , and @xmath4 the gas - mass density .
the pressure profile is related with these quantities through , @xmath5 , and inserting in eq .
( [ boltzmann_equation ] ) leads to @xmath6 , where the prime indicates differentiation with respect to the radial coordinate .
the newtonian gravitational potential @xmath2 can be read from the newtonian limit of the modified einstein field equations ( as addressed later ) , _ i.e. _ a suitably modified poisson equation . however
, this does not prevent us from formally writing the latter as the effect of the non - minimal coupling ascribed to the `` mimicked '' dark matter component , so that    ^2 = 4 g ( + _ dm )  . [ poissoneq ]    this approach shall be discussed in more detail in the following section .    for a gas sphere , with temperature profile @xmath7 ,
the velocity dispersion becomes @xmath8 , where @xmath9 is the boltzmann constant , @xmath10 the mean mass particle , and @xmath11 the proton mass . introducing this expression into eq .
( [ boltzmann_equation ] ) gives    ( kt m_p_g ) & = & -_g ddr  ,    or , equivalently ,    -ddr & = & kt m_pr  .
[ eq : boltzmannpotential ]    therefore , from the models for the gas density and temperature profiles , the total mass enclosed in a sphere with a given radius @xmath12 can be obtained from ,    m ( r ) & = & 4(_g + _ dm ) r^2 dr = r^2 g ddr  , [ dynamical_mass_possion ] leading to the usual expression for the mass profile of a spherical mass distribution in hydrostatic equilibrium m ( r ) & = & -kt rgm_p(d_g dr+dt dr )  .
[ dynamical_mass ]    as already stated , this is the total mass that one should observe in a cluster and that , in general , agrees well with observations from gravitational lensing @xcite .
however , the baryonic component of the cluster mass , galaxies mass in the visible spectrum and gas in x - ray spectrum , where the latter dominates , is much less than the total observed mass leading one to conclude the existence of a dark matter component that in the present work results from the presence of the non - minimal coupling .    hence , the equilibrium equation , eq .
( [ dynamical_mass ] ) , which leads to the total mass within a cluster , can be used to derive the amount of dark matter present in a cluster of galaxies , and its spatial distribution , by the mass difference between this and the gas - mass estimates as provided by x - ray observations .    in the following ,
one tries to reproduce the total mass profiles of a sample of galaxy clusters , where the additional density component emerges not from an exotic form of `` dark '' matter but from a change in geometry due to the presence of the non - minimal coupling in the action .
clusters of galaxies present a velocity dispersion with a relative constant value along the radius , which reveals an implicit matter density profile of the form @xmath13 .
the singular isothermal sphere is the simplest model whose density profile , @xmath14 , leads to a constant and isotropic velocity dispersion @xmath15 along the clusters profile .
however , numerical fits to the observational data show that the density profile of clusters are flatter than the isothermal one at small radii and steeper at larger radii . a general form replicating this behaviour
is given by the cusped density profile ,    _ m&=_0 & ( r r_s ) ^-p(1+r r_s)^p - q  , [ cusped_profile ]    where the parameters @xmath16 and @xmath17 determine the respective slopes of the power - law density profile in the inner and outer regions , and @xmath18 is the so - called core radius , signaling the steepening of the profile .
these parameters usually present values within the ranges @xmath19 and @xmath20 , for models such as the navarro - frenk - white ( nfw ) profile ( @xmath21 and @xmath22 ) , the moore profile ( @xmath23 and @xmath22 ) , and the rasia profile ( @xmath21 and @xmath24 ) .
in fact , one of the important results of hierarchical cdm models is that the density profile of cold dark matter halos is described by the nfw model @xcite .
the more relevant baryonic component is the icm in the form of hot gas , which emits observable x - ray radiation . in general , the observed surface - brightness profiles of clusters that reveal the underlying visible matter profile of this intracluster gas are well described by the so - called @xmath25 model :    _ g(r)=_g0(1+r^2 r^2_c)^-3/2  , [ visible_profile ]    where @xmath26 is the core radius , which has a specific value for each cluster @xcite .      in this cluster , the emission measured profile is well described by the @xmath25 model , eq .
( [ visible_profile ] ) , and , from a least - squares fit of two different sets of observational data , the obtained values for the free parameters are @xmath27 and @xmath28 @xcite .
the icm particle number density profile is given directly by the analytic fit for the projected emission measure profile , eq .
( [ visible_profile ] ) , and is easily converted into the gas density .    the deep potential well of the cluster compresses the icm gas to x - ray emitting temperatures , such that the gas temperature can be inferred from the derived spectrum .
for a polytropic temperature model @xcite ,    t(r)&=&t_0 ( 1+r^2 r_c^2)^-3(-1)/2  , [ temperature_profile ]    where @xmath26 and @xmath25 are parameters obtained in the fit for the visible mass density , and a least - squares fit leads to the values @xmath29 , @xmath30 @xcite .
naturally , setting @xmath31 leads to the isothermal temperature profiles .
the total mass follows from the hydrostatic equilibrium equation , eq .
( [ dynamical_mass ] ) , using the previous expressions for the visible density and temperature profiles    m(r)&= & m_0 ( r r_c)^3(1+r^2 r^2_c ) ^-1 - 32(-1 )  , [ mass_profile ]    with    m_0 = 3 k t_0 r_c g m_p  .    from the relation for a spherical mass distribution , = 4r^2 ( r )  ,    one may approximate the dark matter density by the total density profile using eq .
( [ mass_profile ] ) in the last expression and solving for @xmath32 ,    [ dm_profile ] & & _
dm(r ) _ dm0(1+)(1+r^2r_c^2)^-2  .    where    _ dm0&=&3 k t_0 4 g m_p r_c^2 = m_0 4r_c^3 [ defdm0 ]  ,    and one defines @xmath33\be$ ] , for convenience .
finally , one introduces the virial radius @xmath34 , which signals the distance where the density of the cluster approaches a multiple @xmath35 of the critical density of the universe @xmath36 , where @xmath37 is the hubble expansion rate at the present .
thus , from @xmath38 , one may estimate the virial radius from the relation eq .
( [ dm_profile ] ) ,    & & ( _ dm0_c ) ^1 /(2-+ )  ,    which , for @xmath39 , with @xmath40 , leads to a virial radius of @xmath41 .
recently , a generalization of the @xmath0 modified theories of gravity was considered that further extends the presence of curvature invariant terms in the einstein
hilbert action by non - minimally coupling the scalar curvature with matter @xcite .
the action is obtained by adding an additional @xmath42 term , dependent on the scalar curvature @xmath43 and coupled with the matter lagrangian @xmath44 , to the usual metric theories of gravity , leading to @xcite    s = d^4 x  , [ model ]    where @xmath45 ( with @xmath46 ) are arbitrary functions of the scalar curvature and @xmath47 is the metric determinant . setting @xmath48
one obtains the usual @xmath0 theories of gravity , and gr is recovered if one takes the linear function @xmath49 , where @xmath50 .
variation with respect to the metric , @xmath51 , leads to the field equations ,    & & ( f_1 + 2 f_2_m)r_-1 2f_1g_= + & & = ( _ -g_)(f_1 + 2 f_2_m)+(1+f_2)t _  , [ field_equations ] where @xmath52 is the energy - momentum tensor
. the respective trace equation is given by    & & ( f_1 + 2 f_2_m)r-2f_1= + & & = -3(f_1 + 2 f_2_m)+(1+f_2)t  .
[ field_equations_trace ]      from the covariant derivative of the field equations , eqs .
( [ field_equations ] ) , along with the bianchi identities , @xmath53 , one encounters the covariant non - conservation of the energy - momentum tensor @xcite ,    _ t^&= & ( g^_m - t^ ) _ ( 1+f_2 )  , [ non - cons ]    which can be regarded as an energy - momentum exchange between matter and geometry . the usual covariant conservation for the energy - momentum tensor
is naturally obtained when the non - minimal coupling vanishes , @xmath48 .
this non - conservation implies that a massive particle in the absence of forces will not describe a geodesic curve .
strong variations on the extra term in the last equation can imply a violation of the equivalence principle , leading to a possible way of testing and setting bounds on the coupling functions , @xmath42 .
contrary to the usual jordan  brans  dicke theories , it is impossible to perform a conformal transformation to the einstein frame such that the coupling disappears and the covariant conservation of the energy - momentum tensor is recovered for all matter forms @xcite .
indeed ( following ref .
@xcite ) , using the conformal transformation    g__= f_2 g _  ,    one obtains the covariant conservation law    _
^= 0  ,    only if @xmath54 and @xmath55 . from the energy - momentum tensor of a perfect fluid , @xmath56 , where @xmath57 is the energy density , @xmath16 the pressure and @xmath58 the four - velocity ( with @xmath59 ) , the respective trace is @xmath60 , such that one may write ,    2= 3p - p=(2+)/3  .
however , one obtains different equations of state depending on the chosen lagrangian density @xmath61 .
for instance , if @xmath62 @xcite , one has the state equation @xmath63 , a perfect fluid with negative pressure , even though not a cosmological constant for which @xmath64 ; for @xmath65 , then @xmath66 , which is the equation of state for ultra - stiff matter .
therefore , this non - conservation is , in fact , a fundamental property of the model , eq .
( [ model ] ) , meaning that even under a suitable conformal transformation the energy
momentum tensor is not conserved for all different types of matter .    assuming that matter is described by a perfect fluid , the non - conservation , eq .
( [ non - cons ] ) , results in an extra force of the form    f^&=&1 + p h^ , [ extra - force ]    which has units of an acceleration , and where the projection operator is given by @xmath67 , such that the extra force is orthogonal to the four - velocity , @xmath68 .      in order to test the effect of the non - minimal coupling ,
one focuses on possible deviations from gr due to the presence of the coupling by considering the linear function @xmath69 ; following ref .
@xcite , one adopts a power - law non - minimal coupling    f_2(r)= ( r r_n ) ^n  , [ powerlawcoupling ]    where @xmath70 is a characteristic curvature .    considering a pressureless perfect fluid , @xmath71 , and the respective trace equation , @xmath72 , along with the lagrangian density , @xmath73 ( a choice discussed in depth in ref .
@xcite ) , the field eqs .
( [ field_equations ] ) become    & & r_- 1 2rg_= + & & = n ( g_- _ ) + & & + 1 2u_u _  , [ first_equation ]    and the respective trace , eq .
( [ field_equations_trace ] ) , reads    r=1 2 - 3n  .
[ second_equation ]    if one considers a strong coupling @xmath74 , the above has an implicit solution given by the vanishing of the derivative term in the _ r.h.s . _
,    = _ dm _
n ^1/(1-n )  , [ implicit ]    where one introduces the characteristic density , @xmath75 .
one dubs this a `` static '' solution , not because it is time - independent ( which it is , if one assumes a time  independent visible matter density @xmath32 ) , but in opposition to a more evolved , `` dynamical '' solution that would arise if this second term did not vanish .
indeed , it was shown that ( see ref .
@xcite ) the most general solution of the differential equation corresponds to a dominance of the gradient terms on the _ r.h.s .
_ of this equation : however , it was found that this `` dynamical '' solution exhibits very small oscillations around the `` static '' one ( with a very short - wavelength providing the dominance of the gradient term ) .
thus , the scaling law @xmath76 remains valid even if the gradient term dominates the dynamics .
furthermore , it was shown that , although the above analytical treatment remains valid , this implies that the non - minimal coupling is actually perturbative , _
i.e. _ @xmath77
. this will be used in the following section to avoid the destabilizing effect of the extra
 force on closed orbits .    as already seen in a previous work on the mimicking of the dark matter component in galaxies @xcite , considering the tensor character of the field eqs .
( [ field_equations ] ) leads to the following relation between the energy density and pressure of the mimicked dark matter component for a power - law coupling function @xcite :    [ dm_density ] _ dm & = & 1-n1 - 4n2r  , + p_dm & = & n 1 - 4n2r  ,    which , combined with eq .
( [ implicit ] ) , leads to a relation between the visible and dark matter density profiles ,    = 1-n 1 - 4n^1/(1-n )  .
[ dm_visible2 ]    the above clearly shows that one may obtain a `` dark '' component from the visible matter density profile , with a change of slope given by the scaling exponent @xmath78 , _
i.e. _ , @xmath79 .
conversely , if @xmath80 and @xmath81 , one has    = m^mn=1-mm^  , [ n - estimative ]    which leads to an estimate for the exponent @xmath82 of the power - law coupling function , eq .
( [ powerlawcoupling ] ) based on the knowledge of the visible and dark matter density profiles .    for convenience , eq .
( [ second_equation ] ) can be rewritten in a dimensionless form , by defining the rescaled quantities @xcite ,    y & & r_n r  , [ var1 ] + & & ( _ n)^1/(1-n )  , [ var2 ] + & & 2r _ n=1 r r_n  , [ var3 ]    where one defines the length scale , @xmath83 .
the resulting dimensionless form of the trace eq . ( [ second_equation ] ) is    & = & ^-n+(1 - 2n)^n - 6n |(^n-1 )  , [ rhovar ]    where @xmath84 .
the `` static '' equivalent solution , , is thus , @xmath85 , as may easily be checked .
hence , by considering a known visible matter density profile @xmath57 ( which in a cluster is approximately equal to the icm gas density , @xmath86 ) , one may solve the differential equation eq .
( [ rhovar ] ) and thus read the resulting mimicked dark matter profile through eqs .
( [ var3 ] ) and ( [ dm_density ] ) .
this component can account for the large mass difference in the cluster total density profile , as will be discussed in sec .
[ results ]  while the perturbative nature of the non - minimal coupling ensures that no deviation from the geodesic motion occurs .      before further exploring the mimicking mechanism offered by the scaling law @xmath76 , one addresses an issue mentioned in sec .
[ hydro ] , namely , that the ensued modification of the poisson equation due to the presence of a non - minimal coupling would imply that the derivation of is flawed .
one begins by recalling that the poisson equation is obtained from the newtonian limit of the metric ,    ds^2 = -[1 + 2(r)]dt^2 + _
ijdx^idx^j  .[metric ]    inserting the trace of the einstein field equation ( for non - relativistic particles with @xmath87 ) into the @xmath88 component of the latter , yields the well - known result    r = -t 2= 2 r_00 ^2 = 1 4  .    since the above metric yields @xmath89 , one may instead resort directly to the scalar curvature , to obtain the same result :    r = 2r_00 2 ^2 = 1 2  .
since this study focuses particularly on the trace , the latter form is more suitable  but both approaches are completely equivalent , as they depend only on the assumed metric , not the underlying gravity model .    as before , may be written as @xmath90 , with    _ dm = ( 1 - 2n)(r r_n)^n-6n  .
thus , the above shows that one trivially obtains the same form for the poisson equation , .
one can argue that this rewriting obfuscates the effect of the non - minimal coupling , as its _ r.h.s .
_ also involves the scalar curvature @xmath91 .
this , however , misses the main point of this study : one is not aiming at fully determining the newtonian gravitational potential @xmath2 and the interconnected equilibrium density profile @xmath57  something that would require the derivation of a modified hydrostatic equilibrium equation , as followed in ref .
@xcite . instead
, the focus of this study ( cf . the previous work addressing galactic dark matter @xcite ) is to establish a relation between the observed gas density profile and what is usually interpreted as dark matter , not to explore why the former adopts a particular configuration : the mimicking mechanism proposed here acts as a `` translator '' between visible and dark matter .      before dwelling into the analytical and numerical details of the scenario under scrutiny
, one must first consider a possible issue : the presence of an extra - force , as given by eq .
( [ extra - force ] ) . considering a pressureless perfect fluid along with the lagrangian , @xmath92
, the same equation reduces to    f^&=&-_(1+f_2 ) h^  .
for spherical symmetry , the radial component is    f^r&=&-_r ( 1+f_2 ) = - f_2(r)1+f_2(r ) r^  ,    where the prime denotes a derivative with respect to the radial coordinate .
considering a power - law coupling of the form of , one has    f^r & = & -n ( r / r_n)^n1+(r / r_n)^n r^r  .
[ f_r ]    if the curvature @xmath43 decreases along the radius , a negative exponent @xmath82 will produce an inward binding force , and a positive one an outward force , which could break the stability of an orbit .    as already seen in a previous work @xcite , different solutions for the differential equation , eq .
( [ rhovar ] ) , lead to different behaviours of this extra - force , depending on the relevance of the coupling function .
if it dominates , @xmath93 , then the force becomes    f_r - nr^r  .    on the contrary
, a perturbative regime leads to a suppression of the expression above by a factor @xmath94 ,    f_r - n(r r_n)^nr^r  .    the `` static '' solution given by
is derived from the dominant condition @xmath95 ; taking only the outer slope behaviour of cusped density profiles for the dark matter component , @xmath96 , and a positive exponent @xmath82 ( as will be used in the following section ) , one concludes that the extra force has an outward direction ,    f_r - nr^r ~-n_dm _ dm = nm r  .    for @xmath82 and @xmath97 of order unity ,
one obtains a force much greater than the newtonian counterpart , @xmath98 , thus destabilizing the orbital motion of test particles . in the inner region , @xmath99 ,
one has an approximately constant dark matter density , @xmath100 , yielding a negligible extra force .
however , as discussed in the previous section , one finds that the `` dynamical '' solution to must be considered , as the gradient term on the _ r.h.s .
_ of the latter dominates : this solution is essentially proportional to ( disregarding small , short - wavelength oscillations ) .
thus , for a fixed @xmath70 , it is then clear that @xmath101 increases for lower values of a positive exponent @xmath82 ; thus , one expects that there is a lower bound on @xmath82 , such that the condition @xmath102 fails due to the insufficient suppression of the extra force @xmath103 .    in order to assess the value of this lower bound ,
one first uses eqs .
( [ mass_profile ] ) and ( [ dm_profile ] ) to write the following approximations , valid in the large radius region @xmath104 :    [ radial_scaling_0 ] m(r ) ~m_0 ( r r_c)^1 + -   ,   _ dm~_dm0(r r_c)^-2 + -  .
using the above , one may write    f_r & & -n(r r_n)^nr^r ~-n(_dm _ n)^n_dm _ dm ~ + & & ( 2-+)nr_c ( _ dm0 _ n ) ^n ( r r_c)^(-2 + -)n-1  ,    so that , writing @xmath105 , the perturbative condition @xmath106 reads    ( ) ^n ( r r_c)^(-)n-12-n 1  .    one now assumes that @xmath107 ( as is the case with the abell cluster a586 , discussed in the following section ) .
thus , the relevance of the extra force decreases with the distance from the cluster s center , and one may saturate the _ l.h.s .
_ of the above equation by setting @xmath108 , yielding    b(r_n , n ) ( ) ^n r_0 r_c 1  .
[ lower_bound ]    finally , notice that , for a positive exponent @xmath82 , the perturbative condition @xmath109 translates into    ( r r_n)^n 1 _
dm0 2r_n r_0 r_n  .
thus , a lower exponent @xmath82 will increase the _
_ , as expected . in the following section , this lower bound for @xmath82
will be obtained for the abell cluster a586 .
in this section , one tries to describe the mass difference of the abell cluster a586 .
this particular cluster is a nearby massive strong - lensing cluster and has a high x - ray luminosity , directly proportional to the mass of the intracluster x - ray emitting gas .
these features imply a stringent bound for the two measuring methods for its density profile : strong lensing and hydrostatic equilibrium estimatives @xcite .
the assumption of spherical symmetry is an accurate hypothesis for its morphology @xcite .
the abell cluster a586 is also a relaxed cluster , as pointed out in ref .
.      recalling eqs .
( [ visible_profile ] ) and ( [ radial_scaling_0 ] ) , one writes the approximations valid in the large radius regime @xmath104 ,    _
g~(rr_c)^-3 , _ dm~(rr_c)^-2 - 3(-1 )  .
[ radial_scaling ]    therefore , from the scaling given in eq .
( [ dm_visible2 ] ) , one obtains the relation @xmath110 , allowing for an estimate , through eq .
( [ n - estimative ] ) , for the exponent @xmath82 :    [ expected ] n=2 + 3(-2)2 + 3(-1)0.279  ,    using @xmath111 and @xmath112 @xcite .
thus , for the isothermal case ( @xmath31 )    n=1 - 32(2-)0.223  .    hence , the results arising from the non - minimal coupling are somewhat similar whether one considers an isothermal ( @xmath113 ) or a polytropic temperature profile .      in order to mimic the dark matter component of this cluster ,
one numerically solves the differential eq .
( [ rhovar ] ) , for an assumed gas density profile , as given by eq .
( [ visible_profile ] ) , and varies the parameters @xmath70 and @xmath82 in order to obtain a best fit between the derived mimicked dark matter density @xmath114 and the observed curve .
the considered range for the numerical integration lies between @xmath115 , the scharwzschild radius of the cluster , and a sufficiently large radius , larger than the observed size of the cluster .    as a measure of
the quality one tries to minimize the quantity @xmath116 , given by    _ ^2=_i=1^n ( |_i-_i)^2 ( _ i)^2  , [ sigma2 ]    where the sum is performed over a large number of equal spaced radii @xmath117 , the density @xmath118 corresponds to the dark matter profile obtained from the hydrostatic equilibrium estimatives for the radius @xmath117 , and @xmath119 is the value for the obtained mimicked profile in the same radius .
one also presents the deviation of the final total mass of the dark matter component from the quantity    = |m_mim - m|m  , [ mass_deviation ]    where @xmath120 is the amount of dark matter .
the results of the numerical integration are shown in fig .
( [ densities ] ) and the resulting integrated mass radial profile , @xmath121 , in fig .
( [ masses ] ) , for the values of @xmath122 , @xmath123 , @xmath124 , @xmath125 , and @xmath126 , with a best fit attained for @xmath127 ; @xmath128 is kept constant , as its value does not change the mimicked dark matter density profile ( as discussed below ) .
the best fit is given by the minimization of the relative difference between the obtained mimicked total mass and the expected one from the hydrostatic estimative , eqs .
( [ dynamical_mass ] ) and ( [ mass_profile ] ) , with values shown in table ( [ mass_difference ] ) . the relative mass difference is evaluated at the virial radius of the abell cluster a586 , @xmath129 @xcite .
as one can see , these results show that , in the outer regions of the cluster , the density decreases with the value of @xmath82 along with a decrease in the total mass for the cluster .
the best fit @xmath127 is the one that yields a mimicked profile , which minimizes the mass difference to only @xmath130 at @xmath131 , the virial radius for the density profile derived from eqs .
( [ visible_profile ] ) and ( [ mass_profile ] ) .
in fact , in fig .
( [ masses ] ) , one can see that the values in the range , @xmath132 , lead to a dynamical mass larger than the upper limit set by the model , @xmath133 .
notice that this best fit @xmath134 deviates from the expected value @xmath135 , computed in eq .
( [ expected ] ) .
as the curves for @xmath136 and @xmath137 show , this stems from the fact that the mimicked dark matter density arising from the expected value is unable to follow the observed profile .
indeed , while @xmath135 yields the same slope as the observed total mass density , it leads to a right - shifted curve , giving rise to a much greater total mass for the cluster .
, where the dashing length increases with @xmath82 .
the best fit is obtained for @xmath127 ( dark dashed ) .
the horizontal dotted line corresponds to @xmath138 . ]
( polytropic case ) @xcite . ]    .virial radius and respective mass for @xmath39 .
mass differences for different values of the @xmath82 evaluated at @xmath139 . [ cols="^,^,^,^,^",options="header " , ]      , after the integration of the differential equation , for different values of @xmath82 .
]    in fig .
( [ extra_forces ] ) one presents the plot of the additional radial velocity @xmath140 arising from non - geodesic motion as a function of distance to the cluster s center , for the considered range of values for @xmath82 and a fixed @xmath141 . notice
that , despite the fact that the observable quantity is the velocity dispersion @xmath142 , estimated from the velocities of individual clusters , these curves enable one to see the relevance of the extra force that appears due to the presence of the non - minimal coupling .
since the scalar curvature ( _ i.e. _ the total density ) is decreasing , one concludes that the related force counteracts the gravitational attraction , as discussed in sec .
[ subsectionnon ] .
thus , stability requires that this additional velocity be much smaller than the newtonian radial velocity profile : for a fixed value of @xmath143 , this is attained for higher values of @xmath82 . in particular
, one concludes that a power - law coupling with @xmath127 can effectively mimic the dark matter component of the abell cluster a586 , with a negligible effect on the velocity profile @xmath144 .    as already discussed in sec .
[ nongeodesicalmotion ] , one expects a decrease of the effect of the extra force for increasing exponent @xmath82 ( for fixed @xmath143 ) .
indeed , as fig .
( [ couplings ] ) depicts , the considered range for @xmath82 always leads to a perturbative @xmath42 , for a fixed @xmath145 .
however , values of the exponent in the range @xmath146 , which imply @xmath147 , do not enable a sufficiently strong suppression of the extra force : although the additional ( outward ) radial velocity is indeed smaller than @xmath148 , it may still be above the newtonian counterpart  the converse would require one to compute the condition for @xmath137 , in order to obtain a bound on @xmath149 .
+    +      using and the relevant values obtained from ref .
@xcite , one has @xmath150 and @xmath151 . inserting this into , together with @xmath152 and @xmath153 ,
one finds that the extra force is smaller than the newtonian force when @xmath154 , which is equivalent to the constraint , @xmath155 , and should be compared with the aforementioned lower bound @xmath156 , obtained from fig .
[ extra_forces ] .    the lack of sensitivity of the mimicked dark matter density profile on the value of @xmath157 may be used to establish a more stringent upper bound on this parameter . to do so
, one fixes the best fit value @xmath134 and searches for the value of @xmath143 that fulfills the condition @xmath158 , which is equivalent to the constraint @xmath159 .
one now extends the dark matter mimicking mechanism to a large sample of galaxy clusters , aiming to obtain an encompassing description of their `` dark matter '' content . for this purpose ,
one uses a sample of nearby relaxed clusters derived from the high quality data of _ chandra _ @xcite .    as already seen , in the case of abell cluster a586 , one resorts to the scaling @xmath160 to estimate the exponent @xmath82 of the power - law non - minimal coupling .
however , since the density profiles are now more evolved , alternatively one considers the behaviour of the inner and outer regions , that is one resorts to two different exponents , @xmath161 and @xmath162 , respectively .
indeed , the dark matter density profile follows a nfw profile given by _
t~(r r_s ) ^-1(1+r r_s ) ^-2  ,    such that for small radii , @xmath163 , and for outer regions , @xmath164 ; the gas density profile follows a generalized nfw model @xcite ,    ^2_g~n_0 ^ 2(r / r_c)^-(1+r^2/r_c^2)^3-/2  , [ generalized_beta ]    for which , @xmath165 , in inner regions , and @xmath166 for outer regions .
thus , for small radii one has , from eq .
( [ n - estimative ] ) ,    = 1/2n_in=1 - 2  , [ n1 ]    and , for large radii ,    = 33n_out=1-  .
[ n2 ]    the results of eqs .
( [ n1 ] ) and ( [ n2 ] ) for each cluster are shown in table [ clusters_values ] using the reported values for the @xmath167 and @xmath25 parameters @xcite . since @xmath82 is a parameter of the coupling function eq .
( [ powerlawcoupling ] ) , it should be universal : thus , one performs a simultaneous fit of the `` dark matter '' component of all clusters , with a fixed @xmath82 .
the characteristic length scale @xmath143 is varied individually for each cluster , allowing for the analysis of the perturbative regime of the non - minimal coupling eq .
( [ powerlawcoupling ] ) and the extra force eq .
( [ extra - force ] ) , discussed below  following the approach detailed for the abell a586 cluster in the previous section .
numerically , one finds that the best fit occurs for @xmath168 ; the initial conditions ( _ i.e. _ the initial value for @xmath169 and its derivative ) are varied within one order of magnitude , reflecting the uncertainty on the dark matter distributions reported in ref .
the resulting dark matter profiles yield the dashed curves presented in figs .
( [ mimic1 ] ) and ( [ mimic2 ] ) : these minimize the difference between the observed and mimicked dark matter components , leading to small values of the estimator @xmath170 , eq .
( [ sigma2 ] ) , and mass difference @xmath171 , eq .
( [ mass_deviation ] ) , as can be seen in table [ clusters_values ] .
as in the case of abell cluster a586 , one finds that the mimicked dark matter profile is almost independent of the characteristic length scale @xmath172 , as long as the non - minimal coupling eq .
( [ powerlawcoupling ] ) is perturbative , @xmath173 .
this leads to the bound @xmath174 , which is approximately independent of the cluster .
as before , a more stringent bound , @xmath175 , is obtained from the requirement that the extra force eq .
( [ extra - force ] ) is smaller than the newtonian force .
the independence of the obtained mimicked dark matter distributions on the value of @xmath143 again emerges because a different value for the latter leads to a change in the boundary conditions of eq .
( [ rhovar ] ) , which implies a change of the dimensionless solution @xmath169 ; this is compensated by the presence of @xmath143 in the relation between the latter and the mimicked dark matter density , eq .
( [ var1 ] ) .
+    +    +        notice that the fitting of the abell a586 cluster dark matter profile using the universal exponent @xmath176 ( final graph of fig .
[ mimic2 ] ) leads to the values @xmath177 , and @xmath178  naturally higher than those obtained from the @xmath127 best fit of this cluster alone .
there is a fundamental qualitative difference between the density profiles derived for the a586 cluster and the set of 13 chandra clusters : the later do not exhibit a transition between an inner region with approximately constant density and an outer region with decreasing density of approximately constant slope ( as embodied in the @xmath179-model ) .
this drives the ensuing distinction between the fits for the a586 cluster obtained in this study : indeed , as fig .
[ densities ] shows , the best fit @xmath127 obtained by considering only this cluster yielded an initial plateau for the mimicked dark matter density .
this is to be expected , as most of the mass comes from the integration of the density in this higher , approximately constant region  thus leading to a reduced mass difference @xmath180 .    on the contrary , the mimicked dark matter profile of this cluster obtained in the global fit ( final graph of fig .
[ mimic2 ] ) does not have an approximately flat density for the inner region , but compensates by more closely following the observed slope of the outer regime  hence the proximity between the obtained best fit parameter @xmath176 and the value obtained from the scaling law @xmath76 , @xmath181 .
this occurs because the minimization procedure of the global fit is mostly driven by the remaining thirteen chandra clusters : as it turns out , the best fit for this set ( excluding the a586 cluster ) is @xmath182 ( not shown in this study ) ; a balanced fitting solution with all fourteen clusters is thus obtained by relaxing the constraint of a flattened inner density region for the mimicked dark matter profile of the a586 cluster .
the non - minimal coupling considered here may be regarded as an approximation to a more evolved form for @xmath42 , which would encompass couplings considered in previous studies , which perhaps hints at a laurent series , f_2(r ) = _ n = - ^ ( r r_n)^n  .
[ laurent ]    each term would be valid in a particular regime : early _ vs. _ late time , central _
vs. _ long range , galactic _ vs. _ cluster size , _
etc_. in particular , this series should also include the aforementioned @xmath183 and @xmath184 terms , relevant to mimic dark matter at the galactic level , as well as a linear coupling @xmath185 considered in the context of preheating and inflationary dynamics ; the corresponding parameters are    r_-1 & & ( r_-1)^-1/2 = 21.5  gpc  , + r_-1/3 & & ( r_-1/3)^-1/2 = 1.69 10 ^ 6  gpc  , + r_1 & & r_1 ^ -1/2 = 4.8 10 ^ -29  m  ,    with @xmath186 @xcite .
since the relevant coupling terms at galactic scales have negative exponents , it is trivial to conclude that their effect is negligible during inflation , when the curvature is much higher .
conversely , the effect of the terms involved in preheating and inflation are irrelevant at astrophysical or cosmological scales at present .    similarly , the cosmological effect of the best fit coupling for the considered sample of clusters ( with @xmath176 ) may be disregarded , as it is of order    ( r_cosmor_n)^n ~(r_nr_h)^2n ( 10 ^ -13  kpc4.2  gpc)^0.2 ~10 ^ -4  .    where @xmath187 is the hubble radius .
thus , one is left with the possible conflict between astrophysical relevant terms .
if one considers three power - law couplings with @xmath183 , @xmath188 , and @xmath189 , the resulting dark matter profile will not be a sum of the individual contributions , due to the non - linearity of . even if one disregards the resulting interaction between terms ( a reasonable approximation , as shown in ref .
@xcite ) , the exponents simply yield the slope for each derived dark matter contribution , not for their relative values .
indeed , when numerically integrating , one imposes boundary conditions on the overall dark matter density ( _ i.e. _ the dimensionless function @xmath169 )  the relative contributions of each dark matter component are then determined dynamically .
given the scaling relation , a larger @xmath190 implies a steeper dark matter component @xmath191 .
resorting to the @xmath192 behaviour of the visible mass density eq .
( [ generalized_beta ] ) and using the average value @xmath193 , this translates into    n=-1 :  _ dm ~r^-0.98  & , &  n=-1 3 :  _ dm ~r^-1.5  , + n=0.2 & : &  _ dm ~r^-2.2  .
thus , the @xmath168 merely gives rise to a quickly decaying component : the @xmath183 and @xmath184 contributions eventually dominate at large distances , even if the former dominates at short ones .    in order for the contributions from the @xmath194 and @xmath188 to be neglected , it suffices that they are subdominant at cluster scales , but dominant at the long range of galactic scales .
one does not pursue an explicit computation of the effect of these three power - law terms here , given that the main scope was to describe clusters .
nevertheless , it is not difficult to see that the difference between the galactic and cluster typical densities do allow for a consistent disentangling of the effects of these two sets of terms .
in this work , the possibility of mimicking the dark matter component of a galaxy cluster in the context of a model , , where the scalar curvature is non - minimally coupled with matter is discussed .
for that , one first assesses the viability of the model using the relaxed abell a586 cluster , given the availability of quality data for most of its relevant dynamical parameters .    using a power - law non - minimal coupling @xmath195 ,
one obtains the best fit for @xmath127 , with a strong constraint on the typical length scale @xmath196 arising from the requirement that the induced extra - force is smaller than the newtonian one .
a weaker constraint @xmath197 arises from the requirement that the coupling is perturbative , @xmath173 .
the obtained fit does not correspond to the analytical value for @xmath82 computed from the scaling relation between dark and visible matter , as the latter would give rise to a mimicked dark matter density that does not follow the observational curve ( although its outer behaviour exhibits the same slope ) .
instead , a shift from this analytical value is enforced , whereas the mimicked density is slightly steeper than the observed profile , but with a negligible difference in total mass .
furthermore , the exponent @xmath82 is large enough to keep the non - minimal coupling perturbative , with @xmath198 .
this suppresses the destabilizing effect of the outward extra - force arising from the non - conservation of the energy - momentum tensor , as previously found in the context of galaxy rotation curves ( where two negative exponents @xmath183 and @xmath184 were shown to correspond to the best fit for the rotation curves ) @xcite .
subsequently , this dark matter mimicking mechanism was successfully extended to a large sample of galaxy clusters using the available profiles for the visible matter density and temperature obtained from modifications of the usual @xmath25-model , eq .
( [ visible_profile ] ) , and the polytropic temperature profile , eq .
( [ temperature_profile ] ) @xcite . for this purpose ,
one has used the same power - law coupling function , obtaining an exponent @xmath176 from a simultaneous fit to all considered clusters . as in the detailed case of the abell a586 cluster ,
the fitting procedure is also independent of the length scale @xmath143 , as long as the coupling is perturbative , @xmath173 . by the same token
, a more stringent upper bound of @xmath199 is obtained from the requirement that the extra force due to non - geodesic motion is much smaller than the newtonian force .
the fitting procedure depends on the underlying model , but also on several other assumptions  namely , sphericity , virialization , and thermal equilibrium .
clearly , deviations from any of these hypotheses will introduce an undetermined bias into the best fit scenario : the choice of the abell 586 cluster as a initial test case was motivated by the wish to minimize these sources of error , as this cluster is very well characterized and has a good spherical morphology @xcite .    in galaxy clusters with more evolved structures
, one may observe deviations due to the presence of non - sphericity and non - thermal pressure support , and the density profile may be much rougher than in well - behaved clusters .
however , an adequate description of the mimicked dark matter content was obtained for clusters with more evolved internal structures , namely rxj1159 + 5531 , mkw 4 and usgc s152 , thus supporting the universality and robustness of the proposed mimicking mechanism .
studies addressing different scenarios have resorted to power - law couplings with different exponents than the best fit value @xmath168 here reported  which may be regarded as another term in a putative series expansion of a more general non - minimal coupling .
the specific form of the coupling function can only be obtained from its consideration in different scenarios , in order to reconstitute its full functional dependence on the ricci scalar curvature .
thus , the relevance of each term in the several scenarios considered in previous works was also discussed . by the same reasoning , yet unprobed terms of this series may be charted through the study of other phenomena and environments ( ranging from astrophysical to cosmological scales ) , where distinct curvatures and densities are at play .      d.  n.  spergel _ et al . _ [ wmap collaboration ] , , 377 ( 2007 ) [ http://arxiv.org/abs/astro-ph/0603449[arxiv:astro-ph/0603449 ] ] ; g.  hinshaw _ et al .
_ , , 225 ( 2009 ) [ http://arxiv.org/abs/0803.0732[arxiv:0803.0732 ] ] ; e.  komatsu _ et al .
_ , , 330 ( 2009 ) [ http://arxiv.org/abs/0803.0547[arxiv:0708.2750 ] ] ; , 18 ( 2011 ) [ http://arxiv.org/abs/1001.4538[arxiv:1001.4538 ] ]
. m.  tegmark _ et al . _
[ sdss collaboration ] , , 123507 ( 2006 ) [ http://arxiv.org/abs/astro-ph/0608632[arxiv:astro-ph/0608632 ] ] ; m.  sako _ et al .
_ , , 348 ( 2008 ) [ http://arxiv.org/abs/0708.2750[arxiv:0708.2750 ] ] .
m. cirelli , to appear in the proceedings of the 25th international symposium conference [ http://arxiv.org/abs/1202.1454[arxiv:1202.1454 ] ] ; g. bertone , particle dark matter ( cambridge university press , cambridge , uk , 2010 ) .
o. bertolami , f. gil pedro and m. le delliou , 165 ( 2007 ) [ http://arxiv.org/abs/astro-ph/0703462[arxiv:astro-ph/0703462 ] ] ; , 2839 ( 2009 ) [ http://arxiv.org/abs/0705.3118[arxiv:0705.3118 ] ] ; ( 2011 ) [ http://arxiv.org/abs/arxiv:1105.3033[arxiv:1105.3033 ] ] .


low - energy nonperturbative effects in quantum chromodynamics ( qcd ) , especially at nonzero densities , can only be studied by approximate methods within the framework of various effective models .
it is well known that light meson physics is described by four - fermion models , such as the nambu  jona - lasinio ( njl ) model , which was successfully used to deal with dynamical chiral symmetry breaking ( d@xmath0sb ) both in the vacuum and in hot / dense baryonic matter ( see , e.g. @xcite , as well as the reviews @xcite and references therein ) .
recently , much attention has been paid to the effects of diquark condensation and color superconductivity ( csc ) .
the first studies of the gap equations and the ginzburg - landau free energy for a system of relativistic fermions led to the conclusion that superconductive and color superconductive states may arise in baryonic matter @xcite ( see also the recent papers @xcite ) .
another interesting phenomenon , the condensation of charged pions , which may appear in dense hadronic matter due to an asymmetry of its isospin composition , has been investigated in the framework of qcd - like effective models , including the njl model , as well @xcite .
note that all these phenomena might be inherent to physics of compact stars , where rather strong magnetic as well as gravitational fields are present .
therefore , investigations of the influence of an external ( chromo-)magnetic field on the properties of the d@xmath0sb phase transition @xcite , color superconductivity @xcite and pion condensation @xcite effects are quite appropriate . in particular ,
it was shown in @xcite that external fields significantly change the properties of the chiral and csc phase transitions . in several papers , in the framework of the njl model , the influence of a gravitational field on the d@xmath0sb due to the creation of a finite quark condensate
@xmath1 has been investigated at zero values of temperature and chemical potential  @xcite .
the study of the combined influence of curvature and temperature has been performed in  @xcite .
recently , the dynamical chiral symmetry breaking and its restoration for a uniformly accelerated observer due to the thermalization effect of acceleration was studied in  @xcite at zero chemical potential .
further investigations of the influence of the unruh temperature on the phase transitions in dense quark matter with a finite chemical potential , and especially on the restoration of the broken color symmetry in csc were made in @xcite .
one of the widely used methods of accounting for gravitation is based on the expansion of the fermion propagator in powers of small curvature  @xmath2 @xcite .
for instance , in @xcite , the three - dimensional gross - neveu model in a spacetime with a weakly curved two - dimensional surface was investigated , using an effective potential at finite curvature and nonzero chemical potential . in paper  @xcite , this weak curvature expansion was used in considering the d@xmath0sb at non - vanishing temperature and chemical potential .
it should , however , be mentioned that near the phase transition point , one can not consider the critical curvature @xmath3 to be small and therefore the weak curvature expansion method can not be applied . hence , in the region near the critical regime different nonperturbative methods or exact solutions with finite values of @xmath2 should be used .
this kind of solution with consideration for the chemical potential and temperature in the gravitational background of a static einstein universe has been considered in  @xcite .
there it was demonstrated that chiral symmetry is restored at large values of the space curvature .
analogous studies of diquark condensation and the related color symmetry breaking under the influence of a gravitational field have been performed recently in the model of a static einstein universe in @xcite . recall that this model is widely discussed in literature either as a solution of the einstein equations with a given cosmological constant and a nonvanishing energy - momentum tensor of an ideal fluid as a source , or as an initial state in inflationary cosmology with a scalar field , and the cosmological constant as its vacuum energy ( see , for instance , @xcite ) .
moreover , the einstein universe and other suitably generalized compact curved spacetimes were extensively employed for studying the phenomenon of bose - einstein condensation ( see , e.g. @xcite and references therein ) . as is well known , one of the possible models for the expanding universe is the closed friedmann model . since the formation of quark condensates
is expected to take place considerably faster than the expansion rate of the universe , its radius can be considered in our calculations as constant . in this sense , the chosen model of a static einstein universe
can be considered as a simple cosmological model for a space with positive curvature . on the other hand , the same form of the metric
does also hold for the interior of a ( collapsing ) spherically symmetric star , thereby admitting also a non - cosmological interpretation of the chosen model .    as was mentioned above , in quark matter there
may arise another remarkable phenomenon , pion condensation .
up to now little information about the properties of pion condensation in the presence of external fields has been obtained . in particular , there arises the interesting question , what influence gravitational fields can produce on pion condensation .
the main aim of the present paper is just to give a detailed investigation of this issue .    in order to be able to consider the effects of gravitation on pion condensation in a rigorous nonperturbative way
, we will investigate here the phase structure of isotopically asymmetric quark matter in the framework of an extended njl model with dynamical breaking of chiral and flavor symmetries in the static einstein universe .
in particular , our calculations are performed in the most simple case , restricting ourselves , for simplicity , to the flavor @xmath4 group . by using the mean field approximation
, we will then derive an analytical expression for the thermodynamic potential of the njl model that will enable us to calculate the chiral and pion condensates of quarks moving in a space with constant positive curvature . on this basis
, we study , in particular , the phase portraits of the model . our paper is organized as follows .
section ii , iii contain the nambu - jona - lasinio model in curved spacetime and the general formalism for the derivation of the thermodynamical potential .
numerical calculations and conclusions are given in sections iv , v. finally , some estimates of the role of quantum fluctuations in the case of a closed universe are given in the appendix .
suppose that dense , isotopically asymmetric quark matter ( in this case the densities of @xmath5 and @xmath6 quarks are different ) in curved spacetime is described by an extended njl model with the following action : @xmath7 and the lagrangian @xmath8q+ g\big [ ( \bar qq)^2 + ( \bar qi\gamma^5\vec\tau q)^2\big ] ,   \label{2}\end{aligned}\ ] ] where the quark field @xmath9 is a flavor doublet ( @xmath10 or @xmath11 ) and color triplet ( @xmath12 ) as well as a four - component dirac spinor ( the summation in ( [ 2 ] ) over flavor , color and spinor indices is implied ) ; @xmath13 ( @xmath14 ) are pauli matrices . in 4-dimensional
curved spacetime with signature @xmath15 , the line element is written as @xmath16 the gamma - matrices @xmath17 , the metric @xmath18 and the vierbein @xmath19 , as well as the definitions of the spinor covariant derivative @xmath20 and spin connection @xmath21 are given by the following relations  @xcite : @xmath22 ,   \nonumber \\ & & \omega^{\hat{a}\hat{b}}_{\mu}=\frac{1}{2}e^{\hat{a}\lambda } e^{\hat{b}\rho}[c_{\lambda\rho\mu}-c_{\rho\lambda\mu}-c_{\mu\lambda \rho } ] ,
\quad c_{\lambda\rho\mu}= e^{\hat{a}}_{\lambda}\partial_{[\rho}e_{\mu]\hat{a}}. \label{3_3}\end{aligned}\ ] ] here , the index @xmath23 refers to the flat tangent space defined by the vierbein at spacetime point @xmath24 , and the @xmath25 are the usual dirac gamma - matrices of minkowski spacetime .
moreover @xmath26 , is defined , as usual ( see , e.g. , @xcite ) , i.e. to be the same as in flat spacetime and thus independent of spacetime variables .    in order to describe the quark composition of matter we introduced @xmath27 in ( [ 2 ] ) as the quark number chemical potential .
since the generator @xmath28 of the third component of isospin is equal to @xmath29 , the quantity @xmath30 in ( [ 2 ] ) is half the isospin chemical potential , @xmath31 .
if the bare quark mass @xmath32 is equal to zero , then at @xmath33 , apart from the trivial color su(@xmath34 ) symmetry , the lagrangian ( [ 2 ] ) is invariant under the chiral group @xmath35 transformations .
however , at @xmath36 this symmetry is reduced to the subgroup @xmath37 ( here and above the subscripts @xmath38 imply that the corresponding group acts only on the left , right handed spinors , respectively ) .
it is convenient to present this symmetry as @xmath39 , where @xmath40 and @xmath41 are the isospin and axial isospin subgroup of the chiral @xmath35 group .
quarks are tranformed under these subgroups in the following way @xmath42 where @xmath43 , @xmath44 are independent group parameters . at nonzero
bare quark mass , @xmath45 , and nonzero isotopic chemical potential , i.e. @xmath46 , the lagrangian ( [ 2 ] ) is still invariant under the isospin subgroup @xmath40 , but invariance with respect to @xmath41 does no longer hold in this case .    the linearized version of lagrangian ( [ 2 ] ) , that contains collective bosonic fields @xmath47 and @xmath48 @xmath49 , has the following form @xmath50q   -\frac{1}{4g}\big [ \sigma\sigma+\pi_k\pi_k\big ] .
\label{3}\end{aligned}\ ] ] from the lagrangian ( [ 3 ] ) , one can find the equations of motion for the bosonic fields , @xmath51 it is clear from these relations that under isospin @xmath40 and axial isospin @xmath52 transformations ( [ 2010 ] ) the bosonic fields ( [ 200 ] ) are changed in the following way : @xmath53    in the fermion one - loop ( mean field ) approximation , the effective action for the boson fields is expressed through the path integral over quark fields : @xmath54[dq]\exp\bigl(i\int d^4 x\sqrt{-g } ~\tilde    { \mathcal l}\,\bigr),\ ] ] where @xmath55+\tilde { \cal s}_{\rm { eff } } , \label{4}\ ] ] @xmath56 is a normalization constant .
the quark contribution to the effective action , i.e. the term @xmath57 in ( [ 4 ] ) , is given by : @xmath58[dq]\exp\bigl(i\int d^4x\sqrt{-g } ~\bar q{\cal d}q\,\bigr)=\mbox { det}~{\cal d}. \label{5}\ ] ] in ( [ 5 ] ) , we have used the following notation @xmath59 clearly , @xmath60 is an operator in the coordinate , spinor and flavor spaces . apart from this ,
it is proportional to the unit operator @xmath61 in the color space .
thus , it can be presented in the flavor@xmath62color space in the following matrix form : @xmath63 where operators @xmath64 act in the coordinate and spinor spaces only , and @xmath65 due to the trivial color structure , it follows from ( [ 7 ] ) that @xmath66 now , suppose that @xmath67 and @xmath68 are quantities independent of coordinates .. if the fields @xmath69 do not depend on coordinates , then the effective action is a function of the invariants @xmath70 and @xmath71 only ( this fact is due to the symmetry of the model with respect to the transformations ( [ 201 ] ) ) .
therefore , without loss of generality one can put @xmath67 . ] in what follows , we shall denote the pion condensate @xmath72 by @xmath73 for convenience .
then , using the general formula @xmath74=\det [ \bar aa-\bar ab\bar a^{-1}\bar b ] , \label{10}\end{aligned}\ ] ] we obtain @xmath75 - 3i\ln { \rm det}d \equiv -\omega(\sigma,\delta;\mu,\delta\mu)\int d^4x\sqrt{-g } , \label{11}\ ] ] where we have introduced the thermodynamic potential @xmath76 of the system at zero temperature and @xmath77[i\gamma^\nu \nabla_\nu + ( \mu-\delta\mu)\gamma^0 -\sigma - m ] \big\}. \label{12}\ ] ]
the line element in the static einstein universe is defined by the following relation : @xmath78 where @xmath79 is the radius of the einstein universe ( this quantity is related to the scalar curvature , @xmath80 ) ; @xmath81 , @xmath82 , @xmath83 , @xmath84 .
clearly , @xmath85 in this case anticommutes with all @xmath86 and commutes with all @xmath87 , where @xmath88 and @xmath89 .
moreover , @xmath90 .
so , starting from ( [ 12 ] ) , we have ( @xmath88 ; @xmath91 ) : @xmath92 + \delta\mu^2\big\}. \label{18}\end{aligned}\ ] ] let @xmath93 , @xmath94 .
it is evident that @xmath95 .
hence , we have from ( [ 18 ] ) @xmath96    evidently , @xmath97 is an operator in the hilbert space of functions , depending on the space coordinates @xmath98 . as it is well - known ( see , e.g. @xcite ) ,
each of the eigenvalues @xmath99 of this operator is @xmath100-fold degenerate , @xmath101 thus , one can write @xmath102 where the eigenfunctions @xmath103 ( @xmath104 ) of the operator @xmath97 are also known ( see , e.g. , @xcite ) , and @xmath105 ( see ( [ 17 ] ) ) .
now , let us choose in the hilbert space of functions a basis of the form @xmath106 @xmath107 , where @xmath108 . since each element @xmath109 of this basis is an eigenfunction both of @xmath110 and @xmath97
, one can easily conclude that the operator @xmath111 from ( [ 19 ] ) is diagonal in this basis , i.e. each @xmath112 is an eigenfunction of @xmath113 .
the corresponding eigenvalues @xmath114 of @xmath113 have the following form : @xmath115 it is clear from ( [ 192 ] ) that eigenvalues @xmath116 of the operator @xmath113 do not depend on the quantum number @xmath117 , being @xmath100-fold degenerate .
taking into account in ( [ 11 ] ) the relations ( [ 19 ] ) and det@xmath118tr@xmath119 , as well as the results of appendix [ apb ] , where the quantity tr@xmath120 is calculated ( see ( [ b11 ] ) ) , one finds @xmath121=-3i{\rm tr}\ln\hat o=-3i { \cal t}\sum_{l\eta}\int\frac{dp_0}{2\pi}d_l\ln\big [ \delta^2-(p_0+\mu)^2+(\eta e_l+\delta\mu)^2\big ] . \label{193}\ ] ] here @xmath122 , where @xmath123 stands for an infinite time interval and @xmath124 is the space volume of the einstein universe ( the last relations are due to the fact that @xmath125 depends only on @xmath98 ) .
now , using the definition ( [ 11 ] ) of the thermodynamic potential ( tdp ) and summing in ( [ 193 ] ) over @xmath126 , we have for the zero temperature case @xmath127 + \ln\big [ \delta^2-(p_0+\mu)^2+(e_l+\delta\mu)^2\big ] \big\}. \label{194}\ ] ]    to find the tdp @xmath128 in the case of nonzero temperature @xmath129 , one should use the imaginary time technique , where , after summation over matsubara frequencies ( see , e.g. , @xcite ) , the following expression can be found @xmath130 + \ln\big [ 1+e^{-\beta ( e_l^{(+)}-\mu)}\big ] \big\ } \nonumber\\&-&\frac{3t}{{\cal v } } \sum_{l=0}^\infty d_l\big\{\ln\big [ 1+e^{-\beta ( e_l^{(-)}+\mu)}\big ] + \ln\big [ 1+e^{-\beta ( e_l^{(-)}-\mu)}\big ] \big\ } \label{195}\end{aligned}\ ] ] with @xmath131 and @xmath132 .
it is clear that @xmath128 is an even function with respect to each of the transformations @xmath133 or @xmath134 .
thus , one can deal only with non - negative values of the chemical potentials , @xmath135 , @xmath136 .    from this moment on , we will consider only the case of nonzero isospin chemical potential @xmath137 , whereas the baryon chemical potential is set equal to zero , @xmath138 , since its presence is not of principle importance for us .
so , at @xmath138 two particular cases can be investigated on the basis of the tdp ( [ 195 ] ) .
first , let us choose @xmath139 and @xmath138 , but @xmath46 .
then we obtain from ( [ 195 ] ) the expression : @xmath140 secondly , at @xmath141 , @xmath142 , @xmath46 we obtain : @xmath143 + \ln\big [ 1+e^{-\beta e_l^{(-)}}\big ] \big\}. \label{198}\end{aligned}\ ] ]    next , let us consider the limit of zero curvature or infinitely large radius of the universe .
it is clear that the metric ( [ 17 ] ) never coincides with that of the flat minkowsky spacetime because these two spacetimes have different topologies .
however , in the limit @xmath144 and @xmath145 one can obtain from ( [ 198 ] ) the usual expression for the tdp in flat spacetime ( see , e.g. , @xcite ) by the following substitution : @xmath146 in this case , the tdp looks as follows : @xmath147 + 2 t \ln[1+e^{-\beta e_k^{(-)}}]\right\ } ,   \label{199}\end{aligned}\ ] ] where @xmath148 .
it should be noted that one more particular case , when @xmath141 , @xmath149 , but @xmath33 , can easily be reduced to the investigation of the one - flavored njl model at @xmath141 , @xmath149 in the einstein universe @xcite , and hence we shall not consider it here .
first of all , in order to normalize the tdp , we should subtract a corresponding constant from it , such that @xmath150 .
the thermodynamic potential , normalized in this way , is still divergent at large @xmath151 , and hence , we should introduce a ( soft ) cutoff in the summation over @xmath151 by means of the multiplier @xmath152 @xcite , where @xmath153 is the cutoff parameter . can be specified according to the experimental results . however , in the case of a curved spacetime , in order to fix the cutoff @xmath154 , we need new theoretical / experimental inputs for chiral qcd in ( strong ) gravitational background fields , concerning , for instance , an effective gluon mass , known values of the quark condensate or even experimentally measured characteristics of a pion . due to the lack of experimental knowledge ,
our aim is here to perform only a qualitative study of gravitational effects on the quark and pion condensates by investigating the respective phase portraits of the njl model .
for this reason it is convenient to scale the thermodynamic potential and all relevant quantities like condensates , curvature , chemical potential and temperature by the cutoff @xmath153 . ]    for convenience , we shall multiply all dimensional quantities that enter the thermodynamic potential by the corresponding power of @xmath153 to make them dimensionless , i.e. , @xmath155 , and denote them using the same letters as before : @xmath156 . then the regularized thermodynamic potential is written as @xmath157 + \ln\big [ 1+e^{-\beta e_l^{(-)}}\big ] \big\}. \label{198_1}\end{aligned}\ ] ]    in the following section , we shall perform a numerical calculation of the points of the global minimum of the finite regularized thermodynamic potential @xmath158 ( they should of course coincide with the minima of the potential @xmath159 ) , and with the use of them , consider quark matter phase transitions in the gravitational field of the einstein universe .
in this section , on the basis of the thermodynamic potentials ( [ 197])-([198 ] ) , we will study numerically phase transitions in quark matter and consider only the case of nonzero isospin chemical potential @xmath137 , whereas the baryon chemical potential is set equal to zero , @xmath138 . in order to obtain the values of condensates , one should find the global minimum point ( gmp ) of the thermodynamic potentials over the variables @xmath160 , @xmath73 from the corresponding gap equations @xmath161 formally , there are four different expressions for the gmp : i ) @xmath162 , ii ) @xmath163 , iii ) @xmath164 , iv ) @xmath165 .
the first two gmps correspond to the isotopically invariant phases of the model , whereas the gmps of the form iii ) and iv ) correspond to the phases , in which the ground state is no more @xmath40 invariant . in these phases
the pion condensation phenomenon occurs . for simplicity , we take for the numerical calculations of the gmp the value of the coupling constant @xmath166 ( in our dimensionless choice of parameters ) .
let us first consider phase transitions at zero temperature , @xmath139 and choose the current quark mass to be equal to zero , @xmath167 .
the thermodynamic potential in this case is described by formula ( [ 197 ] ) .
the detailed investigation of the gmp properties vs external parameters @xmath2 and @xmath30 results in the phase portrait shown in fig .
[ pp1 ] . for the points of symmetric phase 1 ,
the gmp is at @xmath168 and @xmath169 . in the phase 3 ,
the minimum is at @xmath168 and @xmath170 , and this indicates that the isospin @xmath40 symmetry of the model is dynamically broken in this phase .
we note that at the same time the chiral @xmath52 symmetry remains unbroken in this phase , since in the gmp we have @xmath168 , as it should be in the case of zero current quark mass , when @xmath170 .
@xmath171    as one can see from fig .
[ pp1 ] , the critical curve , which separates the phases 1 and 3 , has an oscillating character .
this phenomenon is explained by the discreteness of the fermion energy levels ( [ 190 ] ) in compact space . moreover , as it is clearly seen from the curves in [ cond23 ] , the pion condensate @xmath73 vs @xmath2 also oscillates in the phase 3 .
this effect resembles the van alphen - de haas oscillations of different physical quantities in the magnetic field @xmath172 , where fermion levels are also discrete ( the landau levels ) @xcite ( see also @xcite , where a similar influence of a magnetic field on the oscillation behavior of the compton scattering and photoproduction cross - sections was demonstrated ) . indeed , the corresponding magnetic oscillations of the critical curve in the @xmath173-@xmath172 phase portrait of dense cold quark matter with four - fermion interactions were found in papers @xcite .
there , the existence of the standard van alphen - de haas magnetic oscillations of some thermodynamical quantities , including magnetization , pressure and particle density of cold dense quark matter was also demonstrated .
the behavior of the pion condensate @xmath73 as a function of @xmath2 at fixed @xmath30 and as a function of @xmath30 at fixed @xmath2 is shown in fig .
[ cond23 ] ( left and right pictures respectively ) .
the phase portrait at finite current quark mass , @xmath45 , is depicted in fig .
[ pp2 ] . in phase 2 the chiral symmetry
is now broken due to a finite value of the current quark mass , and the global minimum of tdp is at @xmath174 and @xmath175 . in the mixed phase 4 both condensates are nonzero , i.e. @xmath174 and @xmath170 .    the behavior of the chiral condensate @xmath160 in the case of finite quark mass , @xmath45 , is shown in fig .
[ cond5 ] as a function of @xmath2 at @xmath176 ( left picture ) and as a function of @xmath30 at @xmath177 ( right picture ) .
@xmath178    one can see oscillations of @xmath160 on both pictures , although they are rather strong in the left picture for the dependence on @xmath2 , while they are weakly seen in the high @xmath179 tail of the curve in the right picture .
one should note that the right picture resembles , except for these oscillations , the corresponding curve in @xcite ( fig .
1 ) for the flat case .      using formula ( [ 198_1 ] ) for the thermodynamic potential , we can also study the influence of finite temperature on phase transitions . the phase portrait at @xmath180 and zero current quark mass , @xmath181 , is shown in fig .  [ pp3 ] in terms of @xmath182 .
it is seen from this figure that growing temperature leads to a smoothing of oscillations of the phase curve .    for comparision , in fig .
[ pp4 ] , the phase portraits in the @xmath183 and @xmath182 planes ( the latter now at another value of temperature @xmath184 ) are depicted .
first of all , it is clear from fig .
[ pp4 ] that the isospin symmetry is restored due to the vanishing of the pion condensate both at high temperature and high curvature .
the similarity between these two plots leads to the conclusion that curvature and temperature effects play a similar role in the restoration of isospin symmetry .
@xmath185    the phase portrait in fig .
[ pp5 ] for @xmath186 corresponds to the case of flat minkowski spacetime .
it looks very similar to the one obtained , for instance , in @xcite ( see upper panel of their corresponding fig .
let us compare the phase portraits at zero ( fig .
[ pp5 ] ) and finite curvature ( fig .
[ pp4 ] , left picture ) . in the first case the pion condensation appears at arbitrary small nonzero values of the isospin chemical potential , while in the second case the isospin symmetry becomes dynamically broken only at some finite value of the chemical potential @xmath30 .
this phenomenon may be explained by the existence of a gap in the quark spectrum ( [ 190 ] ) , which is proportional to the inverse radius of the einstein universe . in this case
the effect of curvature is similar to the effect of a nonzero current quark mass ( for comparison see lower panel of fig . 12 in @xcite and our fig .
[ pp4 ] ) .
in the framework of an extended nambu  jona - lasinio model , we have studied the influence of a gravitational field on pion condensation in isotopically asymmetric quark matter at finite temperature and isospin chemical potential . as a particular model of a gravitational field configuration
we have taken the static einstein universe .
this particular choice enabled us to investigate phase transitions of the system with an exact consideration of the role of the gravitational field in the formation of the quark and pion condensates and thus to demonstrate its influence on the phase portraits .
in particular , we have found that thermodynamic quantities such as quark and pion condensates as well as the corresponding phase boundaries ( critical curves ) oscillate as functions of curvature .
this oscillating behavior is smoothed out with growing temperature .
there exists also an interesting similarity between the behavior of phase portraits when considered as functions of curvature and as functions of finite temperature ( compare figs .
moreover , we have shown that for massless quarks and for some values of @xmath2 the isospin symmetry becomes in curved spacetime dynamically broken only at some finite value of the chemical potential @xmath30 ( look at fig . 6 ( left picture ) for rather small values of @xmath129 ) .
this is contrary to the case of flat spacetime , where the pion condensation appears at arbitrary small nonzero values of the isospin chemical potential ( see fig .
this effect resembles the pion condensation but at nonzero bare quark mass and may be explained by the presence of a gap in the energy spectrum of quarks in the gravitational field .
finally , let us add some remarks on the possible role of quantum fluctuations of the collective fields and finite size effects in our approach .
in fact , since the volume of the space region modelled by the closed einstein universe adopted in this paper is limited , finite size effects could eventually change the character of phase transitions ( see e.g. discussions in @xcite ) .
eventually , this might even lead to particular situations , where no phase transition can occur .
however , as it is clear from physical considerations , the finite size in itself may not practically forbid the dynamical symmetry breaking , if the characteristic length of the region of space occupied by the system is much greater than the compton wavelength of the excitation responsible for tunneling and restoration of symmetry ( see , e.g. , @xcite ) .
it should be noted that our results , obtained with the use of the mean field approximation in the framework of the njl model are consistent with corresponding estimates of the role of quantum fluctuations ( see appendix b ) .    in conclusion
, we emphasize that the results of this paper are evidently only of a qualitative nature .
they do not allow us to find an exact value of the critical radius , and hence , further studies with more realistic models of gravitational fields should be undertaken .
two of us ( v.zh . and
thank m. mueller - preussker for hospitality during their stay at the institute of physics of humboldt - university , where part of this work has been done , and also daad for financial support .
d.e . is grateful to r. seiler for discussion on the role of quantum fluctuations and finite size effects .
this work was supported in part by the dfg - grant 436 rus 113/477 .
in section [ eu ] we have introduced two operators , @xmath110 and @xmath97 , acting in the hilbert space @xmath187 of all quadratically integrable functions @xmath188 defined on the spacetime manifold of the einstein universe .
now suppose that there is an abstract hilbert space @xmath189 of vectors @xmath190 .
let @xmath191 and @xmath192 be the operators of the space and time coordinates , correspondingly , defined on @xmath193 .
moreover , let @xmath194 be the complete set , or basis , of eigenvectors of @xmath191 and @xmath195 , i.e. @xmath196 , @xmath197 .
the set @xmath198 is usualy called the coordinate basis in @xmath193 .
obviously , the completeness and normalization conditions for the coordinate basis @xmath198 are valid : @xmath199 where @xmath200 is the unit operator in @xmath193 and @xmath201 .
due to ( [ b1 ] ) , it is possible to expand any vector @xmath202 in terms of the basis @xmath198 , namely @xmath203 .
the quantity @xmath204 is called _
@xmath24- ( or coordinate ) representation for the vector @xmath190_. identifying @xmath205 with @xmath188 , we see that the hilbert space @xmath187 of all quadratically integrable functions is simply the coordinate representation of the above introduced abstract hilbert space @xmath193 .
furthermore , if in the hilbert space @xmath193 there is an arbitrary operator @xmath206 , then the matrix @xmath207 , whose matrix elements in the coordinate basis are just the quantities @xmath208 , is called the @xmath24- ( or coordinate ) representation of the operator @xmath206 .
obviously , in order to define an operator @xmath206 in the abstract hilbert space @xmath193 , it is sufficient to define its matrix @xmath209 , which acts in the hilbert space @xmath187 of all quadratically integrable functions , in the coordinate basis .
it is clear that tr@xmath210@xmath211 .
now , let us consider in @xmath193 the two commuting operators @xmath212 and @xmath213 such that in the @xmath24-representation they look like @xmath214 and latexmath:[$\hat{\cal h}\equiv \langle x '    operators @xmath110 and @xmath97 have a common set of eigenfunctions @xmath112 defined in section [ eu ] , from which it follows that @xmath216 the eigenfunctions @xmath112 are the coordinates of the corresponding eigenvectors @xmath217 of the operators @xmath212 and @xmath213 ( recall , @xmath218 ) : @xmath219 ( clearly , @xmath220 and @xmath221 . )
it follows from ( [ b4 ] ) that @xmath222 . using this relation in the normalization condition ( [ b3 ] ) and then integrating there over @xmath24 with the help of ( [ b1 ] )
, we obtain @xmath223 which is the analogue of the normalization condition ( [ b2 ] ) .
it is possible to show that the completness condition for the basis @xmath224 follows from ( [ b5 ] ) : @xmath225 let us construct in @xmath193 the following operator @xmath226 which is diagonal in the basis ( [ b4 ] ) , i.e. each vector @xmath224 is its eigenvector with corresponding eigenvalue @xmath227 ( [ 192 ] ) . in the coordinate representation
its matrix @xmath228 coincides with the operator @xmath113 ( [ 19 ] ) .
so , @xmath229 where the last equality was obtained by employing the completeness relation ( [ b6 ] ) .
now , by using in this formula the eigenvalue condition @xmath230 , the normalization condition ( [ b5 ] ) , and , finally , by performing in the obtained expression the integration and summation over primed indices , we have @xmath231 since in ( [ b9 ] ) the quantities @xmath232 and @xmath233@xmath234 ( the last relation is due to the notations accepted in formula ( [ 191 ] ) and below ) do not depend on the time coordinate , the expression ( [ b9 ] ) is proportional to the infinite time interval @xmath235 .
the remaining @xmath98-integration in ( [ b9 ] ) gives simply unity due to the relation ( [ 191 ] ) .
so we have @xmath236 , \label{b10}\end{aligned}\ ] ] where the fact that each eigenvalue @xmath237 is @xmath100-fold degenerated is taken into account and the notations from ( [ 190])-([192 ] ) are used . in a similar way it is possible to obtain the quantity
tr@xmath120 : @xmath238 . \label{b11}\end{aligned}\ ] ]
it is well known that spontaneous symmetry breaking in low dimensional quantum field theories may become impossible due to strong quantum fluctuations of fields .
the same is also true for systems that occupy a limited space volume .
however , as it is clear from physical considerations , the finite size in itself may not practically forbid the spontaneous symmetry breaking , if the characteristic length of the region of space occupied by the system is much greater than the compton wavelength of the excitation responsible for tunneling and restoration of symmetry .
( indeed , one may recall here well known physical phenomena such as the superfluidity of helium or superconductivity of metals that are observed in samples of finite volume ) .
this idea has been discussed for some scalar field theories in the einstein universe for instance , in @xcite . in this appendix , we shall demonstrate that , under certain conditions , dynamical symmetry breaking in njl - type models is indeed possible in the closed einstein universe .
in particular , we will show that , if the radius of the universe is large enough such that the fluctuations of quantum fields are comparatively small , the symmetry breaking obtained in the mean field approximation is not forbidden .    for illustrations ,
let us confine to the analogous case of the chiral condensate and consider the simplified case of the linearized lagrangian ( [ 3 ] ) with @xmath239 , @xmath240 and the corresponding partition function @xmath241 with @xmath242 now , supposing that @xmath243 , where @xmath244 is the vacuum expectation value of the field @xmath160 and @xmath245 denotes its quantum fluctuation , we obtain : @xmath246 where @xmath247 . thus @xmath248 , where @xmath249 , @xmath250 , and @xmath251 is the effective potential at @xmath252 .
the contribution of quantum fluctuations up to the @xmath253-term to the effective action @xmath254 is given by @xmath255 it is evident , that in the above expansion the term linear in @xmath245 corresponds to the so - called tadpole diagram with one external @xmath245-line and the term quadratic in @xmath245 corresponds to the `` polarization operator '' diagram of the @xmath245 field with one fermion loop . from ( [ v01 ] )
we can write the stationarity condition and find the gap equation , @xmath256 , @xmath257 and the linear terms in @xmath245 , corresponding to the tadpole diagram , cancel out in @xmath258 .
thus , the contribution of fluctuations to the field action is given by @xmath259 next , we shall calculate the contribution of fluctuations @xmath245 to the effective action , taking into account the quark loop in the gravitational field .
( note that this corresponds to the integral ( 11 ) and fig.1a in @xcite . ) for our purpose of making estimates of the role of fluctuations , it is sufficient to limit ourselves to the consideration of fluctuations depending only on time . in this case , we can extract the necessary kinematical factor for the meson fluctuation field and then integrate over the time - component of the loop momentum in the limit of vanishing external momenta .
finally , after going to the basis for the dirac equation in the einstein universe ( see ( [ 19])([191 ] ) ) we obtain a sum over fermion loop quantum numbers @xmath151 instead of an integration over momenta of free quarks made in @xcite .
the sum is divergent and we regularize it by the cut off @xmath153 . in this way we obtain the effective action @xmath260 , \label{phiaction}\ ] ] where @xmath261 and @xmath262 are given in ( [ 190 ] ) .
the summation over @xmath151 in the first term in parenthesis in the above equation cancels out by the term @xmath263 , due to the stationarity condition ( [ gap ] ) @xmath264 after this we obtain @xmath265 or @xmath266 with the lagrange function @xmath267 where @xmath268 is the space volume , and the renormalization @xmath269-factor is defined as @xmath270 for comparision , let us consider the limiting case of flat space which is reached by the replacements @xmath271 in ( [ major ] ) .
then the @xmath269-factor takes the form of the integral in the euclidean spacetime @xmath272 with @xmath273 being the quark energy .
this expression evidently corresponds to the similar formula for the @xmath274-factor in the flat space case considered in @xcite .
next , let us perform the field renormalization @xmath275 in ( [ lagrange ] ) .
quantum fluctuations of the boson field near the ground state @xmath244 can now be estimated , if we consider the renormalized expression ( [ lagrange ] ) as the lagrange function for a harmonic oscillator ( here , we follow the idea of @xcitemodel of a self - interacting scalar field was considered and the scalar loop contribution to the fluctuation lagrangian was calculated . ] ) with the mass @xmath32 and frequency @xmath276 , formally given here by the relations @xmath277 then we can estimate the quantum fluctuations as @xmath278 where m is the mass of the composite @xmath160-meson .
thus , we obtain @xmath279 the estimate ( [ estimate ] ) gives a criterion for the role of quantum fluctuations for a system with finite volume . clearly , quantum fluctuations can be considered negligible , if @xmath280 .
this is in agreement with the physical requirement that quantum fluctuations should be negligible if @xmath281 , i.e. , if the radius of the universe is much greater than the compton wavelength @xmath282 of the @xmath160-meson ( quarks ) ( see also @xcite ) .
, which are of order @xmath283 @xcite .
all these corrections are surely suppressed in the case of large numbers of colors @xmath284 , where the ( induced ) coupling constants become small . ]
the above estimates are certainly of a qualitative nature , and hence they do not allow us to find an exact value of a critical radius such that symmetry breaking for lower values of the curvature radius of the einstein universe is forbidden .          m.k .
volkov , fiz .
atom.yadra * 17 * , 433 ( 1986 ) ; m.k .
volkov and a.e .
radzhabov , phys .
* 49 * , 551 ( 2006 ) ; a.a .
andrianov , d.  espriu , and r.  tarrach , nucl .
b * 533 * , 429 ( 1998 ) ; s.v .
molodtsov , a.n .
sissakian , a.s .
sorin , and g.m .
zinovjev , arxiv : hep - ph/0702178 .
klevansky , rev .
phys .   *
64 * , 649 ( 1992 ) ; t.  hatsuda and t.  kunihiro , phys .
rept .   * 247 * , 221 ( 1994 ) .
m. alford , k. rajagopal , and f. wilczek , phys .
b 422 , 247 ( 1998 ) ; nucl . phys .
b * 537 * , 443 ( 1999 ) ; s.v .
molodtsov and g.m .
zinovjev , phys .  atom .
nucl .   * 66 * , 1349 ( 2003 ) .
g.  nardulli , riv .
nuovo cim .  * 25n3 * , 1 ( 2002 ) ; m. buballa , phys . rep . * 407 * , 205 ( 2005 ) ; i.a .
shovkovy , found .
phys . * 35 * , 1309 ( 2005 ) ; t.  ohsaku , phys .
b * 634 * , 285 ( 2006 ) ; m.g .
alford , a.  schmitt , k.  rajagopal , and t.  schafer , arxiv:0709.4635 .
d.  blaschke , d.  ebert , k.g .
klimenko , m.k .
volkov , and v.l .
yudichev , phys .
d * 70 * , 014006 ( 2004 ) ; d.  ebert , k.g .
klimenko , and v.l .
yudichev , phys .
c * 72 * , 015201 ( 2005 ) ; eur .
j.   c * 53 * , 65 ( 2008 ) ; d.  ebert and k.g .
klimenko , theor .
phys . * 150 * , 82 ( 2007 ) .
d.t .  son and m.a .
stephanov , phys .  atom .
* 64 * , 834 ( 2001 ) ; j.b .
kogut , and d. toublan , phys .
d * 64 * , 034007 ( 2001 ) ; j.b .
kogut , and d.k .
sinclair , phys .
d * 66 * , 014508 ( 2002 ) .
andrianov and d.  espriu , arxiv:0709.0049 ; m.  loewe and c.  villavicencio , braz .
j.  phys .
* 37 * , 520 ( 2007 ) ; j.  erdmenger , m.  kaminski , and f.  rust , arxiv:0710.0334 ; s.  mukherjee , m.g .
mustafa , and r.  ray , phys .
d * 75 * , 094015 ( 2007 ) ; h.  abuki , m.  ciminale , r.  gatto , n.d .
ippolito , g.  nardulli , and m.  ruggieri , arxiv:0801.4254 . k.g .
klimenko , arxiv : hep - ph/9809218 ; d.  ebert , k.g .
klimenko , m.a .
vdovichenko , and a.s .
vshivtsev , phys .  rev .
d * 61 * , 025005 ( 2000 ) ; m.a .
vdovichenko , a.s .
vshivtsev , and k.g .
klimenko , phys .
* 63 * , 470 ( 2000 ) .
vshivtsev , v.ch .
zhukovsky , k.g .
klimenko , and b.v .
magnitsky , phys .
* 29 * , 523 ( 1998 ) ; d. ebert and k.g .
klimenko , nucl . phys . * a 728 * , 203 ( 2003 ) ; t.  inagaki , d.  kimura , and t.  murata , prog .
theor .  phys .
suppl .   * 153 * , 321 ( 2004 ) ; a.a .
osipov , b.  hiller , a.h .
blin , and j.  da providencia , phys .
b * 650 * , 262 ( 2007 ) ; arxiv:0802.3193 ; t.d .
cohen , d.a .
mcgady , and e.s .
werbos , phys .
c * 76 * , 055201 ( 2007 ) ; p.  castelo ferreira and j.  dias de deus , arxiv:0707.4200 ; e.s .
werbos , arxiv:0711.2635 .
d.  ebert , k.g .
klimenko , and h.  toki , phys .
d * 64 * , 014038 ( 2001 ) ; v.ch .
zhukovsky , v.v .
khudyakov , k.g .
klimenko , and d. ebert , pisma zh .
fiz . * 74 * , 595 ( 2001 ) ; d.  ebert et al .
d * 65 * , 054024 ( 2002 ) ; t.  tatsumi , e.  nakano , and k.  nawa , arxiv : hep - ph/0506002 ; e.j .
ferrer , v.  de la incera , and c.  manuel , phys .
lett .   *
95 * , 152002 ( 2005 ) ; nucl .
b * 747 * , 88 ( 2006 ) ; j.l .
noronha and i.a .
shovkovy , phys .
d * 76 * , 105030 ( 2007 ) .
d.  ebert , k.g .
klimenko , v.c .
zhukovsky , and a.m.  fedotov , eur .
j.   c * 49 * , 709 ( 2007 ) ; vestn .
astron .   * 61n2 * , 69 ( 2006 ) . t. inagaki , s.d .
odintsov , and t. muta , prog .
. suppl . * 127 * , 93 ( 1997 ) ( see also further references in this review paper ) .
d. ebert and v.ch .
zhukovsky , phys . lett .
b * 645 * , 267 , ( 2007 ) .
( the paper contains a misprint in the definition of @xmath285 matrices in the charge conjugation operator _
c _ , which should read like this @xmath286 .
the results of the paper , however , were obtained with the above correct expression for this operator and do not depend on this misprint . ) t.s .
bunch and l. parker , phys .
d * 20 * , 2499 ( 1979 ) .
kim and k.g .
klimenko , j. phys .
a * 31 * , 5565 ( 1998 ) .
a. goyal and m. dahiya , j. phys .
g * 27 * , 1827 ( 2001 ) .
x. huang , x. hao , and p. zhuang , astropart .
* 28 * , 472 ( 2007 ) .
d. ebert , a.v .
tyukov , and v.ch .
zhukovsky , phys .
d * 76 * , 064029 ( 2007 ) .
klimenko , theor .
phys . * 70 * , 87 ( 1987 ) .
p.  elmfors , d.  persson , and b.s .
skagerstam , phys .
rev .  lett .
* 71 * , 480 ( 1993 ) ; j.o .
andersen and t.  haugset , phys .
d * 51 * ( 1995 ) 3073 ; a.s .
vshivtsev , k.g .
klimenko , and b.v .
magnitsky , j.  exp .
phys .   * 80 * , 162 ( 1995 ) ; j.  exp .
phys .   *
82 * , 514 ( 1996 ) .
zhukovsky and j.herrmann , yad . fiz .
* 14 * , 150 ( 1971 ) ; * 14 * , 1014 ( 1971 ) .