Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Self-recovery effect of orbital angular momentum mode of circular beam in weak non-Kolmogorov turbulence

Open Access Open Access

Abstract

It is generally true that the orbital angular momentum (OAM) mode persistently degenerate when a vortex beam propagates in the atmospheric turbulence. Here, however, we unveil an interesting self-recovery effect of OAM mode of the circular beam (CiB) in weak non-Kolmogorov turbulence. We show that the CiB displays the self-focusing effect and has clear focus in the weak non-Kolmogorov turbulence if we choose proper complex parameters, and the detection probability of the original OAM mode reaches the maximum at the focus. Our study proposes a method to alleviate the turbulent effects on OAM-based communication.

© 2016 Optical Society of America

1. Introduction

In general, the OAM-carrying beam has the helical phase front with regard to a phase term of exp(i), where l is an arbitrary integer and φ refers to azimuthal coordinate [1]. Different coaxial OAM eigen-modes of vortex beam are orthogonal, which has important application in free space communication [2,3]. While, it is well known that the effect of atmosphere turbulence is one of the key factors that limit the development of OAM-based communication, because the turbulence effects cause the spreading of the spiral spectrum [4–6]. Recently, the effect of Kolmogorov turbulence on single photon and entangled photon pairs has been studied in detail [7,8]. Moreover, the effect of non-Kolmogorov turbulence on some vortex beams, such as hypergeometric-Gaussian Beam [9], Hankel-Bessel beam [10] and Whittaker-Gaussian beam [11], has been discussed as well.

Very recently, some attention is paid to circular beam (CiB) [12,13], which is a very general solution of the paraxial wave equation. The equivalent form, special cases, normalization, Laguerre-Gaussian expansion, free-space divergence and parameter constraints of CiB have been discussed in previous works [12,13]. By assigning particular values to three complex beam parameters, the CiB is reduced to some well-known vortex beams such as standard Laguerre-Gaussian beam, elegant Laguerre-Gaussian beam, generalized Laguerre-Gaussian beam, hypergeometric beams, fractional-order elegant Laguerre-Gaussian beams, hypergeometric-Gaussian Beam and Whittaker-Gaussian beam [12].

The focusing property of beam can be employed in many fields, such as trapping and guiding microparticles [14] and realizing optical micromanipulation [15]. To our knowledge, currently, there is no further research to the influence of focusing property on the propagation of vortex beams in atmosphere turbulence. In this work, we discuss how the beam parameters affect the self-focusing property of CiB. Furthermore, we characterize the influence of non-Kolmogorov turbulence on the OAM modes of self-focusing CiB and reveal its self-recovery effect.

2. Power weight of OAM state for CiB in weak non-Kolmogorov turbulence

The normalized CiB propagating along z axis in cylindrical coordinates reads [13]

CiBp,l0(q0,q1)(r,φ,z)=(i2z0W0)|l0|+1[π|l0|!Ψp,l0(ξ)]121q(z)exp[ikr22q(z)]×[(1+ξ)q˜(z)q(z)]p2[rq(z)]|l0|F11(p2,|l0|+1;r2χ2(z)).×exp(il0φ)

In Eq. (1), there are three complex parameters q0, q1 and p, where q0 is given by q0 = iz0d0. W0 is a real number similar to the waist radii of Gaussian beam and hence z0=kW02/2 is similar to the Rayleigh range of wave number k. d0 is the position of waist of Gaussian envelope. l0 is an integer corresponding to the initial OAM quantum number. Besides, the other four complex parameters in Eq. (1) are defined as q(z) = q0 + z, (z) = q1 + z, ξ=(q1q0)/(q0*q1) and 1/χ2(z) = ik [1/q(z) − 1/(z)]/2. q(z) and (z) are similar to the q-parameter of Gaussian beam and 1/χ2(z) is the scale factor. Ψp,l0(ξ)=F12(p/2,p*/2,|l0|+1;|ξ|2) is the normalization factor. 1F1 and 2F1 denote the confluent hypergeometric function and hypergeometric function, respectively.

