Propagation of femtosecond pulses in carbon nanotubes Текст научной статьи по специальности «Физика»

PROPAGATION OF FEMTOSECOND PULSES IN CARBON NANOTUBES

N.N. Konobeeva1, M.B. Belonenko2

1Volgograd State University, Volgograd, Russia 2Volgograd Institute of Business, Volgograd, Russia

yana_nn@inbox.ru, mbelonenko@yandex.ru PACS 73.63.Fg, 78.67.Ch

The generation of higher harmonics of carbon nanotubes interacting with femtosecond laser pulses was investigated. The analysis was conducted on the basis of quantum kinetic equation for the n-electrons involved in the inside of the band and inter-band transitions. The dynamics of the electromagnetic pulse, depending on the parameters of the problem, were studied.

Keywords: Carbon nanotube, femtosecond pulse, electromagnetic field. 1. Introduction

This paper investigates the dynamics of femtosecond pulses propagating in carbon nanotubes (CNTs). Significant progress in the synthesis of new mesoscopic structures and their great promise for use in nanoelectronics and optics has led to increased research in this area. The unusual manifestation of nonlinear optical properties of these materials when interacting with strong electromagnetic fields was expected. Theoretical and experimental studies predict that the strong nonlinear response is inherent in metal clusters [1-3], fullerenes [4-6], carbon nanotubes and composites based on them [7, 8]. In particular, the generation of higher harmonics, previously observed in noble gases [9,10], has been predicted for mesoscopic structures. Moreover, a distinctive feature of the CNT is a strong selectivity of their high harmonic spectrum. A theoretical study of the nonlinear response of isolated nanotubes by intense laser radiation was studied in this work. A full quantum-theoretical approach to the one-electron approximation in the tight-binding model was used.

2. Statement of a problem 2.1. Kinetic equations

We consider an infinitely long single-walled carbon nanotube oriented along the z axis and irradiated normal to the axis of the ultra-short pulse with a natural frequency w0, polarized to this axis: E(r) = ezEz (x,y) (ez is the unit vector along the z axis). CNTs are considered as a single-layer graphene sheet which is scrolled into a cylinder. We take into account only n-electrons, suggesting that their motion can be described in the tightly-bound approximation [11, 12]. CNT radius is much smaller than the wavelength of the field, which allows us to neglect the spatial inhomogenity of the field in the tubes.

Following [13] in the one-electron approximation, the Schrdinger equation becomes:

db h2

^ = - 2m Ab + [W(r) - e(E ■ r)] b. (1)

Neglecting the harmonic solution of (1) can be written as:

^ = J] C,(pM(p, r), (2)

i, p

where l is the set of quantum numbers characterizing the state of n-electrons with a given quasi-momentum. (p, r) = h-0'5 exp(ipr/h) ui,p (r) are the Bloch functions with an amplitude ui,p (r), which is periodical to an arbitrary vector lattice a: ui,p (r + a) = ui,p (r). Here a = n1a1 + n2a2, n1, n2 are the integers, a1? a2 are the elementary vectors of the graphene hexagonal lattice.

The dispersion law for zig-zag CNTs can be written in the following form [14]:

£s(p) = ±Y0\/1 + 4cos(ap) cos(ns/m) + 4cos2(ns/m), (3)

where s = 1, 2,... m, (m, 0) is type of CNT, y0 is the hopping integral (2.7 eV), a = 3b/2h, b = 0.142 nm is the distance between the adjacent carbon atoms.

From equations (1) and (2) with the spectrum of nanotubes:

Ho =l H,°<p) Hi2o(p) >•

which can be rewritten in the tightly-bound approximation as:

H12 (p) = — Yo it ex^, (5)

j=i

where Tj is the vector connecting the atom to its nearest neighbors. For density matrix elements taking into account F = !R(pcv); $ = pcv); p = pcc, pcv is the density matrix, u is the transition frequency, is the electron charge, we can obtain the following expression:

% + eEz dP- = 2

dt dpz h

dF dF — + e£Z-— = u$, dt dpz

d$ 5$ e (6)

af + eEz ^ = <2P — D — Rii'.p = ih f U. ^«,„ I <Pr.

"'P - 2 J V"''p dpz dpz "'-p

n

The integration domain is the volume Q of two-dimensional unit cell. The initial conditions are following:

Pt=o - Fo (ec (pz, s)) ; F|i=0 - $|i=0 - 0. (7)

This means that at room temperature, the electrons are distributed according to the equilibrium Fermi distribution with zero chemical potential (p — 0).

The boundary conditions reflect the periodicity of solutions in the space of quasi-momentum (similar for F and $):

P |t, - P ft, -. (8)

