Aspects of the pseudo chiral magnetic effect in 2D Weyl-Dirac matter

A connection is established between the continuum limit of the low-energy tight-binding description of graphene immersed in an in-plane magnetic field and the Chiral Magnetic Effect in Quantum Chromodynamics. A combination of mass gaps that explicitly breaks the equivalence of the Dirac cones, favoring an imbalance of pseudo-chiralities, is the essential ingredient to generate a non-dissipative electric current along the external field. Currents, number densities and condensates generated from this setup are investigated for different hierarchies of the energy scales involved.


Introduction
The chiral magnetic effect (CME), predicted in Refs. [1,2] to occur in the quark gluon plasma produced in heavy ion collisions (HIC), has a deep connection with the vacuum structure of quantum chromodynamics (QCD), its topology and symmetries. The non-dissipative current produced by this mechanism points out in the direction of the magnetic field generated in non-central collisions and is a direct consequence of a chiral imbalance. Such imbalance can only occur in certain domains where the gauge field configurations are topologically non-trivial. Interactions of fermions with these fields result in a chirality flip.
The difficulty to extract information about the early stage in HIC -due to unknown electromagnetic properties of the medium, the out-of-equilibrium regime and the lack of effective transport descriptions -makes the CME so far an unique attempt to trace a mechanism that connects some peculiar a e-mail: cvillavicencio@ubiobio.cl quantum properties of QCD to a macroscopic observable. Besides contributing for a more complete picture of QCD, its confirmation would have implications in our understanding about the early universe and baryogenesis. Because of that, it has received considerable attention from the theoretical, experimental and lattice communities in high energy physics, but in spite of all efforts, a conclusive observation of the mechanism is still missing. In effect, results from the CMS collaboration [3] comparing angular correlation between Pb-Pb collisions and p-Pb collisions have challenged previous results from the STAR collaboration [4] that seemed to have observed the CME. In order to clarify this issue, new technology on background analysis is being carried out [5] and new observables have been proposed [5,6]. Remarkably, the 2018 run of RHIC includes isobar collisions in order to disentangle the chiral magnetic effect from background sources [7].
In a different order of ideas, the fast-growing family of Weyl-Dirac materials that have been discovered in the last few years has allowed to test an analog of this mechanism in a condensed matter environment [8]. In this kind of materials, the complex interaction between the charge carriers and the background lattice can be effectively represented considering the former as quasiparticles that obey relativistic-like equations of motion, with the velocity of light replaced by the corresponding Fermi velocity. In this way, it is possible to define chirality for these electrons and to construct the analogy with QCD [9][10][11].
In search for inducing a chiral splitting of charge carriers in Weyl-Dirac materials, an experiment was proposed in Ref. [8] such that by applying a parallel electric and magnetic field to a (3+1)D sample of zirconium pentatelluride (ZrTe 5 ), the observation of a negative magnetoresistence signals the presence of a chiral anomaly [12] and the generation of a nondissipative electric current. This represents yet another novel avenue allowing for interdisciplinary investigation connecting condensed matter and high energy physics. Nevertheless, the Weyl-Dirac behavior of the charge carriers in ZrTe 5 is still controversial. While the authors of Ref. [8] established the ultrarelativistic behavior of charge carriers through angle-resolved photoemission spectroscopy (ARPES), further investigation has opened the possibility that this may not be actually the case after all [13]. A more detailed understanding of the mechanism in the context of condensed matter is still needed and realizations of CME in different materials are welcome. In effect, it has been posteriorly observed in other materials [14][15][16][17][18] and alternative mechanisms have been proposed in order to generate the CME in (3 + 1)D Weyl semimetals [19] but no experimental realization has been achieved so far.
In this paper, we extend the work presented in Refs. [20,21], where some of us proposed an electromagnetic analogue of the CME in Weyl-Dirac materials in two spatial dimensions motivated by graphene and referred to as the pseudo-chiral magnetic effect (PCME).
It is natural to search for this type of analogy because quantum electrodynamics in (2+1)D, dubbed as QED 3 , has been widely used as a toy model for QCD inasmuch as it describes confinement and chiral symmetry breaking [22][23][24][25][26]. This is due to the fact that at high temperature, any field theory can be dimensionally reduced and, on the other hand, a non-Abelian three dimensional gauge field theory abelianizes for a large number of flavors [27]. With the physical realization of graphene and other planar materials containing Dirac electrons, QED 3 was promoted from being a toy model to really describe physical systems and analogues of high energy physics constructed on table-top experiments started to take place.
To describe the PCME, in Refs. [20,21] we constructed a Lagrangian for relativistic fermions in (2+1)D including some effective interactions that simulate a chiral imbalance. The result is an electric current generated in the direction of an external magnetic field applied in-plane. Here we extend this study investigating other relevant quantities: the number densities and condensates that we found to be present and the current in certain limiting cases of magnitude of the magnetic field as compared with temperature, mass gaps and chemical potentials.
The paper is organized as follows: in Sect. 2 we give a brief review of how Dirac fermion emerge from graphene and in this context we present our Lagrangian and describe our model. In Sect. 3 we discuss the propagator of a fermion in the presence of a magnetic field based on the Schwinger proper time method. In Sect. 4 we present our results on the currents, densities and condensates. In Sect. 5 we explain our model in terms of the Fermi liquid model, commonly applied in condensed matter physics. In Sect. 6 we make our final remarks. We complement the article with three appendixes.

