Perfect GMR effect in gapped graphene-based ferromagnetic–normal–ferromagnetic junctions
Karbaschi Hossein†, Reza Rashedi Gholam
Department of Physics, University of Isfahan, Hezar Jerib Ave., Isfahan 81746-73441, Iran

Corresponding author. E-mail: h.karbaschi@sci.ui.ac.ir

Abstract

We investigate the quantum transport property in gapped graphene-based ferromagnetic/normal/ferromagnetic (FG/NG/FG) junctions by using the Dirac–Bogoliubov–de Gennes equation. The graphene is fabricated on SiC and BN substrates separately, so carriers in FG/NG/FG structures are considered as massive relativistic particles. Transmission probability, charge, and spin conductances are studied as a function of exchange energy of ferromagnets ( h), size of graphene gap, and thickness of normal graphene region ( L) respectively. Using the experimental values of Fermi energy in the normal graphene part ( EFN ∼ 400 meV) and energy gap in graphene (260 meV for SiC and 50 meV for BN substrate), it is shown that this structure can be used for both spin-up and spin-down polarized current. The latter case has different behavior of gapped FG/NG/FG from that of gapless FG/NG/FG structures. Also perfect charge giant magnetoresistance is observed in a range of.

Keyword: 73.63.−b; 72.25.−b; 72.80.Vp; graphene; spin-polarized current; giant magnetoresistance; nanoscale structures
1. Introduction

As is well known, graphene is a two-dimensional monolayer of graphite.[1] Its unique properties are indicated to be an appropriate candidate for constructing electronic and spintronic devices.[2, 3] The carriers in the graphene fabricated on SiO2 are massless relativistic particles with a linear dispersion relation near the K point (E (k) = ℏ kvF) where vF (∼ 106 m/s) is the Fermi velocity (Fig. 1(a)).

Recent studies show that despite the extremely high carrier mobility of gapless graphene, it is necessary to open a gap between the conduction and valence bands of graphene in order to realize some device functions, [46] so opening and tailoring a band gap in graphene is probably one of the most important and urgent research topics in graphene research currently.[79]

Since the electronic properties are strongly related to the delocalized electrons in graphene, obviously any modification will influence the electronic properties of the system. To gain the two-dimensional gapped graphene, the possible experimental approaches include the epitaxial growth of a graphene monolayer on an SiC or BN substrate. It is estimated that the graphene sample opens a band gap of about 50 meV on a boron– nitride substrate, [6, 10] and it opens about 260 meV on a silicon– carbide substrate.[11] The carriers in gapped graphene are fermions like massive Dirac particles (see Fig. 1(b)). Unique features of graphene motivate researchers to study electronic and spintronic transport properties of graphene-based systems[1238] as well as photonics and plasmonics.[3942]

Fig. 1. Energy spectra of quasiparticles in (a) gapless and (b) gapped graphene.

In this paper, we consider clean gapped graphene-based ferromagnetic/normal/ferromagnetic (FG/NG/FG) junctions where the graphene is fabricated on the SiC and BN substrates separately. The ferromagnetic state can be induced by depositing a ferromagnetic insulator on top of graphene.[43] The magnitudes of the graphene energy gap for SiC and BN substrates are 260 meV and 50 meV, respectively and the Fermi energy in the normal graphene region (EFN) is 400 meV.

2. Governing equations

Our model is composed of two similar ferromagnetic layers separated by a normal graphene of thickness L (Fig. 2). Direction of the magnetization of the left ferromagnet is fixed to be upward but the right ferromagnet has magnetization parallel (up) or antiparallel (down) to the magnetization direction of the left ferromagnet. The system is considered to be in the xy-plane. Since the considered graphene nanoribbon is short and wide (i.e., W/L ≫ 1, where L is the thickness of the normal graphene region and W is the width of the graphene sheet), microscopic details of edge and chirality are both neglected.[44, 45]

Fig. 2. Schematic illustration of a gapped FG/NG/FG junction.

The wave functions describing the quasiparticles are the solutions of the massive Dirac– Bogoliubov– de Gennes equation defined as

where H0 = − iℏ vF(σ xx ± σ yy) − EFN, with σ x, σ y, and σ z being 2 × 2 Pauli matrices; ε σ is the energy of quasiparticles related to the Fermi level; EFN is the Fermi energy of normal graphene; h is the magnetization exchange energy; vF denotes the Fermi velocity of quasiparticles; σ = ± 1 for up and down spins of quasiparticles; . Both uσ and (electron and hole-like quasiparticle wave functions) have two components of sublattices in the hexagonal lattice of graphene. The solution to Eq. (1) can be clearly seen in different regions. In the ferromagnetic region, wave functions are given by