In weak atmospheric turbulence, the intensity fluctuation is very small that can be neglected, so we only consider the phase aberration on the complex amplitude of CiB. Supposing propagation distance z > 0, the complex amplitude of distorted CiB here can be expressed as

Ψp,q0(q0,q1)(r,φ,z)=CiBp,l0(q0,q1)(r,φ,z)exp[iψ(r,φ,z)],
where ψ(r, φ, z) is the term corresponding to complex phase perturbation caused by the turbulence. To elucidate the weight of all OAM modes, the complex amplitude of distorted CiB in weak turbulence can be expressed as the superposition of spiral harmonics [16]
Ψp,l0(q0,q1)(r,φ,z)=12πl=βp,l0,l(q0,q1)(r,z)exp(ilφ)
with the expansion coefficient
βp,l0,l(q0,q1)(r,z)=12π02πCiBp,l0(q0,q1)(r,φ,z)exp[iψ(r,φ,z)]exp(ilφ)dφ.

Instead of the random variable βp,l0,l(q0,q1)(r,z), we are usually interested in the ensemble average over the turbulence statistics, i.e.

|βp,l0,l(q0,q1)(r,z)|2=12π02π02πCiBp,l0(q0,q1)(r,φ1,z)CiBp,l0(q0,q1)*(r,φ2,z)×exp[il(φ1φ2)]×exp{i[ψ(r,φ1,z)ψ(r,φ2,z)]}dφ1dφ2

The ensemble average in right hand of Eq. (5) reads [17]

exp{i[ψ(r,φ1,z)ψ(r,φ2,z)]}=exp[2r22r2cos(φ1φ2)ρ02],
where ρ0 denotes the spatial coherence radius of a spherical wave propagating in non-Kolmogorov turbulence [18]
ρ0=[8α2Γ(2α2)]12[2(α1)Γ(3α2)πΓ(2α2)k2Cn2z]1α2(3<α<4),
where Cn2 denotes refractive-index structure constant and α denotes non-Kolmogorov turbulence parameter. Equation (7) implies that ρ0 → ∞ if α → 3, and ρ0 → 0 if α → 4. Following Eqs. (5)(7) and using Eq. 8.411.1 in [19], we can obtain
|βp,l0,l(q0,q1)(r,z)|2=12|l0||l0|!Ψp,l0(ξ)(kW0|q(z)|)2|l0|+2|[(1+ξ)q˜(z)q(z)]p|×exp(k2r2W022|q(z)|2)|F11(p2,|l0|+1;r2χ2(z))|2,×r2|l0|exp(2r2ρ02)Ill0(2r2ρ02)
where Ill0 is the modified Bessel function of the first kind with order ll0. By normalizing the detected power Pl=0+|βp,l0,l(q0,q1)(r,z)|2rdr in spiral harmonics with azimuthal number l, the corresponding power weight reads [16]
Cl=Plm=Pm.

Here Cl means the power weight of OAM mode with azimuthal number l and in OAM-based optical communication it’s also the detection probability of OAM state in receiver. Especially, Cl0 is the power weight corresponding to the original OAM mode with azimuthal number l0. Without turbulence and with infinite aperture, Cl0 = 1 but Cll0 = 0. In the presence of turbulence, the transmitted power leak into near modes and therefore Cll0 is called as crosstalk power weight in context.

Because the crosstalk power weight is symmetric about l0 due to In(x) = In(x), it’s convenient to discuss the case of ll0 only. Under the condition of ll0, the infinite series form of Cl can be obtained using Eq. 6.622 in [19] as

Cl=i(1)|l0|+1kW0ρ02π|q(z)||l0|!Ψp,l0(ξ)(η1η+1)|l0|2+14|[(1+ξ)q˜(z)q(z)]p|×m=0n=0(p*2)m(p2)n(|l0|+1)m(|l0|+1)nn!m!×[ρ022(χ2)*(η21)12]m[ρ022χ2(η21)12]n×Qll012|l0|+n+m+12(η),
where η=k2W02ρ024|q(z)|2+1, (*)n denotes the Pochhammer symbol and Qνμ(*) denotes the associated Legendre functions of the second kind.

