Bessel–Gauss photon beams with fractional order vortex propagation in weak non-Kolmogorov turbulence

Jie Gao; Yu Zhu; Donglin Wang; Yixin Zhang; Zhengda Hu; Mingjian Cheng

doi:10.1364/prj.4.000030

Abstract

We model the effects of weak fluctuations on the probability densities and normalized powers of vortex models for the Bessel–Gauss photon beam with fractional topological charge in the paraxial non-Kolmogorov turbulence channel. We find that probability density of signal vortex models is a function of deviation from the center of the photon beam, and the farther away from the beam center it is, the smaller the probability density is. For fractional topological charge, the average probability densities of signal/crosstalk vortex modes oscillate along the beam radius except the half-integer order. As the beam waist of the photon source grows, the average probability density of signal and crosstalk vortex modes grow together. Moreover, the peak of the average probability density of crosstalk vortex modes shifts outward from the beam center as the beam waist gets larger. The results also show that the smaller index of non-Kolmogorov turbulence and the smaller generalized refractive-index structure parameter may lead to the higher average probability densities of signal vortex modes and lower average probability densities of crosstalk vortex modes. Lower-coherence radius or beam waist can give rise to less reduction of the normalized powers of the signal vortex modes, which is opposite to the normalized powers of crosstalk vortex modes.Physics (Grant No. 11447174), and the Fundamental Research Funds for the Central Universities (JUSRP51517).

Keywords

Atmospheric turbulence Photon statistics Quantum communications

1. INTRODUCTION

In recent years, the study of the stable vortex beams whose vortex model per photon can take an arbitrary value within a continuous range, either integer or noninteger, in units of $ℏ$ is a newly burgeoning field ^[1–5]. The propagation property of fractional vortice-imprinted Bessel beams in free space has been given attention both theoretically and experimentally ^[6–8]. It has been verified that vortices with fractional strength do not propagate unaffected in free space. However, neglecting their evolving intensity and vortex structure, a vortex integrated over the whole beam cross section is invariant under free-space propagation ^[7]. The fractional Bessel beams (nondiffracting vortex beams) are of special interest due to their properties of divergence-free propagation and self-repair, and have generated widespread interest in the last decade for applications in optical communications ^[9–13]. Götte proposed a quantum mechanical description of the beams with fractional topological charges ^[13], and a high-order Bessel nonvortex beam of fractional type (HOBNVBs-F) was introduced by Mitri ^[10–12]. Next, we developed the model of average probability densities and the normalized powers of the signal/crosstalk orbital angular momentum modes for the fractional-order Bessel–Gauss (FoBG) photon beams in the turbulent atmosphere of strong irradiance fluctuations ^[14]. On the other hand, we know that Kolmogorov’s power spectral density model is widely used and accepted to describe wave propagation through atmospheric turbulence, but the experimental data in the last several decades shows that atmospheric turbulence generally possesses a structure different from the Kolmogorov model ^[8]. The Kolmogorov turbulence is actually an exception to non-Kolmogorov turbulence. When describing the effect of turbulence on fractional Bessel–Gauss photon beams in different conditions, it is better using non-Kolmogorov turbulence than Kolmogorov turbulence, which is especially embodied in its exponent value. As we know, there is almost no discussion with respect to the effects of non-Kolmogorov turbulence on the probability densities and the normalized powers of the vortex modes for FoBG photon beams in the weak and paraxial atmosphere channel.

In this paper, we analyze the influence of non-Kolmogorov turbulence on the vortex mode probability densities and the normalized powers for FoBG beams along the direction of the beam radius. Our aim is to exploit this changing process to present the effects of the weak non-Kolmogorov turbulence on Bessel–Gauss photon beams. We anticipate that this work will be particularly useful when applied to practical application for optical wireless communication. In Section 2, the models of the probability density for vortex modes and the normalized power of vortex per photon per unit length of a transverse slice of the beam for FoBG beams in non-Kolmogorov turbulence are established. The effects of fractional topological charges, turbulence strength (in weak fluctuation region), propagation distance, the non-Kolmogorov turbulence parameter, and wavelength on the probability densities of signal and crosstalk vortex modes of FoBG photon beams along the direction of the receiving plane are researched in Section 3. Conclusions are presented in Section 4.

2. MODE PROBABILITY DENSITY

