Abstract
We investigate the behavior of the orbital angular momentum (OAM) flux density of partially coherent vortex (PCV) beams in atmospheric turbulence. It is shown that for PCV beams with different spatial coherence structures, the OAM flux density distribution exhibits rich variations along the propagation path. Our findings provide insight into the use of the OAM in free-space optical communications when turbulence effects are significant.
Introduction
During the past few decades, the study of vortex structures in optical fields, and more complicated wavefield singularities, has become a stand-alone area of investigation of modern optics, named singular optics [1–3]. Spatially coherent beams carrying optical vortices, called vortex beams, have attracted much interest due to their potential applications in many areas including coronagraph [4], coherence filtering [5], and free-space optical communication [6]. Since Allen et al. found that Laguerre–Gauss beams carry a well-defined OAM as a consequence of their vortex core [7], vortex beams have also found application in optical tweezing [8], optical spanning [9], design of light-driven machines [10], high-resolution microscopy [11], security holograms [12], quantum information transfer [13], and rotation measurements [14].
Partially coherent beams have been shown to benefit a number of applications, such as free-space optical communication [15], particle trapping [16], and atom cooling [17]. It is therefore of great interest to consider vortex structures in partially coherent beams. Unlike spatially coherent beams, however, partially coherent beams do not have a well-defined phase structure and their behavior can only be described using coherence functions. Some years ago, Schouten et al showed that these coherence functions can possess their own phase singularities (and optical vortices) that are closely related to their coherent counterparts [18], and since then, much attention has been given to these so-called coherence singularities or coherence vortices [19–23]. For partially coherent beams carrying vortex structures, which are called PCV beams, when the coherence decreases, coherence singularities are shown to be more robust than phase singularities [24].
Recently, it was found that the OAM of PCV beams exhibits a number of novel characteristics. For example, different types of PCV beams can have different distributions of the OAM flux density in their cross-sections. The OAM flux density can represent a Rankine vortex, a rigid body rotator, a fluid rotator, or even more rich varieties of OAM distributions; these results indicate that PCV beams can provide additional control over OAM than their coherence counterparts [25–28]. More recently, we found that one can manipulate the source coherence to change the OAM flux density of PCV beams on propagation in free space [29], and this effect joins correlation-induced spectral and polarization changes in a family of source coherence influenced propagation phenomena (for the summary see Ch. 4 of Ref. [30]).
However, despite the fact that free-space optical communication is one of the most significant applications of such beams, the propagation characteristics of the OAM flux density of PCV beams in atmospheric turbulence has yet to be studied. In this paper, we study the propagation of the OAM flux density of PCV beams in atmospheric turbulence. We derive analytic formulas showing how atmospheric turbulence affects the OAM flux density on propagation, and illustrate the changes with a number of model partially coherent beams with structured spatial coherence.
Theory
We consider a scalar, statistically stationary random source that is located in a plane perpendicular to the direction of propagation. The second-order statistical properties of the source can be characterized by the cross-spectral density (CSD) function [31],where represents averaging over a spatial-frequency ensemble of the field , and the asterisk stands for complex conjugate. Here, r is the two-dimensional position vector in the source plane.
To model a variety of partially coherent vortex beams, we consider isotropic Schell-model sources with a definite topological charge, of the formwithrepresenting the normalized Laguerre–Gauss mode of radial order 0 and azimuthal order l, andwhere , is the source spectral degree of coherence, w is the initial beam width, and l is the topological charge.
The CSD function of the PCV beam after propagating in atmospheric turbulence can be calculated with the extended Collins integral [32].where and denote the coordinates of two arbitrary points at the receiver plane, is the wave number of light with being the wavelength. The angle brackets represent an average over an ensemble of turbulence states, which can be expressed as [33]where
Here, is the quantity denoting the effect of turbulence, is the one-dimensional power spectrum of the refractive-index fluctuations of the turbulent medium, and is the spatial frequency.
If turbulence is governed by non-Kolmogorov statistics and the power spectrum has the van Karman form, in which slope 11/3 is generalized to an arbitrary parameter [34], i.e.,withwhere and are the outer and inner scales of the turbulence, respectively; is the power law exponent; is a generalized refractive-index structure parameter with units , and is the gamma function. For the power spectrum in Eq. 7, can be expressed aswhere , is the incomplete Gamma function. For the case of = 11/3, the power spectrum reduces to the van Karman spectrum with Kolmogorov statistics, and Eq. 6 can be expressed aswhere is the coherence length of a spherical wave propagating through the turbulent medium.
In the modeling of a partially coherent vortex beam, we begin with a planar, secondary Gaussian Schell-model vortex (GSMV) source, with a source spectral degree of coherence of the formwhere is the coherence width. Substituting Eqs. 2, 14, 15 into Eq. 5, the cross-spectral density function of the GSMV beam with l = 1 in the plane z can be derived analytically. The result is given by the expressionwith
For a paraxial scalar PCV beam, the average OAM flux density can be shown to be related to the CSD function by the expression [25]where is the free-space permittivity, and represent partial derivatives with respect to and , respectively. It is to be noted that the OAM flux density at a point depends not only on the strength of rotation at the point, but also on the local intensity. To better understand the physics of the OAM distribution, we may consider a normalized OAM flux density , which represents the average OAM flux density per photon,where is the z-component of the Poynting vector. The total OAM per photon can be given by the ratio of integrated and ,
The total OAM is conserved on propagation, while the distribution of OAM flux density may change during propagation. These changes will be a combination of correlation-induced OAM changes and turbulence-induced OAM changes. Substituting Eq. 16 into Eq. 18, one finds that the average OAM flux density of a GSMV beam is given by the expression
With the help of Eqs. 19, 21, we can calculate the distribution of normalized OAM flux density for the GSMV beam propagating in atmospheric turbulence.
To illustrate the possible effects of atmospheric turbulence on the distribution of the OAM flux density of PCV beams with different spatial coherence states, we will involve multi-Gaussian Schell-model vortex (MGSMV) source, for which the source spectral degree of coherence has the formwhere
Here, is the normalization factor, is the Pochhammer symbol (falling factorial) [35], and are the coherence widths of the constituent Gaussian functions. In the far field, the MGSMV beams with represent flat-topped profiles [36], and the MGSMV beams with represent cusped profiles [37]. For , the MGSMV source reduces to the traditional GSMV source in Eq. 15.
We will also consider a source that is the incoherent superposition of GSMV beams with equal and opposite topological charges but different coherence width , for which the CSD function has the form
This source has a net zero OAM flux density but will exhibit OAM flux density changes on propagation due to the different correlation widths of the different OAM modes, even in free space [29].
With the help of Eqs. 16, 19, 21, the normalized OAM flux density of MGSMV beams and the incoherent superposition of GSMV beams can be constructed by superposition of the solutions of GSMV beams with different coherence widths .
Results
To show the effect of atmospheric turbulence on the evolution of OAM of PCV beams on propagation, we calculate the normalized OAM flux density of GSMV beams, MGSMV beams and the incoherent superposition of GSMV beams, respectively. The distribution of spectral density is also given for comparison.
Figure 1 shows the evolution of the normalized OAM flux density of GSMV beams in atmospheric turbulence, for different propagation distances . We take and for the remainder of the paper. It can be seen from Figure 1A that over a short propagation distance , the difference between the OAM distribution near the core for and is small, which is similar with that of the spectral density in Figure 1D, and with the increase of propagation distance the difference becomes larger. However, in the outskirts, even at small propagation distance, the difference between the OAM distribution is obvious, which indicates that the OAM in the outskirts is more affected by atmospheric turbulence than that near the core. This results in an interesting phenomenon that the OAM flux density is larger in atmospheric turbulence than in free space for shorter propagation distances but becomes smaller at longer distances.
FIGURE 1
Figures 2, 3 show the evolution of OAM in atmospheric turbulence for and, respectively. It can be seen that the distribution of OAM in free space mimics that of the spectral density of the respective beams: the beam () with a flat-topped spectral density results in a flat-bottomed “dead zone” in the OAM, and the beam () with a cusped spectral density results in a cusped OAM distribution near the core. As the propagation distance in atmospheric turbulence is increased, the flat-topped spectral density degenerates to a Gaussian profile (see also [38]), while it is interesting to find that the “dead zone” of OAM in the core can be maintained, within which there is no circulation. For the cusped OAM distribution in atmospheric turbulence, with the increase of propagation distance, it disappears as does the cusped spectral density distribution (see also [37]).
FIGURE 2
FIGURE 3
Even more significant variations in the distribution of OAM can be made by the incoherent superposition of GSMV beams with equal and opposite topological charge. We label the positive contribution by the subscript “+” and the negative contribution by the subscript “-”. Figures 4, 5 give the evolution of OAM with the fixed correlation width of the positive vortex beam , and width and for the negative vortex beam, respectively. Counter-rotating regions of OAM flux density are created due to the different local densities of the modes with different correlation widths, creating, for example, in Figure 4 a beam core with a negative flux density and an outer region with a positive flux density.
FIGURE 4
FIGURE 5
It can be seen from the figures that this non-trivial counter-rotating structure can be maintained on propagation in atmospheric turbulence, even though the OAM and the normalized OAM flux density in the source plane. It can also be seen that the distribution of OAM is affected more than the corresponding spectral density.
In comparing Figures 4, 5, it can be seen that the OAM distribution for is affected more than that for . This may be explained by the effective turbulence resistance of lower coherence beams: a beam with lower coherence spreads faster, and therefore turbulence-induced beam spreading is less noticeable. For the constituent GSMV beam with larger coherence width and negative topological charge (), the negative part of the counter-rotating OAM distribution becomes a region which is similar to a “dead zone” as the propagation distance in turbulence increases.
Discussion
In summary, we have studied the evolution of the OAM density flux in atmospheric turbulence for PCV beams with different coherence structures, i.e., GSMV beams, MGSMV beams, and the incoherent superposition of GSMV beams. The analytical expression of the average OAM flux density in atmospheric turbulence is derived, with which we illustrate the changes of OAM distribution on propagation. The results show that after propagating in atmospheric turbulence, the OAM distribution of PCV beams exhibits rich behavior. It is found that the OAM distribution in the outskirts is more affected by atmospheric turbulence than that near the core, which results in an interesting phenomenon that the OAM flux density is larger in atmospheric turbulence than that in free space for shorter propagation distances but becomes smaller at longer distances. Especially, for the MGSMV beams with , the “dead zone” of OAM in the core can be maintained for long propagation distance. For the incoherent superposition of GSMV beams, the non-trivial counter-rotating structure can be maintained on propagation in atmospheric turbulence, even though the OAM and the normalized OAM flux density in the source plane. Meanwhile, for the constituent GSMV beam with large coherence width and negative topological charge, with the increase of propagation distance in atmospheric turbulence, the negative part of the counter-rotating OAM distribution becomes a region which is similar as the “dead zone”. Our findings may be useful in the application of free-space OAM communication in the presence of atmospheric turbulence.
We have illustrated the possible OAM changes in air turbulence only for a very small portion of partially coherent beams; an entirely different behavior may be expected for non-Schell model beams, e.g., non-uniformly correlated, twisted beams or radially accelerating random beams. The calculation of the OAM flux is far from trivial even for general Schell-model beams, even in free space [39], and hence, better methods need to be developed for comprehensive analysis of the OAM flux density of beams with arbitrary coherence states, for propagation in both free space and turbulent media.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding author.
Author contributions
All authors listed have made a substantial, direct, and intellectual contribution to the study and approved it for publication.
Funding
The authors thank the National Natural Science Foundation of China (12174173, 12004218), the Fujian Provincial Natural Science Foundation of China (2022J02047), and the Shandong Provincial Natural Science Foundation of China (ZR2020QA067). GG was supported by the Air Force Office of Scientific Research (FA9550-21-1-0171) and Office of Naval Research (MURI N00014-20-1-2558). OK was supported by the Cooper Fellowship of the University of Miami.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1.
SoskinMSVasnetsovMV. Singular optics. Prog Opt (2001) 42:219–76. 10.1016/s0079-6638(01)80018-4
2.
DennisMRO’HolleranKPadgettMJ. Singular optics: Optical vortices and polarization singularities. Prog Opt (2009) 53:293–363.
3.
GburG. Singular optics. Boca Raton, FL: CRC Press (2017). p. 545.
4.
LeeJHFooGJohnsonEGSwartzlanderGAJr. Experimental verification of an optical vortex coronagraph. Phys Rev Lett (2006) 97:053901. 10.1103/physrevlett.97.053901
5.
PalaciosDRozasDSwartzlanderGAJr. Observed scattering into a dark optical vortex core. Phys Rev Lett (2002) 88:103902. 10.1103/physrevlett.88.103902
6.
GibsonGCourtialJPadgettMJVasnetsovMPas’koVBarnettSMet alFree-space information transfer using light beams carrying orbital angular momentum. Opt Express (2004) 12:5448–56. 10.1364/opex.12.005448
7.
AllenLBeijersbergenMWSpreeuwRJCWoerdmanJP. Orbital angular momentum of light and the transformation of Laguerre-Gaussian laser modes. Phys Rev A (1992) 45:8185–9. 10.1103/physreva.45.8185
8.
GrierDG. A revolution in optical manipulation. Nature (2003) 424:810–6. 10.1038/nature01935
9.
SimpsonNBDholakiaKAllenLPadgettMJ. Mechanical equivalence of spin and orbital angular momentum of light: An optical spanner. Opt Lett (1997) 22:52–4. 10.1364/ol.22.000052
10.
LadavacKGrierDG. Microoptomechanical pumps assembled and driven by holographic optical vortex arrays. Opt Express (2004) 12:1144–9. 10.1364/opex.12.001144
11.
Ritsch-MarteM. Orbital angular momentum light in microscopy. Phil Trans R Soc A (2017) 375:20150437. 10.1098/rsta.2015.0437
12.
RuffatoGRossiRMassariMMafakheriECapaldoPRomanatoF. Design, fabrication and characterization of Computer Generated Holograms for anti-counterfeiting applications using OAM beams as light decoders. Sci Rep (2017) 7:18011. 10.1038/s41598-017-18147-7
13.
ErhardMFicklerRKrennMZeilingerA. Twisted photons: New quantum perspectives in high dimensions. Light Sci Appl (2018) 7:17146. 10.1038/lsa.2017.146
14.
WangLMaJXiaoMZhangY. Application of optical orbital angular momentum to rotation measurements. Results Opt (2021) 5:100158. 10.1016/j.rio.2021.100158
15.
RicklinJCDavidsonFM. Atmospheric turbulence effects on a partially coherent Gaussian beam: Implications for free-space laser communication. J Opt Soc Am A (2002) 19:1794–802. 10.1364/josaa.19.001794
16.
ZhaoCCaiYLuXEyyuboğluHT. Radiation force of coherent and partially coherent flat-topped beams on a Rayleigh particle. Opt Express (2009) 17:1753–65. 10.1364/oe.17.001753
17.
ZhangJWangZChengBWangQWuBShenXet alAtom cooling by partially spatially coherent lasers. Phys Rev A (2013) 88:023416. 10.1103/physreva.88.023416
18.
SchoutenHFGburGVisserTDWolfE. Phase singularities of the coherence functions in Young’s interference pattern. Opt Lett (2003) 28:968–70. 10.1364/ol.28.000968
19.
GburGVisserTD. Coherence vortices in partially coherent beams. Opt Commun (2003) 222:117–25. 10.1016/s0030-4018(03)01606-7
20.
GburGVisserTDWolfE. ‘hidden’ singularities in partially coherent wavefields. J Opt A: Pure Appl Opt (2004) 6:S239–42. 10.1088/1464-4258/6/5/017
21.
WangWTakedaM. Coherence current, coherence vortex, and the conservation law of coherence. Phys Rev Lett (2006) 96:223904. 10.1103/physrevlett.96.223904
22.
RaghunathanSBSchoutenHFVisserTD. Correlation singularities in partially coherent electromagnetic beams. Opt Lett (2012) 37:4179–81. 10.1364/ol.37.004179
23.
StahlCSDGburG. Complete representation of a correlation singularity in a partially coherent beam. Opt Lett (2014) 39:5985–8. 10.1364/ol.39.005985
24.
PalaciosDMMaleevIDMarathayASSwartzlanderGAJr. Spatial correlation singularity of a vortex field. Phys Rev Lett (2004) 92:143905. 10.1103/physrevlett.92.143905
25.
KimSMGburG. Angular momentum conservation in partially coherent wave fields. Phys Rev A (2012) 86:043814. 10.1103/physreva.86.043814
26.
GburG. Partially coherent vortex beams (Invited). Proc SPIE (2018) 10549:1054903.
27.
StahlCSDGburG. Partially coherent vortex beams of arbitrary order. J Opt Soc Am A (2017) 34:1793–9. 10.1364/josaa.34.001793
28.
ZhangYCaiYGburG. Control of orbital angular momentum with partially coherent vortex beams. Opt Lett (2019) 44:3617–20. 10.1364/ol.44.003617
29.
ZhangYKorotkovaOCaiYGburG. Correlation-induced orbital angular momentum changes. Phys Rev A (2020) 102:063513. 10.1103/physreva.102.063513
30.
KorotkovaO. Theoretical statistical optics. London: World Scientific (2021). p. 280.
31.
WolfE. Introduction to the theory of coherence and polarization of light. New York: Cambridge University Press (2007). p. 222.
32.
AndrewsLCPhillipsRL. Laser beam propagation through random media. Bellingham, Washington: SPIE Press (1998). p. 775.
33.
ShiraiTDogariuAWolfE. Mode analysis of spreading of partially coherent beams propagating through atmospheric turbulence. J Opt Soc Am A (2003) 20:1094–102. 10.1364/josaa.20.001094
34.
ShchepakinaEKorotkovaO. Second-order statistics of stochastic electromagnetic beams propagating through non-Kolmogorov turbulence. Opt Express (2010) 18:10650–8. 10.1364/oe.18.010650
35.
KorotkovaOHydeI. Multi-Gaussian random variables for modeling optical phenomena. Opt Express (2021) 29:25771–99. 10.1364/oe.432227
36.
SahinSKorotkovaO. Light sources generating far fields with tunable flat profiles. Opt Lett (2012) 37:2970–2. 10.1364/ol.37.002970
37.
WangFKorotkovaO. Circularly symmetric cusped random beams in free space and atmospheric turbulence. Opt Express (2017) 25:5057. 10.1364/oe.25.005057
38.
KorotkovaOSahinSShchepakinaE. Multi-Gaussian schell-model beams. J Opt Soc Am A (2012) 29:2159–64. 10.1364/josaa.29.002159
39.
WangHYangZLiuLChenYWangFCaiY. Fast calculation of orbital angular momentum flux density of partially coherent Schell-model beams on propagation. Opt Express (2022) 30:16856–72. 10.1364/oe.459089
Summary
Keywords
atmospheric turbulence, propagation, partially coherent, orbital angular momentum, vortex
Citation
Zhang Y, Yu J, Gbur G and Korotkova O (2022) Evolution of the orbital angular momentum flux density of partially coherent vortex beams in atmospheric turbulence. Front. Phys. 10:1073662. doi: 10.3389/fphy.2022.1073662
Received
18 October 2022
Accepted
14 November 2022
Published
30 November 2022
Volume
10 - 2022
Edited by
Carmelo Rosales-Guzmán, Centro de Investigaciones en Optica, Mexico
Reviewed by
Gianluca Ruffato, University of Padua, Italy
Chengliang Zhao, Soochow University, China
Updates
Copyright
© 2022 Zhang, Yu, Gbur and Korotkova.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Olga Korotkova, korotkova@physics.miami.edu
This article was submitted to Optics and Photonics, a section of the journal Frontiers in Physics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.