# Intensity Distribution of Partially Coherent Array Finite Airy Beams Propagating in Atmospheric Turbulence

Authors: Long Jin; Zirui Zhang; Nuo Wang; Zixin Liu; Yuwei Deng; Yanhua Fu
Journal: International Journal of Optics
 (2021)

---

## Abstract

Based on the extended Huygens–Fresnel integral and the Rytov phase structure function, the closed-form expression of a partially coherent array finite Airy beams (PCAFABs) cutting through the Kolmogorov atmospheric turbulence is derived in the space domain under the paraxial approximation. The characteristics of the PCAFABs evoluting in the atmospheric environment are investigated in detail on the basis of the derived wave propagation formulae. We mainly illustrate the intensity profile of this beam changed with the truncation parameter, coherence length, and turbulence factor at several cross sections of the atmospheric space by means of numerical figures. It is convinced that the present concept and derived conclusions will provide useful exploration for learning the optical properties of the PCAFABs transmitting in the atmospheric turbulence, especially for free-space optical communication area.


---

## Body

## 1. Introduction

Atmospheric laser communication is featured with the advantages of large channel capacity, good confidentiality, convenient installation, and low production cost compared with other wireless and wired communication equipment [1]. However, the diffraction and non-self-heal effects of Gaussian beam are very obvious in the transmission process when it meets obstacles or blocks [2]. Moreover, it is difficult to overcome the defects such as low energy utilization rate, large spot divergence angle, slow transmission rate, and high bit error rate for Gaussian beam even with the adoption of partial coherence technology [3]. Therefore, it is becoming the research hot spot for science and engineering experts to find a kind of novel laser light source itself owning the ability to resist atmospheric turbulence in recent decades.The Airy beam with finite energy was experimentally proposed by Siviloglou et al. in 2007, and they mainly used the computer-controlled liquid crystal spatial light modulator to generate the Airy beam model coupled with an exponential term [4]. From then on, it is believed that the self-accelerating, nondiffracting, and self-healing properties of the Airy beam will be of great value in long-range, high-speed atmospheric communication. In 2011, Chu examined the centroid position and skewness of Airy beam intensity in atmospheric turbulence [5]. In 2013, Eyyuboğlu et al. discussed the propagation characteristics of both Airy beam and partially coherent Airy beam in turbulent media, respectively [6, 7]. In 2014, Yang et al. studied the degree of polarization of vector partially coherent Airy beams propagating through atmospheric turbulence [8]. In 2019, Li et al. explored the intensity and phase distributions of a ring Airy Gaussian vortex beam propagating in atmospheric turbulence [9]. However, not only is the amount of energy carried by a single-beam source extremely limited, but also the asymmetric transverse intensity property of the single Airy beam is adverse in some optical application [10], so it is necessary to arrange the multiple Airy beams in the radial symmetry, which is called the array Airy beam [11]. In the past few years, research on the array Airy beam cultivated a new interdisciplinary area. In 2010, Gu et al. investigated the scintillation properties of Airy beam arrays in atmospheric turbulence, and they found that the scintillation of an Airy beam array is significantly reduced and close to the theoretical minimum [12]. In 2014, Chen et al. derived the average intensity of both phase-locked and non-phase-locked Airy array beams propagating through atmospheric turbulence, and they obtained that the differences between the propagation properties of phase-locked Airy array beams and those of non-phase-locked ones become unobservable when the atmospheric turbulence is strong enough [13]. In 2017, Lu et al. investigated the generation of coherent and incoherent Airy beam arrays and the comparisons of their scintillation characteristics in atmospheric turbulence, and they pointed out that the incoherent Airy beam arrays have a smaller scintillation index than the coherent one in the same turbulent condition due to the coherence reduction of the constituent beamlets [14]. In 2020, Suarez et al. presented the experimental generation and analysis of both the intensity and phase of an optical vortex beam originating from a superposition of Airy-vortex beams [15]. However, as far as we know, there are few studies on the evolution rules of the partially coherent array finite Airy beams (PCAFABs) in atmospheric turbulence. We are convinced that the rational use of the evolution properties of PCAFABs in atmospheric turbulence will provide a powerful and convenient way for free-space optical communication.In this paper, we will research systematically the PCAFAB intensity outlines at some cross sections of the atmospheric environment, and the beam transverse self-acceleration parabolic trajectory with the changing of the transmission distance. The rest of the paper is organized as follows: in Section2, the general intensity expression of the PCAFABs passing through the Kolmogorov atmospheric turbulence is derived in the space domain by the extended Huygens–Fresnel integral under the paraxial approximation. In Section 3, the intensity distribution of the PCAFABs changed with the truncation parameter, coherence length, and turbulence factor are explored by numerical examples in detail. In Section 4, a clear, concise summary is outlined at the end of our theses.


