# Waveform Distortion of Gaussian Beam in Atmospheric Turbulence Simulated by Phase Screen Method

Authors: Zhi-Qiang Yang; Li-hong Yang; Lei Gong; Liguo Wang; Xinyi Wang
Journal: Journal of Mathematics
 (2022)

---

## Abstract

The atmospheric turbulence phase screen is generated based on the power spectrum inversion method, and the multiple transmission processes are statistically averaged. The waveform distortion of the Gaussian beam in the atmospheric turbulence is analyzed; the simulation results show that the property of Gaussian beam has been destroyed after its passing through the atmospheric turbulence. When the beam waist radius is close to a certain radius, the degree of change becomes larger. When it approaches the critical value, the wave surface no longer changes sharply, and as the turbulence intensity increases, the phase fluctuation becomes more and more severe, the coherence of the beam is destroyed, and the spot may be split into several pieces. Finally, the relationship of intensity fluctuation, amplitude fluctuation, and bit error rate with distance is analyzed.


---

## Body

## 1. Introduction

The particularity of light propagation in turbulent media lies in the random fluctuation of the refractive index of the medium; the key of numerical simulation for laser atmospheric propagation is to construct a phase screen that can correctly reflect the statistical characteristics of atmospheric turbulence [1–4]. The earliest phase screen simulation turbulence method is proposed by Fleck, who used the “multiphase screen method” in 1976 to simulate the influence of atmospheric turbulence in free space on the transmitted beam. This method is widely used to simulate the atmospheric turbulence effect of the beam [5]. The current methods of generating phase screens are mainly to filter Gaussian white noise by using the refractive index power spectrum of the turbulence, by which the disturbed wavefront phase can be obtained through Fourier transform; this method is called “power spectrum inversion method.” Since Mcglamery proposed the power spectrum inversion method in 1967, this method has been widely used to simulate turbulence [6]. In 1976, Noll proposed to use orthogonal and complete two-dimensional basis functions (such as Zenike polynomials and K-L polynomials) to directly obtain the distorted wavefront of turbulence [7]. In 1983, Walner proposed a method of using the covariance of the phase structure function to generate the Kolmogorov turbulence phase screen [8], but this method is computationally expensive. In 1990, Bowman proposed an algorithm for reconstructing the turbulent wavefront using wavelet transform [9]. Among many methods, the power spectrum inversion method is applicable to a wide range of atmospheric turbulence spectrum models [10–12].Based on the Gaussian beam horizontal propagation theory, this paper uses MATLAB to carry out detailed simulation on the power spectrum inversion method and analyzes the simulation results. The propagation characteristics of the Gaussian beam under the atmospheric turbulence spectrum are discussed.


## 2. Theoretical Model of Laser Transmission in Turbulent Atmosphere

