# Spin-dependent Seebeck effect and spin caloritronics in magnetic graphene

## Abstract

We investigate the spin-dependent thermoelectric effects in magnetic graphene in both diffusive and ballistic regimes. Employing the Boltzmann and Landauer formalisms we calculate the spin and charge Seebeck coefficients (thermopower) in magnetic graphene varying the spin splitting, temperature, and doping of the junction. It is found that while in normal graphene the temperature gradient drive a charge current, in the case of magnetic graphene a significant spin current is also established. In particular we show that in the undoped magnetic graphene in which different spin carriers belong to conduction and valence bands, a pure spin current is driven by the temperature gradient. In addition it is revealed that profound thermoelectric effects can be achieved at intermediate easily accessible temperatures when the thermal energy is comparable with Fermi energy . By further investigation of the spin-dependent Seebeck effect and a significantly large figure of merit for spin thermopower , we suggest magnetic graphene as a promising material for spin-caloritronics studies and applications.

###### pacs:

72.80.Vp, 73.50.Lw, 85.75.-d, 72.25.-b[1] [1] \vekk \vekr [1] [1] [1] H

## I introduction

Thermoelectric effects, although, known for almost two centuries, have received great attention in recent years due to their crucial relevance in meso- and nanoscopic systems Giazotto *et al.* (2006); Dubi and Di Ventra (2011).
Not only can the studies be helpful technologically in managing the generated heat in nanoelectronic devices, however; investigations about thermoelectric effects in mesoscopic regimes are of fundamental interest for condensed matter physicists Giazotto *et al.* (2006).
Starting in the late 1980’s the field of spintronics emerged which focuses on the spin-dependent transport and its coupling to that of charge Zutic *et al.* (2004); Wolf *et al.* (2001); Awschalom and Flatte (2007). Along with the fast growing interest in this field, the pioneering work of Johnson and Silsbee showed that, in spintronic and magnetic systems, heat currents can couple to spin currents as well as charge currents Johnson and Silsbee (1987). In recent years some successive unexpected experimental observations of spin Seebeck effects Uchida *et al.* (2008, 2010); Jaworski *et al.* (2010); Kirihara *et al.* (2012) have attracted a great deal of attention in investigating the thermoelectric and spintronic effects in combination with each other, which has lead to the introduction of a new research field, *spin caloritronics* Bauer *et al.* (2012); Sinova (2010). Besides many promising applications, some fundamental questions have arisen in this field, particularly about the origin of the spin Seebeck effect in different types of materials varying from metals to insulators.

Graphene, as a leading material among recently synthesized two-dimensional atomic monolayers, has received a tremendous amount of interest mostly due to its peculiar electronic structure described by the massless Dirac model Castro Neto *et al.* (2009); Geim and Novoselov (2007). A large number of possible applications in electronics, optics, nanoscale resonators, and even chemistry were suggested and implemented immediately after its discovery a decade ago. One of the main lines of investigation in graphene from the very beginning has been the electronic transport in a variety of regimes from ballistic to diffusive, and also in the extreme regimes of low density or high magnetic fields Novoselov *et al.* (2005); Zhang *et al.* (2005). The experimental observation of
linear dependence of the conductivity on the carrier density initiated a debate in the theoretical community which guided them to include long-ranged charged scatterers for an adequate description of electron transport (Ref. Das Sarma *et al.*, 2011 and references therein provide a thorough review on this topic). Intriguingly the thermoelectric properties of graphene have been also investigated both theoretically and experimentally with special focus on the neutrality or Dirac point Zuev *et al.* (2009); Wei *et al.* (2009); Hwang *et al.* (2009); Xu *et al.* (2014); Alomar and Sánchez (2014).
One of the key findings has been the
sign change of the thermoelectric power across the charge
neutrality point when the carrier type switches from electron to hole, accompanied by the divergent behavior of the Seebeck coefficient Hwang *et al.* (2009).

