## Abstract

We study diffraction of Bessel vortex beams with topological charges of ±1 and ±2 and a wavelength of 130 *µ*m on two-dimensional amplitude periodic gratings. Results of simulations and experiments at the Novosibirsk Free Electron Laser facility show that there appear periodic patterns in the planes corresponding to the classical main and fractional Talbot planes, but instead of self-images of the holes, there are observed periodic lattices of annular vortex microbeams with topological charges corresponding to the charge of the beam illuminating the grating. The ring diameters are the same for beams with different topological charges, but they are proportional to the grating period and inversely proportional to the diameter of the beam illuminating the grating.

© 2018 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

## 1. Introduction

The classical Talbot effect [1–3] is exact self-reproduction of the wavefront behind a periodical grating illuminated by a plane wave in a sequence of so called “main Talbot planes” situated at distances

where*L*= 2

_{T}*p*

^{2}/

*λ*is the Talbot length,

*p*is the grating period, and

*λ*is the wavelength. The diffraction pattern also has a remarkable property of the existence of fractional Talbot planes [4,5] located at distances described by following expression: where

*n*and

*m*are natural numbers with no common factor and

*n < m*. With

*n*/

*m*= 1/2 (“half-Talbot plane”) the image of the grating is also exactly reproduced but shifted along both coordinates by half the period. In planes with

*n*/

*m*= 1/4 and 3/4 (“quarter-Talbot planes”) there are observed grating images of double spatial frequency, whereas at

*n*/

*m*= 1/6, 1/3, 2/3 and 5/6 the patterns have a triple spatial frequency. These remarkable features of the Talbot effect are widely used in image recording devices, interferometers, illuminating systems, metrology systems,

*etc.*, and this effect has kept attracting interest in recent years. As an example we can mention paper [6], in which the self-imaging effect was applied to ablation of substrates with the use of pulsed Nd:YAG laser and a microlens array. To protect the lenses and to increase the spatial frequency of the ablation spots, the substrates were placed in one of the fractional Talbot planes behind the intermediate hole array illuminated by the microlenses. Almost all studies of the Talbot effect have been performed so far in the visible wavelength range. The classical Talbot effect in the terahertz range was first demonstrated in [7], where it was applied to the measurement of the gap between a germanium window and an array of sensors inside a 2D microbolometer camera [8], as well as to the measurement of the wavelength of free electron laser radiation.

Another area of optics still of interest is beams with an orbital angular momentum, or vortex beams [9–11]. However to date, only a few papers have been published in which the Talbot geometry was applied either to study of vortex beams or to formation of vortex-beam lattices. In [12] the two-channel moiré deflectometry was used to measure both the wavefront and the transverse component of the Poynting vector of an optical vortex beam. In [13] a vortex beam with a topological charge of 1 to 3 illuminated a linear amplitude grating, behind which diffraction patterns were observed in the distance *z*′ = 5*L _{T}*/2. The only purpose of the experiment was to demonstrate that distortion of grating slit images reflects the magnitude and sign of the topological charge of the beam. Although this task was achieved, it could be achieved in a more straightforward way, e.g. by diffraction of a vortex beam on a double slit, as it was demonstrated in [14,15]. In recent paper [16], diffraction of vortex beams from amplitude radial grating was investigated.

A greater number of studies have been devoted to the formation of lattices of vortex beams by embedding elements to form vortex beams directly into the openings of the gratings. In work [17], a periodic structure of fork-like holograms was created using a spatial light modulator (SLM). A plane wave passing through the holograms formed square or hexagonal arrays of vortex beams. Similar interferometric techniques using lateral shearing of interfering plane or spherical waves were applied in [18] for generation of a regular network of fork fringes. A lattice of microscopic spiral phase plates fabricated by photopolymerization was used for creation of vortex beam arrays in [19]. In [20], the incident plane wave illuminating a two-dimensional lattice of holographic phase elements, each forming an individual vortex beam. It was shown that in this case a lattice of vortex microbeams appears in the main and fractional Talbot planes. The advantage of this scheme is the possibility of using an SLM to form in each cell a beam with any given topological charge. In [21] two artificially formed quasi-Bessel vortex beams with different transverse wave numbers and individual topological charges interfered in free space, demonstrating that the positions of the Talbot planes do not depend on the magnitude of the topological charges. A similar geometry with annular arrangement of elements to generate vortex beams was used in [22]. For completeness, it is worth mentioning a number of publications [23–28] directly or indirectly related to the problem we are considering.

