# Scintillation index of a Gaussian Schell-model beam on slant atmospheric turbulence.

1. INTRODUCTIONOn propagation through the turbulent atmosphere, an optical wave will experience random irradiance fluctuations, called scintillation, due to perturbations from refractive index fluctuations. Scintillation can lead to power losses at the receiver, resulting in lower signal-to-noise ratios and eventually to fading of the received signal below a prescribed threshold [1]. So, optical scintillation is considered one of the most crucial atmospheric effects and has attracted much attention in practical applications, particularly in laser tracking and ranging systems [2-4] and imaging systems [5-7].

The topic of optical scintillation has been extensively studied for many years. In the early 1960s, Tatarskii and Cherenkov adopted Rytov approximation method [8-10] and obtained the scintillation index of unbounded plane wave and spherical wave, but their scintillation results were limited to weak fluctuations. The limitation of the weak fluctuation theory was clearly demonstrated, and saturation phenomenon of the scintillation was discovered by the experimental data of Gracheva and Gurvich [11] in 1965. A number of theoretical and experimental studies devoted to irradiance fluctuation under strong turbulence regimes then followed. With a theory of saturation developed by Gochelashvili and Shishov [12], several qualitative models of scintillation were subsequently developed and modified by others [13,14]. In 1999, Hopen and Andrews [15] intensity studied fluctuation variance of Gaussian beam propagation in moderate-to-strong fluctuation. Andrews et al. [1,16] made a comprehensive summary of these studies about plane wave, spherical wave, and Gaussian beam. Wu and Wei [17] discussed the scintillation index of a Gaussian beam considering inner-scale and outer-scale on slant path. A few papers reported the scintillation of a partially coherent beam. Many researchers concerned with the spread and average intensity of partially coherent beam [18, 19].

Studies above focused on the horizontal path; however structure constant [C.sup.2.sub.n] in many cases may vary along the propagation path as it does in vertical and slant paths. Beran and Whitman [20] calculated the scintillation index of a plane wave using the two-scale theory and assuming a Kolmogorov structure function with an altitude-dependent structure constant. Baker [21] developed a low-order turbulence (LOT) solution to describe Gaussian beam weak scintillation on ground-to-space propagation. Charnotskii [22] presented a more rigorous theory of beam scintillation on ground-to-space paths using path integrals and analytic techniques, where turbulence is concentrated in a relatively thin layer near the transmitter.

Especially after the introduction of free-space optical atmospheric communication links, the average BER can be significantly reduced by the use of a partially coherent source beam. Many researchers started to study the propagation characteristics of different types of sources in turbulence [23, 24]. Baykal and Eyyubogl [25] developed the scintillation index formula for a flat-topped Gaussian beam source on axis in turbulence. They found that using a flat-topped Gaussian beam in a FSO link will bring advantages in scintillations compared with Gaussian-beam incidence. Then Eyyuboglu and Baykal [26, 27] presented scintillation calculations in weak atmospheric turbulence for partially coherent general beams based on the extended Huygens-Fresnel principle and the scintillation index of cosh-Gaussian beams in a turbulent atmosphere using Rytov method.

In this paper, we use the techniques in [28], the altitude-dependent model of the ITU-R turbulence structure constant model, to develop a scintillation index model for a Gaussian Schell-model beam on the slant path that is applicable under weak-to-strong fluctuations. The result considers the focusing regime and also can be reduced to the result of the horizontal path when atmospheric structure constant is a fixed value. Finally, the numerical results are compared between the fully coherent beam and the partially coherent beam.

2. SCINTILLATION INDEX ON SLANT PATH

2.1. Weak Scintillation Index

The beam is off the optical axis, and the scintillation index is expressed as a sum of radial and longitudinal components in the form at a point.

[[sigma].sup.2.sub.I] (r,L) = [[sigma].sup.2.sub.Ir] (r,L)+ [[sigma].sup.2.sub.I](0,L), (1)

here the first component is radial, and the second is the longitudinal component.

The radial component of the scintillation index in weak fluctuation is given by [16]

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (2)

