Phase dislocations in hollow core waveguides

This paper discusses the basic concepts of phase dislocations and vortex formation in the electric fields of fundamental air core mode of hollow core waveguides with specific types of rotational symmetry of the core cladding boundary. Analysis of the behavior of the electric field phase in the transmission bands shows that the mechanism of light localization in the hollow core waveguides with discrete rotational symmetry of the core cladding boundary cannot be completely described by the ARROW model. For an accurate description of the phase behavior, it is necessary to account for phase jumps of the magnitude of pi when passing through the phase dislocations.


Introduction
The problem we focus on in this paper has been under investigation for ten years in the field of hollow core micro -structured optical fibers [1] and it is related to the explanation of the light localization mechanism in the negative curvature hollow core fibers (NCHCFs) [2,3,4]. These fibers have been alternatively known in literature as inhibited coupling photonic crystal fibers [5], antiresonant hollow core fibers [6], revolver hollow core fibers [7] and so on. In our opinion, this terminological disparity reflects the fact that the mechanism of light localization is not completely understood. Until now, two main concepts can be developed to explain the strong light localization in these fibers. The first concept is called inhibited coupling mechanism and it assumes a low spatial overlap and strong phase mismatch between transverse phases of the air core modes and the cladding modes of the hollow core waveguides [8]. The second one is based on the well -known ARROW (anti -resonant waveguide) mechanism of light localization [9], which assumes that the walls of the cladding elements of the waveguide can be considered as Fabry -Perrot type resonators. In this case, there are wavelengths associated with large losses in the waveguide when the phase incursion in the wall of the cladding element is equal to mπ and, accordingly, with small losses when the phase incursion for the fields of the air core mode is equal to (m+1/2)π, where m is an integer. These hollow core fibers have unique optical properties and various research groups all over the world have achieved outstanding results using these fibers [10,11].
In this paper, we consider two main reasons for the appearance of phase dislocations and vortices [12] in the fields and energy flows of the fundamental air core mode of hollow core waveguides. In particular, we demonstrate that the vortex of the transverse component of the Poynting vector of the fundamental air core mode centered at the origin is caused by leaky behavior of the fundamental air core mode. In this case, due to singularities of the axial components of electric and magnetic fields of the fundamental air core mode at which Re( ) Im( ) 0 On the other hand, hollow core waveguides, except for hollow core PCFs (photonic crystal fibers) [1], can have a complex shape of the corecladding boundary in the form of polygons (Kagome lattice hollow core fibers [14]) or have a cladding consisting of a set of capillaries (negative curvature hollow core fibers). In this case, as will be shown below, the formation of singularities and phase dislocations of the electric fields of the fundamental air core mode occur not only at the origin but also in the cladding. If the core -cladding boundary has a discrete rotational symmetry or the elements of the waveguide cladding are arranged according to a specific type of discrete rotational symmetry, the structure of the radial projection of the Poynting vector has the same type of the rotational symmetry. The phase singularities occur in the cladding elements wall along curves obeying the equation where i represents any of the coordinates   ,, x y z . When passing through the phase dislocation curve, the phase of the field undergoes a jump of the magnitude of π, and, therefore, radiation should be reflected from the walls of the cladding elements inside the transmission bands in a different way than according to the ARROW model.

Role of losses in vortex formation in hollow core waveguides
To demonstrate the physical mechanism of the formation of a central vortex of the transverse component of the Poynting vector located on the z -axis of the hollow core waveguide, we compared the behavior of the transverse component of the Poynting vector of the fundamental core mode for a step -index fiber with no losses and a hollow core waveguide that has an air core in an infinite volume of a dielectric (Fig. 1).
It is known that the axial components of the electric and magnetic fields of hybrid modes of the step -index fiber, localized in the core due to the principle of total internal reflection, can be written as [15]   vector component and, therefore, the optical vortex . They occur due to different mechanisms of light localization in the air core waveguides and in the step -index fiber (Fig. 1). For the air core waveguide, the presence of waveguide losses leads to occurrence of azimuthal and radial projections of the transverse component of the Poynting vector of the fundamental air core mode. We will consider this behavior using the example of a leaky air core waveguide with a core in the form of a round hole in an infinite layer of a dielectric (Fig. 2a).
The axial components of electric and magnetic fields of a fundamental air core mode (HE11 mode) in the air core ( Fig. 2a) are described by the same Bessel functions as in the case of a step -index fiber (1). Assuming the fundamental air core mode has a complex propagation constant

 
as it often happens in practice, the function argument u in (1) can be represented as  magnetic fields of the fundamental air core mode can be represented using an asymptotic of the Bessel function of 1 () Jq at 0 q  : The axial components of electric and magnetic fields of the fundamental air core mode have the same structure (4) and they are proportional to q (4). This means that the real and imaginary parts of the axial components of electric and magnetic fields simultaneously tend to zero and Fig.3a). The transverse component of the time averaged Poynting vector of the fundamental air core mode will have non -zero  -and r -projections at any point of the waveguide cross -section, except for the origin (Fig. 2a). At the origin, both projections of the transverse component of the Poynting vector component are equal to zero. Thus, the leaky fundamental air core mode has the kinetic (Abraham -type) angular momentum density [13].
For the hollow core hexagon (Fig.2b), the axial components of electric and magnetic fields of the fundamental air core mode can be represented as series of Bessel functions of the first kind of different orders with a corresponding asymptotic at the origin as in (4). The propagation constant of the fundamental air core mode is    fundamental air core mode (Fig.2b). Both phase distributions (Fig. 3) are governed by a nonlinear rule [16] and these vortices are called anisotropic. In our work [17], 3. Impact of the core -cladding boundary shape on the vortex formation.

