**Photoinduced Pseudospin Dynamical Effects in Graphene-Like Systems**

Alexander López and Benjamin Santos

Additional information is available at the end of the chapter

http://dx.doi.org/10.5772/67618

#### Abstract

In this chapter, we describe some of our recent results on the laser-induced manipulation of the energy band structure of graphene-like systems. We present numerical results on the quasi-energy spectrum as well as detailed calculations of semi-analytical approximations to other physical quantities of interest. The main message we would like to convey to the interested reader of the chapter is that by properly tuning the perturbation parameters of the radiation field one can control the size and shape of the photoinduced gaps. These in turn would allow the realization of new electronic phases on graphene and its related materials such as silicene.

Keywords: graphene, light-matter interaction, Floquet theory, Dirac fermions, Landau levels, pseudospin polarization

#### 1. Introduction

Since its recent experimental realization [1], graphene has attracted a lot of research interest because of its remarkable transport properties. At low energy, single-layer graphene (SLG) has a linear dispersion spectrum, and charge carriers can be described as massless, chiral, Dirac fermions. Moreover, SLG is a zero-gap semiconductor, has a high mobility of charge carriers and features an unconventional half-integer quantum Hall effect which can be measured at room temperature [2]; therefore, it has been the focus of many experimental [3] and theoretical analyses [4, 5], and it is expected to have several potential applications in carbon-based spintronics. Thus, the study of spin-related transport phenomena is one of the most active research fields in graphene. One important breakthrough was the demonstration of a quantum spin Hall effect in

© 2017 The Author(s). Licensee InTech. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

SLG [6]. This in turn relies on intrinsic spin-orbit coupling (SOC), which although weak compared to other energy scales in the problem [7] opens a gap in the energy spectrum, making SLG into a topological insulator. Another possibility for SOC is due to interactions with a substrate, the presence of electric fields or curvature of the SLG. Therefore, it is called extrinsic or Rashba spin-orbit coupling (RSOC). The RSOC is believed to be responsible for spin polarization [8], and spin relaxation phenomena in graphene [9]. In addition, in the presence of an external magnetic field B, Landau's quantization leads to a ffiffiffi B <sup>p</sup> energy relation which contrasts the linear in <sup>B</sup> dependence of conventional 2DEG, with parabolic energy bands.

Indeed, one interesting aspect of the new graphene-like materials is that they are predicted to support the so-called edge states which are related to the realization of new topological aspects of matter. As it has been acknowledged in the recent literature, since the discovery of the quantum Hall effect [10, 11], the quest for new topological states of matter has attracted a great deal of attention (for a recent review, see Ref. [12]). For instance, some proposals to topologically generate and characterize non-trivial phases in graphene have paved the road for research on how to implement models of non-trivial topological states in different physical systems ranging from semiconducting quantum wells [13], superconductors [14] and neutral [15] and cold atoms [16], just to mention a few.

In Ref. [13], a model for quantum spin Hall (QSH) state in telluride-based quantum wells was analysed. Yet, it was shown afterwards that this model really corresponds to a topological insulator state [17]. In general, a topological insulator is a quantum state of matter where a given material supports edge states that counter-propagate on its boundary. The existence of these edge states allows for non-dissipative transport, which is a most-wanted property in technological implementations, for instance, in quantum information and nanotechnological devices. In principle, the topological-insulating phase is protected by time-reversal symmetry. However, a time-reversal-symmetry-broken QSH state proposal was recently put forward by means of ferromagnetic leads attached to the sample [18].

Upon introduction of time-periodic dynamical modulation, for instance, in semiconducting quantum wells, with a zincblende structure, it has been recently shown that AC driving can induce a topological phase transition leading to the so-called Floquet topological insulator (FTI) phases [19]. This means that a non-trivial topological phase can be induced in a system that in equilibrium behaves trivially [20–22], that is, it does not possess non-dissipative or gapless edge states.

In this chapter, we show some of our results on the photoinduced effects on the Dirac fermions of graphene-like systems presenting two specific scenarios. In Section 2, we describe the Landau levels (LLs) in graphene when a laser is incident perpendicularly to the sample [23]. We show that an exact effective Floquet Hamiltonian can be found from which the dynamics is afterwards described [24–26]. Here, we show that the quasi-energy spectrum presents a leveldependent photoinduced gap and we also show that coherent Landau level states can be synchronized in order to produce Rabi oscillations and quantum revivals [27]. In Section 3, we describe how the energy spectrum and pseudospin polarization in monolayer silicene [28– 33] can be manipulated beyond the so-called off-resonant regime [34, 35] when strong radiation effects are taken into account [36–38]. We find the explicit and realistic parameter regime for the realization of a single-valley polarized state in silicene. In Section 4, an outlook is given whereas in Section 5 we present some technical details of the calculations employed because although we know there are several works describing the dynamical aspects of periodically driven systems, we wanted to be self-contained. Another reason for delving into such details is that we think this explicit technical aspect could be useful for both students and researchers interested in the field.

#### 2. Technical aspects of periodically driven systems

#### 2.1. Floquet Fourier-mode approach

SLG [6]. This in turn relies on intrinsic spin-orbit coupling (SOC), which although weak compared to other energy scales in the problem [7] opens a gap in the energy spectrum, making SLG into a topological insulator. Another possibility for SOC is due to interactions with a substrate, the presence of electric fields or curvature of the SLG. Therefore, it is called extrinsic or Rashba spin-orbit coupling (RSOC). The RSOC is believed to be responsible for spin polarization [8], and spin relaxation phenomena in graphene [9]. In addition, in the presence of an external magnetic

B

Indeed, one interesting aspect of the new graphene-like materials is that they are predicted to support the so-called edge states which are related to the realization of new topological aspects of matter. As it has been acknowledged in the recent literature, since the discovery of the quantum Hall effect [10, 11], the quest for new topological states of matter has attracted a great deal of attention (for a recent review, see Ref. [12]). For instance, some proposals to topologically generate and characterize non-trivial phases in graphene have paved the road for research on how to implement models of non-trivial topological states in different physical systems ranging from semiconducting quantum wells [13], superconductors [14] and neutral

In Ref. [13], a model for quantum spin Hall (QSH) state in telluride-based quantum wells was analysed. Yet, it was shown afterwards that this model really corresponds to a topological insulator state [17]. In general, a topological insulator is a quantum state of matter where a given material supports edge states that counter-propagate on its boundary. The existence of these edge states allows for non-dissipative transport, which is a most-wanted property in technological implementations, for instance, in quantum information and nanotechnological devices. In principle, the topological-insulating phase is protected by time-reversal symmetry. However, a time-reversal-symmetry-broken QSH state proposal was recently put forward by

Upon introduction of time-periodic dynamical modulation, for instance, in semiconducting quantum wells, with a zincblende structure, it has been recently shown that AC driving can induce a topological phase transition leading to the so-called Floquet topological insulator (FTI) phases [19]. This means that a non-trivial topological phase can be induced in a system that in equilibrium behaves trivially [20–22], that is, it does not possess non-dissipative or

In this chapter, we show some of our results on the photoinduced effects on the Dirac fermions of graphene-like systems presenting two specific scenarios. In Section 2, we describe the Landau levels (LLs) in graphene when a laser is incident perpendicularly to the sample [23]. We show that an exact effective Floquet Hamiltonian can be found from which the dynamics is afterwards described [24–26]. Here, we show that the quasi-energy spectrum presents a leveldependent photoinduced gap and we also show that coherent Landau level states can be synchronized in order to produce Rabi oscillations and quantum revivals [27]. In Section 3, we describe how the energy spectrum and pseudospin polarization in monolayer silicene [28– 33] can be manipulated beyond the so-called off-resonant regime [34, 35] when strong radiation effects are taken into account [36–38]. We find the explicit and realistic parameter regime

<sup>p</sup> energy relation which contrasts the linear in <sup>B</sup>

field B, Landau's quantization leads to a ffiffiffi

28 Graphene Materials - Structure, Properties and Modifications

[15] and cold atoms [16], just to mention a few.

means of ferromagnetic leads attached to the sample [18].

gapless edge states.

dependence of conventional 2DEG, with parabolic energy bands.

Before delving into the models of interest, let us present a summary of an important tool in the description of time-dependent Hamiltonian dynamics when the interaction term is periodic in time. These kinds of interactions are ubiquitous in physical systems ranging from cold atoms, cavity QED, superconducting interferometric devices, lasers, and so on. Let us then consider a generic time-dependent periodic Hamiltonian Hðt þ TÞ ¼ HðtÞ. We write it as free Hamiltonian plus a time-dependent interaction

$$H(t) = H\_0 + V(t). \tag{1}$$

The solution of the dynamics for the evolution operator

$$i\hbar \partial\_t U(t) = H(t)U(t),\tag{2}$$

is formally given by

