Abstract
In the aero-optic turbulent boundary layer (TBL), there exist very rich air flow structures that fall into a wide range of scales, with the smallest being roughly of the order of the optical wavelength. However, these fine spatial variations cannot be neglected when one is dealing with light propagation through such structures, since both the amplitude and phase of a light wave undergo modulations. In this study, we studied the influence of TBL on the angular momentum spectrum of light and found that there exists critical point of the azimuthal distribution of the disturbance phase that determines the symmetric properties of the expansion spectrum.
1 Introduction
Allen discovered that Laguerre–Gaussian beams carry orbital angular momentum (OAM) [1], which is associated with their helical wave front. Since then, optical beams with OAM and helical phase structures have garnered increasing research interest. Owing to the infinite number of eigenmodes, beams of this type have found applications in fields such as laser and quantum communications [2–4]. Intense efforts have been put into the measurement of OAM, especially regarding the evolution and the spiral spectrum during light propagation.
For quite a long time, it was believed that there exists a close interconnection between the spiral spectrum and phase vortex structures. However, Charnotskii’s recent research showed that such a connection does not exist [5]. A tricky scenario is propagation in inhomogeneous media. Even if the incident wave consists of only one eigenmode, its spiral spectrum can spread during propagation [6–8], a phenomenon called cross-talk. Light with OAM propagating through inhomogeneous media is an important topic that is employed across many applications.
V. P. Aksenov and Ch. E. Pogutsa gave the mathematical form of the total OAM of a propagating light field, from which they were able to prove that the total OAM for homogeneous media is conserved, while it can change for inhomogeneous media [9]. Such a change is closely related to the spatial distribution of the refractive index (RI).
According to the traditional turbulence theory, provided the fluid is far away from the boundary and has a very high Reynolds number, the physical parameters of the turbulence can be treated as random fields. The properties of turbulence can be described through the statistical moment of fluid parameters [10]. Consequently, by introducing statistical theory and assuming the media to be statistically homogeneous, it is deduced that the ensemble average of total OAM is conserved [5, 9], while the spiral spectrum demonstrates a symmetric distribution with respect to the original [6–8]. Darryl J. Sanchez and Denis W. Oesch tentatively discussed the physical mechanism of the phenomenon: the kinetic energy of free turbulence dissipates as heat, which makes the spatial distribution of RI in turbulence inhomogeneous. According to the cascade theory of turbulence, the RI distribution is statistically homogeneous, so that photons carrying the new OAM must appear in pairs [11].
With the development of fluid mechanics and the progress of numerical techniques, researchers now have a new understanding of fluid motion, especially for fluids close to a definite solid boundary, for example, the turbulence around an airborne plane. Such regions are called the turbulent boundary layer (TBL). E. Jumper [12] pointed out that the power spectrum model cannot be used to describe the characteristics of turbulence in TBL because the model has definite structural properties [13, 14] and is not completely random. Research in this field has evolved to become the so-called aero-optics [15]. The properties of the total OAM and spiral spectrum distribution in TBL are our main focus in this study. In our task, the transmission medium was deterministic and inhomogeneous. Charnotskii calculated the total OAM in turbulence, and his result showed that . This implies that for a single realization of turbulence, the total OAM usually varies. Lingfei Xu [16] calculated the total OAM and the spiral spectrum of OAM modes in a typical TBL by numerical methods and found that the instantaneous total OAM is not conserved and the spiral spectrum of OAM spreads in an asymmetric way. These behaviors are related to the typical structure of RI distribution in TBL.
Although numerical methods can present the evolution characteristics of OAM, it is still challenging to discuss the connection between OAM and the RI spatial distribution. The lack of an analytic form of the spatial structure makes it difficult to study the evolution mechanism. In this study, we do not rely on the concrete form and try to analyze the problem analytically. Hence, our idea is to obtain the asymptotic form of the spiral spectrum of OAM modes. Tow case is concerned: the fluctuation of refractive index is weak or much stronger. By such a mathematical method, we found that the symmetry of the spiral spectrum depends on the distribution of critical points in the scattered light’s phase. Based on such mathematical results, we discuss the mechanism of the influence and propose a method for analyzing the spiral spectrum spread. Our main aim was to analyze the asymptotic behavior of the spiral spectrum.
2 Asymptotic Behavior of the Scattered Light’s Spiral Spectrum
The complex amplitude of light propagating in the turbulent boundary layer can be obtained by the Rytov approximation:where , denotes the diffracted complex amplitude, and denotes the Rytov disturbance. In general, the Rytov disturbance can be written as follows:where is the wavenumber of light. and denote the disturbance on the amplitude and phase, respectively, and the coefficient for each spiral spectrum of OAM modes is given by the following equation: [16]where we write the complex amplitude of LG beam as follows: . denotes the difference between the spiral spectrum of OAM modes and initial beam’s OAM mode. In order to analyze the asymptotic of Eq. 3, we divide the analysis into two cases: 1), the disturbance phase that is much small; and2) the disturbance phase that is very large.
2.1 Asymptotic Behavior of CΔl for
In the weak fluctuating turbulence, the disturbance on the amplitude of beam can be ignored since the turbulent energy at high spatial frequency is really small. So, can be approximately equal to 1. is also small, hence it will be a good approximation if the Rytov complex disturbance is written as follows:substituting Eq. 4 into Eq.3, we can obtainwhere is the Dirac function. Because is a real function, from Eq.5, the following conclusion can be immediately obtained: .
2.2 Asymptotic Behavior for Large P
In this case, the distribution properties of the spiral spectrum are much more complex. It is difficult to analyze and examine the symmetry of the spiral spectrum distribution since the mathematical form of P is hardly obtained in the turbulence boundary layer. Hence, instead of obtaining the precise solution of Eq.3, we tend to discuss the spiral spectrum of OAM modes based on the asymptotic form obtained from Eq. 3. The key point is to discuss the asymptotic behavior of the integral:We can use the stationary phase method [17]. By using the stationary phase, we can obtainwhere , the symbol depends on . is the critical point of the first kind: . Here, we used the fact that and . Thus, the contributions from the critical point of the second kind ( and ) is zero. Hence, there iswhere and . are the subscripts of such critical points of the first kind whose and are the subscripts of the critical points of the first kind whose .
We will discuss about the case wherein the number of the first kind critical points is two or more. Equation 8 indicates that, in general, the spiral spectrum of the OAM modes is asymmetric. It will be symmetric if the terms of in Eq. 8 vanish, for example, (i), the azimuthal difference between each of the two critical points of first kind is , where ; or,(ii) the phase difference between each of the two first kind critical points is equal to for and for .
Equation 8 gives us a way to analyze the symmetric properties of the scattered light’s spiral spectrum. If the scale of coherent structures in the turbulent boundary layer is much larger than the wavelength of light, the disturbance on the phase is approximately equal to the integral of the refractive index difference: , where denotes the difference between the RI in the turbulent boundary and in the free space. can be treated as 1 approximately. Then, by Eq. 8, the symmetric properties can be analyzed by examining the critical points of RI distribution even when the exact scattered complex amplitude is not solved.
3 Verification With a Well-Known Case: Phase Tilt
The influence on the light’s OAM and the spiral spectrum of OAM modes is widely studied in the case in which the phase is tilted. The conclusion is well known: the spiral spectrum of OAM modes is expanded in a symmetric way. Here, we assume that the additional phase is tilted along with one direction, said, , where, a is the tilt degree. Hence, the spiral spectrum of OAM modes is given by the following equation:where is the -th order Bessel function. Equation 9 shows that the spiral spectrum of OAM modes is symmetric with the help of . The symptomatic behavior of Bessel function is given by the following equation:Substituting Eq. 10 and Eq. 11 into Eq. 9, the symptomatic behavior of spiral spectrum can be obtained as follows:Equations 12, 13 give the asymptotic behavior of spiral spectrum from its exact mathematical form (Eq. 9). Now, we will examine the asymptotic behavior of the spiral spectrum from Eq.5 to Eq. 8. By Eq.5, for a small a,Equation 14 shows the following:which coincides with Eq. 12 when . is extremely small for a large value of , and we can say that two kinds of asymptomatic behavior are same in the case of small a.For a large a, according to the method of the stationary phase, we have to calculate critical point of the first kind of . It is extremely easy to find two critical points and . Additionally, , , and . Substituting into Eq. 8, we get the following:therefore, the coefficient of the spiral spectrum is given as follows:which also coincides with Eq. 13. We also give the comparisons between the theoretical results of Eq. 9 and the asymptomatic results, (Figure 1). By comparison between the asymptomatic results and the numerical simulation result, it can be found that the asymptomatic method can describe the behavior of the spiral spectrum with a minor error, which is about . For the small , the errors will be smaller.
FIGURE 1
4 Symmetric Property of Spiral Spectrum in Some Simple Structured Media
Vortex is the basic structures in the turbulence. The refractive index of a vortex flow field can be treated as a Gaussion distribution. In order to avoid solving Navier–Stokes equations, we prefer to use such a simple vortex structure to discuss the symmetric property of the spiral spectrum by the method proposed in Section 2. According to Section 2.1, if the disturbance phase is small, the spiral spectrum of the scattered light is always symmetric. Hence, in this section, the symmetric property of light with strong phase disturbance is discussed. Moreover, the disturbance on the amplitude is ignored, assuming that the scale of vortex is much larger than that of the wavelength. Of course, it is immediately obtained by Eq. 3 that if the vortex flow field is on axis, the spiral spectrum will not expand at all. Hence, we choose a single vortex structure, the center of which is off-axis. The disturbance phase P can be written as follows:The critical points of the first kind of Eq. 18 are and . Thus, the spiral spectrum of the OAM modes is, of course, symmetric. The normalized RI distribution and spiral spectrum are shown in Figure 2; Table 1 gives the numerical results for from to . An interesting example is given by two off-axis vortex structures. The phase disturbance can be written as follows:Such structures are composed by two off-axis Gaussian vortex structure. The center of one vortex is at and another is at . The phase distribution is shown as Figure 3. Following the proposed method in Section 2.2, the critical point of such distribution can be obtained as: , , , and . There are four critical points. Calculating P and for each critical point, we have the following:Hence, the asymptotic behavior of Eq. 6 is as follows:where and . From Eq. 21, the distribution of the spiral spectrum of the OAM modes is symmetric. Equation 21 also demonstrates that if is odd, will be vanished at . The spiral spectrum can be directly calculated by Eq. 3, shown as Figure 4; Table 2 gives the precise numerical results. It expands in a symmetric way and vanished at the position where is odd. Such conclusions coincide with those obtained by the proposed method in Section 2.2. So far, we can conclude that the symmetric property is related to the critical point of the phase distribution. In order to make this conclusion more credible, a special phase is constructed according to the rule of symmetry affected by the critical point proposed in Section 2.2:The critical points can be calculated at , , and . After tedious calculation, we obtain the following:Under the integral for ρ, as shown in Eq.3, the terms will not vanish in general. Thus, the spiral spectrum appears asymmetric. The phase distribution and spiral spectrum calculated by Eq. 3 are shown in Figure 5; Table 3.
FIGURE 2
TABLE 1
| 0 | |||||||
| 0.014 | 0.019 | 0.0192 | 0.85 | 0.023 | 0.019 | 0.014 | |
Coefficients of spiral spectrum modes for single Gaussion vortex structure.
FIGURE 3
FIGURE 4
TABLE 2
| 0 | |||||||||
| 0.019 | 0.0001 | 0.11 | 0.59 | 0.11 | 0.019 | ||||
Coefficients of spiral spectrum modes for two off-axis Gaussian vortex structures.
FIGURE 5
TABLE 3
| 0 | |||||||
| 0.069 | 0.13 | 0.19 | 0.12 | 0.048 | 0.052 | 0.033 | |
Coefficients of spiral spectrum modes for an astigmatism phase.
The critical point distribution makes the spiral spectrum asymmetric. In this case, such critical point distribution is mainly affected by the term , which is known as astigmatism.
5 Conclusion
We studied the symmetric property of light’s spiral spectrum propagating through a non-random inhomogeneous media. Without losing the generality, instead of using the calculated phase aberration, we treated the unknown continuous distribution function of phase disturbance in a formal way. Through asymptotic methods, we found that there exists a critical point that determines the symmetric property of the spiral spectrum. The position and refractive index, along with its second derivative, determine the spectrum distribution.
Statements
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author contributions
LX conceived the research and derived equations about the asymptotic form; JH and ZZ wrote the computational codes, computed the numerical results, and made the figure plots; TR organized the article; and the four authors wrote the article.
Acknowledgments
The authors thank Cunfeng Gu in our department for their helpful discussions.
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.
References
1.
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
2.
GibsonGCourtialJPadgettMJVasnetsovMPas'koVBarnettSMet alFree-space information transfer using light beams carrying orbital angular momentum. Opt Express (2004) 22:5448–56. 10.1364/opex.12.005448
3.
MajumdarARiclkinJC. Free-space laser communications: principles and advances. Berlin, Germany: Springer (2007).
4.
WangJYYangIMFazalNAYanYHuangHRenYet alTerabit free-space data transmission employing orbital angular momentum multiplexing. Nat Photon (2012) 6:488–96. 10.1038/nphoton.2012.138
5.
CharnotskiiM. Transverse linear and orbital angular momenta of beam waves and propagation in random media. J Opt (2017) 20:025602. 10.1088/2040-8986/aa9f50
6.
ChengMGuoLLiJHuangQChengQ. Propagation of an optical vortex carried by a partially coherent Laguerre-Gaussian beam in turbulent ocean. Appl Opt (2016) 55:4642–8. 10.1364/ao.55.004642
7.
YanXZhangPFZhangJHChunHQFanCY. Decoherence of orbital angular momentum tangled photons in non-kolmogorov turbulence. J Opt Soc Am A (2016) 33:1831–5. 10.1364/josaa.33.001831
8.
ChenCYangHTongSLouY. Changes in orbital-angular-momentum modes of a propagated vortex Gaussian beam through weak-to-strong atmospheric turbulence. Opt Express (2016) 24:6959–75. 10.1364/oe.24.006959
9.
AksenovVPPogutsaCE. Fluctuations of the orbital angular momentum of a laser beam, carrying an optical vortex, in the turbulent atmosphere. Quan Electron. (2008) 45:343–8. 10.1070/qe2008v038n04abeh013576
10.
SanchezDJOeschDW. Localization of angular momentum in optical waves propagating through turbulence. Opt Express (2011a) 19:25388–96. 10.1364/oe.19.025388
11.
SanchezDJOeschDW. Orbital angular momentum in optical waves propagating through distributed turbulence. Opt Express (2011b) 19:24596–608. 10.1364/oe.19.024596
12.
SiegenthalerJPJumperEJ. Atmospheric propagation vs. aero-optics. Reston, VA: AIAA (2008).
13.
JumperEJGordeyevS. Physics and measurement of aero-optical effects: past and present. Annu Rev Fluid Mech (2017) 49:419–41. 10.1146/annurev-fluid-010816-060315
14.
ElsingaGEvan OudheusdenBWScaranoF. Evaluation of aero-optical distortion effects in piv. Exp Fluids (2005) 39:246–56. 10.1007/s00348-005-1002-8
15.
KirmseTGardnerADKrombholzC. Density effects in model deformation measurements on a naca 0010 profile in a transonic flow. In: 14th international symposium on flow visualization; 2010 Jun 21–24; Daegu, Korea (2010).
16.
XuLZhouZRenTYuXHanB. Propagation characteristics of orbital angular momentum and its time evolution carried by a laguerre-Gaussian beam in supersonic turbulent boundary layer. Opt Express (2020) 28:4032–47. 10.1364/oe.382421
17.
WolfEBornM. Principles of optics - electromagnetic theory of propagation, interference and diffraction of light. 7th ed. New York, NY: Cambridge University Press (1999).
Summary
Keywords
orbital angular momentum, spiral spectrum of OAM modes, symmetric property, critical point of phase, turbulence boundary layer
Citation
Xu L, Hong J, Zhou Z and Ren T (2021) Study on the Relationship Between Symmetric Property of Optical Orbital Angular Momentum Spectrum and Critical Points of Phase in an Inhomogeneous Non-Random Media. Front. Phys. 9:653849. doi: 10.3389/fphy.2021.653849
Received
15 January 2021
Accepted
10 February 2021
Published
13 April 2021
Volume
9 - 2021
Edited by
Hao Chen, Shanghai Institute of Microsystem and Information Technology (CAS), China
Reviewed by
Baifu Zhang, Nanjing University of Science and Technology, China
Guoquan Zhou, Zhejiang Agriculture and Forestry University, China
Rui-Pin Chen, Zhejiang Sci-Tech University, China
Updates
Copyright
© 2021 Xu, Hong, Zhou and Ren.
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: Lingfei Xu, jsj900@hotmail.com
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.