Besides many other promising applications recently graphene has been suggested for in spintronics devices in particular due to the long spin relaxation
lengths up to a few microns Tombros *et al.* (2007, 2008). Pioneering works of Tombros *et al.* have verified the effective spin injection into graphene via nonlocal magnetoresistance measurements. In addition it has been suggested that spin qubits based on graphene can be used as building blocks for quantum computing Trauzettel *et al.* (2007).
Interestingly a variety of methods have been
suggested to create magnetic graphene,
besides some theoretical predictions about intrinsic ferromagnetism in it Peres *et al.* (2005); Son *et al.* (2006). In practice one can use
an insulating ferromagnetic substrate or, alternatively, add a magnetic material or magnetic impurities on top of the graphene sheet to induce spin imbalance (for a review on magnetism in graphene, see Ref. Yazyev, 2010). In addition very recently the proximity-induced ferromagnetism in graphene/YIG heterostructure has been revealed which indicates a large exchange interaction Wang *et al.* (2015). In contrast to common magnetic materials due
to the gapless excitation spectrum of graphene, and the fine
tunability of its chemical potential , the spin-splitting energy between the up- and down-spin carriers can be even comparable with . So there exists a regime in which majority
and minority spins belong to different bands, conduction and valence. We have called this phase as spin-chiral due to the coupling of the real spin and the chirality and already some of its transport characteristics have been explored Moghaddam and Zareyan (2010); Zareyan *et al.* (2008); Moghaddam and Zareyan (2008); Rameshti *et al.* (2013); Zare Rameshti and Zareyan (2013). We should note that based on Zeeman splitting such
a spin-chiral graphene has been experimentally realized which shows the spin Hall effect without spin-orbit interaction Abanin *et al.* (2011).

In this paper we investigate the combination of charge, heat, and spin transport in graphene in the context of spin caloritronics and spin-dependent thermoelectric phenomena. We consider a magnetic graphene sheet in both ballistic and diffusive regimes when thermal gradients and bias voltages are applied to it. Employing the Landauer-Büttiker scattering method and Boltzmann transport equation for the two regimes we obtain the spin-dependent Seebeck and Peltier coefficients and the spin-dependent figure of merit which is a measure of thermoelectric efficiency. Our findings show that while in the absence of exchange splitting the temperature gradients drive only a charge current, in the case of magnetic graphene a spin current is also established which can be very large in comparison with charge current. In fact our key finding is for the case of undoped spin-chiral magnetic graphene in which different spin carriers are electrons and holes having the same density (see left panel of Fig. 1): a pure spin current (without charge current) is driven by the imposed temperature gradient. An explanation of this effect can be provided noting that the temperature gradient in spin-chiral graphene drives electrons from up-spin subband and holes from down-spin subband. Remembering the fact that holes carry opposite spin and charge of the corresponding electron, both types of carriers (electron and holes) carry the same spin but opposite charges which leads to the pure spin current. However for weakly magnetized and doped graphene when both spins belongs to the conduction or valence bands (right panel in Fig. 1), both spin and charge currents exist while the second dominates the first one. Adding the facts that spin relaxation is very weak in graphene and its electronic properties can be easily tuned, these results show that magnetic graphene could be promising for the spin-caloritronic applications rather than common magnetic metals.

## Ii Theoretical model and basic formalism

In order to study the spin and charge thermoelectric properties of magnetic graphene we work in the linear response regime where the relation between driving forces and the resulting generalized currents are linear. In a thermoelectric device the driving forces are temperature gradient , external electric field , and density gradients where the two last ones can be combined in an effective field with indicating the chemical potential. These fields can drive carriers leading to the charge and heat currents which can differ for different spin channels in general. It is believed that in graphene spin relaxation time is so long that in mesoscopic samples we can treat the two spin channels almost independently. As a result the charge and heat currents carried by spin electrons which can be denoted by and , respectively, are linearly related to the effective electric field and the temperature gradient as

(1) |

By definition the first component is the conductance and the two off-diagonal components are thermoelectric coefficients which are related to each other with Onsager relation (). The last component contributes in the -electron thermal conductivity defined by

(2) |

In the upcoming subsections we will give the explicit relations for the matrix elements in the diffusive and ballistic regimes.

