CC BY
28
NANOSYSTEMS: PHYSICS, CHEMISTRY, MATHEMATICS, 2014, 5 (1), P. 155-159
THE EVOLUTION OF FEW CYCLES OPTICAL PULSES IN A DOUBLE-LAYER GRAPHENE - BORON NITRIDE TAKING INTO ACCOUNT NONLINEARITY
OF A MEDIUM
A. V. Pak1, M.B. Belonenko2
1Volgograd State University, Volgograd, Russia 2Volgograd Institute of Business, Volgograd, Russia
pak.anastasia@gmail.com, mbelonenko@yandex.ru PACS 75.75 -, 73.20 Hb, 73.22 Pr
The propagation of ultra-short optical pulses in a thin film created by graphene grown on a boron nitride was considered, taking into account non-linear medium characteristics. Electron conduction in such a system described with a long-wave effective Hamiltonian for the low temperatures media. The electromagnetic field is taken in the framework of classical Maxwell equations. Dependence of the pulse shape on the initial pulse amplitude and the parameters of the linear and nonlinear polarization is shown.
Keywords: Graphene, Boron Nitride, Optical pulse, Nonlinearity.
1. Introduction
Recently, the number of studies on boron nitride-supported graphene has increased [1-8]. Primarily, this is due to graphene's unusual properties [9]. Note that boron nitride has a similar hexagonal lattice to graphene, as well as the fact that the ionic character of interatomic bonds in the hexagonal boron nitride (h-BN) leads to the absence of surface-based 'dangling' covalent bonds and charge trapping [8]. As found in [8], the roughness of graphene on h-BN is much smaller than that of graphene on SiO2, and charge fluctuations are two orders of magnitude weaker. In general, the electronic characteristics of graphene on h-BN are almost the same as that of free graphene. Thus, to explore graphene on a substrate is much easier and more convenient [8].
The evolution of an ultra-short optical pulse propagating in a double-layer structure of graphene - boron nitride in non-magnetic environments would be revealed taking into account the nonlinear polarization of the medium.
2. Basic equations
We considered a layer of graphene on a substrate of boron nitride. The Hamiltonian we have chosen, in a long-wave approximation, can be written in matrix form as:
H (k)
<0 k* 0 t )
k 0 0 0
0 0 A f *
0 f -A
(1)
where t is the electron overlap integral between the layers of graphene and boron nitride; A is the band gap for boron nitride; k = vFG (kx + iky), vFG - is the Fermi velocity for graphene;
kx, ky - are the electron pulse components: f = vfbn (kx + iky); vfbn - is the Fermi velocity for boron nitride.
The Hamiltonian (1) can be rewritten using a block matrix structure [10]:
H (k)
fo k k 0 0 0 A
* 0 0
t \
0
f *
H11 H12 H21 H22
\t 0 f -a;
In the case of a large band gap in boron nitride compared to the electron's energy, a long-wave approximation can be considered. This makes it possible to write the effective Hamiltonian analogously to bigraphene [10]:
H"> . H11 - H12H-21H21 = -1 _Afk 7fk*kA
(2)
The Hamiltonian (2) is easily diagonalized, which gives the electronic spectrum:
£ (kx, ky) — 2 A
1 -
Jfg
(kx + k,
V
t2
+
\
1 + f (k2 + ky) \2 + (kX + k2)2
t2
A2t2
/
(3)
where vFG, - the Fermi velocity of electrons for graphene and boron nitride, respectively. According to quantum-mechanical laws, the presence of an external electric field E, di-
1 SA
rected along the X axis, the Coulomb gauge can be chosen in the following form: E =---—.
c at
e
It is necessary to replace the momentum with a generalized momentum: p ^ p--A (e is the
c
electron charge). In this case, the effective Hamiltonian (2) can be rewritten as:
H
(p _ eA(t)J apa,
(4)
pa
where 0+., apa -are the creation and annihilation operators of electrons with quasi-momentum p and spin a. The vector-potential A is considered as A = (0,0, A(x, t)).
Maxwell's equations for polarization of the medium can be written as [11]:
<92E ôx2
1 d 2E 4ndj
+ Tôt
ôt2
4n d 2P ^^ ôt2 '
(5)
Here, we neglect the laser beam diffraction in directions perpendicular to the beam propagation axes. In Eqn. (5), E is an electric field of light wave, P = aE + ^ |E| E - is a polarization of the medium, t is a time, c is a speed of light in a vacuum. There is a simple model for the medium nonlinearity when the polarization vector is considered to be parallel to E.
We write the expression for the current density:
Jo
ps
vs(p _ ^A(t)) (a+saps)
(6)
2
ds (p)
where vs(p) = —s—, and brackets denote averaging with the nonequilibrium density matrix
dp
p(t): (B> = Sp(B(0)p(t)).
Further, we consider the case of low temperatures, when the sum (6) contributes only a small area in momentum space near the Fermi level. Therefore, we rewrite the expression for the current density in the form:
A A
e
Jxdpyvy ^r ~ c
AA
j = e / dpxdpyVy (p — ^A(x,t)^ . (7)
The range of pulses integrated in (7) was determined from the particles number equality:
A A A A
dPxdPy = / dPxdPy(^a+x>PyaPx,Py^
A A A A
3. The results of numerical simulation
Equation (5) was solved numerically using a direct finite-difference cross-like scheme [12]. Steps by time and coordinates were determined using standard conditions of stability, strides of finite-difference scheme were halved serially, until the solution did not change in the eighth sign. The initial conditions were chosen as an ultra short optical pulse consisting of a single oscillation of the zero-width field, respectively, that can be specified by setting the potential A as:
A(x,t) = Q exp(-(x — vt)2/y),
f, 2\ 1/2 (8)
Y = (1 — v ) ' ,
where Q is the amplitude, v is the initial ultra-short pulse velocity on the sample input. This initial condition corresponds to the fact that the sample is fed an ultra-short pulse, consisting of a single oscillation of the electric field. The values of energy parameters are expressed in A units. The resulting evolution of the electromagnetic field propagating along the sample is shown in Fig. 1.
The amplitude of the pulse decreased when the linear polarization of the medium gradually increased. This fact can be interpreted as the pulse spending part of its energy on medium polarization.
The pulse behavior presented in Fig. 1 concerns the presence of dispersion, which leads to a broadening of the optical pulse, as well as non-linearity in the same equation, thus leading to a 'narrowing' of the pulse. The competition between these two terms leads to a deformation of the initial pulse shape and rise to its stable form.
More clearly, the effects associated with nonlinear pulse appear at the front and lead to the formation of additional peaks and the broadening of the pulse which can be explained by an imbalance between dispersion and nonlinearity in the system. This is clearly seen in Fig. 3. The dependence of the electric field pulse on the magnitude of the medium's nonlinear polarization is shown in Fig. 3.
As can be seen in Fig. 3, involving non-linear polarization leads to the appearance of a wave front for the second maximum. As in the case of linear polarization, the initial pulse amplitude decreases, which is to be expected.
ISO
t\
.V y \
I''. \\\ \\\i \ V №
I 1 r 1 1 1 1
50 1» 150 200 250 300 350
t, a.u.
Fig. 1. Time dependence of the electric field: a) solid line - a = 0.0; b) dotted line - a = 0.3; c) dashed line - a = 0.5, the nonlinear polarization is absent. Time, a.u. is time 3 ■ 10-16 s.; electric field a.u. is 107 V/m
Fig. 2. Time dependence of the pulse form for different cases of pulse amplitude value: a) solid line - Q = 2; b) dotted line - Q = 3; c) dashed line - Q = 5. Time and electric field a.u. as in Fig. 1
4. Conclusion
As follows from the obtained results, the stable ultra-short optical pulses can undergo propagation in hexagonal boron nitride-supported graphene in a medium with nonlinear polar-
ization.
By increasing the initial pulse amplitude, the wave front is broadened, and there is a second pulse of lower intensity. This effect may be useful in the development of hybrid devices based on the effect of light interaction with graphene electrons.
150
100
JO
3
CO
LU
-50
-100
-150
A* / A1 ■r\ y \
f i i i i i y \ r 1 1 1 1
50
100
150
:oo
:so
Î00
3Í0
t, a.u.
Fig. 3. Time dependence of the electric field on time taking into account the nonlinear polarization of the medium: a) solid line - ft = 0.03; b) dotted line -ft = 0.4; c) dashed line - ft = 0.6 the dispersion is almost zero, a = 0.01. Time and electric field a.u. as in Fig.1
Acknowledgments
This work was supported by Russian Foundation for Basic Research (grant No. 12-0231654).
References
[1] Giovannetti G., Khomyakov P.A., et al. Substrate-induced band gap in graphene on hexagonal boron nitride: Ab initio density functional calculations. Phys. Rev. B, 76, P. 073103 (2007).
[2] Greber T. Graphene and boron nitride single layers. In: The Handbook of Nanophysics, ed. by Sattler K., Taylor and Francis Books, Inc, Oxfordshire, 2010.
[3] Avetisyan A.A., Partoens B., Peeters F.M. Stacking order dependent electric field tuning of the band gap in graphene multilayers. Phys. Rev. B, 79, P. 035421 (2009).
[4] Wakabayashi K., Takane Y., Sigrist M. Perfectly conducting channel and universality crossover in disordered graphene nanoribbons. Phys. Rev. Lett., 99, P. 036601 (2007).
[5] Blase X., Rubio A., Louie S.G., Cohen M.L. Stability and band-gap constancy of boron-nitride nanotubes. Europhys. Lett., 28 (5), P. 355-340 (1994).
[6] Chen Y., Zou J., Campbell S.J., Le Caer G. Boron nitride nanotubes: Pronounced resistance to oxidation. Appl. Phys. Lett., 84, P. 2430-2432 (2004).
[7] Breshenan M.S., Hollander M.J., et al. Integration of hexagonal boron nitride with quasi-freestanding epitaxial graphene: toward wafer-scale, high-perfomance devices. ACS Nano, 6 (6), P. 5234-5241 (2012).
[8] Xue J., Sanchez-Yamagishi J., et al. Scanning tunneling microscopy and spectroscopy of ultra-flat graphene on hexagonal boron nitride. Nature Mater., 10, P. 282-285 (2011).
[9] Rutter G.M., Jung S., et al. Microscopic polarization in bilayer graphene. Nature, 7, P. 649-655 (2011).
[10] Cortijo A., Guinea F., Vozmediano M.A.H., Geometrical and topological aspects of graphene and related materials. arXiv: 1112.2054v1 (2011).
[11] Landau L.D., Lifshitz E.M. Theoretical physics. II. Field theory, Science, Moscow, 1988, 512 p.
[12] Bakhvalov N.S. Calculus of approximations (analysis, algebra, ordinary differential equations), Science, Moscow, 1975, 632 p.
