## Second-harmonic generation in subwavelength graphene waveguides

Daria Smirnova and Yuri S. Kivshar

Nonlinear Physics Center, Australian National University, Canberra ACT 0200, Australia

We suggest a novel approach for generating second-harmonic radiation in subwavelength graphene waveguides. We demonstrate that quadratic phase matching between the plasmonic guided modes of different symmetries can be achieved in a planar double-layer geometry when conductivity of one of the layers becomes spatially modulated. We predict theoretically that, owing to graphene nonlocal conductivity, the second-order nonlinear processes can be actualized for interacting plasmonic modes with an effective grating coupler to allow external pumping of the structure and output of the radiation at the double frequency.

PACS numbers: 78.67.Wj, 42.65.Wi, 73.25.+i, 78.68.+m

## I. INTRODUCTION

Graphene plasmonics is the rapidly developing field of nanophotonics that attracted a lot of attention during the past two years 1-7 . One of the recent novel directions in this field is the study of nonlinear phenomena given the fact that strong confinement of surface p-plasmons in graphene can enhance light-matter interaction 8 and facilitate nonlinear response . 9 In particular, the recent theoretical studies included the analysis of spatial plasmon solitons in graphene layers 10,11 , nonlinear difference frequency generation of THz surface plasmons 12 , and singlephoton operation in a graphene cavity . 9

Double-layer graphene waveguides are known to support symmetric and antisymmetric plasmonic guiding modes classified with respect to the in-plane electric field profile 13-15 , so we may expect the coupling between different modes when the nonlinear response becomes important. In this paper, we suggest and study analytically a novel approach for the second-harmonic generation in double-layer graphene waveguides taking the advantage of the possibility to modulate spatially the conductivity of one of the graphene layers by doping or electrostatic gating 16-18 thus assisting the coupling of freely propagating light to plasmons in graphene 19-22 .

Phased matching between parametrically interacting waves is known to be a crucial requirement for the efficient second-order nonlinear effects 23 , and in our geometry the phase matching can be achieved by tailoring the modal dispersion of the graphene waveguide enabling plasmon-to-plasmon frequency conversion. More specifically, we demonstrate below that the phase matching becomes possible between an antisymmetric fundamental frequency (FF) mode and a symmetric second-harmonic (SH) mode in a double-layer graphene waveguide.

We employ a theoretical analysis based on the perturbation theory and describe the second-order nonlinear process with graphene plasmons in terms of slowly varying modes of a waveguide structure, the approximation frequently used in the nonlinear optics 24 . For definiteness, we consider a planar geometry shown schematically in Fig. 1, where a graphene double-layer waveguide is placed into homogeneous surrounding medium with the dielectric permittivity ε being illuminated by light from the upper half-space x &gt; d/ 2. As a matter of fact, in our analysis we distinguish two parts of the problem and examine them sequentially,

- · Linear scattering. Radiation at the fundamental frequency normally incident onto the structure is scattered by the periodic conductivity grating, and then excites resonantly an antisymmetric mode of the graphene waveguide;
- · Nonlinear plasmon-to-plasmon conversion. Due to the second-order nonlocal nonlinearity of graphene 25-27 , the induced current of the antisymmetric mode serves as a source for the phasematched symmetric mode at the double frequency which eventually radiates into free space.

Similar to the assumptions employed earlier 11,15,28 , namely low dissipation, weak nonlinearity, and small phase mismatch, here we solve Maxwell's equations following the procedure of the asymptotic expansion, and demonstrate that the second-order nonlinear processes

Figure 1. (Color online) Geometry of the problem. A freestanding waveguide created by two graphene layers is illuminated by an external wave being coupled to the phasematched guided modes of the graphene waveguide. Red curve corresponds to an antisymmetric mode at the fundamental frequency, and a blue curve - to a second-harmonic symmetric mode (shown with the tangential electric field component). Conductivity of the upper layer is periodically modulated.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0555v1-figure-946a13ea-9c55-436c-a5fb-77690de8cdd2.jpg)

can be achieved with conductivity modulation to allow external pumping to couple to the guided modes of the structure generating output radiation at the doubled frequency.

## II. LINEAR SCATTERING

We start our derivation by studying the linear resonant excitation of antisymmetric plasmons, and write the corresponding system of Maxwell's equations in the form

