Depolarization of Light in Optical Fibers: Effects of Diffraction and Spin-Orbit Interaction

: Polarization is measured very often to study the interaction of light and matter, so the description of the polarization of light beams is of both practical and fundamental interest. This review discusses the polarization properties of structured light in multimode graded-index optical ﬁbers, with an emphasis on the recent advances in the area of spin-orbit interactions. The basic physical principles and properties of twisted light propagating in a graded index ﬁber are described: rotation of the polarization plane, Laguerre–Gauss vector beams with polarization-orbital angular momentum entanglement, splitting of degenerate modes due to spin-orbit interaction, depolarization of light beams, Berry phase and 2D and 3D degrees of polarizations, etc. Special attention is paid to analytical methods for solving the Maxwell equations of a three-component ﬁeld using perturbation analysis and quantum degrees for the description of strongly focused light beams and their geometrical interpretation are also discussed.


Introduction
The polarization of light must be taken into account in many tasks of fiber optics communication and devices for coupling light in and from optical fibers. The polarization state of light does not change during propagation in a homogeneous, isotropic, nondispersed medium [1,2]. Significant changes in the state and degree of polarization occur during propagation in an inhomogeneous medium and through optical fibers.
Multimode optical fibers have attracted considerable interest in recent years for telecommunications, imaging, fiber lasers and amplifiers, ultrafast photonics, etc. They can operate over a wide range of wavelengths and have high coupling efficiency. Besides this, they degrade weakly under the influence of nonlinear effects. However, imperfections and external perturbations cause mode mixing and random polarization, leading to a strong change in the polarization state. Consequently, it is important to analyze the dynamics of the modes propagating in multimode fibers (MMFs).
Many works have been devoted to the study of depolarization of light in various media. In [3,4] the depolarization of light in randomly inhomogeneous media was considered. Two depolarization mechanisms were shown: geometric and diffraction, due to the Rytov rotation [5,6] of the polarization plane and diffraction, respectively. In single mode optical fibers, the change in the input polarization is usually due to nonlinearity and birefringence of the medium [7][8][9]. Depolarization also occurs in optical fibers without birefringence. The degree of linear polarization of radiation in multimode graded-index fibers decreases with distance [10]. In [11] it is shown that the degree of linear polarization in multimode graded-index fiber decreases with distance according to the quadratic law due to the Rytov rotation of the polarization vector, and the degree of circular polarization is preserved. However, in the experiments [12], the conservation of the degree of polarization in the optical fiber was not observed. In [13], the rotation of the plane of polarization of sagittal rays during propagation in a multimode fiber was experimentally observed. In [14], the rotation of the speckle pattern created by circularly polarized light at the fiber output was calculated, corresponding to the reversal of the sign of circular polarization. In [15], it was experimentally shown that the angle of rotation of the speckle pattern depends on the angle at which the circularly polarized light beam is coupled to the fiber. These effects can be considered as a manifestation of the optical Magnus effect [14,15] and the optical spin-Hall effect resulting from spin-orbit coupling.
In [16,17], the quantum-mechanical formalism of coherent states was used to study the evolution of polarization in a multimode isotropic graded-index medium. It is shown that depolarization of light in an isotropic graded-index medium occurs due to the effects of diffraction and spin-orbit interaction, and this mechanism manifests itself for radiation with linear and circular polarization. The rotation of the polarization vector during propagation in a two-dimensional medium is considered. It is shown that the rotation of the polarization plane depends on the axial displacement and the angle of tilt of the incident beam to the fiber axis. The change in the degree of polarization is closely related to change in the degree of coherence of the radiation. Changes in the degree of polarization of partially coherent beams propagating in optical fibers were considered in [18][19][20].
One of the goals of this review is to present a general picture of the polarizationdependent light transmission through optical fibers associated with the spin-orbit interaction. Recently, the effects of spin-orbit interaction in graded-index (GRIN) media have been of great interest [16,17,[21][22][23][24]. In [25][26][27][28], the propagation of polarized twisted light beams in a cylindrical graded-index optical fiber is investigated by solving the three-component field Maxwell equations. An operator approach to calculating the average values describing propagation of light beams is developed. Spin-dependent effects on twisted light beams propagating in an optical fiber are demonstrated by solving the full three-component field Maxwell equations using perturbation analysis. Polarization and nonparaxial effects are considered in conjunction. Mode splitting in a cylindrical graded-index optical fiber is shown using the perturbation analysis to solve the full Maxwell equations [26].
A detailed description of quantum mechanical methods, such as the coherent states method and the density matrix formalism, for considering the propagation of light in graded-index waveguides, taking into account polarization effects, is given. An operator method for studying the evolution of polarization in a graded-index optical fiber is discussed.
In addition, we emphasize that 2D and 3D degrees of polarization should be used to describe the polarization effects associated with the spin-orbit interaction, depending on whether paraxial or nonparaxial light beams are considered.

2D Polarization
Conventional laser sources emit paraxial light rays that propagate almost parallel to each other. This means that the polarization of the light beam is bounded by a plane perpendicular to the propagation direction, so for circularly polarized light, any resulting spin must be aligned longitudinally along the propagation direction. In conventional optics, the fields are transversal, so the polarization can be well described by a set of Stokes parameters defined from 2 × 2 coherence matrix [1]. In [2], a unified theory of coherence and polarization of random electromagnetic fields is proposed, which can be used to investigate changes in polarization during the propagation of light beams.
The Maxwell equations for the electric field → E exp(−iνt) in an inhomogeneous medium are reduced to: or where k = 2π/λ is the wavenumber and ε = n 2 is the dielectric permittivity of the medium. In the paraxial approximation the Maxwell equations can be reduced to the equivalent time-independent Schrodinger equation [29]. In [16,17], a similar approach is used to obtain a parabolic equation for a two-component vector field. Since the paraxial beam propagates in a direction close to the z-axis, the conditions ∂ 2 E/∂z 2 << k∂E/∂z ≈ ∂ 2 E/∂x 2 ≈ ∂ 2 E/∂y 2 are valid. We have the following equation for the two-component wave function of the field, neglecting the term ∂ 2 E/∂z 2 in Equation (2) [16,17]: where ψ = n ∂x 2 + ∂ 2 ∂y 2 + 1 2n 0 n 2 0 − n 2 Î is the unperturbed Hamiltonian corresponding to the first two terms in Equation (2) and The HamiltonianĤ can be rewritten in terms of annihilationâ and creationâ + operators using the relationsx Thus, we haveĤ 0 = ω kn 0 â + 1â 1 +â + 2â 2 + 1 Î , H 1 = η Âσ z +Bσ x + kM zσy +ĈÎ +Î .
The Pauli matrices satisfy the following relations: The solution to Equation (7) can be expressed as The polarization properties of the light are determined by the elements of the coherency matrix [1,2]: The evolution of the coherency matrix is determined by the expression [16,17]: whereĴ (0) is the coherency matrix operator in the initial plane z = 0.

Rotation of Polarization Plane in a Graded Index Fiber
The polarization plane rotates on propagation of light on the helical trajectory [5,6]. In a single mode optical fiber wound on a cylinder, such rotation was observed experimentally [30] and interpreted in terms of Berry's geometrical phase [31]. The rotation of the polarization plane was observed also in a straight multimode fiber with a step-index type profile [13]. It is also interesting to consider the reverse effect, i.e., the effect of polarization on the trajectory and width of the beam. When a light beam is reflected from an interface, the longitudinal shift of the gravity center of the beam is different for sand p-polarized beams [32], while the transverse shift has reverse signs in the case of a right-and left-hand circularly polarized light beam [33]. In [34], lateral and angular shifts were shown for strongly focused azimuthally and radially polarized beams at the interface of dielectrics. In [14], the rotation of the speckle pattern created by circularly polarized light at the fiber output was shown, corresponding to the reversal of the sign of circular polarization. The dependence of the angle of rotation of the speckle pattern on the angle at which a circularly polarized light beam falls on the fiber input was experimentally demonstrated in [15].
It was shown in [35] that when light propagates along a helical trajectory in a gradedindex fiber, a spin-dependent relative shift occurs between right-hand and left-hand circularly polarized light beams. This effect was observed for a laser beam propagating in a glass cylinder along a smooth helical trajectory [36]. Note that this shift is a manifestation of the optical Magnus effect [37] and the optical spin-Hall effect [38,39], which occurs due to the spin-orbit interaction. Beam and wave optics are used to analyze the propagation of light in graded-index media [40]. The effect of polarization on modes in lens-like media was analyzed in [41]. In [42], the polarization-dependent GH beam shift at the graded-index interface of dielectrics is studied. In [43], the beam shifts with respect to geometric optics caused by spin-orbit interaction and nonparaxiality in a graded-index optical fiber are investigated.
The methods of coherent states can be used to study the evolution of the parameters of a light beam. The operator approach is used to calculate the average values describing the light beam. These methods were used in [16,17] for investigation into the depolarization of light in a graded-index waveguide. A two-dimensional graded-index waveguide was considered: where ω is the gradient parameter, n 0 is the refractive index along the axis of the waveguide, and x and y are the transverse coordinates of the waveguide. In the input plane z = 0, we introduce coherent states represented by Gaussian wave packets: where define the initial coordinates (x 0 , y 0 ) of the center of gravity of the light beam and the tilt angles of the beam p x0 = n 0 sin θ x , p y0 = n 0 sin θ y with respect to the axis of the medium. Coherent states were first introduced by Glauber [44] in 1963 from the consideration of the states of the electromagnetic field oscillators and the study of its statistical properties. In essence, these states are analogous to the Gaussian wave packets, which were constructed and studied by Schrodinger [45] for the quantum harmonic oscillator as part of his study of the connection between the quantum and classical descriptions.
In [46,47], the coherent states method, the integral of motion and the density matrix formalism are used to describe the propagation of a paraxial optical beam and partially coherent radiation in a weakly inhomogeneous medium. The coherent states method is also valuable for considering the diffraction of partially coherent light beams by microlens arrays [48].
Coherent states (10) have the minimal width and diffraction-limited angular divergence on the propagation in a square-law-index medium. The center of such a wave packet propagates along the trajectory of the geometric ray corresponding to the ray optics. Besides, the coherent states (10) are generating functions for the waveguide modes and they form a complete set of functions. This property can be used to decompose an arbitrary beam field | f into a series of coherent states: The trajectory and width of the beam are defined in terms of the matrix elements: The wave functions of the light field linearly polarized along the xand y-axes can be expressed as The wave functions of right-and left-hand circularly polarized light have the form The evolution of the operatorâ 1,2 is determined by the equation whereĤ is the Hamiltonian of the system. For a light beam linearly polarized along the x-axis with initial parameters x 0 = 0, y 0 = 0, p x0 = 0, and p y0 = 0, we obtain the following expressions for the beam trajectory and beam width [35]: It follows from Expressions (16) and (17) that the polarization of the radiation does not affect the path of the beam. However, the polarization effects affect the beam width, which fluctuates with a period equal to the beam oscillation period.
Additional fluctuations in the beam width in the direction of field polarization are shown in Equations (16) and (17). This effect is a first-order correction to the paraxial approximation. In [49,50] the dynamics of beam self-imaging by means of femtosecond laser pulse propagation in nonlinear GRIN multimode optical fibers was studied experimentally. It was shown that the periodical beam self-imaging acts symmetrically on both the x and y axis. The amplitude of the oscillations for the parameters of the graded-index fiber used in [49,50] is only 0.3% of the initial beam width, so it is likely that the effect was not observed in experiments.
For a circularly polarized light beam, the expressions for the trajectory and beam width with the initial parameters x 0 = 0, y 0 = 0, p x0 = 0, and p y0 = 0 have the form [35]: where σ = +1 corresponds to right-hand circularly polarized radiation and σ = −1 corresponds to the left-hand circularly polarized light beam. For circularly polarized light, the plane of propagation of the meridional beam (p y0 = 0) rotates in the direction of circulation of the polarization vector. A similar rotation takes place for the sagittal rays (p y0 = 0). The rotation of the meridional ray and the additional rotation of the sagittal rays increase linearly with distance: The following is a comparison between the Rytov rotation angle or the Berry phase of the polarization plane due to the twisting of the propagation path [17]: the rotation angle of the plane of the meridional ray (19) is due to the circular nature of the polarization of light. It follows from (19) and (20) that the values of these angles are equal to each other for x 2 0 = 1/kω, i.e., in the case when the axial shift of the incident beam is equal to the width of the fundamental waveguide mode w 0 = √ 2/kω. In Figure 1, the angular shift dependence on the beam's angular aperture (diffraction angular divergence) ∆θ = λ/πw 0 for the distance ωz = 2π corresponding to the period of the beam twisting is presented. The angular GH shift increases in proportion to the square of the angular aperture of the beam. A similar effect was experimentally observed for a wave beam when a Gaussian beam falls on the air-glass interface [51]. Substantial positive and negative Goos-Hänchen shifts near the surface plasmon resonance were also shown in subwavelength gratings [52]. The polarization affects the width of the beam propagating in the waveguide. Below, we consider the paraxial propagation of a beam whose initial width is not matched with the fundamental mode of the waveguide. The analysis of the propagation of a strongly focused light beam in a graded-index waveguide is presented in [53,54].
For a linearly polarized axis beam (α = 0) with a field radius a0, the evolution of the beam width is defined by the expressions [43]: The minimum beam width of a circularly polarized beam is given by It follows that the values of the minimum beam width decrease if the width of the incident beam is greater than the radius of the fundamental mode of the waveguide ( 00 aw  ).
Thus, the polarization effects cause an increase in the size of the focal spot. For linearly polarized light, the beam shows a broadening only in the direction of the plane of polarization, and the initially symmetric beam becomes elliptical. Such a broadening of the beam in the direction of the polarization plane was observed in experiments devoted to measuring the size of latex spheres [55]. The polarization affects the width of the beam propagating in the waveguide. Below, we consider the paraxial propagation of a beam whose initial width is not matched with the fundamental mode of the waveguide. The analysis of the propagation of a strongly focused light beam in a graded-index waveguide is presented in [53,54].
For a linearly polarized axis beam (α = 0) with a field radius a 0 , the evolution of the beam width is defined by the expressions [43]: where , and ω + = 2 kw 2 0 is the gradient parameter, w 0 is the radius of the fundamental mode. At z = n 0 π/(2ω), the beam width takes the minimum value: The minimum beam width of a circularly polarized beam is given by It follows that the values of the minimum beam width decrease if the width of the incident beam is greater than the radius of the fundamental mode of the waveguide (a 0 > w 0 ).
Thus, the polarization effects cause an increase in the size of the focal spot. For linearly polarized light, the beam shows a broadening only in the direction of the plane of polarization, and the initially symmetric beam becomes elliptical. Such a broadening of the beam in the direction of the polarization plane was observed in experiments devoted to measuring the size of latex spheres [55].
The effect of polarization effects on the beam parameters is negligible in the case when the ratio of the gradient parameter to the wave number is small. However, when the beam is focused to a region smaller than the radiation wavelength, these effects become significant. This should be taken into account when analyzing experimental data, in particular the results of measurements using near-field microscopes.
Circular polarization leads to rotation of the plane of propagation of the meridional beam, as well as to additional rotation of the sagittal rays, and this rotation depends on the axial displacement of the incident beam. Both circular and linear polarizations lead to fluctuations in the beam width during propagation in a graded-index waveguide. Polarization affects the focus of the light beam in a graded index medium and is a limiting factor. In the case of linear polarization, the cross-section of the light beam in the plane of focus becomes elliptical, elongated in the direction of the plane of polarization.

Depolarization of Light
Depolarization of light propagating in a graded-index fiber was considered in [16,17]. It is shown that the spin-orbit interaction causes an asymmetry effect for the depolarization of right-hand and left-hand circularly polarized light. Depolarization is stronger if the helicity of the trajectory and the photons have the same signs, and weaker if they are not the same.

Linear Polarization
The fields of the linearly polarized beams are given in Equation (13). The coherency matrix is determined by the following matrix elements: The coherence matrix operator in the plane z = 0 is defined by the source. For radiation linearly polarized along the x-axis, we have [1] where I 0 is the total intensity of the incident beam. For linearly polarized radiation along the y axis the coherency matrix operator has the form [ To solve Equation (8) and calculate the matrix elements (24), the commutation relations â i ,â j = δ i,j (i, j = 1, 2), the relationships (5) and the relationŝ can be used.
Solving Equation (8) up to a parameter η 2 and substituting the solutions in Equation (7), we obtain the following expression for the value, determining the depolarization [16,17]: Equation (28) includes a contribution of both sagittal and meridional rays with initial coordinates in the plane z = 0: x 0 = 0, y 0 = 0, p x0 = 0, p y0 = 0.
It follows that the depolarization depends on the wavenumber and disappears at λ → 0. For an axial beam (x 0 = 0, p y0 = 0), we obtain oscillations in the degree of polarization that are purely diffractive in origin. The effect of periodic retrieval of the initial degree of polarization should be observed in a single-mode isotropic optical fiber, when the fundamental mode of the optical waveguide is considered.
The dependence of the degree of polarization of linearly polarized light on the distance z is represented by the curves in Figure 2 for various offsets x 0 , including the meridional rays and the radius-preserving (p y0 = x 0 ω) spirally twisted sagittal rays. It is clear that the degree of polarization decreases inversely with the square of the distance. The decrease in the degree of polarization becomes more pronounced as the displacement or inclination of the beam relative to the fiber axis increases. The depolarization length of the meridional rays is 1.4 times longer than the depolarization length of the sagittal rays. Equation (28) includes a contribution of both sagittal and meridional rays with initial coordinates in the plane z = 0: x0 ≠ 0, y0 = 0, px0 = 0, py0 ≠ 0.
It follows that the depolarization depends on the wavenumber and disappears at λ → 0. For an axial beam (x0 = 0, py0 = 0), we obtain oscillations in the degree of polarization that are purely diffractive in origin. The effect of periodic retrieval of the initial degree of polarization should be observed in a single-mode isotropic optical fiber, when the fundamental mode of the optical waveguide is considered.
The dependence of the degree of polarization of linearly polarized light on the distance z is represented by the curves in Figure 2 for various offsets x0, including the meridional rays and the radius-preserving ( 00 ) spirally twisted sagittal rays. It is clear that the degree of polarization decreases inversely with the square of the distance. The decrease in the degree of polarization becomes more pronounced as the displacement or inclination of the beam relative to the fiber axis increases. The depolarization length of the meridional rays is 1.4 times longer than the depolarization length of the sagittal rays. In Figure 3 the values F = 10 5 4detJ/Sp 2 J as a function of the distance z are presented for meridional and sagittal rays. It can be seen that the increase in this value is accompanied by small fluctuations. The degree of polarization decreases quadratically with the increasing distance. Depolarization increases with an increase in the displacement of the axis or the inclination of the beam to the axis of the waveguide. The decrease in the degree of polarization is accompanied by small fluctuations. Note that Equation (28) is valid only for distances z ≤ z 0 , where the distance z 0 is determined by the accuracy of the solution of Equation (8). In our case, this distance is equal to The degree of polarization can be obtained by substituting Equation (28) into Equation (7): The Equation (29) at distances z >> ω −1 has the form [13,14]:  The degree of polarization can be obtained by substituting Equation (28) into Equation (7): The Equation (29) It follows that the degree of polarization of the axial beam takes its initial value at ωz = mπ (m = 1, 2, 3 …). The depolarization characteristics of the imaging probe with the multi-mode fiber with a length of 6.2 mm were studied in [56]. It is shown that the depolarization coefficient for the modal fields LP01 and LP21 can be estimated at 0.7.

Circular Polarization
The fields for right-handed and left-handed circularly polarized light are defined in Equation (14). The coherency matrix in the initial plane has the form [1] ( ) The Jones vector describing linearly polarized light can be defined in terms of the Jones vector of circularly polarized radiation: The degree of polarization decreases with increasing distance, and depolarization disappears at λ → 0. There are fluctuations in the degree of polarization of purely diffractive origin in the case of the axial beam (x 0 = 0, p y0 = 0): It follows that the degree of polarization of the axial beam takes its initial value at ωz = mπ (m = 1, 2, 3 . . . ). The depolarization characteristics of the imaging probe with the multi-mode fiber with a length of 6.2 mm were studied in [56]. It is shown that the depolarization coefficient for the modal fields LP 01 and LP 21 can be estimated at 0.7.

Circular Polarization
The fields for right-handed and left-handed circularly polarized light are defined in Equation (14). The coherency matrix in the initial plane has the form [1] The Jones vector describing linearly polarized light can be defined in terms of the Jones vector of circularly polarized radiation: The following expression describing the depolarization can be obtained by solving Equation (8) up to a small parameter η 2 [16,17]: 4detJ y0 ω 2 z 2 + sin 2 ωz ± 4kx 0 p y0 ωz sin 2 ωz + 2 sin 2 ωz .
The signs ± correspond to right-handed and left-handed circularly polarized light. The effect of asymmetry occurs with respect to the sign of twisting of the trajectory of sagittal rays (Figure 4). Depolarization is stronger if the helicity of the trajectory of rays and "photons" has the same signs, and less if they do not coincide. The depolarization of the meridional rays is less than that of the sagittal rays.
Equation (8) up to a small parameter η 2 [16,17]: The signs ± correspond to right-handed and left-handed circularly polarized light. The effect of asymmetry occurs with respect to the sign of twisting of the trajectory of sagittal rays (Figure 4). Depolarization is stronger if the helicity of the trajectory of rays and "photons" has the same signs, and less if they do not coincide. The depolarization of the meridional rays is less than that of the sagittal rays. The following relations are used to calculate the coherency matrix elements. Substituting Equation (33) into (7), we obtain an expression for the degree of polarization The degree of polarization at the distances z >> ω −1 may be defined as The following relationŝ are used to calculate the coherency matrix elements. Substituting Equation (33) into (7), we obtain an expression for the degree of polarization The degree of polarization at the distances z >> ω −1 may be defined as The sagittal rays with constant radius of twisting (p y0 = x 0 ω) are considered here. The degree of polarization decreases with increasing distance, similar to the degree of polarization of a linearly polarized beam. A similar result occurs in the case of left-handed polarized light.

Berry Phase and Degree of Polarization
The plane of polarization of a light beam rotates when it propagates along a curved trajectory in an inhomogeneous medium [5,6]. This rotation was experimentally demonstrated [30] in a single-mode fiber wound along a cylinder and interpreted as a geometric Berry phase effect [31]. The effect of a spin-dependent relative shift between right-and left-handed circularly polarized light beams propagating in a glass cylinder along a smooth helical trajectory was also observed [36]. A similar beam shift was shown in a graded index optical fiber [43]. This shift can be considered as a manifestation of the optical Mag- nus effect [14] and the optical spin-Hall effect [38,39], which occurs due to the spin-orbit interaction.
In [57], the evolution of the polarization plane and the Berry phase in a graded index waveguide was investigated using the quantum mechanical formation of coherent states.
The rotation angle of the polarization vector of a beam linearly polarized along the x-axis propagating along a helical trajectory is determined by the expression [57]: The Expression (38) is valid for sagittal rays with initial coordinates For small values of the axial displacement of the beam, the Expression (38) has the form [57]: It can be seen from (39) that the rotation angle of the polarization plane increases quadratically with the angle of tilt θ (p y0 = n 0 sinθ) of the beam to the waveguide axis or with the axis displacement x 0 for sagittal rays with constant radius of twisting (p y0 = x 0 ω) (Figure 5a,b). The calculated results for the angle of rotation of the polarization plane are in good agreement with the experimental data [13]. Note that the cubic dependence of the rotation angle of the polarization vector on the angle between the propagation direction and the axis of the fiber with a step profile of the refractive index was obtained in [37]. The rotation speed of the polarization vector is defined by the expression ( ) It follows that the rotation of the polarization vector is of irregular character. The rotation speed is zero at the points where the polarization vector is parallel to the main normal. Note that the Rytov's rotation [5] or the Berry's phase [31] have a geometrical origin and do not depend on the wavelength of the radiation. However, the variance of the Berry's phase is of wave nature and disappears at λ → 0. Thus the variance of the Berry's phase is equal to [57] It can be seen that the dispersion of the Berry's phase increases nonlinearly at small values of the axis displacement. The Berry's phase dispersion for It follows that the rotation of the polarization vector is of irregular character. The rotation speed is zero at the points where the polarization vector is parallel to the main normal. Note that the Rytov's rotation [5] or the Berry's phase [31] have a geometrical origin and do not depend on the wavelength of the radiation. However, the variance of the Berry's phase is of wave nature and disappears at λ → 0. Thus the variance of the Berry's phase is equal to [57] (δβ) 2  It can be seen that the dispersion of the Berry's phase increases nonlinearly at small values of the axis displacement. The Berry's phase dispersion for x 0 = 20 µm of the order of <δβ> ≈ 1 • occurs at a distance of z ≈ 3.7 mm. This corresponds to the rotation of the polarization plane by an angle of <β> ≈ 4.4 • .
It follows that the dispersion of the Berry's phase is determined by the degree of polarization P. As noted in [29], the depolarization in an optical fiber with a length of z ≈ 7.5 cm is 10-30%, which is consistent with the results obtained ( Figure 6). It can be seen that fluctuations in the Berry's phase dispersion for the axial beam exist, although the phase is zero. In addition, the Berry's phase fluctuations for the meridional rays increase with distance, although the Berry's phase value is zero. Thus, measurements of the degree of polarization will verify experimentally the wave nature of the Berry's phase. The dispersion of the rotation angle of the polarization plane increases with distance and can be determined by measuring the degree of polarization of light. The dispersion of the Berry's phase or the angle of rotation of the polarization plane increases with distance, and there are fluctuations in the Berry's phase dispersion for the axial beam (fundamental mode), which have a wave nature.

3D Polarization
In conventional optics, the fields are transverse, so the polarization can be well described by a set of Stokes parameters determined from the 2 × 2 coherence matrix [1,2]. However, the electromagnetic field in the far zone is transverse and for waves propagating along the z axis, Ez = 0, many problems of near-field optics, optical microscopy, data storage, and scattering need to be analyzed, taking into account the longitudinal component of the field Ez. Strongly focused beams contain a non-zero longitudinal component of the field 0 z E  , so we should consider 3 × 3 coherence matrices to describe them. The Stokes 2 × 2 formalism was extended to 3 × 3 coherence matrices, and the degree of polarization for arbitrary electromagnetic fields was introduced in [58 -61]. The evolution of three-dimensional electromagnetic fields in a graded index medium is studied in [25][26][27][28]62]. The coherence matrix is represented in terms of spherical or rectangular tensors, similar to the Gell-Mann matrices. The evolution of the coherence matrix is defined in terms of the Hamiltonian obtained from the three-component Maxwell field equations. The polarization state of the most general two-dimensional field can be completely described by three parameters. This is referred to as vector polarization. For a three-dimensional field, except for three components of the vector polarization, there are five components of tensor polarization. The complexity of the 3D electromagnetic field  The dispersion of the rotation angle of the polarization plane increases with distance and can be determined by measuring the degree of polarization of light. The dispersion of the Berry's phase or the angle of rotation of the polarization plane increases with distance, and there are fluctuations in the Berry's phase dispersion for the axial beam (fundamental mode), which have a wave nature.

3D Polarization
In conventional optics, the fields are transverse, so the polarization can be well described by a set of Stokes parameters determined from the 2 × 2 coherence matrix [1,2]. However, the electromagnetic field in the far zone is transverse and for waves propagating along the z axis, E z = 0, many problems of near-field optics, optical microscopy, data storage, and scattering need to be analyzed, taking into account the longitudinal component of the field E z . Strongly focused beams contain a non-zero longitudinal component of the field E z = 0, so we should consider 3 × 3 coherence matrices to describe them. The Stokes 2 × 2 formalism was extended to 3 × 3 coherence matrices, and the degree of polarization for arbitrary electromagnetic fields was introduced in [58 -61]. The evolution of threedimensional electromagnetic fields in a graded index medium is studied in [25][26][27][28]62]. The coherence matrix is represented in terms of spherical or rectangular tensors, similar to the Gell-Mann matrices. The evolution of the coherence matrix is defined in terms of the Hamiltonian obtained from the three-component Maxwell field equations.
The polarization state of the most general two-dimensional field can be completely described by three parameters. This is referred to as vector polarization. For a three- dimensional field, except for three components of the vector polarization, there are five components of tensor polarization. The complexity of the 3D electromagnetic field arises from the tensor polarization [59,60]. In contrast to the 2D case, where any beam of light can be represented as the superposition of totally polarized and unpolarized light, in the 3D case, an arbitrary beam cannot be represented as a mixture of an unpolarized beam (one parameter) and a polarized beam (five parameters) because nine independent parameters are needed to characterize the coherence matrix.
Maxwell's equations for the electric field → E exp(−iνt) in a general inhomogeneous medium with dielectric constant ε(x, y) are reduced to: where k = 2π/λ is the wavenumber and n 2 (x, y) is the dielectric constant of the medium, β is the propagation constant.
Here the dependence on time and z is assumed as exp[−i(νt − βz)]. Rewrite the system (42) as [27,62]: where Formally, taking the square root of both sides, we obtain an equation equivalent to the stationary Schrodinger equation for the reduced field Ψ(x, y): where ε = δβ = β − kn 0 and Ψ(x, y) are the eigenvalue and eigenfunction of the Hamiltonian, accordingly, n 0 = n(0, 0) > H (1) is the nonparaxial correction to the operator of paraxial propagation, which has the form The Hamiltonian > H can be rewritten in terms of annihilation and creation operators in cylindrical coordinates [17,27,62]. This allows the averaged values describing the light beams to be calculated analytically.