Assuming that the monochromatic electromagnetic wave is transmitted along theZ direction, the wave field equation is(1)U⇀x,y,z,t=Ux,y,zexp−iωt−kz,where ω is the frequency of monochromatic electromagnetic waves. When it propagates in the atmosphere (nonmagnetic medium μ = 1) without free charge (ρ0 = 0) and free current (J0 = 0), its wave equation can be expressed as(2)∇2Ux,y,z,t+k2nr⟶2Ux,y,z,t=0,where nr⇀ is a function of r⇀ which represents the local refractive index at this point. It is very difficult to accurately solve this basic equation of light transmission in atmosphere.By removing the vector operator, the equation can be turned into a scalar equation:(3)2ik∂Ux,y,z∂z+∇2Ux,y,z+2k2n1rUx,y,z=0,where Ux,y,z is a slowly varying function of z; it can only change at the distance of random medium scale l, so as long as l >> λ; then, ∂2U/∂z2≪2κ∂U/∂z, and 2n1r⇀ approximately replaces 2n1r⇀+n12r⇀ (n1r⇀ represents the fluctuation of atmospheric refractive index at position r⇀). By simplifying formula (3), the following formula can be obtained:(4)2ik∂Ux,y,z∂z+∇⊥2Ux,y,z+2k2n1rUx,y,z=0.Equation (4) is called the parabolic equation or quasi-optical approximation relationship, which is suitable for the propagation of narrow-angle expanded beams. k is equal to 2π/λ which represents wave number in free space, and ∇⊥2=∂2/∂x2+∂2/∂y2. The infinitely extending random medium can be divided into many sections with thickness of Δz; the influence of turbulence in each Δz section on the phase of light waves can be concentrated on a thicknessless phase screen, but this phase screen has no effect on the amplitude of the light wave. The change of light wave amplitude is the cumulative result of free-space diffraction at the Δz distance between many two-phase screens.In order to specify each random phase screen, the correlation function and power spectrum of the phase fluctuation on each phase screen can be calculated. Assuming that there is no diffraction on each phase screen, from the perspective of geometric optics, the correlation function of the phase fluctuation can be written:(5)Γθρ⇀1,ρ⇀2=k2∫0Δz∫0Δznrz1,ρ⇀1n1z2,ρ⇀2dz1dz2=k2ΔzΓnρ⇀1−ρ⇀2,where ρ⟶ is the vector in the coordinate plane perpendicular to the z direction (the xy plane), Γθρ⇀ is the correlation function of the phase fluctuation, and Γnρ⇀ is the correlation function of refractive index fluctuation. Assuming Az>ρ0 (ρ0 is the coherence length of the turbulence), formula (5) can be simplified:(6)Γnρ⇀=2π∫−∞∞∫−∞∞ϕnk⇀⊥,kz=0expik⇀⊥⋅ρ⇀dk⇀⊥.So, the relationship between the power spectrum function and the power spectrum of the turbulent refractive index of each phase screen can be obtained:(7)ϕθk⇀⊥=2πk2Δzϕnk⇀⊥,kz=0.Because the transmission of the light wave field between two adjacent phase screens can be approximated as the transmission in free space, which represents the fluctuation of the refractive indexn1r=0, it can be obtained from equation (4):(8)∂Ux,y,z∂z=i2k∇⊥2Ux,y,z.It can be obtained by Fourier transform of two dimensional on both sides at the same time:(9)∂U˜k⇀⊥,z∂z=−ik⇀⊥2kUk⇀⊥,z,where k⟶⊥ is the wave vector perpendicular to the transmission direction; by integral in the Z direction on both sides of equation (9), solution can be obtained:(10)∂U˜k⇀⊥,z+δz=U˜k⇀⊥,zexp−ik⊥2Δz2k.Equation (10) shows the free-space diffraction between the two-phase screens.As shown in Figure1, after replacing the infinitely extending turbulent medium with a large number of phase screens, the transmission of light waves through the turbulent atmosphere is equivalent to the transmission of light waves between phase screens after passing through the phase screens.(a)
The contribution of the phase screen to the phase of the light wave field is as follows:(11)Uρ⟶,zi+=Uρ⇀,zi−expiθρ⇀,wherezi+ represents the right side of the ith phase screen and zi− represents the left side of the ith phase screen.(b)
Perform Fourier transform onUρ⇀,zi+ to get U˜k⇀⊥,zi+.(c)
Through the free-space transmission with a distance ofΔz,(12)U˜k⇀⊥,zi+1−=U˜k⇀⊥,zi+exp−ik⊥2Δz2k.(d)
By performing Fourier transform onU˜k⇀⊥,zi+1+, Uρ⇀,zi+1− can be obtained:(13)Uρ⇀,n+1Δz−=Ft−1FtUρ⇀,nΔz−×expiθnρ⇀exp−ik⊥22kΔz.Figure 1
Schematic diagram of simulating the propagation of light waves in turbulent atmosphere with the phase screen approximation method.Then, the transmission of light waves in the entire infinitely extending turbulent medium can be realized by repeatedly using equation (13) until the transmission of the last phase screen is completed.


## 3. Numerical Simulation Method