$$\begin{array} {$$

where j (1 2) , are the surface current densities induced in the graphene layers placed at x = ± d/ 2, as indicated by Dirac's delta functions δ . Assuming the field to be p-polarized with the magnetic component H = H x,z,t ( ) y 0 , from Eq. (1) we find

$$& \quad \vdots \psi, \varepsilon, \nu / \Psi$$

where, if we assume the harmonic field dependencies ∼ exp( -iωt ), the currents are given by

$$j ^ { ( 1, 2 ) } ( \omega, z ) & = \sigma ^$$

Here σ (2) ( ω, z ) = σ ω ( ), σ (1) ( ω, z ) = σ ω ( )[1 + f ( z )] are surface conductivities, σ ω ( ) ≡ σ ( R ) ( ω ) + iσ ( I ) ( ω ) is the linear frequency-dependent surface conductivity of graphene, while f ( z ) = f ( z + a ) = f 1 cos ( 2 π a z ) +

f 2 cos ( 4 π a z ) is assumed to be a periodic function of z , i.e. conductivity of one of the graphene layers is spatially modulated.

Assuming the dissipative losses and modulation amplitudes f 1 2 , to be small, we formally introduce the smallness parameter µ

$$\text{never } \mu \\ \mu = \max \left \{ \left | \frac$$

∣ ∣ and adopt the following asymptotic ansatz for the magnetic field

$$\text{and approve one moving any purpose across or one may} \\ \text{etc field}$$

where h x ( ) is the transverse profile of the linear plasmonic mode of the guiding structure, A 1 2 , are slowly varying mode amplitudes, where the subscripts '1' and '2' refer to the forward and backward propagation along the z direction, respectively. Also ˜ H 1 2 , are small corrections to the eigenmode profile, and the term ¯ H 0 does not contain a fast dependence on z .

In the zero-order approximation in µ , for the function h x ( ) Eq. (2) takes the form,

$$\text{nant} \quad & \text{ant} \quad & \frac { d ^ { 2 } h } { d x ^ { 2 } } - ( k _ { \text{sp} } ^ { 2 } - k _ { 0 } ^ { 2 } \varepsilon ) h ( x ) = \\ & \quad = \frac { 4 \pi } { c } i \sigma ^ { ( I ) } ( \omega ) \frac { \partial } { \partial x } \left \{ \left [ \delta \left ( x - \frac { d } { 2 } \right ) + \delta \left ( x + \frac { d } { 2 } \right ) \right ] e _ { z } ( x ) \right \}, \\ \left | \right | \text{$z_{\infty}$} \,. \quad \text{$\dots$} \,. \quad \text{$i$} \,. \quad \text{$\dots$} \,. \quad \text{$\dots$} \,. \quad \text{$\dots$} \,. \quad \text{$\dots$} \,.$$

where e z ( x ) = i k ε dx 0 dh . This equation yields the dispersion relation

$$\left [ 1 - \frac { 2 \pi \sigma ^ { ( I ) } \kappa ^ { ( s, a ) } } { \omega \varepsilon } ( 1 \pm e ^ { - \kappa ^ { ( s, a ) } d } ) \right ] = 0 \,, \quad ( 7 )$$

where κ (s,a) = √ k (s,a)2 sp -k ε 2 0 , k 0 = ω/c , and the modal transverse profiles are

$$h ^ { ( s, \mathbf a ) } ( \omega, x ) = \frac { - i k _ { 0 } \varepsilon } { \kappa ^ { ( s, \mathbf a ) 2 } } \frac { d e _ { z } ^ { ( s, \mathbf a ) } } { d x } \,, \quad \quad ( 8 )$$

where the continuous tangential components of the electric field are expressed as and

$$\begin{array} { c } \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \end{array} \ e _ { z } ^ { ( s ) } ( \omega, x ) = \begin{cases} e ^ { - \kappa ^ { ( s ) } ( x - d / 2 ) }, & x > d / 2 \,, \\ \cosh ( \kappa ^ { ( s ) } x ) & \{ \cosh ( \kappa ^ { ( s ) } d / 2 ) \\ e ^ { \kappa ^ { ( s ) } ( x + d / 2 ) }, & x < - d / 2 \end{cases} \\ \end{array} \ e ^ { ( s ) } ( x + d / 2 ), & x < - d / 2 \end{cases}$$

$$e _ { z } ^ { ( a ) } ( \omega, x ) = \begin{cases} e ^ { - \kappa ^ { ( a ) } ( x - d / 2 ) }, \quad x > d / 2 \,, \\ \sinh ( \kappa ^ { ( a ) } x ) \\ \sinh ( \kappa ^ { ( a ) } d / 2 ) \\ - e ^ { \kappa ^ { ( a ) } ( x + d / 2 ) }, \quad x < - d / 2 \end{cases} \\ \end{cases} \\ \text{for both symmetric and antisymmetric eigenmodes} \begin{cases} e ^ { - \kappa ^ { ( a ) } ( x - d / 2 ) }, \quad x > d / 2 \,, \\ \sinh ( \kappa ^ { ( a ) } x ) \\ \sinh ( \kappa ^ { ( a ) } d / 2 ) \\ - e ^ { \kappa ^ { ( a ) } ( x + d / 2 ) }, \quad x < - d / 2 \end{cases}$$

for both symmetric and antisymmetric eigenmodes guided by a double-layer waveguide 14,29 , as indicated by the superscripts.

We now assume that the period of the conductivity grating is chosen to almost compensate the momentum mismatch with the antisymmetric plasmon, | k (a) sp ( ω ) -2 π/a | ≡ | ∆ k | /lessmuch 2 π/a so that | ∆ k/k (a) sp ( ω ) | /lessorsimilar µ . Substituting the expression of the total field (5) into Eq. (2) and keeping only resonant terms, in the first order

in µ for the field perturbations ˜ H 1 2 , , we obtain

$$& \text{in } \mu \text{ for the field perturbations } \tilde { H } _ { 1, 2 }, \text{ we obtain } \quad & \text{varyin} \\ & \frac { d ^ { 2 } \tilde { H } _ { 1, 2 } } { d x ^ { 2 } } - ( k ^ { ( a ) 2 } _ { \text{sp} } ( \omega ) - k ^ { 2 } _ { 0 ^ { 2 } \varepsilon ) \tilde { H } _ { 1, 2 } - \frac { 4 \pi } { c } i \sigma ^ { ( I ) } ( \omega ) } \quad \pm 2 \\ & \times \frac { \partial } { \partial x } \left \{ \left [ \delta \left ( x - \frac { d } { 2 } \right ) + \delta \left ( x + \frac { d } { 2 } \right ) \right ] \tilde { E } _ { 1, 2 z } \right \} = F _ { 1, 2 } \,, \\ & F _ { 1, 2 } ( x, z ) = \mp 2 i k ^ { ( a ) } _ { \text{sp} } ( \omega ) \frac { \partial \mathcal { A } _ { 1, 2 } } { \partial z } h ^ { ( a ) } ( \omega, x ) + \frac { 4 \pi } { c } \sigma ^ { ( R ) } ( \omega ) \mathcal { A } _ { 1, 2 } } \quad \text{where} \\ & \times \frac { \partial } { \partial x } \left \{ \left [ \delta \left ( x - \frac { d } { 2 } \right ) + \delta \left ( x + \frac { d } { 2 } \right ) \right ] e ^ { ( a ) } _ { z } ( \omega, x ) \right \} } \quad \text{gratin} _ { \text{pling } } \quad \text{form} \\ & + \frac { 4 \pi } { c } i \sigma ^ { ( I ) } ( \omega ) \frac { \partial } { \partial x } \left \{ \delta \left ( x - \frac { d } { 2 } \right ) \left [ \frac { f _ { 1 } } { 2 } \tilde { E } _ { 0 z } e ^ { \mp i \Delta k z } } \\ & + \frac { f _ { 2 } } { 2 } \mathcal { A } _ { 2, 1 } \, e ^ { \mp 2 i \Delta k z } e ^ { ( a ) } _ { z } ( \omega, x ) \right ] \right \}, \\ & \text{whereas } \tilde { H } _ { 0 } \text{ satisfies the equation,} \\ & \quad \varepsilon \quad \cdots \quad \, \quad \, \wedge \, \quad \, \cdots \quad \, \quad \, \cdots \quad \, \quad \, \cdots \quad \,.$$