3. Numerical results and discussion

In this section, we discuss the relation of self-focusing effect to the parameters of CiB at first and then investigate the influence of non-Kolmogorov turbulence on the original OAM mode. W0 relates to the scale of beam and is taken as 0.04m here. The imaginary part of q1, which is similar to that of q0, is taken as 1200m here. Besides, wavelength is taken as λ = 1550nm in the paper. For quantitative analysis of the focusing property of CiB, we use the second order momentum to describe the beam width as

W2(z)=4r2=4r2|CiBp,l0(q0,q1)(r,z)|2d2r(z0),
where r⃗ denotes the point vector on transverse plane and * denotes the ensemble average. The waist, i.e. the minimal beam width, is found at z = ze. ze is also the position where the energy of beam is most concentrated. If ze > 0m, the CiB converges at ze and we call the beam is self-focusing in context. If ze=0m, the CiB is not self-focusing.

In Fig. 1, We investigate the effect of some parameters of the CiB on ze. This curves In Figs. 1(a)–1(d) are non-differentiable at the joint point of ze=0m and ze > 0m due to the mandatory condition z ≥ 0m in Eq. (11) resulting in ze ≥ 0m. In Figs. 1(a) and 1(c), we plot the waist position ze against beam parameters Re(q0) and Re(q1). ze decreases as Re(q1) increases or Re(q0) decreases. We observe an approximate linear relation between Re(q1) and ze, which shows how the waist of Gaussian envelope affects the waist of the CiB. If the value of Re(q1) is greater than about 2000m, ze reaches 0m that means the CiB is not self-focusing. From Figs. 1(a) and 1(c), We suppose that the value of ze can be changed greatly if Re(q1) is varied adequately. In Figs. 1(b) and 1(d), we plot the waist position ze against beam parameter p. As Im(p) increases from 0 to 20, ze increases at first and then decreases. ze increases as Re(p) increases. It is clear that when the value of Im(p) is smaller than about 0, ze is 0m, that is seemingly independent of Re(p). It seems that the value of ze can but be changed into a limited range if only Im(p) is varied. The information contained in Fig. 1 indicate that one can manipulate ze by varying p, Re(q0) and Re(q1).

 figure: Fig. 1

Fig. 1 The waist position ze against (a) Re(q0) and Re(q1) with p=2+20i; (b) p with Re(q0)=Re(q1)=0m; (c) Re(q0) and Re(q1) with p=2+10i; (d) p with Re(q0)=Re(q1)=1000m. Other parameters: l0 = 1, Cn2=1015m3α and α = 3.67, respectively.

Download Full Size | PDF

Figure 2(a) shows the location of waist of one kind of self-focusing CiB in weak non-Kolmogorov turbulence. The location of waist of this CiB is at ze=2264m. The difference between the beam width at ze and the one at z=0m is about 100mm. The phase patterns of distorted CiB in weak non-Kolmogorov turbulence are shown in Figs. 2(b)–2(d). In these phase patterns, The gray scales from black to white correspond to the phase varying from 0 to 2π. Figure 2(b) shows the phase pattern at z=0m. This undistorted azimuthal phase increases counterclockwise while the phase contour lines twine clockwise. Figure 2(c) shows the distorted phase pattern at the waist position z = ze=2264m. A serious distorted azimuthal spiral phase is difficult to recognize at farther position such as z=5000m shown in Fig. 2(d).

 figure: Fig. 2

Fig. 2 (a) beam width W (z) of the CiB against propagation distance z in non-Kolmogorov turbulence. Phase patterns of the CiB on xy plane of (b) z=0m, (c) z=2264m and (d) z=5000m in non-Kolmogorov turbulence. Other parameters: Cn2=1015m3α, α=3.67, q1=1200i m, l0=1 and p=2+20i, respectively.