The effective model
The essence of the CME lies on the fact that the vacuum of the gauge sector in QCD is actually a superposition of vacua including configurations with non-trivial topology, and the interaction of these gauge fields with quarks can cause a flip on their chirality. In that case, an imbalance between chiralities occurs and, inside a certain domain where the topological gauge fields act, one has a difference in the number of leftand right-quarks, whose abundance and sign depend on the details of the topology of the gauge fields. The large magnetic field present in non-central heavy-ion collisions triggers the CME by generating an electric current along its direction. Therefore, the chiral imbalance and the magnetic field constitute the two pillars of the CME.
Planar graphene-like materials provide elements we need for a robust analogy. This kind of materials have the remarkable property that the charge carriers behave as Dirac electrons, which means that they have a linear dispersion relation. Besides, the continuum limit of its tight-binding Hamiltonian (or Lagrangian) exhibit Time-and Space-inversion symmetries such that massless charge carriers are, in fact, chiral. Theoretically, it is possible to break these symmetries if we consider special currents in the Lagrangian, as will be clear in a moment.
Let us briefly describe the low-energy continuum dynamics of graphene-like materials. For a more detailed description one counts on several books and reviews on the subject, including Refs. [9][10][11]28] and references therein. The honeycomb array of these materials does not have the structure of a Bravais lattice required for a tight-binding description. Nevertheless, it can be regarded in terms of two equivalent triangular sub-lattices A and B which are Bravais lattices, in which the nearest neighbor of a carbon atom belonging to sub-lattice A are all atoms in the sub-lattice B and vice versa, as shown in Fig. 1. With the primitive vectors a 1 = a(1/2, √ 3/2), and a 2 = a(−1/2, √ 3/2), where a denotes the interatomic distance, and correspondingly in reciprocal space K 1 = 2π/a(1, √ 3/3) and K 2 = 2π/a(−1, √ 3/3), from Bloch theorem it is straightforward to find that the energy-momentum dispersion relation is where t is a constant that encodes the magnitude of the interaction (the hopping parameter) and k i correspond to the momentum components in the two spatial directions. This function is plotted in Fig. 2. One can notice that the valence and conduction band touch in six points, generating no gap. Those are high symmetry points in the Brillouin zone of which only two are inequivalent, the so-called Dirac points K and K . Expanding the energy around these points, one can verify that the dispersion relation is linear, E =hv F k, and we can think of charge carriers as massless Dirac electrons, with the speed of light replaced by the Fermi velocity. Therefore, the tight-binding description of monolayer graphene corresponds in the continuum to a massless version of the fermion sector of QED 3 , but photons are allowed to move in (3+1) dimensions. The influence of external electromagnetic fields as usual is added via minimal couplings.
The symmetry between the two sublattices allows to represent the charge carriers as four component spinors, with two indices corresponding to the sub-lattice flavor or pseudo-spin and the Dirac point (valley) index, respectively. Working in Weyl representation, the gamma matrices can be explicitly written: Notice that in this case, the free massless Dirac Lagrangian with ∂ = iγ μ ∂ μ and μ = 0, 1, 2, is invariant under the chiral-like transformations Preserving the symmetries in Eq. (4), a couple of terms, bilinear on the fermion fields, can be added, m oψ γ 3 γ 5 ψ and m eψ γ 3 ψ. The first term is known as Haldane mass and in QED 3 it is directly related to a Chern-Simons (CS) effective action upon integrating fermionic degrees of freedom. The second, proportional to m e , corresponds to sublattice symmetry breaking.
Let us now discuss the role of the parity anomaly and how it connects to the aforementioned masses. The ABJ anomaly in (3+1)D can be interpreted as production of Weyl fermions in the presence of external electric and magnetic fields [29]. This was initially proposed in a nuclear physics context and in [30] this interpretation was extended to condensed matter systems. The authors show that the contribution to the anomalous electric current coming from fermions with different chiralities have opposite sign and, since it was previously shown that in any lattice theory with locality, chiral fermions appear in pairs of opposite chirality [31] -fermion doubling -the net current vanishes.
In odd dimensions, it is known that the conservation laws for fermionic currents do not present the chiral anomaly, but it might still emerge abnormal currents for each species of fermions due to the parity anomaly [32]. In [33] the author proposes a (2+1)D analogue of [30] in a graphene-like system. Introducing a gap between the bands induced by a sublattice symmetry breaking, the author show that the contribution from the two inequivalent Dirac points cancel out, similarly to what happens in (3+1)D. The two Dirac points are therefore the analogues for chirality.
In our work we propose something similar, but we look for a Lagrangian structure where the masses for the two valleys are different and therefore for suitable values of the parameters the net current does not vanish. We obtain this by the combination of Haldane mass m o and the mass m e , which carries an identical structure as the one Semenoff uses in [33]. The anomaly is therefore parametrized in m ± in the same way that the chiral anomaly is parametrized in μ 5 in the original chiral magnetic effect. The underlying Hamiltonian is therefore, equivalent to the one in [33] but with different masses for each Dirac point. Further, as we explain below, we promote the interaction of this system with an external magnetic field via minimal coupling and consider thermal effects.
To complete our analogy, we consider an external magnetic field aligned along the graphene membrane described by the vector potential A ext 3 = By, where y represents the second spatial coordinate along the graphene plane and we assume B > 0. This field is assumed to be classical and does not play a role in quantum corrections. Since it is external, in our setup A μ lives in a bulk rather than in a twodimensions sheet. Other field theory approaches for QED in planar condensed matter systems have been presented previously considering an explicit treatment of the gauge sector, see for instance the Pseudo Quantum Electrodynamics (PQED) [34][35][36][37] and the Reduced Quantum Electrodynamics (RQED) [38]. In the present work we do not deal with the Maxwell term and do not consider quantum corrections coming from the gauge sector and therefore we do not make use of PQED/RQED. However, in a work in progress we are currently investigating if within this more complete framework m o can be dynamically generated by the presence of parallel electric and magnetic fields. One of us have checked some formal aspects of the theory in mixed dimensions [39,40] in order to prepare for the numerical calculations and we expect to report soon. In our approach we just consider the fermion sector and the interaction with the classical field. We restrict the dynamical term ∂ψ/∂x 3 = 0 but we preserve the fermion interaction with A 3 .
The scenario described above can be represented by the Lagrangian: where / D = (∂ 0 − iμ, v F ∇), e = −|e| is the fundamental charge, v F is the Fermi velocity of quasiparticles in the crystal which from now onward we set in the natural units of the system (namely, v F = 1), and μ the chemical potential. In the Weyl representation for the gamma matrices, the matter content of the theory can be unveiled by introducing the chiral-like projection operators χ ± = 1 2 (1 ± γ 5 ) which verify χ + + χ − = 1, and χ 2 ± = χ ± , and thus allow us to write where m ± = m e ± m 0 and the "left-" and "right-handed" fields are ψ ± = χ ± ψ. Hence, the Lagrangian represents a system where the two different charge carrier species become non-degenerate in gaps: one of them becomes lighter and the other heavier as m o grows (see Fig. 3). In order to have access to currents, densities and condensates arising in our system, we require the corresponding Green function for fermions from the Lagrangian in Eq. (5), which already takes into account the influence of the in-plane external magnetic field, the chemical potential and temperature. We start by splitting the complete Green function into separate Green functions for each chirality, G ± , such that Since we are interested in exploring this effect at room temperature we must include thermal effects. This is introduced in the usual way by replacing the zeroth component of the momentum with the Matsubara frequencies k 0 → iω n = i(2n + 1)π T and momentum integrals are replaced by sums over these frequencies, Secondly the inclusion of the external magnetic field can be treated within the Schwinger proper time method [42].
Recall, however, that in the presence of finite chemical potential, some care must be taken in the range of integration in the proper-time integrals in order to guarantee their correct convergence [43]. This can be done by the mnemonic rule where the regulator reads Thus, we express the resulting propagator for each chirality as [20]: with where we made use of the shorthand notation M ± = 1 2 y + y eB + m ± with (r 1 , r 2 , r 3 ) ≡ (x, y, z). We also define K = iω n + μ, k 1 , 0 which is the parallel component to the magnetic field, oriented in the x-direction. Notice that the Green function is non-local due to the term M ± .
In the particular case where μ < πT , the proper-time integral can be Wick-rotated in a simple way (see [20] for more details). In this case, the propagator in Eq. (11) can be written as So, for now onward we consider this regime, where the system is diluted enough to fulfill the condition μ < πT . Below we derive expressions for currents, densities and condensates accessible from this Green function.