whereas ¯ H 0 satisfies the equation,

$$\Lambda _ { 0 } ^ { 2 } \bar { H } _ { 0 } & _ { 0 } ^ { 2 } \varepsilon \bar { H } _ { 0 } = \frac { 4 \pi } { c } i \sigma ^ { ( I ) } ( \omega ) \frac { \partial } { \partial x } \Big \{ \Big [ \delta \Big ( x + \frac { d } { 2 } \Big ) + \delta \Big ( x - \frac { d } { 2 } \Big ) \Big ] \quad \text{where in} \\ \times E _ { 0 z } ( \omega, x ) + \frac { f _ { 1 } } { 2 } \Big ( \mathcal { A } _ { 1 } e ^ { i \Delta k z } + \mathcal { A } _ { 2 } e ^ { - i \Delta k z } \Big ) \delta \Big ( x - \frac { d } { 2 } \Big ) \Big \} \,. \quad \text{inition}, \\ \text{Equation} \ (12) \ \text{includes reflection of a plane wave in} \quad q ^ { ( \alpha ) }$$

Equation (12) includes reflection of a plane wave incident from the semi-space x &gt; d/ 2, written as H 0 exp[ - √ ik 0 √ ε x ( -d/ 2)], where we assume H 0 ∼ µ , E 0 = H / 0 ε , E 0 z ( ω, x ) is the tangential electric field distribution at f 1 2 , = 0. Taking into account the boundary conditions at the graphene interface, we find the amplitude of the electric field at x = d/ 2 as follows