## 2. Mathematical Model and Calculation Method

The electric field of the two-dimensional single finite Airy beam at source plane (z = 0) can be described as [16](1)E1r′,z=0=E1x0,0·E1y0,0=Aix0ω0,×expa0x0ω0·Aiy0ω0×expa0y0ω0,where the vector r′=x0,y0 denotes the transverse coordinates at two axial directions, ω0 characters the incident beam waist size, 0 <a0 <1 provides the truncation parameter, and the symbol Ai expresses the Airy function, of which the standard integral function is [17](2)Aix=12π×∫−∞+∞expiu33+iuxdu.The symmetric array finite Airy beams (AFABs) can be constituted through superimposing four parts of single Airy beams. By adjusting each beamlet’s phase parameter, the incident AFABs are given by [18](3)E2r′,z=0=E1x0,0×E1y0,0+E1−x0,0E1y0,0+E1x0,0×E1−y0,0+E1−x0,0×E1−y0,0.The cross-spectral density function (CSDF) is usually used to describe a partially coherent laser beam. Here, we select the Gaussian Schell model, leading to the expression of the incident PCAFABs as [19](4)Wr1′,r2′,0=E2r1′,0×E2∗r2′,0×exp−r1′−r2′2σ2,where σ represents the initial coherence length. The CSDF of PCAFABs evoluting in the atmospheric turbulence can be evaluated by the extended Huygens–Fresnel integral. Under the paraxial approximation, the integral expressions for this wave propagation at arbitrary z-plane are given by [20](5)Wr1,r2,z=k2πz2∫−∞+∞∫−∞+∞∫−∞+∞∫−∞+∞expik2zr1−r1′2−r2−r2′2×Wr1′,r2′,0×<expφr1′,r1+φ∗r2′,r2>mdr1′dr2′,where W (r1, r2, z) and W (r1’, r2’, z) are the CDSF of PCAFABs at source and the output planes, respectively, k = 2π/λ is the wave number with being the λ wavelength, r1 and r2 are two different coordinate vectors, respectively, z is the propagation distance, φ (r’, r) denotes the complex random phase fluctuation of wave transmission in atmospheric turbulence, and <.>m represents the ensemble average and can be written approximately by the Rytov phase structure function [21]:(6)<expφr1′,r1+φ∗r2′,r2>m=exp−r1′−r2′2+r1′−r2′r1−r2+r1−r22ρ02=exp−x1′−x2′2+y1′−y2′2ρ02×exp−x1−x22+y1−y22ρ02×exp−x1′−x2′x1−x2+y1′−y2′y1−y2ρ02,where ρ0=0.545Cn2k2z−3/5 is known to be the coherence length and Cn2 is the structure constant of the refractive index fluctuations of the turbulence. Because the PCAFABs are symmetric in x and y directions, the W (r1, r2, z) can be separated with each other; that is to say,(7)Wr1,r2,z=Wx1,x2,z×Wy1,y2,z.Demanding thatx1 = x2 = x, y1 = y2 = y, the intensity expression for the partially coherent beam can be written by(8)Ix,z=Wx,x,z,Iy,z=Wy,y,z.Finally, the PCAFAB intensity distribution at the receiver plane is obtained as(9)Ir,z=Ix,z×Iy,z+I−x,z×Iy,z+Ix,z×I−y,z+I−x,z×I−y,z,where(10)Ix,z=Aixω0−z24kω04−4a0b0+b02+za0+b0kω02i×Aixω0−z24kω04−4a0b0+b02−za0+b0kω02i×exp2a0+b0xω0−z22k2ω04−4a0b0+b02+3a0b02+b03,(11)Iy,z=Aiyω0−z24kω04−4a0b0+b02+za0+b0kω02i×Aiyω0−z24kω04−4a0b0+b02−za0+b0kω02i×exp2a0+b0yω0−z22k2ω04−4a0b0+b02+3a0b02+b03,with(12)b0=z2ck2ω02,c0=12σ2+1ρ02.Other optical field intensity values in equation (9) can be achieved through changing the plus or minus signs of the transversal coordinates x and y in equations (10) and (11). Each beamlet of the PCAFABs in atmospheric turbulence follows a ballistic trajectory in the x-z and y-z planes; for example, the central beam location with the changes z in the first quadrant is described by the following parabolas:(13)x=y=z24kω04−4a0b0+b02×ω0.Other three beamlets’ trajectory coordinates can be acquired by changing the plus or minus signs of the transversal coordinatesx and y in equation (13).


## 3. Graphical Results and Analysis