Problem of loss reduction
In this section, we would like to consider the correlation between the vortex formation in the transverse component of the Poynting vector and phase dislocations of the fields of the fundamental air core mode, not only at the origin, but also in the waveguide cladding. Here it is important to determine of light localization mechanism in hollow core waveguides with a complex shape of the core -cladding boundary and to reduce respective waveguide losses.
To better understand the difference between the processes of the energy leakage in hollow core waveguides with various rotational symmetries of the core -cladding boundary, we consider two examples. In the first example, we examine the distribution of the radial projection of the transverse component of the Poynting vector of the fundamental air core mode in two waveguides discussed in the previous section. The results of the numerical calculation are shown in Fig. 4. In both cases shown in Fig. 4, the sign of the radial projection of the transverse component of the Poynting vector changes from zero to some positive values, taking no negative values. This is due to the fact that in both cases, the transverse component of the Poynting vector has only one vortex located at the origin. In addition, the radial projection of the transverse component of the Poynting vector in the case of a hollow core hexagon has an inhomogeneous distribution in the waveguide cross -section not related to the discrete rotational symmetry of the core -cladding boundary. The losses in both cases are approximately equal as indicated in the previous section (the imaginary parts of the propagation constants).
However, similar waveguides with a finite wall thickness show a completely different behavior. Let the waveguides considered above be located in empty space and have a wall thickness of 1 µm, so that the  Fig.4a (Fig. 5a). Only the imaginary part of the propagation constant decreases approximately by two orders of magnitude 2 (0.99958 7. 7 7) ie For the hollow core hexagon, the structure of the radial projection of the transverse component of the Poyning vector of the fundamental air core mode has, in contrast to the distribution shown in Fig. 4b, a certain discrete rotational symmetry that coincides with the symmetry of the corecladding boundary. It should be noted that, as it was shown in our work [18], in contrast to the capillary (Fig. 5a) the radial projection of the transverse Poynting vector component for the hollow core polygon has both positive (dark red) and negative (dark blue) values (Fig. 5b). For the hollow core hexagon, the imaginary part of the propagation constant is somewhat lower than for the capillary To qualitatively explain why the hollow hexagon with a finite wall thickness has the structure of the radial projection of the transverse component of the Poynting vector shown in Fig. 5b let us consider refraction and reflection of radiation from the core -cladding boundary. For the waveguides shown in Fig. 4 there is a quasi -standing wave (the air core mode) and refracted waves at the core boundary running away to infinity outside the waveguides. It is known that for an optical vortex to form under the interference of plane waves, it is necessary to have at least three waves in the considered space region [19]. Therefore, vortices in electromagnetic fields and energy flows of the fundamental air core mode for the waveguides shown in Fig.4 can be formed only at the origin, as it was shown in Section 2. A similar reasoning can be applied to the case of the capillary (Fig. 5a), namely, in addition to the air core mode and the wave going to infinity outside the waveguide, there are only two cylindrical waves in the capillary wall. For the axial components of electric and magnetic fields of the fundamental air core mode, they can be expressed as the sum of two Hankel functions:  and don't occur in the core -cladding capillary wall and, therefore, their formation is possible only at the origin. For the hollow core hexagon, the presence of a finite thickness of the core -cladding boundary combined with angular regions of the wall leads to a complex interference inside the waveguide wall, in which the axial components of electric and magnetic fields of the fundamental air core mode are described by series: The phase distributions of the axial component of the electric field of the fundamental air core mode for the capillary and hollow core hexagon are shown in Fig. 6. The Fig. 6a shows that the distribution of the axial component of the electric field of the fundamental air core mode has only one phase singularity located on the axis of the waveguide. The difference between the phases of the axial component of electric field in the hollow core and in the outer space is approximately 1.6 radians at each point of the capillary wall. This value is consistent with the values of the phase (a) (b) Figure 7. Hollow core waveguides formed by six capillaries with an outer diameter of 20 µm and a wall thickness of 1 µm inserted into the corners of the hollow core hexagon (Fig. 2b): (a) radial projection of the transverse component of the Poynting vector of the fundamental air core mode at a wavelength of 3.39 µm; (b) phase distribution of the axial electric field component of the fundamental air core mode at a wavelength of 3.39 µm. Thin red and white lines correspond to the condition that the real (red) and imaginary (white) parts of the axial component of the electric field being equal to zero. The phase changes from -π (dark blue) to π (dark red). The refractive index of the capillaries is 1.45. difference obtained from the ARROW model, when the phase difference should be equal to π/2 for the antiresonant condition in the longest wavelength transmission band. Also, there are no phase inhomogeneities along the perimeter of the capillary. For the hollow core hexagon (Fig. 6b), the phase distribution of the axial component of the electric field of the fundamental air core mode is inhomogeneous along the perimeter of the hexagon. The phase distribution is periodic in accordance with the discrete rotational symmetry of the core -cladding boundary and the phase difference between the field inside the air core and in the outer space varies in adjacent sections of the waveguide wall from exact value of π to arbitrary values of radians lying in the range from 0 to π/2 radians. It can be shown that the phase of transverse components of the electric field of the fundamental air core mode behaves in the same way. At the same time, the losses for the hollow hexagon, as it follows from the values of the imaginary parts of the propagation constants of the fundamental air core mode, are slightly lower than for the capillary.
It can be seen from the above that this a difference in the interaction of the fundamental air core mode with the core -cladding boundary wall for the capillary and the air core hexagon should lead to different types of light reflection from the waveguide wall and the air core mode energy leakage. 0 m  and the phase shift is π/2). This phase shift occurs along the entire corecladding boundary of the capillary. It is well known from singular optics that when passing through the phase dislocation of electromagnetic fields, the phase abruptly changes by π [20]. For the hollow core hexagon, the phase jump of π for the electric field of the air core fundamental mode inside the air core and outside the waveguide occurs due to phase dislocations arising in the core -cladding boundary wall at the core - light reflection from the waveguide walls occurs precisely in this phase regime and is not consistent with the ARROW model. To confirm this finding, let us consider a well -known hollow core waveguide [21] with a cladding consisting of six capillaries of the same thickness as the wall thickness of the waveguides discussed above (Fig. 5). The capillaries will be located at the vertices of the hollow core hexagon shown in Fig. 2b. Figure 7a shows the distribution of the radial projection of the transverse component of the Poynting vector of the fundamental air core mode of the waveguide. As can be seen, the distribution exactly repeats the similar distribution obtained for the hollow core hexagon (Fig. 5b) and does not coincide with the distribution in Fig. 4b. The radial projection of the transverse component of the Poynitng vector has both positive and negative values as in Fig. 5b .This suggests that the presence of thinwalled capillaries at the vertices of the hexagon leads to the same light localization mechanism as for the hollow core hexagon (Fig. 5b). In addition, the effective mode area of the waveguide shown in Fig. 7 is smaller than in the case shown in Fig. 5b To confirm that the radiation is localized in the hollow core waveguide (Fig. 7a) due to the presence of phase dislocations of the fields of the fundamental air core mode, the phase distribution of the axial component of electric field z E was plotted in Fig. 7b. It can be seen from Fig. 7b that the phase distribution of the axial component of the electric field has a vortex structure and a singularity at the origin. Also, the phase experiences a jump equal to π between the field inside the cladding capillary and outside it. This suggests that in this case the mechanism of radiation reflection and its localization in the air core is the same as for the hollow core hexagon (Fig.  6b). The difference is that for the waveguide with a cladding consisting of capillaries the phase jump equal to π occurs only in the capillary wall. Figure 8a shows enlarged image of the phase distribution of the axial component of the electric field (Fig. 7b) near and inside the capillary wall in the longest wavelength transmission band. Figure 8b shows the same distribution for the cladding capillary wall thickness of 2.42 µm (according to the ARROW model, the wavelength of 3.39 µm is antiresonant for 1 m  at this value of the capillary wall thickness). The figures clearly show the phase jumps of the magnitude of π when passing through the phase dislocations of the electric field.