$$& \bar { E } _ { 0 z } ( \omega, x = d / 2 ) = E _ { 0 } ( 1 + r ) \\ & - \frac { 2 \pi } { c \sqrt { \varepsilon } } i \sigma ^ { ( I ) } ( \omega ) \frac { f _ { 1 } } { 2 } \left ( \mathcal { A } _ { 1 } e ^ { i \Delta k z } + \mathcal { A } _ { 2 } e ^ { - i \Delta k z } \right ) \,. \quad \text{(13)} \quad \text{ and } s \\ & \quad. \quad \text{$\text{the a\choose \epsilonf\colon$} }.$$

$$& \text{where } \ r \ = \ - \frac { 4 \pi } { c \sqrt { \varepsilon } } i \sigma ^ { ( I ) } ( \omega ) \left ( 1 + \frac { 4 \pi } { c \sqrt { \varepsilon } } i \sigma ^ { ( I ) } ( \omega ) \right ) ^ { - 1 } \ \approx \quad & \text{coffic} \\ & - \frac { 4 \pi } { } i \sigma ^ { ( I ) } ( \omega ) \text{ stands for the reflection coefficient from} \quad & \text{where}$$

two closely spaced graphene layers, k 0 εd /lessmuch 1.

- √ 4 π c ε iσ ( I ) ( ω ) stands for the reflection coefficient from √

Substituting the expression (13) into Eq. (11), in order for the corrections ˜ H 1 2 , to be nondiverging, in accord with the Fredholm alternative 30 , we have to satisfy the orthogonality condition for the right-hand side of Eq. (11) with the solution of its homogeneous part, or equivalently, with the plasmonic mode itself. This overintegrating is mathematically written as

$$\underset { - \infty } { \overset { + \infty } { \int } } F _ { 1, 2 } ( x, z ) h ^ { ( \mathbf a ) } ^ { * } ( \omega, x ) d x = 0, \quad \ \ ( 1 4 ) \quad \text{ the} \\ - \infty$$

and it leads to the nonlinear equations for the slowly

.

varying amplitudes A 1 2 ,

$$& \cdots _ { \nu } \cdots _ { \circ } \cdots _ { r } \cdots \cdots \cdots _ { a, e } \\ & \quad \pm 2 i k _ { s p } ^ { ( a ) } ( \omega ) \left ( \frac { \partial \mathcal { A } _ { 1, 2 } } { \partial z } \pm ( \gamma ^ { ( a ) } + \gamma _ { r } ^ { ( a ) } ) \mathcal { A } _ { 1, 2 } \right ) \\ = & ( - 2 i k _ { s p } ^ { ( a ) } ( \omega ) \gamma _ { r } ^ { ( a ) } + Q _ { 2 } ^ { ( a ) } ) \mathcal { A } _ { 2, 1 } e ^ { \mp 2 i \Delta k z } + Q _ { 1 } ^ { ( a ) } E _ { 0 } ( 1 + r ) e ^ { \mp i \Delta k z }, \\. &.$$

where linear damping due to losses in graphene γ (a) , grating-induced damping due to radiation γ (a) r , and coupling coefficients Q (a) 1 2 , are derived, respectively, in the form

$$\text{m} \\ \gamma ^ { ( \mathbf a ) } = \frac { 4 \pi } { c } \sigma ^ { ( R ) } ( \omega ) \frac { k _ { 0 } \varepsilon } { k _ { \text{sp} } ^ { ( \mathbf a ) } ( \omega ) q ^ { ( \mathbf a ) } } \,, \\ \gamma ^ { ( \mathbf a ) } _ { r } = \left ( \frac { \pi \sigma ^ { ( I ) } ( \omega ) f _ { 1 } } { c } \right ) ^ { 2 } \frac { k _ { 0 } \varepsilon ^ { 1 / 2 } } { k _ { \text{sp} } ^ { ( \mathbf a ) } ( \omega ) q ^ { ( \mathbf a ) } } \,, \\ Q ^ { ( \mathbf a ) } _ { 1, 2 } = \frac { 2 \pi } { c } \sigma ^ { ( I ) } ( \omega ) f _ { 1, 2 } \frac { k _ { 0 } \varepsilon } { q ^ { ( \mathbf a ) } } \,, \\ \text{row in all the software in a major two for all the time of}$$