Propagation of femtosecond pulses in carbon nanotubes 3. Electric current in nanotube

A quantum-mechanical operator of current density can be written in the following form [13]:

^ ieh ( d d \

(r) = - 2mol dz6 (r - r° +6 (r - r') ôï) •

(9)

Total current density was decomposed into two components: jz = jZ ) + jV, where intraband transitions is responsible for:

4e f dec (p)

(2nh)2 J dpz 1ZB

p(i p)

and inter-band:

jZ2)

8e

(2nh)2 h

ec (p) RcV (pz ,s)$(i, p) d2p.

(10)

(11)

1ZB

We take into account that pvv + pcc =1 and ev = — ec.

1 d a

Choosing the gauge field, E =---—, (where is the light velocity in a vacuum). We

c dt

need to change the momentum to a generalized momentum: p ^ p — eA/c, so:

j ( p - —< F

$ I , $

F

pp

Then, equation (5) takes the following form with the replacement p — 0.5 ^ p:

2

pt = — AtR$, c

< Ft = w$, 2

$t = — AtRp — wF.

c

The matrix of the system coefficients can be written as:

(12)

0 0 A

j = 0 0 w

A —w 0

Where the notation: A = —2c 1AtR and according to the approximation described in [15-19]:

ex = /+j + 4e

(A2 — w2) (A2 — w2)2

3!

+

+ ... l+j2 I XT +

5!

1 (A2 — w2) (A2 — w2)

22

2! 4!

2 / \2\k

+

6!

+

jz1)

(2nh)2

dec (p — eA/c) A ^ ^ _(A (i) — w (p) )

1

1ZB

dpz

k=0

(2k + 2)! (i + exp(—ec (p) /kBT))

d2p,

j-z2) =

8e

(2nh)2 h

1ZB

(1 + exp(—ec (p) /kBT))

2\k

—2A (i) ec (p — eA/c) Rc, (p) ^ (A (i)2 — w (p)2)

£

k=0

(2k + 1)!

-d2p,

where kB is the Boltzmann constant, T is the temperature.

The Maxwell equations for the dielectric non-magnetic medium can be written as following [16]:

d2A 2 No d4A 2—, 4n , , 1, A ,02A 12nn d2A f dA I2

+ —a —--0bA + — (ji + j2) — ^ (1 + 4na) —--— ( — ) =0, (13)

dx2

ôt4

ôt2

c4 ôt2 V dt

where A is the vector-potential, PL = aE is the polarization, PNL = n |E| E is the nonlinear part of polarization, t is the time. We consider a simple model, where the polarization vector is parallel to the vector E is the electric field of light wave, N0, a, a1, ..., b, b1, ... are the empirical constants of medium dispersion [17].

4. Numerical analysis and results

Equation (13) was solved numerically [20]. Boundary conditions took on a Gaussian form with one field- (14a) and two field oscillations (14b):

(A)

(B)

40

20

-20

-40

1 1 ftb ll 1 1 a i i i 1 j "

¡1 11 Ja__J\

- j1' V 1 1 if i i i №' -1 1 1

50

100

150

200

250

300

-13

t,16 s

350

400

Fig. 1. Time dependence of the electric field at various space points a) x - 0.1 ■ 10-5 m; b) x - 0.2 ■ 10-5 m. (A) two field oscillations; (B) one field oscillation

A(0,t) = Q • e-(ut)2/Y, dA(0,t) = 2Qute-(wt) 2/y dx

A(0,t) = Q • e

Y

-(ut)2/Y

sin(ut),

dA(0,t)

dx

2Qute-(ut)2/y sin(-ut) + Q • e-(ut)2/7 cos(ut),

(14a)

(14b)

Y

where Q is the pulse amplitude, u is the pulse velocity, y2 = 1/(1 — u2/A2).

The evolution of the electric field during its propagation in the sample is shown in

Fig. 1.

It can be seen from this dependence that when the propagation distance increases, pulses drift away from each other. Moreover, in the case of two field oscillations, the amplitude of the main peak is stable, and in the case of one oscillation there is a decrease in the main pulse amplitude and an increase in the trailing pulse. This can be attributed to the dispersion characteristics of the medium.

The dependence of the electric field on the initial pulse amplitude is shown in Fig. 2.

As expected, low-amplitude pulses propagate almost without distortion and exhibit only the dispersion spread. However, pulses with high amplitudes exhibit a larger distortion caused by both the interference of the edge pulse and the specific form of the medium's nonlinearity.

Figure 3 demonstrates that the dispersion constants have a great influence on the electric field (and the effect of the b constant is manifested much more strongly compared to the constant a). In our opinion, this effect is due to this type of dispersion reducing the linear response of the system to an external field.