While for the normal graphene region

Here kF, n, σ and kN, σ are the components of electron-like wave vector perpendicular to the interface in ferromagnetic and normal graphene regions, respectively. In our calculations, we have in the two-dimensional plane of graphene. Also, α 1, σ are injection angles of electron-like quasiparticles in the first ferromagnetic region, with respect to the axis normal to the interface. α F, σ and β as the incident angles in the second ferromagnetic and normal graphene regions, are calculated using the conservation of momentum in the y direction.

The refection coefficient r and the transmission coefficient t are calculated by applying the boundary conditions at the interface between ferromagnetic and normal regions. The interfaces between ferromagnets and normal graphene are perpendicular to the x axis at x = 0 and x = L, that is

for x = 0, and

for x = L. Using the boundary conditions (6) and (7), the transmission coefficient for the antiparallel case is obtained as follows:

For the parallel case we only need to set α F, σ = α 1, σ .

Using the approach of Refs. [46]– [48], the differential electronic conductance at zero temperature is

For differential spin conductance at zero temperature, we have

where G0 = 2e2/hNσ (ϵ V) is spin-dependent normal conductance, ρ σ = + 1 (− 1) for σ = ↑ (↓ ), and the density of states is with W being the width of the graphene-based FNF junction.

3. Numerical results and discussion

Figure 3 demonstrates the transmission probability (τ = | t| 2) of quasiparticles in terms of the exchange energy of ferromagnets in parallel and antiparallel structures for a gapless case and three gapped cases.

Fig. 3. Transmission probabilities of electrons versus h/EF of ferromagnets, here, the thickness of the graphene-based normal region is kL= 10. (a) Gapless graphene, (b) gapped graphene fabricated on BN substrate, , (c) gapped graphene fabricated on SiC substrate, , (d) gapped graphene with arbitrary gap size, .

It is observed that for h/EF > 3 the transmission probability is zero for all cases. Also there is an energy interval around the Fermi energy with a size of graphene gap of in which the transmission probabilities of spin-down quasiparticles are negligible in both parallel and antiparallel structures. However, the considered gapless structure can be used for only spin-down polarized current. It is clear that in gapped cases both spin-up and spin-down polarized currents may be produced. Since the transmission probability for the exchange energies is nonzero only for spin-up quasiparticles in parallel setup and in the cases of a large graphene gap, we can use the parallel setup to obtain spin-up polarized current. Like the gapless case, for the exchange energies the spin-down polarized current passes through the antiparallel structure of gapped setups.

Figures 4 and 5 show the charge conductances of parallel and antiparallel structures versus h/EF of ferromagnets, respectively.

Fig. 4. Normalized charge conductances of parallel structure versus h/EF of ferromagnets. Here, the thickness of the graphene-based normal region is kL= 10.

Fig. 5. Normalized charge conductances of antiparallel structure versus h/EF of ferromagnets. Here, the thickness of graphene-based normal region is kL= 10.

Increasing the strength of ferromagnet h, the potential well will be deeper and charge transport suppresses. It can be seen that the conductance diagrams of gapped parallel structures in terms of the exchange-energy have a special behavior at the energies close to 2EF. For the antiparallel structures, the special behavior at energy close to EF can be seen for both gapless and gapped structures. Also in the range of , which is the magnitude of the graphene energy gap, the charge conductances of all gapped structures are zero. Since a giant magnetoresistance (GMR) effect is observed as a significant change in the electrical resistance depending on whether the magnetizations of ferromagnetic layers are in a parallel or an antiparallel alignment, the latter case leads to perfect charge gaint magnetoresistance (GMR) in this energy range for gapped graphene-based FG/NG/FG junctions as shown in Fig. 6.

Fig. 6. Charge GMRs versus h/EF of ferromagnets. Here, the thickness of graphene-based normal region is kL= 10 .

Figures 7 and 8 show spin conductances of parallel and antiparallel structures versus h/EF of ferromagnets. The spin conductance of the gapped parallel structure has different behaviors in energies less than from those in energies greater than . The difference is due to the distinct behavior in the number of available states which decreases with increasing the exchange energy for exchange energies , but for exchange energies with increasing the exchange energy, the number of available states increases. Also it is clear that the spin conductances of gapped parallel structures have considerable values in exchange energies and the size of this interval is equal to graphene energy gap.