Currents, densities and condensates
From the full Green function in Eqs. (7), (10) and (11), we are in position of obtaining several physical observables as functions of the magnetic field B, temperature T , chemical potential μ as well as the charge carrier species gaps m ± with different hierarchies among these scales. The nonlocal nature of the Green function in Eq. (11) indicates that finite size effects are important to be considered. The relevant expectation value of the objects can be defined with respect to a given origin, which we choose, without loss of generality, to be the center of the layer in the y-direction. The setup is symmetric with respect to the x-direction. The appropriate value can be considered then as the average in space, namelȳ Let us comment about the apparent dependence of the current on the coordinate y. A translation in our propagator can be written as A 3 = By → B(y − y 0 ). Note that this is in fact equivalent to shift the mass m e → m e + eBy 0 and therefore different integration limits yield to the same result as far as m e is correctly shifted. In the absence of magnetic field this parameter corresponds to the gap generated by sublattice symmetry breaking obtained for example placing graphene over a boron nitride substrate. If a magnetic field is turned on, m e plays the role of a bare mass that can absorb an eventual constant added to the gauge field and its physical value must be fixed by some observable like the electric current. This is similar for instance to the procedure adopted in [44], where a constant uniform electric field generated from a static potential produces a non-local term, which is considered as a coordinate-dependent chemical potential and must be fixed by a suitable quantity like the number density.
All the relevant quantities can be expressed in a closed analytical form, instead of integrals in proper-time. We explore two regimes for the magnetic field with respect to the temperature and chemical potential. More specifically, the relevant scale to which we compare the field strength is (π T ) 2 − μ 2 , namely, the quantity that needs to be positive to enable for a Wick rotation. Details of the asymptotic low-and highmagnetic field approximations are described in the appendices. The case of pure magnetic contribution for the condensates will be treated separately.