The spin-dependent Seebeck and the Peltier coefficients which, for each spin channel, describe the voltage generation due to the temperature gradient and heat current induction due to the charge current, respectively, are then given by

(3) |

From these relations one can define the charge and spin Seebeck and the Peltier coefficients which are as follows,

(4) |

The ability of a material to efficiently produce thermoelectric power is usually described by a dimensionless figure of merit denoted by . In spin caloritronics we can generalize this concept for the resulting charge and spin currents due to the temperature gradients separately. So the charge and spin figures of merit for a magnetic system can be defined versus Seebeck coefficients as

(5) |

Here, () denotes the charge (spin) conductivity of the system and the the electron thermal conductivity is given by . We concentrate on low enough temperatures where only electrons contribute effectively in thermal transport. In Sec. III.3 based on some estimations, we will discuss how this assumption is verified and what are its limitations.

In the remaining of this section, the theoretical frameworks to calculate the spin and charge thermoelectric coefficients in the diffusive and ballistic transport regimes employing Boltzmann and Landauer formalisms, respectively, will be presented.

### ii.1 Diffusive regime: Boltzmann transport

In this section we give the semiclassical Boltzmann equation to establish the transport coefficients in the diffusive regime. In particular, we take into account two important cases of short-range (SR) impurities with Dirac delta potentials and long-range (LR) Coulomb impurities in our investigation. The spin-dependent thermoelectric properties due to the presence of both electric fields and temperature gradient will be found in the scheme of relaxation time approximation.

In the diffusive regime, the transport coefficients can be calculated from the following general expression for electron current and energy flux density,

(6) |

in which is the semiclassical velocity of the spin carriers. The nonequilibrium distribution function describes the evolution of the electron distribution in the presence of external perturbations. The Boltzmann formalism in the relaxation time approximation scheme and in the linear response leads to the following expression for disturbed function as

(7) |

with denoting the spin-dependent relaxation time and the equilibrium-state Fermi-Dirac distribution at temperature . We note that since graphene has an isotropic dispersion relation, the relaxation time depends only on the energy of electrons .

By invoking the above expression for into Eq. (6) for spin-dependent charge and heat currents,
the matrix coefficients can be expressed in terms of some spin-dependent *kinetic coefficients* as the following,

(8) |

All of the coefficients obey the relation

(9) |

with spin-dependent conductivity given by

(10) | |||||

with spin-dependent density of states (SDOS) . The formalism introduced so far is general for any isotropic magnetic material and all of the thermoelectric properties described by can be found as functions of spin-dependent relaxation time and SDOS . Now we switch to the case of our investigation, magnetic graphene.

A monolayer graphene sheet at the presence of induced spin splitting can be described by a low-energy Dirac Hamiltonian of the form

(11) |

with Fermi velocity , momentum , and exchange splitting . Pauli matrices and () operate on the subspaces of pseudospin (originating from two different trigonal sublattices and of the graphene’s hexagonal structure) and real spin, respectively. The spin-dependent band dispersion then follows with corresponding to the two spin directions and indicating the chirality of states. Since the velocity of carriers in graphene is constant and the SDOS is given by the spin-dependent Boltzmann conductivity of magnetic graphene takes the Drude form as

(12) |

Early investigations of quantum transport in graphene at the presence of impurities have shown that the relaxation time for the SR impurities varies inversely with the DOS as while the LR Coulomb impurities result in Das Sarma *et al.* (2011); Nomura and MacDonald (2007). Therefore in magnetic graphene, the conductivity becomes constant when only SR scatterers are present while the conductivity caused by scattering from LR impurities is proportional to the square of density of states as . In the next section we will use the relations of conductivities to obtain the Seebeck coefficients and corresponding figures of merit.

### ii.2 Ballistic regime: Landauer-Büttiker formula

Within the Landauer-Büttiker approach the electric and thermal currents carried by electrons with spin are obtained from the transmission probabilities integrated over the energy and the angle ,

(13) | |||||

