Abstract
Based on the generalized Huygens Fresnel integral, we derive the analytical formula of the cross-spectral density of a twisted partially coherent array beam propagating in non-Kolmogorov anisotropic turbulence, and investigate the changes in orbital angular momentum (OAM). The results show that the anisotropy of the turbulence causes different effects in horizontal and vertical directions. The spectral density distribution of twisted partially coherent array beam in turbulence presents self-splitting and rotation, which combines the interesting effects of the twist phase and coherent structure. Although OAM is conserved, the spatial distribution of OAM flux density can be changed by changing the propagation distance, power and anisotropy of turbulence, and the modulation of the twist phase affects not only the magnitude of OAM but also its distribution. Our work is helpful for exploring new forms of OAM sources, and promote the application of free-space optical communications and optical field modulation.
© 2022 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement
1. Introduction
In recent years, orbital angular momentum (OAM) has attracted extensive attention because of its great application prospects in the optical communication field [1–5]. OAM has an infinite number of mutually orthogonal modes, and modulating each signal on a different mode will greatly enhance the channel capacity and spectrum utilization of optical communication systems. Gibson et al. encoded the information with OAM to achieve 15m free-space optical communication, and the information encoded in this way has strong confidentiality [1]. Tamburini et al. made OAM multiplexing in the experiment for the first time and realized 442m wireless optical communication [2]. Bozinovic et al. achieved 1.6 Tbit/s information transmission by using a specially designed ring fiber to transmit the OAM beam [3]. Xie et al. designed an ultra-wideband OAM transmitter for 1.2Tbit/s OAM communication [4]. Pang et al. effectively reduce power loss and modal crosstalk in free-space optical communication by simultaneous transmission of four signals through mode orthogonalization [5]. In addition to the communication field, OAM also has broad application prospects in many other fields, such as optical wrenches [6], designing optical micromachines to manipulate particles [7], quantum memory [8], detecting rotating objects [9], etc.
A simple way to make the beam carry OAM is to modulate the phase of the beam. As an unconventional phase structure, the vortex phase has attracted the interest of researchers because of its spiral wavefront property during transmission [10]. The number of OAM that can be carried by each photon of a fully coherent vortex beam is lћ, l is the topological charge number, ћ is the approximate Planck constant. Unlike the vortex phase, in the partially coherent light field, there is also a special phase that can carry OAM called twist phase [11–18]. The twist phase is a position-dependent quadratic function that cannot be simply decomposed into the product of two one-dimensional covariates, therefore, can only exist in partially coherent beams. It is interesting to note, partially coherent beams carrying twist phase can be generated by superposition of vortex beams [12]. In 1993, Simon et al. introduced the concept of twist phase and introduced it into the Gaussian Schell-model beam to obtain the twisted Gaussian Schell-model(TGSM) beam [13], which rotate around the axis during transmission. Then the TGSM beam was generated by Friberg using an elliptical Gaussian Schell-model beam and two sets of cylindrical lenses [14]. When a twisted beam is used as an illumination source, the resulting image resolution can exceed the Rayleigh limit [15]. In addition, twisted beams also have potential applications in ghost imaging [16] and particle capture [17]. It is now well appreciated that twisted partially coherent beams can suppress the negative effects induced by turbulence [18]. However, due to the difficulty of the construction of the new correlation structure, the above studies are confined to the TGSM beam. Many thanks to Gori et al. for proposing a method for constructing special correlated partially coherent beams [19] and extended it to vector beams [20], variety of interesting new special correlated beams had been proposed [21–24]. Multi-Gaussian Schell-model beams can produce a flat-top distribution in the far field [21]. Hermite-Gaussian correlated beam splits itself during propagation [22]. Self-focusing of non-uniformly correlated beams [23]. Arbitrarily manipulable partially coherent array beams [24]. When adding consideration of twist phase, Borghi et al. provided the proof condition of whether the axisymmetric partially coherent beam can carry the twist phase based on the modal analysis method [25]. Gori et al. proposed a method to generate a twisted beam without symmetric constraints [26]. Many special correlated twisted partially coherent beams have also been generated [27–32], such as radially polarized twisted beams [27], twisted Laguerre-Gaussian Schell-model beams [28], twisted vortex Gaussian Schell-model beams [29], twisted array beams [30,31], and twisted elliptical multi-Gaussian Schell-model beams [32].
In the past, the calculation of partially coherent beam OAM mainly used the method of finding second-order moment parameters, which required complex matrix operations [11]. In 2012, Kim et al. proposed another method of calculating OAM, which can be realized only by simple differentiation [33], compared with the second-moment method, the calculation process was greatly simplified. Later, this method was used to calculate the OAM of twisted Laguerre-Gaussian Schell-model beams [34] and twisted Hermite Gaussion Schell-model beams [35] in free space and atmospheric turbulence, and the results showed that the total average OAM per photon was conserved when transmitting in both free space and atmospheric turbulence. Zhang et al. found that different OAM flux density distributions can be controlled by changing the initial coherence of the beam [36]. Due to the important applications of lasers for free-space-optical communication [37–39], their interaction with turbulence has become an important issue that researchers must consider. It is shown that partially coherent beams can well suppress the negative effects of beam propagation in atmospheric turbulence [40–48]. For a long time, the isotropic Kolmogorov power spectrum model has been used to study atmospheric turbulence. However, results have shown that the turbulence above the atmospheric boundary layer (1 to 2km from the earth's surface) does not always follow the Kolmogorov power spectrum model, and the general non-Kolmogorov turbulence statistics should be considered [49–51]. Moreover, turbulence in the free atmosphere has strong inhomogeneity and anisotropy at large scales [52,53]. Therefore, the non-Kolmogorov anisotropic turbulence power spectrum model is a better choice for simulating beam propagating in the atmosphere turbulence [54].
In our previous work, we proposed a class of twisted Schell-model array correlated sources [31], which gradually split the spot into arrays during propagating and can suppress the effects of atmospheric turbulence. In this paper, the OAM variation of twisted array beams have attracted our interest, the OAM evolution characteristics of twisted radially and rectangularly symmetric array beams in anisotropic turbulence are further investigated, and the results have potential applications in the field of laser atmospheric communication.
2. Theoretical model
First, we recall the cross-spectral density (CSD) function of a twisted radial array beam with an elliptical amplitude profile at the source plane can be expressed as [31]
According to the generalized Huygens-Fresnel principle, the CSD function of the beam propagation in atmospheric turbulence can be expressed as [18,38]
In this paper, we study a twisted partially coherent beam which carrying OAM. For scalar partially coherent beams, the OAM flux density represents the spatial distribution of OAM in the beam cross-section, and the OAM flux density along the z-axis can be expressed as [33]
Substituting Eq. (1) and Eq. (3) into Eq. (2), the analytical formula for the CSD function of a twisted radial array beam propagating in non-Kolmogorov anisotropic turbulence is obtained
The expression of the CSD function of the twisted rectangular array beam is obtained by referring to [31]
Except for Eq. (12), other parameter substitutions are similar to the radial array beam, ${C_0} = (2M + 1)(2N + 1)$ is the normalization factor, a denotes row spacing of array. By substituting the Wcir or Wrec into Eqs. (5)-(8), the spatial distribution of OAM flux density, spectral intensity and DOC of twisted array beams in atmospheric turbulence can be calculated.
3. Numerical simulation and analysis
In this section, we simulate the propagation characteristics of twisted array beams in anisotropic turbulence. Unless otherwise specified, the parameters are set to σx = 4 cm, σy = 2 cm, δ0 = 1 cm, λ=632.8 nm, εx = 1, εy = 3, ξ=11/3, μ=-0.5km-1, $C_n^2 = 5 \times {10^{ - 15}}{m^{3 - \xi }}$, ${M_{ra}} = 0$, ${N_{ra}} = 3$, M = N = 1, r0 = 0.15mm-1, a = 5 mm. Figure 1 depict the spatial distribution of the spectral intensity, the modulus of the DOC $\eta ({x,y,0,0} )$ at ${{\mathbf r}_1}\textrm{ = }({x,y} )$ and ${{\mathbf r}_2}\textrm{ = }({0,0} )$, OAM flux density, and normalized OAM flux density at the source plane for the twisted radial array and rectangular array beam, respectively. By observing Figs. 1(a2), (c2) and (d2), we can see that the DOC $\eta ({x,y,0,0} )$ of twisted radial array beam have the same distribution for different σy and µ, and in Fig. 1(b2) the rectangular array beam has different lattice coherence distribution at the source plane, this is attributed to the difference in the structure of the correlation [31]. It is of interest to find from Figs. 1(a3)–1(b3), for circularly symmetrical amplitude profiles, the distribution of OAM flux density at the source plane is circular with a dark hollow in the center, which is very similar to the OAM flux density of a partially coherent vortex beam [33]. This phenomenon further verifying the conclusion [12] that a twisted partially coherent beam can be generated by the superposition of vortex beams. As can be seen from Fig. 1(c), when the beam waist width in x and y directions are different, the distribution of spectral intensity and OAM flux density at the source plane becomes elliptical, while the DOC are not affected. With the addition of the twist factor in Fig. 1(d), we find that both the OAM flux density and the normalized OAM flux density will increase by a fixed multiple, indicating that modulation of the twist phase can control the magnitude of OAM.
Figure 2 illustrates the normalized spectral intensity distribution of the twisted radial array beam propagating in atmospheric turbulence. The spectral intensity gradually splits from a single spot to four spots distributed in a radial array, which is the same as [24]. Compare Figs. 2(a) and 2(b), we find that the intensity distribution in Fig. 2(b) gradually becomes elliptical. The reason for this phenomenon is that when the beam propagating in anisotropic turbulence, the spectral intensity distribution in the near field is mainly influenced by the initial beam parameters. After a certain propagation distance, the effect of turbulence will dominate due to the cumulative effect of turbulence, which will cause anisotropic beam spreading due to the different statistical properties in the horizontal and vertical directions. In Figs. 2(c)–2(d), the beam waist width is set to σx = 4 cm, σy = 2 cm in order to better observe the rotation of the spot. We can clearly observe that the beam is self-splitting while the individual sub spots are accompanied by rotation, which is consistent with the general twisted GSM beam [13].
As we all know, the twisted beam can rotate around the propagation axis during propagating and carrying OAM. Figure 3 shows the DOC $\eta ({u,v,0,0} )$ of the twisted radial array beam propagating in turbulence, the parameter settings correspond to Fig. 2. It is clearly seen that as the propagation distance increases, the initial lattice-like DOC is disrupted. The DOC distribution rotates due to twisted phase, but different source parameters result in different rotation angles. Compared with Figs. 2(a)–2(b) and Figs. 3(a)–3(b), one finds that the effect of anisotropic turbulence on the DOC distribution is more intuitive compared to the intensity, the rotation center of the spectral density splits as the beam splits during propagating, while the rotation center of the DOC remains invariant.
The OAM flux density distribution of the twisted radial array beam at different distances in turbulence is shown in Fig. 4, the evolution of the OAM flux density is similar to that of the spectral intensity, which splits into a multi-petal distribution, rotates and diffuses with the increase of propagation distance. Unlike the uniform distribution of the source field in Figs. 1(a3)-1(d3), the OAM flux density has negative values during transmission, where the OAM rotation direction is opposite to the positive region. The OAM negative region gradually increases with the increase of the distance. This related to the array distribution, which redistributes the OAM over a larger radial distance [34]. Integrating the OAM flux density in the plane yields the total OAM of the beam, which is conserved according to [33]. The spatial distribution of different OAM flux densities can be produced by tuning the anisotropy parameter of turbulence, but it does not affect the magnitude of the total OAM.
Furthermore, Fig. 5 shows the characteristics in the far field, parameters correspond to Fig. 2(b), 3(b) and 4(b). It is worth noting that the twisted beam rotation reaches π/2 as $z \to \infty$ in free space [13,14]. As shown in Figs. 5(a) and 5(b), the rotation of intensity and DOC distribution is no longer apparent in the far-field. Observation of Figs. 5(a1)–5(a4) reveals that when the beam is transmitted to the far field, the cumulative effect of turbulence becomes stronger as the propagation distance increases, and the beam spreading caused by the anisotropy of the turbulence causes the beam to lose its array distribution [31]. Due to the decoherence effect of turbulence, the contour of DOC distribution gradually decreases, and shows elliptical distribution due to the anisotropy of turbulence in Figs. 5(b3)–5(b4). The beam spreading causes the OAM flux density to redistribute and obtains an array OAM flux density distribution in the far field as Fig. 5(c4).
Figure 6 shows the normalized OAM flux density distribution of the twisted radial array beam at different distances in atmospheric turbulence. The normalized OAM flux density represents the average OAM per photon at a specific location in space. It is clear from Eq. (6) that the calculation of the normalized OAM flux density OAM depends on the ratio of OAM flux density to intensity, this results in regions with very low or 0 intensity have a large normalized OAM. In the numerical simulation, we filtered these invalid ranges, these erroneous high values of normalized OAM flux density near the edges/corners are not displayed. In Fig. 6, the distribution of normalized OAM flux density is generally low in the central region and high in the edge region. Which is similar to the distribution of the multi-Gaussian Schell mode vortex beam in [36]. In addition, due to the effects of anisotropic turbulence, negative values of the normalized OAM flux density appear in the central region.
Figure 7 demonstrates the normalized spectral intensity, OAM flux density, and normalized OAM flux density distributions of the twisted radial array beam propagating in anisotropic turbulence at z = 5 km, with different turbulence power exponent $\xi$. For non-Kolmogorov anisotropic turbulence as described by Eq. (3), we restrict it to $3 < \xi < 4$ since beyond this interval the wave structure function is not defined [50] and some of the integrals involved in propagation laws may diverge. Comparing Figs. 7(a1)–7(c1) and Figs. 7(a4)–7(c4), it is found that in the case near the extreme, both the normalized spectral intensity, OAM flux density and normalized OAM flux density distributions have significantly difference. The physical interpretation is on the extreme of $\xi \to 3$, the turbulence tends to vanish. On the other extreme of $\xi \to 4$, the phase effects dominate in the form of random tilts, which cause beam wander [50,54].
Figure 8 demonstrates the normalized intensity, OAM flux density, and normalized OAM flux density distribution for an elliptical twisted radial array beam with different twist factors in anisotropic turbulence at z = 5 km. By observing Figs. 8(a1)–8(a4), it is obtained that as the twist factor increases bringing about the spread of the beam [31], which requires a reasonable choice of the initial optical parameters in practical applications. As can be seen from Figs. 8(b1) and 8(c1), when the distortion factor is 0, the region of OAM flux density and normalized OAM flux density is equal in size and opposite in direction, indicating that the beam does not carry OAM at this time. As the distortion factor increases, the positive region of OAM flux density is larger than the negative region, and the total OAM of the beam is positive and shows a counterclockwise rotation, which can explain the rotation of the light spot phenomenon. As an important conclusion, we found that the modulation of the twist phase affects not only the magnitude of OAM but also its distribution.
In order to visualize the effect of array structure distributions on the evolution of the OAM in anisotropic turbulence. Based on Eq. (11) and (12), Fig. 9 demonstrates the normalized spectral intensity, OAM flux density, and normalized OAM flux density distribution of twisted rectangular array beam in anisotropic turbulence at z = 5 km. Figures 9(a1)–9(c1) N = M = 0, i.e. there is only one spot and no array, Figs. 9(a4)–9(c4) N = M = 1, for 3 × 3 array. As in Figs. 9(a1)–9(a4), the number of intensity array spots gradually increases by modulating the initial coherent structure. Moreover, arbitrary array intensity distributions can be generated. From Figs. 9(b1)–9(b4), due to the effect of anisotropic turbulence, the distributions of the OAM flux density is not a simple summation of the array beam. As we observe Figs. 9(b)–9(c), it is found that the change of the array structure not only affects the distribution structure of the OAM flux density and normalized OAM flux density, but also affects their numerical magnitudes.
4. Conclusion
In this paper, we simulate the propagating of twisted partially coherent array beams in non-Kolmogorov anisotropic turbulence, and analyze the evolution of the spectral intensity and OAM properties in turbulence. The results show when all parameters are ideally isotropic (circularly symmetric), the OAM flux density of the twisted array beam exhibits a circular distribution with a dark hollow in the center as show at the source plane. which is the same as that of the partially coherent vortex beam, further verifying the conclusion that the twisted beam can be generated by the superposition of vortex beams. However, the anisotropy of the turbulence causes anisotropic intensity distribution in horizontal and vertical directions. and the spectral density distribution rotates with increasing distance in turbulence. The twisted array beams have both the rotational characteristics of the TGSM beam and the self-splitting characteristics of the array beam. Although OAM is conserved, the spatial distribution of OAM flux density can be changed by changing the propagation distance, power and anisotropy of turbulence, and the modulation of the twist phase affects not only the magnitude of OAM but also its distribution. Our work is helpful for exploring new forms of OAM sources. The results of this paper have potential applications in free-space optical communication and optical field modulation.
Funding
National Natural Science Foundation of China (11904253); Open Research Fund of Key Laboratory of Atmospheric Optics, Chinese Academy of Sciences (HTAD-JJ-19-02); Taiyuan University of Science and Technology Scientific Research Initial Funding (20202013); Shanxi Provincial Central Government Guides Local Science and Technology Development Fund Project (YDZX20201400001386).
Disclosures
The authors declare no conflicts of interest.
Data availability
Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.
References
1. 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]
2. F. Tamburini, E. Mari, A. Sponselli, B. Thidé, A. Bianchini, and F. Romanato, “Encoding many channels on the same frequency through radio vorticity: first experimental test,” New J. Phys. 14(3), 033001 (2012). [CrossRef]
3. N. Bozinovic, Y. Yue, Y. Ren, M. Tur, P. Kristensen, H. Huang, A. E. Willner, and S. Ramachandran, “Terabit-scale orbital angular momentum mode division multiplexing in fibers,” Science 340(6140), 1545–1548 (2013). [CrossRef]
4. Z. Xie, T. Lei, F. Li, H. Qiu, Z. Zhang, H. Wang, C. Min, L. Du, Z. Li, and X. Yuan, “Ultra-broadband on-chip twisted light emitter for optical communications,” Light: Sci. Appl. 7(4), 18001 (2018). [CrossRef]
5. K. Pang, H. Song, Z. Zhao, R. Zhang, H. Song, G. Xie, L. Li, C. Liu, J. Du, A. F. Molisch, M. Tur, and A. E. Willner, “400-Gbit/s QPSK free-space optical communicationlink based on four-fold multiplexing of Hermite-Gaussian or Laguerre-Gaussian modes by varying both modal indices,” Opt. Lett. 43(16), 3889–3892 (2018). [CrossRef]
6. N. B. Simpson, K. Dholakia, L. Allen, and M. J. Padgett, “Mechanical equivalence of spin and orbital angular momentum of light: an optical spanner,” Opt. Lett. 22(1), 52–54 (1997). [CrossRef]
7. K. Ladavac and D. Grier, “Microoptomechanical pumps assembled and driven by holographic optical vortex arrays,” Opt. Express 12(6), 1144–1149 (2004). [CrossRef]
8. A. Nicolas, L. Veissier, L. Giner, E. Giacobino, D. Maxein, and J. Laurat, “A quantum memory for orbital angular momentum photonic qubits,” Nat. Photonics 8(3), 234–238 (2014). [CrossRef]
9. M. P. Lavery, F. C. Speirits, S. M. Barnett, and M. J. Padgett, “Detection of a spinning object using light's orbital angular momentum,” Science 341(6145), 537–540 (2013). [CrossRef]
10. L. Allen, M. W. Beijersbergen, R. J. Spreeuw, and J. P. Woerdman, “Orbital angular momentum of light and the transformation of Laguerre-Gaussian laser modes,” Phys. Rev. A 45(11), 8185–8189 (1992). [CrossRef]
11. J. Serna and J. M. Movilla, “Orbital angular momentum of partially coherent beams,” Opt. Lett. 26(7), 405–407 (2001). [CrossRef]
12. F. Gori and M. Santarsiero, “Twisted Gaussian Schell-model beams as series of partially coherent modified Bessel-Gauss beams,” Opt. Lett. 40(7), 1587–1590 (2015). [CrossRef]
13. R. Simon and N. Mukunda, “Twisted Gaussian Schell-model beams,” J. Opt. Soc. Am. A 10(1), 95–109 (1993). [CrossRef]
14. A. T. Friberg, E. Tervonen, and J. Turunen, “Interpretation and experimental demonstration of twisted Gaussian Schell-model beams,” J. Opt. Soc. Am. A 11(6), 1818–1826 (1994). [CrossRef]
15. Z. Tong and O. Korotkova, “Beyond the classical Rayleigh limit with twisted light,” Opt. Lett. 37(13), 2595–2597 (2012). [CrossRef]
16. Y. Cai, Q. Lin, and O. Korotkova, “Ghost imaging with twisted Gaussian Schell-model beam,” Opt. Express 17(4), 2453–2464 (2009). [CrossRef]
17. C. Zhao, Y. Cai, and O. Korotkova, “Radiation force of scalar and electromagnetic twisted Gaussian Schell-model beams,” Opt. Express 17(24), 21472–21487 (2009). [CrossRef]
18. Y. Cai and S. He, “Propagation of a partially coherent twisted anisotropic Gaussian Schell-model beam in a turbulent atmosphere,” Appl. Phys. Lett. 89(4), 041117 (2006). [CrossRef]
19. F. Gori and M. Santarsiero, “Devising genuine spatial correlation functions,” Opt. Lett. 32(24), 3531–3533 (2007). [CrossRef]
20. F. Gori, V. Ramírez-Sánchez, M. Santarsiero, and T. Shirai, “On genuine cross-spectral density matrices,” J. Opt. A: Pure Appl. Opt. 11(8), 085706 (2009). [CrossRef]
21. S. Sahin and O. Korotkova, “Light sources generating far fields with tunable flat profiles,” Opt. Lett. 37(14), 2970–2972 (2012). [CrossRef]
22. Y. Chen, J. Gu, F. Wang, and Y. Cai, “Self-splitting properties of a Hermite-Gaussian correlated Schell-model beam,” Phys. Rev. A 91(1), 013823 (2015). [CrossRef]
23. H. Lajunen and T. Saastamoinen, “Propagation characteristics of partially coherent beams with spatially varying correlations,” Opt. Lett. 36(20), 4104–4106 (2011). [CrossRef]
24. S. Zhu, J. Wang, X. Liu, Y. Cai, and Z. Li, “Generation of arbitrary radially polarized array beams by manipulating correlation structure,” Appl. Phys. Lett. 109(16), 161904 (2016). [CrossRef]
25. R. Borghi, F. Gori, G. Guattari, and M. Santarsiero, “Twisted Schell-model beams with axial symmetry,” Opt. Lett. 40(19), 4504–4507 (2015). [CrossRef]
26. F. Gori and M. Santarsiero, “Devising genuine twisted cross-spectral densities,” Opt. Lett. 43(3), 595–598 (2018). [CrossRef]
27. L. Liu, X. Peng, Y. Chen, L. Guo, and Y. Cai, “Statistical properties of a radially polarized twisted Gaussian Schell-model beam in a uniaxial crystal,” J. Mod. Opt. 64(7), 698–708 (2017). [CrossRef]
28. X. Peng, L. Liu, F. Wang, S. Popov, and Y. Cai, “Twisted Laguerre-Gaussian Schell-model beam and its orbital angular moment,” Opt. Express 26(26), 33956–33969 (2018). [CrossRef]
29. C. S. D. Stahl and G. Gbur, “Twisted vortex Gaussian Schell-model beams,” J. Opt. Soc. Am. A 35(11), 1899–1906 (2018). [CrossRef]
30. L. Wan and D. Zhao, “Twisted Gaussian Schell-model array beams,” Opt. Lett. 43(15), 3554–3557 (2018). [CrossRef]
31. J. Wang, H. Huang, Y. Chen, H. Wang, S. Zhu, Z. Li, and Y. Cai, “Twisted partially coherent array sources and their transmission in anisotropic turbulence,” Opt. Express 26(20), 25974–25988 (2018). [CrossRef]
32. H. Wang, X. Peng, L. Liu, F. Wang, and Y. Cai, “Twisted elliptical multi-Gaussian Schell-model beams and their propagation properties,” J. Opt. Soc. Am. A 37(1), 89–97 (2020). [CrossRef]
33. S. M. Kim and G. Gbur, “Angular momentum conservation in partially coherent wave fields,” Phys. Rev. A 86(4), 043814 (2012). [CrossRef]
34. Y. Liu, R. Lin, F. Wang, Y. Cai, and J. Yu, “Propagation properties of Laguerre-Gaussian Schell-model beams with a twist phase,” J. Quant. Spectrosc. Radiat. Transfer 264, 107556 (2021). [CrossRef]
35. R. Lin, H. Yu, X. Zhu, L. Liu, G. Gbur, Y. Cai, and J. Yu, “The evolution of spectral intensity and orbital angular momentum of twisted Hermite Gaussian Schell model beams in turbulence,” Opt. Express 28(5), 7152–7164 (2020). [CrossRef]
36. Y. Zhang, O. Korotkova, Y. Cai, and G. Gbur, “Correlation-induced orbital angular momentum changes,” Phys. Rev. A 102(6), 063513 (2020). [CrossRef]
37. A. Ishimaru, Wave Propagation and Scattering in Random Media (Academic, 1978).
38. L. C. Andrews and R. L. Phillips, Laser Beam Propagation Through Random Media (SPIE, 2005).
39. T. H. Yura, “Mutual Coherence Function of a Finite Cross Section Optical Beam Propagating in a Turbulent Medium,” Appl. Opt. 11(6), 1399–1406 (1972). [CrossRef]
40. J. C. Ricklin and F. M. Davidson, “Atmospheric turbulence effects on a partially coherent Gaussian beam: implications for free-space laser communication,” J. Opt. Soc. Am. A 19(9), 1794–1802 (2002). [CrossRef]
41. T. Shirai, A. Dogariu, and E. Wolf, “Mode analysis of spreading of partially coherent beams propagating through atmospheric turbulence,” J. Opt. Soc. Am. A 20(6), 1094–1102 (2003). [CrossRef]
42. O. Korotkova, “Model for a partially coherent Gaussian beam in atmospheric turbulence with application in Lasercom,” Opt. Eng. 43(2), 330–341 (2004). [CrossRef]
43. G. P. Berman, A. A. Chumak, and V. N. Gorshkov, “Beam wandering in the atmosphere: the effect of partial coherence,” Phys. Rev. E 76(5), 056606 (2007). [CrossRef]
44. Y. Gu and G. Gbur, “Scintillation of nonuniformly correlated beams in atmospheric turbulence,” Opt. Lett. 38(9), 1395–1397 (2013). [CrossRef]
45. Z. Mei and O. Korotkova, “Cosine-Gaussian Schell-model sources,” Opt. Lett. 38(14), 2578–2580 (2013). [CrossRef]
46. R. Chen, L. Liu, S. Zhu, G. Wu, F. Wang, and Y. Cai, “Statistical properties of a Laguerre-Gaussian Schell-model beam in turbulent atmosphere,” Opt. Express 22(2), 1871–1883 (2014). [CrossRef]
47. O. Korotkova and E. Shchepakina, “Rectangular Multi-Gaussian Schell-Model beams in atmospheric turbulence,” J. Opt. 16(4), 045704 (2014). [CrossRef]
48. G. Wu, H. Zhou, T. Zhao, H. Liu, Z. Li, and X. Zhao, “Propagation properties of electromagnetic multi-Gaussian Schell model beams propagating through atmospheric turbulence,” J. Korean Phys. Soc. 64(6), 826–831 (2014). [CrossRef]
49. L. C. Andrews, R. L. Phillips, and R. Crabbs, “Propagation of a Gaussian-beam wave in general anisotropic turbulence,” Proc. SPIE 9224, 922402 (2014). [CrossRef]
50. I. Toselli, O. Korotkova, X. Xiao, and D. G. Voelz, “SLM-based laboratory simulations of Kolmogorov and non-Kolmogorov anisotropic turbulence,” Appl. Opt. 54(15), 4740–4744 (2015). [CrossRef]
51. J. Wang, S. J. Zhu, and Z. H. Li, “Vector properties of a tunable random electromagnetic beam in non-Kolmogrov turbulence,” Chin. Opt. Lett. 14(8), 080101 (2016). [CrossRef]
52. J. Wang, H. Wang, S. Zhu, and Z. Li, “Second-order moments of a twisted Gaussian Schell-model beam in anisotropic turbulence,” J. Opt. Soc. Am. A 35(7), 1173–1179 (2018). [CrossRef]
53. I. Toselli, “Introducing the concept of anisotropy at different scales for modeling optical turbulence,” J. Opt. Soc. Am. A 31(8), 1868–1875 (2014). [CrossRef]
54. I. Toselli and O. Korotkova, “General scale-dependent anisotropic turbulence and its impact on free space optical communication system performance,” J. Opt. Soc. Am. A 32(6), 1017–1025 (2015). [CrossRef]