And the longitudinal scintillation index is defined under the Rytov approximation by [6]

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (3)

here, [I.sub.0](x) is a modified Bessel function, [[PHI].sub.n][kappa] the spatial power spectrum of refractive index fluctuations, and the Gaussian Schell-model source illumination beam in receiver plane without atmospheric effects is characterized by [LAMBDA], [THETA] beam parameters defined by

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (4)

where W represents the beam radius at the receiver. In arriving at Equation (4), here we have replaced the receiver beam parameters for a fully coherent beam in free space with the equivalent free-space receiver parameters for a partially coherent beam [28], where [[??].sub.s] is the source coherence parameter, which is the degree of the partial coherence of the source laser beam at the transmitter. For [[zeta].sub.s] = 1, the effective beam parameters in Equation (4) reduce to parameters of the fully coherent Gaussian beam, i.e., [[LAMBDA].sub.1] = [[LAMBDA].sub.0]/[[[THETA].sup.2.sub.0]+[[ LAMBDA].sup.2.sub.0]], [[THETA].sub.1] = [[THETA].sub.0]/[[[THETA].sup.2.sub.0]+[[ LAMBDA].sup.2.sub.0]]. The parameter [THETA] describes amplitude change in the wave caused by refraction, and [LAMBDA] describes amplitude change caused by diffraction.

2.2. Radial Component

In the presence of optical turbulence, additional diffraction and refraction cause further broadening of the beam spot size and further focusing and defocusing [29,31]. We will follow the work of Ricklin et al. [32] and replace the beam parameter A with the effective beam parameter Ae in the radial component, based on Equation (2). It has been shown that the use of effective beam parameters permits us to formally extend weak fluctuation expressions for the spatial coherence radius into the strong fluctuation regime [29, 30]. In the case where the inner scale of turbulence is smaller than the spatial coherence length of the optical wave, we base our calculations on the conventional Kolmogorov spectrum [[PHI].sub.n]([kappa]) = 0.033[C.sup.2.sub.n] ([xi]H)[[kappa].sup.-11/3], Equation (2) yields

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (5)

Let a = 2[[LAMBDA].sub.e]r[xi], b = [[LAMBDA].sub.e]L[[xi].sup.2]/[kappa], the evaluation of Equation (5) leads

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (6)

By use of the following formula

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (7)

where [sub.1][F.sub.1] (a, b ; c; x) is a hyper geometric function.

[[integral].sup.[infinity].sub.0] exp (-b[[kappa].sup.2])[[kappa].sup.-8/3] d[kappa] = 1/2 [GAMMA](-5/6)[b.sup.5/6]. (8)

In terms of Equation (7), Equation (8), the radial component Equation (6) leads to

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (9)

where [[sigma].sup.2.sub.0] = 1.23[C.sup.2.sub.n0][[kappa].sub.7/6][L.sup.11/6] is the Rytov variance.

The effective wave radius [W.sub.e](L) for the GSM beam through atmospheric turbulence and the global coherence parameter [zeta] are expressed by

[W.sub.e](L) = [W.sub.0][([[THETA].sup.2.sub.0] + [zeta][[ LAMBDA].sup.2.sub.0]).sup.1/2], [zeta] = [[zeta].sub.S] + 2[W.sup.2.sub.0]/[[rho].sup.2.sub.T], (10)

[[rho].sub.T] = [[1.46[[kappa].sup.2]L [[integral].sup.1.sub.0] [C.sup.2.sub.n] ([xi]L)[(1 - [xi]).sup.5/3] d[xi]].sup.-3/5]. (11)

[[LAMBDA].sub.e] = 2L/[kappa][W.sup.2.sub.e]. (12)