$$\mathcal{U}(t) = \mathcal{T}\left[\exp\left(-\frac{\mathrm{i}}{\hbar}\int\_{0}^{t} H(t')dt'\right)\right],\tag{3}$$

with T the time-ordering operator. However, since the Hamiltonian HðtÞ is periodic, one can resort to Floquet theorem [24, 25]. It asserts that the general solution to the dynamics (2) when HðtÞ is periodic can be written as

$$
\mathcal{U}I(t) = P(t)e^{-iH\_\mathcal{F}t/\hbar},\tag{4}
$$

with PðtÞ periodic and HF a constant matrix, respectively. Using Eq. (4), it is easy to verify the stroboscopic property

$$\mathcal{U}(nT) = P(nT)e^{-i\pi T H\_{\mathcal{F}}/\hbar} \tag{5}$$

$$\mathcal{U}(nT) = \left[\mathcal{U}(T)\right]^n,\tag{6}$$

where <sup>U</sup>ðTÞ ¼ <sup>e</sup>�iHFT=<sup>ℏ</sup>. Thus, although the time-dependent Hamiltonian <sup>H</sup>ðt<sup>Þ</sup> is periodic, the corresponding evolution operator UðtÞ is not, that is, Uðt þ TÞ 6¼ UðtÞ. In fact, UðTÞ carries non-trivial information on the dynamics of the periodic system. The eigenvalues of HF give the quasi-energy spectrum of the periodically driven problem. In order to determine these quasienergies, one standard approach (see below) consists of performing an expansion in the (infinite) eigenbasis of time periodic functions <sup>ξ</sup>nðtÞ ¼ <sup>e</sup>inω<sup>t</sup> (Fourier modes), where <sup>ω</sup> <sup>¼</sup> <sup>2</sup>π=T. Then, in order to deal with the infinite eigenvalue problem one resorts to a truncation procedure in order to determine the Floquet exponents.

In the following, we would like to explicitly describe the Floquet-Fourier mode strategy. For this purpose, we assume that we can solve the dynamics of the free part H0. Its eigenbasis is spanned by the spinors jφα〉, where α describes a set of quantum numbers. We now use the eigenstates jφα〉 as expansion basis for the eigenstates of the full Hamiltonian in Eq. (1).

In order to analyse the evolution equation,

$$i\hbar \partial\_t |\Phi(t)\rangle = H(t)|\Phi(t)\rangle\tag{7}$$

we take advantage of the periodicity of the Hamiltonian so we can resort to Floquet's theorem [24, 25]. For this purpose, we define an auxiliary Hermitian Hamiltonian

$$\mathcal{H}(t) = H(t) - i\hbar \partial\_t,\tag{8}$$

along with the so-called Floquet states

$$|\Psi\_a(t)\rangle = \exp\left(i\varepsilon\_a t/\hbar\right)|\Phi(t)\rangle\tag{9}$$

such that

$$\mathcal{H}(t)|\Psi\_a(t)\rangle = \varepsilon\_a|\Psi\_a(t)\rangle\tag{10}$$

which are periodic functions of time, jΨαðt þ TÞ〉 ¼ jΨαðtÞ〉. In addition, the eigenvalues εα form the quasi-energy spectrum, and are the analogous of the quasi-momenta for Bloch electrons in a spatially periodic structure.

We can verify that the states

$$|\Psi\_{an}(t)|\rangle = \exp\left(imat)|\Psi\_a(t)\rangle\right.\tag{11}$$

are also eigenstates of the Hamiltonian HðtÞ but with corresponding eigenvalues εα ! εα þ nℏω. Thus, we can work in the first Brillouin zone �ℏω=2 ≤ εα ≤ ℏω=2.

Using the periodic temporal basis ξnðtÞ ¼ exp ðinωtÞ, which satisfies

$$\frac{1}{T} \int\_{0}^{T} \xi\_{n}^{\*}(t)\xi\_{m}(t) = \delta\_{nm},\tag{12}$$

we write the Fourier-mode expansion

$$\Phi\_a(t) = \exp\left(-i\varepsilon\_a t\right) \sum\_{n=-\infty}^{n=\infty} |\mathsf{C}\_a^{(n)}\rangle \xi\_n(t),\tag{13}$$

Now, we use the expansion <sup>j</sup>Cðn<sup>Þ</sup> <sup>α</sup> 〉 <sup>¼</sup> <sup>X</sup> β Λ<sup>ð</sup>n<sup>Þ</sup> αβ jφβ〉 such that Eq. (10) becomes

#### Photoinduced Pseudospin Dynamical Effects in Graphene-Like Systems http://dx.doi.org/10.5772/67618 31

$$\begin{split}H(t)\sum\_{n=-\infty}^{n=\infty}\sum\_{\beta}\Lambda\_{\alpha\beta}^{(n)}|\phi\_{\beta}\rangle\xi\_{n}(t)+\\\sum\_{n=-\infty}^{n=\infty}\sum\_{\beta}\Lambda\_{\alpha\beta}^{(n)}|\phi\_{\beta}\rangle\xi\_{n}(t)(\varepsilon\_{\alpha}-n\hbar\omega)=0.\end{split}\tag{14}$$

Multiplication by 〈φγjξ� <sup>m</sup>ðtÞ, then average over one temporal period, leads to

$$\sum\_{n=-\infty}^{n=\infty} \sum\_{\beta} [\langle a | H^{(m-n)} | \beta \rangle - (\varepsilon\_a - n\hbar \omega) \delta\_{nm} \delta\_{a\beta}] \Lambda\_{a\beta}^{(n)} = 0. \tag{15}$$

and we have used the simplifying notation <sup>j</sup>α〉 � jφα〉, and <sup>H</sup>ðm�n<sup>Þ</sup> <sup>¼</sup> <sup>1</sup>=<sup>T</sup> Z <sup>T</sup> 0 ξ� <sup>m</sup>ðtÞHðtÞξnðtÞ.

Then, the quasi-energies εα are eigenvalues of the secular equation.

$$\det|H\_F - \varepsilon\_a| = 0 \tag{16}$$

where 〈〈αmjHFjnβ〉〉 <sup>¼</sup> <sup>H</sup><sup>ð</sup>m�n<sup>Þ</sup> αβ þ mωδnmδαβ. Here, we simplify the notation to represent the direct product of orbital and periodic eigenfunctions jαm〉〉 � jα〉 ⊗ ξðtÞ.

For example, for a two-level problem we could take the interaction in the form VðtÞ∝σ ! � <sup>A</sup> ! ðtÞ. Then, if the vector potential given as A ! <sup>ð</sup>tÞ ¼ <sup>A</sup><sup>ð</sup> cos <sup>ω</sup>t; sin <sup>ω</sup>tÞ, we find that <sup>H</sup><sup>m</sup>�<sup>n</sup> has nonvanishing elements only for m � n ¼ 0; � 1 as it is described in the following section.

#### 2.2. Off-resonant approximation

quasi-energy spectrum of the periodically driven problem. In order to determine these quasienergies, one standard approach (see below) consists of performing an expansion in the (infinite) eigenbasis of time periodic functions <sup>ξ</sup>nðtÞ ¼ <sup>e</sup>inω<sup>t</sup> (Fourier modes), where <sup>ω</sup> <sup>¼</sup> <sup>2</sup>π=T. Then, in order to deal with the infinite eigenvalue problem one resorts to a truncation proce-

In the following, we would like to explicitly describe the Floquet-Fourier mode strategy. For this purpose, we assume that we can solve the dynamics of the free part H0. Its eigenbasis is spanned by the spinors jφα〉, where α describes a set of quantum numbers. We now use the eigenstates jφα〉 as expansion basis for the eigenstates of the full Hamiltonian in Eq. (1).

we take advantage of the periodicity of the Hamiltonian so we can resort to Floquet's theorem

which are periodic functions of time, jΨαðt þ TÞ〉 ¼ jΨαðtÞ〉. In addition, the eigenvalues εα form the quasi-energy spectrum, and are the analogous of the quasi-momenta for Bloch

are also eigenstates of the Hamiltonian HðtÞ but with corresponding eigenvalues

Xn¼∞ n¼�∞

<sup>j</sup>Cðn<sup>Þ</sup>

αβ jφβ〉 such that Eq. (10) becomes

εα ! εα þ nℏω. Thus, we can work in the first Brillouin zone �ℏω=2 ≤ εα ≤ ℏω=2.

ΦαðtÞ ¼ exp ð�iεαtÞ

β Λ<sup>ð</sup>n<sup>Þ</sup>

<sup>α</sup> 〉 <sup>¼</sup> <sup>X</sup>

Using the periodic temporal basis ξnðtÞ ¼ exp ðinωtÞ, which satisfies 1 T Z <sup>T</sup> 0 ξ�

[24, 25]. For this purpose, we define an auxiliary Hermitian Hamiltonian

iℏ∂tjΦðtÞ〉 ¼ HðtÞjΦðtÞ〉 (7)

HðtÞ ¼ HðtÞ � iℏ∂t; (8)

jΨαðtÞ〉 ¼ exp ðiεαt=ℏÞjΦðtÞ〉 (9)

HðtÞjΨαðtÞ〉 ¼ εαjΨαðtÞ〉 (10)

jΨα<sup>n</sup>ðtÞj〉 ¼ exp ðinωtÞjΨαðtÞ〉 (11)

<sup>n</sup>ðtÞξmðtÞ ¼ δnm; (12)

<sup>α</sup> 〉ξnðtÞ, (13)

dure in order to determine the Floquet exponents.

30 Graphene Materials - Structure, Properties and Modifications

In order to analyse the evolution equation,

along with the so-called Floquet states

electrons in a spatially periodic structure.

we write the Fourier-mode expansion

Now, we use the expansion <sup>j</sup>Cðn<sup>Þ</sup>

We can verify that the states

such that

In some physical scenarios, the frequency of the radiation field is way much larger than any other energy scale in the problem. Within this regime, one can derive an approximate Floquet Hamiltonian that captures the essence of the photoninduced bandgaps. In order to describe this so-called off-resonant regime, let us begin by explicitly showing the effective Floquet Hamiltonian: For this purpose, let us start from the general periodic Floquet Hamiltonian given in Eq. (1)

$$H(t) = H\_0 + V(t),\tag{17}$$

where H<sup>0</sup> is the static contribution and Vðt þ TÞ ¼ VðtÞ is the time-periodic interaction. Going to Fourier space, we get the general structure of the Floquet Hamiltonian in matrix form as

⋱ ⋯ V�2;�<sup>3</sup> H<sup>0</sup> � 2Ω V�2;�<sup>1</sup> V�2;<sup>0</sup> V�2;<sup>1</sup> V�2;<sup>2</sup> ⋯ ⋯ V�1;�<sup>3</sup> V�1;�<sup>2</sup> H<sup>0</sup> � Ω V�1;<sup>0</sup> V�1;<sup>1</sup> V�1;<sup>2</sup> ⋯ ⋯ V<sup>0</sup>;�<sup>3</sup> V<sup>0</sup>;�<sup>2</sup> V<sup>0</sup>;�<sup>1</sup> H<sup>0</sup> V<sup>0</sup>;<sup>1</sup> V<sup>0</sup>;<sup>2</sup> ⋯ ⋯ V<sup>1</sup>;�<sup>3</sup> V<sup>1</sup>;�<sup>2</sup> V<sup>1</sup>;�<sup>1</sup> V<sup>1</sup>;<sup>0</sup> H<sup>0</sup> þ Ω V<sup>1</sup>;<sup>2</sup> ⋯ ⋮⋮ ⋮ V<sup>2</sup>;<sup>0</sup> V<sup>2</sup>;<sup>1</sup> H<sup>0</sup> þ 2Ω ⋯ ⋱ 0 BBBBBBBB@ 1 CCCCCCCCA ; (18)

where the interaction submatrices are defined as

$$V\_{nn'} = \frac{1}{T} \int\_0^T dt V(t) e^{-i(n-n')\Omega t}. \tag{19}$$