In the weak fluctuation region ^[15] and in the half-space $z > 0$ , the normalized complex amplitude of FoBG beams in a cylindrical coordinate system $(r, φ, z)$ can be expressed as ^[16] where $r = | r |$ , $r = (x, y)$ is the two-dimensional position vector in the source plane; $z$ is propagation distance; $ψ (r, φ, z)$ is the complex phase of waves propagating through turbulence; $l_{0}$ describes the topological charge of the helical structure of the wave front around a wave front singularity; and ${FoBG}_{free} (r, φ, z)$ is the normalized complex amplitude of the FoBG beam at the $z$ plane in the free turbulence channel. In the paraxial region, the normalized complex amplitude ${FoBG}_{free} (r, φ, z)$ of a nondiffracting Bessel–Gauss photon beam with fractional topological charge has the form ^[7,15] where $μ (z) = 1 + i z / z_{R}$ , $z_{R} = k w_{0}^{2} / 2$ , $k = 2 π / λ$ is the wavenumber; $w_{0}$ is the beam waist at the $z = 0$ plane; $λ$ is the wavelength; $z_{R}$ is the Rayleigh range; $J_{| l_{0} |}$ is the Bessel function of integer orders; $k_{r} = ξ k$ is the transverse wavenumber; and $γ$ is the fractional index of any real number ^[7]. To keep it simple, we use $C_{γ}$ to denote the term $\frac{{(\pm i)}^{γ} \sin (\pm π γ)}{π} \cdot \frac{i^{| l_{0} |}}{\pm γ - l_{0}}$ .

Sign up for Photonics Research TOC. Get the latest issue of Photonics Research delivered right to you！Sign up now

As we know, from the view of quantum theory, the atmospheric turbulence fluctuations disturb the complex amplitude of the propagating photon beam, which is no longer guaranteed to be in the original eigenstate of topological charge. The resulting photon beam now can be regarded as a superposition of the waves with new topological charge. Similar to the discussion of ^[7,17] for beams with integer topological charge, the FoBG is expanded as an integral with orthogonal basis $\exp (i l φ)$ which carries topological charge $l$ ; then the complex amplitude $FoBG (r, φ, z)$ of the FoBG beam can be written as where the term $β_{l} (r, z)$ is the mode amplitude of the vortex mode with topological charge $l$ at the position $(r, z)$ given by

As we all know, the mode probability density is equal to ${| β_{l} (r, z) |}^{2}$ . In turbulent media, the mode probability density is associated with the ensemble average over turbulent medium, i.e., ${〈 {| β_{l} (r, z) |}^{2} 〉}_{a t}$ , where $〈 \cdot 〉$ denotes the average over the ensemble of the turbulent atmosphere. By Eqs. (1) and (4) and proceeding with the ensemble average of atmospheric turbulence for $β_{l} (r, z) β_{l}^{*} (r, z)$ , where $*$ denotes complex conjugate, the ensemble averaging mode probability density $D_{l} (r, z)$ of vortex models $l$ of an FoBG beam in paraxial channel $(k_{r} \to \infty)$ is given by

In the far-field region ^[15], ${〈 \exp [ψ (r, φ, z) + ψ^{*} (r^{'}, φ^{'}, z)] 〉}_{a t}$ is given by where $Φ_{n} (\cdot)$ denotes the spatial power spectrum of atmospheric turbulence and $J_{0} (\cdot)$ is the Bessel function of the first kind with 0 order. For non-Kolmogorov turbulence, the spatial power spectrum is commonly given by the expression ^[18] where $α$ is the parameter of non-Kolmogorov turbulence, $Γ (α)$ denotes the Gamma function, and $A (α) = Γ (α - 1) \cos (α π / 2) / 4 π^{2}$ ; $κ_{0} = 2 π / L_{outer}$ and $κ_{m} = {2 π Γ [(5 - α) / 2] A (α) / 3}^{1 / (α - 5)} / L_{inner}$ , where $L_{outer}$ and $L_{inner}$ are the outer and inner scale of turbulence, respectively. $C_{n}^{2}$ is the generalized refractive-index structure parameter with units $m^{3 - α}$ ; for simplicity’s sake, we set $C_{n}^{2}$ for FoGB beams as a constant to discuss the influence of $α$ on the beam. To simplify the numerical discussion in Section 3, we set $L_{outer} = \infty$ without considering the effect of outer scale, and $L_{inner} = 1 mm$ ^[15].