Download Full Size | PDF

In Fig. 3, we plot the received power weight Cl0 against propagation distance z in different cases. Figure 3(a) is plotted for different values of parameter p. All the lines show the similar character that as z increases, Cl0 decreases except for an arch near ze. Since the previous works show that detection probability of initial OAM state decreases as the relative beam width W (z)/r0 increases [7], where r0 is Fried parameter, here the presence of arch can be explained as the case that the ratio W (z)/ρ0 decreases as z increases. Around ze, Cl0 varies slowly. This mechanism proposes a novel method to alleviate the effect of turbulence on free space communication. In other words, one can manipulate ze according to the position of the receiver to obtain maximal Cl0 or alleviate the influence of atmospheric turbulence.

 figure: Fig. 3

Fig. 3 The received power weight Cl0 for the CiB against propagation distance z (a) with different Re(p) where Cn2=1015m3α, α=3.97, Im(p) = 20 and l0=1; (b) with different α where Cn2=1015m3α, l0=1, p=2+20i; (c) with different Cn2 where α=3.67, l0=1, p=2+20i; (d) with different l0 where Cn2=1015m3α, α=3.67, p=2+20i. Other parameter: q1=1200i m.

Download Full Size | PDF

Figure 3(a) also shows that as Re(p) increases, ze increases, but meanwhile the value of Cl0 at ze decreases. For example, Re(p) increases from 2 to 8, ze expands from about 2.25km to about 3.0km while Cl0 at ze decreases from about 0.93 to 0.89. Parameter Re(p) of CiB has similar meaning to the radial index of LGB, and bigger value of Re(p) results in larger beam width that makes the OAM of CiB be more vulnerable to spatial aberrations in atmospheric turbulence.

Figure 3(b) shows that Cl0 decreases as α increases. With different α, the curves also show an arch near ze, But with smaller α the curve is flatter and the arch is not clear. Figure 3(c) shows similar character that Cl0 decreases as Cn2 increases. While turbulence is weak, e.g. Cn2=7.5×1016m3α, Cl0 is almost 1 from 0km to 5km. In Figs. 3(b) and 3(c), it is obvious that the self-recovery effect of self-focusing CiB is more clear in the strong perturbation case i.e. large Cn2 or α than that in the weak perturbation case. Figures 3(b) and 3(c) also show that ze is almost the same and not relevant to the weak atmospheric turbulence. This conclusion is approximative because the weak turbulence is presupposed to affect the phase of the beam only in Eq. (2). But this presupposition is unreasonable in strong turbulence.

Figure 3(d) is plotted for different l0. Cl0 decreases as l0 increases, because bigger l0 means larger beam width W(z) that makes the OAM of beam be more vulnerable to turbulence, which is similar to LGB propagating in atmospheric turbulence in [20]. From Fig. 3(d) we also find that as l0=1 and 2, the arch emerges and the recovery effect of CiB arises. As l0=3,4 and 5, the local arches are not so clear and Cl0 descends very slowly around ze. Therefore, the resistance ability of OAM mode is weaker with larger l0. To obtain better signal quality, it is better to employ lower OAM modes as channels in OAM-based free space communication. For example, the OAM modes with l0=−4,−3...3,4 can be employed for coaxial transmission in condition of Fig. 3(d).

In Fig. 4, we consider the crosstalk power weight Cl of CiB mode. Cl decreases as Δl = |ll0| increases. As z increases, Cl increases except a slight declining near ze=2.264km. The reason for this declining is that received power weight Cl0 reaches its maximum near ze and then less power from initial mode falls into other modes. As Δl = 1, this declining is most obvious. It is clear that the crosstalk power weights for Δl=2,3 or 4 are very small within short propagation distance and they can be omitted comparing with the crosstalk power weights for Δl=1. If some communication systems utilize CiB modes for short distance message transmission, the crosstalk is very small when nonadjacent modes are employed.

 figure: Fig. 4