Thus, based on the results of numerical calculations, one can draw the conclusion that the pulse propagation is determined primarily by the dispersion characteristics of the medium, as well as the process of establishing a balance between the dispersion and interference effects of the wave front.

Fig. 2. Time dependence of the electric field at various initial amplitudes of the pulse (x = 0.2 ■ 10-5 m): a) Q = 1 a.u.; b) Q = 2 a.u. (the case of two field oscillations)

40

^ 20

ь

it

-20

f

-40

0

100

200 -13

300

400

U9 s

Fig. 3. Time dependence of the electric field at various dispersion constants (x = 0.1 ■ 10-5 m): a) b = 0.001 a.u; b) b = 0.005 a.u. (the case of two field oscillations)

Acknowledgments

This study was supported by the Russian Foundation for Basic Research (project No. 12-

[1] Calvayrac F., Reinhard P.-G., Suraud E., Ullrich C.A. Nonlinear electric dynamics in metal clusters. Phys. Rep., 337, P. 493-578 (2000).

[2] Veniard V., Taieb R., Maguet A. Atomic clusters submitted to an intense short laser pulse: A density-functional approach. Phys. Rev. A, 65, P. 013202(1-7) (2001).

[3] Daligault J., Guet C. Kinetic and Coulombic effects in the explosion dynamics of metal clusters in intense femtosecond laser fields. Phys. Rev. A, 64, P. 043203(1-5) (2001).

[4] Hertel I.V., Steger H., et al. Giant plasmon excitation in free C6o and C70 molecules studied by photoioniza-tion. Phys. Rev. Lett., 68, P. 784-787 (1992).

[5] Hunsche S., Starczewski T., et al. Ionization and fragmentation of C60 via multiphoton-multiplasmon excitation. Phys. Rev. Lett., 77, P. 1966-1969 (1996).

[6] Cambell E.E.B., Hansen K., et al. From above threshold ionization to statistical electron emission: The laser pulse-duration dependence of C60 photoelectron spectra. Phys. Rev. Lett., 84, P. 2128-2131 (2000).

[7] Stanciu C., Ehlich R., et al. Experimental and theoretical study of third-order harmonic generation in carbon nanotubes. Appl. Phys. Lett., 81, P. 4064-4066 (2002).

[8] Chen Y.-C., Raravikar N.R., et al. Ultrafast iptical switching properties of single-wall carbon nanotube polymer composites at 1.55 mm. Appl. Phys. Lett., 81, P. 975-977 (2002).

[9] Ferray M., L'Huillier A., et al. Multiple-harmonic conversion of 1064 nm radiation in rare gases. J. Phys. B: At. Mol. Opt. Phys., 21, L31 (1988).

[10] Chang Z., Rundquist A., et al. Generation of coherent soft X rays at 2.7 nm using high harmonics. Phys. Rev. Lett., 79, P. 2967-2970 (1997).

[11] Il'inskii Yu.A., Keldysh L.V. Electromagnetic Response of Material Media. Plenum, New York, 1994, 304 p.

[12] Davydov A.S. Quantum Mechanics. Pergamon, Oxford, 1976, 168 p.

[13] Slepyan G.Ya., Khrutchinski A.A., Nimelentsau A.M., Maksimenko S.A. High-order optical harmonic generation on carbon nanotubes: quantum-mechanical approach. Int. J. Nanoscience, 3, P. 343-354 (2004).

[14] Wallace P.R. The Band Theory of Graphite. Phys. Rev., 71, P. 622-634 (1947).

02-31654).

References

[15] Belonenko M.B., Glazov S.Yu., Meshcheryakova N.E. Dynamics of electromagnetic pulses with wide spectra in semiconductor superlattices. J. Russian Laser Research, 29 (2), P. 114-122 (2008).

[16] Vongradova M.B., Rudenko O.V., Sukhorukov A.P. Wave theory. Nauka, Moscow, 1990, 432 p.

[17] Yanyushkina N.N., Belonenko M.B. The influence of self-nonlinearity on the propagation of ultrashort optical pulses in carbon nanotubes in dispersive non-magnetic dielecrtric media. Zh. Techn. Fiz., 83 (4), P. 155-158 (2013).

[18] Sazonov S.V. The effect of resonance transparency anisotropic medium with a permanent dipole moment. Zh. Eksp. Teor. Fiz., 124, P. 803-819 (2003).

[19] Sazonov S.V. Nonlinear propagation regimes of resonant pulses in multilevel quantum media. Opt. Spektrosk., 95, P. 666-674 (2003).

[20] Kalitkin N.N. Calculation methods. Nauka, Moscow, 1978. 512 p.