In the paraxial approximation of the wave phase structure function, Eq. (6) can be expressed in the following form ^[15]: where $ρ_{0}$ is the spatial coherence radius of a spherical wave propagating in the weak non-Kolmogorov turbulence ^[18,19] and given by based on the integral expression ^[20] where $I_{n} (η)$ is the Bessel function of the second kind with $n$ order. Making use of Eqs. (2), (5), and (7) as well as the Bessel function orthogonal completeness we have the probability density of the vortex modes for FoBG beams in paraxial turbulence channel where $w (z) = w_{0} \sqrt{1 + z_{R}^{2} / z^{2}}$ is the Gaussian spot size and $Δ l = l - l_{0}$ denotes the value of crosstalk. For $Δ l = 0$ , $D_{l_{0}} (r, z)$ is the probability density of the signal vortex mode $l_{0}$ for FoBG beams, and for $Δ l \neq 0$ , $D_{l \neq l_{0}} (r, z)$ is the probability density of the crosstalk vortex mode $l = l_{0} + Δ l$ , which represents the probability density of the part of the energy launched into the signal vortex mode redistributed into other vortex modes by atmospheric turbulence ^[21]. The term $C_{γ}^{2}$ indicates the vortex model distribution for different fractional values $γ$ when $l_{0}$ changes (see Fig. 1).

$Distribution of vortex models with fractional topological charge γ when l0 changes. Every subplot has a peak at the nearest integer to γ. When γ=3.5 (half-integer), this results in two peaks of equal height at the two neighboring integers. The spread in the distribution is determined by the fractional value γ.$

Figure 1.Distribution of vortex models with fractional topological charge $γ$ when $l_{0}$ changes. Every subplot has a peak at the nearest integer to $γ$ . When $γ = 3.5$ (half-integer), this results in two peaks of equal height at the two neighboring integers. The spread in the distribution is determined by the fractional value $γ$ .

The normalized power of vortex per photon per unit length of a transverse slice of the beam is calculated with the expression ^[7] where $B$ is the is the beam width of the channel, $\sum_{l = - \infty}^{\infty} B P_{l}$ is the beam power, and $B P_{l}$ is the power of each vortex mode with topological charge $l$ with and given by ^[7]

Similarly, for $l = l_{0}$ , $L_{z} (l_{0})$ expresses the normalized power of the signal vortex mode for FoBG photon beams and, for $l \neq l_{0}$ , $L_{z} (l)$ is the normalized power of the crosstalk vortex modes.

3. NUMERICAL DISCUSSION

As we all know, the longer wavelength $λ$ is beneficial to the propagation of photon beams for less attenuation and lower crosstalk. Besides, it is entirely reasonable that the average probability density of signal vortex modes and signal normalized powers decrease, and the average crosstalk probability of vortex modes and crosstalk normalized powers increase when the propagation distance $z$ increases. Then, we do not discuss the influences owed to the wavelength $λ$ and the propagation distance $z$ . To simplify, in the following figures, we now define mode probability as average probability density ( $D_{l_{0}} (r, z)$ ) of the signal vortex modes, and crosstalk probability as the average probability density ( $D_{l \neq l_{0}} (r, z)$ ) of the crosstalk vortex modes.

First, we study the effects of the fractional order of topological charge $l$ , the fractional parameter $γ$ , the non-Kolmogorov turbulence coherence radius $ρ_{0}$ (in weak fluctuation region), and beam waist $w_{0}$ on the average probability densities $D_{l} (r, z)$ along the direction of $r$ for FoBG beams. The simulation results are shown in Figs. 2–5. In Fig. 2, we evaluate performance of average probability densities $D_{l} (r, z)$ along the direction of the beam radius $r$ for FoBG beam with $w_{0} = 0.1 m$ , $ξ = 0.01$ , and $Δ l = 0$ , 1, 2, 3, and 4. The parameters for the simulation of atmospheric turbulence channel are as follows: $z = 1 k m$ , $α = 11 / 3$ , and $C_{n}^{2} = 10^{- 15} m^{3 - α}$ . From Fig. 2, we can see that crosstalk probability $D_{l \neq l_{0}} (r, z)$ decreases with the increase of topological charge deviation $Δ l$ and it drops rapidly after the slow growth with the increase of $r$ . When $r$ is large enough (in Fig. 2, $r > 0.2 m$ ), mode probability $D_{l_{0}} (r, z)$ and crosstalk probability $D_{l \neq l_{0}} (r, z)$ tend to be the same. The results show that we can obtain the higher signal crosstalk probability ratio only in the circle region of radius $r < 0.2 m$ , and the center of the circle region is the beam center. So in the next context, we adopt the circle region of radius $r < 0.2 m$ .