For a monochromatic harmonic perturbation, this reduces to a block-tridiagonal matrix

⋱ ⋯ V�<sup>1</sup> H�<sup>2</sup> Vþ<sup>1</sup> 000 ⋯ ⋯ 0 V�<sup>1</sup> H�<sup>1</sup> Vþ<sup>1</sup> 0 0 ⋯ ⋯ 0 0 V�<sup>1</sup> H<sup>0</sup> Vþ<sup>1</sup> 0 ⋯ ⋯ 000 V�<sup>1</sup> H<sup>1</sup> Vþ<sup>1</sup> ⋯ ⋮⋮⋮ 0 V�<sup>1</sup> H<sup>2</sup> ⋯ ⋱ 0 BBBBBBBB@ 1 CCCCCCCCA ; (20)

and to simplify the notation we have set VN ¼ Vn�<sup>m</sup> and Hj ¼ H<sup>0</sup> þ jΩ. If we set out the eigenstate for a given number of Fourier modes n, we will have

$$\Phi = \begin{pmatrix} \phi\_{-n} \\ \phi\_{-n+1} \\ \vdots \\ \phi\_{-2} \\ \phi\_{-1} \\ \phi\_0 \\ \phi\_1 \\ \phi\_2 \\ \vdots \\ \phi\_{n-1} \\ \phi\_n \end{pmatrix}, \tag{21}$$

with each φ<sup>n</sup> a vector of dimensionality determined by H0. For instance, if we approximate the problem in such a way that we only consider one Fourier mode (n ¼ 1), we have to solve the following system of coupled equations:

$$\begin{aligned} H\_{-1}\phi\_{-1} + V\_{+1}\phi\_0 &= E\phi\_{-1} \\ V\_{-1}\phi\_{-1} + H\_0\phi\_0 + V\_{+1}\phi\_{+1} &= E\phi\_0 \\ H\_{+1}\phi\_{+1} + V\_{-1}\phi\_0 &= E\phi\_{+1} \end{aligned} \tag{22}$$

From the first and last equations, we get

$$\begin{aligned} \phi\_{-1} &= (E - H\_{-1})^{-1} V\_{+1} \phi\_0 \\ \phi\_{+1} &= (E - H\_{+1})^{-1} V\_{-1} \phi\_0, \end{aligned} \tag{23}$$

such that we get an effective equation for φ<sup>0</sup>

