NANOSYSTEMS: Rakhmanov S.Z., et al. Nanosystems:
PHYSICS, CHEMISTRY, MATHEMATICS Phys. Chem. Math., 2023,14 (2), 164-171.
http://nanojournal.ifmo.ru
Original article DOI 10.17586/2220-8054-2023-14-2-164-171
Optical high harmonic generation in a quantum graph
SaparboyZ. Rakhmanov1", Ikhvoliddin B. Tursunov2'6, Khikmatjon Sh. Matyokubov3c, Davron U. Matrasulov4,d
1Chirchik State Pedagogical University, Chirchik, Uzbekistan
2 National Universty of Uzbekistan, Vuzgorodok, Tashkent, Uzbekistan
3 Urgench State University, Urgench, Uzbekistan
4 Turin Polytechnic University in Tashkent, Tashkent, Uzbekistan
asaparboy92@gmail.com, 6ixvoliddin.tursunov.1998.04.07@gmail.com, chikmat0188@mail.ru, ddmatrasulov@gmail.com
Corresponding author: Saparboy Z. Rakhmanov, saparboy92@gmail.com
Abstract High ordered harmonic generation in a quantum graph is studied by considering quantum star graph interacting with external monochromatic optical field. Using the numerically obtained solutions of the time-dependent Schrodinger equation on quantum graph, main characteristics of high harmonic generation are computed. In particular, time-dependence of the average dipole moment and high harmonic generation spectra, determined as the generated field intensity as a function of harmonic order are analyzed. Extension of the proposed model to the case of other graph topologies and application to the problem of tunable high harmonic generation are discussed.
Keywords quantum graphs, optical high harmonic generation, average dipole moment, harmonic generation spectra, quantum graph in electromagnetic field.
Acknowledgements This work is supported by the grant of the Ministry of Innovative Development of the Republic of Uzbekistan, World Bank Project "Modernizing Uzbekistan National Innovation System" (Ref. Nr. REP-05032022/235).
For citation Rakhmanov S.Z., Tursunov I.B., Matyokubov Kh.Sh., Matrasulov D.U. Optical high harmonic generation in a quantum graph. Nanosystems: Phys. Chem. Math., 2023,14 (2), 164-171.
1. Introduction
Study of the nonlinear optical phenomena occurring in quantum regime is of importance for many practically important tasks appearing in the intersection such topics as quantum optics, atom optics, and optoelectronics. From the
viewpoint of fundamental research, the importance of such studies is related to attosecond physics and the physics of ultrafast phenomena, while the practical importance is caused by the relevance to high-power laser generation, optical materials design and optoelectronic device fabrication. An interesting aspect of this topic is optical harmonic generation in quantum regime, which attracted much attention recently [1,2] in the context of atomic physics and some confined low-dimensional quantum systems, such as quantum wells and dots. One of the main tasks in this field is achieving
slowly-decaying (as a function of harmonic order) harmonic generation intensity. Solving of such problem is complicated due to the typical features of the harmonic emission spectra of an atom in a strong optical field, which are known as
"the plateau" and "the cutoff". These latter are a wide frequency region with harmonics of comparable intensities, and an abrupt intensity decrease at the end of the plateau. Physical mechanisms of such effects have been explained within the so-called "three-step" model. Existence of such effect makes difficult generating of very high order harmonics and ultrashort pulses, as their intensity becomes very small at high harmonic orders. Therefore, revealing of the high harmonic generation (HHG) regime, where the intensity would slowly decay as a function of harmonic order is of importance for different practical tasks. Study of HHG in low-dimensional quantum systems can be one of the ways providing achieving of such a goal. In this paper, we address the problem of optical high harmonic generation in branched quantum wires. These latter are modeled in terms of the quantum graphs, which are determined as a set of wires connected to each other according to a rule, called topology of the graph [3-5]. Generation of the topical harmonics is assumed as to be caused by the interaction of branched quantum wires with an external monochromatic optical field. Using the eigenfunctions of the quantum graph, we compute numerically the average dipole moment as a function of generated frequencies. Also, as the main characteristics of the high harmonic generation process, we compute the spectrum of the high harmonic generation which is determined as the intensity of the field generated by the interaction of the branched quantum wire with the external optical field. Earlier, quantum graphs have attracted much attention in different contexts for modeling particle and quasiparticle dynamics in branched structures (see Refs. [3-16]).
© Rakhmanov S.Z., Tursunov I.B., Matyokubov Kh.Sh., Matrasulov D.U., 2023
2. Schrodinger equation on graphs
Here we briefly recall basic quantum graphs theory following Refs. [4,10]. Graphs are the systems consisting of vertices which are connected by edges. The edges are connected according to a rule which is called the topology of the graph. Topology of the graph is described in terms of the adjacency matrix [4,5]:
1, if i and j are connected; 0, otherwise;
i, j
1,2,..., N.
Earlier, quantum graphs were extensively studied in the context of quantum chaos theory [4-6]. Strict mathematical formulation of the boundary conditions was given by Kostrykin and Schrader [3]. Inverse problems on quantum graphs have been studied in Refs. [7-9]. An experimental implementation of quantum graphs on (optical) microwave waveguide networks is discussed in the Ref. [8]. Despite the fact that different issues of quantum graphs and their applications have been discussed in the literature, the problem of driven graphs is still remaining as less-studied topic.
Quantum particle dynamics on a graph is described by one-dimensional multi-component Schrodinger equation [4,5] (in the units h = 2m =1):
d2^(x)
dx2
k2^6(x), b = (i,j),
(1)
where b denotes a bond connecting i-th and j-th vertices, and for each bond b, the component of the total wave function ^ is a solution of Eq. (1). In Eq. (1) components are related through the boundary conditions, providing continuity and current conservation [4]:
• Continuity,
|x=0 = y-, lx=Lij = y, for all i < j and C-,, = 0.
• Current conservation,
(2)
- E Ci,.
j<i
dx
+ E Ci,.
(x)
j>i
dx
= A-y-
Here, the parameters Ai are free parameters which determine the type of boundary conditions. The Dirichlet boundary conditions correspond to the case when all the Ai = œ. Solution of Eq. (1) obeying the above boundary conditions can be written as
C- •
¿i j =-i^t— (wi sin k(Li j — x) + y, sin kx),
sin kLi,j
where the quantities yi are the solutions of the algebraic system following from the continuity conditions:
fcC,
sin kLi
-(-y- cos kLi,j + y, ) = A-y-.
The eigenvalues of Eq. (1) can be found from the spectral equation
det(hi,j (k)) = 0
(3)
where
hi,j (k)
^ ^ Ci,m cot kLi,m
m=i
Ci,j sin kLi,j,
k
i = j
i = j
Here we focus on the simplest graph topology, star graph. In this case, since the graph has only single non-boundary vertex, unlike the notations in Eqs. (1) and (2), indices are assigned to the bonds (not to the vertices), i.e. Lj means j-th bond. Thus, in special case of star graph, the boundary conditions can be written as [10]
^l|y=0 = ^2 |y=0
^N |y=0,
^1|y=Li = ^2|y=L2 = ... = ^N |y=LN N
(4)
j=1
Ed
dy^j |y=0 = 0.
The eigenvalues of star graph can be found from the following spectral equation [10]
N
^tan-1(k„LT ) = 0. j=i
L
x
x
0
Fig. 1. Sketch of a metric star graph. Lj is the length of the j th bond with j=1, 2,. ..,N
Bond 1
Bond 2
0 200
400 600 t
0
-5
s"10
TÏ
v -15 -20 -25
800 1000
0 200
Bond 3
0 200
400 600 t
400 600 t
800 1000
800 1000
Fig. 2. Average dipole moment as a function of time on star graph at different values of a parameter for external field's amplitude, e = 0.1 and frequency, w0 = 0.01
Corresponding eigenfunctions are given as [10]
&j,n =
where
Bn
sin( knLj )
sin(k„(Lj - y)),
(6)
Bn =
Lj +sin(2knLj )
j sin2 (knLj )
E
The functions n are orthonormal, i.e., fulfill the following condition:
N j
/ ^*j,m(X)^j,n(X)d'X = Smn. j = 1 0
In the next section, we use these eigenfunctions for computing the high harmonic generation spectrum.
2
0 -10 -20 > -30 -40 -50 -60
Bond 1
Bond 2
200
--20
-30
-40
-20
--40
-60
-80
400 600 800 1000
-
-10
0 200 400 600
-
Bond 3
0 200 400 600 800 1000
-
800
1000
Fig. 3. Average dipole moment as a function of time on star graph at different values of amplitude of the external field for field's frequency, w0 = 0.01, for length of each bond, L1 = 200.31, L2 = 205.53 and L3 = 210.18
3. Driven star graph by external potential
Consider a quantum star graph (Y-junction) interacting with an external linearly polarized monochromatic electromagnetic field. Such system can be described by the following time-dependent (multi-component) Schrodinger equation:
.<9^, (x,t) i dt
d2
——tt + ex cos W0Î
dx2
j = 1, 2, 3,
where e and w0 are the amplitude and the frequency of the optical field, respectively.
The solution of Eq. (8) can be written in terms of the complete set of solutions of Eq. (1)
(x,t) = E Cn(i) jn(x),
(8)
(9)
where ^¿,n(x) are given by Eq. (6).
Substituting (9) into Eq. (8) and using the condition given by Eq. (7), we get a system of first order ordinary differential equations which has the form
iCn(t) = £nCn(t) + E Cm(i)Vnm, (10)
where
Km = Ej j^(i)^j,mdx = e cos ^01 Ej xj^j,mdx
j 0 j 0 and e„ are the eigenvalues of the unperturbed system.
Using numerically found solutions of Eq. (8), one can compute physically observable characteristics of the system "quantum star graph + external optical field".
In numerical calculations, we choose the Gaussian wave packet given on the first bond of the graph as the initial state:
1 (x-M)2
0)=$(x) = e---T-, (11)
with ct being the width of the packet. For other bonds, the initial wave function is assumed to be zero, i.e. 0) 0) =0.
0
0
0
Bond 1
Bond 2
100
10"-
15 16 17 18 19 20
50
10
0 10
10°
I10-5
10"'
0 10
20 30
Harmonic order
Bond 3
20 30
Harmonic order
40 50
40 50
Fig. 4. Spectrum of harmonic generation on star graph at different values of a parameter for external field's amplitude, e = 0.1, and frequency, = 0.01, for duration of interaction, T = 1000
100
10"-
Bond 1
10
20 30
Harmonic order
10°
10
100
10"-
40 50
0 10
Bond 3
Bond 2
20 30
Harmonic order
0 10 20 30 40 50
Harmonic order
40 50
Fig . 5. Spectrum of harmonic generation on star graph at different values of amplitude of external field for field's frequency, w0 = 0.01, for length of each bond, Li = 200.31, L2 = 205.53 and L3 = 210.18, for duration of interaction, T = 1000
For such initial conditions, the expansion coefficients at t = 0 can be written as
11
Cn(0) = J $(x)#n(x)dx.
0
4. High harmonic generation in star graph
Consider high harmonic generation caused by the interaction of the quantum star graph with the external monochromatic optical field. One of the main characteristics of such process is the average (expectation) value of the dipole moment which is determined as [1]
(dj(t)> = -(tfj(x,t)|x|tfj(x,t)>, j = 1, 2, 3. To compute (dj (t)>, one can use the solution of Eq. (8). In terms of the expansion coefficients in Eq. (9), one can write the average dipole moment as
(dj (t)> = - E C*m(t)Cn(t) j
where Vjm„ = j xj^j,mdx.
0
In the following, all the plots are produced by choosing the initial state as Gaussian wave packet determined in Eq. (11) and coming from the first bond. Furthermore, we choose the bond lengths as Lj = aj. The width and the initial position of the packet are fixed as a = 10, ^ = 50, respectively. Fig. 2 presents the plots of the average dipole moment for star graph with the bond lengths, h = 200.31, l2 = 205.53 and l3 = 210.18, as a function of time at different values of a. The amplitude and the frequency of the electromagnetic field are taken as e = 0.1 and w0 = 0.01, respectively. As it can be seen from this plot, (dj(t)> is quasi periodic in time and the period decreases as the parameter a increases.
In Fig. 3, plots show time-dependence of the average dipole moment of each bond at different values of amplitude of the external field strength, e = 0.01 (blue line), e = 0.1 (red line) and e =1 (yellow line). One can see the strong dependence on the amplitude of the electromagnetic field. Period of the quasi periodic average dipole moment decreases as the amplitude of the electromagnetic field grows.
Another important characteristics of the high harmonic generation in quantum regime is its spectrum which is determined in terms of its intensity on each bond (as a function of frequency generated) and given as
1 T
Ij H = |dj HI2 =
e-^(dj (t)>dt
(12)
where T is the total duration of interaction of the quantum graph with the external optical field. The total harmonic generation intensity is defined as the sum of those on each bond. The main practically relevant feature of the high harmonic generation intensity is the appearance of "the plateaus" at certain values of the generated frequency. In other words, the high harmonic generation spectrum, which is determined as the ratio of the external field frequency and the frquency of generated one, N = w/w0, consists of a plateau where the harmonic intensity is nearly constant over many orders and a sharp cutoff. As higher the number of plateaus, as attractive the harmonic generation from the viewpoint of attosecond physics and high power laser generation [17,18]. Fig. 4 presents the plot of I(N) on each bond of the quantum star graph for different values of a choosing the field parameters as e = 0.1, w0 = 0.01 and T = 1000. In Fig. 4, one can see the decreasing of the spectrum as harmonic order increases. Moreover, for large values of the length of bonds, the intensity of high order harmonics increases. Existence of "quasi-plateau" can be seen from the inset. Fig. 5 shows plots of the spectrum of the harmonic generation at different values of the external field strength and field's frequency, w0 = 0.01 for bond length chosen as Lx = 200.31, L2 = 205.53 and L3 = 210.18 for duration of interaction, T = 1000. The plot shows that the HHG intensity is strongly depends on the field strength. For higher values of e, one can observe increasing of the intensity. The above study can be directly extended to any branching topology of a quantum graph (e.g., tree, hexagon, loop, H-graph, etc). Of course, complicated branching topologies should provide wider opportunity for the HHG-process.
5. Conclusion
In this paper, we studied high harmonic generation in quantum star graph caused by its interaction with the external monochromatic optical field. The average dipole moment as a function of time is analyzed as one of the characteristics of HHG. The harmonic generation spectrum described in terms of the generation intensity as a function of the harmonic order given by the ratio of the generated harmonics to that of the external field, is computed at different values of the external field strength, as well as for different bond lengths. Appearance of narrow plateau was shown by analyzing the plots of the HHG spectra. The model considered in this study can be implemented, e.g., in different quantum wire networks fabricated on the basis of superconductors, semiconductors or other low-dimensional quantum materials. Some versions of such quantum wire networks have been studied earlier in Refs. [19-22]. Extension of the above model to the case of other branching topologies and external (e.g., bi-, or poly-chromatic, or nonlinearly polarized) fields can provide an effective tool for tunable high harmonic generation in low-dimensional structures. In addition, the model can be considered as a version of the quantum graph based laser concept. Such concepts have been earlier discussed in [23-43].
References
[1] Boyd R.W. Nonlinear Optics, 2007, 3rd ed, Academic Press.
[2] Tong X., Chu Sh. Theoretical study of multiple high-order harmonic generation by intense ultrashort pulsed laser fields: A new generalized pseudospectral time-dependent method. Chem. Phys. B, 1997, 217(2-3), P. 119-130.
[3] Kostrykin R., Schrader R. Kirchhoff's rule for quantum wires. J. Phys. A: Math. Gen, 1999, 32, P. 595-630.
[4] Kottos T., Smilansky U. Periodic Orbit Theory and Spectral Statistics for Quantum Graphs. Ann. Phys., 1999, 274, P. 76-124.
[5] Gnutzmann S. and Smilansky U. Quantum graphs: Applications to quantum chaos and universal spectral statistics. Adv. Phys., 2006, 55, P. 527625.
[6] Gnutzmann S., Keating J.P., Piotet F. Eigenfunction statistics on quantum graphs. Ann. Phys., 2010, 325, P. 2595-2640.
[7] Boman J., Kurasov P. Symmetries of quantum graphs and the inverse scattering problem. Adv. Appl. Math., 2005, 35, P. 58-70.
[8] Hul O., Bauch S., Pakonski P., Savytskyy N., Zyczkowski K. and Sirko L. Experimental simulation of quantum graphs by microwave networks. Phys. Rev. E., 2004, 69, P. 056205/1-7.
[9] Cheon T., Exner P., Turek O. Approximation of a general singular vertex coupling in quantum graphs. Ann.Phys., 2010, 325, P. 548-578.
[10] Keating J.P. Fluctuation statistics for quantum star graphs. Quantum graphs and their applications. Contemp. Math., 2006, 415., P. 191-200.
[11] Matrasulov D.U., Yusupov J.R., Sabirov K.K., Sobirov Z.A. Time-dependent Quantum Graph. Nanosystems: Physics, Chemistry, Mathematics, 2015, 6(2), P. 173-181.
[12] Yusupov J.R., Sabirov K.K., Ehrhardt M., Matrasulov D.U. Transparent quantum graphs. Phys. Lett. A., 2019, 383, P. 2382-2388.
[13] Matrasulov D.U., Sabirov K.K. and Yusupov J.R. PT-symmetric quantum graphs. J. Phys. A., 2019, 52, P. 155302/1-11.
[14] Nikiforov D.S., Blinova I.V. and Popov I.Y. Schrodinger and Dirac dynamics on time-dependent quantum graph. Indian J. Phys., 2019, 93(7), P. 913-920.
[15] Yusupov J.R., Sabirov K.K., Asadov Q.U., Ehrhardt M. and Matrasulov D.U. Dirac particles in transparent quantum graphs: Tunable transport of relativistic quasiparticles in branched structures. Phys. Rev. E., 2020,101, P. 062208/1-8.
[16] Yusupov J.R., Sabirov K.K. and Matrasulov D.U. Dirac particles on periodic quantum graphs. Phys. Rev. E., 2021, 104, P. 014219/1-7.
[17] Strelkov V.V., Platonenko V.T., Sterzhantov A.F. and Ryabikin M.Yu. Attosecond electromagnetic pulses: generation, measurement, and application. Generation of high-order harmonics of an intense laser field for attosecond pulse production. Phys. Uspekhi, 2016, 59(5), P. 425-445.
[18] Winterfeldt C., Spielmann C., Gerber G. Optimal control of high-harmonic generation. Rev. Mod. Phys., 2008, 80(2), P. 117-140.
[19] Kohmoto M. Quantum wire networks with quantized Hall effect. J. Phys. Soc. Jpn., 1999, 62, P. 4001-4008.
[20] Aharony A., Entin-Wohlman O. Discrete versus Continuous Wires on Quantum Networks. J. Phys. Chem. B, 2009, 113, P. 3676-3680.
[21] Medina J., Green D., Chamon C. Networks of quantum wire junctions: A system with quantized integer Hall resistance without vanishing longitudinal resistivity Phys. Rev. A., B, 2013, 87, P. 045128.
[22] Caudrelier V., Mintchev M., Ragoucy E. Quantum wire network with magnetic flux. Phys. Lett. A., 2013, 377, P. 1788-1793.
[23] Lepri S., Trono C., Giacomelli G. Complex Active Optical Networks as a New Laser Concept. Phys.Rev.Lett., 2017, 118, P. 123901.
[24] Rotter S. Network lasers. Nature Photonics, 2019, 13, P. 140-145.
[25] Brabec T., Krausz F. Intense few-cycle laser fields: Frontiers of nonlinear optics. Rev. Mod. Phys., 2000, 72(2), P. 545-591.
[26] Yousef I., et.al. Relativistic high-power laser matter interactions. Phys. Rep., 2006, 427(2-3), P. 41-155.
[27] Winterfeldt C., Spielmann C., and Gerber G. Colloquium: Optimal control of high-harmonic generation. Rev. Mod. Phys., 2008, 80(1), P. 117-140.
[28] Krausz F., Ivanov M. Attosecond physics. Rev. Mod. Phys., 2009, 81, P. 163-234.
[29] Nisoli M., Sansone G. New frontiers in attosecond science Author links open overlay panel. Prog. Quant. Electr., 2009, 33, P. 17-59.
[30] Kohler M.C., Pfeifer T., Hatsagortsyan K.Z.,Keitel C.H. Frontiers of Atomic High-Harmonic Generation. Advances In Atomic, Molecular, and Optical Physics, 2012, 61, P. 159-208.
[31] Ahn D., Chuang S.L. Optical gain in a strained-layer quantum-well laser. IEEE J. Quantum Electron, 1988, 24, P. 2400-2406.
[32] Rosencher E., Bois Ph. Model system for optical nonlinearities: Asymmetric quantum wells. Phys. Rev. B., 1991, 44, P. 11315-11327.
[33] Milanovic V., Ikonic Z., Indjin D. Optimization of resonant intersubband nonlinear optical susceptibility in semiconductor quantum wells: The coordinate transform approach. Phys. Rev. B., 1996, 53, P. 10887-10893.
[34] Xie W., J.Lumin. The nonlinear optical rectification of a confined exciton in a quantum dot. 131. P. 943-946.
[35] Li B., Guo K-X., Zhang C.-J., Zheng Y.-B. The second-harmonic generation in parabolic quantum dots in the presence of electric and magnetic fields. Phys. Lett. A, 2007, 367, P. 493-497.
[36] Liu J-T., Su F-H., Wang H. Model of the optical Stark effect in semiconductor quantum wells: Evidence for asymmetric dressed exciton bands. Phys. Rev. B., 2009, 80, P. 113302/1-4.
[37] Shao Sh., Guo K-X., Zhang Z-H., Li N., Peng Ch. Studies on the second-harmonic generations in cubical quantum dots with applied electric field. Physica B: Physics of Condensed Matter, 2011, 406, P. 393-396.
[38] Karimi M.J., Keshavarz A., Second harmonic generation in asymmetric double semi-parabolic quantum wells: Effects of electric and magnetic fields, hydrostatic pressure and temperature. Physica E, 2012, 44, P. 1900-1904.
[39] Zhai W. A study of electric-field-induced second-harmonic generation in asymmetrical Gaussian potential quantum wells. Physica B, 2014, 454, P. 50-55.
[40] Bondarenko V., Zauzny M. Intrinsic optical intersubband bistability in quantum well structures: Role of multiple reflections. Phys. Rev. B., 2015, 91, P. 035303/1-20.
[41] Kotova L.V., Platonov A.V., Kats V.N. Optical activity of quantum wells. Phys. Rev. B., 2016, 94, P. 165309/1-5.
[42] Mohammadi S.A., Khordad R., Rezaei G. Optical properties of a semispherical quantum dot placed at the center of a cubic quantum box: Optical rectification, second and third-harmonic generations. Physica E, 2016, 76, P. 203-208.
[43] Rakhmanov S., Matrasulov D., Matveev V.I. Quantum dynamics of a hydrogen-like atom in a time-dependent box: Cooling, compressing and diffusive ionization in non-adiabatic regime. Eur. Phys. J. D., 2018, 72 P. 177/1-8.
Submitted 14 February 2023; revised 28 March 2023; accepted 29 March 2023
Information about the authors:
S. Z. Rakhmanov - Chirchik State Pedagogical University, 104 Amir Temur Str., 111700, Chirchik, Uzbekistan; ORCID 0000-0001-8569-1709; saparboy92@gmail.com
I. B. Tursunov - National Universty of Uzbekistan, Vuzgorodok, 100174, Tashkent, Uzbekistan; ixvoliddin.tursunov.1998.04.07@gmail.com
Kh. Sh. Matyokubov - Urgench State University, 14 H. Olimjon Str., 220100 Urgench, Uzbekistan; ORCID 0000-00026642-5488; hikmat0188@mail.ru
D. U. Matrasulov - Turin Polytechnic University in Tashkent, 17 Niyazov Str., 100095, Tashkent, Uzbekistan; ORCID 0000-0001-8957-0058; dmatrasulov@gmail.com
Conflict of interest: the authors declare no conflict of interest.