In this paper we first study the diffraction pattern that appears behind two-dimensional amplitude periodic gratings of round holes illuminated by a wide vortex beam. Since a vortex beam has more independent parameters than a plane wave does, we intended to find out how the diffraction pattern depends on the beam and grating parameters and to what extent the results agree with the classical Talbot effect. We did not expect obtaining exact self-imaging of the grating, which is a distinctive feature of the Talbot effect, but assumed that we would see some similarities with the properties of the classical Talbot effect. Almost all the experiments were performed using the high power radiation of the Novosibirsk free-electron laser (NovoFEL) [29] at a wavelength of 130 *µ*m, and one experiment was carried out at a wavelength of 49 *µ*m. A large wavelength makes it easy to manufacture the amplitude gratings used in the experiments. An additional advantage in the use of terahertz radiation is the opportunity of reaching maximum possible spatial resolution because the size of the sensitive elements of the microbolometer array [30,31] is equal to 51 *µ*m.

The plan for this article is as follows. In Section 2 we describe the experimental setup. Results of simulations are presented in Section 3. In Section 4 we compare the results of numerical simulations with the experimental results. Finally, all the results and their possible application are discussed in Section 5.

## 2. Experimental configuration

The experimental schematic is shown in Fig. 1. The Gaussian beam of the NovoFEL with an average power of about 100 W was transformed by a silicon phase binary axicon with a diameter of 30 mm and a helical structure of the zones into a vortex Bessel beam with a topological charge of ±1 or ±2

*ρ*and

*θ*are the polar coordinates and

*κ*= 2 mm

^{−1}is the radial wave number. The properties of the axicon used and the characteristics of the beams formed are described in detail in [32]. Application of these beams for formation of surface plasmon polaritons is described in [33]. The diameters of the maxima of the first rings of the beams were 1.9 and 3.0 mm. The beams were expanded by a telescopic system consisting of a silicon lens and a parabolic mirror, the magnification of which was

*M*=

*f*

_{2}/

*f*

_{1}(in the experiments

*M*= 5 or 3.3). The distance from the mirror to the grating was 600 mm. The power of the Bessel beam incident on the grating was approximately 10 W due to large Fresnel losses in the silicon optical elements and radiation loss in the positive diffraction orders on the axicon. The gratings were copper or aluminum plates with a thickness of 200 to 800

*µ*m, in which round holes of diameter

*D*= 1 or 2 mm with a period

*P*= 2 to 6 mm were drilled. The diameters of the gratings were much larger than the diameter of the axicon and the parabolic mirror. Further on, the gratings will be denoted by the abbreviation DdPp, where d and p are the numerical values of the diameter of the holes and the grating period in mm.

The diffraction pattern was recorded by a 320 × 240 microbolometer terahertz camera [30,31], which was moved by a linear translator stage along the optical axis, with a frame rate of 30 frames per second. The length of recorded video films was 750 frames. The physical size of the images was 16.32 × 12.24 mm^{2}. The simulations were performed for an area of 30 × 30 mm^{2}, but the calculated patterns are presented in the figures in the same size as the experimental data. On the one hand, such a representation makes it easy to compare numerical calculations with experimental results. On the another hand, due to geometric limitations, a number of periods on a 30-mm grating in the terahertz range cannot be more than 5 – 15. As is known (see, for example, [34]), in this case the Talbot images are not distorted only in the central part of the pattern.

## 3. Results of simulations

