Quantum Magneto-Optics in Graphene

: The optical conductivity of graphene in quantizing magnetic ﬁelds is studied. Both dynamical conductivities, longitudinal and Hall’s, are analytically evaluated. The conductivity peaks are explained in terms of electron transitions. The optical transitions obey the selection rule with ∆ n = 1 for the Landau number n . The light transmission and Faraday rotation in the quantizing magnetic ﬁelds are calculated.


Introduction
Comprehensive literature on the graphene family is expressed usually in terms of the Dirac gapless fermions. According to this picture, there are two bands at the K hexagon vertexes of the Brillouin zone without any gap between them, and the electron dispersion can be considered as linear in the wide wave-vector region. For the dispersion linearity, this region should be small compared with the size of the graphene Brillouin zone, i.e., less than 10 −8 cm −1 , providing the small carrier concentration n 10 16 cm −1 . Pristine graphene at zero temperature has no carriers, and the Fermi level should divide the conduction and valence bands. However, undoped graphene cannot be really obtained, and so far, the purest graphene contains about n ∼ 10 9 cm −2 of carriers. Then, the following problem appears: how do Coulomb electron-electron interactions renormalize the linear dispersion, and does graphene become an insulator with a gap?
Semiconductors with a gap are needed for electronic applications. Investigations of the graphene bilayer and multilayer are very popular, as the gap appears when the bias is applied. We see how physics made a circle for a half century, returning to graphite studies. Here, Slonczewski, Weiss and McClure (SWMC) should be mentioned, because they have stated the description of a layered matter [1] with interactions strong in the layer and weak between layers.
The most accurate investigation of the band structure of metals and semiconductors is a study of the Landau levels through experiments, such as magneto-optics [2][3][4][5][6][7][8][9][10], magneto-transport [11][12][13][14][15] and magnetic oscillation of the optical phonon [16,17]. In magnetic fields, the classical and quantum Hall effects are observed, as well as the polarization rotation for transmitted (Faraday rotation) or reflected light (Kerr rotation). However, the interpretation of the experimental results involves a significant degree of uncertainty, because it is not clear how the resonances can be identified and to which electron transitions they correspond.
The theoretical solution for the band problem in magnetic fields often cannot be exactly found. A typical example is presented by graphene layers. For bilayer graphene and graphite, the effective Hamiltonian is a 4 × 4 matrix, giving four energy bands. The trigonal warping described by the effective Hamiltonian with a relatively small parameter γ 3 provides an evident effect. Another important parameter is the gate-tunable bandgap U in bilayer graphene. In this situation, the quantization problem cannot be solved within a rigorous method. To overcome this difficulty, several methods have been proposed for approximate [9,[18][19][20][21], numerical [22][23][24][25][26] and semiclassical quantization [27,28].
The present paper is organized as follows. In Section 2, we recall the electron dispersion in graphene. In Section 3, we describe in detail the quantization in magnetic fields. The optical conductivity, light transmission, and Faraday rotation are discussed in Section 4.

Electron Dispersion in Graphene
The symmetry of the K point is C 3v with the three-fold axis and reflection planes. This group has a two-fold representation with the basis functions transforming each in other under reflections and obtaining the factors exp (±2πi/3) in rotations. The linear momentum variations from the K point p ± = ∓ip x − p y transform in a similar way. The effective Hamiltonian is invariant under the group transformations, and we have the unique possibility to construct the invariant Hamiltonian linear in the momentum as: where v is a constant of the velocity units. The same Hamiltonian was written using the tight binding model. The eigenvalues of this matrix give two bands: where the sign ∓ corresponds to holes and electrons. The gapless linear spectrum arises as a consequence of the symmetry, and the chemical potential at zero temperatures coincides with the band crossing due to the carbon valence. The cyclotron mass has the form: and the carrier concentration at zero temperature: n(µ) = µ 2 /πh 2 v 2 is expressed in terms of the chemical potential µ.
Tuning the gate voltage, the linearity of the spectrum has been examined in the Shubnikov-de Haas studies [29] with the help of the connection between the effective mass and the carrier concentration at the Fermi level m(µ)v = ∓h πn(µ). The "constant" parameter v was found to be no longer constant, but at low carrier concentrations n ∼ 10 9 cm −2 , it exceeds its usual value v = 1.05 ± 0.1 × 10 8 cm/s (at concentrations n > 10 11 cm −2 ) by a factor of three.
This is a result of electron-electron interactions, which becomes stronger at low carrier concentrations. The logarithmic renormalization of the velocity was found by Abrikosov and Beneslavsky in [30] for the three-dimensional case and in [31,32] for two-dimensional graphene. Notice that no phase transition was revealed, even at the lowest carrier concentration. We can conclude that the Coulomb interactions do not create any gap in the spectrum.

Graphene in Magnetic Fields
In the presence of the magnetic field B, the momentum projections p + and p − become the operators with the commutation rule {p + ,p − } = −2ehB/c. We use the relations: involving the creation a + and annihilation a operators with ω B = v 2|e|hB/c. We will write only one of two x, y space coordinates, including the corresponding degeneracy proportional to the magnetic field in the final results.
For graphene, we search the eigenfunction of Hamiltonian (1) in the form: where ϕ n (x) are orthonormal Hermitian functions with the Landau number n ≥ 0. Canceling the Hermitian functions from the equations, we obtain a system of linear equations for the eigenvector C sn : giving the eigenvalues: ε sn = ∓ω B √ n with s = 1, 2 and n = 0, 1, 2... For n = 0, there is only one level ε 10 = 0 with C 1 0 = 0, C 2 0 = 1, as follows from Equation (2). The wave function columns are written:

Magneto-Optic Effects in Graphene
An important peculiarity of conductivities in the presence of magnetic fields is an appearance of the Hall component σ xy (ω). The Hall conductivity violates the rotation symmetry of graphene around the major axis. This implies the rotation of the linear polarized electromagnetic wave, i.e., the Faraday and Kerr effects for transmitted and reflected waves, correspondingly. Electron transitions are possible between the levels with the neighboring Landau numbers n and various bands s, and therefore, the resonance denominators ∆ ss n = ε sn − ε s ,n+1 arise in the conductivity tensor.
Calculations [18] give the conductivities for graphite in the collisionless limit when the electron collision frequency Γ is much less than the level splitting: n,s,s ∆f ss n ∆ ss n Here, ∆f ss n = f (ε s n+1 ) − f (ε sn ) is the difference of the Fermi functions and: is the universal graphene conductivity. For low carrier concentrations, e.g., for the low chemical potential, the electron transitions take place between the certain Landau levels obeying the condition ∆n = 1. For instance, the Landau levels at 7 Tin the upper band are ε 2,n = 0, 98, 138, and 169 meV for n = 0 to n = 3 and in the low band ε 1,n = 0, −98, −138, and -169 meV.
The transmission coefficient: and the Faraday rotation angle: calculated for those cases when the chemical potential increases and the carriers fill sequentially the levels from n = 0 to n = 3 are shown in Figures 2 and 3. For instance, if the level with n = 0 is occupied, the transition into the level |2, 1 is possible (at 98 meV in Figures), as well as from the level |1, 2 (at 236 meV, is not shown in the Figures). We use the notation |s, n for the levels with the band symbol s = 1, 2 and the Landau number n.  At high carrier concentrations, the chemical potential µ is much larger than the intraband level splitting ∆ ssn = eB/cm(ε) ≡ Ω, called the cyclotron frequency, with the effective mass m(ε) = ε/v 2 , and the conductivities (3) coincide with the semiclassical Drude conductivities: with ω * = ω + iΓ.

Discussion
Graphene layers affect the transmission and the Faraday rotation in a linear order in the fine structure constant α, whereas the reflected light intensity is quadratic to α. Thus, we discuss here the characteristics of the transmitted light. Expressed linearly in conductivity σ 0 = e 2 /4h, the transmission coefficient T and the Faraday angle for free standing graphene are shown in Figures 2 and 3. The peaks in absorption (Figure 2) correspond to the electron transitions. The series of lines are observed in the 0.1-0.4 eV interval. They are doublets excited by the electron transitions of the type |1, n → |2, n + 1 and |1, n + 1 → |2, n for n > 0. All of these lines obey the selection rule ∆n = 1. Let us note that the transition frequencies in Figure 2 determine the maximum values of the derivative in the Faraday rotation ( Figure 3).

Conflicts of Interest
The author declare no conflict of interest.