Physically speaking, the global coherence parameter is a measure of the global degree of coherence of light across each transverse plane along the propagation. It is used to define a related dimensionless quantity, the source coherence parameter [[zeta].sub.s] as the degree of partial coherence of the source beam at the transmitter. For [[zeta].sub.s] = 1, Equation (10) represents fully coherent beam; [[zeta].sub.s] increasing, source beam becomes less coherent. Here [[rho].sub.T] is the coherence length of a spherical wave propagating on the slant turbulent atmosphere, and [C.sup.2.sub.n]([xi]L) is the model for the refractive-index structure constant in the atmosphere. This effective wave radius for the partially coherent GSM beam through atmospheric turbulences agrees with that by Korotkova [33], where [W.sub.e] = [W.sub.1] [(1 + 4[q.sub.c][[LAMBDA].sub.1] + 1.63[[sigma].sup.2.sub.0][[ LAMBDA].sup.5/6.sub.1]).sup.1/2]. It is noted that Ref. [33] introduces a thin random phase screen to model the diffuser in front of the laser transmitter and characterize the phase screen by a power spectrum function rather than by a correlation function as in the Gaussian Schell-model. It should be noted that [W.sub.1] is Gaussian beam spot radius in receiver plane not partially coherent beam because the effects of partially coherent beam are accounted by [q.sub.c]. We note that the scintillation index for the partially coherent Gaussian beam presented in Equation (9) correctly reduces to the well-known scintillation index for the Gaussian beam on horizontal path when [[zeta].sub.s] is taken as unity, and atmospheric structure constant is a fixed value [C.sup.2.sub.n]([xi]H) = [C.sup.2.sub.n0] in Equation (9). Under weak fluctuations, the result is simplified to

[[sigma].sup.2.sub.I,r](r,L) [congruent to] 4.42[[sigma].sup.2.sub.0][[LAMBDA].sup.5/6.sub.1] [r.sup.2]/[W.sup.2.sub.1], r < [W.sub.1]. [[sigma].sup.2.sub.0] [much less than] 1. (13)

where [W.sub.1] is the free-space spot size radius for Gaussian beam at receiver plane. Here Equation (13) agrees with that in [1].

2.3. Longitudinal Component of Scintillation Index

In weak fluctuation regimes where inner-scale effects are negligible, we base our calculations on the Kolmogorov spectrum. By use of this spectrum model and let a = [LAMBDA]L[[xi].sup.2]/[kappa],b = L[xi](1 - [bar.[[THETA].sub.[xi]]])/[kappa], Equation (6) along the optical axis (r = 0) becomes

[[sigma].sup.2.sub.l] (0,L) = 0.264[[pi].sup.2][[kappa].sup.2]L [[integral].sup.1.sub.0][C.sup.2.sub.n]([xi]H) [[integral].sup.[infinity].sub.0] [[kappa].sup.-8/3] exp(-a[[kappa].sup.2]) [1-cos(b[[kappa].sub.2])] d[kappa]d[xi], (14)

There is a small error a = [LAMBDA]L[xi]/[kappa] in Equation (3.55) in [33]. Using Equations (7), (8), the evaluation of Equation (14) leads to

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (15)

Substituting a, b into Equation (15), we have

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (16)

Equation (16) is restricted to weak fluctuations. Under strong fluctuation conditions, we use the fluctuation theory developed by Andrews et al. [16], which yields

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (17)

where [[sigma].sup.2.sub.l] is defined by Equation (16). We combine the radial component Equation (9) with the longitudinal component Equation (17) to obtain

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (18)

Equation (18) represents our approximation to the scintillation index at any transverse position in the beam at the receiver plane of a general lowest-order Gaussian Schell-model beam wave under all irradiance fluctuation conditions.

3. ATMOSPHERIC TURBULENCE STRUCTURE CONSTANT MODEL

From above derivation, we can see that scintillation index depends largely on height distribution of turbulence atmospheric structure constant. There, we employ ITU-R atmospheric turbulence structure constant model [34] expressed by

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (19)

where [v.sub.RMS] = 21 (m/s) is the rms wind speed in meters per second (m/s). h is in meters (m), [C.sup.2.sub.n0] the structure parameter at the ground in [m.sup.-2/3], (its typical value is 1.7 x [10.sup.-14][m.sup.-2/3]).

