1. Introduction
A phase singularity nested in the cross section of the coherent beam can form a vortex beam with a spiral wavefront [
1], and its Poynting vector propagates along the spiral [
2]. The typical feature of vortex beams is that the phase factor
is ubiquitous in the phase structure. On average, each photon in a vortex beam carries the orbital angular momentum (OAM) of
[
3], where
is the quantum number (also known as topological charge), which is generally an integer or a fraction [
4]. The rapid development of OAMbased coding and multiplexing technology in recent years has made OAM a new degree of freedom for encoding information [
5,
6].
Due to the application value of the vortex beam as an information carrier in improving communication capacity and spectral efficiency [
7], it has become a research hotspot in the field of free space optical communication. However, influenced by turbulent effects including beam wander [
8], scintillation [
9], OAM modes crosstalk [
10], and so on, and due to the divergence of the optical vortex beams themselves [
11], the transmitted information will randomly fluctuate at the communication system receiver. In order to mitigate the divergence of vortex beams, a series of structured beams with non-diffracting and self-healing properties in recent years have been widely used as the background to nest vortices [
12,
13,
14]. Among them, autofocusing Airy beams managed to catch attention due to their unique focusing characteristics, which can effectively weaken the beam spread effect caused by atmospheric turbulence [
15]. In recent years, the propagation characteristics of autofocusing Airy vortex beams (AAVBs) in free space have been extensively studied [
16,
17,
18]. Since not being satisfied with the turbulence-resistant capability inherent in AAVBs, researchers have further regulated the phase structures to obtain better transmission features. For example, AAVBs with nontraditional spiral phases, represented by the power exponential spiral phase [
19,
20,
21], have been successively proposed [
22,
23,
24]. Although the focusing ability of these new vortex beams is significantly enhanced compared with that of the traditional AAVBs, the vortex splitting phenomenon occurs obviously even when the beams transmit in homogeneous medium. According to the research, vortex splitting will lead to the OAM transmit–receive mismatch [
25] in the communication system, resulting in serious distortions of the information transmitted. Therefore, the aim of this paper is to study the way to effectively reduce the vortex splitting in free space when strengthening the focusing ability of AAVBs. In addition, the advantages of the polarization state as the degree of freedom of the light field are not considered in the above literature. More and more studies have shown that the non-uniform modulation of the polarization states of vortex beams can control the energy distribution and spin-orbit angular momentum coupling of light fields in free space propagation [
26].
Through further study of the above nontraditional spiral phase regulation mechanism, we found that the reason for the deepening of the vortex splitting phenomenon is that the rotational symmetry of the traditional spiral phase is broken. In this paper, the phase of the AAVBs with radial-variant polarization is modulated by the Pancharatnam- Berry phase that is closely related to the variation of polarization state, and a new form of nontraditional spiral phase is proposed. The characteristic of this nontraditional spiral phase is that the wavefront isophase line is curved (here, the tailored wavefront isophase line is curved in source plane, not straight in tradition) while maintaining the rotational symmetry of the traditional spiral phase. Based on the generalized Huygens Fresnel principle and multilayer random phase screen theory, a model of radially polarized AAVBs propagation in atmospheric turbulence is established. On this basis, it is proven that the AAVBs with radial polarization variation are more advantageous than the traditional scalar AAVBs in restraining vortex splitting caused by atmospheric turbulence.
2. Methods
The electric field intensity of the AAVBs with radial-variant polarization in the source plane [
27] can be rewritten as
where
is the topological charge,
is the polar coordinate,
represents Airy function,
is the radius of the main ring,
is the cross-section scale factor,
is the apodization factor,
and
are the radial polarization indices, and
is the parameter regulating the overall polarization states of the light field. It should be noted that Equation (1) in this paper differs from that in the literature [
27] in the following three points: first, the spiral phase factor
and the radial polarization regulation index
are added in this paper; second,
is the radius of the light field, not the radius of the main ring in the literature [
27]; third,
in Equation (1) in this paper is a parameter describing the overall polarization state of the light field, while
in the literature [
27] represents the polarization state at
.
By solving Equation (1), the influence of
,
, and
on the radial distribution of the polarization states of the autofocusing Airy beams can be obtained, as shown in
Figure 1. The parameters used for calculation were set as:
, wavelength
, main ring radius
, cross-section scale factor
, apodization factor
, and light field radius
. In
Figure 1, the first row is the polarization state distribu- tions of the light field under different radial polarization indices. The second row is a sim- plified schematic to clearly describe the regularity of polarization state change shown in the first row. The third and fourth rows are the horizontal and vertical components of the polarization distribution shown in the first row, respectively.
Figure 1(a1–a4) shows the polarization distribution of the scalar autofocusing Airy beams, and
Figure 1(b1–b4, c1–c4, d1–d4, e1–e4) shows the polarization distributions of the vector autofocusing Airy beams. As shown in
Figure 1(b1,b2), the polarization states of the light field rotate clockwise along the radius direction from the center of the light field. By comparing
Figure 1(b2) with
Figure 1(c2), it can be found that the polarization states of the light field rotate as a whole with the change of
. By comparing
Figure 1(d2) with
Figure 1(e2), it can be found that, with the increase of
, the radial change frequency of the beam polarization state is obviously accelerated. The reasons are as follows: as can be seen from the expression
, when
,
,
increases by
(one period) in the process of
increasing from 0 to
. As shown in
Figure 1(d2), the light field is vertically polarized at
and deflected with the increase of
. However, when
increases to
, the polarization state of the light field returns to the vertical polarization state, that is, a period of change is completed. With the further increase of
, a period of change is completed as
increases from 0 to
. Therefore, the radial polarization state change frequency of the light field increases with the increase of the radial polarization index
, as shown in
Figure 1(d3,e3,d4,e4). By comparing
Figure 1(d1–d4) with
Figure 1(b1–b4), it can be found that with the increase of
, the radial change frequency of the beam polarization state becomes smaller, because
in this paper, when
,
,
cannot change by
as
increases from
to
. The phenomenon that the radial change frequency becomes smaller caused by the increase of
can also be clearly illustrated by comparing
Figure 1(c3) with
Figure 1(d3) and
Figure 1(c4) with
Figure 1(d4). On the contrary, when
,
, the radial polarization state change frequency of the light field increases with the increase of
. In conclusion, both the radial polarization indices
and
can be used as the control parameters to modulate the distribution of the radial polarization states.
Each polarization state of light can be presented as a point on the surface of the high-order Poincare sphere [
28].
Figure 2a,b are the phase and polarization state for an arbitrary point A on the high-order Poincare sphere. In the seminal article [
29], a quantum system in the eigenstate
returned to the original state after a cyclic adiabatic transformation, and the final wave function
acquired a geometric phase factor depending only on the geometry of the path
. In optics, according to the Pancharatnam-Berry phase theory [
30], as the light’s polarization state undergoes a cyclic transformation, it acquires a corresponding Pancharatnam-Berry phase
purely depending on the closed path (e.g., A–B–C–A) in the parameter space. For example, under the action of phase
, the spiral phase shown in
Figure 2a will be transformed into that illustrated in
Figure 2c.
Therefore, when the scalar light field shown in
Figure 1(a1) changes into
Figure 1(b1–e1) by regulating the distribution of the radial polarization state, different Pancharatnam-Berry phases will be generated. When
, the Pancharatnam-Berry phases generated by radial polarization regulation will affect the spatial structure and transmission characteristics of spiral phase.
Figure 3 shows the phase distribution of the AAVBs in the source plane under different regulations of radial polarization states when
. It should be noted that the values of
and
in
Figure 3 were arbitrarily selected, and other parameters are the same as those in
Figure 1. In general, the wavefront isophase line of scalar vortex beams is straight in the source plane, as shown in
Figure 3a, and the wavefront isophase line will be curved in propagation. However, it can be seen from
Figure 3(b–f) that the wavefront isophase line became curved in the source plane caused by the Pancharatnam-Berry phases generated by the radial polarization regulation, and its curvature change is closely related to the values of the radial polarization indices
and
. According to
Figure 3a,b,e, the curvature of the wavefront isophase line increases with the increase of
. According to
Figure 3b,c,f, the curvature of the wavefront isophase line decreases with the increase of
. The regularity of the wavefront isophase line curvature changing with
and
corresponded to that of the polarization state frequency varying with
and
.
According to the vector angular spectrum theory, the electric field of the vector AAVBs shown in Equation (1) after transmitting a certain distance along the
z-axis can be expressed in the Cartesian coordinate system as [
27]:
where
and
are the angular spectra of the light field distribution in the source plane, which can be obtained by the Fourier transform of
and
after transforming them in the Cartesian coordinate system;
and
are the wavenumbers; and
and
are the spatial frequencies in the
x-axis and
y-axis directions. According to Equation (2), the intensity distributions of vector AAVBs in the focal plane can be obtained, as shown in
Figure 4.
In view of the tiny distinction between the wavefront curvature in the main ring region of the light fields shown in
Figure 3c,d and that in
Figure 3a, which both approximated straight lines, and in order to compare the influences of the wavefront isophase line curvature on the focal field distribution,
Figure 4a–d, respectively, gives the normalized distribution of the light field intensity on the focal planes to which the light fields shown in
Figure 3a,b,e,f transmit.
Figure 4e is the intensity distribution curves of the light field cross-section, as shown in
Figure 4a–d. From
Figure 4e, it can be clearly seen that, compared with the traditional scalar AAVBs, the focal field intensity of AAVBs with radial polarization change is significantly increased. The focal field intensities were
,
,
,
,
(in order of intensity). In addition, the positions of dotted lines with different colors on the x-axis in
Figure 4e represent the focal spot main ring radii of the AAVBs under the corresponding wavefront curvatures. It can be seen that the main ring radii of the focal field under different wavefront curvatures were
,
,
,
,
(in order of radius size). The reason for this phenomenon is that when light wave propagates, its phase represents its local direction, and the normal direction of the isophase surface (i.e., wavefront) is consistent with the energy flow direction of the light wave. As the parameters
n and
m change, the curvature of the wavefront isophase line and phase distribution change, and then the intensities and main ring radii in the focal plane change accordingly.
As illustrated in
Figure 5, when implementing this scheme in practice, similar to the prior demonstration [
22], we can generate AAVBs by using an expanded Gaussian beam to illuminate a common programmable spatial light modulator (SLM) on the surface of which the specific phase hologram (e.g.,
Figure 3) is displayed. In our scheme, the Pancharatnam-Berry phase induced by the variation of polarization states of AAVBs can be modulated by controlling the parameters
and
, with no need for complex optical devices such as a metasurface [
31].
When the proposed AAVBs propagate in atmosphere, according to the generalized Huygens-Fresnel principle, the disturbance of atmospheric turbulence on the beam can be equivalent to a complex phase, which directly acts on the Huygens Fresnel paraxial expression of the light field propagating in vacuum. When the fluctuation of atmospheric refraction index is slight, the change of the light field caused by the inhomogeneity of the atmospheric refraction index is thought to be reflected only in the phase. Therefore, the process of beam propagation in atmospheric turbulence can be decomposed into two processes: one is that the beam propagates in a vacuum; the other is that the phase is disturbed by the atmospheric refraction index fluctuation. These two independent and simultaneous processes are multiple phase screen models of atmospheric turbulence, which can be expressed mathematically as [
32]:
where
is the light field distribution in the source plane,
or
,
represents the
j-th phase screen, and parameters such as
,
,
,
,
,
,
,
, and
are the same as those in reference [
32]. According to the Kolmogorov turbulence theory, although the turbulence is random and non-isotropic on the whole, it can be considered as homogeneous and isotropic in a given local area. This is the basic theory for studying optical propagation in turbulent atmosphere. Based on this theory, all components of the vector beams will also follow the wave equation in the same scalar form [
33]. Therefore, all components of the vectorial AAVBs after passing through atmospheric turbulence can be obtained by substituting Equation (2), respectively, into Equation (3), and then the final light field distribution can be obtained by vector synthesis.
4. Discussion
The reason for the phenomena shown in
Figure 6 and
Figure 7 is that, in the weak turbulence, the intensity of the individual vortex wandering in the cross-section of the light field is lower, allowing for the weak main ring energy to be adequate enough to bound the individual vortex. At this moment, the larger focal main ring radius will make the wandering effect caused by atmospheric turbulence relatively small [
32]. When the turbulence intensity increases, the weak intensity of the focal main ring cannot confine the wandering individual vortices to the area of focal main ring enclosed. Therefore, as shown in
Figure 4e and
Table 1, the stronger the intensity of the focal main ring is, the smaller the vortex splitting rate is.
Although it seems that the amount of reduction of splitting ratio achieved is a little shown in
Table 1, it should be made clear that the main aim in this paper was that a new spiral phase structure is proposed to control the focusing property of AAVBs. The most salient feature of this new spiral phase is that its wavefront isophase line is curved. Through the numerical experiment, we can find that the focusing property of AAVBs is not a monotonic function of the curvature of the isophase line. It is therefore reasonable to infer that there must be optimal wavefront curvatures for different turbulences. As long as we can find the optimal wavefront curvatures, the amount of reduction of splitting ratio could be enlarged enough to satisfy the typical applicative scenarios. This paper provides valuable insights for the study of vector vortex beam dynamics.
It should be noted that, due to complexity of the function of autofocusing Airy vortex beam with the radial polarization indices, it is very difficult to obtain analytic solutions by solving the paraxial wave equation. Numerical computations thus have to play the role of an efficient and important tool to quantitatively investigate the propagation effects [
35]. With the help of the split-step propagation method and angular-spectrum theory,
Figure 8 shows the trajectories of AAVBs with
m = 0, 0.5, 1, and 2. On the one hand, in comparison with the propagation distance (z = 3.5 km), the size of AAVBs (
,
) is very small, which satisfies the paraxial approximation condition. On the other hand, according to the research [
36], in the nonparaxial regime, Airy beams can bend at large angles. As illustrated in the figure, comparing the trajectory of AAVB (
m = 0) with that of AAVBs (
m = 0.5, 1 and 2), the introduction of
m does not make the trajectories of AAVBs bend at large angles. Therefore, Equation (1) with parameter
m satisfies the paraxial wave equation.