In this section, as numerical examples, we investigate the influence of beam parameters and turbulence factors on the PCAFABs’ evolution characteristics at different cross sections of atmospheric environment. Firstly, the variations of beam intensity at source plane and beam transverse acceleration parabolic trajectory witha0 are demonstrated in Figure 1, and other parameters are listed as follows: λ = 1.55um, σ = 100 mm, ω0 = 0.2 m, and Cn2 = 10−11m−2/3. It shows that the side lobes of incident PCAFABs gradually diminish with the increasing a0, and the beam evolves into a partially coherent array Gaussian beams (PCAGBs) when a0 = 0.5, as learned in Figure 1(c). Moreover, one can observe from Figure 1(d) that three lines of beam transverse acceleration parabolic trajectory overlap each other, manifesting that the central main lobe location is rarely affected by the changes of a0.Figure 1
Influence ofa0 on PCAFABs’ evolution characteristics; other parameters are fixed as λ = 1.55um, σ = 100 mm, ω0 = 0.2 m, and Cn2 = 10−13m−2/3: (a) incident PCAFABs, a0 = 0.01; (b) incident PCAFABs, a0 = 0.1; (c) incident PCAFABs, a0 = 0.5; (d) beam transverse acceleration parabolic trajectory for three a0.
(a)(b)(c)(d)The influence ofσ on incident PCAFABs’ silhouette is depicted in Figure 2, and other parameters are the same as those in Figure 1 except for a0 = 0.1 fixed. It is gotten that the initial coherence length plays an insignificant role on incident beam intensity distribution by comparing Figure 2(a) and Figure 2(b). As the transverse distance increases, the exponential term(14)expa0r1′+r2′ω0≫exp−r1′−r2′2σ2.Figure 2
Incident PCAFABs changed with differentσ; other parameters are the same as those in Figure 1 except for a0 = 0.1: (a) σ = 0.5 mm; (b) σ⟶∞.
(a)(b)Thus, two drawings below are much of a muchness.In order to better understand the effect of the coherence length on beam’s evolution characteristics, Figure3 describes the contour graphs of PCAFABs with different σ at several intersecting surfaces during evolution in atmospheric turbulence, and other parameters are the same as those in Figure 1 except for a0 = 0 fixed. It is recognized from these figures that the side loses of PCAFABs fade away with the σ decreases and the transmission distance increases. There are great deals of side lobs spreading over each coordinate axis when σ = 1 mm and z = 100 m, as shown in Figure 3(c); however, by viewing Figure 3(d), we discover that lots of these side lobes disappear when beam transmits to z = 500m; hence, we elucidate that the larger the σ can result in the stronger antiturbulence ability of the PCAFABs when it propagates in atmospheric environment. Meanwhile, we find that the beam’s self-bending property is increasingly obvious along with the changing of σ during the wave transmission process by viewing in Figure 3(f) and Figure 3(g), but these lines are nonmonotonic versus the increases of distance z. A mass of simulations manifest that the PCAFABs covert into the PCAGBs at the first maximum point (inflexion A); hence, equation (9) is no longer valid to describe the intensity of PCAFABs when the transmission distance z is larger than the inflexion A. Observed in Figure 3(h), it is indicated that the bigger the σ is, the more severe the bending is formed at the receiver plane, and the inflexion point A increases with the increasing of σ. As we know, the PCAFABs convert into a AFABs when σ⟶∞ and the beam self-bending acceleration at this point is the fastest than that of the other σ. However, on the other hand, the PCAFABs can result in better directive performance than that of the AFABs while transmission in the same atmospheric circumstances. Hence, we can regulate the path of the PCAFABs in atmospheric environment by choosing the different value of the σ to meet tangible needs.Figure 3
Influence ofσ on PCAFABs’ evolution characteristics; A expresses the inflexion point, and other parameters are the same as those in Figure 1 except for a0 = 0: (a) σ = 0.1 mm, (z) = 100m; (b) σ = 0.3 mm, z = 100m; (c) σ = 1 mm, z = 100m; (d) σ = 1 mm, z = 500m; (e) σ = 10 mm, z = 500m; (f) parabolic trajectory for σ = 0.3 mm; (g) parabolic trajectory for σ = 1 mm; (h) beam transverse self-acceleration parabolic trajectory for different σ.
(a)(b)(c)(d)(e)(f)(g)(h)Next, we delve into the effect ofCn2 on the PCAFABs’ deviation at some observation planes, the results are delineated in Figure 4, and we set other parameters to be the same as those in Figure 1(b). The acceleration direction of each beamlet is labeled by the blue arrow head, as shown in Figure 4(a). By viewing the beam intensity evolution rules at the circumstances Cn2 = 10−13m−2/3 in Figure 4(b) and Figure 4(c), it is found that each beamlet of the PCFAABs gradually converges to the original central location at the beginning of the propagation process until the smallest main and side spots are formed. Then, they deviate farther and farther away from the original central location with the increasing distance z in the rest of atmospheric turbulence. When z = 500m in Figure 4(c), the PCAFABs separate into four independent portions because of the self-bending peculiarity, which confirms the conclusion obtained in [22]. Moreover, we find that the side lobes of PCAFABs gradually fade away during transmission in atmospheric environment, and the larger the Cn2 is, the shorter the transmission distance z is required to achieve the disappearance of beam side petals. While Cn2 = 10−11m−2/3, the transmission distance is z = 500 where the PCAFABs gradually degrade into the PCAGBs, as can be seen in Figure 4(f). The singular properties of the PCAFABs cease to exist provided that its side lobes disappear entirely; therefore, while Cn2 = 10−8m−2/3 in Figure 4(j), the PCAFABs will never separate from each other in the rest of the atmospheric turbulence, and the spot size of these array beams increases in proportion to the Gaussian function with the increasing transmission distance until the power eventually burns out.Figure 4
Influence ofCn2 on PCAFAB evolution characteristics; other parameters are the same as those in Figure 1(b), and the blue arrow head denotes each beamlet’s acceleration direction: (a) Cn2 = 10−13m−2/3, (z) = 50m; (b) Cn2 = 10−13m−2/3, (z) = 100m; (c) Cn2 = 10−13m−2/3, (z) = 500m; (d) Cn2 = 5 × 10−12m−2/3, (z) = 1000m; (e) zoom in on Figure 4(d) in the first quadrant; (f) Cn2 = 10−11m−2/3, (z) = 500m; (g) Cn2 = 10−9m−2/3, (z) = 100m; (h) Cn2 = 10−9m−2/3, (z) = 250m; (i) Cn2 = 10−8m−2/3, (z) = 50m; (j) Cn2 = 10−8m−2/3, (z) = 100m.
(a)(b)(c)(d)(e)(f)(g)(h)(i)(j)We end up with this section by discussing the PCAFABs through atmospheric turbulence which contains different numbers of the single Airy beam so that the conclusions obtained are universal. The results are presented in Figure5, where r0 and θ0 are the initial transverse displacement and angle displacement between two incoming Airy beams, respectively, n denotes the beam number, and other parameters are in accord with that in Figure 4(a). It reveals from these figures that the incident PCAFAB intensity profile is determined by the a0, ω0, r0, θ0, and n, whereas the beam transmission properties in atmospheric space are characterized by the Cn2, σ, and z; in other words, the derived analytical formulae and conclusions can be extended to the study on different numbers of array Airy beams’ transmission in atmospheric turbulence.Figure 5
Transverse intensity distribution of PCAFABs in atmospheric space for different value of (n); other parameters are the same as those in Figure 4(a), and two blue arrow heads denote the initial transverse displacement and angle displacement between two incoming Airy beams, respectively: (a) (n) = 2, (z) = 0 m; (b) (n) = 2, (z) = 500 m; (c) (n) = 3, (z) = 0 m; (d) (n) = 3, (z) = 500 m; (e) (n) = 6, (z) = 0 m; (f) (n) = 6, (z) = 500 m; (g) (n) = 8, (z) = 0 m, first kind; (h) (n) = 8, (z) = 500 m, first kind; (i) (n) = 8, (z) = 0 m, second kind; (j) (n) = 8, (z) = 500 m, second kind.
(a)(b)(c)(d)(e)(f)(g)(h)(i)(j)