where in all the expressions we employ the following definition,

$$\int \cdot \quad \cdot \quad \cdot \quad \cdot q ^ { ( a ) } = k _ { 0 } ^ { 2 } \frac { \varepsilon ^ { 2 } } { \kappa ^ { ( a ) 3 } } \left [ 1 + \frac { \sinh \left ( \kappa ^ { ( a ) } d \right ) + \kappa ^ { ( a ) } d } { 2 \sinh ^ { 2 } \left ( \kappa ^ { ( a ) } d / 2 \right ) } \right ]. \quad \ \ ( 1 7 )$$

When we assume a simple homogeneous case, ∂ ∂z = 0, and also the exact momentum matching, ∆ k = 0, from Eq. (15) we obtain the slowly varying amplitudes,

$$\mathcal { A } _ { 1 } = \mathcal { A } _ { 2 } \equiv \mathcal { A } = \frac { Q _ { 1 } ^ { ( \mathbf a ) } E _ { 0 } ( 1 + r ) } { 2 i k _ { \text{sp} } ^ { ( \mathbf a ) } ( \omega ) [ \gamma ^ { ( \mathbf a ) } ( \omega ) + 2 \gamma _ { r } ^ { ( \mathbf a ) } ] - Q _ { 2 } ^ { ( \mathbf a ) } } \ ( 1 8 ) \\ \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \end{cases}$$

and substitute them into Eq. (13). Thus, the excitation of the antisymmetric mode results in changing the reflection coefficient compared to the case without a conductivity grating, r A ≈ r (1 + f 1 A / 2 E 0 ).

## III. NONLINEAR PLASMON-TO-PLASMON CONVERSION

First we notice that, for the case of TM-polarized waves with the tangential electric field of the form E τω = z 0 Ee iqz (the monochromatic time dependence ∼ exp( -iωt ) is omitted), the induced second-harmonic current in graphene depends on the in-plane momentum q manifesting its nonlocal nature. It is expressed as 26,31

$$\mathfrak { j } _ { 2 \omega } = - \mathbf z _ { 0 } \frac { 3 } { 8 } \frac { e ^ { 3 } v _ { F } ^ { 2 } } { \pi \hslash ^ { 2 } \omega ^ { 3 } } q E ^ { 2 } e ^ { i 2 q z } \,, \quad \quad ( 1 9 )$$

where v F ≈ c/ 300 is the Fermi velocity. This result is derived in the semiclassical limit from the Boltzmann kinetic equation written for Dirac electrons with a linear energy dispersion under the approximations /planckover2pi1 ω ≤ E F ,

k B T /lessmuch E F and qv F /ω /lessmuch 1, where E F is the Fermi energy, k B is the Boltzmann constant, T is the temperature. Thereby, a nonlinear source associated with the antisymmetric mode will generate a synchronous symmetric mode. Importantly, in general, when discussing the parametric interaction that involves the modes of different symmetries, the phase matching is not a sufficient condition, and the overlaps of the modes with the respective nonlinear sources should be examined. Since a quadratic nonlinear source localized at the graphene layers is symmetric, the conversion between the FF antisymmetric mode and SH symmetric mode becomes possible.

Similar to the analysis of the fundamental frequency field, we employ here the perturbational method and derive the equations for the envelopes of the secondharmonic fields B 1 2 , ,

$$& \text{harmonic fields $\mathcal{B}_{1,2}$,} \\ & \pm 2 i k _ { \text{sp} } ^ { ( s ) } ( 2 \omega ) \left ( \frac { \partial \mathcal { B } _ { 1, 2 } } { \partial z } \pm ( \gamma ^ { ( s ) } + \gamma _ { r } ^ { ( s ) } ) \mathcal { B } _ { 1, 2 } \right ) } \quad \begin{array} { c } \alpha \ u e \\ \varepsilon \ = \ 1 \\ \text{metric and } r \\ \text{and } h \\ ( 2 0 ) \quad \text{dots in} \\ \text{harmoc} \end{array}$$