Fig. 7. Normalized spin conductances of parallel structure versus h/EF of ferromagnets. Here, the thickness of the graphene-based normal region is kL= 10.

Fig. 8. Normalized spin conductances of antiparallel structure versus h/EF of ferromagnets. Here, the thickness of the graphene-based normal region is kL= 10.

Unlike the parallel structure case, the only effect of increasing the graphene energy gap is to reduce the spin conductance of the antiparallel structure, which has significant magnitude in energy around h ∼ 2EF.

Since parallel and antiparallel spin conductances have significant values in different ranges of energies (except EF < h < 2EF), the spin GMR is almost 1 in most exchange energies as illustrated in Fig. 9.

Fig. 9. Spin GMRs versus h/EF of ferromagnets.

Our calculations show oscillating dependences of transmission probability, charge, and spin conductance on the thickness of the normal graphene region for all exchange energies except (Figs. 10 and 11). This is because of the interference between incident and reflected quasiparticles from the two different FN interfaces inside the normal region.[4952] In the range of , there is a damping factor for all quasiparticles except spin-up quasiparticles in the parallel setup.

Fig. 10. Transmission probabilities of electrons versus thickness of the normal region, h = 4EF.

Fig. 11. Charge and spin conductance of electrons versus thickness of the normal region, h = 4EF.

4. Summery and conclusion

We studied the electronic and spintronic conductances in parallel and antiparallel clean gapped graphene-based FG/NG/FG junctions by taking into account the effect of the exchange energy of magnetic layers and thickness of the normal graphene region.

In our calculations, the experimental data given in Refs. [6], [10], and [11] are used. For gapped graphene grown on SiC and BN substrates, the values of energy gap are 260 meV and 50 meV, respectively and the Fermi energy of the normal graphene region (EFN) is 400 meV.

It is illustrated that despite gapless FG/NG/FG structures, which can be used to acquire spin-up (or spin-down) polarized current (see Ref. [16]), the gapped FG/NG/FG structures can be used for both spin-up and spin-down polarized currents. It means that they can be used as a spin-polarized current switch. Also in the gapped FG/NG/FG structure, perfect charge GMR is observed in the range of while in the case of gapless FG/NG/FG structure, perfect charge GMR is obtained only in h = EF. The transmission property in gapped structure has an oscillating dependence on the thickness of the normal graphene layer, which is similar to that in the gapless FG/NG/FG structure.

