# Intrinsic Spin-Orbit Coupling in Zigzag and Armchair Graphene Nanoribbons

Authors: Ying Li; Erhu Zhang; Baihua Gong; Shengli Zhang
Journal: Journal of Nanomaterials
 (2011)

---

## Abstract

Starting from a tight-binding model, we derive the energy gaps induced by intrinsic spin-orbit
(ISO) coupling in the low-energy band structures of graphene nanoribbons. The armchair graphene
nanoribbons may be either semiconducting or metallic, depending on their widths in the absence of
ISO interactions. For the metallic ones, the gaps induced by ISO coupling decrease with increasing
ribbon widths. For the ISO interactions, we find that zigzag graphene nanoribbons with odd
chains still have no band gaps while those with even chains have gaps with a monotonic decreasing
dependence on the widths. First-principles calculations have also been carried out, verifying the
results of the tight-binding approximation. Our paper reveals that the ISO interaction of graphene
nanoribbons is governed by their geometrical parameters.


---

## Body

## 1. Introduction

Graphene, first obtained in 2004 [1], has been explored for spintronics due to its unique and intriguing properties [2]. In spintronics, one of the main interactions that could affect the spin behavior is intrinsic spin-orbit (ISO) coupling. In the last few years, a volume of research has been performed on the ISO coupling in graphene [3–8] and shows that the ISO coupling can cause a gap in the band structure. Recently, it is reported that the geometric curvature [5, 9, 10] and impurity [11] can enhance the effective strength of the ISO by orders of magnitude.When graphene is patterned into a narrow ribbon, one obtains a graphene nanoribbon (GNR), whose carriers are confined to the one-dimensional (1D) system. There are two basic shapes of GNRs, called armchair GNRs (AGNRs) and zigzag GNRs (ZGNRs). Both results by tight-binding(TB) method [12–15] and first-principles calculations [16–19] show that spectra of GNRs depend strongly on the shapes of their edges. GNRs, offering a system to investigate confinement effects, could be key elements in future electronics and spin-transport applications. In the case of GNRs with ISO coupling, on one hand, based on TB model, Zarea et al. [20–22] have studied the effects of ISO interactions on low-energy properties of AGNRs and ZGNRs. However, these works were based on the effective π-orbital TB approaches and used the model from graphene plane by imposing hard-wall boundary conditions. They did not study ISO from that of each atom. They also did not give the exact data of the gap induced by ISO coupling. On the other hand, based on first principles, few investigations of the ISO interaction on GNRs have so far been performed. Thus, the ISO coupling is still significant and challenging in GNR.In this paper, we study the ISO effects on ZGNRs and AGNRs by thesp3 TB method from a graphene ribbon model. The first-principles calculations are also performed to verify the results. Our investigation reveals how the geometrical parameters, such as width and the edge shape, affect the strength of the ISO.


## 2. Tight-Binding Model