We plot the logarithm of the ITU-R model Equation (19) in Figure 1 as a function of altitude for three values of the parameter [C.sup.2.sub.n0] and three values of high-altitude wind speed. [C.sup.2.sub.n0] affects the profile model only up to roughly 1 km, and high-altitude wind is the element most around altitudes of 10 km. It shows that the value of atmospheric structure constant [C.sup.2.sub.n] is the biggest near ground. From 1 km to 4 km, the atmospheric structure constant does not depend on sub aerial atmospheric structure constant and wind velocity due to atmospheric absorption. The ITU-R atmospheric structure constant mode is closer to the actual turbulence model than negative exponent index model.

4. NUMERICAL ANALYSIS CONCLUSION

In this section, we have analyzed the characteristics of the scintillation index with the variations of propagation length and root of the Rytov variance. Our results are shown in Figures 2-6. We have chosen the wavelength of operation to be [lambda] = 1.55 [micro]m, since this wavelength is the most widely used in the current FSO links. In computation, since the atmospheric structure constant is variable with height on slant propagation path, the computing formulas have to reserve [C.sup.2.sub.n](h) path integral calculus form; H is the height between transmitter and receiver. Then we deal with a Gaussian Schell-model collimated beam scintillation index.

Figure 2 shows the radial scintillation index as a function of the square root of the Rytov variance [[sigma].sub.0] for different values of the source coherence parameter [[zeta].sub.s] representing beams from the coherent to the partially coherent. The illustrated GSM beam is assumed collimated beam ([[THETA].sub.0] = 1), and the propagation distance is allowed to vary. At first, in weak fluctuation the radial scintillation index increases quickly with propagation distance, then it eliminates quickly once the beam has passed the focusing regime. In all irradiance fluctuation conditions, as the source beam becomes less coherent, the radial scintillation index decreases as expected. When [[sigma].sub.0] is very large, radial component of the scintillation index will vanish, independent of source coherence parameter.

In Figure 3, we compare behavior of the coherent beam with that of the partially coherent beam and plot the radial component of the scintillation index Equation (9) versus r for different values of the source coherence parameter [[zeta].sub.s]. With increasing off-axis parameter r, the radial component of the scintillation index will increase, but the radial scintillation index for the incoherent source ([[zeta].sub.s] = 1000) is zero independent of the off-axis parameter r. The Gaussian beam scintillation index is bigger than the partially coherent beam. This agrees with the result that the partially coherent beam is less sensitive to the effects of turbulence than fully coherence one.

In Figure 4, we show Equation (17) as a function of the propagation distance. As the source beam becomes less coherent, the scintillation index decreases at a given distance. With the increase of distance, longitudinal scintillations index increases and tends to saturation. For short ranges (up to 6-8 km), the partially coherent and incoherent beams assume lower levels of scintillation index than a fully coherent beam; therefore the advantage of the partially coherent beam can be expected in weak or focusing atmospheric regimes.

The effect of receiver altitude H on the longitudinal scintillation index is illustrated in Figure 5. Before focusing atmospheric regimes (location of peak scintillation), the higher the receiver is, the smaller the scintillation index of the partially coherent GSM beam is. In Figure 6, we show the total scintillation Equation (15) as a function of [[sigma].sub.0] with several values of height. We set [C.sup.2.sub.n0] = [10.sup.-13] [m.sup.-2/3] and allow the propagation distance to vary. The result, as in Figure 5, shows that the total scintillation index is largely dependent on the longitudinal scintillation index. When the atmospheric structure constant is a fixed value, our result is consistent with the result of horizontal path [30].

5. SUMMARY

With the motivation of understanding whether employing a GSM beam may bring some advantages in reducing the scintillation in long atmospheric optics telecommunications links, we evaluate the scintillation index of a GSM incidence on the slant path for all atmospheric turbulence conditions. With the change of altitude and source coherence parameters, the longitudinal and radial components of scintillation index are discussed, respectively. In general, the larger the source coherence parameter ss and the higher receiver H are, the smaller the scintillation index is. The scintillation index of the GSM is smaller on the slant path than in the horizontal path in optical turbulence. However, at relatively large [[sigma].sub.0], the scintillation advantage of the slant path tends to disappear. For incoherent beam, off axis radial scintillation index is zero. Combining the results in current work with previous result, we can assert that the GSM beam on slant path reduces scintillation index more than the fully coherent Gaussian beam on horizontal path. This scintillation index model of the GSM beam on slant path for laser beam propagation in atmosphere turbulence is significant for the visible and infrared imaging, laser tracking, controlling and guiding, and LIDAR systems.