Vector Laguerre-Gauss (LG) Beams in a Graded-Index Fiber
Diffraction and refraction forces govern a light beam propagating in an inhomogeneous medium in the scalar approximation. However, the vector Maxwell equations are the basis of a consistent electromagnetic theory of light propagation. To analyze the vector wave fields, the influence of additional effective forces associated with polarization (spin angular momentum (SAM)) and orbital angular momentum (OAM) should be taken into account. It is known that when the linearly polarized optical beam propagates along a spiral trajectory, the polarization vector undergoes a Rytov rotation [5,6].
In [28], the full three-component field Maxwell equations are solved to investigate the effect of polarization (spin) and OAM on the shape of a vortex light beam with a nonzero radial number propagating through a cylindrical graded-index fiber. Remote focusing and highly efficient transfer of a strongly focused spot in a graded-index waveguide are demonstrated in [54]. It is known that different paraxial modes or rays are periodically in phase as they propagate in graded-index optical waveguides, and periodic focusing can be observed [63]. However, nonparaxial effects destroy this periodic focusing. Surprisingly, the restoration of the initial beam shape and strong focusing occur at extremely large distances due to the revival effect caused by interference between propagating modes [28,54]. Expressions of the path and beam width describing nonparaxial propagation in gradedindex (GRIN) media were obtained earlier in [64]. However, the resulting expressions do not include the revival effect. In [65,66], the correct expressions for the beam trajectories and widths are obtained, which contain the long-term period oscillations of the beam trajectory and its width. In [67], the revival effect in GRIN media was considered up to a second-order approximation. In [53,54,[68][69][70], exact analytical expressions for the trajectory and width of nonparaxial wave beams are obtained. It should be noted that the revival effect in a cylindrical GRIN rod, which occurs during nonparaxial propagation, was also recently studied in [71]. However, the problem was considered in the framework of the scalar Helmholtz equation, which is not the case of nonparaxial propagation of a light beam in a cylindrical waveguide. In addition, Hermite-Gauss beams were considered as modal solutions, whereas for a cylindrical GRIN waveguide they may not be modal solutions.
Analytical solutions can be obtained for a rotationally symmetric graded-index optical fiber: where n 0 is the refractive index on the waveguide axis, ω is the gradient parameter, r = x 2 + y 2 , a is the waveguide radius. The unperturbed Equation (44) has the form where the HamiltonianĤ 0 is given in (45). The solution of the unperturbed equation is described by the radially symmetric Laguerre-Gauss functions (LG) ψ vl (r, φ) = |v, l (Figure 7) [62]: where v = 2p + l is the principal quantum number, p and l are the radial and azimuthal indices, accordingly, and l = v, v − 2, v − 4, ...1 or 0, ω = 2/(kw 2 0 ), w 0 is the radius of the fundamental mode. The numbers v and l express the eigenvalues of the unperturbed Hamiltonian > H 0 |v, l = ω/kn 2 0 (v + 1)|v, l , and eigenvalues L = l/k of the angular momentum operator L z |v, l = (l/k)|v, l .
The hybrid wave functions consisting of transverse and longitudinal components are solutions of Equation (43): where σ = +1 and σ = −1 correspond to right-handed and left-handed circularly polarized beams, accordingly, and σ = 0 corresponds to the linear polarization.
where 1  =+ and 1  =− correspond to right-handed and left-handed circularly polarized beams, accordingly, and 0  = corresponds to the linear polarization. There is no mode conversion at propagation if the incident beam is expressed by (48). Note that the hybrid wave function (48) It follows from (49) that the group velocity There is no mode conversion at propagation if the incident beam is expressed by (48). Note that the hybrid wave function (48) cannot be factorized into the product of spin and orbital parts since the mixing of OAM and SAM exists. Thus, the modal solutions of the Maxwell equations in a GRIN media are the hybrid vector Laguerre-Gauss modes with the spin-orbit entanglement. The longitudinal field component can be expressed through the transverse field components, i.e., |e z = (i/kn 0 )∇ ⊥ e ⊥ . Note that besides SAM and OAM, the LG beams possess the radial number associated to the intensity distribution in the transverse plane [72][73][74].
The propagation constant that is correct to the first-order nonparaxial term has the form [23]: where η = ω/kn 2 0 , j = l + σ. It follows from (49) that the group velocity v g = 1 c dβ dk −1 depends on the OAM and SAM of the propagating pulsed beam [26]. Analogically, the speed of the vortex light beam in the free space also depends on the OAM value [75].
Below, we consider the vector vortex beams with right-circular and left-circular polarizations, accordingly: Ψ + 0 = ( v l|, i v l|, e z ) and Ψ − 0 = ( v l|, i v l|, e z ), where |v l is given by (47), and ω = 2/(ka 2 0 ), a 0 is the radius of a beam, which is different from the radius of the fundamental mode of the medium w 0 = 2/(kω).
The incident beam can be decomposed into modal solutions, so the evolution of the beam in the medium (46) can be represented as [27]: where a vlσ are the coupling coefficients. Only propagating modes are considered below, i.e., all the propagating power is inside a waveguide, the evanescent waves do not reach the far-field zone [76].

Simulation Results
The intensity distribution changes with distance are determined by the functions I ⊥ (r, ϕ, z) = |Ψ ⊥ (r, ϕ, z)| 2 and I z (r, ϕ, z) = |e z (r, ϕ, z)| 2 . Figure 8 shows the intensity distributions of the linearly polarized beams with the radial number p = 1 and different OAMs in the focal plane. The gradient parameter ω = 7 · 10 −3 µm −1 of the medium (46) and the refractive index n 0 = 1.5 are considered. The beams with wavelength λ = 0.63 µm are considered. The numerical aperture is determined by N A = aω = n 0 √ 2∆, where a is the waveguide radius and ∆ ≈ (n 0 − n(a))/n 0 . The initial beam width is a 0 = 45 µm. As can be seen, the intensity distributions depend on the SAM and OAM of the incident beam.
The focused spot in the longitudinal field is split into two equal parts if l = 0 (Figure 8b). An asymmetry occurs between the intensity distributions of the longitudinal field component for incident beams with opposite OAM (Figure 8f,h). There is no such asymmetry for the transverse field components.
Note that the modes exhibiting hybrid entanglement between SAM and OAM may be useful for classical implementations of quantum communication and computational tasks [77].

The Splitting of the Degenerate Modes
In [23], the splitting of modes in a graded-index optical fiber due to the spin-orbit interaction is shown by solving the Maxwell equations using perturbation analysis. Due to the axial symmetry of the fiber, the modes are (v + 1)-fold degenerated, and for a given value of principal number v = 2p + l, where p and l are the radial and azimuthal indices, the modes are degenerate by l. The modes are 2(v + 1)-fold degenerated if the polarization degeneracy of the modes is taken into account.
The polarization degeneracy is eliminated by an additional term in the Hamiltonian describing the spin-orbit interaction. It is shown in [26] that the first order correction to the eigenvalue ε (1) = vlσ|H (1) |vlσ depends on the sign of the helicity and the total angular momentum where vlσ| = vl|, −iσ vl|, e z |, η = ω/kn 2 0 is the small parameter, j = l + s is the total angular momentum, σ is the helicity, s is the spin angular momentum, and the longitudinal field component can be expressed through the transverse field components, i.e., |e z = (i/kn 0 )∇ ⊥ e ⊥ .
The last term in (52) refers to the spin-orbit interaction and is analogous to the correction obtained for the modes of a ring resonator from a circular dielectric waveguide [78]. However, in their calculation results, the terms corresponding to the nonparaxiality were omitted. In fact, the first two terms in Equation (52) do not contribute to splitting a level with the same principal number and total angular momentum.
Using (49) the group delays of the modes are expressed by where z is the length of the fiber, c is the velocity of light in free space. The group velocities v g = z/τ of the vortex modes with right-and left-handed polarizations differ from each other, so the effective anisotropy is induced due to the spin-orbit interaction. In the case of zero orbital momentum, there is no such asymmetry.

Simulation of Delay Time and Group Delay Splitting
In Figure 9a the dependences of the propagation delay time of modes in comparison with the fundamental mode on the radial mode number for different fundamental mode radii w 0 = (2/kω) 1/2 are presented. In Figure 9b the group delay splitting of the azimuthal modes of fixed radial indices (p = 0) is presented. It follows that the group delay splitting is an order of magnitude smaller than the relative propagation delay of modes. Laguerre-Gauss beams are also modes of the Fabry-Perot resonator and fiber ring resonators made of a graded-index fiber. It follows from the resonance condition that the frequency shift and the corresponding wavelength shift are expressed by [26]: 8 where 0 / f cn  = .
Note that the frequency shift of the same order of magnitude was obtained in [79,80]. In Figure 10a the wavelength shifts  In addition to the polarization degeneracy, there is the degeneracy of modes with different orbital angular momentum and polarization, but with the same total angular momentum. In this case, corrections to the eigenvalues cannot be obtained from (52). The eigenvalues of the modes can be obtained analytically using the degenerate perturbation Laguerre-Gauss beams are also modes of the Fabry-Perot resonator and fiber ring resonators made of a graded-index fiber. It follows from the resonance condition that the frequency shift and the corresponding wavelength shift are expressed by [26]: where f = cn 0 /λ. Note that the frequency shift of the same order of magnitude was obtained in [79,80]. In Figure 10a the wavelength shifts ∆λ l = λ l+1 − λ l as a function of the azimuthal mode number for different beam radii are presented. The splitting is very sensitive to the size w 0 of the fundamental mode and increases with the azimuthal index. A similar result was obtained in almost hemispherical micro-cavities experimentally [81]. Laguerre-Gauss beams are also modes of the Fabry-Perot resonator and fiber ring resonators made of a graded-index fiber. It follows from the resonance condition that the frequency shift and the corresponding wavelength shift are expressed by [26]: where 0 / f cn  = .
Note that the frequency shift of the same order of magnitude was obtained in [79,80]. In Figure 10a the wavelength shifts   In addition to the polarization degeneracy, there is the degeneracy of modes with different orbital angular momentum and polarization, but with the same total angular momentum. In this case, corrections to the eigenvalues cannot be obtained from (52). The eigenvalues of the modes can be obtained analytically using the degenerate perturbation  In addition to the polarization degeneracy, there is the degeneracy of modes with different orbital angular momentum and polarization, but with the same total angular momentum. In this case, corrections to the eigenvalues cannot be obtained from (52). The eigenvalues of the modes can be obtained analytically using the degenerate perturbation theory [26]: There is a degeneracy lifting δε = (5/8)η 2 (v + 1) (p + 1)(p + l) of the levels with the same total angular momentum j. The hybrid modes, which are solutions of Equation (44), correspond to these levels The vector modes of a graded-index fiber are hybrid polarization-OAM-entangled states. Note that a similar OAM mixing of LG modes was also found in micro-cavities [82]. Because of the fact that the shift also takes place for the level with zero OAM, its mechanism is not only connected with the spin-orbit interaction, but also with nonparaxiality and tensor interaction [61]. The degeneracy lifting can be considered as an optical analogue of the Lamb shift, in which there is a separation of levels between degenerate states with the same total angular momentum.
The relative degeneracy lifting δλ l /λ for different beam radii w 0 is shown in Figure 10b. The degeneracy lifting affects the quality Q-factor, which characterizes the resonant bandwidth of the modes. For resonators with the radius of the fundamental mode of the order of the wavelength of light, this splitting is about 1 nm, which can be observed experimentally.

Vector and Tensor Polarization Degrees of Light
The polarization state of the most general two-dimensional field can be fully described by three parameters. This is related to as vector polarization. For a three-dimensional field, in addition to the three components of the vector polarization, there are five components of tensor polarization.
Many works have been devoted to the description of the polarization of a threedimensional field, in particular, several measures of the 3D degree of the polarization P have recently been proposed [83][84][85][86][87][88][89][90]. There is no unambiguous mathematical definition of the degree of polarization, that is, different measures of P are not equivalent, which indicates that a single parameter is not enough to fully describe all the properties of a threedimensional field. The problem is that all measures of the 3D degree of polarization cannot distinguish between the two fields with different polarization matrices [84,85]. It is shown in [59] that two measures of the degrees of polarization (vector and tensor) are sufficient for a complete description of three-dimensional polarization, i.e., for distinguishing two fields with different polarization matrices. Examples illustrating various possible polarizing states are considered.
The elements of the Stokes vector are the intensity, vector and tensor polarizations which are defined in terms of the density matrix [61,62] (coherence matrix) ρ: whereT KQ are the spherical tensor operators, K = 1, 2; Q = 0, ±1, ±2.
Using the Equation (56) the density matrix may be written in terms of the expectation values t KQ : where the scalar, vector and tensor components corresponding to monopole, dipole and quadrupole radiations are separately presented. The spherical tensor operators in terms of components of the spin operatorσ are [61,62]:τ The advantages of this representation are that these irreducible operators are transformed in the polarization space in the same way that spherical harmonics are transformed in the coordinate space [91], and that the measured light beam intensities can simply be expressed in terms of their mathematical expectations. The correspondence of the tensor elements t KQ to the field correlation components E * i E j is easy to find from the expansion (57). Note that spherical operators can also be defined in terms of the Gell-Mann operators [58,59,83]. However, unlike the decomposition of three-dimensional fields onto Gell-Mann matrices, this decomposition has a physically clear interpretation.
The degree of polarization is described by the expression It is reasonable to consider the vector and tensor polarization components in (59) separately, i.e., representing the total degree of polarization as [61] where are the vector and tensor polarization degrees, accordingly. Vector polarization is defined by three parameters t 10 , t 1−1 , t 11 and the tensor polarization is defined by five parameters t 20 , t 2±1 , t 2±2 , and the normalization Trρ is not considered as a parameter. The values of D v and D t are determined by the possible values of the components t 1Q and t 2Q , which can be determined from the condition of positivity property of the diagonal elements of the density matrix (ρ 11 ≥ 0, ρ 22 ≥ 0, ρ 33 ≥ 0), and the normalization condition for the total intensity.
Geometrically, each possible vector polarization state can be represented by a point P in a sphere with radius of √ 3/2 whose coordinates are t 10 , √ 2Ret 11 , and √ 2Imt 11 , respectively ( Figure 11).
Acceptable values of the degree of vector polarization are 0 where R v is the length of the vector OP ( Figure 11). If we determine the vector polarization in the direction of the z-axis, the components t 1±1 disappear. In the case of the direction of the polarization vector along the principal axis of the symmetric second-rank tensor formed by T 2Q , the components t 2±1 also disappear. The possible states of tensor polarization cannot be represented by a spherical area; these points form a cone inside a sphere with a radius of √ 2, and the permissible values of the degree of tensor polarization are 0 ≤ D t ≤ 1 (Figure 12).   The degree of polarization is determined by the length of the vector OP: D = R/ √ 2, R = OP. Thus, the degree of polarization D < 1 for a pure vector polarized beam (D t = 0), i.e., the vector polarized beam is only partially polarized ( Figure 11). This indicates that it cannot be represented by a pure polarization state. However, there is a completely polarized beam that does not have a preferred direction of the polarization vector (t 11 = t 1−1 = t 10 = 0), and this light beam is a pure tensor polarized beam (D = D t = 1). This corresponds to the point t 20 = − √ 2 at the top of the cone (point A in Figure 12). As can be seen from Figure 12, there are completely polarized states consisting of vector and tensor components of polarization on the circle abcd of the cone basis (for example, points a and c with t 10 = ± √ 3/2 and t 20 = √ 1/2). The points on the circle abcd correspond to circularly and linearly polarized beams and elliptically polarized states.
The unpolarized state of the beam is located in the center of the sphere. Note that in the 2D case, an unpolarized beam can be represented as the sum of two orthogonally (linearly or circularly) polarized beams, or as an incoherent mixture of equal amounts of two oppositely polarized beams. In the 3D case, an incoherent mixture of opposite polarization states χ 11 = (1, 0, 0) T and χ 1−1 = (0, 0, 1) T does not give a density matrix equal to the density matrix for an unpolarized beam [61]. The degree of polarization of this mixture is equal to D = 1 2 |t 20 | 2 1/2 = 1 2 . Moreover, the degrees of vector polarization for the pure state χ 0 = (0, 1, 0) T and the incoherent mixture of beams in the states χ 11 and χ 1−1 are equal to zero [61]. This means that vector polarization is not sufficient to completely determine the polarization state of a 3D light beam. Additional parameters (tensor polarization components) should be taken into account.
In the classical domain, there are nonzero intrinsic fluctuations of the Stokes parameters ∆σ 2 i = σ 2 i − σ i 2 associated with tensor polarization components. Here, It can be shown that two parameters corresponding to the vector and tensor degrees of polarization are sufficient for a complete description of the characteristics of the electromagnetic field, i.e., to distinguish between the two fields with different polarization matrices. Let us consider two diagonalized polarization matrices ρ and ρ with different eigenvalues (λ 1 , λ 2 , λ 3 ) and ( λ 1 , λ 2 , λ 3 ), respectively, having the same degrees of polarization D = D. It is easy to demonstrate that D v = D v and D t = D t if only the matrices ρ and ρ are identical, i.e., λ 1 = λ 1 , λ 2 = λ 2 , and λ 3 = λ 3 . Indeed, the degrees of polarization for these fields are determined by the expressions [61]: The condition of equality of the vector and tensor polarizations (D =D, D v =D v , D t =D t ) leads to the equalities λ 1 − λ 3 =λ 1 −λ 3 and λ 1 − λ 2 + λ 3 =λ 1 −λ 2 +λ 3 . It can be seen that these equations are valid only in the case of λ 1 =λ 1 , λ 2 =λ 2 , and λ 3 =λ 3 , i.e., two different density matrices cannot obey the condition D v =D v , and D t =D t , simultaneously. Thus, an unambiguous description of the polarization state of a light beam can be carried out using two parameters corresponding to the vector and tensor polarization degrees. Indeed, for the polarization matrices considered in [85], we have different vector and tensor polarization degrees for two different fields. For example, two density matrices with the eigenvalues λ 1 = 0.7, λ 2 = 0.2, λ 3 = 0.1 andλ 1 = 0.65,λ 2 = 0.343,λ 3 = 0.007 have the same values of the total degree of polarization D = 0.56. However, these matrices have different degrees of vector and tensor polarization: D v = 0.52, D t = 0.2, andD v = 0.56, D t = 0.0145, respectively.
To obtain a simple physical interpretation of the tensor polarization, it is useful to use the analogy between the polarization tensors and the multipole moments of the electrostatic potential generated by the distributed charge. It can be noted that the decomposition into irreducible representations of the rotation group SO(3) used above is actually a multipole decomposition [91]. The total intensity corresponds to an electric monopole, the vector polarization to an electric dipole, and the tensor polarization components to an electric quadrupole. However, there are limitations to the analogy with electric moments. If the electric dipole is a polar vector, then the polarization vector is an axial vector. Geometrically, tensor polarization can be represented by an ellipsoid of charge or an ellipsoid surface. In the case where one of the expectation values σ 2 x , σ 2 y , and σ 2 z disappears, the polarization tensor is represented by a degenerate ellipsoid, namely, a cylinder.
Note that the polarization ellipsoid associated with the polarization tensor is also analogous to the moment of inertia ellipsoid. The geometrical interpretation of the 3D polarization matrix in terms of the inertia ellipsoid and the orientation vector was considered in [88].
Since all the polarization components are associated with field correlations, the vector and tensor degrees of coherence corresponding to the vector and tensor parts of the density matrix decomposition can be introduced. In the far zone E z → 0, and all tensor components of the second-rank are transformed into vector components. This means that the 2D formalism can be applied to describe the field in the far zone.
Thus, to describe an arbitrarily polarized light beam, the vector and tensor polarization degrees can be introduced. The representation in terms of vector and tensor operators allows us to obtain a simple physical insight into the proposed parameters for the polarization of an arbitrary field. It is shown that the tensor polarization components can be associated with the intrinsic fluctuations of the Stokes parameters. It is shown that for a complete description of three-dimensional polarization, two measures of the degrees of polarization (vector and tensor) are sufficient, i.e., to distinguish two different polarized fields. The proposed definition of polarization can be useful for describing electromagnetic fields in near field optics and analyzing propagation of a focused light beam.

Discussion
The polarization of electromagnetic waves play an important role in the interaction of light and matter, optical transmission, photonics, nanoplasmonics, high-power lasers and amplifiers, optical coherence tomography and imaging.
Multimode fibers are becoming important for improving the capabilities of nextgeneration telecommunication systems due to the spatial division multiplexing technique [92]. In recent years, they have attracted considerable attention in optical communication [93,94], soliton dynamics [95,96], ultrafast photonics [97,98], and other areas. In [99], the formation and manipulation of breathing solitons in graded-index multi-mode fibers was numerically demonstrated. Multimode fibers and nonlinear propagation within them can also be used to improve the efficiency of all optical data processing devices [100].
It was shown in [70] that classical beam propagation in optical waveguides can simulate quantum phenomena such as squeezing, Schrodinger cat states, collapse and revival. The possibility of measuring extremely small displacements, which can be used in laser gravitational-wave experiments and in ultra-high resolution spectroscopy, is shown. The results obtained can be useful in various fields, such as gravitational wave experiments, high-precision spectroscopy, and quantum information processing and metrology.
Due to the high mode dispersion, it is important to continuously analyze and monitor the dynamics of the modes propagating in multimode fibers. In [101], the complete polarization control of coherent light transmitted through an MMF with a strong mode and polarization coupling is experimentally demonstrated. It is necessary to propose a new measure of the degree of depolarization for simultaneous monitoring of the evolution of polarization. Recently, a new definition of the depolarization degree of partially polarized beams-the dispersion degree of polarization-has been proposed to measure the performance of fiber depolarizers [102].
It is of interest to determine the length of depolarization of light in real optical fibers. It is known [8,9,12] that when coherent light propagates in optical fibers, significant changes in the state and degree of polarization occur. As can be seen from Formulas (30) and (37), the degree of polarization depends on the gradient parameter ω of the fiber, the axial displacement x 0 of the beam and the wavelength λ. From Equations (30) and (37), it is possible to determine that the depolarization length corresponds to a 50% decrease in the degree of polarization: Substituting here the values for the gradient parameter ω = 7 · 10 −3 µm −1 (a fiber with a radius r 0 = x 0 = 25 µm and a relative difference in the refractive index ∆ ≈ 6.9 · 10 −3 ) and the wavelength λ = 0.63 µm, we get the resulting distance l d ≈ 8.3 cm. Such depolarization was experimentally observed in an isotropic multimode optical fiber [12]. Depolarization is enhanced by increasing the axial displacement of the incident beam or the angle of its inclination to the waveguide axis. Since rays with a large axial displacement x 0 correspond to high-order fiber modes, the depolarization of modes with large numbers will be stronger compared to low-order modes. Therefore, in coherent communication systems, it is better to use single-mode fibers without birefringence, where the degree of polarization is retrieved with a period z = π/ω on propagation. This length in graded-index optical fibers is equal to z ≈ 0.5 mm.
Note that the depolarization of light has a diffraction origin and the consideration of the Rytov rotation as a "geometric" depolarization mechanism is not correct. The Rytov rotation leads only to the rotation of the plane of polarization and not to a change in the degree of polarization. In addition, the angle of rotation of the polarization plane does not decrease with decreasing wavelength, whereas depolarization disappears at λ → 0. As shown above, depolarization also occurs for meridional rays for which there is no Rytov rotation.
Physically, depolarization corresponds to a decrease in the degree of correlation between different components of the field or to the appearance of a field component that is not correlated with the original component. In real situations, depolarization is most often caused by the scattering of light on small-scale inhomogeneities of the medium. The physical cause of depolarization in our case is the scattering of light on large-scale inhomogeneities of the medium. Since the influence of large-scale inhomogeneities on the polarization characteristics is small, it is usually neglected. However, in waveguide problems, small corrections accumulate with distance and can lead to noticeable effects.
Thus, the decrease in the degree of polarization in multimode optical fibers has a diffraction origin and can only be explained using the wave approach. Since depolarization disappears at both λ → 0 and x 0 → 0, it can be interpreted as the result of the interaction of polarization (spin) and trajectory (orbital angular momentum). The degree of polarization of both linearly and circularly polarized beams decreases with distance according to the quadratic law. The depolarization of the meridional rays is less than that of the sagittal rays. For right-and left-handed circularly polarized light, the effect of asymmetry with respect to the sign of the twist of the sagittal ray trajectory is observed. Depolarization is enhanced by increasing the axial displacement of the beam, the gradient parameter of the fiber and the radiation wavelength.
For three-dimensional polarization, the coherence matrix of three-dimensional electromagnetic fields is represented as spherical or rectangular tensors, similar to the Gell-Mann matrices. The evolution of the coherency matrix can be defined in terms of the Hamiltonian obtained from the three-component field Maxwell equations. It is shown that two measures of the degrees of polarization (vector and tensor) are sufficient for a complete description of the three-dimensional polarization, that is, for distinguishing two fields with different polarization matrices [61]. The paper [59] shows the possibility of an unambiguous description of three-dimensional polarization using two measures, namely, the vector and tensor degrees of polarization. Note that the polarization matrix of any three-dimensional field can be directly measured using nanoprobe scatterers [103,104], which enables an experimental determination of the vector and tensor degrees of polarization of the light under consideration.
Finally, we present the fundamental effect of spin symmetry breaking via spin-orbit and tensor interactions, which occurs even in rotationally symmetric fibers. Due to this effect, the conversion of the incident linear polarized light into a particularly left-or righthanded circularly polarized light can be achieved. Although the effect manifests itself in a graded-index fiber, it manifests itself in general cylindrical potentials, in particular the Coulomb cylindrical potential. The phenomenon can be used in polarization control applications, such as new circular polarizers, polarization transformers and modulators.
Recently, new developments in the phenomenon known as transverse spin, which occurs in laterally confined electromagnetic waves, have been of great interest [105]. The spin-pulse blocking effect associated with transverse spin has already found promising applications for efficient spin-directed couplers [106,107].
More recently, interesting optical phenomena related to the angular momentum of light have been demonstrated theoretically and experimentally [108][109][110]. In [109], vectorial optomechanical effects mediated by the anisotropy of a nematic liquid crystal material and the longitudinal field component were discovered. It was shown in [108] that transverse spin is essentially a polarization-independent phenomenon that persists even in fields generated by completely unpolarized sources. The transverse spin appears in any paraxialto-nonparaxial transformation, even without a change in the degree of polarization. This result establishes an important connection between three-dimensional polarization in the transverse spin [103,111] and non-paraxial fields [58][59][60][61]112,113]. The origin of this phenomenon lies in the intrinsic spin-orbit interaction of light. These results provide a new opportunity to use polarization-independent spin from unpolarized sources, such as in nanoelectronics, telecommunications, photovoltaics, optical imaging, etc. Transverse spin can also be useful for unidirectional coupling of light in and from optical fibers or waveguides, which may have implications for future quantum networks [114,115].

