<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.0 20120330//EN" "JATS-archivearticle1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink">
  <front>
    <journal-meta />
    <article-meta>
      <article-id pub-id-type="doi">10.18287/1613-0073-2015-1490-97-104</article-id>
      <title-group>
        <article-title>Simulation of linear gradient lenses for subwavelength focusing of Gaussian beams</article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <string-name>Savelyev D.A.</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <aff id="aff0">
          <label>0</label>
          <institution>Samara State Aerospace University, Image Processing Systems Institute, Russian Academy of Sciences</institution>
        </aff>
      </contrib-group>
      <pub-date>
        <year>2015</year>
      </pub-date>
      <fpage>97</fpage>
      <lpage>104</lpage>
      <abstract>
        <p>This paper demonstrates the use of linear gradient lenses (a diverging lens and a converging lens) for subwavelength focusing of Gaussian beams. Numerical simulations have shown that increasing the length lens provides a more compact focal region in both the transverse and longitudinal directions.</p>
      </abstract>
      <kwd-group>
        <kwd>subwavelength focusing</kwd>
        <kwd>linear gradient lense</kwd>
        <kwd>FDTD</kwd>
        <kwd>Gaussian beams</kwd>
      </kwd-group>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>Introduction</title>
      <p>
        The subject of the study gradient optics is phenomena associated with optical
effects in media with a gradual change of the refractive index [
        <xref ref-type="bibr" rid="ref1">1</xref>
        ]. Medium consisting
of materials which are of the type GRIN (GRadient INdex), have distributed on the
gradient refractive index material [
        <xref ref-type="bibr" rid="ref2">2</xref>
        ]. Thus, light rays are bent in a curve passing
through them, which allows using them for the best focus [
        <xref ref-type="bibr" rid="ref3 ref4 ref5 ref6">3-6</xref>
        ]. In [
        <xref ref-type="bibr" rid="ref3">3</xref>
        ] used
hyperbolic planar lens, which allowed to focus the light in the focal spot with a
diameter of the half-width at half intensity (FWHM) equal 0,131λ, in [
        <xref ref-type="bibr" rid="ref4">4</xref>
        ] dealt with a
2-D gradient microlenses Mikaelian, in [
        <xref ref-type="bibr" rid="ref5 ref6">5-6</xref>
        ] - two-dimensional photonic crystals.
Photonic crystals due to their properties are used for several applications [
        <xref ref-type="bibr" rid="ref7">7</xref>
        ], among
which their application to optical waveguides [
        <xref ref-type="bibr" rid="ref8 ref9">8-9</xref>
        ].
      </p>
      <p>
        Currently in the integrated and fiber optics there is a large variety of optical
waveguides with different properties. Planar waveguides and fibers along the profile
of the spatial distribution of the refractive index can be divided into two groups: with
step refractive index profile and with gradient profile [
        <xref ref-type="bibr" rid="ref10 ref11">10-11</xref>
        ]. In the case of a gradient
profile refractive index varies smoothly from the center of the waveguide to the
boundary defined by the law [
        <xref ref-type="bibr" rid="ref11">11</xref>
        ].
      </p>
      <p>One of the major advantages of the lenses belonging to the gradient optics is that
the optical surfaces of gradient lenses may be flat, which is important for the
inputoutput of radiation of the fiber.</p>
      <p>
        It is known that high-order mode exhibit greater divergence, however, the use of
multimode method is one of the possible ways to increase the capacity of modern
communication systems [
        <xref ref-type="bibr" rid="ref12">12</xref>
        ]. And the most promising is the multiplexing of modes
with different orders of the vortex phase singularity [
        <xref ref-type="bibr" rid="ref13 ref14">13-14</xref>
        ].
      </p>
      <p>
        Distribution and generation of free-space laser vortices studied quite well [
        <xref ref-type="bibr" rid="ref15 ref16 ref17 ref18">15-18</xref>
        ],
in [
        <xref ref-type="bibr" rid="ref19 ref20 ref21">19-21</xref>
        ] shows how to use them to focus the simple micro-elements of square and
