Nonlinear electromagnetic response of few-layer graphene: A nonperturbative description

Nonperturbative approach based on exact solution of Boltzmann kinetic equation in the relaxation time approximation is developed for the study of nonlinear response of electron-doped few-layer graphene to a high-frequency electromagnetic field. It is shown that nonperturbative approach can be applied to a two-dimensional conductor with an arbitrary isotropic spectrum of carries. The cases of ABC stacked three-layer and ABC staked four-layer graphene are considered. The low-energy electron spectrum of such graphenes is characterized by the third and the fourth power in the momentum dependence, correspondingly. The transmission, reflection and absorption coefficients are calculated. It is shown that the transmission coefficient T of three- and four-layer graphene decreases under increase in the intensity of the incident wave I by the asymptotic law T I   , respectively. It is found that three-layer and four-layer graphenes demonstrate power-induced reflectance, in contrast to power-induced transparency in monolayer graphene. It is found that in the four-layer graphene irradiated by a monochromatic wave the induced electrical current contains only the first and the third harmonics. The analytical expression for the efficiency of the third-harmonic generation G TGH (the ratio of the intensity of third-harmonic radiation to the intensity of the incident wave) of the four-layer graphene is obtained. It is shown that only at rather small incident wave intensity the efficiency G TGH is proportional to the second power of I and at large incident wave intensity the efficiency G TGH approaches the constant quantity. Saturation of the efficiency of the third-order harmonic generation is caused by the increase of the reflection. In contrast, the efficiency of the third-order harmonic generation of monolayer graphene depends nonmonotonically on the intensity of the incident wave. The maximum is reached at rather small intensity I аnd the efficiency in the maximum is of order of 10 -4 that is two orders of magnitude smaller