Conclusions
In this review, we focused on studying the effects of spin-orbit interactions in the propagation of vortex light through multimode graded-index fibers. Several detailed examples discussed in this review demonstrate that optical fibers supporting vortex beams can be considered using an operator approach well-developed in quantum mechanics. With this approach, all the dynamics of the system are transferred to the operators. This makes it possible to study the evolution of the beam parameters using purely algebraic procedures, i.e., without using explicit expressions for the wave functions of the field and without calculating the corresponding integrals. This model is simple and provides a physical clarification of the phenomenon. Note that the operator approach can be used to study the evolution of polarization in optical fibers with birefringence, absorption, or amplification, as well as with random inhomogeneities.
Understanding the relationship between the degree of polarization, the spin-orbit interaction and the Berry phase is important for the design of optical devices, such as new polarizers, polarization transformers, and modulators. Fluctuations of the Berry's phase, which have a wave nature, during the propagation of light in an optical fiber were found in [57]. It was shown that the variance of the rotation angle of the polarization plane increases with distance and can be determined by measuring the degree of polarization.
Spin-dependent effects on vortex light beams propagating in an optical fiber are demonstrated by solving the full three-component field Maxwell equations using perturbation analysis. The splitting of degenerate modes due to the spin-orbit interaction is shown. Vector and tensor polarization degrees are introduced to characterize the polarized beam. It is shown that the degree of the pure vector polarized light beam is less than 1, that is, the vector polarized light beam is only partially polarized. The geometric interpretation of vector and tensor polarizations is presented.