Fig. 4 The crosstalk power weight Cl for the CiB against propagation distance z with Δl=|ll0|=1,2,3,4. Other parameters: l0=1, p=2+20i, q1=1200i m, α=3.67 and Cn2=1015m3α.

Download Full Size | PDF

4. Conclusions

In this paper, we have revealed the self-focusing property and the self-recovery effect of the CiB in weak non-Kolmogorov turbulence. This work is carried out theoretically and numerically. Our results show that the appearance of self-focusing property of the CiB strongly depends on Re(q1) and Im(p), and the location of waist of the CiB can be changed by varying Re(q0),Re(q1), Re(p) and Im(p). The key conclusion shows that the self-recovery effect may arise while the CiB has self-focusing property, which implies that self-focusing CiB may have more excellent resistance ability to turbulence effect on the position of waist than vortex beam with no self-focusing property. It is noticeable that the self-recovery effect is only available for the CiB with small azimuthal number l0. Anyway, it proposes an opportunity to alleviate the turbulent effects on OAM-based optical communication.

Funding

We acknowledge the support of the National Natural Science Foundation of China (NSFC) under Grant nos. 60908034, 61205122 and 11474048, and Fundamental Research Funds for the Central Universities under grant nos. ZYGX2013J052 and ZYGX2015J042.

References and links

1. L. Allen, M. Beijersbergen, R. Spreeuw, and J. Woerdman, “Orbital angular momentum of light and the transformation of Laguerre-Gaussian laser modes,” Phys. Rev. A 45(11), 8185–8189 (1992). [CrossRef]   [PubMed]  

2. A. Willner, H. Huang, Y. Yan, Y. Ren, N. Ahmed, G. Xie, C. Bao, L. Li, Y. Cao, Z. Zhao, J. Wang, M. Lavery, M. Tur, S. Ramachandran, A. Molisch, N. Ashrafi, and S. Ashrafi, “Optical communications using orbital angular momentum beams,” Adv. Opt. Photon. 7, 66–106 (2015).

3. G. Gibson, J. Courtial, M. Padgett, M. Vasnetsov, V. Pas’ko, S. Barnett, and S. Franke-Arnold, “Free-space information transfer using light beams carrying orbital angular momentum,” Opt. Express 12(22), 5448–5456 (2004). [CrossRef]   [PubMed]  

4. B. Rodenburg, M. P. J. Lavery, M. Malik, M. N. O’Sullivan, M. Mirhosseini, D. J. Robertson, M. Padgett, and R. W. Boyd, “Influence of atmospheric turbulence on states of light carrying orbital angular momentum,” Opt. Lett. 37(17), 3735–3737 (2012). [CrossRef]   [PubMed]  

5. G. A. Tyler and R. W. Boyd, “Influence of atmospheric turbulence on the propagation of quantum states of light carrying orbital angular momentum,” Opt. Lett. 34(2), 142–144 (2009). [CrossRef]   [PubMed]  

6. J. A. Anguita, M. Neifeld, and B. Vasic, “Turbulence-induced channel crosstalk in an orbital angular momentummultiplexed free-space optical link,” Appl. Opt. 47(13), 2414–2429 (2008). [CrossRef]   [PubMed]  

7. C. Paterson, “Atmospheric turbulence and orbital angular momentum of single photons for optical communication,” Phys. Rev. Lett. 94(15), 153901 (2005). [CrossRef]   [PubMed]  

8. C. Gopaul and R. Andrews, “The effect of atmospheric turbulence on entangled orbital angular momentum states,” New J. Phys. 9, 94 (2007). [CrossRef]  

9. Y. Zhu, L. C. Zhang, Z. D. Hu, and Y. X. Zhang, “Effects of non-Kolmogorov turbulence on the spiral spectrum of Hypergeometric-Gaussian laser beams,” Opt. Express 23(7), 9137–9146(2015). [CrossRef]   [PubMed]  

10. Y. Zhu, X. J. Liu, J. Gao, Y. X. Zhang, and F. S. Zhao, “Probability density of the orbital angular momentum mode of Hankel-Bessel beams in an atmospheric turbulence,” Opt. Express 22(7), 7765–7772 (2014). [CrossRef]   [PubMed]  