ACKNOWLEDGMENT

This work is supported by the National Natural Science Foundation of China (60971079, 6127710), and Xianyang Normal School (07XSYK281).

REFERENCES

[1.] Andrews, L. C. and R. L. Phillips, Laser Beam Propagation through Random Media, SPIE Press, Bellingham, 1998.

[2.] Banakh, V. A. and V. L. Mironov, Lidar in a Turbulent Atmosphere, Artech House, Boston, 1987.

[3.] Bufton, J., "Scintillation statistics measured in an earth-space-earth retroreflector link," Appl. Opt., Vol. 16, 2654-2660, 1977.

[4.] Bufton, J., R. Iyer, and L. Taylor, "Scintillation statistics caused by atmospheric turbulence and speckle in satellite laser ranging," Appl. Opt., Vol. 16, 2408-2412, 1977.

[5.] Goodman, J. W., Statistical Optics, Wiley, New York, 1985.

[6.] Lawrence, T. W., D. M. Goodman, E. M. Johansson, and J. P. Fitch, "Speckle imaging of satellites at the U.S. Air Force Maui Optical Station," Appl. Opt., Vol. 31, 6307-6321, 1992.

[7.] Labreyrie, A., "Attainment of diffraction-limited resolution in large telescopes by Fourier analysing speckle patterns in star images," Astron. Astrophys., Vol. 6, 85-87, 1970.

[8.] Tatarskii, V. I., Wave Propagation in a Turbulent Medium, McGraw-Hill, New York, 1961.

[9.] Chernov, L. A., Wave Propagation in a Random Medium, McGraw-Hill, New York, 1960.

[10.] Ishimaru, A., Wave Propagation and Scattering in Random Media, Academic Press, New York, 1978.

[11.] Gracheva, M. E. and A. S. Gurvich, "Strong fluctuations in the intensity of light propagated through the atmosphere close to the earth," Izvestiya VUZ Radiofizika, Vol. 8, 717-724, 1965.

[12.] Gochelashvili, K. S. and V. I. Shishov, "Saturated fluctuations in the laser radiation intensity in a turbulent medium," Sov. Phys. JETP, Vol. 39, 605-609, 1974.

[13.] Fante, R. L., "Inner-scale size effect on the scintillations of light in the turbulent atmosphere," J. Opt. Soc. Am., Vol. 73, 277-281, 1983

[14.] Frehlich, R. G., "Intensity covariance of a point source in a random medium with a Kolmogorov spectrum and an inner scale of turbulence," J. Opt. Soc. Am. A, Vol. 4, 360-365, 1987.

[15.] Hopen, C. Y. and L. C. Andrews, "Optical scintillation of a Gaussian beam in moderate-to-strong irradiance fluctuations," Proc. SPIE, Vol. 4, 142-150, 1999.

[16.] Andrews, L. C., R. L. Phillips, and C. Y. Hopen, Laser Beam Scintillation with Applications, SPIE, 2001.

[17.] Wu, Z. S. and H. Y. Wei, "Study on scintillation considering inner- and out-scales for laser beam propagation on the slant path through the atmospheric turbulence," Progress In Electromagnetics Research, Vol. 80, 277-293, 2008.

[18.] Li, J., Y. Chen, S. Xu, Y. Wang, M. Zhou, Q. Zhao, Y. Xin, and F. Chen, "Average intensity and spreading of partially coherent four-petal Gaussian beams in turbulent atmosphere," Progress In Electromagnetics B, Vol. 24, 241-261, 2010.

[19.] Wang, F., Y. Cai, H. T. Eyyuboglu, and Y. K. Baykal, "Average intensity and spreading of partially coherent standard and elegant Laguerre-Gaussian beams in turbulent atmosphere," Progress In Electromagnetics Research, Vol. 103, 33-56, 2010.

