Control of electronic transport in graphene by electromagnetic dressing K. Kristinsson,1 O. V. Kibis,2, 1, ∗ S. Morina,1, 3 and I. A. Shelykh1, 3, 4 1
Division of Physics and Applied Physics, Nanyang Technological University 637371, Singapore 2 Department of Applied and Theoretical Physics, Novosibirsk State Technical University, Karl Marx Avenue 20, Novosibirsk 630073, Russia 3 Science Institute, University of Iceland, Dunhagi-3, IS-107, Reykjavik, Iceland 4 ITMO University, St. Petersburg 197101, Russia
SUPPLEMENTARY INFORMATION: Full and consistent derivation of the energy spectrum and wave functions of dressed electrons in graphene
Let us consider a graphene sheet which lies in the plane (x, y) at z = 0 and is subjected to an electromagnetic wave propagating along the z axis (dressing electromagnetic field). Then electronic properties of the graphene are described by the Hamiltonian ˆ = vσ · (~k − eA), H
(1)
where σ = (σx , σy ) is the Pauli matrix vector, k = (kx , ky ) is the electron wave vector in the graphene plane, v is the electron velocity in graphene near the Dirac point, e is the electron charge, and A = (Ax , Ay ) is the vector potential of the electromagnetic wave in the graphene plane. Solving the Schr¨odinger equation with the Hamiltonian (1), we can find the energy spectrum of dressed electrons and their wave functions as follows. I. Circularly polarized dressing field
For the case of circularly polarized electromagnetic wave, its vector potential A = (Ax , Ay ) can be written as E0 E0 A= cos ωt, sin ωt , (2) ω ω where E0 is the electric field amplitude of the wave, and ω is the wave frequency. Then the Hamiltonian (1) is ˆ=H ˆ0 + H ˆk, H where ˆ0 = H and ˆk = H
−
veE0 ω
0
eiωt
(3) e−iωt , 0
0 v(~kx − i~ky ) . v(~kx + i~ky ) 0
(4)
(5)
The nonstationary Schr¨ odinger equation with the Hamiltonian (4), i~
∂ψ0 ˆ 0 ψ0 , =H ∂t
(6)
describes the time evolution of electron states in the Dirac point (k = 0). The exact solutions of the equation (6) can be sought in the form ψ0 = e−iαt/~ [AΦ′1 (r)e−iωt/2 + BΦ′2 (r)eiωt/2 ], Φ′1,2 (r)
(7)
where are the basic functions of the 2 × 2 matrix Hamiltonian (1), and α, A and B are the sought parameters. Substituting the wave function (7) into the Schr¨odinger equation (6), we arrive at the system of algebraic equations veE0 ~ω +B = 0 A α+ 2 ω ~ω veE0 A +B α− = 0. (8) ω 2
2 The condition of nontrivial solution of the system (8), α + ~ω veE0 ω veE 2 ~ω = 0, 0 α − ω 2 gives the two different parameters, α = ±Ω/2, where s Ω=
2
(~ω) +
2veE0 ω
2
.
(9)
Therefore, there are two sets of solutions of the system (27), which correspond to these two parameters and satisfy the normalization condition, |A|2 + |B|2 = 1. As a result, there are two wave functions (7), "r # r Ω ± ~ω ′ Ω ∓ ~ω ′ e ± ±iΩt/2~ −iωt/2 iωt/2 ψ0 = e Φ1 (r)e ± Φ2 (r)e , (10) 2Ω |e| 2Ω which exactly describe electron states of irradiated graphene in the Dirac point (k = 0). Since the two wave functions (10) are the complete function system for any time t, we can seek the solution of the Schr¨odinger equation with the full Hamiltonian (3) as an expansion ψk = a+ (t)ψ0+ + a− (t)ψ0− .
(11)
Substituting the expansion (11) into the Schr¨odinger equation with the full Hamiltonian (3), i~
∂ψk ˆ k, = Hψ ∂t
we arrive at the system of two differential equations for the coefficients a± (t), Ω + ~ω e W0 e−i(Ω/~−ω)t a− (t) (kx cos ωt + ky sin ωt)a+ (t) − (kx − iky ) ia˙ + (t) = v |e| Ω 2Ω Ω − ~ω + (kx + iky ) e−i(Ω/~+ω)t a− (t) , 2Ω e W Ω + ~ω 0 − − ia˙ (t) = −v ei(Ω/~−ω)t a+ (t) (kx cos ωt + ky sin ωt)a (t) + (kx + iky ) |e| Ω 2Ω Ω − ~ω i(Ω/~+ω)t + − (kx − iky ) e a (t) , 2Ω
(12)
(13)
where W0 =
2veE0 ω
is the characteristic kinetic energy of rotational electron motion induced by the circularly polarized field. In what follows, we will assume that the field frequency ω is high enough to satisfy the condition W0 /~ω ≪ 1.
(14)
W0 W0 Ω + ~ω Ω − ~ω ≈ ≈ 0, ≈ 1, ≈ 0, Ω ~ω 2Ω 2Ω
(15)
Then we have
and, correspondingly, Eqs. (13) take the form e (kx − iky )e−i(Ω−~ω)t/~ a− (t), |e| e ia˙ − (t) = −v (kx + iky )ei(Ω−~ω)t/~ a+ (t). |e| ia˙ + (t) = −v
(16)
3 ± −i(±Ω/2∓~ω/2+εk )t/~ , where C ± and εk are the sought parameters. Let us seek solutions of Eqs. (16) as a± c (t) = C e Solving the system of equations (16), we arrive at q εk = ± (εg /2)2 + (~vk)2 , (17)
and
s
2
(~ω) +
εg =
2veE0 ω
2
− ~ω.
(18)
Following the conventional terminology, Eq. (17) describes the spectrum of quasienergies of electrons in graphene dressed by a circularly polarized field, where signs “+” and “−” correspond to conduction band and valence band, respectively. The two wave functions (11), which correspond to the two energy branches (17), can be written as " s ! r r |εk | ∓ εg /2 −iθ/2 Ω + ~ω ′ Ω − ~ω ′ e −iεk t/~ iωt ψk = e ± e Φ1 (r) + Φ2 (r)e 2|εk | |e| 2Ω 2Ω s !# r r |εk | ± εg /2 iθ/2 Ω − ~ω ′ Ω + ~ω ′ e −iωt − e Φ1 (r)e − Φ2 (r) , (19) 2|εk | 2Ω |e| 2Ω where θ is the azimuth angle of the wave vector, k = (k cos θ, k sin θ). Let us write the basis functions Φ′1,2 (r) in the conventional Bloch form, Φ′1,2 (r) = Φ1,2 (r)ϕk (r), where Φ1,2 (r) are the periodical functions arisen from atomic π√ orbitals of the two crystal sublattices of graphene, ϕk (r) = eik·r / S is the plane electron wave, and S is the graphene area. Keeping in mind inequalities (15), we arrive from (19) to the sought wave functions of dressed electrons in the final form, s "s # |ε | ∓ ε /2 |εk | ± εg /2 iθ/2 k g −iεk t/~ −iθ/2 ψk = ϕk (r)e (20) e Φ1 (r) ± e Φ2 (r) . 2|εk | 2|εk | II. Linearly polarized dressing field
For the case of electromagnetic wave linearly polarized along the x axis, its vector potential A = (Ax , Ay ) can be written as E0 A= cos ωt, 0 . (21) ω Then the Hamiltonian (1) is ˆ=H ˆ0 + H ˆk, H
(22)
0 −1 veE0 cos ωt. −1 0 ω
(23)
0 v(~kx − i~ky ) . v(~kx + i~ky ) 0
(24)
where ˆ0 = H and ˆk = H
The nonstationary Schr¨ odinger equation with the Hamiltonian (23), i~
∂ψ0 ˆ 0 ψ0 , =H ∂t
(25)
describes the time evolution of electron states in the Dirac point (k = 0). The exact solutions of this Schr¨odinger equation can be sought in the form veE0 ′ ′ ψ0 = [AΦ1 (r) + BΦ2 (r)] exp −iα sin ωt , (26) ~ω 2
4 where α, A and B are the sought parameters. Substituting the wave function (26) into the Schr¨odinger equation (25) with the Hamiltonian (23), we arrive at the system of algebraic equations Aα + B = 0 A + Bα = 0.
(27)
The condition of nontrivial solution of the system (27), α 1 1 α = 0,
gives the two different parameters, α = ±1. Therefore, there are two sets of solutions of the system (27), which correspond to these two parameters and satisfy the normalization condition, |A|2 + |B|2 = 1: 1 A=B= √ , 2 1 A = −B = √ . 2
(28)
As a result, there are two wave functions (26), ψ0±
veE0 1 ′ ′ = √ [Φ1 (r) ± Φ2 (r)] exp ±i sin ωt , ~ω 2 2
(29)
which exactly describe electron states of irradiated graphene in the Dirac point (k = 0). Since the two wave functions (29) are the complete function system for any time t, we can seek the solution of the nonstationary Schr¨odinger equation with the full Hamiltonian (22) as an expansion ψk = a+ (t)ψ0+ + a− (t)ψ0− .
(30)
Substituting the expansion (30) into this Schr¨odinger equation with the full Hamiltonian (22), i~
∂ψk ˆ k, = Hψ ∂t
we arrive at the system of two differential equations for the coefficients a± (t), 2veE0 + + − ia˙ (t) = vkx a (t) + ivky a (t) exp −i sin ωt , ~ω 2 2veE0 sin ωt − vkx a− (t). ia˙ − (t) = −ivky a+ (t) exp i ~ω 2
(31)
(32)
It follows from the Floquet’s theorem that the functions a± (t) can be written as a± (t) = e−iεk t/~ e a± (t),
(33)
where εk is the electron quasienergy in the irradiated graphene (the energy spectrum of dressed electrons), and e a± (t) are the periodical functions with the period T = 2π/ω. Let us apply the Jacobi-Anger expansion, eiz sin γ =
∞ X
Jn (z)einγ ,
n=−∞
to the exponents in the right side of Eqs. (32) and expand the periodical functions e a± (t) into the Fourier series ±
e a (t) =
∞ X
n=−∞
inωt c± . ne
5 Then the differential equations (32) can be transformed to the algebraic equations ∞ X εk 2veE0 − = 0, − − nω c+ + ivk c J y n n−m m ~ ~ω 2 m=−∞ ∞ X εk 2veE0 + − vkx + = 0, + nω cn + ivky cn−m Jm ~ ~ω 2 m=−∞
vkx −
(34)
where Jm (z) is the Bessel function of the first kind. Equations (34) can be easily solved in the case of high-frequency field satisfying the conditions ~ω ≫ εk , ω ≫ vk.
(35)
In this simplest case, Eqs. (34) for n 6= 0 can be reduced to the equation c± n ≈i
vky nω
2veE0 . ∓ c∓ J n−m m ~ω 2 m=−∞ ∞ X
(36)
± Keeping in mind that |c± n | ≤ 1 and |Jn (z)| ≤ 1, Eq. (36) leads to the evident solution cn6=0 ≈ 0. After substitution of this solution into Eqs. (34), the expressions (34) are reduced to the two equations εk + 2veE0 vkx − J c0 + ivky c− = 0, 0 0 ~ ~ω 2 2veE0 εk − J c0 + ivky c+ vkx + = 0. 0 0 ~ ~ω 2 (37)
Solving the system of the two algebraic equations (37) accurately, we can easily obtain both the coefficients c± 0 and the energy spectrum of dressed electrons, εk = ±~vkf (θ),
(38)
where f (θ) =
s
cos2 θ + J02
2veE0 ~ω 2
sin2 θ.
Correspondingly, the sought wave functions of dressed electrons (30) take the final form s 2 cos θ + f (θ) −iεk t/~ [Φ1 (r) ± Φ2 (r)] e±i(veE0 /~ω ) sin ωt ψk = ϕk (r)e 4f (θ) ! 2veE0 sin θ ∓i(veE0 /~ω 2 ) sin ωt [Φ1 (r) ∓ Φ2 (r)] e . J0 − i cos θ + f (θ) ~ω 2
(39)
Taking into account Eq. (38), the two conditions (35) can be reduced to the solely condition, ω ≫ vk, which describes borders of applicability of Eqs. (38) and (39).
∗
Electronic address:
[email protected]
(40)