the current understanding of topological matter @xcite is mostly based on a single - particle description .
a paradigmatic example is the integer quantum hall effect : in the absence of interactions , the hall conductivity has a deep topological interpretation@xcite , which explains its quantization and stability .
a more recent example is provided by the classification of time - reversal invariant insulators @xcite , which , again , relies on the properties of the noninteracting bloch functions .
understanding the effect of interactions on topological matter has become a very active area of research @xcite .
a natural model in which to explore such issues is the _ haldane - hubbard model_. the haldane - hubbard model describes spin-1/2 electrons on the honeycomb lattice , interacting via a local hubbard interaction of strength @xmath0 .
the electrons hop between nearest neighbor sites with hopping strength @xmath1 , and between next - to - nearest neighbor sites with alternating hopping parameters @xmath2 : the phases @xmath3 describe a transverse magnetic field , with zero net flux through the honeycomb plaquette .
finally , the system is also exposed to a staggered chemical potential , with strength @xmath4 on the two triangular sublattices . in the absence of interactions
@xcite this model shows , depending on the value of its parameters , a trivial insulating phase with vanishing transverse conductivity @xmath5 , or a quantum hall phase with @xmath6 .
these topological phases are separated by two critical curves in the @xmath7 plane , intersecting at the crossing points @xmath8 and @xmath9 . along the critical curves ,
the energy bands touch at a conical intersection ; at the crossing points , there are two such conical intersections , as in standard graphene . indeed if @xmath10 the system describes graphene with short range interaction .    from a theoretical viewpoint ,
the haldane topological phases have been argued to emerge in pure graphene sheets by spontaneous mass generation , due to the strong , unscreened coulomb repulsion @xcite . from an experimental viewpoint , the haldane model has been realized in ref.@xcite , and the topological phase transition has been observed .
the inclusion of a tunable hubbard interaction seems to be accessible by the present technology .
therefore , studying its effects on the transport coefficients is of fundamental importance for the next generation of cold atom experiments .
so far , the properties of the haldane - hubbard model have been investigated mostly via mean - field , variational and numerical analyses@xcite .    concerning the _ transverse conductivity _ , topological arguments for interacting systems@xcite
ensure that , away from the critical curves , @xmath11 can only take integer values , in units of @xmath12 ( here @xmath13 are the elements of the conductivity matrix , in the limit of zero frequency and zero temperature ) .
however , its specific value at a given point in the phase diagram can be different from the corresponding non - interacting value , in particular in the vicinity of the critical curves ( at weak coupling , far from the critical lines , the conductivity is known to be independent of the interaction @xcite ) .
the relevant question here is to distinguish between two scenarios : the first , in which small interactions are not able to generate new phases and their main effect is a shift of the critical curves , as found in certain 3d topological insulators @xcite ; and a second one , characterized by the emergence of a novel , interaction - induced , topological phase , like the one corresponding to @xmath14 , predicted for the haldane - hubbard model in ref.@xcite .
regarding the _ longitudinal conductivity _
, in the absence of interactions it is equal to @xmath15 , for all the values of @xmath16 _ on _ the critical lines , with the exception of the crossing points @xmath17 , where it is equal to @xmath18 ( of course , away from the critical lines @xmath19 ) .
there are no topological arguments ensuring that the _ critical _ longitudinal conductivity should remain quantized when the interaction is switched on : therefore , the relevant question here is whether the interaction introduces corrections breaking this exact quantization or not .
this question is related to a similar one discussed in the context of graphene , in which recent experiments @xcite showed that the optical longitudinal conductivity is essentially universal , and in excellent agreement with the value computed for the non - interacting model @xcite ; on the contrary , the interaction produces dramatic effects on other physical quantities , such as the fermi velocity @xcite . on the theoretical side ,
the universal behavior in graphene is still not completely understood @xcite , see , in particular , ref.@xcite for a recent review .
a rigorous result for short range interaction @xcite showed that , in order to get exact universality of graphene s longitudinal conductivity , one needs to fully take into account the non linear correction to the bands , even if such terms are irrelevant in the renormalization group ( rg ) sense .    in this paper
we compute the conductivity matrix of the haldane - hubbard model via exact rg methods , close to and at the critical lines , for weak interactions .
we take lattice effects into account , and we exploit _ lattice _ symmetries in order to reduce the number of independent running couplings , in a way similar to ref.@xcite for graphene .
the use of exact rg methods is motivated by the fact that the computation of conductivity is extremely sensitive to the choice of the regularization scheme@xcite . even though they are irrelevant in the rg sense , lattice and interaction effects produce , in general , _ finite corrections _ to the physical observables , and they must be taken into account in order to prove or disprove the emergence of new interaction - induced topological phases , as well as to address the issue of universality of the critical longitudinal conductivity .    by choosing the chemical potential @xmath20 so to fix the fermi energy half - way between the valence and conduction bands , the band gap can only close at the two fermi points @xmath21 , where @xmath22 is the _ valley _ index .
we prove that , close to criticality , the interacting euclidean two - point function is : & & s_2(k_0,p_f^+k ) = + & & = -    ik_0 z_1,r - m_r & v_r(-ik_1+k_2 ) + v_r(ik_1+k_2 ) & ik_0
z_2,r+m_r    ^-1-.1truecm(1+r(k_0,k))where the error term @xmath23 is subleading in the effective mass @xmath24 , in the matsubara frequency @xmath25 and in the quasi - momentum @xmath26 .
the parameters @xmath27 , depend non - trivially on the valley index @xmath28 and on the interaction .
in particular , the renormalized mass @xmath29 reads : [ fpm ] m_r,=w 3t_2-f_(u , w , ) where @xmath30 is the valley index , and @xmath31 is expressed in the form of a _ convergent _ renormalized series , whose first non - trivial order is given by eq .
the dressed critical lines , defined by the condition that the renormalized mass vanishes , are also modified by the interaction , see fig.[figaltezza ]
. similarly , the fermi velocity @xmath32 and the wave function renormalizations @xmath33 have non - trivial interaction corrections and , remarkably , @xmath34 and @xmath35 are _ different _ , as shown in fig.[fig1 - 2 ] .
all these non - universal renormalizations are _ absent _ in effective relativistic descriptions : by neglecting the ( irrelevant ) non - linear corrections to the energy bands , one would obtain a nambu - jona lasinio model , in which lorentz and chiral symmetry would imply the invariance of @xmath32 , the invariance of the critical lines and @xmath36 .
however , these extra symmetries are broken by the lattice , and the renormalization of the effective parameters is a physical signature of many body interaction that should be visible in real systems , e.g. , in cold atom experiments .
these non - universal parameters also enter the computation of the conductivity : remarkably , they are related by exact _
lattice ward identities _ , which induce non - trivial cancellations and imply subtle universality properties , as stated in the following theorem .
we recall that @xmath13 are the elements of the kubo conductivity matrix , in the limit of zero frequency and zero temperature .
we also denote by @xmath37 their values on the renormalized critical curves .    , @xmath38 , and different values of @xmath0 .
@xmath39 , where the values of @xmath40 are reported in the figure . for @xmath41 the system is a topological insulator ( ti ) , while for @xmath42 the system is a trivial , normal , insulator ( ni ) . ]    * theorem . *
_ there exists @xmath43 such that for @xmath44 , the system is massless if and only if the right side of eq .
this condition defines two renormalized critical curves intersecting at @xmath17 , separating two non - trivial topological phases , characterized by transverse conductivity @xmath45 , from two standard insulating phases , see fig.[figaltezza ] . on the renormalized critical curves ,
the critical longitudinal conductivity @xmath46 , @xmath47 , is quantized : if @xmath48 , _ ii^ = , while @xmath49 at @xmath17 . _
thus , the critical lines acquire non - universal , interaction - dependent corrections , but they still separate topological regions labelled by @xmath50 from the trivial ones , labelled by @xmath51 , see fig.[figaltezza ] .
new intermediate phases characterized by the quantum number @xmath52 are _ rigorously excluded _ at weak coupling : the universality class of the topological transition remains unchanged .
the effect of the repulsive interaction is to _ enlarge _ the topologically non - trivial region , see fig.[figaltezza ] .
this enhancement agrees with the numerical findings of ref.@xcite and is presumably a sign that repulsive interactions in graphene - like systems can favor the spontaneous generation of the topological insulating phase @xcite .
even if not protected by any topological argument , the critical longitudinal conductivity @xmath53 is _ exactly universal _ and equal to half the one of graphene , on the whole critical line , with the exception of the special crossing points @xmath8 and @xmath9 , at which the value of @xmath53 is the same obtained for interacting graphene@xcite , namely @xmath54 : each dirac cone contributes with a universal quantity @xmath55 to the critical longitudinal conductivity . of course , away from the critical curves , the longitudinal conductivity is exactly zero .
our results are in agreement with a low energy description in terms of an effective action that includes a non - trivial chern - simons term , whose coefficient ( the hall conductivity ) is proportional to the difference of the signs of the renormalized masses , @xmath56 , rather than the bare ones , as one would get in the relativisitic approximation @xcite .    the theory that we develop is _ non perturbative _ , in the sense that it allows us to express all the correlations and transport coefficients in terms of _ convergent _ series .
as it will appear from the analysis , our non - perturbative bounds on the correlation functions , once combined with ward identities , allow us to conclude the universality of the conductivity , without exploiting explicit cancellations at all orders .
we have not tried to optimize the estimate for the radius of the convergence domain , which , therefore , is expected to be far from the values of @xmath0 where interaction - induced phase transitions might take place .
however , we believe that the range of validity of our convergent expansions could be improved by combining our analysis with numerical techniques , as it is done , for instance , for the stability of kam tori in classical mechanics .
finally , we stress that our analysis only requires the interaction to be short - ranged , we considered the hubbard interaction just for the sake of definiteness .
the paper is organized as follows . in section [ sec2 ]
we define the haldane - hubbard model , and derive the exact lattice ward identities for its correlation functions . in section [ sec3 ]
we perform an exact renormalization group analysis for the correlations , we classify the allowed running coupling constants by the exact lattice symmetries of the system , and we compute the decay of the correlations at large distances , as well as the renormalized critical line . in section [ sec4 ]
we prove the quantization of the hall conductivity across the critical line , and the universality of the critical longitudinal conductivity .
the haldane - hubbard model describes interacting fermions on the honeycomb lattice @xmath57 , which can be understood as the superposition of two triangular sublattices @xmath58 and @xmath59 ; see fig .
the triangular sublattice @xmath58 is generated by the basis vectors _ 1 = ( 3 , -),_2 = ( 3 , ) . with each sublattice
, we introduce fermionic creation and annihilation operators @xmath60 , @xmath61 , where @xmath62 is the spin degree of freedom , @xmath63 .
the hamiltonian is : [ eq : h ] h = h_0 +
uv - n , where : @xmath64 is the noninteracting hamiltonian , @xmath65 is the hubbard interaction and @xmath66 fixes the chemical potential .
the noninteracting hamiltonian is@xcite : [ eq : h0 ] h_0 & = & -_=_x , y t_1 + & & - _ = _ x , y   + & & + w_= ; the first sum is over nearest - neighbours on @xmath57 , while the second is over next - to - nearest neighbours .
each site on @xmath58 is connected to its three nearest - neighbours on @xmath59 by the vectors : _ 1 = ( 1 , 0),_2 = ( -1 , ) , _ 3 = ( -1 , - ) .
the next - to - nearest neighbour hopping parameter @xmath67 is defined as : t_2(x , x + _ i ) = e^it_2,t_2(x , x - _ i ) = e^-it_2 , for @xmath68 .
explicitely ( see fig .
[ fig3 ] ) : _ 1 = _ 1 - _ 2,_2 =
_ 2,_3 = -_1 .
the hubbard interaction term is , as usual : [ eq : v ] v = _ x , where the sum ranges over the full honeycomb lattice ; the density operator @xmath69 is : n_x , = \ {    cc a^+_x , a^-_x , & + b^+_x , b^-_x , &    .
, in terms of which we also have @xmath70 . the factors @xmath71 in eq .
( [ eq : v ] ) amount to a redefinition of @xmath72 , and simplify the functional integral representation of the model ( see section [ sec : funcint ] ) .
we denote the finite volume version of @xmath73 by @xmath74 , with periodic boundary conditions .
the finite volume and finite temperature gibbs state is : _ ,
l=\{e^-h_l}/\{e^-h_l } , and we let _ = _
l_,l , = _ _ , l .
_ correlations , current and conductivity .
_ it is convenient to define ^+_x,= ( a^+_x , , b^+_x+_1,),^-_x , = ( ^+_x,)^ ; also , for any inverse temperature @xmath75 , we let @xmath76 be their evolution at ` imaginary time ' @xmath77 . for general @xmath78 ,
we extend @xmath79 anti - periodically ( of anti - period @xmath75 ) beyond the basic interval @xmath80 .
the fourier transform of the fields is defined as @xmath81 , where @xmath82 is the brillouin zone @xcite .
the 2-point function is s_2(,)&= & + & = & _ re^-i(-)s_2(),where @xmath83 , @xmath84 , @xmath85 is the fermionic time - ordering operator ( which orders imaginary times in decreasing order @xcite ) , and @xmath86 , where @xmath25 is the matsubara frequency .
note that @xmath87 is a @xmath88 matrix ( with indices in the ` sublattice ' space ) and its definition is independent of the choice of @xmath89 .
the current is defined via the peierls substitution ( see app.[app.ver ] ) , and is equal to j_p(x_0)=_= _ ^+_k+p,(x_0)m(k , p)^-_k,(x_0).[eq : curr]the two components @xmath90 , @xmath47 , of @xmath91 are the bare vertex functions , which are @xmath88 matrices , with elements labelled by the spinor indices . for the explicit expression of the bare vertex functions , see app.[app.ver ] .
the current - current and the vertex correlations are defined , respectively , as & & k_()=_rdx_0 e^-i p_0x_0 , + & & g_(,)=_rdx_0_rdy_0 e^-i p_0x_0+i(k_0+p_0)y_0 + & & 1.2truecm , [ eq : kmn ] where @xmath92 , @xmath93 with m_0(p ) =    1 & 0 + 0 & e^-ip_1    ; the labels @xmath94 refer to the components of the current defined in eq .
( [ eq : curr ] ) .
moreover , @xmath95 is the trace per unit volume , and the semi - colon indicates that the expectation is truncated . for later reference , we also introduce the vertex function : _ ( , ) = s_2 ^ -1(+)g_(,)s_2 ^ -1(),[eq : vf]where @xmath96 is the inverse of the 2-point function , thought of as a @xmath88 matrix .    finally , the d.c . kubo conductivity is defined in terms of the current - current correlation , in units such that @xmath97 , as : _ ij = - _ p_00^+ , [ eq : sijdef]where @xmath98 and @xmath99 is the area of the fundamental cell .    _
ward identities . _ the continuity equation for the lattice current eq .
, when averaged against an arbitrary number of field operators , implies exact identities among correlation functions ( ward identities ) , valid for any value of the interaction @xmath0 . in particular , the one relating the 2-point and the vertex functions , which will play an important role in the following , reads as follows : _
= 0 ^ 2 ( i)^_,0p_g_(,)=s_2(+)m_0(p)- m_0(p ) s_2().[1.14 ] if we derive this equation with respect to @xmath100 , compute the result at @xmath101 and recall the definition of the vertex function , we find : _ ( , * 0*)=(-i)^_,0_s_2 ^ -1()+.[eq : wi-2]in the following , @xmath102 will be denoted simply by @xmath103 .    _ the non - interacting case . _
if @xmath104 , the band structure and the phase diagram can be computed explicitly : the bloch hamiltonian is @xcite & & h_0(k ) = [ hbloch ] + & & =    -2t_2_1(k ) + m(k ) & -t_1^*(k ) + -t_1(k ) & -.3truecm -2t_2_1(k ) - m(k )    where @xmath105 and & & _
i=1 ^ 3(k_i),_2(k ) = _
i=1 ^ 3 ( k_i ) + & & m(k ) = w - 2t_2 _ 2(k ) .
the corresponding energy bands are @xmath106 to make sure that the energy bands do not overlap , we assume that @xmath107 .
the two bands can only touch at the _ fermi points _
@xmath108 , which are the two zeros of @xmath109 , around which @xmath110 .
the condition that the two bands touch at @xmath111 , with @xmath112 , is that @xmath113 , with m_ m(p_f^ ) = w + 3t_2.therefore , the unperturbed critical curves are given by the values of @xmath7 such that : w=33 t_2 , which correspond to the dotted curves in fig.[figaltezza ] .
fixing the chemical potential in such a way that the fermi energy lies in between the two bands , = -2t_2_1(p_f^)=-3 t_2 , the system passes from a semi - metallic behavior , when @xmath16 is on the critical line , to an insulating behavior , characterized by the exponential decay of correlations , when @xmath114 .
the insulating phase consists of four disconnected regions in the @xmath16 plane , two of which are ` topologically trivial ' , while the other two have non - zero hall conductivity , see fig.[figaltezza ] : more precisely , if @xmath114 , _ 12 = ,
= ( m_- ) - ( m_+ ) .
we now construct the interacting correlations and phase diagram , by using a convergent renormalized expansion , in the spirit of ref.@xcite . in this section
, we introduce the functional integral formulation of the model , discuss the exact lattice symmetries of the fermionic action , and describe the infrared integration , including the study of the flow of the running coupling constants .
one of the main results of this section is the equation for the interacting critical line .
we are interested in the semi - metallic and insulating regimes of the interacting system .
we , therefore , set the chemical potential accordingly ( its value will be different , in general , from the unperturbed one ) : = -2 t_2_1(p_f^)- , where @xmath115 ( the shift of the chemical potential ) must be chosen as a function of @xmath116 so that the renormalized propagator either has a linear , ` conical ' , infrared singularity ( along the interacting critical line ) , or is gapped ( in the insulating phase ) .    the generating function @xmath117 for correlations , in which @xmath118 is the external field conjugated to the fermionic fields , and @xmath119 is the external field conjugated to the current , can be written as the following grassmann integral : e^(f , a)=,[eq.11 ] where : @xmath120 , with @xmath121 and @xmath122 , is a two - component grassmann spinor ( it is the grassmann counterpart of @xmath79 ) , whose components will be denoted by @xmath123 , with @xmath124 ; @xmath125 is the fermionic gaussian integration with propagator g(,)=_r e^-i(- ) g(),[prop]where , letting @xmath126 , g()=    & -ik_0+r(k)+ m(k ) & -t_1 ^*(k ) + & -t_1 (k ) & -i k_0+r(k)- m(k )    ^-1and , at contact , @xmath127 should be interpreted as @xmath128 $ ] ; v()=_r dx_0_x_a _
= 1,2(un^_,n^_,+_= , n^_,),where @xmath129 ; and , finally , & & ( , f)=_rdx_0_x_a_=(^+_,f^-_,+f^+_,^- _ , ) , + & & ( j , a)=_r^3 a_,_,,where @xmath130 , in which @xmath131 are the _ bare _ vertex functions , namely : @xmath132 , and , if @xmath47 , @xmath133 are the two components of the ( matrix - valued ) vector @xmath91 defined in and following lines . in terms of these definitions , the correlations can be re - expressed as s_2 ( , ) & = & ( 0,0 ) , + k _ ( , ) & = & ( 0,0 ) , and of suitable linear combinations of g_2,1;(,,)=(0,0 ) .
we now compute the generating function eq . via a renormalized expansion , which is convergent uniformly close to ( and even on ) the critical line .
note that , on this line , the grassmann integral has an infrared problem . in order to resolve and re - sum the corresponding singularities ,
we proceed in a multi - scale fashion .
first of all , we distinguish the ultraviolet modes , corresponding to large values of the matsubara frequency , from the infrared ones , by introducing two compactly supported cut - off functions , @xmath134 , supported in the vicinity of the fermi points ( more precisely , we let @xmath135 , where @xmath136 is a smoothed out characteristic function of the ball of radius @xmath137 , with @xmath137 equal to , say , @xmath138 , and @xmath139 ) and by letting @xmath140 .
we correspondingly split the propagator in its ultraviolet and infrared components : g(,)=g^(1)(,)+_= e^-ip_f^(x - y)g_^(0)(,)[eq:4 ] where @xmath141 and @xmath142 are defined in a similar way as eq .
, with @xmath143 replaced by @xmath144 and by @xmath145 , respectively .
we then split the grassmann field as a sum of two independent fields , with propagators @xmath146 and @xmath147 : @xmath148 and we rewrite the grassmann gaussian integration as the product of two independent gaussians : @xmath149 . by construction , the integration of the ` ultraviolet ' field @xmath150 does not have any infrared singularity and , therefore , can be performed in a straightforward manner , thus allowing us to rewrite the generating function @xmath151 as the logarithm of p(d^(0))e^-v^(0)(^(0))+b^(0)(^(0 ) , f , a),[eq : v0 ] where @xmath152 and @xmath153 are , respectively , the effective potential and the effective source ( which depend explicitly on , respectively , @xmath154 and @xmath155 ) , @xmath156 is independent of @xmath154 ( and depends explicitly on @xmath157 ) , and @xmath158 . both @xmath152 and @xmath153
are expressed as series of monomials in the @xmath159 fields , whose kernels ( given by the sum of all possible feynman diagrams with fixed number and fixed location of the external legs ) are _ analytic functions _ of the interaction strength , for @xmath0 sufficiently small .
the proof of their analyticity is based on a determinant expansion and on a systematic use of the gram - hadamard bounds , see ref.@xcite .      before tackling the multi - scale integration of the infrared modes ,
we make a digression about the symmetry structure of the effective potential , and in particular of its local parts : the purpose is to classify the possible relevant and marginal coupling constants . in the case
@xmath160 ( standard graphene model ) the lattice symmetries severely constrain the form of the quadratic terms in the effective potential : in particular , the interaction does not shift the chemical potential , nor does it generate a mass@xcite . in the general case ( @xmath161 ) the model is invariant under the following symmetry transformations ( since they do not mix the spin indices , for notational convenience we temporarily drop the spin labels from the formulas ) .
we discuss the symmetries in the _ absence _ of external fields , since we will use them only to infer the structure of the relevant and marginal contributions to the effective potential @xmath152 .
once the structure of these terms is known , the structure of the marginal contributions to the effective source @xmath153 can be computed by using the ward identity .
\(1 ) _ discrete rotations : _ _ ^-e^ik(_3-_1)n_-_t^- , _ ^+_t^+e^-ik(_3-_1)n_-
where , denoting the pauli matrices by @xmath162 , we defined n_-=(1-_3)/2,t=(k_0,e^-i_2k ) ; that is , @xmath163 is the spatial rotation by @xmath164 in the counter - clockwise direction .
\(2 ) _ complex conjugation _ : ^_^_- , combined with cc^*,- , where @xmath165 is a generic constant appearing in @xmath125 or in @xmath166 .
\(3 ) _ horizontal reflections _ : ^-__1 ^ -_r_h,^+_^+_r_h_1 , with r_h=(k_0,-k_1,k_2),(w,)(-w,- )    \(4 ) _ vertical reflections _ : ^_^_r_v , with r_v=(k_0,k_1,-k_2),- .
\(5 ) _ particle - hole _
: ^-_i^+,t_p,^+_i^-,t_p , with p=(k_0,-k_1,-k_2),- .
note that , at fixed @xmath167 , the theory is invariant under the transformations ( 1 ) , ( 2)+(4 ) , and ( 2)+(5 ) . in particular
, these transformations leave the quadratic part @xmath168 of the effective potential @xmath169 invariant .
this means that : w_2()&=&e^-ik(_1-_2)n_- w_2(t^-1)e^ik(_1-_2)n_- + & = & w_2^*(-k_0,-k_1,k_2)[symm ] + & = & w_2^(-k_0,k_1,k_2).as we will see in the next section , the values of @xmath170 and of its derivatives at the fermi points define the _ effective coupling constants_. by , we find , for @xmath171 , -.2truecm w_2(_f^)&=&e^-i3n_-w_2(_f^)e^i3n_- + & = & w_2^*(_f^)= w_2^(_f^),which implies that w_2(_f^)=_+__3,[eq : w2_loc0]for two _ real _ constants @xmath172 and @xmath173 .
if we derive with respect to @xmath174 and compute the result at @xmath175 , we find : _ w_2(_f^)&=&e^-i3n_-t_w_2(_f^)e^i3n_- + & = & ( -r_v)_w_2^*(_f^)[symm_der ] + & = & ( -p)_w_2^(_f^),where @xmath176 ( resp .
@xmath177 ) is the diagonal matrix with diagonal elements @xmath178 ( resp .
@xmath179 ) . by using , it is straightforward to check that _w(_f^ ) =    -i z_1, k_0 & -u_(-i k_1
+ k_2 ) + -u_(i k_1 + k_2 ) & -i z_2, k_0    , [ eq:6]where @xmath180 , and @xmath181 are _ real _ constants . in conclusion , for general values of @xmath167 , the linearization of @xmath170 at @xmath175 is parametrized by 5 real constants , namely @xmath182 and @xmath183 , the first two are relevant coupling constants , and the other three are marginal .
note that , in general , the values of these constants depend on @xmath184 ( therefore , there are 5 of them at @xmath185 and 5 more at @xmath186 ) .
note also that , in general , @xmath187 , i.e. , the wave function renormalization depends explicitly on the spinor index , an effect that can be checked explicitly at second order in perturbation theory ( see below ) , and can not be explained purely in terms of the relativistic approximation of the model around the fermi points .
note that there are special points in the @xmath188 plane , for which the model has more symmetries , and where the number of independent couplings is smaller than in the general case .
for instance , if @xmath189 , the model is invariant under all the 5 symmetry transformations listed above , in which case it is straightforward to see that & & _ =_-,_=0,u_ = u_- , + & & z_1, = z_2, = z_1,- = z_2,- .
a similar discussion applies to the case @xmath190 , @xmath191 .    finally , if @xmath192 , the model is invariant under the following additional symmetry transformation ( see also ref.@xcite ) : ^-_,-i _ 1_3 ^ -_-r_v,,^+_,-i^+_-r_v,_3_1 , which implies that @xmath193 so that , in particular , _ =0,z_1, = z_2, .
a similar discussion applies to @xmath194 .
let us now describe the integration of the infrared fields .
we shall focus on the semi - metallic behavior of the system at , or very close to , a generic point of the critical line .
moreover , since we are interested in the behavior of the current - current correlations around @xmath101 , we shall assume that the external field @xmath195 is supported in the vicinity of the origin ( in particular , we assume that it vanishes in the vicinity of @xmath196 , @xmath171 ) .    by dimensional considerations , the quadratic terms in the effective action are _ relevant _ , and , the ones corresponding to the renormalization of the mass are of particular importance .
the flow of the effective mass tends to diverge linearly under the rg iterations , which signals that , in general , the location of the critical lines is changed by the interaction . in order to construct a convergent expansion , we need to dress the mass , after which we determine the location of the renormalized critical lines , which is given by the condition that the dressed mass vanishes .
more in detail , we proceed as follows .
we perform the integration of the infrared modes in iteratively , by decomposing the fermionic fields as @xmath197 as @xmath198 , where @xmath199 is a grassmann field whose propagator is supported on the momenta @xmath174 such that @xmath200 , and by integrating the fields @xmath201 step by step . after the integration of the modes on scales @xmath202 , we rewrite the generating function @xmath117 as the logarithm of p(d^(h))e^-v^(h)(^(h))+b^(h)(^(h ) , f , a),[eq : vh ] where @xmath203 and @xmath204 are , respectively , the effective potential and source terms , to be defined inductively in the following . moreover , @xmath205 is the grassmann gaussian integration with propagator ( diagonal with respect to the @xmath62 and @xmath184 indices ) g^(h)_(,)&=&p(d^(h))^-(h)_,,^+(h)_,, + & = & e^-i(- ) g_^(h)( ) , where @xmath206 and , letting @xmath207 , @xmath208 $ ] and @xmath209 ( here @xmath136 is the cutoff function defined a few lines before ) , g_^(h)()=_h( )    a_1, , h( ) & b^*_ , h( ) + b_ , h( ) & a_2, , h( )    ^-1,[eq : prop.2]with & & a_, , h()=-ik_0z_, , h+r_(k)+(-1)^-1 m_ , h(k ) , + & & b_ , h()=-v_ , h ( ik_1+k_2)+s_(k).[eq : prop.3 ] in which @xmath210 , @xmath211 and @xmath212 are , respectively , the wave function renormalizations , the effective mass and effective velocities , to be defined inductively in the following .
their initial values are : z_j , ,0=1,m_,0(k)=m(k+p_f^),v_,0=32t_1 .    in order to clarify the inductive definition of the effective potential , source , etc , we now describe the integration step at scale @xmath213 .
we start from , where @xmath214 is a sum of even monomials in the @xmath215 fields , whose kernels of order @xmath216 are denoted by @xmath217 ( for notational simplicity , we temporarily drop the space - time , spin , spinor and valley indices of the fermionic fields ) .
similarly , we denote the kernels of @xmath204 of order @xmath216 in @xmath215 , @xmath218 in @xmath118 and @xmath219 in @xmath119 , by @xmath220 .
the scaling dimension of the kernels @xmath221 and @xmath220 is ( see ref.@xcite ) d=3-n - m - q,[eq : scdim]with the convention that @xmath222 corresponds to relevant , @xmath223 to marginal , and @xmath224 to irrelevant operators .
note that the only relevant terms are those with @xmath225 , and the only marginal terms are those with @xmath225 and @xmath226 ( note that , by construction , @xmath227 is positive and even ) . in particular , the effective electron - electron interaction , corresponding to the case @xmath228 and @xmath229 , is _
irrelevant_.    in order to define a convergent renormalized expansion , we need to re - sum the relevant and marginal terms . for this purpose
, we split @xmath203 and @xmath204 into their _
local _ and _ irrelevant _ parts ( here , for simplicity , we spell out the definitions only in the @xmath230 case , the general case is treatable analogously , along the lines of , e.g. , sect .
12 of ref.@xcite , or ref.@xcite ) : @xmath231 and @xmath232 , where , denoting the quadratic part of @xmath203 by @xmath233 and the part of @xmath204 of order @xmath234 in @xmath235 by @xmath236 we let : & & v^(h)()=_=_= + & & ^+_,, [ w^(h)_2;(*0*)+_w^(h)_2;(*0*)]^-_,, , and & & b^(h)(,0,a)=_=_=_=0 ^ 2 + & & a_,^+_+,, w^(h)_2,1;,(*0*,*0*)^-_,, . by the symmetries
discussed in the previous section ( see , in particular , and ) & & v^(h)()=_=,where @xmath237 are real constants .
moreover , by using the ward identity , we find that & & b^(h)(,0,a)=_=_=_=0 ^ 2 + & & a_,^+_+,, _ ,  , h^-_,,,[eq : lbh]where _ 0, , h&=&-_=1 ^ 2(z_, , h+z_, , h)n _ + _ 1, , h&=&-(v_ , h+u_ , h)_2 + _ 2, , h&=&-(v_ , h+u_ , h)_1 in which @xmath238 and @xmath239 are the standard pauli matrices .
once the effective potential and source have been split into local and irrelevant parts , we combine the part of @xmath240 in the second line of with the gaussian integration @xmath205 , thus defining a dressed measure @xmath241 whose propagator @xmath242 is analogous to @xmath243 , with the only difference that the functions @xmath244 , @xmath245 in - are replaced by a_, , h-1()&=&-ik_0z_, , h-1()+r_(k ) + & & + ( -1)^-1 m_ , h-1( ) , + b_ , h-1()&=&-v_ , h-1( ) ( ik_1+k_2)+s_(k ) , with & & z_, , h-1()=z_, , h+z_, , h_h( ) , + & & m_ , h-1()=m_ , h(k)+_ , h_h( ) , + & & v_ , h-1()=v_ , h+u_ , h_h().now , by rewriting the support function @xmath246 in the definition of @xmath242 as @xmath247 , we correspondingly rewrite : @xmath248 , where @xmath249 is defined exactly as in - , with @xmath213 replaced by @xmath250 , and @xmath251 defined by the flow equations : & & z_, , h-1=z_, , h+z_, , h , + & & m_ , h-1(k)=m_ , h(k)+_ , h,[eq : flow.1 ] + & & v_ , h-1=v_ , h+u_ , h.at this point , we integrate the fields on scale @xmath213 , and define : & & e^-v^(h-1)()+b^(h-1)(,f , a)+w^(h)(f , a)=c_hp(d^(h ) ) + & & e^-f_^(h)(^(h)+)+v^(h)(^(h)+)+b^(h)(^(h)+ , f , a),where @xmath252 is the gaussian integration with propagator @xmath253 , @xmath254 , and @xmath255 . finally ,
letting @xmath256 , we obtain the same expression as , with @xmath213 replaced by @xmath250 .
this concludes the proof of the inductive step , corresponding to the integration of the fields on scale @xmath213 .
the integration procedure goes on like this , as long as the two effective masses @xmath257 are small , as compared to @xmath258 .
if we are not exactly at the ` graphene point ' @xmath189 , i.e. , if we are close to , or at , any other point on the critical line but the origin , then after a while we reach a scale @xmath259 at which @xmath260 ( possibly , @xmath261 , in the case that @xmath262 is of order 1 , i.e. , if @xmath167 are far enough from the graphene point ) .
note that , once we reach scale @xmath259 , the field @xmath263 is massive ` on the right scale ' @xmath264 . at that point , we integrate out the field @xmath263 in a single step , and we are left with a ( chiral ) theory , whose only dynamical degree of freedom is @xmath265 , with @xmath266 .    from that scale on ,
we integrate @xmath267 in a multi - scale fashion , analogous to the one discussed above , with the important difference that only the running coupling constants corresponding to the valley index @xmath268 continue to flow .
the multi - scale integration goes on until we reach a scale @xmath269 such that @xmath270 , at which point we can integrate out the remaining degrees of freedom in a single step .
the criticality condition , i.e. , the condition that the system is on the ( renormalized ) critical line , corresponds to the condition that @xmath271 .
the multi - scale integration described in the previous section defines a flow for the effective chemical potential @xmath272 , the effective mass @xmath273 , the effective wave function renormalization @xmath274 , and the effective fermi velocity @xmath212 .
the flow of @xmath275 and @xmath212 is driven by eqs.([eq : flow.1 ] ) , while @xmath276 where @xmath277 is the ( @xmath115-component of the ) beta function , which is defined in terms of the sum of all the local quadratic contributions in renormalized perturbation theory , and should be thought of as a function of @xmath0 and of the sequence of the effective coupling constants .
remember that the flow drives the effective couplings with both @xmath278 and @xmath279 , up to the scale @xmath259 ; then the flow of the couplings with @xmath280 is stopped , and only the couplings with @xmath268 continue to flow until scale @xmath269 ( possibly @xmath271 )
.    the multi - scale procedure is well defined , and the effective potentials are , step by step , given by _ convergent _ expansions , provided : ( i ) @xmath0 is small enough , ( ii ) @xmath281 remain small for all scales , and ( iii ) @xmath282 remain close to their initial ( bare ) values , for all scales .
note that , in order for condition ( ii ) to be valid , we need to properly fix the initial condition on the chemical potential , as discussed in the following .
in addition , note that , once that the flows of @xmath274 and @xmath212 are controlled , then the marginal contributions to the effective source term @xmath283 are automatically under control , thanks to and following lines .
the key fact , which allows us to control the flow of the effective couplings , is that , since the electron - electron interaction is _ irrelevant _ , with scaling dimension @xmath284 ( cf . with ) , then the scaling dimensions of all diagrams with at least one interaction vertex can be effectively improved by one , see ref.@xcite .
in particular , @xmath285 , for any @xmath286 and a suitable constant @xmath287 , and similarly for the beta functions of @xmath274 and @xmath212 . [ the reason why we lose , in general , an @xmath288 in the decay exponent as @xmath289 , is that we need to use a little bit of decay @xmath290 in order to sum over all diagrams and scales , see ref.@xcite for details . ]    in order to guarantee that the flow of the chemical potential remains bounded , we fix the initial data ( via a fixed point theorem , such as the contraction mapping theorem ) so that @xmath291 , in the limit as @xmath292 . thanks to the dimensional gain of @xmath293 , due to the irrelevance of the interaction
, we actually find that @xmath294 tends to zero , as @xmath289 , exponentially fast : @xmath295(const.)@xmath296 .
once we imposed that @xmath294 remains bounded for all scales @xmath297 , we can a posteriori check that @xmath298 is also bounded for all scales @xmath299 : in fact , the beta function @xmath300 , for @xmath301 , can be rewritten as @xmath302 $ ] , where the difference in square brackets can be straightforwardly shown to be proportional to @xmath303 [ if all the masses @xmath304 were zero , then the model would be symmetric under the exchange of @xmath184 in @xmath305 , as in ref.@xcite , see also section [ sec.symm ] above ; therefore , the difference @xmath306 between the contributions with different valley indices must be proportional to a mass term @xmath307 , which is smaller than ( const.)@xmath308 .
therefore , the flow of @xmath298 , for @xmath301 , remains close to the one of @xmath294 ( which is uniformly bounded for all scales ) , up to terms that are proportional to @xmath303 and , therefore , are bounded by ( const.)@xmath309 ( here @xmath310 is the dimensional amplification factor arising from the scaling dimension @xmath311 of the chemical potential terms , while @xmath293 is the dimensional gain coming from the irrelevance of the interaction ) .
recalling that @xmath312 , we find that @xmath295(const.)@xmath296 , for all scales @xmath301 .    finally , once the chemical potential is fixed so that @xmath313(const.)@xmath296 , we immediately infer that the beta functions of @xmath274 and @xmath212 are bounded by ( const.)@xmath296 , as well : therefore , their flows converge exponentially fast , and the dressed values of @xmath274 and @xmath212 are analytic functions of @xmath0 , analytically close to their bare values .
the discussion in the previous section guarantees that , once the chemical potential is properly fixed , then the flows of the chemical potential , wave function renormalizations , and fermi velocity converge exponentially fast .
the values of the chemical potential , as well as of the dressed wave functional renormalizations , dressed fermi velocity , and dressed critical lines are expressed in terms of convergent expansions ( they are _ analytic _ functions of @xmath0 ) , which are dominated by the first non trivial order in perturbation theory , provided @xmath0 is not too large ( note that the condition of convergence of the renormalized expansion is _ uniform _ in the gap , and is valid , in particular , _ on the critical line _ ) .
the explicit lowest order contributions to the chemical potential @xmath115 , to the renormalized fermi velocity @xmath314 and the wave function renormalizations @xmath315 on the renormalized critical line @xmath271 are the following :    1 .
chemical potential : = -_=1 ^ 2g_(+_f^_2)g_()g_(+ ) ; 2 .
fermi velocity : v_r=32 t_1-iu^2_k_1g_12(+_f^_2)g_12()g_21(+);[eq : vr ] 3 .
wave function renormalizations : z_,r=1+iu^2_k_0g_(+_f^_2)g_()g_(+).[eq : zrr ]    moreover , the equation for the critical line @xmath271 reads : m__2=2,where @xmath316 , @xmath317 and @xmath318 where defined after .
this is a fixed point equation for @xmath319 , whose solution leads to the plot in fig.[figaltezza ] .
note that , as discussed in sect.[sec.symm ] , there is no symmetry reason why @xmath34 should be equal to @xmath35 .
actually , an explicit computation shows that @xmath320 is different from zero along the critical line , unless we are at one of the highly symmetric points @xmath321 or @xmath192 , see fig.[fig1 - 2 ] , where we plot the value of @xmath320 on the critical line at second order in @xmath0 , for two different values of @xmath0 .
in this section we compute the jump discontinuity of the hall conductivity across the critical line , as well as the value of the longitudinal conductivity on the same line , and prove a universality result for both of them , i.e. , we prove that their values are quantized and exactly independent of the interaction strength @xmath0 .
note that this fact is highly non - trivial , due to the unusual renormalization of the fermi velocity and of the wave function renormalizations , which depends explicitly on the spinor index and break the asymptotic relativistic invariance of the propagator : the cancellations behind universality need to take lattice ( and , therefore , rg - irrelevant ) effects into account , and do not follow from asymptotic relativistic computations .
we stress that our result is exact at all orders of the ( convergent , renormalized ) expansion for the conductivity .
one key ingredient used in the proof is the lattice ward identity , which is rigorously valid ( without any sub - leading correction ) , thanks to the exact lattice symmetries and the fact that the correlations appearing at both sides can be computed in terms of convergent expansions , following from the multi - scale construction described above .
here we compute the universal jump discontinuity of the hall conductivity across the renormalized critical line . for the moment , we assume not to be at the graphene points @xmath322 and @xmath190 , @xmath323 ; we shall discuss later the ( straightforward ) adaptation to these special cases .
therefore , the goal is to compute : = _ m_r0^+ _ 12 - _ m_r0 ^ -_12,where @xmath324 is the mass gap of the dressed propagator .
the condition that we are not at a graphene point means that @xmath325 should be kept finite as @xmath326 . using the definition , as well as the fact that @xmath327 is differentiable in @xmath100 outside the critical line , we can rewrite @xmath328.\ ] ]
the interacting current - current correlation can be computed via the multiscale renormalized expansion discussed in sect .
[ sec.rg ] : in particular , proceeding as in ref.@xcite , among the contributions to @xmath329 we can distinguish the dominant contribution , coming from the ` dressed bubble ' , from the sub - dominant one , which is the sum over all the renormalized diagrams with at least one interaction term .
thanks to the irrelevance of the interaction , these sub - dominant diagrams have a dimensional gain ( of order @xmath258 on scale @xmath213 ) , which makes the corresponding contribution to @xmath327 _ differentiable _ at @xmath101 , in the limit @xmath326 . in particular , they give zero contribution to @xmath330 .    the dominant contribution to @xmath327 ( i.e. , the ` dressed bubble ' ) is k_ij^()&=&-2\{s_2()_i ( , ) + & & s_2(+)_j(+,-)},where @xmath331 is the vertex function defined in , and the factor 2 in front of the integral takes into account the spin degrees of freedom .
both @xmath332 and @xmath333 are given by convergent renormalized series , which depend on the details of the microscopic model .
the finite contribution to the jump - discontinuity of @xmath334 across @xmath335 comes from the integration over @xmath174 in the vicinity of @xmath336 , since the rest is continuous as @xmath326 . for the same reason , for the purpose of computing @xmath330
, we can replace @xmath337 by @xmath338 , and @xmath332 by its linearization @xmath339 at @xmath336 , |s( ) =    -ik_0 z_1,r+m_r & -v_r(-ik_1+_2 k_2 ) + -v_r(ik_1+_2 k_2 ) & -ik_0 z_2,r - m_r    ^-1,[eq : linprop ] where @xmath340 and @xmath32 are _ analytic _ functions of @xmath0 , for @xmath0 small , whose expansions at second order in @xmath0 are given explicitly by - .
recall that , a priori , @xmath341 are complicated infinite series in @xmath0 . thus , a direct computation of the jump - discontinuity , starting from the expression of the dressed bubble and from the feynman rules for the generic term in the renormalized expansions for @xmath340 , @xmath32 and @xmath342 , would be hopeless .
the key fact is that , thanks to the ward identity , [ eq : wivertex ] _ i(_f^_2)=_k_i|s^-1(*0 * ) that is , [ eq : wivertex2 ] _ 1(_f^_2)=-v_r_2,_2(_f^_2)=-_2v_r_1 .
therefore , & & = ( _ m_r0^+-_m_r0 ^ -)_|k| _ r + & & \ { |s( ) _
1|s^-1(*0 * ) _ 0|s()_2|s^-1(*0*)},[eq : delta]where we used that @xmath343 , and we denoted by @xmath288 a small , arbitrary , positive constant . using the identity _ 0 |s( )    @xmath344 by @xmath345 ( which is allowed , for the purpose of computing @xmath330 , simply because the difference is continuous at @xmath335 ) , we can further rewrite @xmath346 as & & = -(_m_r0^+-_m_r0 ^ -)_|k| _ r + & & \ { |s( ) _ 1|s^-1(*0 * ) |s()_0|s^-1(*0*)|s()_2|s^-1(*0*)}.the integral over @xmath25 can be evaluated explicitly and , after a straightforward computation , we get & & = _ m_r0^+m_r + & & _
|k| dk^-3/2.thus , introducing [ eq : vmtilde ] v_r = , m_r = m_r , we see that @xmath330 can be rewritten as , performing the change of variables @xmath347 : & = & _
m_r0^+ m_r _
|k|v_r dk^-3/2 + & = & _ m_r0^+ _ |k|v_r / m_r dk^-3/2 + & = & , where we recall that the result is expressed in units such that @xmath97 . therefore , the cancellation between the parameters @xmath32 , @xmath34
, @xmath35 gives a universal result .
finally , at the graphene points , the analogous computation gives twice the same value , because of an extra factor 2 coming from the valley degeneracy .
a similar discussion as the one in the previous subsection can be repeated for the _ longitudinal _ conductivity on the renormalized critical line .
the point here , as compared to the computation of @xmath330 in the previous subsection , is to take first the limit @xmath326 , and then @xmath348 ( recall the definition of conductivity , eq . ) .
once again , we assume for definiteness not to be exactly at the graphene point ( a similar discussion applies there , too ) .
note that , by the very definition of current - current correlations , @xmath349 is _ even _ in @xmath350 .
therefore , all the contributions to @xmath349 that are differentiable in @xmath350 give zero contribution to the longitudinal conductivity on the critical line . by repeating a strategy analogous to the one that led us to , for the purpose of computing the longitudinal conductivity on the critical line , we can : ( i )
replace the full current - current correlation by its dominant contribution ( from the ` dressed bubble ' ) ; ( ii ) restrict the integration over the loop momenta in the vicinity of @xmath336 ; ( iii ) linearize the propagators and vertex functions around @xmath336 ; ( iv ) use the ward identity eq .
( [ eq : wivertex ] ) to replace the vertex functions by the derivatives of the inverse two - point function .
after these replacements , we get ( denoting the value of the longitudinal conductivity on the critical line by @xmath351 ) : @xmath352,\ ] ] with f(,p_0)&=&\ { |s_0( ) _
i |s_0 ^ -1(*0 * )    & = & v_r^2\ { |s_0( ) _
i |s_0(k_0+p_0,k ) _ i } where @xmath353 is the linearized propagator , computed at @xmath335 , and the last step follows from ( [ eq : wivertex ] ) , ( [ eq : wivertex2 ] ) . by evaluating the integral over @xmath25 explicitly , and setting @xmath354 as in eq .
( [ eq : vmtilde ] ) , the computation of @xmath351 reduces to the contribution of just one dirac cone to the longitudinal conductivity of noninteracting graphene @xcite , with fermi velocity @xmath355 .
thus , proceeding as in ref.@xcite , we get , in units such that @xmath97 : _
ii^ = 12_p_00^+_0^v_r dx=18 .
notice that , as for graphene , the fermi velocity ( in general a nontrivial function of the hubbard interaction strength @xmath0 ) disappears , thus yielding a universal result .
the analogous computation performed at the graphene points gives twice the same value , in agreement with the result of ref.@xcite .
we studied the haldane - hubbard model by rigorous renormalization group techniques .
our analysis predicts that the critical lines separating the distinct topological phases are modified non - trivially by the hubbard interaction , in particular that the non - trivial topological phase , characterized by the topological quantum number @xmath41 , is enlarged by weak repulsive interactions .
moreover , our results rule out the presence of new interaction - induced topological phases in the vicinity of the phase boundaries .
such predictions may be verified experimentally in optical lattice realizations of the system@xcite , where the on - site interaction can be produced and tuned by means of feshbach resonances . concerning numerical simulations ,
our results agree with those of ref.@xcite .
the interaction affects the relativistic structure of the two - point function by non - universal renormalizatized coefficients , which differ from those obtained by approximate treatments of the system based on the effective dirac theory .
in particular , we find that there are two different wave function renormalizations , one for each pseudo - spin index . despite the non - universal renormalization of the two - point function and of the vertex functions , lattice ward identities guarantee the quantization and the universality of the conductivity matrix at the critical line .
concerning the transverse conductivity @xmath11 , its quantization follows from topological arguments ; however , these arguments do not provide any information regarding which values @xmath11 might take .
for instance , numerical and mean - field analyses predict that , at intermediate coupling strengths , new topological phases might appear , corresponding to the values @xmath356 , which are not present in the noninteracting theory .
our exact analysis rules out such new phases at small coupling .
the second part of our result focuses on the critical longitudinal conductivity @xmath357 ( away from criticality @xmath358 is trivially zero ) . in constrast to @xmath11 ,
this quantity is not protected by any topological argument .
nevertheless , we show that it is universal : all interaction and lattice corrections disappear .
each dirac cone contributes with a universal quantum of conductivity @xmath55 ; in particular , at the doubly critical points where the two critical curves cross ( see fig . [ figaltezza ] ) , the critical longitudinal conductivity is @xmath54 , which is the same value measured in graphene @xcite .
our results require the interaction to be weak and short - range ; instead , different features are expected in the presence of long - range interactions .
for instance , it is known that , at the graphene point , long - range interactions have dramatic effects on several physical properties @xcite , and their role on the renormalization of the optical conductivity is still actively debated @xcite .
we expect such effects to have profound implications for the haldane - hubbard model , especially in the proximity of the critical lines separating the different topological phases .
we plan to investigate this issue in future work .
the work of a.g . has been carried out thanks to the support of the a*midex project hypathie ( no anr-11-idex-0001 - 02 ) funded by the `` investissements davenir '' 25 french government program , managed by the french national research agency ( anr ) , and by a c.n.r.s .
visiting professorship spent at the university of lyon-1 .
the work of m.p . has been carried out thanks to the support of the nccr swissmap .
in order to define the current , we couple the electron gas to an external vector potential @xmath359 , by multiplying the hopping strength from @xmath360 to @xmath361 by an extra phase factor @xmath362 ( peierls substitution ) .
we denote by @xmath363 the modified hamiltonian , and let the ( paramagnetic ) current be @xmath364 , where @xmath47 label the two ( orthogonal ) coordinate directions @xmath365 and @xmath366 . an explicit computation leads to , with @xmath367 and , defining @xmath368 , & & m_11(k ,
p)=-it_2_j=1 ^ 3_=_j_p_je^i(-k_j ) , + & & m_12(k ,
p)=-it_1_j=1 ^ 3_j_p_je^-ik(_j-_1 ) , + @xmath369 and @xmath370 .
in this appendix , we discuss some of the details of the numerical computations from which figs.[figaltezza]-[fig1 - 2 ] were produced .
the program used to carry them out is available online @xcite , has been named hhtop , and is released under an apache license .
the source code includes a documentation file , in which the computations are described in greater detail .
the numerical computations carried out in this work involve numerical evaluations of integrals .
the algorithm that was used to carry these out is based on _ gauss - legendre _ quadratures , by which , given an integer @xmath371 , an integral is approximated by a discrete sum with @xmath372 terms : @xmath373 where @xmath374 are the roots of the @xmath372-th legendre polynomial @xmath375 , and @xmath376    if @xmath118 is an analytic function , then one can show that the remainder @xmath377 decays exponentially in @xmath372 .
however , in order to compute the difference of the wave - function renormalizations , we need to compute the integral of an integrand that , instead of being analytic , is a class-2 gevrey function ( a class-@xmath378 gevrey function is a @xmath379 function whose @xmath216-th derivative is bounded by @xmath380 , so that analytic functions are class-1 gevrey functions ) .
the remainder @xmath377 can be shown to be bounded , if @xmath118 is a class-@xmath378 gevrey function with @xmath381 and @xmath372 is large enough ( independently of @xmath118 and @xmath378 ) , by @xmath382 for some @xmath383 , that only depend on @xmath118 . for a proof of this statement ,
see lemma  a3.1 in the documentation of hhtop@xcite . in short , this estimate is obtained by expanding @xmath118 in chebyshev polynomials , and using a theorem of a.c .
curtis and p. rabinowitz@xcite that shows that , if @xmath118 is the @xmath384-th chebyshev polynomial , then @xmath377 is bounded uniformly in @xmath384 .
the decay of the coefficients of the chebyshev expansion of class-@xmath378 gevrey polynomials allows us to conclude .
there is a single , minor , pitfall in the numerical evaluation of @xmath387 : we wish to use gauss - legendre quadratures ( see app.[app.gauss_legendre ] ) to carry out the computation , but the integrand in ( [ eq.fnumeric ] ) is not smooth : indeed , if @xmath388 , then its second derivative diverges at @xmath389 due to the divergence of the derivative of @xmath390 .
however , by switching to polar coordinates @xmath391 , this singularities is regularized , that is , the integrand becomes a smooth function of @xmath392 and @xmath393 . at this point
, there is yet another danger to avoid : while the integrand is smooth , the upper bound of the integral over @xmath392 is a function of @xmath393 , which is , due to the rhombic shape of @xmath394 , only smooth by parts .
the integral over @xmath393 must , therefore , be split into parts in which the bounds of the integral over @xmath392 are smooth .
this can be done very easily using the @xmath395 rotation symmetry .
once both of these traps have been thwarted , gauss - legendre quadratures yield very accurate results .    in order to compute the correction to the critical line
, we solve @xmath396 for @xmath397 and @xmath398 . for the sake of clarity ,
we have made the @xmath16 dependence of @xmath387 explicit . to solve ( [ eq.mass1 ] ) , we fix @xmath398 , and use a newton algorithm to compute the critical value of @xmath397 : we set @xmath399 , and compute @xmath400 provided @xmath401 is not too far from the solution of ( [ eq.mass1 ] ) , @xmath402 converges _ quadratically _
( i.e. @xmath403 , in which the constant depends on the supremum of @xmath404 , which is bounded ) to the solution of ( [ eq.mass1 ] ) .      at second order in @xmath0
, @xmath320 is @xmath405 where @xmath406 the computation is carried out on the critical line , that is , when @xmath407 .
the integrals over @xmath350 and @xmath408 can be carried out explicitly : @xmath409\hskip15pt\cdot      \left(\frac{(\xi_p+\xi_q+\xi_f)(\frac{m_p}{\xi_p}+\frac{m_q}{\xi_q}-\frac{m_f}{\xi_f}-\frac{m_pm_qm_f}{\xi_p\xi_q\xi_f})z}{(z^2-(\xi_p+\xi_q+\xi_f)^2)^2}\right )    \end{array }    \label{eq.z1-z2}\ ] ] where , using the definitions of @xmath317 , @xmath410 and @xmath109 after and after , @xmath411 , @xmath412 , @xmath413 , @xmath414 , @xmath415 and @xmath416 , @xmath417 , @xmath418 .
the numerical evaluation of the integral in ( [ eq.z1-z2 ] ) involves a similar difficulty to that in ( [ eq.fnumeric ] ) : the integrand has divergent derivatives if any of the following conditions hold : @xmath419 , @xmath420 or @xmath421 .
these singularities can not be regularized by changing @xmath422 and @xmath423 to polar coordinates , since @xmath424 is a singular function of the polar coordinates of @xmath422 and @xmath423 ( due to the fact that it behaves , asymptotically , as @xmath425 approaches @xmath426 , as @xmath427 , which has divergent second derivatives ) . however , there are coordinates , which we call _ sunrise coordinates _ ( since @xmath428 is the value of the so - called _ sunrise _ feynman diagram ) , which regularize these singularities .
their expression is rather long , and will not be expounded here ; the interested reader is invited to consult the documentation file bundled with the source code of hhtop@xcite . once written in terms of the sunrise
coordinates , the integral in ( [ eq.z1-z2 ] ) can be computed using gauss - legendre quadratures very accurately .        c. l. kane and e. j. mele , phys .
95 , 146802 ( 2005 ) .
j. e. moore and l. balents , phys .
b 75 , 121306(r ) ( 2007 ) .
r. roy , phys .
b 79 , 195321 ( 2009 ) .
l. fu and c. l. kane , phys . rev .
b 74 , 195312 ( 2006 ) .
d. carpentier , p. delplace , m. fruchart and k. gawedzki , phys .
114 , 106806 ( 2015 ) .
d. carpentier et al . , nucl .
b 896 , 779 - 834 ( 2015 ) .
a. kitaev , in advances in theoretical physics : landau memorial conference , edited by v. lebedev and m. feigelaman , aip conf . proc.,vol .
1134 ( aip , melville , ny , 2009 ) , pp .
a. p. schnyder , s. ryu , a. furusaki , and a. w. w. ludwig , phys . rev .
b 78 , 195125 ( 2008 ) .
s. raghu , x. qi , c. honerkamp , s. zhang , phys .
100 , 156401 ( 2008 ) . c. weeks , m. franz , phys
b 81 , 085105 ( 2010 ) .
p. ghaemi , j. cayssol , d. n. sheng , a. vishwanath , phys . rev .
108 , 266801 ( 2012 ) .
d. a. abanin and d. a. pesin , phys .
109 , 066802 ( 2012 ) .
j. gonzlez , j. high energ .
07 , 175 ( 2013 ) .
b. roy and i. herbut , phys .
b 88 , 045425 ( 2013 ) .
m. daghofer and m. hohenadler , phys .
b 89 , 035103 ( 2014 ) .
j. he et al .
, phys . rev .
b 84 , 035127 ( 2011 ) .
j. he , y. liang , s .-
kou , phys .
b 85 , 205107 ( 2012 ) .
zhu et al .
, j. phys .
: cond . mat .
26 , 175601 ( 2014 ) .
wu , n. li , and s .- p .
kou , eur .
j. b 88 , 255 ( 2015 ) . w. zheng , h. shen , z. wang and h. zhai , phys .
b 91 , 161107(r ) .
gu , k. li and j .- x .
li , arxiv:1512.05118
. j. wu , j. p. l. faye , d. senechal and j. maciejko , phys .
b 93 , 075131 ( 2016 ) .
j. he , s .-
kou , y. liang and s. feng , phys . rev .
b 83 , 205116 ( 2011 ) .
j. maciejko and a. ruegg , phys . rev .
b 88 , 241101(r ) ( 2013 ) . c. n. varney , k. sun , m. rigol and v. galitski , phys
b 82 , 115125 ( 2010 ) ; phys . rev .
b 84 , 241105(r ) ( 2011 ) .
wang , z .- c .
gu , c .- d . gong and d. n. sheng , phys .
107 , 146803 ( 2011 ) .
i. vasic , a. petrescu , k. le hur and w. hofstetter , phys .
b 91 , 094502 ( 2015 ) .
a. amaricci et al . , phys .
114 , 185701 ( 2015 ) ; phys .
b , in press ( arxiv:1603.04263 ) .
p. goswami and s. chakravarty , phys .
104 , 196803 ( 2011 ) ; arxiv:1603.03763 .
k. kobayashi , t. ohtsuki , k .-
imura and i. f. herbut , phys .
112 , 016402 ( 2014 ) .
b. roy , p. goswami and j. d. sau , phys .
b 94 , 041101 ( 2016 ) ; b. roy , y. alavirad and jay d. sau arxiv:1604.01390 .
e. g. mishchenko , europhys .
83 , 17005 ( 2008 ) .
i. herbut , v. jurii and o. vafek , phys .
100 , 046403 ( 2008 ) .
v. jurii , o. vafek and i. herbut , phys . rev .
b 82 , 235402 ( 2010 ) .
d. e. sheehy and j. schmalian , phys .
b 80 , 193411 ( 2009 ) . a. giuliani and v. mastropietro , phys .
b 85 , 045420 ( 2012 ) .
i. sodemann and m. m. fogler , phys .
b 86 , 115408 ( 2012 ) .
b. rosenstein , m. lewkowicz , t. maniv , phys .
110 , 066602 ( 2013 ) .
d. l. boyda , v. v. braguta , m. i. katsnelson , m. v. ulybyshev , phys . rev .
b 94 , 085421 ( 2016 ) .
j. m. link , p. p.
orth , d. e. sheehy and j. schmalian , phys . rev .
b 93 , 235447 .
i. herbut , v. jurii , o. vafek and m. j. case , arxiv:0809.0725 .
the definition of @xmath433 is meant as the limit as @xmath434 of @xmath435 , which is the anti - periodic extension in @xmath436 and @xmath437 ( of anti - period @xmath75 in both variables ) of @xmath438 , with @xmath439 . at equal times
, we let @xmath440 .
the same conventions hold for the mixed current - current and current - field correlations introduced below .