11. Y. X. Zhang, M. J. Cheng, Y. Zhu, J. Gao, W. Y. Dan, Z. D. Hu, and F. S. Zhao, “Influence of atmospheric turbulence on the transmission of orbital angular momentum for Whittaker-Gaussian laser beams,” Opt. Express 22(18), 22101–22110 (2014). [CrossRef]   [PubMed]  

12. M. A. Bandres and J. C. Gutierrez-Vega, “Circular beams,” Opt. Lett. 33(2), 177–179 (2008). [CrossRef]   [PubMed]  

13. G. Vallone, “On the properties of circular beams: normalization, Laguerre-Gauss expansion, and free-space divergence,” Opt. Lett. 40(8), 1717–1720 (2015). [CrossRef]   [PubMed]  

14. P. Zhang, J. Prakash, Z. Zhang, M. S. Mills, N. K. Efremidis, D. N. Christodoulides, and Z. Chen, “Trapping and guiding microparticles with morphing autofocusing Airy beams,” Opt. Lett. 36(15), 2883–2885 (2011). [CrossRef]   [PubMed]  

15. N. K. Efremidis and D. N. Christodoulides, “Abruptly autofocusing waves,” Opt. Lett. 35(23), 4045–4047 (2010). [CrossRef]   [PubMed]  

16. L. Torner, J. P. Torres, and S. Carrasco, “Digital spiral imaging,” Opt. Express 13(3), 873–881 (2005). [CrossRef]   [PubMed]  

17. L. C. Andrews and R. L. Phillips, Laser Beam Propagation Through Random Media (SPIE, 2005). [CrossRef]  

18. C. H. Rao, W. H. Jiang, and N. Ling, “Spatial and temporal characterization of phase fluctuations in non-Kolmogorov atmospheric turbulence,” J. Mod. Opt. 47(6), 1111–1126 (2000). [CrossRef]  

19. I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series and Products, 7th ed. (Academic, 2007).

20. Y.-D. Liu, C. Q. Gao, M. W. Gao, and F. Li, “Coherent-mode representation and orbital angular momentum spectrum of partially coherent beam,” Opt. Commun. 281(8), 1968–1975 (2008). [CrossRef]  

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (4)

Fig. 1
Fig. 1 The waist position ze against (a) Re(q0) and Re(q1) with p=2+20i; (b) p with Re(q0)=Re(q1)=0m; (c) Re(q0) and Re(q1) with p=2+10i; (d) p with Re(q0)=Re(q1)=1000m. Other parameters: l0 = 1, C n 2 = 10 15 m 3 α and α = 3.67, respectively.
Fig. 2
Fig. 2 (a) beam width W (z) of the CiB against propagation distance z in non-Kolmogorov turbulence. Phase patterns of the CiB on xy plane of (b) z=0m, (c) z=2264m and (d) z=5000m in non-Kolmogorov turbulence. Other parameters: C n 2 = 10 15 m 3 α , α=3.67, q1=1200i m, l0=1 and p=2+20i, respectively.
Fig. 3
Fig. 3 The received power weight C l 0 for the CiB against propagation distance z (a) with different Re(p) where C n 2 = 10 15 m 3 α , α=3.97, Im(p) = 20 and l0=1; (b) with different α where C n 2 = 10 15 m 3 α , l0=1, p=2+20i; (c) with different C n 2 where α=3.67, l0=1, p=2+20i; (d) with different l0 where C n 2 = 10 15 m 3 α , α=3.67, p=2+20i. Other parameter: q1=1200i m.
Fig. 4
Fig. 4 The crosstalk power weight Cl for the CiB against propagation distance z with Δl=|ll0|=1,2,3,4. Other parameters: l0=1, p=2+20i, q1=1200i m, α=3.67 and C n 2 = 10 15 m 3 α .

Equations (11)

Equations on this page are rendered with MathJax. Learn more.