Currents
First we obtain the electric currents for each sub-lattice and the chiral current according to with i = 1, 2. Tracing over gamma matrices, one can notice that the only non-vanishing components of the currents are j 1 and j 15 along the magnetic field direction. This is the essence of the PCME. We express these currents as with y ± = m ± /eB and the function j (η) defined as The Matsubara sums can be performed straightforwardly using the definition of the Jacobi Θ 3 -function and then, using the inversion relation for the Jacobi function (see Appendix A) we can write The approximation of low and high magnetic field is obtained from Eq. (19). For |eB| where A comparison of the weak and strong field expansion and the full numerical result with all the parameters scaled with |eB| is shown in Fig. 4. We can observe that the approximations are in good accordance with the numerical expression. It is interesting that in fact for strong magnetic field (low temperature) the dependence of the currents on the temperature is negligible.

Number densities
We proceed the same as in the case of electric and chiral currents. The number density and the chiral number density are defined as and we can express these quantities as where the auxiliary function ν (η) is defined as Performing the Matsubara sums, it can also be expressed as The approximations for weak and strong limits of the magnetic field are obtained from Eq. (27). For |eB| (π T ) 2 − μ 2 , where Li 2 (x) is the dilogarithm function.
Notice that a very similar expression was derived in Ref. [45] Eq.(74), replacing eBη by a mass gap . This is totally expected since the relevant value to calculate the currents, as can be checked in Eq. (25), is η = y + y ± , where y ± = m ± eB . It means that the argument in the exponentials of the equation above is a function of m ± only, which plays precisely the same role as : breaking sublattice symmetry.
For |eB| (π T ) 2 − μ 2 , This is a remarkable result. Comparing with Eq. (22), we can see that | j| = |eν|. In fact, integrating over the plane, we get where J 1 and N 5 are j 1 and n 5 integrated over the plane. This is exactly the same result obtained in [2], for the CME in QCD. The comparison between the strong and weak magnetic field limits and the full function is described in Fig. 5 for all parameters scaled with |eB|.
Similarly, the chiral separation effect is naturally obtained in our prescription. It is defined as a chiral charge separation that appears in the case a (non-chiral) chemical potential is present [46]. In the context of the quark gluon plasma it can be intuitively understood as follows: right-handed particles align with the magnetic field while their right-handed antiparticle anti-align with the field. If there is a chemical potential there will be a net current of right-handed particles in the direction of the field. In the same way, a current of left-handed particles is generated in the opposite direction, contributing with the same magnitude to the chiral current. Observing Eqs. (30) and (22), we verify that:

Condensates
Finally, we proceed with a similar reasoning to obtain condensates. The non-vanishing condensates are generated through the gaps m e and m o . The reason to refer to these objects as condensates and not as a currents is because there is no dynamics along the z-direction. In this sense, σ 3 and σ 35 behave like order parameters for the symmetry breaking generated by the mass gaps in the sense that an expansion around these values naturally give a finite gap in the corresponding Lagrangian. It is convenient to write these condensates as where the function σ in this case is We must be careful because the σ function is UV divergent. Therefore we separate the condensate in three contributions: the thermal-magnetic, magnetic, and vacuum part: We start by analyzing the thermal-magnetic contribution, which is defined as where we have subtracted the T, μ = 0 part, as can be seen from the Θ 3 − 1 term.
In the same way we proceeded with the currents and number densities, we can obtain approximations for weak and strong magnetic field from equation (37), now taking care that we have to subtract the T, μ = 0 contribution. For eB For |eB| where we also consider the fact that T √ |eB|, using the non-relativistic gas approximation. Here, Li 1 2 (x) is the Polylogarithmic function.
The thermo-magnetic contribution to the condensates is described in Fig. 6, compared with the weak and strong magnetic field approximations. The strong magnetic field (low temperature) approximation seems to be zero, but it is slowly growing with the temperature. This is not appreciable in this graph. Now, the pure magnetic contribution is defined as where η = y + m χ /eB. The last term removes the vacuum part, so this function vanishes for B = 0. Until now, all the approximations have been made in terms of the magnetic field compared with the temperature. For the pure magnetic contribution, however, we have available two more energy scales to compare against, the masses (gaps) and the (inverse) of the length of the sample. Thus, to perform the proper-time integration and get a better comprehension about the scales for which the magnetic field strength can be considered weak or strong, we calculate the average of the condensate in space, defining it in the same way as in Eq. (14), Then, for |eB L y /2m χ | 1, and for |eB L/2m χ | 1, Finally, the third contribution is the vacuum part, which is divergent. Thus, we introduce a UV cutoff Λ to regulate the integrals, namely, assuming in the last step that Λ m χ .