[20.] Beran, M. J. and A. M. Whitman, "Scintillation index calculations using an altitude-dependent structure constant," Appl. Opt., Vol. 27, 2178-2182, 1988.

[21.] Baker, G. J., "Gaussian beam weak scintillation: Low-order turbulence effects and applicability of the Rytov method," J. Opt. Soc. Am. A, Vol. 23, 395-417, 2006.

[22.] Charnotskii, M., "Beam scintillations for ground-to-space propagation. Part I: Path integrals and analytic techniques," J. Opt. Soc. Am. A, Vol. 27, 2169-2179, 2010.

[23.] Chu, X., "Evolution of beam quality and shape of Hermit Gaussian beam in non-Kolmogorov turbulence," Progress In Electromagnetics Research, Vol. 120, 339-353, 2011.

[24.] Li, Y.-Q., Z.-S. Wu, and L.-G. Wang, "Polarization characteristics of a partially coherent Gaussian Schell-model beam in slant atmospheric turbulence," Progress In Electromagnetics Research, Vol. 121, 453-468, 2011.

[25.] Baykal, Y. and H. T. Eyyuboglu, "Scintillation index of flat-topped Gaussian beams,"Appl. Opt., Vol. 45, 7066-7073, 2006.

[26.] Eyyuboglu, H. T. and Y. Baykal, "Scintillation characteristics of cosh-Gaussian beams," Appl. Opt., Vol. 46, 1099-1106, 2007.

[27.] Eyyuboglu, H. T., Y. Baykal, and Y. Cai, "Scintillation calculations for partially coherent general beams via extended Huygens-Fresnel integral and self-designed Matlab function," Appl. Phys. B, Vol. 100, 597-609, 2010.

[28.] Miller, W. B., J. C. Ricklin, and L. C. Andrews, "Scintillation of initially convergent Gaussian beams in the vicinity of the geometric focus," Appl. Opt., Vol. 34, 7066-7073, 1995.

[29.] Rickin, J. C. and F. M. Davidson, "Atmospheric optical communication with a Gaussian Schell beam," J. Opt. Soc. Am., Vol. 20, 856-863, 2003.

[30.] Andrews, L. C., W. B. Miller, and J. C. Ricklin, "Spatial coherence of a Gaussian beam in weak and strong optical turbulence," J. Opt. Soc. Am. A., Vol. 11, 1653-1660, 1994.

[31.] Young, C. Y. and L. C. Andrews, "Effects of a modified spectral model on the spatial coherence of a laser beam," Waves Random Media, Vol. 4, 385-397, 1994.

[32.] Ricklin, J. C., W. B. Miller, and L. C. Andrews, "Effective beam parameters and the turbulent beam waist for initially convergent Gaussian beams," Appl. Opt., Vol. 34, 7059-7065, 1995.

[33.] Korotkova, O., "A model for a partially coherent Gaussian beam in atmospheric turbulence with applications for lasercom and lidar systems," M.S. University of Central Florida, 2003.

[34.] ITU-R Document 3J/31-E, "On propagation data and prediction methods required for the design of space-to-earth and earthto-space optical communication systems," Radiocommunication Study Group Meeting, International Telecommunication Union, 2001

Ning-Jing Xiang * and Zhen-Sen Wu

School of Science, Xidian University, P. O. Box 48, Xi'an 710071, China

Received 30 December 2012, Accepted 8 April 2013, Scheduled 18 April 2013

* Corresponding author: Ning-Jing Xiang (xiangningjing@sohu.com).

Printer friendly Cite/link Email Feedback | |

Author: | Xiang, Ning-Jing; Wu, Zhen-Sen |
---|---|

Publication: | Progress In Electromagnetics Research M |

Article Type: | Report |

Geographic Code: | 9CHIN |

Date: | Mar 1, 2013 |

Words: | 3791 |

Previous Article: | Design of metaparticles as sharp frequency-selective obscurant aerosols. |

Next Article: | Analytical investigation of periodic coplanar waveguides. |

Topics: |