Reference
1 Novoselov K S, Geim A K, Morozov S V, Jiang D, Zhang Y, Dubonos S V, Grigorieva I V and Firsov A A 2004 Science 306 666 DOI:10.1126/science.1102896 [Cited within:1]
2 Heersche H B, Herrero P J, Oostinga J B, Vand ersypen L M K and Morpurgo A F 2007 Nature 446 56 DOI:10.1038/nature05555 [Cited within:1]
3 Schedin F, Geim A K, Morozov S V, Hill E W, Blake P, Katsnelson M I and Novoselov K S 2007 Nat. Mater. 6 652 DOI:10.1038/nmat1967 [Cited within:1] [JCR: 35.749]
4 Novoselov K 2007 Nat. Mater. 6 720 DOI:10.1038/nmat2006 [Cited within:1] [JCR: 35.749]
5 Harada N, Ohfuti M and Awano Y 2008 Appl. Phys. Express 1 024002 DOI:10.1143/APEX.1.024002 [Cited within:1] [JCR: 2.731]
6 Semenoff G W, Semenoff V and Zhou F 2008 Phys. Rev. Lett. 101 087204 DOI:10.1103/PhysRevLett.101.087204 [Cited within:3] [JCR: 7.943]
7 Zhu W, Wang Z, Shi Q, Szeto K Y, Chen J and Hou J G 2009 Phys. Rev. B 79 155430 DOI:10.1103/PhysRevB.79.155430 [Cited within:1]
8 Jiang L W, Zheng Y S, Li H D and Shen H H 2010 Nanotechnology 21 145703 DOI:10.1088/0957-4484/21/14/145703 [Cited within:1] [JCR: 3.842]
9 Mousavi H 2013 Physica B 414 78 DOI:10.1016/j.physb.2013.01.025 [Cited within:1] [JCR: 1.327]
10 Giovannetti G, Khomyakov P A, Brocks G, Kelly P J and van den Brink J 2007 Phys. Rev. B 76 073103 DOI:10.1103/PhysRevB.76.073103 [Cited within:2]
11 Zhou S Y, Gweon G H, Fedorov A V, First P N, de Heer W A, Lee D H, Guinea F, Castro Neto A H and Lanzara A 2007 Nat. Mater. 6 770 DOI:10.1038/nmat2003 [Cited within:2] [JCR: 35.749]
12 Soodchomshom B, Tang I M and Hoonsawat R 2009 Phys. Lett. A 373 [Cited within:1] [JCR: 1.11]
13 Kotov V N, Pereire V M and Uchoa B 2008 Phys. Rev. B 78 075433 DOI:10.1103/PhysRevB.78.075433 [Cited within:1]
14 Zareyan M, Mohammadpour H and Moghaddam A G 2008 Phys. Rev. B 78 193406 DOI:10.1103/PhysRevB.78.193406 [Cited within:1]
15 Linder J, Zareyan M and Sudbø A 2009 Phys. Rev. B 80 014513 DOI:10.1103/PhysRevB.80.014513 [Cited within:1]
16 Mir Mojarabian F Rashedi G 2011 Physica E 44 647 DOI:10.1016/j.physe.2011.10.019 [Cited within:1] [JCR: 1.522]
17 Salehi M and Rashedi G 2010 Physica C 470 703 DOI:10.1016/j.physc.2010.07.003 [Cited within:1] [JCR: 0.718]
18 Gao W Z, Sun L and Zheng Y S 2007 Chin. Phys. Lett. 24 1693 [Cited within:1] [JCR: 0.811] [CJCR: 0.4541]
19 Ma Y N, Ren J F, Zhang Y B, Liu D S and Xie S J 2007 Chin. Phys. Lett. 24 1697 [Cited within:1] [JCR: 0.811] [CJCR: 0.4541]
20 Wang T X, Wei H X, Han X F, Langford R M, Thornton M, Bari M A and Coey J M D 2006 Chin. Phys. Lett. 23 2852 [Cited within:1] [JCR: 0.811] [CJCR: 0.4541]
21 Zhang G B, Wang S J and Li L 2006 Chin. Phys. Lett. 23 1570 [Cited within:1] [JCR: 0.811] [CJCR: 0.4541]
22 Yuan S L, Liu L, Zhong Q H, Cao H, Xiao X, Chen W, Miao J H, Niu L Y, Zhang G H, Xia Z C and Liu S 2004 Chin. Phys. Lett. 21 1352 [Cited within:1] [JCR: 0.811] [CJCR: 0.4541]
23 Xiu M X, Ren J F, Wang Y M, Yuan X B and Hu G C 2010 Acta Phys. Sin. 59 8856(in Chinese) [Cited within:1] [JCR: 1.016] [CJCR: 1.691]
24 Qiao S Z, Zhao J Q, Jia Z F, Zhang N Y, Wang F X, Fu G and Ji Y J 2010 Acta Phys. Sin. 59 8856(in Chinese) [Cited within:1] [JCR: 1.016] [CJCR: 1.691]
25 Li P, Jin C, Mi W B and Bai H L 2013 Chin. Phys. B 22 047505 DOI:10.1088/1674-1056/22/4/047505 [Cited within:1] [JCR: 1.148] [CJCR: 1.2429]
26 Yuan S P, Shen C, Zheng H Z, Liu Q, Wang L G, Meng K K and Zhao J H 2013 Chin. Phys. B 22 047202 DOI:10.1088/1674-1056/22/4/047202 [Cited within:1] [JCR: 1.148] [CJCR: 1.2429]
27 Wang Y M, Ren J F, Yuan X B, Dou Z T and Hu G C 2012 Chin. Phys. B 21 108508 DOI:10.1088/1674-1056/21/10/108508 [Cited within:1] [JCR: 1.148] [CJCR: 1.2429]
28 Wang J and Zhang L 2011 Chin. Phys. B 20 127203 DOI:10.1088/1674-1056/20/12/127203 [Cited within:1] [JCR: 1.148] [CJCR: 1.2429]
29 Duan Z G, Wang H Y, Liao W H and Zhou G H 2010 Chin. Phys. B 19 037301 DOI:10.1088/1674-1056/19/3/037301 [Cited within:1] [JCR: 1.148] [CJCR: 1.2429]
30 Guo Y, Guo W and Chen C 2008 Appl. Phys. Lett. 92 243101 DOI:10.1063/1.2943414 [Cited within:1] [JCR: 3.794]
31 Behpour M, Meshki M and Masoum S 2013 J. Nanostructures 3 243 DOI:10.1016/j.ab.2014.11.011 [Cited within:1]
32 Guo Y, Guo W and Chen C 2009 Phys. Rev. B 80 085424 DOI:10.1103/PhysRevB.80.085424 [Cited within:1]
33 Guo Y, Zhang Z and Guo W 2010 J. Phys. Chem. C 114 14729 [Cited within:1]
34 Ghoreishi F S, Ahmadi V and Samadpour M 2013 J. Nanostructures 3 453 [Cited within:1]
35 Iranmanesh A and Alizadeh Y 2011 Int. J. Nanosci. Nanotechnol. 7 28 [Cited within:1]
36 Guo Y and Guo W 2012 J. Phys. Chem. C 117 692 [Cited within:1] [JCR: 4.814]
37 Monemzadeh M, Taie M and Khoshnevisan B 2013 J. Nanostructures 3 315 [Cited within:1]
38 Guo Y and Guo W 2012 J. Appl. Phys. 111 074317 DOI:10.1063/1.3702877 [Cited within:1] [JCR: 0.71]
39 Fang Z, Liu Z, Wang Y, Ajayan P M, Nordland er P and Halas N 2012 Nano Lett. 12 3808 DOI:10.1021/nl301774e [Cited within:1] [JCR: 13.025]
40 Fang Z, Wang Y, Liu Z, Schlather A, Ajayan P M, Koppens F H L, Nordland er P and Halas N 2012 ACS Nano 6 10222 DOI:10.1021/nn304028b [Cited within:1] [JCR: 12.062]
41 Fang Z, Thongrattanasiri S, Schlather A, Liu Z, Ma L, Wang Y, Ajayan P M, Nordland er P, Halas N and de Abajo F J G 2013 ACS Nano 7 2388 DOI:10.1021/nn3055835 [Cited within:1] [JCR: 12.062]
42 Fang Z, Wang Y, Schlather A, Liu Z, Ajayan P M, de Abajo F J G, Nordland er P, Zhu X and Halas N 2014 Nano. Lett. 14 299 DOI:10.1021/nl404042h [Cited within:1] [JCR: 13.025]
43 Peres N M R, Araújo M A N and Bozi D 2004 Phys. Rev. B 70 195122 DOI:10.1103/PhysRevB.70.195122 [Cited within:1]
44 Katsnelson M I 2006 Eur. Phys. J. B 51 157 DOI:10.1140/epjb/e2006-00203-1 [Cited within:1] [JCR: 1.282]
45 Tworzydlo J, Trauzettel B, Titov M, Rycerz A and Beenakker C W J 2006 Phys. Rev. Lett. 96 246802 [Cited within:1] [JCR: 7.943]
46 Berry M V and Mondragon R J 1987 Proc. Roy. Soc. Lond. A 412 53 DOI:10.1098/rspa.1987.0080 [Cited within:1]
47 Yokoyama T, Linder J and Sudbø A 2008 Phys. Rev. B 77 132503 DOI:10.1103/PhysRevB.77.132503 [Cited within:1]
48 Bozovic M and Radovic Z 2002 Phys. Rev. B 66 134524 DOI:10.1103/PhysRevB.66.134524 [Cited within:1]
49 Annunziata G, Enoksen H, Linder J, Cuoco M, Noce C and Sudbø A 2011 Phys. Rev. B 83 144520 DOI:10.1103/PhysRevB.83.144520 [Cited within:1]
50 Soodchomshom B, Tang I M and Hoonsawat R 2009 Physica E 41 1475 DOI:10.1016/j.physe.2009.04.017 [Cited within:1] [JCR: 1.522]
51 Soodchomshom B, Tang I M and Hoonsawat R 2008 Phys. Lett. A 372 5054 DOI:10.1016/j.physleta.2008.05.051 [Cited within:1] [JCR: 1.11]
52 Zheng Z M, Qi Y N, Xing D Y and Dong J M 1999 Phys. Rev. B 59 14505 DOI:10.1103/PhysRevB.59.14505 [Cited within:1]