Currents, densities and condensates at zero temperature
Although graphene is usually studied at room temperature, it is interesting to check the zero temperature limit. The calculation in this regime presents a technical issue whenever chemical potential and magnetic field are present in the Schwinger formalism: the Wick rotation in proper time is not trivial in this case, and therefore the propagator is highly oscillatory. However, we can explore what happens at small magnetic field. In this case, in Eq. (11) the trigonometric functions can be approximated as tan(eBs) ≈ eBs and sec(eBs) ≈ 1.
where in the last equation we consider the chemical potential dependent term only and do not consider the vacuum part of the condensate. It is interesting to see in this case that the magnetic conductivity is finite if only one pseudochirallity is present. For example if |m − | < μ < |m + |, the electric current density is j 1 = e 2 B/2π .

Chiral chemical potential
In QCD, the CME is described by a single parameter, the so-called chiral chemical potential, which is proportional to the time derivative of the CS term. To establish our analogy, we have described a 2D Weyl-Dirac material where the parity anomaly is manifested through the mass term m ± in the Lagrangian. But because we need the system to be filled with charge carriers, the chemical potential also plays a role in the dynamics. If we consider the fermion system of our material degenerated with respect to the energy gaps m ± , we can explore what happens with the charge carriers, or holes, moving near the Fermi surface. There are two scenarios in this case, which we consider separately. Let us turn-off the external magnetic field for a while. The Lagrangian (6) can now be written as From the equations of motion of each chiral field, we can obtain the dispersion relations for particles Thus, particles with opposite chiralities propagate differently. If the chemical potential is larger than the mass gaps, we can consider each species propagating near their Fermi surface (Fig. 7). In this case, antiparticles can be omitted in the discussion because those states are already filled by assumption. Now, the Fermi momentum for each chirality is: Expanding around it, we get where now + and − label quasi-particles and their antiparticles (quasi-holes), respectively. We have restored the velocity units to make clear the dynamics. Then, we observe that we must define a new Fermi velocity for each chirality, namely, Thus, quasi-particles and holes propagate at different velocities under this perspective.
We can see that the Lagrangian that describes the motion of our Fermi liquid system reads which describes a gapless fermion system, with quasi-particles and holes moving differently. Moreover, the corresponding fields in momentum space are also different, This can be seen as if the Dirac points described in Fig. 1 were displaced as K → K + p + and K → K + p − .
Let us explore this approximation in another way, now by assuming that all the quasi-particles and holes propagate with the same Fermi velocity v F . Thus, the Fermi momentum for each chirality can be defined as Expanding Eq. (51) around p χ we obtain In this way, we obtain the above dispersion relations from the effective Lagrangian with the new fields defined in Eq. (56), and where the chemical potentials are defined as At this moment, we have not specified the orientation of the Fermi momentum because the Fermi velocity is unknown. Certainly, the velocity of propagation of the charge carriers we consider is near the original average Fermi energy μ. Therefore, we can set μ = 0 and what we get is an explicit relation of the Fermi velocity and chiral chemical potential, From the above expressions, we observe that mass gaps can be parametrized in terms of this single chiral chemical potential, just as in the CME (Fig. 8). This effective Lagrangian corresponds to remove the masses and add a chemical potential ±μ 5 to each chirality.). In particular, in the case of zero temperature and small magnetic field, we have Let us emphasize once more that our construction is based on the symmetry breaking of Dirac Points which in those planar systems with relativistic fermions play the role of pseudochirality. As such, μ 5 encodes this asymmetry rather than being related to actual chirality.