## 4. Conclusion

In this paper, we have explored thoroughly the properties of the PCAFABs evoluting in the Kolmogorov atmospheric turbulence based on the extended Huygens–Fresnel integral and the Rytov phase structure function under the paraxial approximation. The influence of beam parameters and turbulence factors on the PCAFABs normalized intensity distribution at different cross sections of atmospheric environment was discussed in detail. We found that the side lobes of incident PCAFABs gradually diminish with the increasinga0, but the central main lobe location is rarely affected by the changes of a0. In stark contrast, σ plays an insignificant role on the incident PCAFABs’ intensity distribution, but the side loses of PCAFABs fade away during the transmission process with the σ decreases, and the beam self-bending property is increasingly obvious along with the changing of σ; the bigger the σ is, the more severe the bending is formed at the receiver plane, but these lines are nonmonotonic versus the increases of distance z. A mass of simulations manifest that the PCAFABs covert into the PCAGBs at the first inflexion; hence, equation (9) is no longer valid to describe the intensity of PCAFABs when the transmission distance z is larger than the inflexion A. Moreover, it is found that the side lobes of PCAFABs gradually diminish during transmission in atmospheric space, and the larger the Cn2 is, the shorter the transmission distance z is required to achieve complete disappearance of the side petals. We look forward to that the derived analytical formulae and the corresponding conclusions can be useful for learning the optical properties of the PCAFABs which contains different numbers of single Airy beam transmitting in the atmospheric turbulence, especially for free-space optical communication area.



---

*Source: 6649144-2021-02-08.xml*