where ˜ ∆ k = k (s) sp (2 ω ) -2 k (a) sp ( ω ) is a small detuning, and other parameters are

$$\text{where} \, \Delta k = k _ { s p } ^ { ( s ) } ( 2 \omega ) - 2 k _ { s p } ^ { ( a ) } ( \omega ) \, \text{is} \, a \, \text{mal} \, \text{detuning}, \, \text{and} \\ \text{other parameters are} \\ q ^ { ( s ) } = 4 k _ { 0 } ^ { 2 } \frac { \varepsilon ^ { 2 } } { \kappa ^ { ( s ) 3 } } \left [ 1 + \frac { \sinh \left ( \kappa ^ { ( s ) } d \right ) - \kappa ^ { ( s ) } d } { 2 \cosh ^ { 2 } \left ( \kappa ^ { ( s ) } d / 2 \right ) } \right ] \,, \\ \gamma ^ { ( s ) } = \frac { 8 \pi } { c } \sigma ^ { ( R ) } ( 2 \omega ) \frac { k _ { 0 } \varepsilon } { k _ { s p } ^ { ( s ) } ( 2 \omega ) q ^ { ( s ) } } \,, \\ \gamma ^ { ( s ) } _ { r } = \left ( \frac { \pi \sigma ^ { ( I ) } ( 2 \omega ) f _ { 2 } } { c } \right ) ^ { 2 } \frac { 2 k _ { 0 } \varepsilon ^ { 1 / 2 } } { k _ { s p } ^ { ( s ) } ( 2 \omega ) q ^ { ( s ) } } \,, \\ \Phi ^ { ( s ) } _ { \Phi } \Phi ^ { ( s ) } _ { \Phi } \\ g = - i \frac { 1 6 \pi } { c } \sigma _ { 2 } ( \omega ) \frac { k _ { 0 } \varepsilon } { q ^ { ( s ) } } \,, \quad ( 2 1 ) \, \right ) \, \text{is} \, \quad \\ \text{with} \, \sigma _ { 2 } ( \omega ) \, \text{ denoting} \, \text{ the quadratic} \, \text{conductivity} \, \text{ of} \\ \text{random} \, \text{do} \, \text{mod} \, \text{eq}$$

with σ 2 ( ω ) denoting the quadratic conductivity of graphene defined as

$$\sigma _ { 2 } ( \omega ) = - \frac { 3 } { 8 } \frac { e ^ { 3 } v _ { F } ^ { 2 } } { \pi \hslash ^ { 2 } \omega ^ { 3 } } k _ { \text{sp} } ^ { ( \text{a} ) } ( \omega ). \\.. \. \.$$

Within the undepleted pump approximation, when the amplitude at fundamental frequency is not affected by nonlinearity, for the exact phase matching, ˜ ∆ k = ∆ k = 0, and homogeneous case discussed above, we have

$$\mathcal { B } _ { 1 } = \mathcal { B } _ { 2 } \equiv \mathcal { B } = \frac { g \mathcal { A } ^ { 2 } } { 2 i k _ { s p } ^ { ( s ) } ( 2 \omega ) ( \gamma ^ { ( s ) } + 2 \gamma _ { r } ^ { ( s ) } ) }. \quad ( 2 2 )$$

Since the radiation associated with SH field is emitted equally to both directions, the normalized conversion efficiency can be introduced as a ratio of the doubled SH energy flux density in the upper half-space to the energy flux of the incident fundamental wave. Equivalently, this definition takes the following form

$$\eta = \frac { 2 | \bar { E } _ { 0 z } ( 2 \omega, x = d / 2 ) | ^ { 2 } } { | E _ { 0 } | ^ { 2 } }, \quad \begin{array} { c c c } & \text{inc} \\ ( 2 3 ) & \text{are} \\ \text{no} \\ \text{We} \end{array}$$

$$\text{where } \bar { E } _ { 0 z } ( 2 \omega, x = d / 2 ) = - \frac { 2 \pi } { c \sqrt { \varepsilon } } i \sigma ^ { ( I ) } ( 2 \omega ) f _ { 2 } \mathcal { B }.$$

Figure 2. (Color online) Plasmonic guide index vs. frequency at the fixed separation between the graphene layers, d = 123 nm, calculated for the doping level of E F = 0 3 eV, . N = 5, ε = 1. Blue and red solid curves correspond to the symmetric and antisymmetric modes, respectively. Dashed gray and red lines show dispersion of an isolated graphene sheet and half-frequency antisymmetric mode, for reference. Color dots in the horizontal axis mark the fundamental and secondharmonic frequencies.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0555v1-figure-bc3b0c8a-4635-4976-8d71-ba3d86bbff4c.jpg)

