i l St. Petersburg Polytechnic University Journal. Physics and Mathematics. 2022 Vol. 15, No. 3.2 Научно-технические ведомости СПбГПУ. Физико-математические науки. 15 (3.2) 2022
Conference materials UDC 535.42
DOI: https://doi.org/10.18721/JPM.153.232
The comparison of the optical vortices focusing by silicon diffraction axicons and ring gratings with variable relief heights using high-performance computer systems
D. A. Savelyev , 2e 1 Samara National Research University, Samara, Russia
2 Image Processing Systems Institute of RAS - Branch of the FSRC "Crystallography and Photonics" RAS, Samara, Russia
H dmitrey.savelyev@yandex.ru
Abstract. A study was made of the spatial distribution of the intensity of the Laguerre-super-Gaussian (1,0) modes with circular, radial and azimuthal polarization depending on the change in the height of silicon subwavelength optical elements, the height of which varied from 0.2 to 3 wavelengths. Simulation by the finite difference method in the time domain showed that a change in the height of the considered optical elements significantly affects the diffraction pattern in the near zone. The smallest focal spot size was obtained for "-" circular polarization at an element height equal to two wavelengths.
Keywords: optical vortices, sub-wavelength structures, FDTD, high performance computer systems, ring gratings
Funding: This study was funded by the Ministry of Science and Higher Education, as part of the implementation of the State assignment (project no. 0777-2020-0017) in the part «Introduction», by the Samara University Development Program for 2021-2030 as part of the Priority 2030 program with the support of the Government of the Samara Region in the part «Materials and Methods», scholarship of the President of the Russian Federation (SP-1173.2022.5) in the part «Results and Discussion» and was performed under the «ERA.Net RUS plus» program and funded by RFBR, project number 20-52-76021 in the part «Conclusion».
Citation: Savelyev D. A., The comparison of the optical vortices focusing by silicon diffraction axicons and ring gratings with variable relief heights using high-performance computer systems, St. Petersburg State Polytechnical University Journal. Physics and Mathematics. 15 (3.2) (2022) 172-177. DOI: https://doi.org/10.18721/JPM.153.232
This is an open access article under the CC BY-NC 4.0 license (https://creativecommons. org/licenses/by-nc/4.0/)
Материалы конференции УДК 535.42
DOI: https://doi.org/10.18721/JPM.153.232
Сравнение фокусировки оптических вихрей кремниевыми дифракционными аксиконами и кольцевыми решетками с переменной высотой рельефа с использованием высокопроизводительных компьютерных систем
Д. А. Савельев 12е
1 Самарский национальный исследовательский университет им. академика С.П. Королева г. Самара, Россия;
2 Институт систем обработки изображений - филиал ФНИЦ «Кристаллография и фотоника» РАН, г. Самара, Россия н dmitrey.savelyev@yandex.ru
Аннотация. В данной работе исследуется пространственное распределение интенсивности мод Лагерра-супер Гаусса (1,0) с круговой, радиальной и азимутальной поляризацией в зависимости от изменения высоты кремниевых субволновых оптических элементов, высота которых менялась от 0.2 до 3 длин волн. Моделирование методом конечных разностей во временной области показало, что изменение высоты рассма-
© Savelyev D. A., 2022. Published by Peter the Great St.Petersburg Polytechnic University.
триваемых оптических элементов существенным образом влияет на дифракционную картину в ближней зоне. Наименьший размер фокального пятна был получен для «-» круговой поляризации при высоте элемента равной двум длинам волн.
Ключевые слова: оптические вихри, субволновые структуры, FDTD, высокопроизводительные компьютерные системы, кольцевые решетки
Финансирование: Работа выполнена при поддержке Министерства науки и высшего образования, в рамках выполнения Государственного задания (проект № 0777-20200017) в части «Введение», в рамках реализации Программы развития Самарского университета на 2021-2030 годы в рамках программы "Приоритет-2030" при поддержке Правительства Самарской области в части «Материалы и методы», Стипендии Президента Российской Федерации (СП-1173.2022.5) в части «Результаты и обсуждение», а также под эгидой программы «ERA.Net RUS plus» при финансовой поддержке РФФИ в рамках научного проекта № 20-52-76021 в части «Заключение».
Ссылка при цитировании: Савельев Д. А. Сравнение фокусировки оптических вихрей кремниевыми дифракционными аксиконами и кольцевыми решетками с переменной высотой рельефа с использованием высокопроизводительных компьютерных систем // Научно-технические ведомости СПбГПУ. Физико-математические науки. 2022. Т. 15. № 3.2. С. 172-177. DOI: https://doi.org/10.18721/ JPM.153.232
Статья открытого доступа, распространяемая по лицензии CC BY-NC 4.0 (https:// creativecommons.org/licenses/by-nc/4.0/)
Introduction
The structured laser beams [1], including those in the field of metamaterials and metasurfaces [2], are being actively studied in various fields of optics and photonics at present. In particular, it is known to use metasurfaces to generate an optical vortex in order to increase the throughput of optical communication systems [3]. Among structured beams, a special place is occupied by singular light beams [4], among which we should note beams with screw dislocations (optical vortices), which are currently being actively studied [5]. Optical vortices are known to be used for a number of applications, including the transmission of information over an optical fiber [6], in quantum informatics [7], and in wireless communication systems [8].
The introduction of a vortex phase singularity into the incident beam makes it possible to amplify the longitudinal component of uniformly polarized laser beams on the optical axis in the focal region [9], which makes it possible to change the diffraction pattern due to energy redistribution between the components of the electromagnetic field [10].
The presence of a strong longitudinal component in the focus region makes it possible to improve the optical resolution and is used for a number of applications, among which one can note optical manipulation and material processing [11].
Currently, the silicon and its compounds have a variety of applications [12, 13]. In particular, they are used for the manufacture of semiconductor devices [14], solar batteries [15], in medicine [16], as well as for solving photonics problems [17]. Due to the high refractive index of silicon, it is possible to obtain a high integration density [13], as well as to reduce the size of the focal spot during sharp focusing with optical microelements, including diffractive axicons [18].
The diffraction of optical vortices on silicon subwavelength mycoelements is studied in this paper: diffractive axicons and ring gratings of variable height. The influence of a change in the height of such elements on the subwavelength focusing of the considered laser beams is considered. Numerical calculations (3D) of laser radiation propagation were performed using the finite difference time domain method (FDTD) using high-performance computer systems.
Materials and Methods
This section presents the simulation parameters, types of elements, input beam: wavelength X = 1.55 ^m, the size of the computational domain for the selected radiation wave x, y, z was in the range [-5.8X; 5.8X]. The spatial sampling step is X/30, the time sampling step is X/(60c), where c is the speed of light. The thickness of the absorbing PML layer surrounding the compu-
© Савельев Д. А., 2022. Издатель: Санкт-Петербургский политехнический университет Петра Великого.
tational domain on all sides is 1.16X. The source is located inside the substrate, which occupies the entire space below the relief and is immersed in the absorbing PML layer. The refractive index of the element is n = 3.47.
Let us consider the action of the simplest implementation of an axicon in the form of a binary diffractive axicon with a numerical aperture NA = 0.95 (lattice period 1.05X), whose phase takes values of 0 and n radians. The height of the relief of a binary element, corresponding to the phase n radians, with a refractive index of the material of the element n = 3.47 is equal to:
h = -
n
0.21,
(1)
k (n -1)
where X is the wavelength of the considered laser radiation, n is the refractive index, k = 2n/X is the wave number.
An optical vortex as an input beam significantly changes the focal pattern compared to a conventional one, for example, a Gaussian beam [9], and the direction of rotation of circular polarization becomes important. Let us call "-" circular polarization the polarization in which the sign of the circular polarization is opposite to the sign of the introduced vortex phase singularity. The case when the sign of the circular polarization is co-directed with the introduced vortex phase singularity will be called "+" circular polarization.
Previously, it was shown that at the second order of the optical vortex and higher for "-" circular polarization, a shadow round light spot was formed [9, 19]. In this paper, the first order of an optical vortex in an incident beam for this type of polarization is considered.
The beams with an amplitude distribution approximated by a super-Gaussian function are often used [20, 21] for the cases where a uniform intensity distribution over the beam cross section is required.
The Laguerre-superGauss (1,0) modes, mainly with "-" circular polarization, were considered as input laser radiation (cases for "+" circular, radial and azimuthal polarizations are also given) in this paper. The amplitude of the Laguerre-superGaussian mode (1,0) is given by:
A(r, 9) = r ■ exp where g is the beam size in ^m, p = 6.
2ap
(2)
Results and Discussion
This section presents the study's results of the changing in the height relief effect of silicon microelements on the diffraction pattern in the near zone.
Figure 1 shows the results of propagation of the considered laser radiation through a diffractive axicon with a relief height h from 0.2X to 3X. The total intensity and intensity of the longitudinal
h = 0.2 Î,
fWHM = 0.57,1
FWHM-Q4U
=
FWHM = 0.46,\
FWHM = 0.42.1
h - 2X
FWHM = Ü.44.1
FWHM = 0.45,1
FWHM = W.44,1
Fig. 1. The longitudinal cross section (xz) of the Laguerre-superGauss mode (1,0) propagation («-»circular polarization, intensity): total intensity (a, b, c, d), longitudinal component intensity (e, f, g, h)
component of the electric field are shown. The size of the focal spot on the optical axis was estimated from the full width at half maximum (FWHM) of the all intensity.
All main maximums are formed outside the optical element. It should be noted that, at first, the size of the focal spot at the maximum decreases up to the height h = 2X, and then an increase is observed. At a height h = 3X, a significant shift of the maximum peak to the edge of the element relief is observed with a broadening of the focal spot size.
If estimate the value of the focal spot in the immediate vicinity of the element (at a distance of 0.1X), then an increase in intensity is observed at all heights considered (from 40% at h = 0.2X to 90% at h = 3X). The focal spot width was smaller than at the maximum and was FWHM = 0.36X for heights from 0.2X to 2X and FWHM = 0.46X for h = 3X. The width of the focal spot for the longitudinal component of the electric field was less than the total width. It should also be noted that as the height of the relief increases, a reduction in the length of the light needle is observed.
The best value for the focal spot width was obtained for the height h = 2X (total intensity). For this case, we will also consider other polarizations ("+" circular, radial and azimuthal). The results are shown in figure 2.
Fig. 2. The longitudinal cross section (xz) of the Laguerre-superGauss mode (1,0) propagation, intensity: «+ »circular polarization (a, d), radial polarization (b, e), azimuthal polarization (c)
As previously shown [18], the following is true for optical vortices of the first order: a zero value of the central focal spot indicates "+" circular polarization, and a nonzero one indicates "-" circular polarization. Similarly, for radial polarization there will be a zero value at the central focal point, and a non-zero value means azimuthal polarization.
Actually, you can see confirmation of this fact (Fig. 2). The rings obtained in the case of "+" circular polarization are smaller in length along the axis and the central shadow spot is larger than in the case of radial polarization (the length along the optical axis is greater). For the case of azimuthal polarization, the formation of a light needle on the optical axis is observed, but the intensity value is only 57% of the maximum. The size of the focal spot at the maximum on the optical axis is 0.57X and is formed at a distance of 3.66X from the edge of the relief.
Previous studies have shown that under certain conditions, a conventional round protrusion can be used to recognize polarization and focusing, and the focusing will be better than for a diffractive axicon with a similar grating size [18]. But the relief height also plays a key role [21].
In further studies, the central part of the relief changed, as shown in figure 3 ("-" circular polarization). Let's start from the height of the central cylinder hl = X, the height of the remaining zones is the same and equal to 0.2X. Next, the height h has been increased to 2X, the height of the next zone will be equal to X, and the rest by 0.2X (Fig. 3, b, e). And finally, hx = 3X, the height of the first ring is 2X, the height of the second ring from the center is X, and the rest are 0.2X each (Fig. 3, c, f).
hi - x
FWHM = 0.44,\
I
hi -2X
FWHM. = 0.43,\
FWHM = 0.45,1
hi = 31
FWHM »3.5 A
FWHM=0.38h
FWHM = 0.4 A
Fig. 3. The longitudinal cross section (xz) of the Laguerre-super-Gauss mode (1,0) propagation on ring gratings with variable relief height: hx = X (a, d), hx = 2X (b, e), hx = 3X (c, f)
It should be noted that only for the case hl = X, the main maximum is formed outside of the element (at a distance of 0.76X from the edge of the element), so all FWHM values (Fig. 3) are given on the 0.1X from the edge of the central cylinder. That is, increasing the height in the selected format leads to the formation of an intensity peak inside the element. In the general case, a broadening of the focal spot size is observed outside the element; nevertheless, at hl = 3X for the longitudinal component of the electric field, the result improved by 13.6% at a comparable intensity value.
Conclusion
The simulation by the finite difference time domain method showed that a change in the height of subwavelength silicon elements significantly affects the diffraction pattern in the near zone. Diffractive axicons and annular gratings with different heights of individual grating rings were considered as such elements. The relief height of optical elements varied from 0.2X to 3X.
The Laguerre-superGaussian (1,0) modes with circular, radial, and azimuthal polarizations were considered as input laser radiation.
An analysis of the electric field intensity pattern showed that the smallest focal spot size was obtained for the diffractive axicon at h = 2X: FWHM = 0.44X (for the maximum) and 0.36X (in the immediate vicinity of the element, at 65% of the total intensity). All main maxima for axicons were formed outside the optical elements. It should also be noted that as the height of the relief increases, a reduction in the length of the light needle of radiation is observed.
The results of numerical simulation for ring gratings with different heights showed that, in the general case, an increase in height leads to the formation of a maximum inside the element, however, at hl = 3X for the longitudinal component of the electric field, the result improved by 13.6% at a comparable intensity value (outside, near the edge of the element).
REFERENCES
1. Yan L., Kristensen P., Ramachandran S., Vortex fibers for STED microscopy, Apl Photonics. 4 (2) (2019) 022903.
2. Zhao Y., Askarpour A. N., Sun L., Shi J., Li X., Alu A., Chirality detection of enantiomers using twisted optical metamaterials, Nature communications. 8 (1) (2017) 1—8.
3. Ma X. et al., A planar chiral meta-surface for optical vortex generation and focusing, Scientific Reports. 5 (1) (2015) 1-7.
4. Cozzolino D., et al., Orbital angular momentum states enabling fiber-based high-dimensional quantum communication, Physical Review Applied. 11 (6) (2019) 064058.
5. Tkachenko G., Chen M., Dholakia K., Mazilu M., Is it possible to create a perfect fractional vortex beam, Optica. 4 (3) (2017) 330-333.
6. Pryamikov A., Alagashev G., Falkovich G, Turitsyn S., Light transport and vortex-supported wave-guiding in micro-structured optical fibres, Scientific Reports. (2020) 10 (1) 1-12.
7. Li S., Pan X., Ren Y., Liu H., Yu S., Jing J., Deterministic generation of orbital-angular-momentum multiplexed tripartite entanglement, Physical Review Letters. 124 (8) (2020) 083605.
8. Li S., et al., Efficient optical angular momentum manipulation for compact multiplexing and demultiplexing using a dielectric metasurface, Advanced Optical Materials. 8 (8) (2020) 1901666.
9. Savelyev D. A., Khonina S. N., Characteristics of sharp focusing of vortex Laguerre-Gaussian beams, Computer Optics. 39 (5) (2015) 654-662.
10. Khonina S. N., Alferov S. V., Karpeev S. V., Strengthening the longitudinal component of the sharply focused electric field by means of higher-order laser beams, Optics Letters. 38 (17) (2013) 3223.
11. Zhan, Q., Cylindrical vector beams: from mathematical concepts to applications, Advances in Optics and Photonics. 1 (2009) 1-57.
12. Cheng L., Mao S., Li Z., Han Y., Fu H. Y., Grating couplers on silicon photonics: design principles, emerging trends and practical issues, Micromachines. 11 (7) (2020) 666.
13. Mu X., Wu S., Cheng L., Fu H. Y., Edge couplers in silicon photonic integrated circuits: A review, Applied Sciences. 10 (4) (2020) 1538.
14. Myny K., The development of flexible integrated circuits based on thin-film transistors, Nat. Electron. 1 (2018) 30-39.
15. Garnett E., Yang P., Light trapping in silicon nanowire solar cells, Nano Lett. 10 (2010) 1082-1087.
16. Xu Y., et al., Silicon-based sensors for biomedical applications: a review, Sensors. 19(13) (2019) 2908.
17. Wu S., Mu X., Cheng L., Mao S., Fu H. Y., State-of-the-art and perspectives on silicon waveguide crossings: A review, Micromachines. 11 (3) (2020) 326.
18. Savelyev D., Kazanskiy N., Near-field vortex beams diffraction on surface micro-defects and diffractive axicons for polarization state recognition, Sensors. 21 (6) (2021) 1973.
19. Savelyev D. A., Khonina S. N., Golub I., Tight focusing of higher orders Laguerre-Gaussian modes, AIP Conference Proceedings. 1724 (2016) 020021.
20. Ding X., Ren Y. and Lu R., Shaping super-Gaussian beam through digital micro-mirror device, Science China Physics, Mechanics & Astronomy. 58(3) (2015) 1-6.
21. Savelyev D. A., The investigation of the features of focusing vortex super-Gaussian beams with a variable-height diffractive axicon, Computer Optics. 45 (2) (2021) 214-221.
THE AUTHORS
SAVELYEV Dmitry A.
dmitrey.savelyev@yandex.ru
ORCID: 0000-0003-2282-3895
Received 18.07.2022. Approved after reviewing 19.07.2022. Accepted 20.07.2022.
© Peter the Great St. Petersburg Polytechnic University, 2022