CC BY
6
www.volsu.ru
ФИЗИКА И АСТРОНОМИЯ
DOI: https://doi.Org/10.15688/mpcm.jvolsu.2020.3.4
UDC 535.3, 54.03 Submitted: 14.04.2020
LBC 22.343.4, 24.5 Accepted: 13.07.2020
MODELING OF TWO-DIMENSIONAL LIGHT BULLETS PROPAGATION IN AN ARRAY OF CARBON NANOTUBES TAKING INTO ACCOUNT THE MECHANICAL TENSION AND MAGNETIC FIELD1
Natalia N. Konobeeva
Doctor of Physical and Mathematical Sciences, Associate Professor, Department of Information Systems and Computer Modeling, Volgograd State University yana'nn@volsu.ru
https://orcid.org/0000-0002-6043-9555
Prosp. Universitetsky, 100, 400062 Volgograd, Russian Federation
о
CM
о
CS
> CO
с
о ©
Dmitry S. Skvortsov
Assistant Lecturer, Department of Information Systems and Computer Modeling,
Volgograd State University
dmitry.skvortsov@volsu.ru
https://orcid.org/0000-0001-9467-4177
Prosp. Universitetsky, 100, 400062 Volgograd, Russian Federation
CO
Q
§ Abstract. In this paper, we study the influence of acoustic and magnetic
fields on the propagation of the two-dimensional light bullet in an array of carbon nanotubes. The acoustic field is taken into account in the framework of the gauge theory. A magnetic field is applied along the nanotube axis and leads to a change in the electronic spectrum of n-electrons. It is shown that the pulse stably propagates in the medium, taking into account both of these factors. In this case, the magnetic field and tension slow down the pulse, as well as change its amplitude.
Key words: localized pulses, gauge theory, magnetic flux, acoustic field, carbon nanotubes.
Introduction
Extremely short optical pulses of femtosecond duration are of great technological and fundamental interest. They are also the subject of many recent studies of nonlinear optics [2-5].
Such short light pulses have potential applications in the spectroscopy of dielectrics, semiconductors and transient chemical processes, light propagation in the atmosphere and at the interface of materials [1].
Nonlinear optics usually works with almost harmonic electromagnetic waves, modulated by the envelope, under which there are about 100-1 000 oscillations. In contrast, an extremely short pulse usually contain only a few oscillations of the electromagnetic field. Moreover, light bullets are of great interest among such pulses. They are localized in all directions and propagating in a nonlinear medium at a constant speed [6; 7].
It should be noted, that the successes of modern laser technologies in terms of the formation of powerful electromagnetic radiation with given characteristics, including light bullets, are an incentive for a comprehensive study of such pulses in various media, including carbon nanotubes [8]. An important issue is the study of the pulse behavior in a medium under the action of strong external fields. The authors of [9] considered the influence of external strains in the one-dimensional case.
In this paper, we study the dynamics of two-dimensional light bullet in a medium with CNTs at the presence of an acoustic and a constant magnetic field.
1. Model and basic equations
The electron spectrum for zig-zag carbon nanotubes (n, 0) in the presence of a magnetic field directed along the CNT axis has the form:
z(vz, s) = Y 1 + 4cos(ap„) cos ( П (s + ф) ) + 4cos2 ( П (s + ф) ), (1)
у \n Ф0 / \n Фо у
where y « 2.7 eV; a = 3b/2h; b = 0,142 nm is the distance between the adjacent carbon atoms; (py, s) is the quasi-momentum, where py is the pulse component along the CNT axis; s = 0, ..,n; $ is the magnetic flux through the CNT cross section; $0 = ^f is the magnetic flux quantum [10], e is the electron charge, c is the light velocity.
We consider the interaction between carbon nanotubes to be weak. And it is not taken into account here. The electric field E is directed along the nanotubes axis (OY), and the electromagnetic pulse moves in the transverse direction (OZ).
The acoustic field is taken into account in the framework of the gauge theory and is associated with the appearance of a stress field caused by the strain field. This field is determined by the corresponding vector potential A', which changes the momentum of electrons in the medium containing the array of carbon nanotubes. In deformed CNTs, all interatomic bonds turn out to be nonequivalent, and the three hopping integrals may not be equal to each other. Since CNTs are a rolled sheet of graphene (the ZOY plane) in the cylinder along the OY axis, we write all the equations for graphene, and then take into account the periodic boundary conditions. In accordance with this, the calibration vector potential has the form: A' = (A'z,A'y), whose components can be written in the form [11]:
\/3 1
A'z = ^ (Y3 = Y2), K = 2 (Y3 + Y2 - 2Yi), (2)
where Yi, y2, Y3 are the hopping integrals.
In the case of weak deformation of the crystal lattice, each hopping integral can be expanded in a series up to the second term [12]:
Yi = Y + r (u¿ - uo), (3)
a2
here y is the unperturbed hopping integral; r¿ is the radius vector of nearest neighbors; Uj is the displacement vector of the z-th atom; u0 is the displacement vector of the central atom; |3 is the Gruneisen electronic parameter, which expresses a change in the frequency of lattice vibrations depending on changes in the volume of the system [13]. For carbon structures it is equal 2.
In the continuum limit, one can expand the displacement of carbon atoms:
Uj - uo a (r¿V)u, (4)
u is the displacement field.
Substituting (3) into (2) taking into account (4), we obtain:
A'z = X pY(Uzz — Uyy), A'y = -x ^^ Uzy, (5)
where x is the parameter depending on the characteristics of the chemical bond in the substance, which can be set for CNTs equal to 1 [14].
Thus, the strain gauge field is proportional to the strain tensor which determines the local strain of the lattice. By setting the deformation field u, we obtain that the strain tensor in the linear approximation has a well-known form [15]:
9aUp + dpUa ... =-2-• (6)
The calibration potential changes the band structure by calibrating the pulse by —qA/c in the same way as the external electromagnetic field. The contributions of the electromagnetic field and lattice deformations are reduced to the sum of the corresponding vector potentials. The model does not take into account the rotations of interatomic bonds. Based on these assumptions, only one correction Az is taken into account in the effective equation for the vector potential.
The longitudinal and transverse components of the strain tensor are related by the known relation [15]: uzz = —v-uyy, h = 0,19 is the Poisson's ratio for CNTs [16].
In the two-dimensional case, the equation describing the propagation of the light bullet in the array of carbon nanotubes takes the form:
d2A d2A £ d2A A ... n
+ M — ^ ^ — + A'/e) = °, (7)
* = x ^ (1 + h),
here £ is the dielectric constant of the medium; the current density j is determined by the formula:
n
3 = 2£ v(pz,s) • f (pz ,s) dp, (8)
S=1BZ
where v(p,s) = de(py,s)/dpy is the electron velocity; f (py,s) is the Fermi distribution function; z(py, s) is determined by the formula (1), integration is carried out in the first Brillouin zone (BZ). We omit the calculation of current density in this paper, because it repeats the standard calculations for CNTs, and we present only the final form of the effective equation:
д2 А д2 A e д2 A 4nen0ya
■X „ . t a^A + A') 4
g 8m I-с-)
дх2 dz2 с2 dt2
0,
(9)
where
¡ d(py a) • cos(q ару) • exp - cos(qapy) • bsq/квT
Fq = -q
s=1 — п
9=1
-)
n п ix
^ / d(py a) • exM - £ cos(qapy) • bsq/квT
s=1 — п \ q=1
)
п/а
п
d(py a) • e(py, s) • (qapy);
—п/а
n0 is the electron concentration; ks is the Boltzmann constant; T is the temperature; bsq are the coefficients in the expansion of the electron dispersion law (1) in a Fourier series. We make the dimensionlessness of the quantities included in the equation (9):
~ eaA aA x z ~ c • t c
A =-, =-, x = —, z = —, i =-, a0 = —-— .
c c a0 a0 a0 nn0y0
The initial condition is chosen in the form:
À(x, z, 0) = Q • ex^-- t")
dÄ(x,z, 0) 2z vQ
dt
II
• exp
f t2 \ V ^x IV
(10)
vz/с.
Where Q is the amplitude of the electromagnetic pulse at the entrance to the environment with CNT; lx, lz are determined the pulse width along the x and z axis, vz is the pulse velocity along z-axis.
Next, numerical simulation is carried out using an explicit second-order finite-difference scheme [17] with zero Neumann boundary conditions for the following system parameters: CNTs of the zig-zag type (7, 0), e = 4, n0 = equation (9) are calculated at a temperature T
10 m 3, the coefficients Fq determined by
77 K [18], v = 0.9. The relaxation time in CNTs is trei ~ 10-12 s, the pulse duration ipuise « 2 ■ 10-15 s, thus ipuise ^ irei.
The evolution of the electromagnetic pulse is presented in figure intensity I, which can be determined as I = c-2(dA/dt)2.
We draw the pulse
a
Fig. 1. Evolution of light bullet (uzz = 0,1, = n/2): a) i = 0; b) i = 3, 5; c) i = 6, 5; d) i = 9, 5.
Imax is the maximum of the intensity for time value i
From figure 1 it can be concluded that the light bullet is stable in the CNT array under the influence of acoustic and magnetic fields, since the pulse travels about 10 of its wavelengths. Note that the pulse also exibits a slight diffraction spreading. Let us study the influence of the magnetic field on the light bullet propagation. A two-dimensional view is shown in figure 2, a longitudinal section (x = 0) is in figure 3.
Fig. 2. The intensity of the electric field of the light bullet (i = 9, 5): a) uzz = 0, $/$0 = 0. Figures (b-e) for uzz = 0,1: b) $/$0 = 0; c) $/$o = n/6; d) $/$o = n/3; e) $/$o = n/2. !max is the maximum of the intensity for time value i
Figure 2a shows the pulse shape without taking into account the magnetic field and mechanical deformation. Figures 2b-e are plotted in the presence of mechanical tension and for different values of magnetic flux. It cab be seen when the magnetic field increases, the pulse slows down. Also for figure 2d the appearance of a "tail" behind the main pulse is observed as before [8].
The given dependence shows that, controlling the magnitude of the magnetic field, we can change not only the shape, but also the amplitude of the pulse.
The dependence of the intensity of the electric field on the acoustic field magnitude is shown in figure 4.
It can be seen, that in the presence of a magnetic field, but at zero mechanical tension, the pulse retains its localization region, but its amplitude increases. When the tension increases (fig. 3c, d) we observe that the amplitude decreases. Moreover, the larger the
value of uzz (fig. 4d), the closer the pulse shape to the initial one (fig. 4a).
Fig. 3. The intensity of the electric field of the light bullet, a longitudinal section at x = 0 (uzz = 0,1): a) pulse at the initial time moment;
Others for i = 9, 5: b) = 0; c) $/$0 = n/6, for clarity, the curve is raised by 1 unit along the vertical axis;
d) $/$0 = n/3 the curve is raised by 2 units along the vertical axis;
e) $/$0 = n/2 the curve is raised by 3 units along the vertical axis.
I0 is the maximum of the intensity for time value i = 0
I/Ir
1
0.8 0.6 0.4 0.2 0
1 1 1 - — ■ — . — а Ь 1
.........с - H
J й ' J " Л
--1 hî Wi Ч—
12
13
14
15
16
Fig. 4. The intensity of the electric field of the light bullet (a longitudinal section at x = 0, i = 9, 5): a) ùzz = 0, = 0. Figures (b-d) correspond to $/$0 = n/3: b) ùzz == 0; c) uzz = 0, 05; d) uzz = 0,1. I0 is the maximum of the intensity for time value i = 0 at ùzz = 0 and $/$0 = 0
Conclusion
Key results of this work are summarized as follows:
1) The numerical modeling showed that the acoustic and magnetic fields significantly affect the dynamics of the light bullet in the carbon nanotubes array.
2) It is observed, that using the magnetic flux, we can change the pulse shape, its region of the localization and the amplitude. Thus, the magnetic field is the controlling parameter for the electromagnetic field of the light bullet.
3) The effect of deformation is manifested both in a change in the shape of the pulse and its amplitude.
NOTE
1 Konobeeva N.N. expresses gratitude to the Ministry of Science and Higher Education of the Russian Federation for funding numerical modeling under the State Assignment task (0633-2020-0003).
REFERENCES
1. Bloembergen N., Mlejnek M., Moloney J., Wright E. Femtosecond Pulses at the Boundary of a Nonlinear Dielectric. Advances in Laser Physics. Boca Raton, CRC Press, 2000, pp. 15-23.
2. Fermann M.E., Andrejco M.J., Silberberg Y., Weiner A.M. Generation of Pulses Shorter than 200 fs from a Passively Mode-Locked Er Fiber Laser. Indian Journal of Pure and Applied Mathematics, 1993, vol. 18, iss. 1, pp. 48-50. DOI: 10.1364/0L.18.000048.
3. Goorjian P.M., Agrawal G.P. Computational Modeling of Ultrashort Pulse Propagation in Semiconductor Materials. Quantum Optoelectronics, 1997, vol. 9, pp. 95-97.
4. Kaplan A.E., Shkolnikov P.L. Electromagnetic "Bubbles" and Shock Waves: Unipolar, Nonoscillating EM Solitons. Phys. Rev. Lett., 1995, vol. 75, iss. 12, pp. 2316-2319. DOI: 10.1103/PhysRevLett.75.2316.
5. Kozlov S.A., Samartsev V.V. Fundamentals of Femtosecond Optics. Cambridge, Woodhead Publishing, 2013. 273 p.
6. Silberberg Ya. Collapse of Optical Pulses. Optics Letters, 1990, vol. 15, iss. 22, pp. 12821284. DOI: 10.1364/0L.15.001282.
7. Sazonov S.V. Light bullets. Journal of Physics: Conference Series. Bristol, IOP Publishing, 2016, vol. 714, pp. 15-23. DOI: 10.1088/1742-6596/714/1/012015.
8. Belonenko M.B., Lebedev N.G., Popov A.S. Two-Dimensional Light Bullets in an Array of Carbon Nanotubes. JETP Letters, 2010, vol. 91, iss. 9, pp. 461-465. DOI: 10.1134/S0021364010090067.
9. Lyapkosova O.S., Lebedev N.G., Belonenko M.B. Extremely Short Optical Pulses in Strained Graphene in Terms of the Gauge Theory. Physics of the Solid State, 2013, vol. 55, pp. 2602-2607. DOI: 10.1134/S1063783413120226.
10. Ovchinnikov A.A., Atrazhev V.V. Magnetic Susceptibility of Multilayered Carbon Nanotubes. Physics of the Solid State, 1998, vol. 40, pp. 1769-1773. DOI: 10.1134/1.1130653.
11. Vozmediano M.A.H., Katsnelson M.I., Guinea F. Gauge Fields in Graphenes. Physics Reports, 2010, vol. 496, iss. 4-5, pp. 109-148. DOI: 10.1016/j.physrep.2010.07.003.
12. Suzuura H., Ando T. Phonons and Electron-Phonon Scattering in Carbon Nanotubes. Physical Review B, 2002, vol. 65, iss. 23, pp. 235412-235431. DOI: 10.1103/PhysRevB.65.235412.
13. Sanditov D.S., Mantatov V.V., Darmaev M.V., Sanditov B.D. On the Gruneisen Parameter for Crystals and Glasses. Technical Physics, 2009, vol. 54, pp. 385-388. DOI: 10.1134/S1063784209030098.
14. Manes J.L. Symmetry-Based Approach to Electron-Phonon Interactions in Graphene. Physical Review B, 2007, vol. 76, iss. 4, pp. 045430-045440. DOI: 10.1103/PhysRevB.76.045430.
15. Landau L.D., Lifshitz E.M. Theory of Elasticity. Oxford, Pergamon Press, 1970. 248 p.
16. Rafii-Tabar H. Computational Physics of Carbon Nanotubes. Contemporary Physics, 2010, vol. 51, iss. 1, pp. 92-93. DOI: 10.1080/00107510903020519.
17. Strikwerda J. Finite Difference Schemes and PDE. Pacific Grove, Wadsworth & Brooks/Cole Advanced Books & Software, 1989. 435 p.
18. Zhukov A.V., Bouffanais R., Malomed B.A., Leblond H., Mihalache D., Fedorov E.G., Rosanov N.N., Belonenko M.B. Collisions of Three-Dimensional Bipolar Optical Solitons in an Array of Carbon Nanotubes. Phys. Rev. A, 2016, vol. 94, iss. 5, pp. 053823-053836. DOI: 10.1103/PhysRevA.94.053823.
МОДЕЛИРОВАНИЕ РАСПРОСТРАНЕНИЯ ДВУМЕРНОЙ СВЕТОВОЙ ПУЛИ В МАССИВЕ УГЛЕРОДНЫХ НАНОТРУБОК С УЧЕТОМ МЕХАНИЧЕСКОЙ НАГРУЗКИ И МАГНИТНОГО ПОЛЯ
Наталия Николаевна Конобеева
Доктор физико-математических наук,
доцент кафедры информационных систем и компьютерного моделирования,
Волгоградский государственный университет
yana_nn@volsu.ru
https://orcid.org/0000-0002-6043-9555
просп. Университетский, 100, 400062 г. Волгоград, Российская Федерация
Дмитрий Сергеевич Скворцов
Ассистент кафедры информационных систем и компьютерного моделирования, Волгоградский государственный университет dmitry.skvortsov@volsu.ru https://orcid.org/0000-0001-9467-4177
просп. Университетский, 100, 400062 г. Волгоград, Российская Федерация
Аннотация. В данной работе мы исследуем влияние акустического и магнитного полей на процесс распространения локализованного в двух направлениях предельно короткого оптического импульса в среде, содержащей массив полупроводниковых углеродных нанотрубок zig-zag типа. Мы считаем, что расстояние между соседними нанотрубками существенно превышает их диаметр, то есть взаимодействие между углеродными нанотрубками считается слабым и можно им пренебречь. Учет акустического поля проводился в рамках калибровочной теории. Появление этого поля обусловлено появлением поля напряжений, которое вызывается механическими нагрузками на кристаллическую решетку углеродных нанотрубок. Магнитное поле прикладывается вдоль оси нанотрубок и приводит к изменению электронного спектра п-электронов. Показано, что импульс стабильно распространяется в среде с
учетом обоих этих факторов. При этом магнитное поле и деформация приводят к замедлению импульса, а также к изменению его амплитуды. Также наблюдается появление «хвоста» за основным импульсом
Ключевые слова: локализованные импульсы, калибровочная инвариантность, магнитный поток, акустическое поле, углеродные нанотрубки.