where is the sample width and and are the Fermi-Dirac distribution functions of the left and right electronic leads, respectively. Assuming the linear response regime, we can expand the difference of Fermi-Dirac functions in the above formula up to linear terms in a small bias voltage and temperature difference between two reservoirs. This results in a relation very similar to Eq. (1) with thermoelectric conductances related to the kinetic coefficients according to Eq. (8). These coefficients for the ballistic transport regime are obtained after some straightforward algebra,

(14) | |||||

in which denotes the ideal conductance of the junction. It is easy to note that this relation is in complete accordance with Eq. (9) in which only the conductivity is replaced with the conductance of the ballistic system given by the Landauer-Büttiker formula. So the only thing we need is to determine the transmission probabilities .

In the ballistic regime we consider a magnetic graphene sheet between two electrodes in which their chemical potentials despite a very small bias are kept at while the central region’s doping can be varied. Diagonalizing the Hamiltonian (11) results in eigenstates , where indicates the band index and specifies the propagation direction. The total scattering wave function, corresponding to an incoming excitation coming from the left reservoir, inside each region, can be written as

(15) |

with and (for the sake of simplicity the overall factor is dropped in all components of wave functions). The transmission and reflection amplitudes are determined by matching the wavefunctions at the interfaces and and subsequently used to calculate the total transmission probability as

(16) |

We recall that there is no mechanism that couples states with opposite spin indices so no off-diagonal spin channel mixing terms will appear in our calculations.

## Iii Results and Discussion

In this section our numerical and analytical results will be presented. We focus on the charge and spin Seebeck coefficients (, ) and their corresponding figures of merit and . We can also find the Peltier coefficients and in both diffusive and ballistic regimes. However according to the so-called Thomson relation , which is originated from the symmetry properties of the coefficients demanded by Onsager reciprocity, there is no new information on the Peltier coefficients. We divide this section into two parts concentrating on the diffusive and ballistic systems, respectively.

### iii.1 Diffusive transport

It is well known that electron-hole asymmetry around the Fermi level in the band structure or transport properties is responsible for the thermoelectric effects. In fact the key role in thermoelectric effects is played by which according to Eq. (9) vanishes when is a symmetric (even) function of . In the case of graphene, at very first glance, the Dirac dispersion relation and linear energy dependence of DOS suggest a possible source of asymmetry in away from the neutrality point which can lead to thermoelectric phenomena. However as we have seen in the previous section when only short range scatterers are present the conductivities become constant. Therefore diffusive transport caused by SR impurities leaves magnetic graphene with no thermoelectric effects with vanishing charge and spin Seebeck and Peltier coefficients.

In contrast at the presence of long range Coulomb impurities, which are in fact the dominant scatterers in most graphene samples, the spin-dependent conductivities have explicit energy dependence. Invoking the quadratic energy dependence of conductivities in Eq. (9) and performing the integrations over energies we find a simple form for the Seebeck coefficient of spin- carriers,

(17) |

in which depends explicitly on temperature and exchange splitting. The method of the calculation of and the behavior of chemical potential as a function of exchange and temperature will be presented in the Appendix. As one can see from Eq. (17) two spin-dependent Seebeck coefficients reach their maximum absolute values at temperatures , respectively. In addition we see that each of the coefficients passes from zero and changes sign when the Fermi level of the corresponding spin subband lies at the Dirac point . This is similar to the well-known effect in semiconductors in which the thermopower for and types has opposite sign and based on this effect devices made of junctions are used for electronic cooling. However a big advantage in the case of graphene is provided by the fine-tunability of doping in it. So in real experimental situations one can play with , exchange splitting, and also temperature to cover a wide range of parameter space.

Now we turn the discussion to which are more feasible quantities in real experiments. Using Eq. (17) charge and spin Seebeck coefficients can be easily obtained,

(18) | |||||

(19) |

Inserting the numerically calculated in the above relations, the variation of thermopowers with temperature and spin splitting is obtained as shown in Fig. 2. As we expect at very low temperatures, , charge and spin Seebeck effects are very weak and go to zero linearly with . On the other hand at some intermediate temperatures when the thermal energy is comparable with the spin-dependent Fermi levels measured from the Dirac points () profound Seebeck effects can be observed.