Two basic kinds of GNRs are shown in Figure1. The AGNRs are classified by number of dimer lines while ZGNRs by number of chains. To simplify, we use N-AGNR and N-ZGNR to represent the AGNR with N dimer lines and the ZGNR with N chains, respectively. As seen in Figure 1, a GNR is described by a rectangular unit cell containing 2N carbon atoms. The 2N atoms can be classified to N subsystems. Each subsystem contains two atoms A and B. The unit vector a⃗A(a⃗Z) denotes the translational axis of an (a) AGNR (ZGNR). The nearest-neighbor translation vectors of the atoms are
(1)N⃗A={N⃗1A,N⃗2A,N⃗3A}={a(13,0),a(-123,-12),a(-123,12)},N⃗Z={N⃗1Z,N⃗2Z,N⃗3Z}={a(-12,-123),a(12,-123),a(0,13)},
where a=2.46 Å. We employ the tight-binding method within a two-center Slater-Koster approximation for nearest-neighbor hopping, and orthogonality between Wannier functions centered on different atoms is assumed. This gives the elements of the Hamiltonian matrix [4]  (2)hA,μ;A,μ′(k⃗)=hB,μ;B,μ′(k⃗)=tμδμ,μ′,hA,μ;B,μ′(k⃗)=hB,μ′;A,μ*(k⃗)=∑i=nearesteik⃗⋅N⃗itμ,μ′(N⃗i),
where μ(μ′)=s,px,py,pz, A and B represent the two distinct sites in the honeycomb lattice unit cell. tμ is the atomic energies of μ orbital, and tμ,μ′ denotes the hopping energies between the nearest-neighbor sites. In Table 1, we reproduce for completeness the relationship between tμ,μ′ and the four independent Slater-Koster parameters. To simplify, we let θ=(ssσ),  α=(spσ), β=(ppσ), and η=(ppπ). The wave vector k is parallel to x.Table 1
Two-center matrix elements for hopping betweens and p orbital along a direction specified by the unit vector (nx,ny,nz), taken from [4]. s=-8.868 eV,  p=0, (ssσ)=-6.769 eV, (spσ)=5.580 eV, (ppσ)=5.037 eV, and (ppπ)=-3.033 eV.
tsstpx,pxnx2(ppσ)+(1-nx2)(ppπ)tpptpy,pyny2(ppσ)+(1-ny2)(ppπ)ts,s(ssσ)tpz,pznz2(ppσ)+(1-nz2)(ppπ)ts,pxnx(spσ)tpx,pynxny(ppσ)-nxny(ppπ)ts,pyny(spσ)tpx,pznxnz(ppσ)-nxnz(ppπ)ts,pznz(spσ)tpy,pznynz(ppσ)-nynz(ppπ)The geometry of (a)N-AGNR and (b) N-ZGNR.
(a)(b)First, we consider theπ electrons. The 2×2 hamiltonian matrix Hπmn is introduced to describe the interaction of m and n subsystems. The 2N×2N hamiltonian matrix for N-GNR can be written as follows:(3)Hπ=(Hπ1,1Hπ1,20⋯000Hπ2,1Hπ2,2Hπ2,3⋯000⋮⋮⋮⋮⋮⋮000⋯HA,N-2HA,N-1HA,N000⋯0HπN,N-1HπN,N),
where 𝒜 denotes  πN-1.We can diagonalize the above hamiltonian and gain the eigenvalues for GNRs. The results are the same with those of [12]. The energy spectra of AGRNs for different widths can be observed in Figure 2. The system is metallic for AGNRs with widths N=3M-1, where M is an integer, and insulating otherwise. For the metallic AGNRs, the top of valence band and the bottom of conduction band are degenerate at KA=0 in the Brillouin zone (BZ). For the insulating ones, the gaps decrease with an increase in ribbon widths and tend to zero in very large N. The calculated band structures of ZGNRs are depicted in Figure 3 for four different ribbon widths. It is observed that for ZGNRs, the valence band and the conduction band are degenerate at KZ=π/a and KZ′=-π/a in the BZ.Calculated band structures of AGNRs with (a)N=4, (b) N=5, (c) N=6, and (d) N=60.
(a)(b)(c)(d)Calculated band structures of ZGNRs with (a)N=4, (b) N=5, (c) N=6, and (d) N=60.
(a)(b)(c)(d)


## 3. Intrinsic Spin-Orbit Coupling

HSO induces σ-π coupling. So, we discuss the hamiltonian including both π and σ electrons. The hamiltonian for N-GNR is an 8N×8N matrix:(4)Hπσ=(H1,1H1,20⋯000H2,1H2,2H2,3⋯000⋮⋮⋮⋮⋮⋮000⋯HN-1,N-2HN-1,N-1HN-1,N000⋯0HN,N-1HN,N).We consider the ISO interaction as an additional term toHπσ. The ISO interaction HSO of the system comes from each atomic ISO hamiltonian hSO. hSO is in the form of hSO=ξL⃗·S⃗, where ξ=6meV denotes the ISO coupling constant. L⃗ and S⃗ represent the angular momentum operator and the spin operator, respectively. The hamiltonian hSO can be expressed as follows [3–5]:(5)hSO=(0000000000iξ20000ξ20-iξ200000iξ200000ξ2iξ2000000000000ξ200-iξ20000-iξ20iξ2000ξ2-iξ200000),
which acts on the basis states |s↑〉, |px↑〉, |py↑〉, |pz↑〉, |s↓〉, |px↓〉, |py↓〉, and |pz↓〉. Then, the total hamiltonian including ISO interaction is H=Hπσ+HSO. Finally, we solve the total hamiltonian for GNRs and obtain the eigenvalues.Since it has been shown that the ISO coupling in graphene just causes a small gap [3–8], thus the spectra near the Fermi surface attract our interest. Hence, we focus on the metallic GNRs. If ISO interaction is taken into account, the band structures are lifted. In metallic AGNRs and the ZGNRs with even chains, a gap is opened in the spectra by ISO coupling. In contrast, for the ZGNRs with odd chains, conduction and valence edge bands are still degenerated at KZ=π/a. The difference between the influence on odd ZGNRs and even ZGNRs is likely due to their configuration. The odd ZGNRs are in “zigzag/antizigzag” configuration while even ZGNRs are in “zigzag/zigzag” configuration. We display the energy spectrum near Fermi surface of 5-AGNR with ISO interaction near KA point in Figure 4(a). The band gap caused at k=0 is 1.7×10-3 meV. The band structures of 5-ZGNR and 4-ZGNR with ISO interaction near KZ point are presented in Figures 4(b) and 4(c). For 4-ZGNR, the highest valence band state and the lowest conduction band state are located at k=0.9973π/a rather than k=π/a point. There is a band gap of 1×10-5 meV induced by ISO. The ISO interaction is shown to have a stronger effect in 5-AGNR than in 4-ZGNR. It is possible because that the proportion of edge atoms in unit cell atoms which is 4/10 for 5-AGNR is larger than 2/8 for 4-ZGNR. The relation between the gaps and the widths will be discussed in the next section to compare the results with the first-principles calculations data.The band structures of the GNRs with ISO interaction for (a) 5-AGNR, (b) 5-ZGNR, and (c) 4-ZGNR. The band gaps are1.7×10-3 meV, 0, and 1×10-5 meV, respectively.
(a)(b)(c)