The parabolic equation is used to express the light field when the beam propagates in the turbulent atmosphere. The propagation path of the beam can be divided into two parts, the free space and the thin phase screen dispersed in it. When the light beam is transmitted to the phase screen, the wave front of the light is added with the corresponding disturbance part, and it passes through a section of free space to reach the next phase screen. When the beam reaches the receiving end, the entire simulation process ends. From equation (13), the light field at any adjacent phase screen can be expressed by the following formula [13, 14]:(14)Uρ⟶,zj+1=Ft−1FtUρ⟶,zjexpiθnρ⟶exp−ikx22k+ky22kΔzj+1.Equation (14) represents the turbulent phase screen, and kx and ky represent the wave number. The light wave transmission distance is divided into NZ sections, and the jth phase screen is in zj, Δzj is the phase screen spacing, and each phase screen is divided into N × N square grids whose grid spacing is Δx.Power spectrum of the atmospheric turbulence is used to filter the complex Gauss random number matrix; by inverse Fourier transform, the atmospheric disturbance phase will be obtained, which is the core of the power spectrum inversion method.MATLAB is used to generate a complex Gaussian random number matrix which has a mean value of 0 and a variance of 1, and then, the atmospheric refractive index power spectrum is introduced; the modified Von Karman spectrum is as follows:(15)Φnk=0.033Cn2exp−k2km2k2+k02−11/6,0<k<∞.The Kolmogorov spectrum is as follows:(16)Φnκ,z=0.033Cn2zκ−11/3,which is used to produce a continuous phase screen by filtering matrix a with the atmospheric refraction power spectrum [14]:(17)ϕx,y=C∑kx∑kyRκx,κyFϕκx,κyejκxx+κyy,where R represents the spatial spectrum of random function with normal distribution with mean value of 0 and variance of 1 and F represents atmospheric turbulence power spectrum. Δκx and Δκy are the sampling interval, and the constant C comes from the scale factor, which is used to adjust the phase screen variance. Fresnel diffraction is adopted in the vacuum transmission part, and the light field after the beam transmission distance z can be expressed as [15](18)Ux,y=−iλzexpikzexpikx2+y22z∬U0x0+y0×expikx2+y22z×exp−i2πfx0+fy0dx0dy0,where f is the spatial frequency; by Fourier transform of U0x0+y0expikx02+y02/2, the transmission in the vacuum part can be realized.


## 4. Numerical Simulation Results

The simulation parameters are set as follows: wavelength is 1.55μm, number of sampling points is 256 × 256, phase screen width is 0.5 m, phase screen spacing is 300 m, and 1000 phase screens are constructed.For the modified Von Karman spectrum, atmospheric turbulence phase screen for the different atmospheric refractive index structure constant and distortion of Gaussian beam after passing through the phase screen are simulated, and the results are as follows.Figure2 is the atmospheric turbulence phase screen generated by the power spectrum inversion method. Figures 2(a) to 2(f) are the atmospheric turbulence phase screens when the atmospheric refractive index structure constant Cn2 takes different values. From figures, it is can be seen that the turbulence intensity increases gradually with the increase of atmospheric refractive index structure constant Cn2.Figure 2
Atmospheric turbulence phase screen for the different atmospheric refractive index structure constantsCn2. (a) Cn2=2×10−17m−2/3, (b) Cn2=2×10−17m−2/3, (c) Cn2=2×10−15m−2/3, (d) Cn2=2×10−15m−2/3, (e) Cn2=2×10−13m−2/3, and (f) Cn2=2×10−13m−2/3.
(a)(b)(c)(d)(e)(f)Figure3 shows the light intensity distribution diagram of the initial beam with the beam waist radius w0=4cm. Figure 3 shows Gaussian beam after passing through turbulence simulated by the phase screen theory, where Figures 4(a)–4(f) represent the phase fluctuation and light intensity distribution of the distorted Gaussian beam after passing through phase screens of different intensities.Figure 3
Initial Gaussian beam. (a) Spot pattern of initial Gaussian beam. (b) Stereograph of initial Gaussian beam.
(a)(b)Comparing Figures4(a), 4(c), and 4(e) with Figure 3(a), after the Gaussian beam passes through atmospheric turbulence, the Gaussian property of the beam is destroyed. Compared with the initial Gaussian beam, the wave surface becomes rugged, the faster the change within the beam waist radius, but when approaching a certain distance from the center, the wave surface no longer changes sharply but up and down changes. As the atmospheric refractive index structure constant increases, that is, the turbulence intensity increases, in fact, the three images correspond to weak turbulence, medium turbulence, and strong turbulence. When the beam waist radius is close to a certain radius, its degree of change keeps getting bigger. When it approaches the critical value, the wave surface no longer changes sharply, but as the distance is close to the center point, the greater the turbulence intensity, the greater the ups and downs of the wave surface. By comparing Figures 4(b), 4(d), and 4(f) with Figure 3(b), it can be obtained after the Gaussian beam passing through the atmospheric turbulence phase screen which is generated by the spectrum inversion method; the spot becomes blurred and scattered, and the focusing ability becomes worse, which indicates the phase of the Gaussian beam fluctuates. With the increase of Cn2, the intensity of atmospheric turbulence continues to increase, the phase fluctuation becomes more and more severe, the coherence of the beam is destroyed, and the light spot may be split into several pieces. Compared with the initial Gaussian beam, the intensity of the beam attenuates significantly, and the larger the Cn2, the greater the attenuation of the light intensity.Figure 4
Gaussian beam after passing through the phase screen. (a)Cn2=2×10−17m−2/3, (b) Cn2=2×10−17m−2/3, (c) Cn2=2×10−15m−2/3, (d) Cn2=2×10−15m−2/3, (e) Cn2=2×10−13m−2/3, and (f) Cn2=2×10−13m−2/3.
(a)(b)(c)(d)(e)(f)