When the spin splitting is small () both up- and down-spin Fermi levels lie in the conduction band and thermally activated electrons of both spins move along temperature gradient which result in a charge accumulation gradient in the opposite direction due to the negative charge of the electrons. Therefore a negative charge thermopower is obtained for (). By further increase in the exchange splitting then one of the spin subbands’ Fermi level goes to the valence band and then the holes from the spin-down subband will be thermally activated. Such excitations carry positive charge current and as a result their contribution in the charge Seebeck effect has positive sign while spin up electrons from the conduction band still have a negative contribution, which means the excitations from two spin subbands compete with each other. By further increase of the exchange the contribution of minority spin carriers from the hole band dominates and as one can clearly see from Eq. (18) at (indicated by the dashed line in Fig. 2 ) the charge thermopower changes its sign.

The spin Seebeck effect behaves in a somehow opposite way with the variation of exchange field. At low electrons carrying different spins compete with each other to result in a spin accumulation caused by the temperature gradient. As we see in the right panel of Fig. 2 for not so high temperatures the minority spins are dominant and as a result unlike a positive spin Seebeck effect is observed. But at higher temperatures becomes all negative dominated by majority up spins. Upon increasing the exchange splitting when the down spins’ Fermi level goes to the valence band both up-spin electrons and down-spin holes, which carry the same intrinsic angular momentum, accompany each other to give a strong spin signal. In fact as one can immediately see from Eq. (19) for when the sign of the spin thermopower changes and for higher exchanges becomes negative.

Now the key finding of our work is the fact that by moving along the curve , we can completely turn off the charge Seebeck effect, while a spin Seebeck effect can be observed. This is clear if we compare two plots in Fig. 2 and notice that along the line of a large negative spin Seebeck coefficient is obtained which is given by

(20) |

As a result a maximum value can be reached in the absence of corresponding charge signal.

In order to see the exchange dependence of the thermopowers more clearly, Fig. 3 shows the charge and spin Seebeck coefficients and the corresponding figures of merit as functions of for some different temperatures. We see that although charge thermopower always shows sign change at the vicinity of , spin thermopower becomes all negative at higher temperatures, irrespective of spin splitting strength . This is again a clear manifestation of the possibility of pure spin current caused by temperature gradient. In fact at higher this is easily understood from the fact that conduction band spin-up electrons and valence band spin-down holes accompany each other to give rise to a negative . At lower as we mentioned above by increasing temperature majority up spins dominate the thermoelectric effect and since they carry negative current, the spin Seebeck coefficient remains still negative. In addition as expected the spin (charge) figure of merit reaches its maximum value for some splitting above (below) the chemical potential . The figures of merit for both spin and charge Seebeck effects becomes large (of the order of 1) at some intermediate temperatures where the thermoelectric effect is very strong while the heat transport is not. On the other hand at higher temperatures the thermopowers decrease as the inverse of and subsequently the figures of merit show decline with temperature. These effects can be seen from Fig. 4 where the variations of and are shown with temperature . These results again clearly show that the strong thermoelectric effects can be seen at the intermediate temperatures when is comparable with Fermi levels measured from neutrality point. It is worth noting again that in graphene the doping can be varied easily and as a result one can even reach the regimes in which the Fermi energy is comparable to thermal energy . So unlike conventional metals with very large Fermi energy in comparison with , it is reasonable to reach the most efficient values for thermoelectric responses.

### iii.2 Ballistic transport

Now we turn to the spin-dependent thermoelectric properties of ballistic graphene. Unlike the diffusive case due to the complicated energy dependence of transmission coefficients we cannot obtain simple analytic relations. Figure 5 shows the numerically obtained results for spin and charge thermopowers in ballistic regimes. Fascinatingly the overall behavior is almost the same as the diffusive regime in the presence of long-range Coulomb impurities. In fact, comparing diffusive and ballistic results, we only see that they are only slightly different in quantitative manner. For instance the possible maxima of thermopowers and also the lines in which they vanish are different for two cases which is related to details originating from scattering mechanisms of two regimes. In other words, our results show that despite the details of scattering phenomena, the band structure and dispersion of graphene play a main role in the spin-dependent thermoelectric effects. Of course we know that in the case of diffusive transport the presence of long-range impurities is crucial for thermoelectric effects. Nevertheless when the Seebeck effect does exist, the dependence on the temperature and spin splitting is more or less universal and despite the transport regime we see the same features. It should be noted that such universal behavior is partly related to the definition of Seebeck coefficients themselves and generally we find .