The calculations of the electromagnetic field behind the grating were performed by Matlab R2016b using the Rayleigh-Sommerfeld integral, the applicability of which for our geometry follows from Fig. 5 in [35]. The diffraction patterns are presented as they are seen when an observer looks along the z axis. In this case, when the beam has a positive topological charge, its phase increases clockwise. Figures 2 and 3 show the calculated diffraction patterns of a plane wave and a vortex wave with a topological charge *l* = +1 behind the gratings with openings of 1 mm in diameter and a period increasing from 2 to 6 mm for the Talbot planes with *Z* = *z*/*L _{T}* = 1/6 to

*Z*= 1.

The diffraction patterns for a plane wave (Fig. 2) completely correspond to the expectations and indicate that our simulation methods are correct. First, in the *Z* = 1 and 1/2 planes, the images of the grating are precisely reproduced (the blurring of the circle boundaries at large P is explained by the decrease in the number of interfering waves). Second, in the *Z* = 1/4 and 1/6 planes, the images of the double and tripple spatial frequencies are clearly distinguishable (they disappear when the spatial orders overlap). Third, the calculated phase distributions correspond to the expected distributions. For example, in the groups of four points in the 1/4 and 3/4 planes the phase difference, in exact agreement with the theory [4,5], changes clockwise as follows: 0, *π*/2, *π*, and *π*/2. For a given wavelength, the diffraction patterns in the case of the classical Talbot effect depend only on the grating period, and the effect of self-imaging is perfectly demonstrated.