Figure 2.Average mode probability $D_{l} (r, z)$ of FoBG beam along the direction of the beam radius $r$ for different values of $Δ l$ .

$Average probability densities Dl(r,z) of vortex modes for fractional FoBG beam along the direction of the beam radius r with different values of γ. (a) Δl=0, mode probability; (b) Δl=1, crosstalk probability.$

Figure 3.Average probability densities $D_{l} (r, z)$ of vortex modes for fractional FoBG beam along the direction of the beam radius $r$ with different values of $γ$ . (a) $Δ l = 0$ , mode probability; (b) $Δ l = 1$ , crosstalk probability.

$Average probability densities Dl(r,z) of vortex modes for fractional FoBG beam along the direction of the beam radius r with different values of α. (a) Δl=0, mode probability; (b) Δl=1, crosstalk probability.$

Figure 4.Average probability densities $D_{l} (r, z)$ of vortex modes for fractional FoBG beam along the direction of the beam radius $r$ with different values of $α$ . (a) $Δ l = 0$ , mode probability; (b) $Δ l = 1$ , crosstalk probability.

$Average probability densities Dl(r,z) of vortex modes for fractional FoBG beam along the direction of the beam radius r with different values of Cn2. (a) Δl=0, mode probability; (b) Δl=1, crosstalk probability.$

Figure 5.Average probability densities $D_{l} (r, z)$ of vortex modes for fractional FoBG beam along the direction of the beam radius $r$ with different values of $C_{n}^{2}$ . (a) $Δ l = 0$ , mode probability; (b) $Δ l = 1$ , crosstalk probability.

In Fig. 3, the average probability density $D_{l} (r, z)$ is plotted as a function of the beam radius $r$ for topological charge deviation $Δ l = 0$ and 1, with $w_{0} = 0.10 m$ , $ξ = 0.01$ , $z = 1 k m$ , $α = 11 / 3$ , and $C_{n}^{2} = 10^{- 15} m^{3 - α}$ . Figures 3(a) and 3(b) are for $γ = 3.1$ , 3.3, 3.5, 3.7, and 3.9, respectively. As is seen from Fig. 3, the $D_{l} (r, z)$ for $γ = 3.1$ , 3.3, 3.7, and 3.9 are the vibration curves which are around the curve of $γ = 3.5$ . Figure 3(a) indicates that mode probability decreases with the increase of $r$ , while the crosstalk probability decays after the increase in the area close to the circle area in Fig. 3(b). Figures 3(a) and 3(b) show that the curve is without vibration when the value of $γ$ is half-integer. And the more close the value of $γ$ is to half-integer, the closer its curve tends to the curve with $γ = 3.5$ . In order to facilitate the beam properties, we set $γ = 3.5$ .

Figures 4 and 5 present how the various parameters of atmospheric turbulence affect the average probability densities $D_{l} (r, z)$ . Figure 4(a) is a plot of the average mode probability $D_{l_{0}} (r, z)$ versus the varied beam radius $r$ for the different non-Kolmogorov turbulence parameters $α$ . The curve associated with $α = 3.67$ corresponds to Kolmogorov turbulence. The curves in Fig. 4(a) are consistent with ^[18] in that the smaller $α$ comes with the larger $D_{l_{0}} (r, z)$ in the receiving plane for any given propagation distance; the reason is that the spatial coherence radius $ρ_{0}$ tends to infinity when $α$ approaches 3. That is, the interruption of the turbulence on the photon beams is very weak, and a high $D_{l_{0}} (r, z)$ is acquired. However, when $α$ is close to 4, the wave aberration caused by turbulence is a pure wavefront tilt, which shifts the beam off-axis and results in the beam center being away from the receiving plane ^[18]. Thus, $D_{l_{0}} (r, z)$ decreases as $α$ approaches 4. The relationship between the average crosstalk probability $D_{l \neq l_{0}} (r, z)$ and the $α$ of non-Kolmogorov turbulence is revealed by Fig. 4(b). When the propagation distance $z$ is given, the $D_{l \neq l_{0}} (r, z)$ of FoBG photon beams increases rapidly when $α$ changes from 3.07 to 3.97 near the beam center regions (about $r < 0.04 m$ ). When detection position is far away from the beam center in the receiving plane, all the curves in Fig. 4(b) tend to be same.