cylindrical shape. However, obtaining an optical fiber separate vortex modes and their
superpositions is a big problem [
        <xref ref-type="bibr" rid="ref14 ref22 ref23 ref24 ref25 ref26">14, 22-26</xref>
        ]. In contrast to the classic LP-mode (mode
optical fibers), angular harmonics are invariant to the scale at the input and output of
fiber with the help of diffraction microstructures. This gives greater freedom in
choosing the parameters of the optical circuit. In [
        <xref ref-type="bibr" rid="ref27">27</xref>
        ] have been identified application
features sharp focusing of laser modes for introduction into a fiber of smaller diameter
with using the binary micro-relief applied at the output end of the optical fiber.
Overlaying certain conditions can generate laser light having a property of
reproducible [
        <xref ref-type="bibr" rid="ref28">28</xref>
        ], while the phase shift between the modes at certain distances to
approximate the desired cross-sectional intensity distribution of the laser beam [
        <xref ref-type="bibr" rid="ref29">29</xref>
        ].
      </p>
      <p>
        In this paper, numerically investigate the diffraction of Gaussian beams on a linear
gradient lens with different length. For the numerical simulation of diffraction of the
considered laser beams using finite-difference time-domain method (FDTD),
implemented in the software package Meep [
        <xref ref-type="bibr" rid="ref30">30</xref>
        ]. Numerical simulations were
performed using the computing power of a cluster 775 GFlops. Features a cluster: the
number of cores – 116, computes nodes: 7 twin servers HP ProLiant 2xBL220c, the
amount of RAM 112 GB.
      </p>
    </sec>
    <sec id="sec-2">
      <title>Investigation of diffraction of Gaussian beams</title>
      <p>In the paper considers two types of laser beams, which retain their structure during
propagation in the free space in the circular polarization: a fundamental Gaussian
mode, mode Gauss-Laguerre (0,1).</p>
      <p>Figure 1 shows the transverse structure of the considered optical trace element
(diverging lens and the converging lens). Later in the linear gradient lenses, in which
the refractive index decreases from the central part to the edges of the lens will be
called converging lens. The refractive index gradient lenses varies linearly as: all
considered 10 rings (width ), each of which had a their refractive index (from n =
3.47 on the center, until n = 1.5 on the edge of the lens).
0.8
0.6
0.4
0.2
0
0.5
1.0
1.5
2.0
2.5
3.0
3.5
4.0
4.5
5.0</p>
      <p>Simulation parameters: the wavelength  = 1.55 microns, the size of the
computational domain x, y,z[–4.5; 4.5]. The thickness of the absorbing layer
PML ~ 0.65 (1 micron), the sampling step of space – /31, the sampling step of time
– /(62c), where c is the velocity of light. Simulation time - 40 tacts (one for that
meant the propagation time over a single wavelength), and was chosen so that the
kind of laser radiation did not change with increasing simulation time (Figure 2).</p>
      <p>The size of the focal spot on the FWHM considered a global maximum (maximum
intensity). Table 1 shows the diffraction of a fundamental Gaussian mode and mode
Gauss-Laguerre (0,1) on the converging lens when you change the length L.
Considered the overall intensity.</p>
      <p>The table shows both the diffraction pattern and the size of the focal spot, that
with increasing length of the focusing lens has been increasing. If L = 0.25λ presence
lens has almost no effect, then L = 2λ diffraction pattern has changed significantly:
the size of the focal spot on the half-width at half intensity decreased by 2.67 times
(from FWHMmax = 1.95λ for L = 0.25λ, to FWHMmax = 0.73λ for L = 2λ).
Comparing 4 and 5 columns of Table 1, it can be seen that the depth of focus,
decreased significantly with a slight reduction of the transverse dimension of the focal
spot in both these cases. Increasing the length of the lens results in better focus, it is
clearly seen when comparing the first and last column of Table 1 for the case of
fashion Gauss - Laguerre (0.1). Note that FWHMmax also decreases with increasing
length of the lens. Also the global maximum (max) also considered at the edge section
of the optical element (out).
FWHMout = 3.92λ FWHMout=0.48λ
FWHMmax = 1.41λ FWHMmax = 0.76λ</p>
      <p>Table 2 shows research of changes in the length L of the converging linear