## 4. First-Principles Calculations

In order to verify the results from TB model, we study the metallic GNRs by first-principles calculations. The relativisticab initio electronic structure calculations for graphene nanoribbons have been performed using the Vienna ab initiosimulation package (VASP) [23]. We choose the projector augmented wave (PAW) method [24] with a Perdew-Burke-Ernzerhof (PBE) generalized gradient exchange correlation potential (GGA) to describe the electron interactions [25]. A vacuum separation of 10 Å is used to eliminate interactions between graphene layers. An energy cutoff of 950 eV and Monkhorst-Pack k-point meshes are employed to deal with the many atoms in a unit cell of the GNRs for various widths. The ISO interaction is calculated considering only the spherical part of the potential inside muffin tins surrounding the carbon nuclei. The edge carbon atoms of our GNRs are passivated by hydrogen atoms. Due to the smallness of ISO-derived gaps in GNRs, we artificially increase the relativistic potential in VASP code in our calculations.First, we study the band structures of AGNRs with different widthsN=3M-1 in the range of N=5 to N=35. We have known that the AGNRs with N=3M-1 show the metallic behaviors in TB method without ISO interactions, and gaps are opened by the ISO coupling between the valence band and the conduction band. While for the first principles calculations, we find that there are always gaps in the band spectra, but the ISO coupling can reduce the gaps. The first-principles band structure near Fermi surface of 5-AGNR is shown in Figure 5. Without ISO interaction, the band gap is 0.5379 eV. When the ISO coupling is considered, the band gap becomes smaller. As Figure 6 depicts, the change of the gap ΔE decreases when N increases from 8 to 35, which is in good agreement with the results of TB method except N=5. In 5-AGNR, the gap change is smaller in TB method than that expected. In TB method, we do not consider the change of the carbon-carbon bond length on the edge. As the smallest width metallic AGNR, 5-AGNR, is affected most by the edge. So, the model may not exactly express the ISO interaction in 5-AGNR.Figure 5
(Color online) Band structure of the 5-AGNR without (black squares) and with (red triangles) ISO coupling obtained by first-principles calculations.Figure 6
(Color online) The relationship between the gaps induced by ISO and N-AGNRs from tight-binding approximation (red circles) and from First-principles calculations (black diamonds; energies are multiplied by 50).Next, we consider the ZGNRs with the widths forN=3,4,5,6. Without ISO coupling, although the bands are self-avoiding, the ZGNRs are metallic. In the presence of ISO interaction, the bands split, and a gap is opened. Figure 7 displays the band structures near Fermi surface of 5-ZGNR and 4-ZGNR in the presence and absence of the ISO interaction. Increasing the relativistic potential by 300 times, the band gap changes from 1.9×10-3 eV to 1.0×10-5 eV. Since the ISO-opened gap is the second power of the potential expected from the perturbation theory, final results which are obtained through dividing the gaps gotten from VASP by 90000 are 2.0×10-6 meV and 1.0×10-7 meV, respectively. This method was verified of being effectual by the first-principles calculations in [4]. Both first-principles calculations and TB approximation demonstrate that with the increase of the widths for even chains, the gaps decrease rapidly.(Color online) Band structures of the (a) 5-ZGNR and (b) 4-ZGNR without (black squares) and with (red triangles) ISO coupling obtained by first-principles calculations. The insets are a zoom in of the band spectra close to the Fermi surface.
(a)(b)


## 5. Summary

In summary, we have investigated the ISO effects in ZGNRs and AGNRs within four orbital TB approximation and first-principles calculations. In TB method, for metallic AGNRs and the ZGNRs with even chains, a gap is opened in the spectra by ISO coupling while for the ZGNRs with odd chains, conduction and valence edge bands are still degenerated atKZ=π/a. By first-principles calculations, the ISO coupling reduces the gaps for AGNRs with N=3M-1 and causes a gap in ZGNRs. The results show that the edge configuration and width affect the strength of the ISO coupling. The ISO interaction has a stronger effect for AGNRs than for ZGNRs. Besides, the influence of ISO coupling is weaker on GNRs with larger width, which the edge affects less.Note Added
Recently, a related preprint which addresses spin-orbit interactions in GNRs with more than 50 chains by four orbital TB method has been posted [26].



---

*Source: 364897-2011-04-28.xml*