Figure 3. (Color online) Operational fundamental frequency as a function of the separation between the layers for the exact phase-matching at E F = 0 3 eV, . N = 5, ε = 1. The point corresponds to the parameters used in Fig. 2.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0555v1-figure-a61e2d7c-0250-4716-b1e2-d18ff3f84b33.jpg)

## IV. DISCUSSIONS

In the framework of the nonlinear amplitude equations, we can now analyze the parametric effects and estimate the frequency conversion efficiency. To make a period of the conductivity modulation realistic, in our calculations we employ multilayer graphene that effectively increases the equivalent surface conductivity leading to a reduction of the wavenumber of the p-polarized plasmons supported by multilayer graphene structures 11,29 . We take σ ω ( ) = Nσ s ( ω ) for a randomly stacked multilayer graphene film consisting of N layers, with each

layer characterized by a surface conductivity 32

$$\text{layer cancerize} \, \text{by a surface concuativity} \stackrel { \text{so } \text{na} } { \text{$\sigma$} } ( \omega ) = \frac { i e ^ { 2 } } { \pi \hslash } \left \{ \frac { \mathcal { E } _ { F } } { \hslash \left ( \omega + i \tau _ { \text{intra} } ^ { - 1 } \right ) } + \frac { 1 } { 4 } \ln \left | \frac { 2 \mathcal { E } _ { F } - \hslash \omega } { 2 \mathcal { E } _ { F } + \hslash \omega } \right | } \right \}, \, ( 2 4 ) \\ \text{where for doped graphene we assume $\hslash \omega < 1.67\mathcal{E}_{F}$ and} \quad \text{$\text{by}$} \\ k _ { B } T \ll \mathcal { E } _ { F }, \, \text{and } \tau _ { \text{intra} } ^ { - 1 } \, \text{is the relaxation rate}.$$

To reveal a possibility of plasmonic frequency conversion in the graphene waveguide, in Fig. 2 we plot the frequency-dependence of the plasmonic guide indices for the symmetric and antisymmetric modes. We observe that the first harmonic antisymmetric and secondharmonic symmetric modes of the waveguide can be matched at /planckover2pi1 ω = 0 12 eV ( . λ = 2 π/k 0 ≈ 10 6 . µ m, corresponding to one of the radiation lines of high power CO 2 lasers), k (a) sp ≈ 10 5 . k 0 .

where for doped graphene we assume /planckover2pi1 ω &lt; 1 67 . E F and k B T /lessmuch E F , and τ -1 intra is the relaxation rate.

Figure 3 shows the dependence of the fundamental frequency on the separation between the layers. Note while discussing the dispersion relations and the possibility of phase-matching conditions we take σ = iσ ( I ) and do not take losses into account, however these losses are included into our considerations perturbatively. Calculated for f 1 = f 2 = 0 1, . a ≈ 1 µ m, τ intra = 0 3 ps, the incident . wave intensity of 1 MW/cm 2 and the other parameters of Fig. 2, efficiency is η ∼ 8 7 . × 10 -7 . Overall, for this structure the smallness parameter is found to be µ ∼ 0 1, .

- 1 A. N. Grigorenko, M. Polini, and K. S. Novoselov, Nature Photon. 6 , 749 (2012).
- 2 Q. Q. Bao and K. P. Loh, ACS Nano 6 , 3677 (2012).
- 3 M. Jablan, M. Soljacic, and H. Buljan, Proc. IEEE 101 , 1689 (2013).
- 4 X. Luo, T. Qiu, W. Lu, and Z. Ni, Mater. Sci. Eng., R 74 , 351 (2013).
- 5 T. Stauber, J. Phys.: Cond. Mat. 26 , 123201 (2014).
- 6 T. Low and P. Avouris, ACS Nano 8 , 1086 (2014).
- 7 F. J. Garcia de Abajo, ACS Photonics 1 , 135 (2014).
- 8 F. H. L. Koppens, D. E. Chang, and F. J. Garcia de Abajo, Nano Lett. 11 , 3370 (2011).
- 9 M. Gullans, D. E. Chang, F. H. L. Koppens, F. J. Garcia de Abajo, and M. D. Lukin, Phys. Rev. Lett. 111 , 247401 (2013).
- 10 M. L. Nesterov, J. Bravo-Abad, A. Yu. Nikitin, F. J. Garcia-Vidal, and L. Martin-Moreno, Laser Phot. Rev. 7 , L7 (2013).
- 11 D.A. Smirnova, I.V. Shadrivov, A. I. Smirnov, and Y.S. Kivshar, Laser Phot. Rev. 8 , 291 (2014).
- 12 X. Yao, M. Tokman, and A. Belyanin, Phys. Rev. Lett. 112 , 055501 (2014).
- 13 G. W. Hanson, J. Appl. Phys. 104 , 084314 (2008).
- 14 P. I. Buslaev, I. V. Iorsh, I. V. Shadrivov, P. A. Belov, and Yu.S. Kivshar, JETP Lett. 97 , 535 (2013).
- 15 D.A. Smirnova, A. V. Gorbach, I. V. Iorsh, I. V. Shadrivov, and Y.S. Kivshar, Phys. Rev. B 88 , 045443 (2013).
- 16 F. Wang, Y. Zhang, C. Tian, C. Girit, A. Zettl, M. Crommie, and Y. R. Shen, Science 320 , 206 (2008).