In the case of ballistic devices we also investigate the effect of gate voltage in the middle region. This could be of great importance in real applications since the gate voltages can be easily tuned. Subsequently one can control the spin-dependent thermoelectric properties by changing the chemical potential of middle graphene between the two electrodes and can tune the spin-caloritronic properties. The dependence of spin and charge thermopowers and corresponding figures of merit on the exchange splitting scaled by the leads’ chemical potential at zero temperature are shown in Fig. 6 for a variety of . First of all we see that changing results in shifts in the dependencies of and which can be easily understood due to the fact that the middle magnetic graphene doping plays the main role in transport properties rather than the leads’ doping . Second it is clear that the overall amplitude of the Seebeck coefficients and figures of merit also vary by changing the gate voltages. This is related to the fact that when the chemical potential of electrodes and middle graphene are different the energy-dependent transmission coefficient changes. However unless these changes do not affect the magnitude of thermopowers since both and scale with overall transparency of the scattering region (middle graphene) and only the shift as a function of is observable.

### iii.3 On the experimental reliability

We close our discussion with commenting on the possible experimental realization of the results we find. First of all we should recall that in all our models we have ignored spin relaxation and spin-flip scattering which is verified until the device size is smaller than spin relaxation length. So in order to have a strong spin Seebeck effect and usage for spin-caloritronics application, we need devices of length which is easily accessible in current experimental devices Guimarães *et al.* (2014). On the other side an important step is to combine already existing experimental spintronic and thermoelectric setups based on graphene. This is apparently an easy task since on one side nonlocal magnetoresistance measurements are proven to be very useful to detect spin injection and spin currents in graphene-based spintronic devices. On the other hand the thermoelectric effects themselves have been already observed with significant precision Zuev *et al.* (2009); Wei *et al.* (2009) which suggests that spin-dependent thermoelectric properties, in principle, could be detected with high feasibility.

In order to reveal the relevance of our finding to
the experimental situations, in Fig. 7 the variations of
charge and spin thermopowers with chemical potential ()
for different temperatures are shown in the case of the diffusive regime. In this figure unlike
previous ones instead of scaled dimensionless parameters
we use reliable numerical values of parameters in electronvolts, kelvins, etc., and in particular the temperatures
are exactly the same as in Ref. Zuev *et al.*, 2009. In addition we assume the predicted value for exchange splitting Haugen *et al.* (2008). We see that the numerical values of thermopowers reach values on the order of a hundred which is consistent with previous experimental results for nonmagnetic graphene. In fact when the temperature is large enough in comparison with exchange splitting () the qualitative behavior of is very close to that obtained by Zuev et al. Zuev *et al.* (2009). But very interestingly close to the Dirac point () the spin Seebeck coefficient becomes very large especially for intermediate temperatures () as discussed before.

Finally, we should comment on the possible influences of phonons in our result which we have not considered. The main effect of phonons is their contribution in the thermal conductivity and the charge and spin thermopower are not affected with the presence of phonons. Therefore it is clear that the thermal conductivity of phonons can only affect the figures of merit in our results and since it does not depend on chemical potential or exchange splitting it will only increase depending on temperature. This will decrease but dependence on and will not be changed qualitatively, whatsoever.
Second, the thermal conductivity of phonons will decrease according to some power-law behavior with upon decreasing temperature Xu *et al.* (2014), while the electrons contribution () in the thermal conductivity varies linearly with at low temperatures. This can be easily seen from the exact formula which can be obtained in the diffusive regime,

(21) | |||||