## 5. The Effect of Turbulent Atmosphere on Free-Space Optical Communications

In Free-Space Optical Communications (FSO) [16, 17], the laser from the transmitting end can be treated as a plane wave after being collimated by an optical lens [18, 19]; for the plane wave,(19)x2=σx2=0.31Cn2k7/6L11/6.From this, the relationship between bit error rate and light intensity fluctuation can be obtained [20, 21]:(20)BER=12erfc42×1.24Cn2k7/6L11/6.Under different fluctuation conditions, the variation of logarithmic amplitude fluctuation and logarithmic intensity fluctuation with transmission distance are shown in Figures5 and 6, and the relationship between bit error rate and logarithmic light intensity fluctuation variance is shown in Figure 7. In addition, the influence of transmission distance on bit error rate is also given, as shown in Figure Figure 8.Figure 5
Variation of amplitude fluctuation with transmission distance.Figure 6
Variation of intensity fluctuation with transmission distance.Figure 7
Variation of bit error rate with light intensity fluctuation.Figure 8
Variation of bit error rate with transmission distance.Figures5 and 6 show the variation of logarithmic amplitude fluctuation and logarithmic intensity fluctuation with transmission distance under different intensity turbulence conditions. It can be clearly seen from the figures that when the transmission distance is constant, the stronger the turbulence intensity, the greater the fluctuation of logarithmic amplitude and logarithmic intensity.In Free-Space Optical communication (FSO) system with a rate below 2.5 Gbps, the system’s bit error rate is generally required to be below 10−9. It can be seen from Figure 7 that, under weak fluctuation conditions, for the atmospheric laser communication system’s bit error rate less than the requirement, the light intensity fluctuation should be less than 0.67. It can be seen from Figures 7 and 8 that when the transmission distance is constant, with the increase of turbulence intensity, the bit error rate rises quickly. For the bit error rate 10−9, when Cn2 is equal to 2×10−15m−2/3, the effective distance of communication is close to 4.6 km. With the increase of Cn2, the effective distance of transmission is getting smaller and smaller; when it reaches 2×10−13m−2/3, the effective distance of transmission reduced to about 0.37 km.


## 6. Conclusion

By comparing the phase screens under atmospheric turbulence of different intensities, it is found that atmospheric turbulence of different intensities affects the propagation of Gaussian beams in terms of light intensity and phase fluctuations. The greater the intensity of the atmospheric turbulence is, the greater the light intensity of the Gaussian beam is affected, the more obvious the fluctuation of the light intensity, and the more obvious the tendency of the light spot to split into several pieces. In terms of phase fluctuations, the greater the intensity of the turbulence, the more intense the phase fluctuations of the Gaussian beam, the larger the difference of the light spot in the figure, and the more obvious the chromatic aberration, which puts forward higher and more specific requirements in the accuracy and reliability of laser transmission under atmospheric turbulence of different intensities. With the increase of turbulence intensity, the bit error rate rises quickly. With the increase ofCn2, the effective distance of transmission is getting smaller and smaller; when it reaches 2×10−13m−2/3, the effective distance of transmission reduced to about 0.37 km.



---

*Source: 5939293-2022-02-21.xml*