^{1}

^{2}

^{*}

We investigate the emerging consequences of an applied strong in-plane electric field on a macroscopically large graphene sheet subjected to a perpendicular magnetic field, by determining in exact analytical form various many-body thermodynamic properties and the Hall coefficient. The results suggest exotic possibilities that necessitate very careful experimental investigation. In this alternate form of Quantum Hall Effect, non-linear phenomena related to the global magnetization, energy and Hall conductivity (the latter depending on the strengths of magnetic B- and electric E-fields) emerge without using perturbation methods, to all orders of E-field and B-field strengths. Interestingly enough, when the value of the electric field is sufficiently strong, fractional quantization also emerges, whose topological stability has to be verified.

Dirac-type materials, such as Topological Insulators, monolayer graphene, and three-dimensional (3D) Dirac and Weyl semimetals, appear nowadays as stable (actually very robust) topological phases of matter, displaying behavioral patterns that produce new physics at a very fundamental level and at the same time give the possibility of exotic future applications [

Let us start with the graphene energy spectrum when a monolayer is subjected to an in-plane electric E (taken along the x-direction) and a perpendicular magnetic field B in the z-direction (and let us focus on the positive branch, and also ignore the Zeeman interaction term), and take the Landau gauge A = (0, xB, 0) in which the energy spectrum turns out to be (through a Dirac equation procedure similar to the one in [

ε n , k y = 2 n ℏ e B u f ( 1 − β 2 ) 3 / 4 + ℏ u f β k y , (1.1)

with β = E / u f B a dimensionless parameter (always supposed to be lower than unity, β < 1 ), n = 0, 1, 2, 3. The Landau Level index for the positive branch, u f is the Fermi velocity, and k y is the wave vector along the y-direction. We also find that the guiding center operator’s eigenvalue (projected on the x-axis) X_{0} reads

X 0 = l B 2 k y − sgn ( n ) β l B 2 | n | ( 1 − β 2 ) 1 / 4 , (1.2)

with l B = ℏ c / e B the magnetic length and sgn ( n ) is the sign function. Due to the spatial confinement in the x-direction, the guiding center operator X 0 may acquire any value in the following range:

− L x 2 ≤ X 0 ≤ L x 2 (1.3)

with L x being the x-direction size of the system, which is here supposed to be macroscopically large. Each Landau Level, defined by different values of index n, contains Φ / Φ 0 independent quantum states, with Φ = B S the magnetic flux penetrating the 2D graphene sheet (S is the area) and Φ 0 = h c / e is the flux quantum. Now, because of the spin and valley degeneracies that are present in graphene, each Landau Level (L.L.) may accommodate up to 4 Φ / Φ 0 spinful electrons, according to the Pauli exclusion principle excluding the lowest level n = 0, which can accommodate only up to 2 Φ / Φ 0 spinful electrons, due to mutual sharing with the holes.

As a next step it is convenient (for the thermodynamic calculations followed below), to express (1.1) as a function of the guiding center operator, namely

ε n , X 0 = 2 n ℏ e B u f ( 1 − β 2 ) 1 / 4 + e E X 0 . (1.4)

A few remarks are then in order about the energy gap (the inter-L.L. gap, for a constant guiding center) determined by:

δ ε n = 2 ℏ e B u f ( 1 − β 2 ) 1 / 4 ( n + 1 − n ) . (1.5)

Unlike conventional semiconductors, the energy gap has an E-field dependence. As can easily be seen from (1.5), the larger the E-field gets, the larger the energy gap becomes. On the other hand, the larger the L.L. index, the smaller the energy gap. This interplay will play a major role later on, when we consider the thermodynamic occupations of the energy levels. One can always prove that, for a given L.L. index and a value of E-field (such that β < 1 ) there will always be an unavoidable overlap (states of greater n values have lower energy than states with lower n values). We can set conditions (for arbitrary E and n) for which this kind of overlap is avoided as:

ε n , L x / 2 ≤ ε n + 1 , − L x / 2 (1.6)

To picture this more properly, we provide examples in the graphs shown in

Condition (1.6) along with (1.4) results in:

2 ℏ e B u f ( 1 − β 2 ) 1 / 4 ( n + 1 − n ) ≥ e E L x or δ ε n ≥ e E L x (1.7)

That is, in words, when the work performed by the electric field is smaller than the energy gap at a certain X 0 , no overlap is observed between the L.L.s n and n + 1. Generally, for strong enough electric field or for a small enough L.L. index, the above inequality will be true (see

We define n = i F to be the topmost L.L., which for a given value of E and B maintains an energy gap with its adjacent L.L.s: ( n = i F − 1 , n = i F + 1 ) Clearly, L.L. index n = i F + 1 no longer separates from i F + 2 with an energy gap, and in this case, Inequality (1.7) reverses direction: δ ε i F + 1 ≤ e E L x . In other words, L.L.s with n = 0 , 1 , 2 , ⋯ , i F do not overlap, while L.L.s with higher quantum numbers ( n = i F , i F + 1 , i F + 2 , ⋯ ) do overlap.

The above is a generic case for an arbitrary value of E and n. Of course there might be cases where overlaps start from n = 0 (for a low enough E-field), in which case i F = 0 , and all L.L.s with n > 0 overlap. In this case, with no energy gap present at all, graphene will gain a metallic character. Equation (1.7) with input n = i F + 1 becomes then:

2 ℏ e B u f ( 1 − β 2 ) 1 / 4 ( i F + 2 − i F + 1 ) ≤ e E L x (1.8)

After the above discussion and definitions, we proceed to thermodynamic occupations of the graphene’s energy levels at zero temperature. For this purpose, we consider a collection of N electrons at T = 0, which fill the lowest energy levels until the Fermi energy denoted by ε F . In reality, the Fermi energy is not constant when there is an electric field running through the system; what we then mean by Fermi energy is actually a Fermi point, which is the topmost occupied state in the energy diagram. We also make the supposition that this Fermi point is located on a L.L. indexed with n = ρ − 1 ( ρ = 1 , 2 , 3 , ⋯ ), so that there are always ρ L.L.s occupied at any time (the last level n = ρ − 1 being generally partially occupied). First, we will focus on the special case where all ρ Landau Levels are not overlapping, and can be occupied independently by the N electrons. In this case the following relation must be satisfied:

ρ − 1 ≤ i F (2.1)

and for strong enough magnetic fields, it is guaranteed that Equation (2.1) will always be satisfied, and no overlap between L.L.s with different quantum numbers will be observed. In what follows, we will consider a constant, strong E-field, while the magnetic field may vary, but always in a way that satisfies Equation (2.1).

Considering that the L.L. with n = 0 only has a capacity for 2Φ/Φ_{0} electrons, and that all the other n ≠ 0 L.L.s may host up to 4Φ/Φ_{0} electrons, we find that in order to have ρ L.L.s occupied, the following inequality must hold:

( 2 ρ − 3 ) 2 Φ Φ 0 ≤ N ≤ ( 2 ρ − 1 ) 2 Φ Φ 0 (2.2)

Note that when ρ = 1 then 0 ≤ N ≤ 2 Φ Φ 0 , as N is always a positive number.

Treating N as a constant, we can solve (2.2) with respect to magnetic field B:

1 2 ( 2 ρ − 3 ) n A Φ 0 ≥ B ≥ 1 2 ( 2 ρ − 1 ) n A Φ 0 (2.3)

with n A the electronic surface density and Φ 0 the flux quantum. Also note that we have considered the special case where ρ ≤ i F + 1 , where all L.L.s are well separated with an energy gap, and no inter-L.L. overlapping occurs. (The special case of nonzero overlaps will be considered in the next Section). Note that Equation (2.3) can also describe the well-known unconventional Quantum Hall Effect in graphene, with the Hall conductivity given by the well-known relation:

σ H = e n A B = ( ρ − 1 2 ) 4 e 2 h (2.4)

If one replaces B by the value B = 1 2 ( 2 ρ − 1 ) n A Φ 0 ; σ H is then quantized in half-integer multiples of 4 e 2 h .

In the case we are considering the total energy of the system (minimized at T = 0) is given as a sum over all states lower than “the Fermi state” or Fermi point, namely

E = ∑ n , X 0 ε n , X 0 = ∑ n , X 0 2 n ℏ e B u f ( 1 − β 2 ) 1 / 4 + e E X 0 (2.5)

In the thermodynamic limit L x → ∞ , we may approximate the sum with respect to X 0 as follows:

∑ X 0 → 4 B L y Φ 0 ∫ − L x 2 L x 2 d X 0 for n > 0 and ∑ X 0 → 2 B L y Φ 0 ∫ − L x 2 L x 2 d X 0 for n = 0 . (2.6)

The total internal energy of the system can then be separated in the energy of the fully occupied bands plus the energy of the partially occupied last L.L. (that contains the Fermi point, n = ρ − 1 ):

E = E full + E part (2.7)

In what follows, we will consider the case ρ > 1 , where the contribution of the lowest L.L. n = 0 is negligible. Given that all the bands up to n = ρ − 2 are fully occupied, we may write:

E full = 4 B L y Φ 0 ∑ i = 1 ρ − 2 ∫ − L x 2 L x 2 d X 0 2 n ℏ e B u f ( 1 − β 2 ) 1 / 4 + 4 B L y e E Φ 0 ∑ i = 1 ρ − 2 ∫ − L x 2 L x 2 X 0 d X 0 = 4 Φ Φ 0 2 ℏ e B u f ( 1 − β 2 ) 1 / 4 ∑ n = 1 ρ − 2 n (2.8)

To determine E part , we must first determine the X 0 value at the Fermi point (limit of the X 0 integration at the n = ρ − 1 L.L.). From Equation (1.2) we get for k y = 2 π l / L y and n = ρ − 1 :

X 0 = l B 2 2 π l L y − β l B 2 ( ρ − 1 ) ( 1 − β 2 ) 1 / 4 (2.9)

The quantum number l, appearing in (2.9) has a starting value l 0 that can be

determined by setting X 0 = − L x 2 as follows:

l 0 = − Φ 2 Φ 0 + L y 2 π l B β 2 ( ρ − 1 ) ( 1 − β 2 ) 1 / 4 (2.10)

Now, using the fact that the last L.L. contains N − ( 2 ρ − 3 ) 2 Φ Φ 0 electrons, we

have all the necessary information to determine the guiding center value of the electron placed at the Fermi point:

X 0 , F = 2 π L y l B 2 [ N 4 − ( ρ − 1 ) Φ Φ 0 ] (2.11)

For example, when N = 2 ( 2 ρ − 3 ) Φ Φ 0 , (meaning that the last L.L. is empty of

electrons), we get from (2.11):

X 0 , F = 2 π L y l B 2 [ 2 ( 2 ρ − 3 ) 4 − ( ρ − 1 ) ] Φ Φ 0 = − π L y l B 2 Φ Φ 0 = − L x / 2 , (2.12)

whereas for N = 2 ( 2 ρ − 1 ) Φ Φ 0 we have that X 0 , F = L x / 2 , explicitly demonstrating

the correctness of our results. The remaining task to carry out is to calculate the energy of the partially occupied L.L. n = ρ − 1

E part = 4 B L y Φ 0 ∫ − L x 2 X 0 , F d X 0 2 ( ρ − 1 ) ℏ e B u f ( 1 − β 2 ) 1 / 4 + e E 4 B L y Φ 0 ∫ − L x 2 X 0 , F X 0 d X 0 = 4 B L y Φ 0 2 ( ρ − 1 ) ℏ e B u f ( 1 − β 2 ) 1 / 4 ( X 0 , F + L x 2 ) + e E 4 B L y Φ 0 ( X 0 , F 2 2 − L x 2 8 ) (2.13)

Substituting Equation (2.11) into (2.13) we obtain (in units of Fermi energy (in the absence of the electromagnetic field) ε F = ℏ u F π n A , per particle)

E part N = ε F ( 1 − β 2 ) 1 / 4 ( ( B n A Φ 0 ) 1 / 2 4 ( ρ − 1 ) − 2 4 ( ρ − 1 ) ( B n A Φ 0 ) 3 / 2 [ 2 ρ − 3 ] ) + ( N E 8 L y h B − e S E L y ( ρ − 1 ) + e E B L x 2 n A Φ 0 [ 4 ρ 2 − 8 ρ + 3 ] ) (2.14)

Finally, adding Equations (2.14) and (2.8) we arrive at the following result:

E N = ε F ( 1 − β 2 ) 1 / 4 ( 4 ( B n A Φ 0 ) 3 / 2 [ 2 ∑ n = 1 ρ − 2 n − [ 2 ρ − 3 ] ( ρ − 1 ) ] + ( B n A Φ 0 ) 1 / 2 4 ( ρ − 1 ) ) + ( N E 8 L y h B − e S E L y ( ρ − 1 ) + e E B L x 2 n A Φ 0 [ 4 ρ 2 − 8 ρ + 3 ] ) (2.15)

with ∑ n = 1 ρ − 2 n = − ζ ( − 1 2 , ρ − 1 ) − ζ ( 3 2 ) 4 π where ζ ( − 1 2 , ρ − 1 ) is the Hurwitz Zeta function and ζ ( 3 2 ) is the Riemman function. Note the interesting fact

that, terms proportional to electric field strength E result in the following magnetization:

M E N = − ∂ E N ∂ B = N E 8 L y h B 2 − e E L x 2 n A Φ 0 [ 4 ρ 2 − 8 ρ + 3 ] , (2.16)

which, when considered at full band occupation (meaning that B = 1 2 ( 2 ρ − 1 ) n A Φ 0 ) yields:

M E N = N E 8 L y 4 h ( 2 ρ − 1 ) 2 n A 2 Φ 0 2 − e E L x 2 n A Φ 0 [ 4 ρ 2 − 8 ρ + 3 ] = E L x n A ( ρ − 1 2 ) 2 e 2 h = L x n A σ H 2 E (2.17)

i.e. the proportional constant is equal to half Hall conductivity, similar to the corresponding result in a conventional semiconductor case [

We now proceed to a considerably more difficult case: the low E-field regime. In this regime, as the electric field becomes weaker, the energy gap gets smaller. As a result, there will be an unavoidable point where some L.L.s will overlap (

We fix the Fermi point at n = ρ − 1 , with energy given by:

ε F = 2 ( ρ − 1 ) ℏ e B u f ( 1 − β 2 ) 1 / 4 + e E X 0 F (3.1)

with X 0 F the guiding center position of the last energetically highest electron, occupying the Fermi point state. For the example case shown in

ε F = 8 ℏ e B u f ( 1 − β 2 ) 1 / 4 − e E L x 2 (3.2)

Now, from the

6 ℏ e B u f ( 1 − β 2 ) 1 / 4 + e E X 0 F = 8 ℏ e B u f ( 1 − β 2 ) 1 / 4 − e E L x 2 (3.3)

which yields for X 0 F :

X 0 F = − L x 2 + ℏ e B u f e E ( 1 − β 2 ) 1 / 4 ( 8 − 6 ) (3.4)

From the above, we may determine the initial and final values of the k x (that is, l 0 and l f ), namely

l 0 = − Φ 2 Φ 0 + L y 2 π l B 6 β ( 1 − β 2 ) 1 / 4 (3.5)

l f = − Φ 2 Φ 0 + B L y E h ℏ e B u F ( 1 − β 2 ) 1 / 4 ( 8 − 6 ) + E L y h u F e B e ℏ 6 ( 1 − β 2 ) 1 / 4 (3.6)

The number of states in L.L. n = 3 is then given by the difference of the initial and final values of l, namely

l f − l 0 = B L y E h ℏ e B u F ( 1 − β 2 ) 1 / 4 ( 8 − 6 ) (3.7)

and the total number of states below the Fermi point reads:

g = Φ Φ 0 | n = 0 + 2 Φ Φ 0 | n = 1 , 2 + B L y E h ℏ e B u F ( 1 − β 2 ) 1 / 4 ( 8 − 6 ) (3.8)

Because all these states are filled with electrons (4 electrons in each state for n > 0 and 2 electrons for n = 0) we have that:

N = 2 Φ Φ 0 + 8 Φ Φ 0 + 4 B L y E h ℏ e B u F ( 1 − β 2 ) 1 / 4 ( 8 − 6 ) = 10 Φ Φ 0 + 4 B L y E h ℏ e B u F ( 1 − β 2 ) 1 / 4 ( 8 − 6 ) (3.9)

Equation (3.9) yields for the Hall conductivity:

σ H = e n A B = 10 e 2 h + 4 e L y S E h ℏ e B u F ( 1 − β 2 ) 1 / 4 ( 8 − 6 ) (3.10)

Although (3.10) does not necessarily imply that irrational quantization in graphene is possible, it is quite interesting to notice how the integer quantization (the first term) is destroyed as long as L.L.s become intermixed due to the extra overlaps―which in turn are induced by the E-field. The least one can gain out of these calculations is the importance of the energy gap, and the possible formation of localized and extended states that lie inside this gap. When the gap is destroyed, the system becomes metallic; it may thus not be unreasonable to find a Hall conductivity that is electric/magnetic field-dependent, destroying the plateaux formation. This is something that needs to be further investigated experimentally.

In this work, a thermodynamic study has been conducted with respect to a 2D Graphene monolayer subjected to crossed electric and magnetic fields. Thermodynamic quantities like the global energy and magnetization have been exactly and analytically determined, in combination with transport properties, i.e. the Hall conductivity. The results suggest exotic possibilities that are here pointed out as a natural outcome of an exact and careful calculation in the noninteracting electrons many-body framework, and these necessitate exceedingly careful experimental investigation. Finally, with respect to the range of validity, although our results involve no approximations whatsoever, the general role of disorder in combination with the above physics of the strong fields is certainly something that needs further study.

Konstantinou, G. and Moulopoulos, K. (2018) Strong Electric Field in 2D Graphene: The Integer Quantum Hall Regime from a Different (Many-Body) Perspective. Advances in Materials Physics and Chemistry, 8, 32-43. https://doi.org/10.4236/ampc.2018.81003