this paper addresses the problem of interacting electrons in a band of a two - dimensional crystal .
this problem uncovers all the macroscopic properties of the materials of which the maximum interest is centered at the present on the high-@xmath1 superconductors , in particular , on the nature of their normal state . to fix ideas and
some notation let us review the highlights of the problem .
start with a free theory of , say , electrons in a metal , i. e. , a free fermi gas whose ground state is known and whose low energy elementary excitations are electrons and holes obeying a free dispersion relation . then switch on an interaction among the electrons and ask the question of what the new ground state is and what is the nature of the low energy excitations of the full interacting theory .
this program has been successfully addressed by landau in three spatial dimensions @xmath2 where it leads to the famous theory of landau s fermi liquid @xcite and the concept of quasiparticles .
things are different in @xmath3 where an exact solution can be found and the behavior of the system is always of the luttinger type@xcite .    the case @xmath4 is physically very interesting both in practice ( most of the high-@xmath1 materials are organized in two dimensional layers and experimental data are available ) and under a theoretical point of view . despite the great effort devoted to it
, no exact solution is known at the present and the situation remains controversial superconductivity . ] .
the ultimate origin of the trouble lies on the very complex nature of the many - body vacuum and the corresponding difficulties of doing perturbations around the fermi surface .
the theory of the fermi liquid and the perturbations of it that lead to the bcs and charge density wave instabilities have recently been rephrased as effective field theories @xcite to which renormalization group techniques can be applied @xcite . the main idea of this approach is that , although we are dealing with problems that naturally have an ultraviolet cutoff @xmath5 ( the width of the band ) , we are only interested in the physics at a scale much lower than it , namely , a very narrow strip above and below the fermi surface . by integrating out in the path integral all the momenta above a certain cutoff @xmath6 , we end up with an effective field theory whose couplings are cutoff dependent .
the renormalization group transformations are rescalings of @xmath7 . in this way one can visualize the fermi liquid as a fixed point of this rg transformations ( it is gapless ) that flows towards a bcs state when a weak attraction is introduced between the quasiparticles .    when comparing the effective field theories that arise in this context with standard quantum field theory , we encounter two main differences .
the first one is the nonrelativistic nature of the many body physics where there is a parameter , the fermi velocity , whose typical value is much smaller than the speed of light .
this prevents the use of the general theorems concerning renormalizability .
the more severe difference lies in the above mentioned nature of the many body vacuum . when treating the interacting system as a perturbation of the  free " system constituted by the fixed point , perturbations are to be described around a fermi surface in @xmath8 , a fermi line in @xmath9 , and fermi points in @xmath3 .
we are then faced with couplings that are functions of the shape of the fermi surface which is also modified by the interaction . only in the @xmath3 case
we deal , as in quantum field theory , with a true vacuum around which one can expand and get a bunch of coupling _
constants_.    in ref .
@xcite , we developed a model for the electronic structure of free electrons in a two - dimensional honeycomb lattice .
this lattice constitutes the basic building block of periodic systems made of identical atoms with planar , threefold coordination .
that is the case of graphite , built of carbon atoms .
the rapidly growing family of fullerene compounds @xcite can be thought of as closed , curved surfaces derived from the planar honeycomb lattice .
these compounds show a variety of unexpected features , the most striking of them being the superconductivity of doped crystals of c@xmath10 molecules @xcite .
moreover , this superconductivity can be considered a form of high - t@xmath11 superconductivity , as the number of electrons available to form cooper pairs is unusually low .
the main feature that makes the honeycomb lattice the object of our attention in the present context is that its fermi level at half - filling consists of exactly two fermi _
points_. this peculiarity opens the possibility of performing a complete analysis of the problems mentioned above in analogy with the @xmath3 case .
in particular it makes of this lattice an ideal laboratory for investigating departures from fermi liquid behavior in @xmath9 which have shown very ellusive in more conventional approaches .
the major interest of this investigation relies on the fact , made clear by now from the experiments , that the normal state of high-@xmath1 materials is not of fermi liquid type . besides
, it can be directly applied to the physics of the fullerenes that , as mentioned above , are two - dimensional graphitic structures at half - filling @xcite .
the paper is organized as follows . in section 2
we review results on the free electrons on the honeycomb lattice that lead to the low energy effective action used as a starting point of our analysis .
we emphasize the point that , unlike what usually occurs in many body physics , our effective action constitutes a genuine quantum field theory with a single , well - defined , vacuum state .
next the free theory is completed with the addition of the coulomb interaction . in the field theory , this is mediated by an electromagnetic field whose propagator reproduces , in the nonrelativistic limit , the usual four - fermion interaction considered in the many - body approach .
we end up with a standard local quantum field theory .
a discussion is done on the differences to be expected by the use of the relativistic versus the instantaneous electromagnetic field propagator illustrated by the computation of the density of states in both cases .
the results announce an infrared instability which is clarified in section 4 with the help of the renormalization group .
section 3 is devoted to the study of the renormalizability of the model .
it is shown that , despite its nonrelativistic nature and due to its gauge invariance , it is renormalizable .
the fermi velocity and the electron wave function get renormalized while the electric charge does not .
sections 4 and 5 contain the analysis of the rg flow of the model and its physical implications .
the first section is devoted to the nonrelativistic limit of the theory .
we see that the perturbative expansion has an effective coupling given by the ratio of the square of the electric charge and the fermi velocity .
the rg shows that this coupling decreases to zero in the infrared although there is no fixed point in the nonrelativistic regime .
we then investigate the possibility of having a natural length scale in the theory and make some comments on the system away from half - filling .
finally we find an anomalous dimension for the electron wave function at two - loop order that points towards a non - fermi liquid behavior of the system .
the previous analysis is completed in section 5 where we study the relativistic regime which is the relevant regime in the infrared .
we find a fixed point of the rg flow obtained when the fermi velocity equals the speed of light .
next we compare our model with standard rg analysis of the fermi liquid theory and find significant differences both from the technical and from the physical point of view .
the most significant one is the anomalous dimension found for the electron green function .
we discuss the non renormalization of the electric charge and establish the non - fermi liquid nature of our fixed point .
we set our main conclusions and a summary of the more relevant points of the article in section 6 .
in this section we first review the one - particle electronic states in the two - dimensional honeycomb lattice , which serve afterwards as a starting point for the many - body problem .
the present one - particle description has been worked out with complete detail in ref .
@xcite , regarding the conduction band of a two - dimensional graphite layer .
the honeycomb lattice has two different atoms per primitive cell , so that two degrees of freedom arise in a variational computation of the energy eigenstates , corresponding to the respective orbitals for the two atoms . in the tight - binding approximation
the problem reduces to the diagonalization of the one - particle hamiltonian @xmath12 where the sum is over pairs of nearest neighbors atoms @xmath13 on the lattice and @xmath14 , @xmath15 are canonically anticommuting operators @xmath16 the states which are simultaneously eigenvectors of the hamiltonian and of the lattice generators have to be of the form @xmath17 with coefficients @xmath18 and @xmath19 for black and blank points , respectively , as depicted in figure 1 .
obviously , under the action of ( [ 6 ] ) black points are mapped into blank points and viceversa .
we have , indeed , @xmath20 where @xmath21 is a triad of vectors connecting a @xmath22 atom with its nearest neighbors , and @xmath23 the triad made of their respective opposites ( see figure 1 ) .
therefore , the state ( [ 8 ] ) is an eigenvector of @xmath24 provided that the coefficients @xmath18 and @xmath19 are solutions of the eigenvalue problem @xmath25 the parameter @xmath26 measures the hopping between nearest neighbors and provides the natural energy scale in the one - particle description .
the diagonalization of ( [ 9 ] ) gives the band of levels @xmath27 the parameter @xmath28 being the lattice spacing .
we will be interested in the consideration of this band at half - filling , which is the pertinent instance for the carbon - based materials and , in particular , graphite layers and undoped fullerenes .
the dispersion relation ( [ 12 ] ) has the amazing property that the fermi level at @xmath29 is only reached by six isolated points in momentum space ( see figure 2 ) .
these are the six corners of the first brillouin zone and , due to the periodicity of the dual lattice , they correspond to only two independent states .
the existence of a finite number of fermi points is quite unusual in more than one spatial dimension .
it has important consequences since it makes possible to encode the low - lying excitations about the fermi level into a simple field theory .
when considering an arbitrarily large sample , it is appropriate to amplify the region about any of the two independent fermi points by taking the limit @xmath30 . in either case
we obtain for the operator @xmath31 in equation ( [ 9 ] ) @xmath32 where is the displacement about the fermi point and @xmath33 stand for the two pauli matrices .
the conclusion is that , in the long - wavelength limit @xmath30 , the one - particle ( energy ) eigenstates of the electrons in the honeycomb lattice are given by the field theory of two massless dirac spinors ( one for each independent fermi point ) in two spatial dimensions .
this is all about the free theory .
we now turn to introduce the electronic interaction by applying the method of second quantization .
we remark that , proceeding in this way with the above set of one - particle states , we end up with a formalism which falls more into the framework of quantum field theory than into that of quantum statistical theory .
the difference is subtle but significant . in quantum statistical theory
the ground state of the system is , in general , a condensate with an extended fermi surface ( in more than one spatial dimension ) and the elementary excitations ( quasiparticles ) are only stable right at the fermi level .
the second quantized theory built from the dispersion relation ( [ deq ] ) has , however , a definite vacuum in the sense of quantum field theory .
its excitations are similar to the particle - hole excitations of the dirac sea .
we may consider separately , for instance , each of the two dirac spinors and work out the many - body electron propagator of the free theory @xmath34 we denote by @xmath35 the electron operator @xmath36 in the interaction representation while @xmath37 . in terms of the respective modes @xmath38 for unoccupied and occupied levels we have @xmath39 where @xmath40 is the fermi velocity and the operators @xmath41 satisfy anticommutation relations @xmath42 taking the statistical average with the system at half - filling , we get for the fourier transform of the propagator @xmath43 @xmath44 is a standard set of @xmath45-matrices satisfying @xmath46 and related to the above quoted pauli matrices by @xmath47 .
the above expression coincides precisely with the feynman propagator for a massless spinor in 2 + 1 dimensions , which confirms our assertion that the system at half - filling is equivalent to the dirac sea in quantum field theory .
we must notice however the presence of the fermi velocity @xmath48 in front of the spatial part of the scalar product , which is a feature of the many - body theory .    in the nonrelativistic approximation ,
the introduction of the electronic interaction leads , in first instance , to the second quantized hamiltonian @xmath49 we attach to this hamiltonian the name of `` coulomb '' since it describes the instantaneous interaction between electric charges . the model given in terms of the coulomb interaction is , however , a highly nonlocal field theory , what makes very awkward its investigation from a formal point of view . both for the sake of studying the properties of the quantum theory as well as for computational purposes , a complete description of the interaction with the electromagnetic field @xmath50
is more desirable .
one of the points that we support in what follows is that , actually , the model of interacting electrons is sensitive to retardation effects of the electromagnetic propagation .
the reason for such unconventional behaviour is the massless character of the spinor @xmath51 .
this property is not accidental in the description of the two - dimensional layer , as long as the expansion ( [ deq ] ) gives only one marginal operator dictating the long - distance behaviour of the free theory .
we propose , therefore , a quantum field theory description based on the second quantized hamiltonian .
] @xmath52 the interaction of the electromagnetic field @xmath53 and the electrons in the layer is described in the standard fashion , by coupling to the conserved current @xmath54 this poses some technical problems since the electromagnetic field propagates in three - dimensional space while we want the dynamics of the electrons to be confined to the two - dimensional layer .
although this may not be achieved in general , it turns out to be possible by specializing to the feynman gauge@xcite , which enforces the constraint @xmath55 and has also the property of placing the @xmath56 field in the same direction as the electronic current . in this gauge
we have the propagator @xmath57 and the coupling to the 2 + 1 dimensional current is perfectly consistent . in the computation of electronic properties
we have just to take care of placing the points @xmath58 on the plane of the two - dimensional layer . with this description
we certainly get more than we would need in a nonrelativistic theory , but we can recover at any time the nonrelativistic limit of all quantities by expanding in powers of @xmath59 .    it is important to stress that we have to expect , in general , different answers in the computation of a given observable taking ( [ hcoul ] ) or ( [ hq ] ) as starting point .
the dimensionality of the system is low enough so that the massless condition of the spinor reflects in the appearance of infrared divergences in perturbation theory .
thus , although we should not expect any discrepancy in the computation of local quantities ( like cutoff dependent quantities ) in the nonrelativistic model ( [ hcoul ] ) and in the model given by ( [ hq ] ) in the limit @xmath60 , the response of nonlocal quantities to the interaction turns out to be quite different in the infrared regime . from a technical point of view
, perturbation theory looses predictive power because of the mentioned infrared instabilities , and one has to resort to more sophisticated methods in order to obtain information from the quantum theory .
let us illustrate the above fact with the computation of the quantum corrections to the density of states @xmath61 near the fermi level ( @xmath62 ) .
this observable is given by @xmath63      \label{density}\ ] ] in the theory with the instantaneous coulomb interaction , the density of states is not modified to first order in perturbation theory
. the instantaneous interaction can be obtained from ( [ em ] ) , for instance , by considering formally its expression for @xmath64 . in this limit
we have , computing always at points @xmath65 in the two - dimensional layer @xmath66 then , it becomes obvious that this kind of interaction can not modify the @xmath67 dependence of the propagator @xmath68 , to the first perturbative order ( see also the discussion at the end of section 4 ) .
however , as mentioned above , the propagator ( [ propa ] ) can not be safely used in the computation of nonlocal quantities .
this remark is appropriate to the case of the density of states , whose determination in the quantum field theory given by ( [ hq ] ) goes as follows .
the inverse propagator may be decomposed , as usual , into a free part and the electron self - energy @xmath69 @xmath70 @xmath69 has a perturbative expansion in terms of one - particle - irreducible diagrams@xcite , which are in general ultraviolet divergent .
since our model can be treated as a genuine quantum field theory we expect these divergences to be local and susceptible of being absorbed into @xmath71 . for the moment we are only interested in finite corrections and will leave the issue of renormalization to the next section .
to first order in perturbation theory , the relevant diagram is shown in figure 3 and gives ( in the limit @xmath60 ) the finite contribution to @xmath72 @xmath73 when inserted in equation ( [ density ] ) the above expression gives the first quantum corrections to the density of levels @xmath74 it is already appreciated that the perturbative expansion of @xmath61 is not well - defined , since the last integral in ( [ 21 ] ) is clearly divergent . even if we could manage to regularize in some way this infrared divergence , we still would be left with a more serious problem . at the point
near which we want to measure the density of levels , i.e. near @xmath62 , the first perturbative correction dominates over the zeroth order term , and one can presume that higher order terms can produce increasingly large powers of @xmath75 .
perturbation theory is clearly not reliable in the infrared regime , but the kind of infrared instability described points at new physics near the fermi level .
henceforth we will pay attention to the local quantum field theory ( [ hq ] ) , which offers the possibility of analyzing the infrared behaviour by application of renormalization group methods .
we study first in this section the renormalization properties of the model . though the many - body theory of the 2d layer has a natural ultraviolet cutoff , we insist in absorbing ultraviolet divergent contributions into bare parameters of the theory since
this is a way of extracting relevant information about the scaling properties in the infrared@xcite . as we have mentioned before , the renormalizability of a quantum many - body theory is , in general , a nontrivial issue , since the couplings in the effective theory may depend on the point chosen on the fermi surface . in our case things become much simpler as our low - energy effective theory is a genuine quantum field theory . regarding the usual considerations for renormalizability in field theory ,
the only condition that we are missing is relativistic invariance since the fermi velocity of the electrons @xmath76 does not match the speed of light @xmath77 that appears in the dispersion relation of photons .
for this reason it is still necessary to check that the ultraviolet divergences of the theory can be absorbed into a redefinition of the scale of the fields and of the parameters in the hamiltonian .
we present here the analysis at the one - loop level , and conclude with a compelling argument showing that the model is renormalizable to all orders in perturbation theory .
the conditions which make the theory renormalizable can be stated more clearly in the lagrangian formalism . corresponding to the hamiltonian ( [ hq ] ) , we have the action @xmath78 in the quantum theory the coefficients of the terms in the action have to be adjusted so as to cancel the ultraviolet divergences .
we should start therefore with the bare action @xmath79 the renormalization coefficient @xmath80 may be determined from the electron self - energy @xmath69 , as well as the renormalization coefficient for the fermi velocity @xmath81 .
the coefficient @xmath82 is determined from the renormalization of the time component of the interaction vertex .
it turns out therefore that the spatial components have necessarily to be renormalized by the product @xmath83 .
this is a nontrivial check of the renormalizability of the theory , which arises in the absence of relativistic invariance .
provided that it is fulfilled , the renormalized theory may be made finite in terms of renormalized quantities @xmath84 @xmath85    in the previous section we already quoted the one - loop renormalized contribution to the electron self - energy @xmath86 .
we deal here with the ultraviolet divergent part of this object , regularizing it by working in analytical continuation to dimension @xmath87@xcite . from the diagram in figure 3
we have @xmath88 in order to ensure the finiteness of the electron propagator we take the renormalization coefficients @xmath89 these coefficients have a more complicated structure than those of a relativistic theory , as long as @xmath48 does not match the speed of light which , in our units , means that @xmath90 .
@xmath81 has , for instance , an infinite power series expansion in @xmath48 @xmath91 we will later use this expression in the consideration of the nonrelativistic approximation @xmath92 .
now we come to the renormalization of the interaction term in the action ( [ action ] ) .
the first quantum corrections to the vertex are given by the diagram in figure 4 , and we denote them symbolically by @xmath93 . the divergent contribution to the time component of
the vertex is actually different than that to the spatial components .
we find that @xmath94 we may think of @xmath93 as a correction to the bare vertex in ( [ action ] ) .
we have then @xmath95 by taking the quotient between ( [ dos ] ) and ( [ uno ] ) , we get an independent determination of the renormalization coefficient @xmath81 @xmath96 it is not obvious that this expression coincides with the alternative form ( [ alt ] ) , though it has to for the theory to be renormalizable .
in fact , after using the recursion relations which hold for hypergeometric functions@xcite one can show that the two functions in ( [ one ] ) and ( [ alt ] ) are identical .
we then have the nice situation of dealing with a nonrelativistic quantum field theory in which the fermi velocity of the electrons is renormalized exactly in the same way in the different terms of the action .
this precise coincidence is not casual .
one is led to suspect that it has to be consequence of some symmetry operating in the theory , which turns out to be nothing but gauge invariance .
actually , this should automatically lead in the quantum theory to the same renormalization of the kinetic part and the interaction term in the action , @xmath97 .
it is not difficult to check that this is satisfied in our model , to the perturbative order in which we are working . from expressions ( [ zpsi ] ) and ( [ uno ] ) we have @xmath98 which become identical after use again of the relations among hypergeometric functions .    the gauge symmetry @xmath99 is manifest at the classical level .
provided that it holds for the bare action ( [ action ] ) , it ensures that @xmath97 .
it is easily seen , by taking a gauge parameter @xmath100 independent of @xmath26 , that it also enforces the same renormalization of the fermi velocity @xmath48 in the different terms of the action .
as long as one is able to preserve the gauge invariance in the quantum theory , this turns out to guarantee the renormalizability to all orders of the perturbative expansion .
the rest of the renormalization coefficients in ( [ ren ] ) can be obtained from the following remark . the fact that @xmath101 has to equal @xmath82 to all orders in perturbation theory implies , as in standard quantum electrodynamics , that @xmath102 .
the renormalization coefficient of the electromagnetic field @xmath103 is determined from the polarization tensor , whose first perturbative contribution shown diagrammatically in figure 5 is a finite quantity @xmath104 we conclude , therefore , that @xmath105 remarking the absence of renormalization of the electric charge to first order in perturbation theory .
in this section we seek the connection of our model to the real physics of electrons in a two - dimensional layer . a single graphite sheet or any fullerene aggregate of sufficiently large size is the appropriate physical instance to look at .
the first question concerns the magnitude of the parameters pertinent to the real layer .
a standard value for the hopping parameter in graphite is @xmath106 . for the nearest neighbor separation in the lattice we can take @xmath107 . from these values and restoring for a moment
the dependence on @xmath77 and @xmath108 , the fermi velocity in our model turns out to be @xmath109 .
this is certainly much smaller than the speed of light , for which to extract relevant quantities to the electron physics in the layer we will have to take the limit @xmath60 and keep the leading order in previous formulas .    for a clearer description of our approach to renormalization group we will quote the renormalization coefficients in this section as functions of the ultraviolet momentum cutoff @xmath110 rather than in terms of the dimensional regularization poles
@xmath111 . at the one - loop level
this just amounts to replace @xmath111 by @xmath112 in all formulas .
the renormalization coefficient of the fermi velocity @xmath113 becomes , for instance , in the nonrelativistic limit @xmath114 it is not difficult to see that at each level of the perturbative expansion the leading order in the nonrelativistic approximation enters with a power @xmath115 .
we define , therefore , a consistent perturbative expansion in the nonrelativistic limit by taking @xmath92 ( or , equivalently , @xmath64 ) and @xmath116 this latter quantity becomes then the effective strength of the electronic interaction . in this limit
we obtain from ( [ zkin ] ) @xmath117 so that we do not find wavefunction renormalization at the one - loop level in the nonrelativistic limit .
we recall also the former result ( [ charge ] ) @xmath118    we have a nonrelativistic quantum theory with no renormalization of the electric charge , but in which the effective strength of the interaction @xmath119 is renormalized .
we rely on the original interpretation of renormalization adopted by wilson@xcite to determine the scale dependence of the effective interaction .
since we have shown that the model is renormalizable , we may impose the independence of the renormalized theory on the momentum cutoff @xmath110 . by means of renormalization group transformations we are able therefore to relate bare theories at different scales ( different cutoff ) which correspond to the same renormalized theory@xcite .
the essential information is encoded in the renormalization group equation , which for the electron propagator , for instance , states that under a change in the scale of the cutoff @xmath110 @xmath120 where @xmath121 the well - known solution to this equation takes the form @xmath122 where the effective parameters are given by @xmath123    increasing the scale @xmath124 of the cutoff is equivalent to measure the observables of the theory at large distance scale . from the one - loop order results of the previous section we see that @xmath125 is constant at this level , while @xmath126 is not .
we expect , therefore , the effective coupling of the electronic interaction @xmath119 to have a nontrivial renormalization group flow in the infrared regime . actually , from equations ( [ zetav ] ) and ( [ betav ] ) @xmath127 so that the asymptotic behaviour of the coupling is @xmath128    we have to bear in mind that , when integrating equation ( [ flow ] ) , we are relying heavily on perturbation theory .
if we recall our original estimate of the fermi velocity @xmath48 , we see that the coupling @xmath119 of the nonrelativistic theory is not a small parameter but , in fact , larger than one . nevertheless , we know that when dealing with a renormalized quantum theory the statement that a coupling like @xmath119 is small or large applies to a particular energy scale
. we then have to be very precise in establishing the regime in which the integration of the renormalization group equation becomes consistent .
let us remark that the solution ( [ sol ] ) of the renormalization group equation is equivalent to performing the sum of leading logarithm terms of the perturbative expansion and curing , therefore , the problem of the infrared instability of the density of states addressed at the end of section 2 .
taking the limit @xmath129 in the bare theory corresponds to measure at large distances ( or small energies with respect to the fermi level ) and in this regime perturbation theory becomes more and more accurate .
this is another instance of a physical system which has strong coupling at certain energies and an asymptotic regime at which perturbation theory becomes reliable , similar to what is believed to be true ( as supported by experiment ) in the quantum field theory of the strong interactions .
we are providing the example of a nonrelativistic theory in which this also happens , exchanging the ultraviolet regime in the case of the strong interactions by the infrared regime here .
as for the practical application of the flow ( [ sol ] ) , we remark that the nonrelativistic theory appears to be formally incomplete , since it does not give any hint about what the value of @xmath130 has to be .
this is related to the unlimited growth of the fermi velocity , and to the fact that , in the nonrelativistic approximation , there is no signal of a renormalization group fixed point .
later we will see that in the complete quantum field theory the fixed point is given by @xmath131 , but the flow in the graphite layer may not be necessarily controlled by this value of the renormalized coupling .
we believe , anyhow , that the observation of the phenomenon of scaling in the asymptotic region should be possible in a single graphite sheet .
regarding the possible existence of a natural infrared cutoff , it can be shown that in our model the coulomb interaction is not screened , in the usual metallic manner , by quantum effects .
this is due to the vanishing of the density of states at the fermi level .
thus , a finite screening length is only generated by giving to the system an excess of charge , as a consequence of doping . in order to study this effect
, we may place the system away from half - filling by introducing a chemical potential @xmath132 for the conserved charge @xmath133 the correct expression of the free electron propagator at finite @xmath132 is@xcite @xmath134 obviously , the @xmath135-dependent part can not change the ultraviolet behaviour of the loop - integrals and the renormalization coefficients of the quantum field theory are not modified by the nonvanishing chemical potential .
this may lead , however , to finite corrections .
the fermion loop diagram is now nonvanishing , opposite to what happened before at @xmath136 .
the computation of ( [ ncharge ] ) gives @xmath137 which , after multiplying by 2 ( the two independent fermi points ) gives the correct number of particles .
regarding the polarization tensor , a nonvanishing value of @xmath138 signals the appearance of a photon effective mass . at small
@xmath132 one finds actually the estimate for the inverse of the screening length @xmath139 which gives a measure of the screening effects away from half - filling .
we now turn , finally , to the issue of whether the electron observables may get anomalous dimensions in the nonrelativistic limit .
we remark that the absence of one - loop wavefunction renormalization in the limit @xmath92 is just an accident .
the signal of wavefunction renormalization is given by an ultraviolet divergent self - energy contribution of the form @xmath140 to establish the degree of divergence of the local contributions one can safely implement the nonrelativistic limit in the photon propagator , which leads to the expression ( [ propa ] ) .
the @xmath67-dependence in the vertex diagram in figure 3 can be absorbed then by a simple change of variables in the loop integral , so that the divergent part of the diagram can not bear any dependence on it .
at the two - loop level , the same independence of the local contributions on @xmath67 holds except for the vertex diagrams in figure 6 . in these diagrams
the loop integrals overlap in such a way that the previous argument does not apply . by using the nonrelativistic approximation ( [ propa ] ) and after a very lengthy calculation , one obtains from the sum of the two - loop diagrams @xmath141 where @xmath142 is a nonvanishing quantity , @xmath143 .
we see therefore that the electron observables acquire anomalous dimensions in the infrared regime .
according to equation ( [ rge ] ) , the electron propagator transforms under a change of scale @xmath124 by the factor @xmath144    we have therefore an electronic system which deviates slightly from fermi liquid behaviour .
in particular , we have seen that the vanishing of the density of states at the fermi level suppresses the usual screening of the coulomb interaction .
this property is maintained in the quantum theory , at least at the one - loop level , since the fermi level is not shifted by quantum corrections . also a gap does not open at the fermi points to the one - loop order ( as it would be if an effective mass for the spinor had been generated in perturbation theory ) .
the system is not an insulator , though it does not exhibit usual metallic properties either .
in the previous section we have mentioned that the nonrelativistic theory is formally incomplete since it does not know of any bound in the effective fermi velocity @xmath126 of electrons .
as this parameter grows steadily at large distances the nonrelativistic approximation becomes less accurate , up to a point in which the power series expansion in @xmath59 looses its meaning . in this relativistic regime
, we have to use the full expression of the beta function @xmath145 to analyze the behaviour of @xmath126 .
in particular , we are interested in knowing if a possible fixed point of the renormalization group flow exists , characterized by the condition @xmath146 by quoting ( [ alt ] ) in terms of elementary functions , this equation reads , at the one - loop level , @xmath147 for positive @xmath48 , the only solution to ( [ fix ] ) is @xmath148 , ) as @xmath149 . ] which in our units means @xmath48 equal to @xmath77 .
thus , the full quantum field theory has a built - in upper bound of the renormalization group flow , and the speed of light turns out to be the fixed point for the fermi velocity of electrons .    in fact
, it is not difficult to check that the one - loop beta function has a simple zero at @xmath150 , and to compute the flow in the neighborhood of the fixed point .
the renormalization group equation for @xmath126 becomes @xmath151 the asymptotic behaviour of the solution is , taking for @xmath130 the fixed point value , @xmath152 in figure 7 we show the plot of this scaling behaviour , which is practically indistinguishable in the range considered from that obtained with the nonrelativistic approximation .    recalling the vanishing of the electric charge beta function @xmath153 at the one - loop level
, we may assert that to this perturbative order we have identified an infrared fixed point given by @xmath131 and arbitrary ( though small ) value of @xmath125 .
the one - loop result is not satisfactory , however , since the independence of @xmath125 on the flow may be destabilized by higher - loop effects . in the renormalization group approach to fermi liquid theory@xcite , the consideration of the one - loop graph in figure 5 is enough since it encodes , via the renormalization group , information of the series of ladder diagrams . in the usual analysis of the many - body theory with an extended fermi surface
these are which dominate in the infrared regime .
because of the existence of fermi points in the honeycomb lattice at half - filling , though , there is no argument here to conclude the relevance of one set of diagrams against the rest . in our quantum field theory description
it seems that no diagram can be discarded a priori .
it is also for this reason that there are diagrams which are irrelevant in fermi liquid theory but give rise in our model to wavefunction renormalization and anomalous dimensions in the infrared regime .    regarding the renormalization of the electric charge
, we have analyzed the polarization tensor @xmath154 to second order in perturbation theory .
as stated above , the renormalization of the photon field @xmath56 is directly related by gauge invariance to the renormalization of the charge @xmath155 .
computing in the framework of dimensional regularization we have not found any divergent contribution of order @xmath156 to @xmath154 .
the reason is the same operating in the one - loop polarization tensor .
a nave dimensional analysis of the diagram shows that it should diverge linearly in momentum space . in the dimensional regularization scheme
this leads to the appearance of gamma functions at half - integer values of the argument , instead of poles at @xmath157 .
we believe that this scheme is the best suited for the regularization of the gauge theory it certainly is not going to break the gauge invariance in any event .
if , however , one insists on understanding the regularization in terms of a cutoff dependence , it can be checked that this does not lead either to a renormalization of the scale of the field @xmath50 .
one would have to look for divergent contributions of @xmath158 type accompanying the gauge invariant expression ( [ pol ] ) for @xmath154 .
these are not present at the two - loop order and one may conjecture , following the above dimensional argument , that they do not appear to any order in perturbation theory .    at or near the fixed point , the two - dimensional electron liquid we are describing is definitely not a fermi liquid .
the main reason for this assertion is the anomalous dimension acquired by the electron field @xmath51 , which may be obtained from ( [ zpsi ] ) evaluated now at @xmath150 .
at the fixed point the behaviour of the green function under a change of scale is @xmath159 or , equivalently , after introducing the precise value of @xmath160 @xmath161 we arrive in this way to solve the problem of the proliferation of @xmath162 terms in the perturbative computation of the density of states .
as is well - known , the renormalization group approach succeeds in summing up all these increasingly divergent contributions , in the form of an anomalous dependence on the electron frequency .
we have found that this sum is consistently achieved in the relativistic regime , where the fixed point lies and the electron field has a simple , anomalous scaling behaviour .
we have studied a system of electrons in the two - dimensional honeycomb lattice interacting via a coulomb potential .
the fact that the fermi surface of this lattice at half - filling consists of two isolated points renders the problem similar to the @xmath3 case and allows the formulation of the low energy effective theory as a standard quantum field theory .
we have been able to investigate the infrared behavior of the system by using renormalization group techniques and have seen that , indeed , it resembles more that of one - dimensional electrons than that of a two - dimensional fermi liquid .    before entering into the analysis of the phenomenological consequences of our work , there are a few points that we would like to emphasize .
let us first comment on the use made through the paper of the electromagnetic field propagator .
this use is unconventional both from the condensed matter point of view where the photon is mostly absent , and from the quantum field theory since , as we will strength now , we are not dealing with standard massless qed .
expression ( [ em ] ) shows the usual three - dimensional photon of qed .
this is what we use in our relativistic computations .
the constraint that , in the electronic interactions , the two points lay on the lattice plane , allows us to integrate out the perpendicular dimension .
the resulting expression read as a propagator on the plane certainly differs from the standard propagator in qed@xmath0 .
in particular its infrared behavior is quite different , it goes like @xmath163 instead of the standard @xmath164 .
this is the origin of , for instance , the finiteness of the photon self - energy up to the two - loop level and , hence , of the nonrenormalization of the electric charge .
let us now comment on the comparison between the instantaneous versus the retarded propagator .
the effect of retardation on the infrared behavior of nonlocal quantities has been exemplified with the computation of the quantum corrections to the density of states at the end of section 2 . in the theory given by the coulomb interaction ( [ hcoul ] ) , no one - loop corrections were found to the free value of the density of states . in the theory with the retarded interaction
our rg analysis shows that , despite the infrared divergence found at the one - loop level , the approach is consistent in the sense that the perturbative expansion of the nonrelativistic approximation does not have to make sense in the infrared when the fermi velocity is still far from the fixed - point value .
hence , we must resort to the full relativistic theory when computing the infrared behavior of nonlocal observables .
there , the rg allows to interpret the one - loop infrared instability shown in section 2 as the manifestation , in perturbation theory , of the anomalous dimension .
we have then found that our original system flows in the infrared towards a fixed point of the renormalization group that is a relativistic quantum field theory with anomalous dimension of the fermionic propagator .
our analysis has important consequences in condensed matter physics .
one of the more interesting is the anomalous screening of the coulomb interactions of the model . in a fermi liquid , formation of virtual electron - hole pairs ( i. e. , quantum corrections to the photon self - energy ) generates a new scale , the screening length .
this is the basis of the absence of long - range interactions in the ultralocal approach to interacting electrons ( hubbard model ) .
what we have seen is that , due to the vanishing of the density of states at the fermi level and to the absence of renormalization of the photon self - energy , there is no metallic screening in our materials @xcite .
the robustness of this result stands on the non - renormalization of the chemical potential ( the fermi level remains a set of two points after the corrections due to the interactions ) .
moreover we have seen that , although the low energy density of states goes to zero , no finite gap appears in the spectrum of charge excitations .
our results rule out the validity of descriptions in terms of on - site interactions only .
screening effects are not strong enough to reduce the couplings to simple on - site terms .
the one - particle properties are very different from those of a conventional fermi system .
the wavefunction renormalization leads to the disappearance of the quasiparticle peak at low energies .
thus , no effective description in terms of independent , and coherent , quasiparticles is possible .
the situation resembles the 1d luttinger liquid , where any amount of electron - electron interaction suffices to wipe out the quasiparticle pole , leading to incoherent one - particle excitations .
finally , the existence of nontrivial scaling laws implies that anomalous exponents should show up in a variety of properties , like the specific heat , susceptibility , etc .
the possibility of having an experimental test of these results remains open .    to conclude this summary and to point out possible future work
, we would like to remark the interest of investigating the connections of our work with the infrared catastrophe mentioned in the introduction as well as with other unconventional fermi liquid behavior discussed in the literature @xcite . on this respect , this work can be seen as a rigorous realization of the ideas settled in @xcite .
* acknowledgements *    this work has been partially supported by the cicyt .    99    a. a. abrikosov , l. p. gorkov and i. e. dzyaloshinski , _ methods of quantum field theory in statistical physics _
( dover , new york , 1975 ) .
e. m. lifshitz and l. p. pitaevskii , _ statistical physics _ , vol .
2 ( pergamon press , oxford , 1980 ) .
f. d. m. haldane , _ j. phys _ * c14 * ( 1981 ) 2585 .
p. w. anderson , _ phys .
lett . _ * 64 * ( 1990 ) 1839 ; * 65 * 2306 ; * 66 * 3226 .    for a recent review on the infrared catastrophe see the talk of p. w. anderson in _  the physics and the mathematical physics of the hubbard model " _ , nato a.r.w . , san sebastian , october 1993 .    for a theoretical physicist s vision of the problem
see j. polchinski , _  effective field theory and the fermi surface " _ , hep - th-9210046 .
g. benfatto and g. gallavotti , _ j. stat .
* 59 * ( 1990 ) 541 ; _ phys .
rev . _ * b42 * ( 1990 ) 9967 .
r. shankar , _ physica _ * a177 * ( 1991 ) 530 . for an extense review see r. shankar , _
 renormalization group approach to interacting fermions " _