so that our asymptotic approach is well justified.

## V. CONCLUSIONS

By employing the approximation of slowly-varying field amplitudes, we have derived nonlinear equations for describing the second-harmonic generation from a double-layer graphene structure with modulated conductivity. We have revealed that the quadratic phase matching between the plasmonic modes of different symmetries becomes possible in a planar waveguide geometry with conductivity modulation playing a role of an effective grating that couples the external pumping radiation to the waveguiding modes with the subsequent parametric conversion into radiation at the double frequency.

## ACKNOWLEDGMENTS

The authors are grateful to I.V. Shadrivov and A.I. Smirnov for in-depth comments, and also thank S.A. Mikhailov for useful discussions at an early stage of this project. This work was supported by the Australian National University.

- 17 A. Vakil and N. Engheta, Science 332 , 1291 (2011).
- 18 H. Liu, Y. Liu, and D. Zhu, J. Mater. Chem. 21 , 33353345 (2011).
- 19 Y. V. Bludov, N. M. R. Peres, and M. I. Vasilevskiy, Phys. Rev. B 85 , 245409 (2012).
- 20 N. M. R. Peres, A. Ferreira, Y. V. Bludov, and M. I. Vasilevskiy, J. Phys.: Cond. Mat. 24 , 245303 (2012).
- 21 T. M. Slipchenko, M. L. Nesterov, L. Martin-Moreno, and A. Yu. Nikitin, J. Opt. 15 , 114008 (2013).
- 22 Yu. V. Bludov, A. Ferreira, N. M. R. Peres, and M. I. Vasilevskiy, Int. J. Mod. Phys. B 27 , 1341001 (2013).
- 23 R. W. Boyd, Nonlinear Optics (Elsevier Science, USA, 2003).
- 24 Yu. S. Kivshar and G. A. Agrawal, Optical Solitons: From Fibers to Photonic Crystals (Academic, San Diego, 2003).
- 25 M.M. Glazov, JETP Lett. 93 , 366 (2011).
- 26 S. A. Mikhailov, Phys. Rev. B 84 , 045432 (2011).
- 27 M.M. Glazov and S.D. Ganichev, Phys. Rep. 535 , 101 (2014).
- 28 A.V. Gorbach, Phys. Rev. A 87 , 013830 (2013).
- 29 D.A. Smirnova, I. V. Iorsh, I. V. Shadrivov, and Y.S. Kivshar, JETP Lett. 99 , 456 (2014).
- 30 G.A. Korn and T. M. Korn, Mathematical Handbook for Scientists and Engineers (McGraw-Hill, New York, 1961).
- 31 D.A. Smirnova, I. V. Shadrivov, A.E. Miroshnichenko, A.I. Smirnov, and Y. S. Kivshar, Phys. Rev. B 90 , 035412 (2014).
- 32 L. A. Falkovsky and A. A. Varlamov, The Eur. Ph. J. B, 56 , 281 (2007).
Daria Smirnova Yuri S. Kivshar
Second-harmonic generation in subwavelength graphene waveguides
We suggest a novel approach for generating second-harmonic radiation in
subwavelength graphene waveguides. We demonstrate that quadratic phase matching
between the plasmonic guided modes of different symmetries can be achieved in a
planar double-layer geometry when conductivity of one of the layers becomes
spatially modulated. We predict theoretically that, owing to graphene nonlocal
conductivity, the second-order nonlinear processes can be actualized for
interacting plasmonic modes with an effective grating coupler to allow external
pumping of the structure and output of the radiation at the double frequency.