Figure 5 is constructed to assess how the non-Kolmogorov turbulence coherence radius $C_{n}^{2}$ affects the average probability densities $D_{l} (r, z)$ with varied beam radius $r$ . The curves in Fig. 5(a) show that increasing $r$ makes mode probability $D_{l_{0}} (r, z)$ decrease sooner, while increasing $C_{n}^{2}$ doesn’t make the signal probability density change apparent. It is important to note that crosstalk probability $D_{l \neq l_{0}} (r, z)$ increases with the lower $C_{n}^{2}$ near the beam center (about $r < 0.04 m$ ), and the main energy still remains at the transmit signal channel. However, crosstalk probability decays and the energy of the adjacent channel decreases immediately as the detection position moves away from the beam center. For instance, all the curves slowly tend to overlap; namely, the effects of different $C_{n}^{2}$ value on crosstalk probability densities become less distinct when $r > 0.15 m$ .

Figure 6 investigates the average probability $D_{l} (r, z)$ as a function of beam radius $r$ for a FoBG beam with topological charge deviation $Δ l = 0$ and 1, $γ = 3.5$ , and $ξ = 0.01$ for $w_{0} = 0.02 m$ , 0.05 m, 0.08 m, and 0.10 m, respectively. The parameters of the atmospheric turbulence channel for the simulation are as follows: $z = 1 k m$ , $α = 11 / 3$ , and $C_{n}^{2} = 10^{- 15} m^{3 - α}$ . Figure 6(a) shows that increasing $w_{0}$ makes mode probability decay slowly. The results of Fig. 6(b) reveal that higher $w_{0}$ corresponds to the higher crosstalk probability densities and makes the curve peak away from the beam center.

$Average probability densities Dl(r,z) of vortex modes for fractional FoBG beam along the direction of the beam radius r with different values of beam waist w0. (a) Δl=0, mode probability; (b) Δl=1, crosstalk probability.$

Figure 6.Average probability densities $D_{l} (r, z)$ of vortex modes for fractional FoBG beam along the direction of the beam radius $r$ with different values of beam waist $w_{0}$ . (a) $Δ l = 0$ , mode probability; (b) $Δ l = 1$ , crosstalk probability.

For the sake of simplicity, we will discuss the impact of $ρ_{0}$ on the mode probability densities and normalized powers instead of the basic parameters $α$ and $C_{n}^{2}$ . The normalized powers $L_{z} (l)$ of signal and crosstalk vortex modes for the parameters $α = 11 / 3$ , $C_{n}^{2} = 10^{- 15} m^{3 - α}$ , $w_{0} = 0.10 m$ , and $Δ l = 0$ , 1, 2, 3, and 4 are shown in Fig. 7. Here we see that the crosstalk-normalized powers $L_{z} (l \neq l_{0})$ between vortex modes decrease with the increasing of $Δ l$ .

$Normalized powers Lz(l) of fractional vortex modes for FoBG beam along the direction of the transmission distance z with different values of Δl.$

Figure 7.Normalized powers $L_{z} (l)$ of fractional vortex modes for FoBG beam along the direction of the transmission distance $z$ with different values of $Δ l$ .

Finally, the effects of $ρ_{0}$ and $w_{0}$ on the normalized powers $L_{z} (l)$ of signal or crosstalk vortex modes per photon per unit length of a transverse slice for FoBG beams are shown in Fig. 8. To do that, we fix the parameters $γ = 3.5$ and $ξ = 0.01$ . Lower $ρ_{0}$ can give rise to increase of the normalized power of the signal vortex mode $L_{z} (l_{0})$ and decrease of the normalized powers of the crosstalk vortex modes $L_{z} (l \neq l_{0})$ . Obviously, greater $w_{0}$ can enhance the normalized power of the signal vortex mode and reduce the normalized powers of crosstalk vortex modes.