Discussion
In this work, the mechanisms of light localization of the fundamental air core mode in hollow core waveguides were considered from the perspective of the phase dislocations of the mode fields and their impact on vortex formation in the transverse component of the Poynting vector. According to the ARROW model, the most efficient reflection of radiation of the air core mode from the waveguide wall can be achieved at an antiresonant condition when phase incursion in the waveguide wall is A completely different behavior occurs in hollow core waveguides with a finite wall thickness comparable to the wavelength and with a specific discrete rotational symmetry of the core -cladding boundary. In this case, the lines of zero values of the real and imaginary parts of the axial component of the electric field can intersect and overlap not only at the origin but also in the cladding forming phase dislocations of different types. The phase changes abruptly to π when passing through this phase dislocation of the air core mode electric fields and its behavior cannot be explained by the ARROW model. This behavior of the fundamental air core mode fields in the cladding can generate the corresponding vortices in the transverse component of the Poynting vector [17]. In our opinion, all of the above consideration should be taken into account for a correct understanding of light localization mechanism in hollow core waveguides and application of the ARROW model.
Funding: This research was funded by Russian Science Foundation, grant number 19 -12 -00134.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author. The data are not publicly available due to the fact that the most of the data was obtained using a commercial package COMSOL.

Conflicts of Interest:
The authors declare no conflict of interest.