$$\left[\left(V\_{-1}(E - H\_{-1})^{-1}V\_{+1} + H\_0 + V\_{+1}(E - H\_{+1})^{-1}V\_{-1}\right)\phi\_0 = E\phi\_0\tag{24}$$

which explicitly reads

$$\left[\left(V\_{-1}(E - H\_0 + \Omega)\right)^{-1}V\_{+1} + H\_0 + V\_{+1}(E - H\_0 - \Omega)^{-1}V\_{-1}\right]\phi\_0 = E\phi\_0\tag{25}$$

For Ω≫jjH0jj, that is, frequencies much larger than the typical energy scales of the static problem, we can approximate this as

$$
\phi \left( H\_0 + \frac{V\_{-1}V\_{+1}}{\Omega} - \frac{V\_{+1}V\_{-1}}{\Omega} \right) \phi\_0 \approx E\phi\_0,\tag{26}
$$

so we get the effective approximate Floquet Hamiltonian, valid for large frequencies

$$H\_{\rm F} \approx H\_0 + \frac{[V\_{-1}, V\_{+1}]}{\Omega}.\tag{27}$$

With a similar procedure, one can show that for n ¼ 2 one gets the approximate Floquet Hamiltonian

$$H\_{\rm F} \approx H\_0 + \frac{[V\_{-1}, V\_{+1}]}{\Omega} - \frac{1}{2\Omega} \frac{[V\_{-1}^2, V\_{+1}^2]}{\Omega^2}.\tag{28}$$

These results essentially imply that a photoinduced energy bandgap can be induced by means of the dressed Floquet states that emerge from the off-resonant condition. Needless to say that sometimes the photoinduced energy bandgaps can be tiny; yet, in some circumstances one only needs to make sure that there is a gap, however, small and the topological properties of the driven system can be qualitatively different as those of the undriven (gapless) system.

#### 3. Photoinduced effects of Landau levels

Vnn<sup>0</sup> <sup>¼</sup> <sup>1</sup> T Z <sup>T</sup> 0

eigenstate for a given number of Fourier modes n, we will have

⋱

32 Graphene Materials - Structure, Properties and Modifications

0

BBBBBBBB@

following system of coupled equations:

From the first and last equations, we get

such that we get an effective equation for φ<sup>0</sup>

which explicitly reads

½V�<sup>1</sup>ðE � H�<sup>1</sup>Þ

�1

dtVðtÞe

and to simplify the notation we have set VN ¼ Vn�<sup>m</sup> and Hj ¼ H<sup>0</sup> þ jΩ. If we set out the

0

BBBBBBBBBBBBBBBB@

with each φ<sup>n</sup> a vector of dimensionality determined by H0. For instance, if we approximate the problem in such a way that we only consider one Fourier mode (n ¼ 1), we have to solve the

<sup>V</sup>�<sup>1</sup>φ�<sup>1</sup> <sup>þ</sup> <sup>H</sup>0φ<sup>0</sup> <sup>þ</sup> <sup>V</sup>þ<sup>1</sup>φþ<sup>1</sup> <sup>¼</sup> <sup>E</sup>φ<sup>0</sup>

Vþ<sup>1</sup> þ H<sup>0</sup> þ Vþ<sup>1</sup>ðE � Hþ<sup>1</sup>Þ

�1 Vþ<sup>1</sup>φ<sup>0</sup>

�1

�1

<sup>H</sup>�<sup>1</sup>φ�<sup>1</sup> <sup>þ</sup> <sup>V</sup>þ<sup>1</sup>φ<sup>0</sup> <sup>¼</sup> <sup>E</sup>φ�<sup>1</sup>

<sup>H</sup>þ<sup>1</sup>φþ<sup>1</sup> <sup>þ</sup> <sup>V</sup>�<sup>1</sup>φ<sup>0</sup> <sup>¼</sup> <sup>E</sup>φþ<sup>1</sup>:

<sup>φ</sup>�<sup>1</sup> ¼ ð<sup>E</sup> � <sup>H</sup>�<sup>1</sup><sup>Þ</sup>

<sup>φ</sup>þ<sup>1</sup> ¼ ð<sup>E</sup> � <sup>H</sup>þ<sup>1</sup><sup>Þ</sup>

<sup>φ</sup>�<sup>n</sup> <sup>φ</sup>�nþ<sup>1</sup> ⋮ <sup>φ</sup>�<sup>2</sup> <sup>φ</sup>�<sup>1</sup> φ0 φ1 φ2 ⋮ <sup>φ</sup><sup>n</sup>�<sup>1</sup> φn

1

CCCCCCCCCCCCCCCCA

For a monochromatic harmonic perturbation, this reduces to a block-tridiagonal matrix

⋯ V�<sup>1</sup> H�<sup>2</sup> Vþ<sup>1</sup> 000 ⋯ ⋯ 0 V�<sup>1</sup> H�<sup>1</sup> Vþ<sup>1</sup> 0 0 ⋯ ⋯ 0 0 V�<sup>1</sup> H<sup>0</sup> Vþ<sup>1</sup> 0 ⋯ ⋯ 000 V�<sup>1</sup> H<sup>1</sup> Vþ<sup>1</sup> ⋯ ⋮⋮⋮ 0 V�<sup>1</sup> H<sup>2</sup> ⋯

Φ ¼

�iðn�n<sup>0</sup> ÞΩt

⋱

: (19)

1

CCCCCCCCA

; (21)

<sup>V</sup>�<sup>1</sup>φ0; (23)

V�<sup>1</sup>�φ<sup>0</sup> ¼ Eφ<sup>0</sup> (24)

; (20)

(22)

In this section, we summarize our results reported in Ref. [23] where we theoretically analyse the dynamical manipulation of the LL structure of charge carriers on suspended monolayer graphene when a periodically driving radiation field is applied perpendicular to the sample. For this purpose, we focus on the low-energy properties of non-interacting spinless charge carriers in a suspended monolayer graphene subject to a perpendicular, uniform and constant magnetic field B ¼ B^z. The dynamics is governed by Dirac's Hamiltonian. In coordinate representation, it reads

$$\mathcal{H}\_{\eta}(\mathbf{r}) = \upsilon\_{\mathbf{F}}(\eta \pi\_{\mathbf{x}} \sigma\_{\mathbf{x}} + \pi\_{y} \sigma\_{y}),\tag{29}$$

where vF ∼ 106 m/s is the Fermi velocity in graphene. In addition, the canonical momenta π<sup>j</sup> ¼ pj þ eAj (j ¼ x; y) contain the vector potential (∇ · A ¼ B), �e is the electronic charge (e > 0), and η ¼ �1 describes the valley degree of freedom. Using the definition of the magnetic length l �2 <sup>B</sup> <sup>¼</sup> eB=<sup>ℏ</sup> and the cyclotron energy <sup>ℏ</sup>ω<sup>c</sup> <sup>¼</sup> ffiffiffi 2 <sup>p</sup> vFℏ=lB, the Hamiltonian Eq. (29) at each K (K') Dirac point, which corresponds to η ¼ þ1 (η ¼ �1), can be written in the form

$$H\_{+1} = \hbar \omega\_{\mathfrak{c}} \begin{pmatrix} 0 & a \\ a^{\sharp} & 0 \end{pmatrix} \tag{30}$$

$$H\_{-1} = -\hbar\omega\_c \begin{pmatrix} 0 & a^\dagger \\ a & 0 \end{pmatrix} \tag{31}$$

where the annihilation and creation operators are defined by standard relations as

$$a = \frac{l\_{\mathcal{B}}(\pi\_{\mathfrak{x}} - i\pi\_{\mathfrak{y}})}{\sqrt{2}} \qquad \text{and} \ a^{\dagger} = \frac{l\_{\mathcal{B}}(\pi\_{\mathfrak{x}} + i\pi\_{\mathfrak{y}})}{\sqrt{2}}.\tag{32}$$

The eigenenergies of the Hamiltonian (29) are then

$$E\_n^{s,\eta} = s\eta\sqrt{n}\hbar\omega\_c\tag{33}$$

with s ¼ �1. Positive (negative) values of sη represents the conduction (valence) band at each Dirac point. In addition, the integer quantum number n ¼ 0; 1; 2… labels the Landau-level structure of monolayer graphene. Using the eigenstates <sup>j</sup>n〉 of the operator <sup>a</sup>†a, the corresponding eigenstates <sup>j</sup>ϕ<sup>s</sup>;<sup>η</sup> <sup>n</sup> 〉 read

$$|\varphi\_n^{s,+1}\rangle = \frac{1}{\sqrt{2}}\binom{s|n-1\rangle}{|n\rangle}\tag{34}$$

$$|\varphi\_n^{s,-1}\rangle = \frac{1}{\sqrt{2}} \begin{pmatrix} -s|n\rangle\\|n-1\rangle \end{pmatrix} \tag{35}$$

for n 6¼ 0. The zero-energy eigenstate (n ¼ 0) is given in each case by

$$|\phi\_0^{+1}\rangle = \begin{pmatrix} 0\\|0\rangle \end{pmatrix} \tag{36}$$

$$
\langle \varphi\_0^{-1} \rangle = \begin{pmatrix} |0\rangle \\ 0 \end{pmatrix}. \tag{37}
$$

Due to time-reversal symmetry, we have T Hþ<sup>1</sup>T ¼ H�1. Let us now consider the effect of intense circularly polarized terahertz electromagnetic radiation, incident perpendicularly to the sample. We assume that the beam radiation spot is large enough compared to the lattice spacing so we can neglect any spatial variation. According to the standard light-matter interaction formulation, the dynamical effects of a monochromatic radiation field incident perpendicular to the sample can be described by means of a time-dependent vector potential

$$\mathbf{A}(t) = \frac{\mathcal{E}}{\omega} (\cos \omega t, \delta \sin \omega t), \tag{38}$$

where E and ω are, respectively, the amplitude and frequency of the electric field given in turn by the standard relation Eðt޼�∂tAðtÞ. In addition, we are using δ ¼ þ1 (δ ¼ �1) for right (left) circular polarization. We are using circular polarization because it has been shown to provide the maximal photoinduced bandgap. Starting from the ordinary dipolar interaction term �ep � AðtÞ, introduced to the tight-binding Hamiltonian via the Peierls substitution, we can evaluate the effects of the driving at each Dirac point as

$$V\_{\eta} = e\upsilon\_{\mathcal{F}}[\eta \sigma\_x A\_x(t) + \delta \sigma\_y A\_y(t)],\tag{39}$$

which explicitly reads

Hþ<sup>1</sup> ¼ ℏω<sup>c</sup>

H�<sup>1</sup> ¼ �ℏω<sup>c</sup>

where the annihilation and creation operators are defined by standard relations as

Es;<sup>η</sup>

<sup>n</sup> <sup>¼</sup> <sup>s</sup><sup>η</sup> ffiffiffi

with s ¼ �1. Positive (negative) values of sη represents the conduction (valence) band at each Dirac point. In addition, the integer quantum number n ¼ 0; 1; 2… labels the Landau-level structure of monolayer graphene. Using the eigenstates <sup>j</sup>n〉 of the operator <sup>a</sup>†a, the

> ffiffiffi 2 p

> > ffiffiffi 2 <sup>p</sup> �sjn〉 jn � 1〉 � �

<sup>0</sup> 〉 <sup>¼</sup> <sup>0</sup>

<sup>0</sup> 〉 <sup>¼</sup> <sup>j</sup>0〉 0 � �

Due to time-reversal symmetry, we have T Hþ<sup>1</sup>T ¼ H�1. Let us now consider the effect of intense circularly polarized terahertz electromagnetic radiation, incident perpendicularly to the sample. We assume that the beam radiation spot is large enough compared to the lattice spacing so we can neglect any spatial variation. According to the standard light-matter interaction formulation, the dynamical effects of a monochromatic radiation field incident perpen-

where E and ω are, respectively, the amplitude and frequency of the electric field given in turn by the standard relation Eðt޼�∂tAðtÞ. In addition, we are using δ ¼ þ1 (δ ¼ �1) for right (left) circular polarization. We are using circular polarization because it has been shown to

j0〉 � �

sjn � 1〉 jn〉 � �

<sup>a</sup> <sup>¼</sup> lBðπ<sup>x</sup> � <sup>i</sup>πy<sup>Þ</sup> ffiffiffi 2

<sup>n</sup> 〉 read

<sup>j</sup>ϕ<sup>s</sup>;þ<sup>1</sup> <sup>n</sup> 〉 <sup>¼</sup> <sup>1</sup>

<sup>j</sup>ϕ<sup>s</sup>;�<sup>1</sup> <sup>n</sup> 〉 <sup>¼</sup> <sup>1</sup>

<sup>j</sup>ϕþ<sup>1</sup>

<sup>j</sup>ϕ�<sup>1</sup>

dicular to the sample can be described by means of a time-dependent vector potential

<sup>A</sup>ðtÞ ¼ <sup>E</sup>

for n 6¼ 0. The zero-energy eigenstate (n ¼ 0) is given in each case by

The eigenenergies of the Hamiltonian (29) are then

34 Graphene Materials - Structure, Properties and Modifications

corresponding eigenstates <sup>j</sup>ϕ<sup>s</sup>;<sup>η</sup>

0 a a† 0 � �

0 a† a 0 � �

<sup>p</sup> and <sup>a</sup>† <sup>¼</sup> lBðπ<sup>x</sup> <sup>þ</sup> <sup>i</sup>πy<sup>Þ</sup> ffiffiffi

2

p : (32)

<sup>n</sup> <sup>p</sup> <sup>ℏ</sup>ω<sup>c</sup> (33)

: (37)

<sup>ω</sup> <sup>ð</sup> cos <sup>ω</sup>t; <sup>δ</sup> sin <sup>ω</sup>tÞ; (38)

(30)

(31)

(34)

(35)

(36)

$$V\_{\eta} = \xi \eta (\sigma\_x \cos \omega t + \eta \delta \sigma\_y \sin \omega t). \tag{40}$$

with the effective coupling constant ξ ¼ evFE=ω. This makes the total Hamiltonian

$$H\_{\eta}(t) = H\_{\eta} + V\_{\eta}(t),\tag{41}$$

periodic in time Hηðt þ TÞ ¼ HηðtÞ, with T ¼ 2π=ω the period of oscillation of the driving field. Therefore, if we focus on the K Dirac point (η ¼ 1), the physics at the K' Dirac point (η ¼ �1) can be easily found by the substitutions ξ ! �ξ and ω ! �ω.

Thus, let us focus on the K point physics and afterwards, we can make the necessary substitutions. In order to simplify the notation, we set Hþ<sup>1</sup> ¼ H<sup>0</sup> and Vþ<sup>1</sup>ðtÞ ¼ VðtÞ. Hence, defining rising σ<sup>þ</sup> and lowering σ� pseudospin operators by the standard formulas

$$
\sigma\_{\pm} = \frac{\sigma\_{\times} \pm i\sigma\_{y}}{2},
\tag{42}
$$

the time-dependent interaction potential can be rewritten as

$$V(t) = \xi(e^{-i\delta\omega t}\sigma\_+ + e^{i\delta\omega t}\sigma\_-),\tag{43}$$

Now, we invoke Floquet's theorem which states that the time evolution operator of the system induced by a periodic Hamiltonian can be written in the form [24]

$$
\mathcal{U}I(t) = P(t)e^{-iH\_\mathbb{F}t/\hbar},\tag{44}
$$

with PðtÞ a periodic unitary matrix and HF a time-independent dynamical generator referred to as the Floquet Hamiltonian. The eigenvalues of the Floquet Hamiltonian HF represent the quasi-energy spectrum of the periodically driven system. Typically, in order to solve for the quasi-energy spectrum, one can expand each term of the time-dependent Schrödinger equation in Fourier space and numerically solve an infinite eigenvalue problem. Instead, we will take a perturbative approach as discussed below.

Accordingly, for our problem we can find approximate solutions to the dynamics by modifying slightly the analytical strategy presented in Ref. [39]. Then, one finds that the excitation number operator, Na, defined as

$$N\_a = \left(a^\dagger a + \frac{1}{2}\right)\mathbb{I} + \frac{\sigma\_z}{2},\tag{45}$$

which commutes with the Hamiltonian H<sup>0</sup> and satisfies the eigenvalue equation

$$N\_a|\varphi\_n^s\rangle = n|\varphi\_n^s\rangle.\tag{46}$$

Na generates a time-dependent unitary transformation jΨðtÞ〉 ¼ PðtÞjΦðtÞ〉 given as

$$P(t) = \exp\left(-i\mathcal{N}\_d \delta \omega t\right),\tag{47}$$

such that the time-dependent Schrödinger equation

$$i\hbar \partial\_t |\Psi(t)\rangle = H(t)|\Psi(t)\rangle\tag{48}$$

can be transformed with a time-independent operator HF governing the dynamics of the problem

$$i\hbar \partial\_t |\Phi(t)\rangle = H\_F |\Phi(t)\rangle,\tag{49}$$

where HF and jΦðtÞ〉 are the Floquet Hamiltonian and Floquet eigenstate, respectively. Conducting the explicit calculation, HF is found to be given by

$$H\_{\rm F} = H\_0 - N\_a \delta \hbar \omega + \xi \sigma\_x,\tag{50}$$

where the effective coupling to the radiation field is ξ ¼ evFE=ω. This Hamiltonian can be treated via non-degenerate perturbation theory to find the approximate quasi-energy spectrum (see Ref. [23] for details on the derivation), which read

$$
\epsilon\_m^s = s \sqrt{m(\hbar \omega\_c)^2 + \left(m\xi\right)^2}, \text{ mod}\,\hbar\omega. \tag{51}
$$

The associated mean energies are in turn found as

$$
\overline{\epsilon}\_{m}^{s} = \epsilon\_{m}^{s} - \omega \frac{\partial \epsilon\_{m}^{s}}{\partial \omega},
\tag{52}
$$

which are invariant under e<sup>s</sup> <sup>m</sup> ! <sup>e</sup><sup>s</sup> <sup>m</sup> þ lℏω, for l being an integer. Conducting the explicit calculation, the mean energies are found to be given by the expression

$$\overline{\mathbf{e}}\_{m}^{s} = s \left( \frac{m(\hbar \omega\_{c})^{2} + 2m^{2} \xi^{2}}{\sqrt{m(\hbar \omega\_{c})^{2} + m^{2} \xi^{2}}} \right). \tag{53}$$

As can be seen in Figure 1, these mean energies are plotted as function of the quantizing magnetic field B, for different values of the Landau-level index changing the effective coupling ξ. We notice that, at intermediate light-coupling strength, the energy resolution of these levels becomes much better and could experimentally be tested for not so large

Figure 1. Approximate mean energies for the driven scenario (continuous lines) as function of the quantizing magnetic field B: The dotted lines represent the undriven spectra for the corresponding LL. We notice a level-dependent energy gapped that leads to non-trivial behaviour of physical quantities, as discussed below (see main text). Adapted from Ref. [23] with permission from APS.

quantizing magnetic fields B. Moreover, we find that to this order of approximation the LL becomes gapped, with the striking feature that the photoinduced gap is level-dependent. These gap openings appear except for the m ¼ 0 level which, as discussed before, remains insensitive to the radiation field.

#### 3.1. Discussion

Na <sup>¼</sup> <sup>a</sup>†

which commutes with the Hamiltonian H<sup>0</sup> and satisfies the eigenvalue equation

Najϕ<sup>s</sup>

Na generates a time-dependent unitary transformation jΨðtÞ〉 ¼ PðtÞjΦðtÞ〉 given as

such that the time-dependent Schrödinger equation

36 Graphene Materials - Structure, Properties and Modifications

Conducting the explicit calculation, HF is found to be given by

trum (see Ref. [23] for details on the derivation), which read

es <sup>m</sup> ¼ s

<sup>m</sup> ! <sup>e</sup><sup>s</sup>

e s <sup>m</sup> ¼ s

calculation, the mean energies are found to be given by the expression

The associated mean energies are in turn found as

which are invariant under e<sup>s</sup>

problem

a þ 1 2 � �

<sup>n</sup>〉 <sup>¼</sup> <sup>n</sup>jϕ<sup>s</sup>

can be transformed with a time-independent operator HF governing the dynamics of the

where HF and jΦðtÞ〉 are the Floquet Hamiltonian and Floquet eigenstate, respectively.

where the effective coupling to the radiation field is ξ ¼ evFE=ω. This Hamiltonian can be treated via non-degenerate perturbation theory to find the approximate quasi-energy spec-

ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi

<sup>m</sup> � ω

mðℏωcÞ

As can be seen in Figure 1, these mean energies are plotted as function of the quantizing magnetic field B, for different values of the Landau-level index changing the effective coupling ξ. We notice that, at intermediate light-coupling strength, the energy resolution of these levels becomes much better and could experimentally be tested for not so large

q

0 B@ <sup>2</sup> þ ðmξ<sup>Þ</sup> 2

> ∂ε<sup>s</sup> m

<sup>2</sup> <sup>þ</sup> <sup>2</sup>m<sup>2</sup>ξ<sup>2</sup> ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi mðℏωcÞ

<sup>2</sup> <sup>þ</sup> <sup>m</sup><sup>2</sup>ξ<sup>2</sup>

mðℏωcÞ

q

es <sup>m</sup> <sup>¼</sup> <sup>e</sup><sup>s</sup> 1 þ σz

<sup>2</sup> ; (45)

<sup>n</sup>〉: (46)

PðtÞ ¼ exp ð�iNaδωtÞ, (47)

iℏ∂tjΨðtÞ〉 ¼ HðtÞjΨðtÞ〉 (48)

iℏ∂tjΦðtÞ〉 ¼ HFjΦðtÞ〉; (49)

HF ¼ H<sup>0</sup> � Naδℏω þ ξσx; (50)

<sup>m</sup> þ lℏω, for l being an integer. Conducting the explicit

1

; modℏω: (51)

<sup>∂</sup><sup>ω</sup> ; (52)

CA: (53)

So far, we have shown that upon introduction of a perpendicularly radiation field in the Terahertz frequency domain, the Landau-level-quantized scenario can be turned into a leveldependent-gapped system. As shown in reference, this energy bandgap effects can be traced via the oscillations of the pseudospin polarization as well as the temporal evolution of the autocorrelation function for an initially prepared coherent superposition of the static Landau-level configuration. Thus, due to the natural connection between coherent-state superpositions and applications in quantum optics, one could expect that the experimental consequences of such radiation effects on the Landau-level structure of graphene would have some applications within experimentally accessible parameter regimes in the quantum optics realm.

#### 4. Irradiated silicene

As expected, the family of new two-dimensional materials with graphene-like properties has grown in the recent years. However, we show in this section that, although similar in lattice structure, both materials can be experimentally found to have different physical properties. In this section, we focus our attention on silicene which consists of a two-dimensional honeycomb lattice structure of silicon atoms analogous to that of graphene. Some works have reported the synthesis of silicene [30–32]. Silicene has a corrugated or buckled lattice structure that makes the silicon atoms in one sublattice to be perpendicularly displaced with respect to the other sublattice. For this reason, when a perpendicular electric field Ez is applied to silicene, the atoms belonging to each sublattice respond differently to Ez, giving rise to a staggered potential. Due to this peculiar pseudospin response to applied electric fields, and despite their similarities, the electronic properties of silicene are predicted to considerably differ from those of graphene.

In particular, since its intrinsic spin-orbit coupling is much larger than that of pristine graphene, an interesting interplay among intrinsic spin-orbit and electric field effects was predicted to appear because the bandgap can be electrically controlled. Moreover, the addition of an exchange potential term (which physically could represent the proximity effect due to coupling of ferromagnetic leads) allows for topological quantum phase transitions in the static regime [33]. Furthermore, in the presence of circularly polarized electromagnetic radiation, the realization of the so-called single Dirac cone phase in silicene has been recently proposed. At this topological phase, it is found that well-defined spin-polarized states are supported at every Dirac point. Moreover, within this configuration different spin components propagate in opposite directions giving rise to a pure spin current [34]. Yet, these photoinduced topological phase changes [20, 21] reported by Ezawa [34] were derived under the off-resonant assumption, that is, dynamical processes such that the frequency (coupling strength) of the radiation field is much larger (smaller) than any other energy scale in the problem. Under these assumptions, it is possible to derive an effective time-independent Floquet Hamiltonian [24, 25] with a tiny photoinduced bandgap correction that stems from virtual photon absorption and emission processes. Since the sign of the bandgap term (i.e., the effective bandgap) determines important topological properties of the material, it is vital both for potential practical implementations, for instance, in technological realizations of silicene-based devices, and from a fundamental point of view, to effectively achieve manipulation of this quantity.

In this section, we show via an exactly solvable model, where in order to detect physically relevant photoinduced effects in the energy band structure of silicene under strong circularly polarized electromagnetic radiation in the terahertz (frequency) domain one needs to go beyond the aforementioned off-resonant approximation. Indeed, we find that a zero momentum, the obtained exactly solvable time-dependent Hamiltonian, suggests the range of control parameters that physically might lead to experimentally feasible realization of new topological phases in silicene.

#### 4.1. Model

Let us consider the Dirac cone approximation to describe the dynamics of non-interacting charge carriers in silicene subject to a perpendicular, uniform and constant electric field E ¼ Ez^z. This is given by the 8 · 8 Hamiltonian [33] (ℏ ¼ e ¼ 1, with e being the electron's charge)

4. Irradiated silicene

38 Graphene Materials - Structure, Properties and Modifications

4.1. Model

As expected, the family of new two-dimensional materials with graphene-like properties has grown in the recent years. However, we show in this section that, although similar in lattice structure, both materials can be experimentally found to have different physical properties. In this section, we focus our attention on silicene which consists of a two-dimensional honeycomb lattice structure of silicon atoms analogous to that of graphene. Some works have reported the synthesis of silicene [30–32]. Silicene has a corrugated or buckled lattice structure that makes the silicon atoms in one sublattice to be perpendicularly displaced with respect to the other sublattice. For this reason, when a perpendicular electric field Ez is applied to silicene, the atoms belonging to each sublattice respond differently to Ez, giving rise to a staggered potential. Due to this peculiar pseudospin response to applied electric fields, and despite their similarities, the electronic properties of silicene are predicted to considerably differ from those of graphene.

In particular, since its intrinsic spin-orbit coupling is much larger than that of pristine graphene, an interesting interplay among intrinsic spin-orbit and electric field effects was predicted to appear because the bandgap can be electrically controlled. Moreover, the addition of an exchange potential term (which physically could represent the proximity effect due to coupling of ferromagnetic leads) allows for topological quantum phase transitions in the static regime [33]. Furthermore, in the presence of circularly polarized electromagnetic radiation, the realization of the so-called single Dirac cone phase in silicene has been recently proposed. At this topological phase, it is found that well-defined spin-polarized states are supported at every Dirac point. Moreover, within this configuration different spin components propagate in opposite directions giving rise to a pure spin current [34]. Yet, these photoinduced topological phase changes [20, 21] reported by Ezawa [34] were derived under the off-resonant assumption, that is, dynamical processes such that the frequency (coupling strength) of the radiation field is much larger (smaller) than any other energy scale in the problem. Under these assumptions, it is possible to derive an effective time-independent Floquet Hamiltonian [24, 25] with a tiny photoinduced bandgap correction that stems from virtual photon absorption and emission processes. Since the sign of the bandgap term (i.e., the effective bandgap) determines important topological properties of the material, it is vital both for potential practical implementations, for instance, in technological realizations of silicene-based devices, and from a fundamental point of view, to effectively achieve manipulation of this quantity.

In this section, we show via an exactly solvable model, where in order to detect physically relevant photoinduced effects in the energy band structure of silicene under strong circularly polarized electromagnetic radiation in the terahertz (frequency) domain one needs to go beyond the aforementioned off-resonant approximation. Indeed, we find that a zero momentum, the obtained exactly solvable time-dependent Hamiltonian, suggests the range of control parameters that physically might lead to experimentally feasible realization of new topological phases in silicene.

Let us consider the Dirac cone approximation to describe the dynamics of non-interacting charge carriers in silicene subject to a perpendicular, uniform and constant electric field

$$\begin{array}{l} \mathcal{H}^{\eta} = \upsilon\_{\mathcal{F}} (k\_{x}\sigma\_{x} + \eta k\_{y}\sigma\_{y}) + \sigma\_{z}(\eta \mathbf{s}\_{z}\lambda\_{s\circ} - l\mathbf{E}\_{z})\\ \qquad + \eta \sigma\_{z}h\_{11} + h\_{22} \end{array} \tag{54}$$

where vF <sup>¼</sup> ffiffi 3 <sup>p</sup> atb <sup>2</sup> <sup>≈</sup> <sup>8</sup>:<sup>1</sup> · 105 <sup>m</sup>=s is the Fermi velocity for charge carriers in silicene, with <sup>a</sup> <sup>¼</sup> <sup>3</sup>:<sup>86</sup> Å the lattice constant and tb ¼ 1:6 eV the hopping parameter within a tight-binding formulation, whereas l ¼ 0:23 Å measures half the separation among the two sublattice planes. In addition, η ¼ �1 describes the Dirac point, σ<sup>i</sup> and si (i ¼ x; y; z) are Pauli matrices describing pseudo- and real-spin degrees of freedom, respectively, whereas k ! ¼ ðkx; ηkyÞ is the momentum measured from the Dirac point η ¼ �1. The parameter λso ¼ 3:9 meV represents the strength of the intrinsic spin-orbit contribution. Moreover, the two contributions given by the terms

$$h\_{11} = a\lambda\_{R2}(k\_y s\_x - k\_x s\_y),\tag{55}$$

$$h\_{22} = \lambda\_{R1} (\eta \sigma\_x s\_y - \sigma\_y s\_x)/2,\tag{56}$$

describe the spin-orbit coupling associated to the next nearest neighbour hopping and nearest neighbour tight-binding formulation, respectively.

The term h<sup>11</sup> has its origin in the buckled structure of silicone, whereas h<sup>22</sup> is induced by the application of an external static electric field Ez. Using first-principle calculations, the authors of Ref. [40] found that λR<sup>2</sup> ¼ 0:7 meV and typically h<sup>22</sup> are of order μeV and thus much smaller than the other energy scales in the problem. Therefore, these two non-conserving contributions will be neglected in the following, although in the appendix we show that the largest contribution h<sup>11</sup> can be easily incorporated in the solution to the dynamical evolution presented below. Yet, we have verified that our results do not qualitatively change by the introduction of these two small corrections.

Within the approximation h<sup>22</sup> ¼ 0, let us now consider the pseudospin dynamics under an intense radiation field represented by the time-dependent vector potential

$$\mathbf{A}(t) = A(\cos \Omega t, \sin \Omega t), \tag{57}$$

with A ¼ E=Ω and Ω its amplitude and frequency, respectively. It describes a monochromatic electromagnetic wave incident perpendicular to the sample. This vector potential can in turn be derived from the corresponding electric field by means of Eðt޼�∂tAðtÞ, where E is the amplitude of the time-dependent electric field.

Using the standard minimal coupling prescription given as k ! ! k ! þ A ! , we get the dynamical generator

$$\mathcal{H}^{\eta}(\overrightarrow{k},t) = \upsilon\_{\text{F}}(k\_{\text{x}}\sigma\_{\text{x}} + \eta k\_{\text{y}}\sigma\_{\text{y}}) + \sigma\_{z}(\eta \text{s}\_{z}\Lambda\_{\text{so}} - l\text{E}\_{z})\tag{58}$$

$$\upsilon\_{\text{F}}A[e^{i\eta\Omega t}\sigma\_{-} + e^{-i\eta\Omega t}\sigma\_{+}].\tag{59}$$

In the following, we explore the emerging photoinduced dynamical features at zero momentum since this scenario allows for an exact analytical solution to the dynamical evolution equations. Given the fact that we have an exact analytical solution, we can explore the low-, intermediate-, and strong-coupling regimes of the charge carriers in silicene under the radiation field. We then discuss this exact solution and argue the need to explore either the intermediate or strong lightmatter-coupling regimes in order to obtain experimentally observable modifications in the physical properties of the system within the irradiation configuration.

#### 4.2. Physics at k = 0

In this subsection, we explicitly analyse another exactly solvable model for a graphene-like system. Let us then focus on the zero-momentum scenario for which the extrinsic spin-orbit term h<sup>11</sup> vanishes and the z-component of spin sz ¼ �1 is a good quantum number. Therefore, the following analysis is independent of taking into account the aforementioned spin-orbit contribution. Setting for notational convenience α ¼ vFA and Vz ¼ lEz, the physics at zero momentum k ! ¼ 0 is described by the dynamical generator

$$\mathcal{H}^{\eta}(0,t) = (\eta \text{s} \lambda\_{\text{so}} - V\_z)\sigma\_z + a[e^{i\eta \Omega t}\sigma\_- + e^{-i\eta \Omega t}\sigma\_+]. \tag{59}$$

Using the unitary transformation

$$\mathcal{P}^{\eta}(t) = \mathfrak{e}^{-i\eta(\mathbb{I} + \sigma\_{\overline{\mathbf{z}}})\Omega t/2} \tag{60}$$

we get the effective time-independent Floquet Hamiltonian <sup>H</sup>Fð<sup>k</sup> <sup>¼</sup> <sup>0</sup>Þ¼ðP<sup>η</sup><sup>Þ</sup> † <sup>ð</sup>tÞH<sup>η</sup>ð0;tÞP<sup>η</sup>ðt<sup>Þ</sup> �iðP† Þ η <sup>ð</sup>tÞP\_ <sup>η</sup> ðtÞ

$$\mathcal{H}\_{\rm F}(k=0) = -\frac{\eta \Omega}{2} \mathbb{I} + \left[ \eta \left( s\_z \lambda\_{\rm so} - \frac{\Omega}{2} \right) - V\_z \right] \sigma\_z + a \sigma\_x. \tag{61}$$

Thus, the zero-momentum quasi-energy spectrum is given as

$$
\varepsilon\_{\rm ss}^{\eta}(\mathbf{k}=\mathbf{0}) = -\frac{\eta\Omega}{2} + \sigma\sqrt{a^2 + \left(\Lambda\_s^{\eta}\right)^2},\tag{62}
$$

where s; σ ¼ �1 represent the real- and pseudospin degrees of freedom, respectively. In addition, we have defined the effective gap

$$
\Delta\_s^\eta = \eta \left( s \lambda\_{\text{so}} - \frac{\Omega}{2} \right) - V\_z. \tag{63}
$$

On the other hand, the zero-momentum exact Floquet eigenstates are

$$|\psi\_{s\sigma}^{\eta}(t)\rangle = \frac{e^{-ie\_{s\eta}^{\eta}t}}{\sqrt{2\Gamma}} \begin{pmatrix} e^{-i\eta\Omega t} \sqrt{\Gamma + \sigma\Delta\_{s}^{\eta}} \\ \sigma\sqrt{\Gamma - \sigma\Delta\_{s}^{\eta}} \end{pmatrix},\tag{64}$$

with Γ ¼ ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi <sup>α</sup><sup>2</sup> þ ðΔ<sup>η</sup> s Þ 2 q . Some comments are in order at this point. The exact quasi-energy spectrum resembles the solution for the Rabi problem and one could expect that Rabi oscillations should appear in the dynamical evolution of this zero-momentum solution. Moreover, since we are interested in analysing the behaviour of a topological quantity not in terms of the momentum variable but assuming as a toy model that one explores the Vz parameter space along with the 'torus' for

the time domain 0 <sup>≤</sup> <sup>t</sup> <sup>≤</sup> <sup>T</sup>, one could evaluate the components of the pseudospin operator ^ S ! α defined within each subband as

$$\dot{\overrightarrow{\mathbf{S}}}^{\eta} = \langle \psi\_{s\sigma}^{\eta}(t) | \overrightarrow{\mathbf{s}} | \psi\_{s\sigma}^{\eta}(t) \rangle \tag{65}$$

After some algebra, we get

$$\hat{\vec{S}}^{\eta} = \eta (-\sin\beta\sin\Omega t, \sin\beta\cos\Omega t, -\cos\beta) \tag{66}$$

where

In the following, we explore the emerging photoinduced dynamical features at zero momentum since this scenario allows for an exact analytical solution to the dynamical evolution equations. Given the fact that we have an exact analytical solution, we can explore the low-, intermediate-, and strong-coupling regimes of the charge carriers in silicene under the radiation field. We then discuss this exact solution and argue the need to explore either the intermediate or strong lightmatter-coupling regimes in order to obtain experimentally observable modifications in the phys-

In this subsection, we explicitly analyse another exactly solvable model for a graphene-like system. Let us then focus on the zero-momentum scenario for which the extrinsic spin-orbit term h<sup>11</sup> vanishes and the z-component of spin sz ¼ �1 is a good quantum number. Therefore, the following analysis is independent of taking into account the aforementioned spin-orbit contribution. Setting for notational convenience α ¼ vFA and Vz ¼ lEz, the physics at zero

iηΩt

2

� Vz

ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi <sup>α</sup><sup>2</sup> þ ðΔ<sup>η</sup> s Þ 2

� �

q

� �

σ� þ e

�iηΩt

�iηð1þσ<sup>z</sup> <sup>Þ</sup>Ωt=<sup>2</sup> (60)

σþ�: (59)

†

σ<sup>z</sup> þ ασx: (61)

; (62)

; (64)

� Vz: (63)

<sup>ð</sup>tÞH<sup>η</sup>ð0;tÞP<sup>η</sup>ðt<sup>Þ</sup>

ical properties of the system within the irradiation configuration.

40 Graphene Materials - Structure, Properties and Modifications

¼ 0 is described by the dynamical generator

<sup>H</sup><sup>η</sup>ð0;tÞ¼ðηsλso � VzÞσ<sup>z</sup> <sup>þ</sup> <sup>α</sup>½<sup>e</sup>

<sup>P</sup><sup>η</sup>ðtÞ ¼ <sup>e</sup>

<sup>2</sup> <sup>1</sup> <sup>þ</sup> <sup>η</sup> szλso � <sup>Ω</sup>

<sup>2</sup> <sup>þ</sup> <sup>σ</sup>

where s; σ ¼ �1 represent the real- and pseudospin degrees of freedom, respectively. In addi-

� �

2

e�iηΩ<sup>t</sup> ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi <sup>Γ</sup> <sup>þ</sup> <sup>σ</sup>Δ<sup>η</sup> s

p σ ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi <sup>Γ</sup> � <sup>σ</sup>Δ<sup>η</sup> s

p � �

<sup>s</sup> <sup>¼</sup> <sup>η</sup> <sup>s</sup>λso � <sup>Ω</sup>

η sσt ffiffiffiffiffi <sup>2</sup><sup>Γ</sup> <sup>p</sup>

we get the effective time-independent Floquet Hamiltonian <sup>H</sup>Fð<sup>k</sup> <sup>¼</sup> <sup>0</sup>Þ¼ðP<sup>η</sup><sup>Þ</sup>

<sup>s</sup>σð<sup>k</sup> <sup>¼</sup> <sup>0</sup>޼� ηΩ

Δη

On the other hand, the zero-momentum exact Floquet eigenstates are

<sup>s</sup>σðtÞ〉 <sup>¼</sup> <sup>e</sup>�i<sup>ε</sup>

jψη

<sup>H</sup>Fð<sup>k</sup> <sup>¼</sup> <sup>0</sup>޼� ηΩ

Thus, the zero-momentum quasi-energy spectrum is given as

εη

4.2. Physics at k = 0

momentum k

�iðP† Þ η <sup>ð</sup>tÞP\_ <sup>η</sup> ðtÞ

with Γ ¼

!

Using the unitary transformation

tion, we have defined the effective gap

ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi <sup>α</sup><sup>2</sup> þ ðΔ<sup>η</sup> s Þ 2

.

q

$$
\tan \beta = \frac{2\alpha}{|\Delta\_s^\eta|}. \tag{67}
$$

Thus, within this limit, the average out-of-plane pseudospin polarization S z <sup>η</sup> is the only finite component and it reads

$$\overline{S}^{\overline{z}}\_{\eta} = -\frac{\eta |\Delta^{\eta}\_{s}|}{\Gamma} \tag{68}$$

Thus, the main advantage of finding semi-analytical solutions to the dynamical evolution is that closed expressions for the physical quantities of interest can be found in such a way that one can provide further insight into the nature of the physical mechanisms involved and how their interplay leads to a given behaviour of the polarization, charge current and so on.

Let us now assume that the system is initially prepared in the arbitrary state

$$|\Phi(0)\rangle = \begin{pmatrix} \cos\frac{\theta}{2}e^{i\phi/2} \\ \sin\frac{\theta}{2}e^{-i\phi/2} \end{pmatrix},\tag{69}$$

with 0 ≤ θ ≤ π and 0 < φ ≤ 2π being spherical coordinates over the Bloch sphere. Thus, the evolution of the pseudospin polarization is given by the standard relation <sup>σ</sup>zðtÞ ¼ 〈Φð0ÞjU† <sup>F</sup>ðtÞ <sup>σ</sup>zUFðtÞjΦð0Þ〉, with UFðt<sup>Þ</sup> being the unitary Floquet evolution operator UFðtÞ ¼ <sup>P</sup><sup>η</sup>ðtÞe�iHFt (note that <sup>σ</sup><sup>z</sup> and <sup>P</sup><sup>η</sup>ðt<sup>Þ</sup> commute with each other). The initial polarization in the state (65) is given by σzð0Þ ¼ cos θ. After some algebra, we find

$$\begin{split} \sigma\_z(t) &= \frac{2\alpha}{\Gamma} \sin \theta \sin \Gamma t \Big( \frac{\Delta\_s^\eta}{\Gamma} \sin \Gamma t \cos \phi - \cos \Gamma t \sin \phi \Big) \\ &\quad + \cos \theta \Big( 1 - \frac{2a^2}{\Gamma^2} \sin^2 \Gamma t \Big). \end{split} \tag{70}$$

In addition, the one-period mean-value pseudospin polarization

$$
\langle \sigma\_z \rangle = \frac{1}{T} \int\_0^T \sigma\_z(t) dt,\tag{71}
$$

with T ¼ 2π=Ω being the period of oscillations of the driving field, gives then

$$\begin{split} \langle \sigma\_z \rangle &= \alpha \sin \theta \Big[ \frac{\Lambda\_s^\eta}{\Gamma^2} \cos \phi (1 - \text{sinc}(2\Gamma T)) - T \sin \phi \text{sinc}^2(\Gamma T) \Big] \\ &+ \cos \theta \Big[ 1 - \frac{\alpha^2}{\Gamma^2} (1 - \text{sinc}(2\Gamma T)) \Big], \end{split} \tag{72}$$

where sincðxÞ ¼ sin <sup>ð</sup>x<sup>Þ</sup> <sup>x</sup> .

In particular, for initial states that have zero polarization (θ ¼ π=2), we get the simplified expressions

Figure 2. Zero-momentum pseudospin mean polarization for different combinations of the product sη. The vertical lines correspond to α ¼ 0:1Ω (continuous), α ¼ 0:25Ω (large dashed) α ¼ 0:5Ω (short dashed) and α ¼ 0:75Ω (large dots), respectively. See discussion in the main text. Adapted from Ref. [35] with permission from APS.

$$\langle \sigma\_z \rangle = a \left[ \frac{\Delta\_s^\eta}{\Gamma^2} \cos \phi (1 - \text{sinc}(2\Gamma T)) - T \sin \phi \text{sinc}^2(\Gamma T) \right]. \tag{73}$$

Setting the value φ ¼ π=4, we plot in Figure 2 the mean pseudospin polarization for the different spin and valley sη product combinations.

From this figure, we find that within the low coupling regime (α ≤ 0:1Ω), it is in general not possible to induce appreciable changes of the pseudospin polarization and this is related to the fact that the quasi-energy behaviour is essentially controlled by the parameters Vz and λso which determine the gap behaviour in the static regime. On the other hand, for intermediate (α ¼ 0:5Ω) and large (α ¼ 0:75Ω) values of the coupling to the driving field, that is, beyond the off-resonant condition, effective pseudospin inversion is achievable and therefore a qualitatively different behaviour emerges within this coupling regime.

#### 5. Conclusions

〈σz〉 <sup>¼</sup> <sup>1</sup> T Z <sup>T</sup> 0

<sup>Γ</sup><sup>2</sup> cos <sup>φ</sup>ð<sup>1</sup> � sincð2ΓTÞÞ � <sup>T</sup> sin <sup>φ</sup>sinc2

In particular, for initial states that have zero polarization (θ ¼ π=2), we get the simplified

Figure 2. Zero-momentum pseudospin mean polarization for different combinations of the product sη. The vertical lines correspond to α ¼ 0:1Ω (continuous), α ¼ 0:25Ω (large dashed) α ¼ 0:5Ω (short dashed) and α ¼ 0:75Ω (large dots),

respectively. See discussion in the main text. Adapted from Ref. [35] with permission from APS.

<sup>Γ</sup><sup>2</sup> <sup>ð</sup><sup>1</sup> � sincð2ΓTÞÞ � �

� �

with T ¼ 2π=Ω being the period of oscillations of the driving field, gives then

<sup>þ</sup> cos <sup>θ</sup> <sup>1</sup> � <sup>α</sup><sup>2</sup>

s

〈σz〉 <sup>¼</sup> <sup>α</sup> sin <sup>θ</sup> <sup>Δ</sup><sup>η</sup>

42 Graphene Materials - Structure, Properties and Modifications

where sincðxÞ ¼ sin <sup>ð</sup>x<sup>Þ</sup>

expressions

<sup>x</sup> .

σzðtÞdt; (71)

ðΓTÞ

(72)

;

In this chapter, we have described some photoinduced consequences of using periodically driven interactions in graphene-like systems. In essence, these examples provide new insights on the non-trivial behaviour of systems taken out of equilibrium and given the fact that Floquet's theorem allows for a dynamical analysis in terms of an equivalent static description of the physics via the Floquet states, one can infer new physical aspects that can emerge on these systems when subject to such periodic interactions. We consider that these approaches do provide an arena for analysing some interesting phenomena beyond the static limit. The extension of the results presents more involved scenarios, where numerical tools can profit from the simple models described in this chapter and we hope the interested reader could profit from the material presented. One could further explore the dynamical features at finite momentum but refer the reader to Ref. [35] where it is proven that the intermediate lightmatter-coupling regime is the physically correct scenario for describing the emergent topological features such as the realization of the single-valley polarized state in silicene. Needless to say that other quantities of interest such as charge and spin currents can be treated within the approximate scenarios discussed in the chapter, they go beyond the main interest of the chapter and are thus not further discussed here.

#### Appendix

For completeness, in this appendix we include one approximation strategy that can be useful in the semi-analytical treatment of periodically driven systems. This is called the rotating wave approximation concerns with the dynamical evolution of a periodically driven system when the parameters are closed to a resonance. These resonant effects have been shown to be relevant in several experimental situations such as within cavity QED as well as within the context of cold atoms, among others.

#### Approximate RWA dynamics

In some physical systems, the most important physical effects happen around the first resonance located at around ω<sup>0</sup> ¼ ω, even a low coupling ξ ¼ 0:25ω. This motivates the search of an analytical solution by means of the rotating wave approximation (RWA) approximation. We now find an approximate analytical solution to the dynamics generated by

$$H(t) = \frac{\omega\_0}{2}\sigma\_z + A\sin\omega t \sigma\_x. \tag{74}$$

Then, in the low coupling regime where the amplitude of the driving field satisfies <sup>A</sup> <sup>ω</sup> ≪ 1, we can use the rotating wave approximation (RWA) which amounts to keep those contributions in the driving that are close to resonance, that is, ω~ω0.

Changing to the interaction representation, the dynamics is dictated by VIðtÞ ¼ <sup>e</sup><sup>i</sup>ω0tσz=<sup>2</sup>Vðt<sup>Þ</sup> <sup>e</sup>�iω0tσz<sup>=</sup>2, where <sup>V</sup>ðt<sup>Þ</sup> is the time-dependent contribution to Eq. (74). Using the algebra of the Pauli matrices and the explicit time dependence of the driving field, we get

$$V\_I(t) = A \sin \omega t (\cos \alpha\_0 t \sigma\_{\tilde{x}} - \sin \alpha\_0 t \sigma\_{\tilde{y}}) \tag{75}$$

At this point, we invoke the RWA to get the approximate interaction contribution

$$V\_{RWA}(k,t) = i\frac{A}{2}(e^{-i\mu t}\sigma\_+ - e^{i\mu t}\sigma\_-),\tag{76}$$

where the frequency detuning that characterizes near-resonance contributions is defined as μ ¼ ω � ω0. To get Eq. (76) from Eq. (75), we have assumed ω � ω<sup>0</sup> ~0 and thus have disregarded the quickly oscillating terms proportional to ω þ ω<sup>0</sup> (which are known in the literature as the secular contributions).

Switching back to the Schrödinger representation, we get the total RWA Hamiltonian

$$H\_{RWA}(t) = \frac{\omega\_0}{2}\sigma\_z + i\lg(e^{-i\omega t}\sigma\_+ - e^{i\omega t}\sigma\_-),\tag{77}$$

where we have introduced the effective coupling constant g ¼ <sup>A</sup> <sup>2</sup> The Hamiltonian (77) is taking to a time-independent form HF by means of a periodic unitary transformation <sup>P</sup>ðtÞ ¼ <sup>e</sup>�iðσ0þσ<sup>z</sup> <sup>Þ</sup>ωt=<sup>2</sup> by means of the standard transformation rule HF <sup>¼</sup> <sup>P</sup>† <sup>ð</sup>tÞHRWAPðtÞ � iP† ðtÞ∂tPðtÞ. Explicitly, we have

$$H\_F = -\frac{\omega}{2}\sigma\_0 - \frac{\mu}{2}\sigma\_z - g\,\sigma\_y. \tag{78}$$

Defining the Rabi frequency <sup>Ω</sup><sup>R</sup> <sup>¼</sup> ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi <sup>μ</sup><sup>2</sup> <sup>þ</sup> <sup>4</sup> <sup>g</sup><sup>2</sup> <sup>p</sup> , the approximate quasi-energies are found to be given as

Photoinduced Pseudospin Dynamical Effects in Graphene-Like Systems http://dx.doi.org/10.5772/67618 45

$$
\varepsilon\_{\alpha} = -\frac{(\omega + \alpha \Omega\_{\mathbb{R}})}{2} \tag{79}
$$

where α ¼ �1 is the energy subband index. The corresponding approximate Floquet eigenstates are

$$|\Phi\_a(k,t)\rangle = \frac{e^{-ie\_at}}{\sqrt{2\Omega\_{\mathbb{R}}}} \begin{pmatrix} e^{-i\omega t} \sqrt{\Omega\_{\mathbb{R}} - \alpha\mu} \\ i\alpha \sqrt{\Omega\_{\mathbb{R}} + \alpha\mu} \end{pmatrix}. \tag{80}$$

As before, once the Floquet quasi-energies and eigenstates are found, one can evaluate any physical quantity of interest. For instance, one can evaluate the spin (or any other related quantity that satisfies the Pauli matrices algebra) polarization. In addition, the evaluation of the charge as well as spin currents can be given analytical expressions within the RWA. The relevance of resonant phenomena in nature has already been highlighted in the introduction, but we would like to remember that within this scheme, photon-assisted transport is one of the many resonant processes that can be properly described via the RWA. However, we must mention that one of the drawbacks of the RWA is that it does not take into account a systematic shift of the resonant energies which occurs at higher coupling strengths to the driving field. These are the so-called Bloch-Siegert shifts which emerge from the neglected secular terms. Thus, in order to go to higher coupling strengths and preserve a more realistic physical picture, this secular term should be included. In this context, the Van Vleck quasi-degenerate perturbation theory as well as Magnus expansion approach can be useful semi-analytical tools and we refer the interested reader to Ref. [41] for a detailed treatment of the Magnus expansion approach and to Ref. [42] for a thorough description of the Van Vleck approach.

#### Author details

Approximate RWA dynamics

44 Graphene Materials - Structure, Properties and Modifications

by

have

given as

In some physical systems, the most important physical effects happen around the first resonance located at around ω<sup>0</sup> ¼ ω, even a low coupling ξ ¼ 0:25ω. This motivates the search of an analytical solution by means of the rotating wave approximation (RWA) approximation. We now find an approximate analytical solution to the dynamics generated

<sup>2</sup> <sup>σ</sup><sup>z</sup> <sup>þ</sup> <sup>A</sup> sin <sup>ω</sup>tσx: (74)

σ�Þ, (76)

σ�Þ, (77)

<sup>2</sup> The Hamiltonian (77) is taking to

<sup>2</sup> <sup>σ</sup><sup>z</sup> � <sup>g</sup> <sup>σ</sup>y: (78)

ðtÞ∂tPðtÞ. Explicitly, we

VIðtÞ ¼ A sin ωtð cos ω0tσ<sup>x</sup> � sin ω0tσyÞ (75)

<sup>ω</sup> ≪ 1, we

<sup>H</sup>ðtÞ ¼ <sup>ω</sup><sup>0</sup>

Pauli matrices and the explicit time dependence of the driving field, we get

VRWAðk; tÞ ¼ i

HRWAðtÞ ¼ <sup>ω</sup><sup>0</sup>

HF ¼ � <sup>ω</sup>

where we have introduced the effective coupling constant g ¼ <sup>A</sup>

by means of the standard transformation rule HF <sup>¼</sup> <sup>P</sup>†

Defining the Rabi frequency <sup>Ω</sup><sup>R</sup> <sup>¼</sup> ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi

At this point, we invoke the RWA to get the approximate interaction contribution

A 2 ðe �iμt

Switching back to the Schrödinger representation, we get the total RWA Hamiltonian

<sup>2</sup> <sup>σ</sup><sup>z</sup> <sup>þ</sup> igð<sup>e</sup>

a time-independent form HF by means of a periodic unitary transformation <sup>P</sup>ðtÞ ¼ <sup>e</sup>�iðσ0þσ<sup>z</sup> <sup>Þ</sup>ωt=<sup>2</sup>

<sup>2</sup> <sup>σ</sup><sup>0</sup> � <sup>μ</sup>

�iωt

σ<sup>þ</sup> � e iωt

<sup>ð</sup>tÞHRWAPðtÞ � iP†

<sup>μ</sup><sup>2</sup> <sup>þ</sup> <sup>4</sup> <sup>g</sup><sup>2</sup> <sup>p</sup> , the approximate quasi-energies are found to be

where the frequency detuning that characterizes near-resonance contributions is defined as μ ¼ ω � ω0. To get Eq. (76) from Eq. (75), we have assumed ω � ω<sup>0</sup> ~0 and thus have disregarded the quickly oscillating terms proportional to ω þ ω<sup>0</sup> (which are known in the

σ<sup>þ</sup> � e iμt

the driving that are close to resonance, that is, ω~ω0.

literature as the secular contributions).

Then, in the low coupling regime where the amplitude of the driving field satisfies <sup>A</sup>

can use the rotating wave approximation (RWA) which amounts to keep those contributions in

Changing to the interaction representation, the dynamics is dictated by VIðtÞ ¼ <sup>e</sup><sup>i</sup>ω0tσz=<sup>2</sup>Vðt<sup>Þ</sup> <sup>e</sup>�iω0tσz<sup>=</sup>2, where <sup>V</sup>ðt<sup>Þ</sup> is the time-dependent contribution to Eq. (74). Using the algebra of the

> Alexander López<sup>1</sup> \* and Benjamin Santos<sup>2</sup>

\*Address all correspondence to: alexander.lopez@physik.uni-regensburg.de

1 Escuela Superior Politécnica del Litoral, ESPOL, Facultad de Ingeniería Mecánica y Ciencias de la Producción, Campus Gustavo Galindo Km 30.5 Vía Perimetral, Guayaquil, Ecuador

2 INRS-EMT, Université du Québec, Varennes, Québec, Canada

#### References


[3] Geim A. K. and Novoselov K. S., Nat. Mater. 6, 183 (2007).

[6] Kane C. L. and Mele E. J., Phys. Rev. Lett. 95, 226801 (2005), ibid, Phys. Rev. Lett. 95,

[7] Huertas-Hernando D., Guinea F., and Brataas A., Phys. Rev. B 74, 155426 (2006).

[10] Ando T., Matsumoto Y. and Uemura Y., J. Phys. Soc. Jpn. 39, 279 (1975). [11] Klitzing K. v, Dorda G. and Pepper M., Phys. Rev. Lett. 45, 494 (1980).

[13] Bernevig B. A., Hughes T. L. and Zhang S. C., Science 314, 1757 (2006).

[14] Fu L., Kane C. L. and Mele E. J., Phys. Rev. Lett. 98, 106803 (2007).

[19] Lindner N. H., Refael G. and Galinski V., Nat. Phys. 7, 490 (2011).

[22] Cayssol J., Dóra B., Simon F., and Moessner R., Phys. Status Solidi RRL 7, 101 (2013).

[29] Guzmán-Verri G. G. and Lew Yan Voon L. C., Phys. Rev. B 76, 075131 (2007).

[23] López A., Di Teodoro A., Schliemann J., Berche B., Santos B., Phys. Rev. B 92, 235411

[12] Hasan M. Z. and Kane C. L., Rev. Mod. Phys. 82, 3045 (2010).

[15] Goldman N. et al, Phys. Rev. Lett. 105, 255302 (2010).

[16] Schneider U. et al., Nat. Phys. 8, 213 (2012).

[18] Yang Y. et al, Phys. Rev. Lett. 107, 066602 (2011)

[20] Kitagawa T. et al., Phys. Rev. B 82, 235114 (2010). [21] Kitagawa T. et al., Phys. Rev. B 84, 235114 (2011).

[24] Grifoni M. and Hänggi P., Phys. Rep. 304, 229 (1998). [25] Chu S.-I. and Telnov D. A., Phys. Rep. 390, 1 (2004).

[27] Krueckl V. and Kramer T., New J. Phys. 11, 093010 (2009). [28] Takeda K. and Shiraishi K., Phys. Rev. B 50, 14916 (1994).

[30] Vogt P. et al., Phys. Rev. Lett. 108, 155501 (2012).

[26] Shirley J. H. Phys. Rev. 138, B979 (1965).

[17] König M. et al., Science 318, 766 (2008).

(2015).

[4] Castro Neto A. H. et al., Rev. Mod. Phys. 81, 109 (2009). [5] Das Sarma S. et al., Cond-mat.mes-hall, 1003.4701 (2010).

[8] Rashba E., Phys. Rev. B 79, 161409(R) (2009).

46 Graphene Materials - Structure, Properties and Modifications

[9] Ertler C. et al., Phys. Rev. B 80, 041405(R) (2009).

146802 (2005).