Ключевые слова: графен, возрастание отражательной способности в мощном поле, генерация третьей гармоники Introduction Linear dependence of electron spectrum on the momentum causes many unusual properties of graphene. One of them is a strong nonlinear electromagnetic response. Nonlinearity of the response can be seen from the dependence of the velocity on the momentum of charged quasiparticles with the spectrum ( where n is the density of carriers. In [1,2] nonlinear electromagnetic properties of graphene were studied within the quasiclassical approach based on the Boltzmann kinetic equation. The density matrix approach was developed in [3][4][5][6]. It was predicted in [3][4][5][6] that the third-harmonic generation (THG) intensity in electron-doped graphene has the main peak at were obtained within the diagrammatic approach [7].
Strong enhancement of THG in a system of two graphenes, one of which is the electron-doped, and the other, the hole-doped, was predicted in [8]. The approaches used in [1][2][3][4][5][6][7][8] are perturbative ones. In [9] a nonperturbative theory of nonlinear electromagnetic response of graphene was developed. The theory is based on the exact solution of the kinetic Boltzmann equation within the relaxation-time approximation. One of interesting results of [9] is the absence of optical bistability in graphene predicted in [10] but not confirmed in [9]. Note that in [10] the exact solution of the kinetic Boltzmann equation was also used but the answer was expanded in series in the amplitude of the electric field. The difference between the semiclassical [1,2,9,10] and quantum [3][4][5][6][7][8] approaches is that the former ones take into account only the electron band. The graphene spectrum near Dirac points contains the electron and hole bands that touch each other in Dirac points. A nonperturbative approach that takes into account two bands was developed in [11,12]. The results of [12] are based on a heuristic solution of the time-dependent Dirac equation for the two-component wave function (graphene Bloch equation). It was found in [12] that the contribution of the hole band into the nonlinear response of electrondoped graphene is essential only at F   . While the approach [9] is not valid at high frequencies (for instance, it cannot reproduce peaks in THG intensity [3][4][5][6][7][8]), at low frequencies it gives adequate description. In particular, in [9] the power-induced transparency in graphene and third-, fifth-and seventh-harmonic generation at large input power were described quantitatively. The electron spectrum of bilayer graphene differs from one of monolayer graphene. The bilayer graphene has AB stacking of layers and the spectrum of the electron band can be approximated by free electron spectrum , where m is the effective mass. Therefore at F  and low temperature the electron-doped bilayer graphene should not demonstrate any nonlinear electromagnetic response. In the general case nonlinear response of bilayer graphene can be comparable with one of monolayer graphene [8,13]. The spectrum of few-layer graphene depends on stacking of layers. Quite interesting situation is realized in ABC stacked N-layer graphenes. Their low-energy electron spectrum can be approximated as [14,15] 1 0 cm/s is the Fermi velocity in the monolayer graphene, and 0.4 t   eV is the nearestneighbor interlayer hopping energy. Note that Eq. (1) gives the spectrum of the monolayer and bilayer graphene as well. One can expect that the spectrum (1) with 2 N  reveals itself in an unusual nonlinear electromagnetic response. In this paper we study nonlinear response of graphene with the spectrum (1) using the semiclassical nonperpurbative approach. The main attention is given to the four-layer graphene ( 4 N  ) for which an analytical expression for the nonlinear part of the electrical current can be obtained. For 4 N  we calculate the dependence of the transmission, reflection and adsorption coefficients and the efficiency of the third-harmonic generation on the intensity of the incident wave. The results are compared with ones for monolayer graphene.

Nonperturbative expression for the electrical current
We start from the Boltzmann kinetic equation in the relaxation time approximation where 00 ( , ) The electric field is directed along the x -axis.
The exact solution of Eq. (2) with the initial condition 00 ( , , ) | x y t f p p t f   has the form [16] .
The electrical current induced by the field () Et is calculated as   (7) in the system with the spectrum (1) is proportional to 1 0 N E  . In the dirty limit it is equal to Up to the factor ( 1)! N  Eq. (9) coincides with the

Transmission, reflection and absorption of highintensity incident wave in a four-layer graphene
At even N one can obtain from Eq. (7) the explicit dependence of the current on the electric field. In particular, for N=4 The current (10) can be presented as a sum of a linear and a cubic in 0 E terms: j j j . The linear term is equal to ( The obtained dependences (11) and (12) allows to calculate the transmission, reflection and absorption coefficients for a monochromatic wave as a function of the incident intensity. We consider the normal incidence. We are interested in a frequency range in which the wavelength is much larger than the thickness of four-layer graphene. In this case the graphene can be treated as a zero-thickness boundary between the upper and the lower half-spaces. We take the electric field of the transmitted wave in the form ( ) sin( ),

Third-harmonic generation in four-layer graphene
Nonlinear response causes generation of harmonics. In systems with inversion symmetry only odd harmonics are generated. The harmonics are generated by the component the electrical current that oscillates with the corresponding frequency. Usually the intensity of the third harmonic is the largest one. To calculate the efficiency of THG (the ratio of the intensity of the third harmonic (3) I to the input intensity inc I ) we should obtain the equation for the electric field of the wave emitted at the frequency 3 . This field is taken in the form (3) [ sin(3 3 ) where the upper/lower sign corresponds to the field in the upper(lower) half-space. The boundary conditions yield , One can see that even at very large intensity of the incident wave relative losses caused by THG do not exceed 0.02 in the pure limit and 0.2 in the dirty limit (we take into account the emission into the upper and the lower half-spaces).
It is instructive to compare the efficiency of THG in monolayer and four-layer graphene. Let us do it in the dirty limit. In this limit the current (7)  In Fig.1 we present the dependence of the transmission and reflection coefficients on the intensity of the incident wave for four-layer graphene in the pure limit calculated for the parameters given above.
In Fig.2 the transmission, reflection and absorption coefficients in a dirty four-layer graphene as functions of the intensity of the incident waves are shown. Note that in physical units (MW/cm 2 ) the ranges of input intensities in Fig. 1 and Fig. 2 are almost the same. One can see that the main difference between the pure and dirty limits is nonzero absorption in the latter case. The absorption coefficient depends non-monotonically on the intensity of the incident wave and decreases at large input power. The transmission decreases and the reflection increases at large input power in the pure and in the dirty limits. In Fig. 3 the efficiency of THG in pure and dirty fourlayer graphenes is presented. One can see that quadratic dependence of the efficiency on the intensity of the incident wave survives only at very small intensities. For the parameters considered the maximum efficiency does not exceed few percents. At the same input power the efficiency of THG in a dirty graphene is larger than in a pure one. The efficiency of THG is the increase function of the intensity of the incident wave.
To compare the behavior of four-layer and monolayer graphene we calculate the transmission, reflection, and absorption coefficients in a dirty monolayer graphene for the same  and for the same density of carries per layer.
The result is displayed in Fig. 4. Note that for the parameters specified the reference intensity I  for monolayer graphene is in four times smaller than for four-layer graphene. In physical units the range of intensities in Fig. 4 is in two times smaller than in Figs. 1-3. Figs. 2 and 4 illustrate the difference in the nonlinear response of few-layer and monolayer graphene. The former one demonstrates the power-induced reflectance while the latter one, the power-induced transparency. We emphasize that this difference emerges if the frequency of the incident wave satisfies the condition / F   that corresponds to the frequency range up to several terahertz.
The efficiency of THG of monolayer graphene is shown in Fig. 5. One can see that it is in two orders smaller than one for the four-layer graphene (Fig. 3b). Another difference is that the efficiency Fig. 5 is a nonmonotonic function of the intensity of the incident wave with the maximum at rather small intensities. The differences are connected with that the amplitude of the third harmonic of the electrical current (24) saturates at large electric field, while the amplitude of the third harmonic of the current (12) is proportional to the third power of the electric field. The saturation of THG efficiency in Fig. 3 is solely due to decrease of the transmission T. Fig. 5. The efficiency of third-harmonic generation in the dirty limit in monolayer graphene versus the intensity of the incident wave. Fig. 4. Transmission, reflection and absorption coefficients in monolayer graphene in the dirty limit versus the intensity of the incident wave. Fig. 3. The efficiency of third-harmonic generation in the pure (a) and the dirty (b) limit in four-layer graphene versus the intensity of the incident wave. Fig. 2. Transmission, reflection and absorption coefficients in four-layer graphene in the dirty limit versus the intensity of the incident wave.
One should note that in this study we neglect mixing of the first harmonic with generated third harmonic. This effect results in fifth-harmonic generation, а reduction of reflection and lowering of the efficiency of thirdharmonic generation at very large input intensities. In the range of I inc considered two latter effects are small, less than one percent, and an error in dependences presented in Figs. 1-5 does not exceed 0.01.

Conclusion
In conclusion, we have shown that nonperturbative theory of nonlinear electromagnetic response developed for monolayer graphene can be generalized to the system with an arbitraty isotropic spectrum of carriers. The nonperturbative approach is applied to the study of nonlinear electromagnetic properties of few-layer graphenes with ABC stacking of layers. It is established that nonlinear behavior of electron-doped N-layer graphenes with N>2 differs significantly from one of monolayer graphene. We predict that the transmission coefficient T of three-and four-layer graphene decreases under increase in the intensity of the incident wave inc I by the asymptotic law contrast to the case of monolayer graphene where this quantity reaches its maximum at certain inc I and then decreases under further increase of inc I . The difference is connected with saturation of the amplitude of third harmonic of an electrical current induced in monolayer graphene by a monochromatic field and the absence of such saturation in few-layer graphene. The efficiency of third-harmonic generation in a four-layer graphene can reach several percents which is two orders in magnitude larger than the maximum efficiency in monolayer graphene.