_ , to appear .
wilson , _ phys .
rev . _ * b4 * ( 1971 ) 3174 , 3184 .
j. gonzlez , f. guinea , and m. a. h. vozmediano , _ nucl .
_ * b406 * ( 1993 ) 771 .
h. w. kroto , j. r. heath , s. c. obrien , r. f. curl and r. e. smalley , _ nature _ * 318 * ( 1985 ) 162 .
w. krtschmer , l. d. lamb , k. fostiropoulos and d. r. huffman , _ nature _ * 347 * ( 1990 ) 354 .
a. f. hebard , m. j. rosseinsky , r. c. haddon , d. w. murphy , s. h. glarum , t. t. m. palstra , a. p. ramirez , and a. r. kortan , _ nature _ * 350 * ( 1991 ) 600 .
m. j. rosseinsky _
_ , _ phys .
_ * 66 * ( 1991 ) 2830 .
k. holczer , o. klein , s .-
huang , r. b. kaner , k .- j .
fu , r. l. whetten and f. diederich , _ science _ * 252 * ( 1991 ) 1154 .    for a review on the properties of fullerenes ,
see k. holczer and r. l. whetten , _ carbon _ * 30 * ( 1992 ) 1261 .
j. gonzlez , f. guinea and m. a. h. vozmediano , _  electrostatic screening in fullerene molecules " _ , preprint cond - mat/9303007 ( to be published in _ mod . phys
s. sorella and e. tosatti , _ europhys .
* 19 * ( 1992 ) 699 , g. santoro , m. airoldi , s. sorella and e. tosatti , _ phys .
b _ * 47 * ( 1993 ) 16216 .
p. ramond , _ field theory . a modern primer_. ( benjamin / cummings , london , 1981 ) .
j. zinn - justin , lectures delivered at the cargse summer school 1973 ( unpublished ) .
k. symanzik , in _
particles , quantum fields and statistical mechanics _
( edited by m. alexanian and a. zepeda ) ( springer verlag , berlin , 1975 ) .
d. j. amit , _ field theory , renormalization group and critical phenomena _ , chap .
8 ( mcgraw - hill , new york , 1978 ) .
i. s. gradshteyn and i. m. ryzhik , _ table of integrals , series and products _ , p. 1044
( academic press , new york , 1980 ) .
k. g. wilson , _ phys .
* 28 * ( 1972 ) 548 .
s. a. chin , _ ann .
* 108 * ( 1977 ) 301 .
x. g. wen , _ phys .
_ * b42 * ( 1990 ) 6623 .    c. m. varma _ et al .
_ , _ phys .
* 63 * ( 1989 ) 1996 .
* figure 1 : * the planar honeycomb lattice .