with impurity concentration .
From the experimental results on phonon thermal conductivity in suspended graphene Seol *et al.* (2010), we can estimate values on the order of at temperatures , respectively. Then from Eq. (21) and assuming typical values , , we see that for thermal conductivity of electrons decreases linearly with and for instance at , we get which is the same as the phonon contribution. So we can conclude that for low temperatures the electrons dominate the thermal conductivity in graphene. At higher temperatures phonons becomes important but as we mentioned before it only results in the overall decline of predicted figures of merit, without affecting their qualitative behavior.

## Iv Conclusions

In this study we reveal that magnetic graphene could be very promising for spin-caloritronics studies and applications. Employing Boltzmann and Landauer formalisms, the spin-dependent thermoelectric properties of graphene in both diffusive and ballistic regimes are obtained. The main finding is that while in the absence of spin splitting, the temperature gradient drives a charge current in graphene, by imposing spin splitting a significant spin current is established, too. Very intriguingly when we consider an undoped magnetic graphene in which different spin carriers belong to conduction and valence bands, we will have a pure spin thermopower without charge thermopower. This pure spin current generation by temperature gradient can be achieved in the temperature and spin splitting of the order of the unpolarized state Fermi energy which is accessible in current experiments. So based on this study, we believe that besides the suggested applications of graphene for spintronic devices due to long spin relaxation, magnetic graphene can be used as a base material to investigate spin-thermoelectric phenomena.

## Acknowledgments

Authors would like to thank Gerrit Bauer for fruitful discussion. A.G.M. acknowledges warm the hospitality of Michele Govenale and Ulrich Zülicke and their useful comments on the manuscript, during his visit at Victoria University of Wellington where part of this work was done. B.Z.R. thanks the CMSP of ICTP in Trieste for hospitality and support during his visit to this institute where part of this work was done. We are also grateful to Saeed Abedinpour for his useful comments.

## Appendix A Temperature dependence of the chemical potential

In this appendix we will present the results of chemical potential variations with temperature and spin splitting . The conventional way to obtain the temperature dependence of chemical potential is to enforce the following quantity to be constant,

(22) |

with Fermi distribution function

(23) |

which is nothing but the total number of electrons in the system. However in the case of the massless Dirac model for graphene in which there is no lower band for energy the integration over energy diverges. However we can easily overcome this difficulty by subtracting the infinite number of negative energy states. So we define the excess number of charge carriers instead of all electrons,

(24) |

which must be a constant irrespective of temperature variations. Inserting the density of states of magnetic graphene and equating the finite-temperature value of the above expression with its zero-temperature correspondence, we will have

(25) |

### References