gradient lenses separately dedicated to the longitudinal component of the electric
field. In addition to the global maximum (max) is considered the first maximum of the
optical element (fmx).
l
ta e
en od
m m
a n
nd ia</p>
      <p>s
fu s</p>
      <p>u
e a
hT G
)
1
,
0
(
e
r
r
e
u
g
a
L
s
s
u
a
G
e
d
o
M
FWHMfmx = 1.22λ FWHMfmx = 0.33λ FWHMfmx = 0.36λ FWHMout = 0.44λ FWHMout = 0.4λ
FWHMmax = 1.16λ FWHMmax = 0.9λ FWHMmax = 0.77λ FWHMmax = 0.58λ FWHMmax = 0.46λ</p>
      <p>The longitudinal electric field component (component z) in the case of the
fundamental Gaussian mode also becomes smaller with increasing length of the lens.
Comparing Tables 1 and 2 for the case of mode Laguerre - (0,1), it should be noted
that increasing the length of the lens, there is a significant enhancement z components,
and the size of the focal spot begins to form mainly due to the longitudinal component
of the electric field. This is particularly noticeable when analyzing the last column of
table 1 and 2: for mode Gauss - Laguerre (0,1) size of the focal spot next to the central
part of which contains the longitudinal component of the electric field is FWHMout =
0.4λ, while the total intensity of the light Spot contains a cross-polarized side lobes
that broadens the spot size to FWHM = 0.48λ. In the plane of maximum intensity, the
effect of cross-polarized components of the electric field is more significant
(broadened to 0.76λ). Of special note is the case of L = 0.5λ, when the longitudinal
component of the electric field near the plane of the first element in the maximum was
obtained by the minimum size of the focal spot: FWHMfm = 0.33λ, but as you can see
on the the graphs - by increasing the intensity of the side lobes.</p>
      <p>Table 3 shows the result of diffraction of Gaussian beams at the diverging linear
gradient lens when changing the length L, that is, on the lens, whose refractive index
increases from the center to the edges. Such a lens can be used for the introduction of
laser radiation into a hollow (annular) optical waveguide. The refractive index in this
case, i have a similar case considered earlier: also examined 10 rings (width ), each
0.5
0.0
of which had a their refractive index. But they are changed in the reverse order of n =
1.5 in the middle, to n = 3.47 at the edge of the lens (Figure 1b).</p>
      <p>Table 3 shows that the longitudinal component of the electric field using
GaussLaguerre modes (0.1) in the plane of maximum intensity able to obtain a focal spot
size smaller than the collecting lens, and in the case of L = 1.55λ size of the focal spot
is the focal spot in the case of L = 2λ in Table 2. However, the intensity of the spots
materially ~ 5 times less.
0.8
0.6</p>
      <p>FWHMmax = 0.8λ
FWHMmax = 0.63λ</p>
      <p>FWHMmax = 0.46λ</p>
      <p>Figure 3 is a graph showing the cross section of the components of z to
LaguerreGaussian modes (0,1) (bottom row of Table 3). As seen from the graph, a reduction
intensity of the focal spot accompanied by decrease intensity of side lobes.</p>
    </sec>
    <sec id="sec-3">
      <title>Conclusion</title>
      <p>In this paper using the FDTD method, was numerically research of circularly
polarized Gaussian beams passing through a linear gradient lens, considered two
options of lenses - converging and diverging.</p>
      <p>It has been shown that increasing the length lens provides a more compact size of
the focal region in all directions, both for total intensity and the longitudinal electric
field component.</p>
      <p>The minimum size of the focal spot was not obtained at the maximum intensity,
but at the first maximum out of optical element. For modes Gauss-Laguerre (0,1) the
minimum size of the focal spot on the FWHM for the longitudinal electric field
component was 0.33 wavelength for L = 0.5λ; in the plane of maximum intensity
using a diverging lens for z components able to obtain a focal spot size less than when
using a converging lens (at falling intensity ~ 5 times).</p>
    </sec>
    <sec id="sec-4">
      <title>Acknowledgements</title>
      <p>The work was partially funded by the Russian Foundation of Basic Research
Grant (№ 14-07-31079 mol_a) and Russian Federation Ministry of Education and
Science.</p>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          1.
          <string-name>
            <surname>Tomlinson</surname>
            <given-names>WJ</given-names>
          </string-name>
          .
          <article-title>Applications of GRIN-rod lenses in optical fiber communication systems</article-title>
          .
          <source>Applied Optics</source>
          ,
          <year>1980</year>
          ;
          <volume>19</volume>
          (
          <issue>7</issue>
          ):
          <fpage>1127</fpage>
          -
          <lpage>1138</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          2.
          <string-name>
            <surname>Emkey</surname>
            <given-names>W</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Jack</surname>
            <given-names>C</given-names>
          </string-name>
          .
          <article-title>Analysis and evaluation of graded-index fiber lenses</article-title>
          .
          <source>Journal of Lightwave Technology</source>
          ,
          <year>1987</year>
          ;
          <volume>5</volume>
          (
          <issue>9</issue>
          ):
          <fpage>1156</fpage>
          -
          <lpage>1164</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          3.
          <string-name>
            <surname>Kotlyar</surname>
            <given-names>VV</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Kovalev</surname>
            <given-names>АА</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Triandafilov</surname>
            <given-names>YaR</given-names>
          </string-name>
          , Nalimov AG.
          <article-title>Modes of planar gradientindex hyperbolic secant waveguide</article-title>
          .
          <source>Computer optics</source>
          ,
          <year>2010</year>
          ;
          <volume>34</volume>
          (
          <issue>2</issue>
          ):
          <fpage>146</fpage>
          -
          <lpage>154</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          4.
          <string-name>
            <surname>Kotlyar</surname>
            <given-names>VV</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Kovalev</surname>
            <given-names>АА</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Nalimov</surname>
            <given-names>AG</given-names>
          </string-name>
          .
          <article-title>Gradient-index element of micro-optics for super-resolution</article-title>
          .
          <source>Computer Optics</source>
          ,
          <year>2009</year>
          ;
          <volume>33</volume>
          (
          <issue>4</issue>
          ):
          <fpage>369</fpage>
          -
          <lpage>378</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          5.
          <string-name>
            <surname>Matsumoto</surname>
            <given-names>T.</given-names>
          </string-name>
          <article-title>Focusing of light by negative refraction in a photonic crystal slab superlens on silicon-on-insulator substrate</article-title>
          .
          <source>Optics letters</source>
          ,
          <year>2006</year>
          ;
          <volume>31</volume>
          (
          <issue>18</issue>
          ):
          <fpage>2786</fpage>
          -
          <lpage>2788</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref6">
        <mixed-citation>
          6.
          <string-name>
            <surname>Kurt</surname>
            <given-names>H.</given-names>
          </string-name>
          <article-title>The focusing effect of graded index photonic crystals</article-title>
          .
          <source>Applied Physics Letters</source>
          ,
          <year>2008</year>
          ;
          <volume>93</volume>
          (
          <issue>17</issue>
          ):
          <fpage>171108</fpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref7">
        <mixed-citation>
          7.
          <string-name>
            <surname>Zhu</surname>
            <given-names>Q</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Jin</surname>
            <given-names>L</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Fu</surname>
            <given-names>Y.</given-names>
          </string-name>
          <article-title>Graded index photonic crystals: A review</article-title>
          .
          <source>Annalen der Physik</source>
          ,
          <year>2015</year>
          ;
          <volume>527</volume>
          (
          <issue>3-4</issue>
          ):
          <fpage>205</fpage>
          -
          <lpage>218</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref8">
        <mixed-citation>
          8.
          <string-name>
            <surname>Xiong</surname>
            <given-names>C</given-names>
          </string-name>
          .
          <article-title>Slow-light enhanced correlated photon pair generation in a silicon photonic crystal waveguide</article-title>
          .
          <source>Optics letters</source>
          ,
          <year>2011</year>
          ;
          <volume>36</volume>
          (
          <issue>17</issue>
          ):
          <fpage>3413</fpage>
          -
          <lpage>3415</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref9">
        <mixed-citation>
          9.
          <string-name>
            <surname>Liang</surname>
            <given-names>J</given-names>
          </string-name>
          .
          <article-title>Wideband ultraflat slow light with large group index in a W1 photonic crystal waveguide</article-title>
          .
          <source>Journal of Applied Physics</source>
          ,
          <year>2011</year>
          ;
          <volume>110</volume>
          (
          <issue>6</issue>
          ):
          <fpage>063103</fpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref10">
        <mixed-citation>
          10.
          <string-name>
            <surname>Goncharenko</surname>
            <given-names>А</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Karpenko</surname>
            <given-names>V</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Goncharenko</surname>
            <given-names>I.</given-names>
          </string-name>
          <article-title>Fundamentals of the theory of optical waveguides</article-title>
          . Minsk: Belarusian science;
          <year>2009</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref11">
        <mixed-citation>
          11.
          <string-name>
            <surname>Snyder</surname>
            <given-names>AW</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Love</surname>
            <given-names>J</given-names>
          </string-name>
          .
          <source>Optical Waveguide Theory. Springer Science &amp; Business Media</source>
          ,
          <year>1983</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref12">
        <mixed-citation>
          12.
          <string-name>
            <surname>Bozinovic</surname>
            <given-names>N.</given-names>
          </string-name>
          <article-title>Terabit-scale orbital angular momentum mode division multiplexing in fibers</article-title>
          .
          <source>Science</source>
          ,
          <year>2013</year>
          ;
          <volume>340</volume>
          (
          <issue>6140</issue>
          ):
          <fpage>1545</fpage>
          -
          <lpage>1548</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref13">
        <mixed-citation>
          13.
          <string-name>
            <surname>Essiambre</surname>
            <given-names>R</given-names>
          </string-name>
          .
          <article-title>Capacity limits of optical fiber networks</article-title>
          .
          <source>Journal of Lightwave Technology</source>
          ,
          <year>2010</year>
          ;
          <volume>28</volume>
          (
          <issue>4</issue>
          ):
          <fpage>662</fpage>
          -
          <lpage>701</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref14">
        <mixed-citation>
          14.
          <string-name>
            <surname>Khonina</surname>
            <given-names>SN</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Kazanskiy</surname>
            <given-names>NL</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Soifer</surname>
            <given-names>VA</given-names>
          </string-name>
          .
          <article-title>Optical vortices in a fiber: Mode division multiplexing and multimode self-imaging</article-title>
          .
          <source>Editer Dr Moh. Yasin. Recent Progress in Optical Fiber Research. InTech</source>
          ,
          <year>2012</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref15">
        <mixed-citation>
          15.
          <string-name>
            <surname>Savelyev</surname>
            <given-names>DA</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Khonina</surname>
            <given-names>SN</given-names>
          </string-name>
          .
          <article-title>Maximizing the longitudinal electric component at diffraction on a binary axicon linearly polarized radiation</article-title>
          .
          <source>Computer Optics</source>
          ,
          <year>2012</year>
          ;
          <volume>36</volume>
          (
          <issue>4</issue>
          ):
          <fpage>511</fpage>
          -
          <lpage>517</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref16">
        <mixed-citation>
          16.
          <string-name>
            <surname>Khonina</surname>
            <given-names>SN</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Savelyev</surname>
            <given-names>DA</given-names>
          </string-name>
          .
          <article-title>High-aperture binary axicons for the formation of the longitudinal electric field component on the optical axis for linear and circular polarizations of the illuminating beam</article-title>
          .
          <source>Journal of Experimental and Theoretical Physics</source>
          ,
          <year>2013</year>
          ;
          <volume>117</volume>
          (
          <issue>4</issue>
          ):
          <fpage>623</fpage>
          -
          <lpage>630</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref17">
        <mixed-citation>
          17.
          <string-name>
            <surname>Khonina</surname>
            <given-names>SN</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Karpeev</surname>
            <given-names>SV</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Alferov</surname>
            <given-names>SV</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Savelyev</surname>
            <given-names>DA</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Laukkanen</surname>
            <given-names>J</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Turunen</surname>
            <given-names>J</given-names>
          </string-name>
          .
          <article-title>Experimental demonstration of the generation of the longitudinal E-field component on the optical axis with high-numerical-aperture binary axicons illuminated by linearly and circularly polarized beams</article-title>
          .
          <source>Journal of Optics</source>
          ,
          <year>2013</year>
          ;
          <volume>15</volume>
          :
          <fpage>085704</fpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref18">
        <mixed-citation>
          18.
          <string-name>
            <surname>Khonina</surname>
            <given-names>SN</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Savelyev</surname>
            <given-names>DA</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Kazanskiy</surname>
            <given-names>NL</given-names>
          </string-name>
          .
          <article-title>Vortex phase elements as detectors of polarization state</article-title>
          .
          <source>Optics Express</source>
          ,
          <year>2015</year>
          ;
          <volume>23</volume>
          (
          <issue>14</issue>
          ):
          <fpage>17845</fpage>
          -
          <lpage>17859</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref19">
        <mixed-citation>
          19.
          <string-name>
            <surname>Khonina</surname>
            <given-names>SN.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Savelyev</surname>
            <given-names>DA.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Ustinov</surname>
            <given-names>АV</given-names>
          </string-name>
          .
          <article-title>Diffraction of laser beam on a two-zone cylindrical microelement</article-title>
          .
          <source>Computer Optics</source>
          ,
          <year>2013</year>
          ;
          <volume>37</volume>
          (
          <issue>2</issue>
          ):
          <fpage>160</fpage>
          -
          <lpage>169</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref20">
        <mixed-citation>
          20.
          <string-name>
            <surname>Kotlyar</surname>
            <given-names>VV</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Stafeev</surname>
            <given-names>SS</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Feldman АYu</surname>
          </string-name>
          .
          <article-title>Photonic nanojets formed by square microsteps</article-title>
          .
          <source>Computer Optics</source>
          ,
          <year>2014</year>
          ;
          <volume>38</volume>
          (
          <issue>1</issue>
          ):
          <fpage>72</fpage>
          -
          <lpage>80</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref21">
        <mixed-citation>
          21.
          <string-name>
            <surname>Savelyev</surname>
            <given-names>DA</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Khonina</surname>
            <given-names>SN</given-names>
          </string-name>
          .
          <article-title>Numerical analysis of subwavelength focusing using a silicon cylinder</article-title>
          .
          <source>Computer Optics</source>
          ,
          <year>2014</year>
          ;
          <volume>38</volume>
          (
          <issue>4</issue>
          ):
          <fpage>638</fpage>
          -
          <lpage>642</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref22">
        <mixed-citation>
          22.
          <string-name>
            <surname>Berdague</surname>
            <given-names>S</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Facq</surname>
            <given-names>P</given-names>
          </string-name>
          .
          <article-title>Mode division multiplexing in optical fibers</article-title>
          .
          <source>Applied Optics</source>
          ,
          <year>1982</year>
          ;
          <volume>21</volume>
          :
          <fpage>1950</fpage>
          -
          <lpage>1955</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref23">
        <mixed-citation>
          23.
          <string-name>
            <surname>Mikaelian</surname>
            <given-names>AL</given-names>
          </string-name>
          .
          <article-title>Optical Methods in Informatics</article-title>
          . Moscow: Nauka (Science) Publishers,
          <year>1990</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref24">
        <mixed-citation>
          24.
          <string-name>
            <surname>Soifer</surname>
            <given-names>VA</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Golub</surname>
            <given-names>MA</given-names>
          </string-name>
          .
          <article-title>Laser beam mode selection by computer-generated holograms</article-title>
          .
          <source>Boca Raton</source>
          , CRC Press,
          <year>1994</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref25">
        <mixed-citation>
          25.
          <string-name>
            <surname>Karpeev</surname>
            <given-names>SV</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Khonina</surname>
            <given-names>SN</given-names>
          </string-name>
          .
          <article-title>Experimental excitation and detection of angular harmonics in a step-index optical fiber</article-title>
          .
          <source>Optical Memory &amp; Neural Networks (Information Optics)</source>
          ,
          <year>2007</year>
          ;
          <volume>16</volume>
          (
          <issue>4</issue>
          ):
          <fpage>295</fpage>
          -
          <lpage>300</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref26">
        <mixed-citation>
          26.
          <string-name>
            <surname>Khonina</surname>
            <given-names>SN</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Striletz</surname>
            <given-names>AS</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Kovalev</surname>
            <given-names>AA</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Kotlyar</surname>
            <given-names>VV</given-names>
          </string-name>
          .
          <article-title>Propagation of laser vortex beams in a parabolic optical fiber</article-title>
          .
          <source>Proceedings SPIE</source>
          ,
          <year>2010</year>
          ;
          <volume>7523</volume>
          :
          <fpage>75230B</fpage>
          -
          <fpage>1</fpage>
          -12.
        </mixed-citation>
      </ref>
      <ref id="ref27">
        <mixed-citation>
          27.
          <string-name>
            <surname>Savelyev</surname>
            <given-names>DA</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Khonina</surname>
            <given-names>SN</given-names>
          </string-name>
          .
          <article-title>Sharp focusing by means of binary relief at the end of the optical fiber</article-title>
          .
          <source>Proceedings SPIE</source>
          ,
          <year>2013</year>
          ;
          <volume>9156</volume>
          :
          <fpage>915609</fpage>
          -1-6.
        </mixed-citation>
      </ref>
      <ref id="ref28">
        <mixed-citation>
          28.
          <string-name>
            <surname>Khonina</surname>
            <given-names>SN</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Kotlyar</surname>
            <given-names>VV</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Soifer</surname>
            <given-names>VA</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Paakkonen</surname>
            <given-names>P</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Turunen</surname>
            <given-names>J</given-names>
          </string-name>
          .
          <article-title>Measuring the light field orbital angular momentum using DOE</article-title>
          .
          <source>Optical Memory and Neural Networks (Allerton Press)</source>
          ,
          <year>2001</year>
          ;
          <volume>10</volume>
          (
          <issue>4</issue>
          ):
          <fpage>241</fpage>
          -
          <lpage>255</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref29">
        <mixed-citation>
          29.
          <string-name>
            <surname>Almazov</surname>
            <given-names>AA</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Khonina</surname>
            <given-names>SN</given-names>
          </string-name>
          .
          <article-title>Periodic self-reproduction of multi-mode laser beams in graded-index optical fibers</article-title>
          .
          <source>Optical Memory and Neural Networks (Information Optics)</source>
          ,
          <year>2004</year>
          ;
          <volume>13</volume>
          (
          <issue>1</issue>
          ):
          <fpage>63</fpage>
          -
          <lpage>70</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref30">
        <mixed-citation>
          30.
          <string-name>
            <surname>Oskooi</surname>
            <given-names>AF</given-names>
          </string-name>
          .
          <article-title>Meep: A flexible free-software package for electromagnetic simulations by the FDTD method</article-title>
          .
          <source>Computer Physics Communications</source>
          ,
          <year>2010</year>
          ;
          <volume>181</volume>
          :
          <fpage>687</fpage>
          -
          <lpage>702</lpage>
          .
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>