mass gap generation and ( chiral ) symmetry breaking in relativistic fermionic systems can arise from a variety of mechanisms which are often related to certain couplings or interaction channels becoming dominant .
an apparent counter - example is catalyzed symmetry breaking , first studied in the context of magnetic catalysis @xcite , where mass gap generation is triggered by the presence of a magnetic field even for arbitrarily small values of the interaction strength .
this phenomenon can be understood in various ways , the essence being that the long - range fluctuations driving the symmetry - breaking transitions are parametrically enhanced , see @xcite for a recent review .
more concretely , a magnetic field induces a fermionic fluctuation spectrum with landau - levels containing a zero mode that leads to an enhancement of the density of states in the ir and to an effective dimensional reduction favoring symmetry breaking .
magnetic catalysis has found a rich variety of applications both in particle physics ( chiral phases of qcd ) @xcite and condensed matter physics @xcite .
another simple picture for magnetic catalysis has recently been developed within the framework of the functional renormalization group in the context of the @xmath0 gross - neveu model @xcite . in line with the fact that symmetry - breaking phase transitions are often related to fixed points of renormalization group transformations
, also magnetic catalysis can be related to the behavior of rg fixed points as a function of the magnetic field .
this rg picture has already successfully been applied in the context of qcd @xcite .    in the present work ,
we verify the underlying rg picture of catalyzed symmetry breaking in the context of curved spacetimes .
the fact that symmetry breaking and mass generation in fermionic systems can be influenced by negative curvature of the spacetime has been realized early @xcite , and is meanwhile reviewed in textbooks @xcite .
the phenomenon is typically studied at mean - field level and occurs in many different fermionic models @xcite . in @xcite ,
the similarity to magnetic catalysis was realized in terms of an effective dimensional reduction mechanism of the spectral properties of the dirac operator .
this justifies the use of the terminology `` gravitational catalysis '' @xcite .    in fact , we find in the present work that the effective dimensional reduction and the corresponding enhancement of the density of states in the ir is directly related to the fixed point structure as identified below . from this rg viewpoint ,
symmetry breaking arises as a consequence of the fact that the coupling value required for criticality becomes arbitrarily small as a function of the curvature ( the catalyzer ) .
hence , any finite value of the fermionic interactions ultimately becomes supercritical , typically driving the system towards the ordered phase .
we investigate this rg mechanism within the simple @xmath0 gross - neveu model in the present work .
we concentrate on two different curved backgrounds with constant negative curvature : a maximally symmetric spacetime ( anti de sitter ) and a purely spatially curved case ( lobachevsky plane ) . for both cases ,
mean - field studies are already available , see @xcite and @xcite , respectively . whereas the former allows for an analytic treatment in terms of simple functions ,
the latter is potentially relevant for curved layered condensed matter systems .
for instance , the exitonic or anti - ferromagnetic instabilities in graphite and graphene have been associated with quantum phase transitions falling into the @xmath0 gross - neveu universality class @xcite . as catalyzed symmetry breaking
is manifestly driven by the long - range modes , the rg analysis allows us to estimate the required curvature in relation to the length scale of the sample .
this paper is organized as follows : in sect .
[ sec : gnmodel ] , we briefly introduce the model and its formulation in curved space .
section [ sec : frg ] is devoted to an evaluation of the rg flow in its simplest formulation .
the manifestation of gravitational catalysis is discussed in sect .
[ sec : gc ] .
we estimate the influence of finite probe length on gravitational catalysis in sect . [
sec : pseudo ] by means of a pseudo - critical coupling .
conclusions are drawn in sect .
[ sec : conc ] .
relevant technical details are deferred to the appendices .
we aim at investigating the @xmath0 gross - neveu model @xcite in curved spacetime with metric @xmath1 ( greek indices running from @xmath2 to @xmath3 ) and signature @xmath4 using functional rg methods .
the microscopic action functional @xmath5 at some ultraviolet ( uv ) scale @xmath6 depends on the bare coupling constant @xmath7 , the @xmath8 grassmann - valued fields @xmath9 and the @xmath8 conjugated fields @xmath10 , @xmath11 = { } & \!\ ! { \int \!\!\ ! { } _ { { } _ { { } _ { { } _ { { } _ \text{\small{\ensuremath{x } } } } } } } } \ ! \left [   \sum\limits_{i = 1}^{{n } } \bar{\psi}^{i } \slashed{\nabla } \psi^{i } +   \frac{\bar{\lambda}_{\lambda}}{2 { n } } \left ( \sum\limits_{i = 1}^{{n } } \bar{\psi}^{i } \psi^{i } \right)^{\!\!\!2 } \right ]   \notag \\   = { } & \!\ ! { \int \!\!\ ! { } _ { { } _ { { } _ { { } _ { { } _ \text{\small{\ensuremath{x } } } } } } } } \left [ \bar{\psi } \slashed{\nabla } \psi +     \frac{\bar{\lambda}_{\lambda}}{2 { n } } ( \bar{\psi } \psi ) ^2 \right ] \ !   \text { , } \label{eq : gnaction}\end{aligned}\ ] ] where @xmath12 is
a shorthand for the integral over the @xmath13-dimensional spacetime and @xmath14 is the determinant of the spacetime metric .
the differential operator @xmath15 is composed from a set of @xmath16-dimensional gamma matrices satisfying the clifford algebra @xmath17 where @xmath18 is the identity in spinor space . in the present work , we use the irreducible representation , such that @xmath19 specifies the dimension of the @xmath20 matrices as well as the number of dirac components of the fermions . for our computations with fermions in curved space
, we have actually used the weldon formalism @xcite which can be viewed as a generalization of the conventional vierbein formalism . however , a mapping to standard vierbein formulas is straightforward .
the covariant derivative in the weldon formalism reads @xmath21 accounting for covariance with respect to the spinor and spacetime structure . here , @xmath22 is the affine spin connection , implicitly defined by @xmath23,\\   { \ensuremath{(\mathrm{ii})}}\quad { } & 0 = { \operatorname{tr}}\gamma_{\mu } \text { , } \end{aligned}\end{aligned}\ ] ] with @xmath24 the christoffel symbols .
can additionally accommodate a @xmath25 gauge field @xmath26 @xcite . in this case , the right - hand side of @xmath27 in ( [ eq : defspincon ] ) would be given by @xmath28 , where @xmath29 is the charge under the @xmath25 gauge group . in the present work ,
we ignore such a @xmath25 gauge sector , setting @xmath30 . ] in the present work , we are interested in a discrete `` chiral '' @xmath31 symmetry , where the nontrivial transformation is defined by @xcite @xmath32 this symmetry acts simultaneously on all flavors .
it can spontaneously be broken by a chiral condensate @xmath33 , which for finite interactions goes along with a mass gap generation . incidentally
, the @xmath0 gross - neveu model actually has a much larger continuous u(@xmath8 ) flavor symmetry also allowing for more complicated breaking patterns @xcite .
condensed matter systems where the gross - neveu model is considered as an effective theory , the low - energy degrees of freedom can be arranged into @xmath34 4-component dirac spinors , corresponding to a reducible representation of the dirac algebra .
this reducible representation can be constructed from a suitable combination of 2-component spinors such that @xmath35 in terms of the counting of fermions of the present work , see , e.g. , @xcite for a review .
note that the gross - neveu interaction term considered in this work @xmath36 corresponds to @xmath37 in the reducible 4-component notation of @xcite ( or to @xmath38 in the notation of @xcite ) for even @xmath8 .
the critical properties of the discrete chiral transition are , however , identical to a 4-component gross - neveu model with a @xmath39 interaction as considered in @xcite . ]
in the following , we use the functional renormalization group to compute the rg flow of the gross - neveu coupling as a function of the ( negative ) curvature .
we employ the wetterich equation @xcite , describing the flow of a scale - dependent effective action functional @xmath40 as a function of an ir regulator scale @xmath41 , @xmath42 = \frac{{{\ensuremath{\mathrm{i}}}}}{2 } { \operatorname{str}}\left [ \left ( \gamma_{k}^{\ensuremath{\text{\tiny{$(2)$ } } } } + r_{k } \right)^{-1 } \partial_{k } r_{k } \right ] \ !
\text{. } \label{eq : wetteq}\end{aligned}\ ] ] the effective average action @xmath40 is related to the standard generating functional for 1pi correlation functions in the limit @xmath43 . towards the uv cutoff @xmath44 ,
@xmath40 approaches the microscopic action .
the regularization is encoded in the regulator function @xmath45 , see below . for reviews of the functional rg adapted to the present context ,
see refs .
@xcite .    in the present work ,
we evaluate the flow within a rather simple approximation for the effective action . for this , we truncate the effective action to @xmath46 =   { \int \!\!\ ! { } _ { { } _ { { } _ { { } _ { { } _ \text{\small{\ensuremath{x } } } } } } } } \left [ \bar{\psi } \slashed{\nabla } \psi + \frac{\bar{\lambda}_{k}}{2 { n } } ( \bar{\psi } \psi ) ^2 \right]\text{,}\end{aligned}\ ] ] where the only scale dependence lies in the four fermion coupling @xmath47 .
furthermore , @xmath48 parametrically depends on the curvature of the background manifold .
the ir regularization is ensured by a chirally symmetric regulator of the form @xmath49 where the superscript @xmath50 denotes transposition in dirac space , and @xmath51 represents a spin - valued delta distribution , keeping track of the spinor or conjugate - spinor transformation properties associated with the spacetime arguments . in app .
[ app : beta_function ] , we briefly summarize our conventions .    for practical computations
, we use a callan - symanzik type regulator , that facilitates the use of proper time representations , @xmath52 within our investigations we restrict ourselves to negative curvature , giving rise to gravitational catalysis .
it is intuitively clear , that positive curvature ( e.g. , a sphere ) generically suppresses ir modes and thus reduces the density of states of low lying modes @xcite .
we also consider the system mainly in the large-@xmath8 limit . at finite @xmath8 ,
further pointlike fermionic self - interactions are generated which correspond to operators with an explicit curvature dependence .
we expect no qualitative modifications from these operators and hence ignore them in the following .
this approximation becomes exact in the limit @xmath53 .
it is straightforward ( cf .
[ app : beta_function ] ) to calculate the flow of the coupling as an implicit functional of the choice of the manifold , which enters via the spectrum of the dirac operator @xmath54 .
\end{aligned}\ ] ] for this calculation , it suffices to project the flow onto constant fields @xmath55 , with @xmath56 . here
, @xmath57 denotes the spacetime volume .
the operator occurring in the trace is related to the square of the regularized fermionic green s function in curved spacetime .
this has a direct correspondence to a feynman diagram representation , see fig .
[ fig : feynman_graph ] , as the flow in the present simple truncation is driven by a single fermion bubble ( and rg - improved resummations thereof ) . in the following we distinguish between the cases of a maximally symmetric spacetime in sect .
[ sec : maximally_symmetric ] which can be treated fully analytically , and a negatively curved space in sect .
[ sec : flat_time ] , which is a more interesting case in view of two - dimensional condensed matter systems .
the use of a callan - symanzik regulator shape function , cf .
( [ eq : csregulator ] ) , facilitates to rewrite the right - hand side of eq .
( [ eq : betafunction1 ] ) in terms of a simple proper time representation .
other shape functions would still permit to use a proper time representation but would lead to more intricate @xmath41 dependencies .
the laplace transform of the operator of ( [ eq : betafunction1 ] ) reads @xmath58 which upon insertion into ( [ eq : betafunction1 ] ) yields @xmath59 \text{. } \label{eq : lambdabar}\end{aligned}\ ] ]    the expression inside the super trace is known as the heat kernel @xmath60 of the ( squared ) dirac operator .
it satisfies @xmath61 and was calculated in @xcite for any maximally symmetric space , with euclidean signature . for our case
the solution to this equation can be obtained in an easier way with a special ansatz ( cf .
[ app : heat_kernel ] ) : @xmath62 where @xmath63 , and @xmath64 is the geodesic distance between the points @xmath65 and @xmath66 .
the parallel transporter @xmath67 ( wegner - wilson line ) is defined by @xmath68 where @xmath69 denotes the path ordering prescription and @xmath70 is the geodesic between @xmath65 and @xmath66 with @xmath71 and @xmath72 .
now we are able to calculate the supertrace as @xmath73 \ !
= - \frac{2 \omega { n}}{(4 \pi)^{\frac{3}{2 } } \sqrt{s } } \ !
\left ( \frac{1}{s } + { { \ensuremath{\mathrm{i}}}}\frac{{\ensuremath{\left\vertr\right\vert}}}{12 } \right ) \ ! { \mathrm{e}}^{- { { \ensuremath{\mathrm{i}}}}\frac{\pi}{4 } } \text{,}\end{aligned}\ ] ] and finally get @xmath74 in terms of the dimensionless coupling , @xmath75 we obtain the beta function @xmath76 @xmath77 this is an ordinary differential equation that parametrically depends on the curvature
. it can be solved by straightforward integration : @xmath78 the initial value is given by the dimensionless coupling @xmath79 which in terms of the initial bare gross - neveu coupling reads @xmath80 .      for the case of a manifold where the spatial part has a constant negative curvature
, we choose a special set of coordinates such that the metric can be expressed as @xmath81 where @xmath82 ( latin indices running from @xmath83 to @xmath3 ) represents the metric of a two dimensional maximally symmetric space and therefore only depends on the spatial coordinates .
hence the gamma matrices are time independent as well .
the christoffel symbols vanish for every time component @xmath84 and also the curvature tensor vanishes if any index is zero @xmath85 . from eq .
( [ eq : defspincon ] ) @xmath86 , we infer that the spin connection @xmath87 has to be proportional to @xmath18 , @xmath88.\end{aligned}\ ] ] moreover , from eq .
@xmath27 , we conclude that @xmath87 even has to vanish completely , @xmath89 this implies that the operator @xmath90 is separable into @xmath91 where the curvature is only induced by the spatial components .    for a simpler calculation of the beta function , we perform a wick rotation @xmath92 and perform again a laplace transformation of the operator occurring in eq .
( [ eq : betafunction1 ] ) , @xmath93 we arrive at the euclidean analogue of eq .   @xmath94 \text{.}\end{aligned}\ ] ] here , we calculate the super trace again with the aid of the heat kernel using that the differential operator @xmath90 is separable and the delta distribution factorizes in a time like and a spatial part @xmath95 where @xmath96 .
the quantity @xmath97 denotes the spatial coordinates of @xmath65 and should be treated as a set of coordinates and not as a vector .
both factors satisfy a heat - kernel equation and can be solved analytically , cf .
@xcite , @xmath98 where @xmath99 , @xmath100 is the parallel transporter for the spatial part and @xmath101 is the nonnegative spatial geodesic distance between @xmath97 and @xmath102 with @xmath103 .
plugging these relations into eq .
( [ eq : betafunction2 ] ) gives @xmath104 the @xmath105 integral is an integral representation of the modified bessel function of the second kind @xmath106 @xcite .
therefore , we have @xmath107 with these results we are again able to derive the beta function for the dimensionless coupling @xmath108 , @xmath109 the integration of this ordinary differential equation depending parametrically on the curvature can be cast into an integral representation , @xmath110 this result is qualitatively similar to the maximally symmetric case of eq .  .
let us now analyze the consequences of the rg flows for the long - range properties of the gross - neveu model .
for both background manifolds , @xmath76 considered as a function of @xmath111 is a parabola where the prefactor of the quadratic part is scale and curvature dependent , see fig .
[ fig : betafunction ] . for vanishing curvature
, the @xmath76 function vanishes at the two fixed points @xmath112 ( gauian ) and @xmath113 which corresponds to the well - known critical coupling of the gross - neveu model in flat space in this regularization scheme @xcite .
this critical coupling separates the symmetric phase for @xmath114 where the long range behavior is controlled by the non - interacting gauian fixed point from the chiral symmetry broken phase for @xmath115 . in the latter case
, @xmath111 runs to large values towards the infrared . in the present simple truncation
, @xmath111 in fact diverges at a finite scale @xmath116 signaling the transition into the ordered regime .
the scale @xmath116 is thus characteristic for the physical scales in the ordered phase . in @xcite
it has been shown that @xmath116 actually agrees with the value of the dynamically generated fermion mass @xmath117 as obtained in mean - field approximation .
since we are working in the large-@xmath8 limit anyway , we will use this mean - field identification in the following : @xmath118 .
the existence of the non - gauian fixed point @xmath119 can be attributed to the competition between the power - counting scaling ( the linear coupling term in @xmath76 ) and the interaction terms @xmath120 . in our rg picture ,
the interaction terms are enhanced by negative curvature as soon as the wavelength of the fluctuations becomes of the order of the curvature scale . as a consequence ,
the interacting second zero of the @xmath121 function no longer is a true fixed point but becomes scale dependent .
this `` pseudo - critical coupling '' @xmath122 moves towards the gauian fixed point for decreasing scale @xmath41 , see fig .
[ fig : betafunction ] .
any finite initial coupling strength @xmath79 will eventually become larger than @xmath123 for small rg scales @xmath124 . by this mechanism ,
the system is forced into the symmetry - broken phase even at the weakest initial coupling .
we observe this mechanism in both cases of negative curvature , the maximally symmetric as well as the purely spatial curvature case .
function of the coupling @xmath125 for different values of the scale dependent negative curvature @xmath126 ( from top to bottom : @xmath127 ) .
the black lines depict the @xmath76 functions for the maximally symmetric spacetime ( ads ) , cf . eq .
( [ eq : betafunction_maximally_symmetric ] ) ; arrows indicate the flow towards the ir .
the green dotted graphs show the flows for the case of purely spatial curvature ( lobachevsky plane ) , cf .
( [ eq : betafunction_flat_time ] ) .
in addition to the gauian fixed point , there exists a non gauian fixed point ( full red circle at @xmath128 for the present regulator scheme ) , separating the symmetric phase for @xmath129 from the broken phase for @xmath130 for zero curvature . for finite curvature , this critical point becomes scale - dependent and moves towards the gauian fixed point for increasing scale dependent curvature , i.e. , with decreasing ir scale for fixed curvature . in the case of vanishing curvature ,
the symmetry is preserved and no mass is generated for initial values @xmath131 in the blue dashed region.,title="fig:",scaledwidth=48.0% ]    as we see in fig .
[ fig : betafunction ] the influence of the curvature is somewhat stronger in the maximally symmetric case .    as discussed above
, we calculate the symmetry breaking scale @xmath116 by searching for a zero of the inverse coupling .
the fermion mass @xmath117 corresponding to this scale where the rg flow enters the symmetry - broken regime can thus be computed from the criterion @xmath132 upon partial bosonization ( hubbard - stratonovich transformation ) , @xmath133 is related to the mass parameter of a composite bosonic field .
hence the divergence of the fermionic self - interaction simply corresponds to onset of the order parameter @xcite .
let us now analyze the two different backgrounds under consideration in detail .      in the maximally symmetric case , the fermion mass defined by the criterion eq .
can be straightforwardly computed from the running coupling eq .  , yielding @xmath134 plots of this gravitationally catalyzed fermion mass are shown as a function of the curvature as solid lines in fig .
[ fig : comparison_generated_mass ] .
let us discuss this result in various limits .
in the zero - curvature limit @xmath135 , we find @xmath136 for @xmath137 . for super - critical couplings
@xmath138 , we rediscover the standard mean - field result in @xmath0 , @xmath139 this is in perfect agreement with the known behavior in flat spacetime .    provided the fermion system is initially weakly coupled , @xmath140 , a leading order expansion can be performed for any value of the curvature , resulting in @xmath141 } |r| \bar{\lambda}_\lambda } } , \quad \text{for}\,\ , \bar{\lambda}_\lambda   \lambda \ll 1,\ ] ] where we have reinserted the dimensionful initial coupling @xmath142 . if we additionally consider the limit of small curvature , we find a linear dependence of the fermion mass on both the curvature as well as the coupling , @xmath143 by contrast , in the limit of large curvature , @xmath144 , we find that @xmath145 . in other words ,
large curvature induces immediate chiral symmetry breaking , such that the induced mass becomes of the order of the cutoff .
incidentally , this result is similar for the large - coupling limit : for @xmath146 , we again find that @xmath147 to leading order independently of the curvature .    the above results display explicit uv cutoff and regularization - scheme dependencies .
since the @xmath0 gross - neveu model is asymptotically safe and thus non - perturbatively renormalizable @xcite , we can remove the uv cutoff by keeping an ir obvservable fixed while sending @xmath148 .
this `` line of constant physics '' defines a renormalized trajectory .
this can most conveniently be done in the super - critical regime where @xmath149 such that the fermion mass in flat - space @xmath150 of eq .
defines a natural ir renormalization point . at a suitable renormalization point @xmath151 to a specific value . ] in this case , the generated fermion mass in the limit @xmath152 can be written as @xmath153 where the first line holds for arbitrary curvature , and the second line represents a weak - curvature expansion being in perfect agreement with @xcite .
we emphasize that the fermions acquire a mass @xmath154 for any given @xmath79 as long as the curvature is nonvanishing .
while we expect the tendency to drive the fermion system towards the broken phase through gravitational catalysis to remain also beyond our truncation , fluctuations of bosonic composites entering beyond the large-@xmath8 limit typically provide for an opposite tendency .
hence , the status of gravitational catalysis beyond mean - field remains an interesting question .
analogously , the effects of beyond - mean - field fluctuations on magnetic catalysis are under active current investigation @xcite .      in the case of pure spatial curvature , the criterion eq .
can not be resolved analytically , but we can give an implicit equation for the induced fermion mass @xmath117 @xmath155 which can be solved numerically .
though the basic picture does not differ much from the maximally symmetric case , there are some interesting differences .
as can already be inferred from the beta function in fig .
[ fig : betafunction ] ( dotted lines ) , the curvature induced mass in the spatially curved case is smaller than in the maximally symmetric case , cf . fig . [
fig : comparison_generated_mass ] .      from top
to bottom ( @xmath128 ) .
as long as the background manifold is negatively curved , @xmath156 , a finite fermion mass is generated.,title="fig:",scaledwidth=48.0% ]    several limits can be discussed in an analytic fashion , using the series representation of @xmath157 developed in app .
[ app : curvature_expansion ] .
furthermore , we need the integral @xmath158 defined by @xmath159 where @xmath27 fixes the constant of integration . the explicit calculation is done in app .
[ app : curvature_expansion ] .
the defining equation for the fermion mass in terms of @xmath160 is @xmath161 this representation provides us with some interesting insight .
first , we can show that there exists a unique solution to this equation with @xmath162 for any given negative curvature @xmath163 and @xmath164 .
this can be seen in two steps .
the uniqueness is because the function @xmath165 for @xmath166 is bijective .
this property holds , because @xmath157 is positive and therefore @xmath160 is strictly monotonically increasing . owing to @xmath167 for small @xmath168 and @xmath169 for large @xmath170 , the function @xmath160 is not bounded . since by assumption @xmath171
, it follows from eq .   that @xmath172 has to hold , which implies that @xmath173 because of the monotonic behavior
this demonstrates that @xmath174 as claimed above .
let us first check the flat spacetime limit @xmath175 . in complete agreement with eq .
, we find again that @xmath176 = 1 - \frac{m_{\text{f},0}}{\lambda},\end{aligned}\ ] ] where we have used eq .  .    in the weak coupling regime ,
we expect the fermion mass to be small compared to the curvature scale , cf .
. hence , we need @xmath158 for large argument as provided by eq .  , @xmath177 using this asymptotic behavior for @xmath178 , we infer from eq .   that @xmath179 , \,\ ,    \alpha_\lambda= \sqrt{\frac{|r|}{2\lambda^2}}.   \end{aligned}\ ] ] if in addition , the curvature is small , i.e. , @xmath180
, we can use the corresponding expansion of eq .   and replace @xmath181 . taking differences arising from the dirac representation into account ,
the exponential inverse - coupling dependence is in perfect agreement with the results of @xcite .    despite the overall similarities to the maximally symmetric case
, we observe that the weak - coupling and weak - curvature limit of the case with pure spatial curvature shows some distinct differences .
in particular , there is an exponential non - analytic dependence of the fermion mass on the coupling as well as on the curvature .
this difference is reminiscent to magnetic catalysis in @xmath182 and @xmath183 @xcite , where the fermion gap is analytic in @xmath182 , but shows an essential singularity in the coupling in @xmath183 .
also in the present case , such a singularity shows up similar to bcs gap formation , as a consequence of the effective dimensional reduction of the fermionic fluctuation spectrum to @xmath184 @xcite .
also in this respect , our functional rg picture is in perfect agreement with previous studies @xcite .
at zero curvature , the @xmath0 gross - neveu model exhibits a critical coupling strength corresponding to a quantum critical point above which chiral symmetry is broken at large length scales .
this critical coupling manifests itself as a non - gauian fixed point of the rg flow . in the present regularization scheme , we identified @xmath185 , or @xmath128 in dimensionless conventions . as illustrated in fig .
[ fig : betafunction ] , the fixed point strictly speaking no longer exists at finite negative curvature .
the nontrivial zero of the @xmath121 function becomes scale dependent and eventually merges with the gauian fixed point in the deep ir for @xmath186 , such that only the chirally broken branch of the @xmath121 function remains .
let us call the nontrivial zero of @xmath121 a pseudo - critical coupling @xmath187 .
we find @xmath188 for the maximally symmetric case , cf .
.  , and @xmath189^{a } } \right . \!\ ! \sqrt{\frac{|r|}{2k^2 } }   \left . \vphantom{\left[\frac{a}{b}\right]^{a } } \right ) } \label{eq : lpsc}\ ] ] for the spatially curved case , cf .
. since @xmath190 is a monotonically decreasing function of scale @xmath41 , the coupling @xmath125 can eventually exceed @xmath187 such that the system becomes critical and runs towards the ordered phase . in this sense
, @xmath191 can be viewed as a criticality condition @xcite , defining a scale @xmath192 , where the system becomes critical . for lower scale ,
the system is driven towards the symmetry broken regime which is ultimately entered at @xmath193 defined above .
the value of @xmath192 depends on the curvature as well as the initial coupling @xmath79 ( the latter is considered as initially subcritical here and in the following ) .
the preceding discussion implicitly assumed that @xmath41 can run over all scales down to @xmath194 , such that the criticality condition can eventually always be satisfied .
however , if , for instance , the system has a finite volume characterized by a finite length scale @xmath195 , also the fluctuation momenta are restricted , typically leading to an ir cutoff @xmath196 .
one may think of a finite probe length , such as , e.g. , the size of a layer of graphene .
this finite probe length @xmath195 can lead to a screening of the gravitationally catalyzed ordered regime if @xmath197 is larger than the would - be critical scale @xmath192 .
hence , @xmath198 can be viewed as a lower bound for the coupling strength required for symmetry breaking in a real system of finite length .
it thus generalizes the critical coupling at infinite volume and zero curvature to the situation of finite volume and finite curvature .
it is instructive to study @xmath187 as a function of the length scale @xmath195 measured in units of a typical curvature length scale which we define by @xmath199 . for finite systems , this gives an estimate for how strongly a probe has to be curved in order to exhibit gravitational catalysis . in turn , for a given curvature of the probe , @xmath187 provides an estimate for the initial coupling strength required for symmetry breaking in the finite system .
for instance , for the maximally symmetric case , we read off from eq .
that @xmath200 for the spatially curved case , we incidentally find the same result in the limit of small curvature , i.e. , @xmath201 , using the expansion which is accurate even for values of @xmath202 .
by contrast , the large curvature limit for the spatially curved case is different , cf .
: @xmath203 from eq .
, it is obvious that probe length to curvature ratios up to @xmath204 lead to pseudo - critical couplings @xmath187 which deviate from the zero - curvature critical coupling of the gross - neveu model only below the 1% level .
significant deviations only occur for @xmath195 being an order of magnitude larger than the curvature scale @xmath205 . from the viewpoint of curvature - deformed condensed matter systems ,
a ratio of @xmath204 appears to be `` large '' in the sense that
loosely speaking  a spatially curved @xmath206 planar probe embedded in @xmath0 euclidean space would rather look like a @xmath0 object .
another way to interpret these results is the following : consider a finite - probe system with a subcritical bare coupling , @xmath207 , thus being in the symmetric ( e.g. semimetal ) phase . in order to gravitationally
catalyze a transition to a broken ( gapped or insulating ) phase , the criticality condition has to be met for a sufficiently large @xmath208 . in view of eq .
this requires comparatively strong curvature , i.e. , a small curvature length scale @xmath205 compared with the probe length @xmath195 .
we have investigated the phenomenon of gravitational catalysis in the @xmath0 gross - neveu model on specific manifolds with constant negative curvature .
while the mechanism had already been studied frequently with mean - field methods as well as from the viewpoint of the fluctuation spectrum of the dirac operator , we have added a new renormalization group picture to the comprehensive understanding of this phenomenon .
the essence of this picture is that the critical coupling of the fermionic system , corresponding to a quantum phase transition in flat spacetime , is transmuted into a scale - dependent pseudo - critical coupling that flows to zero as a consequence of long - wavelength fluctuations ( compared to the curvature scale ) . in this manner ,
the infinite - volume fermion system becomes critical for any arbitrarily weak coupling .
we have identified the rg ( pseudo- ) fixed point mechanism for two example manifolds of constant negative curvature : a maximally symmetric spacetime ( anti de sitter ) and a purely spatially curved case ( lobachevsky plane ) .
both manifolds support the mechanism of gravitational catalysis , but exhibit a rather different behavior as far as the dependence of chiral symmetry breaking on the coupling and the curvature are concerned .
the maximally symmetric case shows a linear dependence ( to leading order ) on both quantities which makes clear that the phenomenon is essentially perturbative , being reminiscent of the `` quantum anomaly '' for fermions in a magnetic field @xcite .
by contrast , order parameters indicating the symmetry broken state such as the induced fermion mass exhibit an essential singularity in both the coupling and the curvature for the case of the purely spatially curved case .
this is in many respects similar to bcs - type gap formation .
again , our renormalization group picture goes hand in hand here with properties of the fermionic fluctuation spectrum , as analyzed in @xcite .    as a benefit
, the functional renormalization group also gives a simple access to systems of finite extent by identifying the rg infrared cutoff with an inverse length scale . in this manner
, we can estimate the fate of gravitational catalysis in finite systems .
in fact , the phenomenon only occurs , as long as the curvature radius is sufficiently small compared to the systems length scale .
we have been able to phrase this statement quantitatively by introducing a pseudo - critical coupling .
thinking in terms of curved layered condensed matter systems , rather large curvatures are needed compared to a realizable probe length in order to drive a sub - critical system into a phase dominated by gravitational catalysis .
we would like to emphasize that an immediate application of our results to condensed matter systems would only be possible for reparametrization invariant systems such as fluid membranes @xcite , curvature effects of which can be mapped onto the language of riemannian geometry . for tethered membranes or general lattice systems , further phenomena connected to extrinsic curvature or curvature related defects can be become relevant .
in this context , it is interesting to note that an external strain exerted on a graphene sheet in flat space induces a pseudo - magnetic field @xcite , that may also support ( pseudo-)magnetic catalysis . for such systems , we hence expect an interesting interplay between these various effects if we expose them to negative curvature inducing strain .
finally , gravitational catalysis may become relevant in the context of asymptotically safe quantum gravity @xcite . in conjunction with fermionic degrees of freedom @xcite , the uv fixed point determining the shape of the universe at highest energies might go along with a negative ( though scale - dependent ) curvature @xcite . whether or not gravitational catalysis in connection with gravitionally modified critical fermion interactions @xcite could become active and impose constraints on the matter content of the universe
then is a highly involved question that deserves to be investigated in greater depth .
the authors thank astrid eichhorn , lukas janssen , daniel scherer , ren sondenheimer , and andreas wipf for valuable discussions and acknowledge support by the dfg under grants gi  328/5 - 2 ( heisenberg program ) , grk1523 and for 723 .
let us first introduce our conventions . in the weldon formalism @xcite ,
the dirac conjugate spinor @xmath209 is related to the hermitean conjugate of the spinor @xmath210 via the spin metric @xmath211 @xmath212 which is implicitly defined by @xmath213 this gives rise to the following properties of the building blocks of the gross - neveu action : @xmath214 choosing the spin metric @xmath211 to be anti - hermitean in eq .
( i ) is convenient for our metric and clifford algebra conventions , but differs from @xcite . also , the absence of any imaginary factor of `` @xmath215 '' in the fermion kinetic term is due to these conventions .
the flow equation uses a rather condensed notation .
more concretely , we work in field space parameterized by the collective fields @xmath216 representing grassmann - valued functions on the manifold , reminiscent to nambu - gorkov spinors . for instance , the classical action amended by the regulator term which is used for deriving the flow equation from the functional integral ( cf .
@xcite ) , reads in these conventions @xmath217 = { } & s[\phi ] + \frac{1}{2 } { \int \!\!\ ! { } _ { { } _ { { } _ { { } _ { { } _ \text{\small{\ensuremath{x } } } } } } } } \!\!\ ! { \int \!\!\ ! { } _ { { } _ { { } _ { { } _ { { } _ \text{\small{\ensuremath{y } } } } } } } } \bar{\phi}(x ) r_{k}(x , y ) \phi(y).\end{aligned}\ ] ] these conventions differ slightly from those commonly used , see , e.g. , @xcite , but turn out to be advantageous for coordinate space computations on curved manifolds .
for instance , the representation of the unit element in field space becomes rather intuitive , @xmath218 by @xmath51 , we denote the spin - valued delta distribution , which fulfills @xmath219 therefore , @xmath220 has to hold . with this notation
, we indicate that @xmath51 transforms as a spinor in @xmath65 and a dirac - conjugated spinor in @xmath66 , i.e. @xmath51 can be interpreted as @xmath221 , with the standard scalar delta distribution @xmath222 .
we remark that there is a difference between @xmath51 and @xmath223 in the spinor structure , since @xmath224 for the evaluation of the flow equation eq .
, we proceed in a standard fashion .
we decompose @xmath225 into a field - dependent part @xmath226 and a field - independent part @xmath227 in order to expand the wetterich equation in powers of the fields , @xmath228 \text{.}\end{aligned}\ ] ] since @xmath229 & \psi \psi^{{\ensuremath{\mathrm{t } } } } \\
\bar{\psi}^{{\ensuremath{\mathrm{t } } } } \bar{\psi } & { } \!\!\ ! [ ( \bar{\psi } \psi ) { \ensuremath{\mathrm{i}}}+ \psi \bar{\psi}]^{{\ensuremath{\mathrm{t } } } } \ ! { }   \end{pmatrix } \ !
\mathbb{1},\\   { \ensuremath{\mathcal{p}}}_{k } = { } & { { \ensuremath{\mathrm{i}}}}k \ !
\begin{pmatrix }      \sqrt{\tau}\big ( { \ensuremath{\mathrm{i}}}+ r(\tau ) \big ) & 0 \\      0 & \sqrt{\tau^{{\ensuremath{\mathrm{t } } } } } \big ( { \ensuremath{\mathrm{i}}}+ r(\tau^{{\ensuremath{\mathrm{t } } } } ) \big )   \end{pmatrix } \ !
\mathbb{1},\end{aligned}\ ] ] we observe that only the term @xmath230 in eq
.   can contribute to the flow of @xmath47 .
we need @xmath231 with @xmath232 , and therefore only the diagonal of @xmath233 is required . in the limit @xmath234 ,
only the following terms remain : @xmath235_{11}(x , y ) = - \frac{\frac{\bar{\lambda}_{k}^2}{{n}^{2 } k^{2 } } ( \bar{\psi } \psi)^2}{\tau \big ( { \ensuremath{\mathrm{i}}}+ r(\tau ) \big)^{2 } } \ , \frac{\delta(x , y)}{\sqrt{-g}},\\   & \left[({\ensuremath{\mathcal{p}}}_{k}^{-1 } { \ensuremath{\mathcal{f}}}_{k})^2\right]_{22}(x , y ) = \ ! \left[({\ensuremath{\mathcal{p}}}_{k}^{-1 } { \ensuremath{\mathcal{f}}}_{k})^2\right]^{{\ensuremath{\mathrm{t}}}}_{11}(y , x)\text{. } \ ] ] the lhs of eq
.   boils down to @xmath236 and the rhs yields @xmath237 \notag \\   & { } = { { \ensuremath{\mathrm{i}}}}\frac{2 \bar{\lambda}_{k}^2}{{n}^2 k^3 } ( \bar{\psi } \psi)^2 { \operatorname{str}}\ ! \left [ \frac { r'(\tau ) } { \big ( { \ensuremath{\mathrm{i}}}+ r(\tau ) \big)^{3 } } \frac{\delta(x , y)}{\sqrt{-g } } \right].\end{aligned}\ ] ] inserting the callan - symanzik regulator , we end up with eq .   of the main text .
following @xcite , we choose for the heat kernel satisfying eq .
the ansatz @xmath238 where @xmath239 is a scalar function of @xmath240 the geodesic distance and the proper time @xmath105 . plugging this into equation ( [ eq : heatkernel ] ) using @xmath241 , @xmath242 and @xmath243 n^{\nu } u \text { , } \end{aligned}\end{aligned}\ ] ] cf .
@xcite , we get @xmath244 with @xmath245 , @xmath246 . because of the boundary condition for @xmath247 and the regularity of @xmath248 we observe that @xmath249 has to hold , with @xmath222 representing the standard delta distribution . one suitable representation for the delta distribution in 3@xmath13 maximally symmetric space with negative curvature in the limit @xmath250 is @xmath251 next , we factorize @xmath239 into the delta part and an auxiliary function @xmath252 of @xmath240 and @xmath105
, @xmath253 expanding @xmath252 in powers of @xmath254 and using the boundary conditions , @xmath255 leads to @xmath256 with the boundary condition @xmath257 . from equation @xmath86 and the boundary condition ,
we see @xmath258 plugging this into @xmath27 gives @xmath259 where we have eliminated the constant of integration using equation @xmath260 . from this , we finally get @xmath261
for the detailed analysis of the spatially curved case , the integral representation of the running coupling eq .   can be studied in various limits analytically .
more specifically , @xmath157 as defined in eq
.   can be expanded for small and for large values of @xmath262 .
let us first consider an expansion of this function about @xmath263 , starting with an expansion of the integrand , @xmath264 using the standard integral @xmath265 the small @xmath262 expansion of @xmath157 can be computed to any order . to order @xmath266 , we find @xmath267 due to the factorial growth of the coefficients , cf . eq .
( [ eq : bessel_integral ] ) , the expansion is an asymptotic series . by comparison with the numerical result , the accuracy of eq .
( [ eq : mfri_small_alpha ] ) turns out to be above 99% up to @xmath268 .
a similar approximation can be done for large @xmath262 by expanding @xmath269 which holds for any @xmath270 .
next , we employ @xmath271 which can be used for @xmath272 . for @xmath273 , an analytic continuation into the complex
is implicitly understood , leading to a replacement of the term @xmath274 by @xmath275 \big ( 1 - ( n \alpha)^2 \big)^{-3/2}$ ] here and in the following .
this leads to the convergent series @xmath276\end{aligned}\ ] ] for any @xmath277 . neglecting orders higher than @xmath278 and using @xmath279 for @xmath280
, we arrive at @xmath281 where @xmath282 is the euler - mascheroni constant and @xmath283 is the glaisher - kinkelin constant .
the accuracy of this result is above 99% for @xmath284 .    with the series ( [ eq : mfri_sum_rep ] )
it is possible to find a series representation for the required integral ( [ eq : def_mfrf ] ) @xmath285 \ ! \text{.}\end{aligned}\ ] ] this series is convergent for any @xmath286 .
the expansions for small and large @xmath262 are @xmath287 where @xmath288 is defined as the limit @xmath289      v.  p.  gusynin , v.  a.  miransky and i.  a.  shovkovy , phys .
d * 52 * , 4718 ( 1995 ) [ hep - th/9407168 ] .
v.  p.  gusynin , v.  a.  miransky and i.  a.  shovkovy , phys .
b * 349 * , 477 ( 1995 ) [ hep - ph/9412257 ] .
v.  p.  gusynin , v.  a.  miransky and i.  a.  shovkovy , nucl .
b * 462 * , 249 ( 1996 ) [ hep - ph/9509320 ] .
i.  a.  shovkovy , arxiv:1207.5081 [ hep - ph ] .
i.  a.  shushpanov and a.  v.  smilga , phys .
b * 402 * , 351 ( 1997 ) [ hep - ph/9703201 ] .
t.  d.  cohen , d.  a.  mcgady and e.  s.  werbos , phys .
c * 76 * , 055201 ( 2007 ) [ arxiv:0706.3208 [ hep - ph ] ] .
a.  v.  zayakin , jhep * 0807 * , 116 ( 2008 ) [ arxiv:0807.2917 [ hep - th ] ] .
a.  j.  mizher , m.  n.  chernodub and e.  s.  fraga , phys .
d * 82 * , 105016 ( 2010 ) [ arxiv:1004.2712 [ hep - ph ] ] .
r.  gatto and m.  ruggieri , phys .
d * 82 * , 054027 ( 2010 ) [ arxiv:1007.0790 [ hep - ph ] ] .
j.  k.  boomsma and d.  boer , phys .
d * 81 * , 074005 ( 2010 ) [ arxiv:0911.2164 [ hep - ph ] ] .
g.  s.  bali , f.  bruckmann , g.  endrodi , z.  fodor , s.  d.  katz and a.  schafer , phys .
d * 86 * , 071502 ( 2012 ) [ arxiv:1206.4205 [ hep - lat ] ] .
g.  w.  semenoff , i.  a.  shovkovy and l.  c.  r.  wijewardhana , mod .
lett .  a * 13 * , 1143 ( 1998 ) [ hep - ph/9803371 ] .
k.  krishana , n.  p.  ong , y.  zhang , _ et al .
_ , phys .
lett . * 82 * , 5108 ( 1999 ) [ cond - mat/9904138 ] .
d.  v.  khveshchenko , phys .
lett .   * 87 * , 206401 ( 2001 ) [ cond - mat/0106261 ] .
h.  leal and d.  v.  khveshchenko , nucl .
b * 687 * , 323 ( 2004 ) [ cond - mat/0302164 ] .
i.  f.  herbut and b.  roy , phys .
b * 77 * , 245438 ( 2008 ) [ arxiv:0802.2546 [ cond-mat.mes-hall ] ] ; b.  roy and i.  f.  herbut , phys .
b * 83 * , 195422 ( 2011 ) [ arxiv:1102.3481 [ cond-mat.mes-hall ] ] .
e.  j.  ferrer and v.  de la incera , phys .
* 102 * , 050402 ( 2009 ) [ arxiv:0807.4744 [ hep - ph ] ] ; nucl .
b * 824 * , 217 ( 2010 ) [ arxiv:0905.1733 [ hep - ph ] ] ; e.  j.  ferrer , v.  de la incera and a.  sanchez , phys .
lett .   * 107 * , 041602 ( 2011 ) [ arxiv:1103.5152 [ hep - ph ] ] .
d.  d.  scherer and h.  gies , phys .
b * 85 * , 195417 ( 2012 ) [ arxiv:1201.3746 [ cond-mat.str-el ] ] .
k.  fukushima and j.  m.  pawlowski , phys .
d * 86 * , 076013 ( 2012 ) [ arxiv:1203.4330 [ hep - ph ] ] .
i.  l.  buchbinder and e.  n.  kirillova , sov .
j.   * 32 * , 446 ( 1989 ) ; int .
j.  mod .
phys .  a * 4 * , 143 ( 1989 ) .
i.  l.  buchbinder , s.  d.  odintsov and i.  l.  shapiro , `` effective action in quantum gravity , '' bristol , uk : iop ( 1992 ) .
t.  inagaki , t.  muta and s.  d.  odintsov , mod .
lett .  a * 8 * , 2117 ( 1993 ) [ hep - th/9306023 ] .
i.  sachs and a.  wipf , phys .
b * 326 * , 105 ( 1994 ) [ hep - th/9310085 ] .
e.  elizalde , s.  leseduarte , s.  d.  odintsov and y.  i.  shilnov , phys .
d * 53 * , 1917 ( 1996 ) [ hep - th/9505065 ] .
s.  kanemura and h.  -t .
sato , mod .
a * 11 * , 785 ( 1996 ) [ hep - th/9511059 ] .
t.  inagaki , int .
j.  mod .
a * 11 * , 4561 ( 1996 ) [ hep - th/9512200 ] .
t.  inagaki and k.  -i .
ishikawa , phys .
d * 56 * , 5097 ( 1997 ) .
b.  geyer and s.  d.  odintsov , phys .
b * 376 * , 260 ( 1996 ) [ hep - th/9603172 ] ; phys .
d * 53 * , 7321 ( 1996 ) [ hep - th/9602110 ] .
g.  miele and p.  vitale , nucl .
b * 494 * , 365 ( 1997 ) [ hep - th/9612168 ] ; p.  vitale , nucl .
b * 551 * , 490 ( 1999 ) [ hep - th/9812076 ] .
t.  inagaki , t.  muta and s.  d.  odintsov , prog .
suppl .   * 127 * , 93 ( 1997 ) [ hep - th/9711084 ] .
j.  hashida , s.  mukaigawa , t.  muta , k.  ohkura and k.  yamamoto , phys .
d * 61 * , 044015 ( 2000 ) [ gr - qc/9907014 ] .
e.  v.  gorbar and v.  p.  gusynin , annals phys .   * 323 * , 2132 ( 2008 ) [ arxiv:0710.2292 [ hep - ph ] ] .
m.  hayashi , t.  inagaki and h.  takata , arxiv:0812.0900 [ hep - ph ] ; t.  inagaki and m.  hayashi , arxiv:1003.1173 [ hep - ph ] .
s.  sasagawa and h.  tanaka , prog .
.   phys .   *
128 * , 925 ( 2012 ) [ arxiv:1209.2782 [ hep - ph ] ] .
e.  v.  gorbar , phys .
d * 61 * , 024013 ( 2000 ) [ hep - th/9904180 ] ; ukr .
j.  phys .
* 54 * , 541 ( 2009 ) [ arxiv:0809.2558 [ hep - th ] ] .
d.  ebert , a.  v.  tyukov and v.  c.  zhukovsky , phys .
d * 80 * , 085019 ( 2009 ) [ arxiv:0808.2961 [ hep - th ] ] .
d.  v.  khveshchenko , phys .
lett .   * 87 * , 246802 ( 2001 ) [ cond - mat/0101306 ] .
i.  f.  herbut , phys .
lett .   * 97 * , 146401 ( 2006 ) [ cond - mat/0606195 ] .
d.  j.  gross and a.  neveu , phys .
d * 10 * , 3235 ( 1974 ) .
h.  a.  weldon , phys .
d * 63 * , 104010 ( 2001 ) [ gr - qc/0009086 ]
. f.  hofling , c.  nowak and c.  wetterich , phys .
b * 66 * , 205111 ( 2002 ) [ cond - mat/0203588 ] .
l.  janssen and h.  gies , phys .
d * 86 * , 105007 ( 2012 ) [ arxiv:1208.3327 [ hep - th ] ] ; h.  gies and l.  janssen , phys .
d * 82 * , 085018 ( 2010 ) [ arxiv:1006.3747 [ hep - th ] ] .
l.  janssen , `` critical phenomena in ( 2 + 1)-dimensional relativistic fermion systems , '' phd thesis , fsu jena , http://www.db-thueringen.de/servlets/document + servlet?id=20856 ( 2012 )
. j.  braun , h.  gies and d.  d.  scherer , phys .
d * 83 * , 085012 ( 2011 ) [ arxiv:1011.1456 [ hep - th ] ] . c.  wetterich , phys .
b * 301 * , 90 ( 1993 ) .
j.  berges , n.  tetradis and c.  wetterich , phys .
rept .  * 363 * , 223 ( 2002 ) [ hep - ph/0005122 ] .
k.  aoki , int .
j.  mod .
b * 14 * , 1249 ( 2000 ) .
j.  m.  pawlowski , ann .
* 322 * , 2831 ( 2007 ) [ hep - th/0512261 ] .
h.  gies , lect .
notes phys .   *
852 * , 287 ( 2012 ) [ hep - ph/0611146 ] .
b.  delamotte , lect .
notes phys .
* 852 * , 49 ( 2012 ) [ cond - mat/0702365 [ cond - mat ] ] .
p.  kopietz , l.  bartosch and f.  schutz , lect .
notes phys .
* 798 * , 1 ( 2010 ) .
w.  metzner , m.  salmhofer , c.  honerkamp , v.  meden and k.  schonhammer , rev .
phys .   * 84 * , 299 ( 2012 ) [ arxiv:1105.5289 [ cond-mat.str-el ] ] .
j.  braun , j.  phys .
g * 39 * , 033001 ( 2012 ) [ arxiv:1108.4449 ] .
r.  camporesi , commun .
phys .   * 148 * , 283 ( 1992 ) .
i.  s.  gradshteyn and i.  m.  ryzhik , `` table of integrals , series , and products , '' new york , us : academic press ( 2000 ) .
j.  jaeckel and c.  wetterich , phys .
d * 68 * , 025020 ( 2003 ) [ hep - ph/0207094 ] .
u.  ellwanger and c.  wetterich , nucl .
b * 423 * , 137 ( 1994 ) [ hep - ph/9402221 ] .
h.  gies and c.  wetterich , phys .
d * 65 * , 065001 ( 2002 ) [ hep - th/0107221 ] .
v.  skokov , phys .
d * 85 * , 034026 ( 2012 ) [ arxiv:1112.5137 [ hep - ph ] ] .
j.  braun , c.  s.  fischer and h.  gies , phys .  rev .
d * 84 * , 034045 ( 2011 ) [ arxiv:1012.4279 [ hep - ph ] ] .
v.  p.  gusynin and s.  g.  sharapov , phys .
lett .   * 95 * , 146801 ( 2005 ) [ cond - mat/0506575 ] .
d.  nelson , t.  piran and s.  weinberg , `` statistical mechanics of membranes and surfaces , '' singapore , sg : world scientific ( 1989 ) .
f.  guinea , m.  i.  katsnelson and a.  k.  geim , nature phys .
* 6 * , 30 ( 2010 ) [ arxiv:0909.1787 [ cond-mat.mes-hall ] ] ; f.  guinea , a.  k.  geim , m.  i.  katsnelson , k.  s.  novoselov , phys .
b * 81 * , 035408 ( 2010 ) ; n. levy _ et al .
_ , science * 329 * , 544 ( 2010 ) .
s.  weinberg , _ in * hawking , s.w .
, israel , w. : general relativity * , 790 - 831 _ ; m.  reuter , phys .
d * 57 * , 971 ( 1998 ) [ arxiv : hep - th/9605030 ] ; m.  niedermaier and m.  reuter , living rev .
* 9 * , 5 ( 2006 ) ; r.  percacci , in * oriti , d. ( ed . ) : approaches to quantum gravity * 111 - 128 [ arxiv:0709.3851 [ hep - th ] ] ; m.  reuter and f.  saueressig , new j.  phys .   *
14 * , 055022 ( 2012 ) [ arxiv:1202.2274 [ hep - th ] ] .
dona and r.  percacci , arxiv:1209.3649 [ hep - th ] ; u.  harst and m.  reuter , jhep * 1205 * , 005 ( 2012 ) [ arxiv:1203.2158 [ hep - th ] ] ; r.  percacci and d.  perini , phys .
d * 67 * , 081503 ( 2003 ) [ hep - th/0207033 ] .
o.  lauscher and m.  reuter , jhep * 0510 * , 050 ( 2005 ) [ hep - th/0508202 ] .
a.  eichhorn and h.  gies , new j.  phys .   *
13 * , 125012 ( 2011 ) [ arxiv:1104.5366 [ hep - th ] ] .