$Normalized powers Lz(l) of fractional vortex modes for FoBG beam along coherence radius r with different values of beam waist w0. (a) Δl=0, signal normalized powers; (b) Δl=1, crosstalk normalized powers.$

Figure 8.Normalized powers $L_{z} (l)$ of fractional vortex modes for FoBG beam along coherence radius $r$ with different values of beam waist $w_{0}$ . (a) $Δ l = 0$ , signal normalized powers; (b) $Δ l = 1$ , crosstalk normalized powers.

4. CONCLUSION

In conclusion, the model of probability densities and the normalized powers of signal or crosstalk vortex modes for FoBG photon beams in weak fluctuation turbulent atmosphere have been developed. Our results show that average mode probability $D_{l_{0}} (r, z)$ and average crosstalk probability $D_{l \neq l_{0}} (r, z)$ are the same when beam radius $r$ is large enough. For fractional order Bessel–Gauss beams, when $r$ varies, $D_{l} (r, z)$ oscillates except $γ$ is half-integer. For the half-integer beams, the increases of $α$ and $C_{n}^{2}$ result in the degradation of $D_{l_{0}} (r, z)$ , and average crosstalk probability $D_{l \neq l_{0}} (r, z)$ drops rapidly after the slow growth as $r$ increases, namely average crosstalk probability curves have peaks. Near the central region of the beam, enlarging $Δ l$ or minimizing $α$ , $C_{n}^{2}$ , and $w_{0}$ can make the $D_{l \neq l_{0}} (r, z)$ decrease. Furthermore, lower $ρ_{0}$ and $w_{0}$ can give rise to less reduction of the turbulence effect on the normalized power of signal vortex mode. Compared with signal normalized powers, the opposite changes in tendency of crosstalk-normalized powers are observed.

ACKNOWLEDGMENT

Acknowledgment. This work is supported by the Natural Science Foundation of Jiangsu Province of China (Grant No. BK20140128), the National Natural Science Foundation of Special Theoretical Physics (Grant No. 11447174), and the Fundamental Research Funds for the Central Universities (JUSRP51517).

References

[1] J. Leach, E. Yao, M. J. Padgett. Observation of the vortex structure of a non-integer vortex beam. New J. Phys., 6, 71(2004).

[2] W. N. Plick, M. Krenn, R. Fickler, S. Ramelow. Quantum orbital angular momentum of elliptically symmetric light. Phys. Rev. A, 87, 033806(2013).

[3] M. V. Berry. Optical vortices evolving from helicoidal integer and fractional phase steps. J. Opt. A, 6, 259-268(2004).

[4] J. C. Gutiérrez-Vega. Fractionalization of optical beams: II. Elegant Laguerre-Gaussian modes. Opt. Express, 15, 6300-6313(2007).

[5] J. B. Götte, K. O’Holleran, D. Preece, F. Flossmann, S. Franke-Arnold, S. M. Barnett. Light beams with fractional orbital angular momentum and their vortex structure. Opt. Express, 16, 993-1006(2008).

[6] S. H. Tao, X. Yuan. Self-reconstruction property of fractional Bessel beams. J. Opt. Soc. Am. A, 21, 1192-1197(2004).

[7] J. C. Gutiérrez-Vega, C. López-Mariscal. Nondiffracting vortex beams with continuous orbital angular momentum order dependence. J. Opt. A, 10, 015009(2008).

[8] E. Golbraikh, H. Branover, N. S. Kopeika, A. Zilberman. Non-Kolmogorov atmospheric turbulence and optical signal propagation. Nonlinear Processes Geophys., 13, 297-301(2006).

[9] J. C. Gutirrez-Vega, C. Lopez-Mariscal. Nondiffracting vortex beams with continuous orbital angular momentum order dependence. J. Opt. A, 10, 015009(2008).

[10] F. G. Mitri. Vector wave analysis of an electromagnetic high-order Bessel vortex beam of fractional type. Opt. Lett., 36, 606-608(2011).

[11] F. G. Mitri. High-order Bessel nonvortex beam of fractional type. Phys. Rev. A, 85, 025801(2012).