Discussion and conclusions
In this article we have established an analogy of the CME for 2D Dirac-Weyl matter. We have considered a system of fermions with dynamics restricted to a plane which have two different mass gaps. One of those mass gaps is known to be related with the CS term in QED 3 . It is under investigation if the same holds for PQED or RQED. The other mass gap m e breaks the equivalence between the two triangular sublattices of the underlying honeycomb array. When an external magnetic field is applied parallel to this system, the generation of an electric current along the direction of the magnetic field is formed. This corresponds to an analogue of the CME, but in terms of the pseudo-spin of the charge carriers of our system. We have further explored the formation of number density, chiral number density, axial current and condensates associated to the mass gaps.
The most remarkable relation and close connection of the CME and the effect discussed here is represented in Eq. (31), which was obtained in the strong magnetic field (low temperature) regime and is functionally the same obtained for the ordinary CME in QCD [2]. Let us stress that in the latter case, the phenomenon is believed to occur at high-temperature. The CME is described through a chiral chemical potential, while in our case we need two mass gaps. This can be understood considering that in graphene-like system first we need a mechanism that distinguishes the left and right pseudochiralities and further another one that breaks the balance. Considering our Fermi system to be degenerate, it is possible to describe it as a Fermi liquid effective model, and the chiral imbalance can be seen in two ways: as quasiparticles with different pseudo-chiralities propagating at different Fermi velocities, or like a fluid whose imbalance is produced by a chiral chemical potential that depends on the difference of the mass gaps, as can be seen from Eq. (60). In both cases, the consequence is a shift in the Dirac points, hinting a deformation of the crystal structure. This observation suggests the fact that strained graphene is a good candidate to reproduce this effect.
It is interesting to remark that different arrangements of a graphene membrane can mimic many QCD situations. We want to explore in more detail the mass gap generation, including the possibility of generating it through external fields, which requires a proper treatment of the gauge sector. This idea is under scrutiny and results shall be presented elsewhere.
A useful trick used in this article is to include a momentum integral in order to remove some denominator in the Schwinger proper time, which allow us to integrate the proper analytically. For this reason, in some approximations, the proper-time denominator is exponentiated using the formula 1 s n/2 = 2 Γ n 2 ∞ 0 dp p n−1 e −sp 2 . (A.4)

Appendix B: Weak Field Approximation
For the weak field approximation of j (η), ν(η) and σ T,B (η), we first start from Eqs. (19), (27) and (37), respectively. First, we notice that for eB (π T ) 2 − μ 2 , we can set sech(eBs) 1 and tanh (eBs) eBs, allowing the exact integration in proper-time in the case of j 1 (η). For ν(η) and σ (η) there is a s −1 factor that can replaced by the integral in a momentum variable described in Eq. (A.4) with n = 2. After integration in proper time, the Matsubara frequencies can be summed exactly and then we perform the integral in terms of the auxiliary p in the case of number densities and condensates.
In the case of the pure magnetic contribution to the condensates, we average first on the space as is indicated in Eq. (43). Then we expand in powers of eB and finally integrate on the proper time. The expansion itself shows the scale of the expansion, which is m 2 χ /L y .

Appendix C: Strong Field Approximation
For the strong field approximation of j (η), ν(η) and σ T,B (η), we start as in the previous case from Eqs. (19), (27) and (37), respectively. Now, we consider the case of extremely intense field, eB (π T ) 2 − μ 2 . In this regime, we notice that tanh(|eB|s) ≈ 1, and in the case of the condensate, sech(eBs) ≈ 4e −2|eB| . This generates an denominator s −1/2 in the three cases, and we replace it by an integral in momentum described in Eq. (A.4) with n = 1. Now, the integral in proper-time can be solved exactly and also the sum in Matsubara frequencies. Then we integrate the auxiliary momentum term.
In the case of the thermomagnetic contribution for the condensates, the non-relativistic approximation was used. This can be applied when the mass is larger or of the same order than the temperature. Here, the mass term comes from the coth term and is √ 2|eB|. Then, the condition T |2eB| 1/2 is enough to make this approximation.
For the case of the pure magnetic contribution to the condensates, the procedure is a bit more complicate. We first change the proper-time integration variable time with one scaled with the magnetic field,s = |eB|s. Then, we sepa-rate the integral, one with 0 <s < 1 and the other with 1 <s < ∞. In the smalls integral, we can approximate tanh(s) ≈s and sech 2 (s) ≈ 1. In the highs integral, we can set tahh(s) ≈ 1 and sech 2 (s) ≈ 4e −2s . Now, the integrals in proper-time can be readily done, with the tricks of the momentum integral in Eq. (A.4), average with respect to the space as Eq. (43) and expand in terms of |eB| −1 .