CiB p , l 0 ( q 0 , q 1 ) ( r , φ , z ) = ( i 2 z 0 W 0 ) | l 0 | + 1 [ π | l 0 | ! Ψ p , l 0 ( ξ ) ] 1 2 1 q ( z ) exp [ i k r 2 2 q ( z ) ] × [ ( 1 + ξ ) q ˜ ( z ) q ( z ) ] p 2 [ r q ( z ) ] | l 0 | F 1 1 ( p 2 , | l 0 | + 1 ; r 2 χ 2 ( z ) ) . × exp ( i l 0 φ )
Ψ p , q 0 ( q 0 , q 1 ) ( r , φ , z ) = CiB p , l 0 ( q 0 , q 1 ) ( r , φ , z ) exp [ i ψ ( r , φ , z ) ] ,
Ψ p , l 0 ( q 0 , q 1 ) ( r , φ , z ) = 1 2 π l = β p , l 0 , l ( q 0 , q 1 ) ( r , z ) exp ( i l φ )
β p , l 0 , l ( q 0 , q 1 ) ( r , z ) = 1 2 π 0 2 π CiB p , l 0 ( q 0 , q 1 ) ( r , φ , z ) exp [ i ψ ( r , φ , z ) ] exp ( i l φ ) d φ .
| β p , l 0 , l ( q 0 , q 1 ) ( r , z ) | 2 = 1 2 π 0 2 π 0 2 π CiB p , l 0 ( q 0 , q 1 ) ( r , φ 1 , z ) CiB p , l 0 ( q 0 , q 1 ) * ( r , φ 2 , z ) × exp [ i l ( φ 1 φ 2 ) ] × exp { i [ ψ ( r , φ 1 , z ) ψ ( r , φ 2 , z ) ] } d φ 1 d φ 2
exp { i [ ψ ( r , φ 1 , z ) ψ ( r , φ 2 , z ) ] } = exp [ 2 r 2 2 r 2 cos ( φ 1 φ 2 ) ρ 0 2 ] ,
ρ 0 = [ 8 α 2 Γ ( 2 α 2 ) ] 1 2 [ 2 ( α 1 ) Γ ( 3 α 2 ) π Γ ( 2 α 2 ) k 2 C n 2 z ] 1 α 2 ( 3 < α < 4 ) ,
| β p , l 0 , l ( q 0 , q 1 ) ( r , z ) | 2 = 1 2 | l 0 | | l 0 | ! Ψ p , l 0 ( ξ ) ( k W 0 | q ( z ) | ) 2 | l 0 | + 2 | [ ( 1 + ξ ) q ˜ ( z ) q ( z ) ] p | × exp ( k 2 r 2 W 0 2 2 | q ( z ) | 2 ) | F 1 1 ( p 2 , | l 0 | + 1 ; r 2 χ 2 ( z ) ) | 2 , × r 2 | l 0 | exp ( 2 r 2 ρ 0 2 ) I l l 0 ( 2 r 2 ρ 0 2 )
C l = P l m = P m .
C l = i ( 1 ) | l 0 | + 1 k W 0 ρ 0 2 π | q ( z ) | | l 0 | ! Ψ p , l 0 ( ξ ) ( η 1 η + 1 ) | l 0 | 2 + 1 4 | [ ( 1 + ξ ) q ˜ ( z ) q ( z ) ] p | × m = 0 n = 0 ( p * 2 ) m ( p 2 ) n ( | l 0 | + 1 ) m ( | l 0 | + 1 ) n n ! m ! × [ ρ 0 2 2 ( χ 2 ) * ( η 2 1 ) 1 2 ] m [ ρ 0 2 2 χ 2 ( η 2 1 ) 1 2 ] n × Q l l 0 1 2 | l 0 | + n + m + 1 2 ( η ) ,
W 2 ( z ) = 4 r 2 = 4 r 2 | CiB p , l 0 ( q 0 , q 1 ) ( r , z ) | 2 d 2 r ( z 0 ) ,
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.