In the case of vortex beams (Fig. 3), the patterns obey more complex regularities. The first conclusion that we can make is that the illumination of a grating by a vortex beam produces some ordered pictures with Talbot spatial periods. Remarkably, arrays of ring-like vortex beamlets can be formed in the “classical” Talbot distances (Fig. 3). This can be a reason to call this phenomenon “the quasi-Talbot effect”. Despite the variety of observed diffraction patterns, even the first glance at Fig. 3 suggests a discernible regularity in their appearance. This is clearly seen if we compare the columns of figures for the Talbot distances, equal to 1/4, 1/2, and 1, in which arrays of annular vortex beams appear (hereinafter, we normalize the distances by the Talbot length *Z* = *z*/*L _{T}* ). Summarizing the results presented in Fig. 3, as well as similar compilations (which we do not show here) of the diffraction of vortex beams with different parameters, we find that five parameters should be considered in order to describe this regularity: diameter of grating openings and distance from the grating to the observation plane, grating period, size of the ring-like illuminating beam, and its topological charge.

The diameter of the holes determines the angle of diffraction of radiation behind the grating. At *D* = 1 mm and *λ* = 130 *µ*m, the diffraction angle on a single hole is equal roughly to 0.13. For a beam of ∼ 13 mm in diameter, waves diverging from the left and right edges of the grating intersect in a distance of approximately *z*_{0} = 100 mm. Thus, the first two parameters determine the “zone of limited diffraction” (see Fig. 4) within which waves with opposite phases do not interfere. This distance corresponds to the frames located in the upper left corner of Fig. 3. In these frames, the “local” classic Talbot effect is clearly visible. This is irrefutably confirmed, for example, by the phase distribution in the above-mentioned groups of four images in the *Z* = 1/4 plane for gratings D1P2 and D1P3. The helicoidal nature of the illuminating beam, however, manifests itself in these frames as a general phase dependence exp (+*ilθ*). In distances greater than 100 mm, a periodic lattice of perfect vortex microbeams with an annular intensity distribution resembling the Fourier transform of the Bessel beam arises in the *Z* = 1/4, *Z* = 1/2, and *Z* = 1 planes. The topological charge of each microbeam corresponds to the topological charge of the illuminating beam, but the phase distribution Φ = exp [*i*(*θ* + ∆*θ _{sm}*)] is shifted depending on the microbeam position (

*x*). The diameter of the rings in the corresponding Talbot planes increases with the grating period. For the same grating, the diameter grows up with the number of the plane. When the rings begin to overlap, complex distributions of intensity and phase arise in place of the system of rings. This corresponds to the upper region in Fig. 4.

_{s}, y_{m}Dependences on the topological charge, the diameter of the holes, and the diameter of the illuminating beam are clearly visible from the collections of diffraction patterns shown in Fig. 5. The diameter of the rings [Fig. 5(a)] remains constant when the sign and the magnitude of the topological charge change, which may be useful for applications. Increase in the diameter of the holes with the period unchanged does not alter the diameter of the rings [Fig. 5(b)], but the pattern is distorted when the rings begin to overlap. Another dependence [Fig. 5(c)], which should be specially noted, is the decrease in the diameter of the rings with increase in the illuminating beam diameter (in the experiment this is increase in the telescope magnification *M*). The dependencies shown in Fig. 5 can be represented graphically. The grating period is uniquely related to the Talbot length via relation (1), so we can present the dependence of the diameter of the rings as a function of the Talbot length [see Fig. 6(a)]. With good precision, this dependence turned out to be a linear one. On the contrary, the dependence of the diameter of the rings of the beamlets on the diameter of the beam illuminating the grating [Fig. 6(b)] is inversely proportional to the size of the latter.

In addition, we give some more numerical results. We investigated the diffraction of vortex Laguerre-Gaussian beams on the same gratings. In this case, lattices of annular microbeams are also observed in the Talbot planes. When gratings with asymmetrical holes are illuminated by Bessel beams, oval-shaped rings are observed in the Talbot planes.

## 4. Comparison with experiments

A set of frames extracted from terahertz videos recorded by the microbolometer array is presented in Fig. 7 for distances expressed in the Talbot length. Comparison of the frames shows that the patterns recorded in the experiments are in good agreement with the results of the numerical calculations. This is clearly seen from comparison of rows **e–f** and **g–h**. Comparing series **a–b** and **c–d**, we see that the experimental frames corresponding to the calculated ones are situated a little bit farther from the grating. This can be easily explained if we assume that in the experiments with *M* = 5 a slightly divergent wave was formed after the telescope, as a result of which the position of the Talbot planes shifted.

Comparing the results of the experiments and calculations, we see that in both cases bright spots are observed in the center of the rings formed in the Talbot *L* = 1/2 plane, and in the experimentally recorded frames they are somewhat brighter. The emergence of this spots are easily explained by the fact that, in addition to the main Bessel beam, a weaker beam with a plane wave front is formed in our experiments behind the binary silicon axicon. The appearance of the spots is caused by two reasons. The binary axicon used by us was calculated for a wavelength of 141 *µ*m, whereas the experiments were carried out at a wavelength of 130 *µ*m, and therefore a small part of the radiation propagated as a zero-order wave of diffraction. This portion of the radiation was also taken into account in the calculations. However, one more source of the plane wave was not taken into account in the calculations: the Fresnel reflections inside the axicon. In our case, the interference of these waves turned out to be constructive (for details, see Sec. III(b) in [32]), which additionally increased the intensity of the undifferentiated wave in the experiment. Thus, these waves are responsible for the appearance of the classical Talbot pattern in the *Z* = 1/2 plane.

Fig. 7(b,d) shows that in the experiment the rings are transformed in spirals. It is known that such patterns appear when the vortex beam interferes with a spherical wave. Our numerical calculations of Talbot carpets, which we do not present here, have shown that, because of the finite dimensions of the incident beam, the electromagnetic wave behind the grating is divergent. The interference of the vortex waves with the spherical wavefront of the non-diffracted portion of the beam forms helical (a) and double-helical (b) patterns. In the figure they are not seen well enough due to poor spatial resolution, but a similar diffraction pattern with an excellent resolution is shown in Fig. 26 in a review [11].

In addition to the experiments described above, we carried out an experiment on the diffraction of vortex beams on gratings at the wavelength of the NovoFEL equal to 49 *µ*m. At this wavelength, we also observed the formation of annular beams.

## 5. Discussion

The Talbot effect is revisited for the case of vortex Bessel beams illuminating two-dimensional gratings. Regular patterns with spatial frequencies corresponding to the main and fractional Talbot planes are observed in respective distances. The diffraction patterns are more complex than in the case of the classical Talbot effect and depend on a larger number of parameters. Instead of images of the holes, there are ring-shaped microbeam arrays formed in the Talbot planes within a certain distance interval. For each of these rings, as it was shown in the calculations, the topological charge reproduces that of the illuminating beam. The greater is the grating period, the larger are the ring diameters. Expansion of the illuminating beam leads to decrease in the ring diameters. Keeping the diameter of the illuminating beam unchanged and increasing the lattice period, we find that the diameter of the rings increases linearly. “Images” of holes in the Talbot planes are ideal single rings. Recall that the spatial Fourier spectrum of Bessel beams is a ring and that the diameter of this ring for beams created by a binary axicon does not depend on the topological charge, and if the beam is expanded by a telescopic system, its angular spectrum decreases linearly. Recall also that the Talbot effect can be described by means of a fractional Fourier transform [28]. Taking into account what was said above, one can conclude from the dependences obtained in this paper that the properties of the “ideal” vortex beams observed in the Talbot planes (the beams whose diameter does not depend on the topological charge) are determined by the Fourier-transform properties of the Bessel beams.” More detailed analytical calculations will be carried out later.

The phenomenon described in the paper can be used in various applications, for example, to create an array of optical traps, in micromachining, and in microfluidic devices. This method has a useful property of the production in the planes with *Z* = 1/4 of arrays of beamlets with a double spatial frequency, which densely fill the space. The formation of a vortex beam array with a simple metal grating allows using high power laser beams. This is especially important in the terahertz optics, where metals reflect radiation practically completely, whereas dielectric optical elements have higher absorption as compared with elements in visible range, and it is difficult to create, for example, terahertz light modulators commonly used for transformation of wavefronts.

## Funding

RFBR grant (15-02-06444); Russian Science Foundation (14-50-00080); The Ministry of Education and Science of the Russian Federation (RFMEFI62117X0012, 16.7894.2017/6/7).

## Acknowledgments

The authors are grateful to G. N. Kulipanov, V. G. Serbo and V. A. Soifer for stimulating discussions, B. O. Volodkin for technical support, and Ya. V.Getmanov, V. V. Kubarev, T. V.Salikova, M. A. Scheglov, O. A.Shevchenko, D. A. Skorokhod, and other members of the NovoFEL team for the invaluable support of the experiments.

## References and links

**1. **H.F. Talbot, “Facts relating to optical science,” Philos. Mag. Ser. **3**, 401–407 (1836).

**2. **J. T. Winthrop and C. R. Worthihgton, “Theory of Fresnel images. I. Plane periodic objects in monochromatic light,” J. Opt. Soc. Am. **55**, 373–381 (1965). [CrossRef]

**3. **J. M. Wen, Y. Zhang, and M. Xiao, “The Talbot effect recent advances in classical optics, nonlinear optics, and quantum optics,” Adv. Opt. Photon. **5**, 83–130 (2013). [CrossRef]

**4. **J. Westerholm, J. Turunen, and J. Huttunen, “Fresnel diffraction in fractional Talbot planes - a new formulation,” J. Opt. Soc. Am. A **11**, 1283–1290 (1994). [CrossRef]

**5. **L. Martínez-León, M. Araiza-E, B. Javidi, P. Andrés, V. Climent, J. Lancis, and E. Tajahuerce, “Single-shot digital holography by use of the fractional Talbot effect,” Opt. Express **17**, 12900–12909 (2009). [CrossRef] [PubMed]

**6. **M. Aymerich, D. Nieto, E. Álvarez, and M. T. Flores-Arias, “Laser surface microstructuring of biocompatible materials using a microlens array and the Talbot effect: evaluation of the cell adhesion,” Materials **10**, 10020214 (2017). [CrossRef]

**7. **B. A. Knyazev, V. S. Cherkassky, Y. Y. Choporova, V. V. Gerasimov, and M. G. Vlasenko, “The Talbot effect in the terahertz spectral range,” in Proceedings of the 35th International Conference on Infrared Millimeter and Terahertz Waves (IRMMW-THz), Rome, Italy, 5–10 Sept 2010, 5612558.

**8. **M. A. Demyanenko, D. G. Esaev, B. A. Knyazev, G. N. Kulipanov, and N. A. Vinokurov, “Imaging with a 90 fps microbolometer focal plane array and high-power terahertz free electron laser,” Appl. Phys. Lett. , **92**, 131116 (2008). [CrossRef]

**9. **L. Allen, M. W. C. Beijersbergen, R.J. Spreeuw, and J. P. Woerdman, “Orbital angular momentum of light and the transformation of Laguerre-Gaussian laser modes,” Phys. Rev. A **45**, 8185 (1992). [CrossRef] [PubMed]

**10. **D. L. Andrews and M. Babiker, *The angular momentum of light*. (Cambrige University Press, 2013).

**11. **B. A. Knyazev and V. G. Serbo, “Photon beams with non-zero projection of orbital momentum – new results,” Physics Uspekhi (to be published in May 2018).

**12. **M. Yeganeh, S. Rasouli, M. Dashti, S. Slussarenko, E. Santamato, and E. Karimi, “Reconstructing the Poynting vector skew angle and wavefront of optical vortex beams via two-channel moiré deflectometery,” Opt. Lett. **38**, 887–889 (2013). [CrossRef] [PubMed]

**13. **P. Panthong, S. Srisuphaphon, A. Pattanaporkratana, S. Chiangga, and S. A. Deachapunya, “Study of optical vortices with the Talbot effect,” J. Opt. **18**, 035602 (2016). [CrossRef]

**14. **H. I. Sztul and R. R. Alfano, “Double-slit interference with Laguerre-Gaussian beams,” Opt. Lett. **31**, 999–1001 (2006). [CrossRef] [PubMed]

**15. **O. Emile and J. Emile, “Young’s double-slit interference pattern from a twisted beam,” Appl. Phys. B **117**, 487–491 (2014). [CrossRef]

**16. **D. Herbi, S. Rasouli, and M. Yaganeh, “Intensity based measurement of the topological charge alteration by the diffraction of vortex beams from amplitude sinusoidal radial grating,” J. Opt. Soc. Am. B **35**, 724–730 (2018). [CrossRef]

**17. **A. Dreischuh, S. Chervenkov, D. Neshev, G. G. Paulus, and H. Walther, “Generation of lattice structures of optical vortices,” J. Opt. Soc. Am. B **19**, 550–556 (2002). [CrossRef]

**18. **D. P. Ghai, S. Vyas, P. Senthilkumaran, and R. S. Sirohi, “Vortex lattice generation using interferometric techniques based on lateral shearing,” Opt. Commun. **282**, 2692–2698 (2009). [CrossRef]

**19. **E. Brasselet, M. Malinauskas, A. Žukauskas, and S. Juodkazis, “Photopolymerized microscopic vortex beam generators: Precise delivery of optical orbital angular momentum,” Appl. Phys. Lett. **97**, 211108 (2010) [CrossRef]

**20. **Y. Lu, B. Jiang, S. Lü, Y. Liu, S. Li, Z. Cao, and X. Qi, “Arrays of Gaussian vortex, Bessel and Airy beams by computer-generated hologram,” Opt. Commun. **363**, 85–90 (2016). [CrossRef]

**21. **J. S. Rodrigues, C. V. Mendes, E. J. Fonseca, and A. J. Jesus-Silva, “Talbot effect in optical lattices with topological charge,” Opt. Lett. **42**, 3944–39472017. [CrossRef] [PubMed]

**22. **H. Gao, Y. Li, L. Chen, J. Jin, M. Pu, X. Li, P. Gao, C. Wang, X. Luo, and M. Hong, “Quasi-Talbot effect of orbital angular momentum beams for generation of optical vortex arrays by multiplexing metasurface design,” Nanoscale , **10**, 666–671 (2018). [CrossRef]

**23. **J. Courtial, R. Zambrini, M. R. Dennis, and M. Vasnetsov, “Angular momentum of optical vortex arrays,” Opt.Expr. , **14**, 938–949 (2006).

**24. **A. Bezryadina, D. N. Neshev, A. S. Desyatnikov, J. Young, Z. Chen, and Y. S. Kivshar, “Observation of topological transformations of optical vortices in two-dimensional photonic lattices,” Opt. Expr. , **14**, 8317–8327 (2006). [CrossRef]

**25. **N. Gao and C. Xie, “Experimental demonstration of free-space optical vortex transmutation with polygonal lenses,” Opti. Lett. , **37**, 3255–3257 (2012). [CrossRef]

**26. **A. Ferrando and M. A. García-March, “Theory for the control of dark rays by means of discrete symmetry diffractive elements,” J. Opt. , **15**, 044014 (2013). [CrossRef]

**27. **A. Sabatyan and S. M. T. Balanoji, “Square array of optical vortices generated by multiregion spiral square zone plate,” JOSAA **33**, 1793–1797(2016). [CrossRef]

**28. **C. Schnebelin and H. G. de Chatellus, “Fractional Fourier transform-based description of the Talbot effect: application to analog signal processing,” Appl. Opt. , **56**, A62–A68 (2017). [CrossRef]

**29. **G. N. Kulipanov, E. G. Bagryanskaya, E. N. Chesnokov, Y. Y. Choporova, V. V. Gerasimov, Y. V. Getmanov, S. L. Kiselev, B. A. Knyazev, V. V. Kubarev, S. E. Peltek, V. M. Popik, T. V. Salikova, M. A. Scheglov, S. S. Seredniakov, O. A. Shevchenko, A. N. Skrinsky, S. L. Veber, and N. A. Vinokurov, “Novosibirsk free electron laser: facility description and recent experiments,” IEEE Trans. THz Sci. Technol. **5**, 798–809 (2015). [CrossRef]

**30. **M. A. Demýanenko, D. G. Esaev, V. N. Ovsyuk, B. I. Fomin, A. L. Aseev, B. A. Knyazev, G. N. Kulipanov, and N. A. Vinokurov, “Microbolometer detector arrays for the infrared and terahertz ranges,” J. Opt. Technol. **76**, 739–743 (2009). [CrossRef]

**31. **B. A. Knyazev, V. S. Cherkassky, Y. Choporova Y., V. Gerasimov V., M. G. Vlasenko, M. A. Demýanenko, and D. G. Esaev, “Real-time imaging using a high-power monochromatic terahertz source: comparative description of imaging techniques with examples of application,” J. Infrared Millim. THz Waves **32**, 1207–1222 (2011). [CrossRef]

**32. **Yu. Yu. Choporova, B. A. Knyazev, G. N. Kulipanov, V. S. Pavelyev, M. A. Scheglov, N. A. Vinokurov, B. O. Volodkin, and V. N. Zhabin, “High-power Bessel beams with orbital angular momentum in the terahertz range,” Phys. Rev. A **96**, 023846 (2017). [CrossRef]

**33. **B. A. Knyazev, Yu. Yu. Choporova, M. S. Mitkov, V. S. Pavelyev, and B. O. Volodkin, “Generation of terahertz surface plasmon polaritons using nondiffractive Bessel beams with orbital angular momentum,” Phys. Rev. Lett. **115**, 163901 (2015). [CrossRef] [PubMed]

**34. **V. Arrizón and G. Rojo-Velázquez, “Fractional Talbot field of finite gratings: compact analytical formulation,” J. Opt. Soc. Am. A **18**, 1252–1256 (2001). [CrossRef]

**35. **Yu. Yu. Choporova, B. A. Knyazev, and M. S. Mitkov, “Classical holography in the terahertz range: recording and reconstruction techniques,” IEEE Trans. THz Sci. Technol. **5**, 836–844 (2015). [CrossRef]