- F. Giazotto, T. T. Heikkilä, A. Luukanen, A. M. Savin, and J. P. Pekola, Rev. Mod. Phys. 78, 217 (2006).
- Y. Dubi and M. Di Ventra, Rev. Mod. Phys. 83, 131 (2011).
- I. Zutic, J. Fabian, and S. D. Sarma, Reviews of Modern Physics 76, 323410 (2004).
- S. A. Wolf, D. D. Awschalom, R. A. Buhrman, J. M. Daughton, S. von Molnár, M. L. Roukes, A. Y. Chtchelkanova, and D. M. Treger, Science 294, 1488 (2001).
- D. D. Awschalom and M. E. Flatte, Nat. Phys. 3, 153 (2007).
- M. Johnson and R. H. Silsbee, Phys. Rev. B 35, 4959 (1987).
- K. Uchida, S. Takahashi, K. Harii, J. Ieda, W. Koshibae, K. Ando, S. Maekawa, and E. Saitoh, Nature 445, 778 (2008).
- K. Uchida, J. Xiao, H. Adachi, J. Ohe, S. Takahashi, J. Ieda, T. Ota, Y. Kajiwara, H. Umezawa, H. Kawai, G. E. W. Bauer, S. Maekawa, and Saitoh, Nature Mater. 9, 894 (2010).
- C. M. Jaworski, J. Yang, S. Mack, D. D. Awschalom, J. P. Heremans, and R. C. Myers, Nature Mater. 9, 898 (2010).
- A. Kirihara, K. Uchida, Y. Kajiwara, M. Ishida, Y. Nakamura, T. Manako, E. Saitoh, and S. Yorozu, Nature Mater. 11, 686 (2012).
- G. E. W. Bauer, E. Saitoh, and B. J. van Wees, Nature Mater. 11, 391 (2012).
- J. Sinova, Nature Mater. 9, 880 (2010).
- A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov, and A. K. Geim, Rev. Mod. Phys. 81, 109 (2009).
- A. K. Geim and K. S. Novoselov, Nature Mater. 6, 183 (2007).
- K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, M. I. Katsnelson, I. V. Grigorieva, S. V. Dubonos, and A. A. Firsov, Nature 438, 197 (2005).
- Y. Zhang, Y.-W. Tan, H. L. Stormer, and P. Kim, Nature 438, 201 (2005).
- S. Das Sarma, S. Adam, E. H. Hwang, and E. Rossi, Rev. Mod. Phys. 83, 407 (2011).
- Y. M. Zuev, W. Chang, and P. Kim, Phys. Rev. Lett. 102, 096807 (2009).
- P. Wei, W. Bao, Y. Pu, C. N. Lau, and J. Shi, Phys. Rev. Lett. 102, 166808 (2009).
- E. H. Hwang, E. Rossi, and S. Das Sarma, Phys. Rev. B 80, 235415 (2009).
- Y. Xu, Z. Li, and W. Duan, Small 10, 2182 (2014).
- M. I. Alomar and D. Sánchez, Phys. Rev. B 89, 115422 (2014).
- N. Tombros, C. Jozsa, M. Popinciuc, H. T. Jonkman, and B. J. van Wees, Nature 448, 571 (2007), 0706.1948 .
- N. Tombros, S. Tanabe, A. Veligura, C. Jozsa, M. Popinciuc, H. T. Jonkman, and B. J. van Wees, Phys. Rev. Lett. 101, 046601 (2008).
- B. Trauzettel, D. V. Bulaev, D. Loss, and G. Burkard, Nat. Phys. 3, 192 (2007).
- N. M. R. Peres, F. Guinea, and A. H. Castro Neto, Phys. Rev. B 72, 174406 (2005).
- Y. W. Son, M. L. Cohen, and S. G. Louie, Nature 444, 347 (2006).
- O. V. Yazyev, Rep. Prog. Phys. 73, 056501 (2010).
- Z. Wang, C. Tang, R. Sachs, Y. Barlas, and J. Shi, Phys. Rev. Lett. 114, 016603 (2015).
- A. G. Moghaddam and M. Zareyan, Phys. Rev. Lett. 105, 146803 (2010).
- M. Zareyan, H. Mohammadpour, and A. G. Moghaddam, Phys. Rev. B 78, 193406 (2008).
- A. G. Moghaddam and M. Zareyan, Phys. Rev. B 78, 115413 (2008).
- B. Z. Rameshti, A. G. Moghaddam, S. H. Abedinpour, S. Abdizadeh, and M. Zareyan, Phys. Rev. B 88, 245424 (2013).
- B. Zare Rameshti and M. Zareyan, Appl. Phys. Lett. 103, 132409 (2013).
- D. A. Abanin, S. V. Morozov, L. A. Ponomarenko, R. V. Gorbachev, A. S. Mayorov, M. I. Katsnelson, K. Watanabe, T. Taniguchi, K. S. Novoselov, L. S. Levitov, and A. K. Geim, Science 332, 328 (2011).
- K. Nomura and A. H. MacDonald, Phys. Rev. Lett. 98, 076602 (2007).
- M. H. D. Guimarães, J. J. van den Berg, I. J. Vera-Marun, P. J. Zomer, and B. J. van Wees, Phys. Rev. B 90, 235428 (2014).
- H. Haugen, D. Huertas-Hernando, and A. Brataas, Phys. Rev. B 77, 115406 (2008).
- J. H. Seol, I. Jo, A. L. Moore, L. Lindsay, Z. H. Aitken, M. T. Pettes, X. Li, Z. Yao, R. Huang, D. Broido, N. Mingo, R. S. Ruoff, and L. Shi, Science 328, 213 (2010).