## Abstract

The source and receiver plane characteristics of flat topped (FT) beam propagating in turbulent atmosphere are investigated. To this end, source size, beam power and *M*
^{2} factor of source plane FT beam are derived. For a turbulent propagation medium, via Huygens Fresnel diffraction integral, the receiver plane intensity is found. Power captured within an area on the receiver plane is calculated. Kurtosis parameter and beam size variation along the propagation axis are formulated. Graphical outputs are provided displaying the variations of the derived source and receiver plane parameters against the order of flatness and propagation length. Analogous to free space behavior, when propagating in turbulence, the FT beam first will form a circular ring in the center. As the propagation length increases, the circumference of this ring will become narrower, giving rise to a downward peak emerging from the center of the beam, eventually turning the intensity profile into a pure Gaussian shape.

©2006 Optical Society of America

## 1. Introduction

Flat topped (FT) beams are derived from the fundamental Gaussian beam by introducing an order of flatness parameter. They possess a number of attractive features. Amongst these are less spreading in propagation and medical applications [1]. It is the former that is of interest to us in this article. From propagation point of view, there exists a wide selection of source beam formulations. Broadly considered, four basic classifications are possible, namely, super-Gaussian, flattened-Gaussian, Fermi-Dirac and super-Lorentzian [2, 3]. Among these, our model is based on flattened-Gaussian beam of Li [3, 4] which consists of taking the original flattened-Gaussian case proposed by Gori [5] and including the Gaussian exponential inside the summation after applying the conditions of flatness. In all cases however, there is an order parameter whose role is to control the flatness of the beam.

Recently we have analyzed the scintillation properties of flattened-Gaussian type beams, concluding that it will have favorable characteristics at certain propagation ranges, source sizes and operating wavelengths [6].

Up to now, some studies have been conducted by Cai, Lü and other researchers, mainly exploiting the free space propagation aspects of FT beams [7–14]. Very few studies concerning propagation in atmospheric turbulence have appeared in the literature [15]. On the other hand a comprehensive analysis of the source and the receiver plane parameters of FT beams in the presence of atmospheric turbulence have not so far been undertaken. In this respect, our current study has the following novelties

1. The inclusion of turbulence in the propagation of FT beams

2. A comprehensive analysis of source and receiver plane parameters of FT beams.

## 2. Formulation

#### 2.1. Formulation of the source field and its related parameters

Our source field expression is adopted from Eq. (8) of [3], by dropping the order parameter in the denominator. Hence the field expression, *u*_{s}
(*s*) becomes

in Eq. (1), (*s*_{x}*,s*_{y}
) designates the decomposition of the vector **s** into *x* and *y* components of the transverse source plane, *N* is the order parameter for flatness, such that at *N*=1, *u*_{s}
(*s*) will yield the fundamental mode (Gaussian) source field, correspondingly *α*_{sx}
and *α*_{sy}
becoming the Gaussian source sizes along *s*_{x}
and *s*_{y}
directions. When *N*≠1, on the other hand, the beam sizes will depart from Gaussian source sizes, and according to the definition developed by Carter [16], will be given by

where, *I*_{s}
(*s*) denotes the intensity of the source plane to be obtained from Eq. (1) as follows

Hence upon substitution into Eq. (2) from Eqs. (1) and (4), and expanding the terms in the square bracket as binomial series, thus performing the integrations, beam size of the *s*_{x}
direction will read

where the forms of $\left(\begin{array}{c}{C}_{1}\\ {C}_{2}\end{array}\right)$ refer to binomial coefficients such that, $\left(\begin{array}{c}{C}_{1}\\ {C}_{2}\end{array}\right)={C}_{1}!\u2044[\left({C}_{1}-{C}_{2}\right)!{C}_{2}!]$. In Eq. (5), the expression for *α*_{syN}
is attained by simply interchanging *α*_{sx}
with *α*_{sy}
. From Eq. (5), it is clear that at *N*=1, αsxN will appropriately reduce to *α*_{sx}
.

The power contained in the FT beam of flatness order *N*, denoted as *P*_{sN}
, can be calculated by integrating the intensity over the entire transverse plane, which means that it is equivalent to the denominator in Eqs. (2) or (3). This way *P*_{sN}
will become

*M*
^{2} is reported to be a measure of beam quality [13]. Its formulation, again separated into *s*_{x}
and *s*_{y}
components, is given by

where the pairs ${\sigma}_{\mathit{\text{sx}}}^{2}$ and ${\sigma}_{\mathit{\text{sy}}}^{2}$, ${\sigma}_{\mathit{\text{fx}}}^{2}$ and ${\sigma}_{\mathit{\text{fy}}}^{2}$ are respectively, the second moments associated with spatial and frequency spatial domains that are described, for the *x* direction, by the following equations

Here the *f*_{x}
and *f*_{y}
are the variables of frequency spatial domain. *I*_{s}
(*s*_{x}
) and *I*_{f}
(*f*_{x}
) are attained by respectively decoupling the *x* and *y* components in the manner *I*_{s}
(*s*_{x}*,s*_{y}
)=*I*_{s}
(*s*_{x}
)*I*_{s}
(*s*_{y}
), *I*_{f}
(*f*_{x}*,f*_{y}
)=*I*_{f}
(*f*_{x}
)*I*_{s}
(*f*_{y}
) where the function, *I*_{f}
(*f*_{x}*,f*_{y}
), is retrieved from the two dimensional spatial frequency Fourier transform of *I*_{s}
(*s*_{x}*,s*_{y}
). Identical expressions will hold for the *y* component. By applying the same steps as before, that is resorting to the use of binomial expansion, the ${M}_{x}^{2}$ of Eq. (7) will develop into

Note that ${M}_{y}^{2}$ will be identical to ${M}_{x}^{2}$. Thus, both ${M}_{x}^{2}$ and ${M}_{y}^{2}$ carry no dependency on the source parameters involving the transverse directions. Furthermore, the definition of *M*
^{2} is confined to the source plane, therefore for an FT beam considered in this paper, Eq. (11) will be valid regardless of propagation medium.

Before ending this section, we would like to note that, although it is theoretically possible to have selected *N* to be non-integer as pointed out by Li [4], we limit the present study to *N* being (positive) integer, since, the summations would otherwise turn into infinite series, thus preventing us gaining a proper insight into the subject.

#### 2.2. Formulation of the average receiver intensity and its related parameters

After the FT beam has propagated in the turbulent atmosphere, the average intensity on a receiver plane located at axial distance of *L* away from the source plane, will be supplied via the Huygens Fresnel integral in the following way

$$\times <\mathrm{exp}\left[\psi ({\mathbf{s}}_{1},\mathbf{p})+{\psi}^{*}({\mathbf{s}}_{2},\mathbf{p})\right]>,$$

where (**P**)=(*p*_{x}*,p*_{y}
) are the x and y components of the transverse receiver plane, *k*=2*π/λ* with *λ* being the wavelength of propagation. The ensemble average term appearing on the second line of Eq. (12) is [17]

where *D*_{ψ}
(s_{1}-s_{2}) represents the wave structure function, and *ρ*
_{0}=(0.545${C}_{n}^{2}$
*k*
^{2}
*L*)^{-3/5} is the coherence length of a spherical wave propagating in the turbulent medium, ${C}_{n}^{2}$ indicates the structure constant. As expressed in Eq. (13), in our derivation we employ the quadratic approximation of Rytov’s phase structure function.

When the source field expression of Eq. (1) is substituted into Eq. (12), thereby benefiting from the term by term integration facility invoked by the use of binomial expansion, the average receiver intensity takes the form of

$$+{\displaystyle \sum _{n=0}^{N}}{\displaystyle \sum _{{n}_{1}=0}^{n}}{\displaystyle \sum _{{n}_{2}=0}^{{n}_{1}}}{\left(-1\right)}^{{n}_{1}}\left(\begin{array}{c}N\\ n\end{array}\right)\left(\begin{array}{c}n\\ {n}_{1}\end{array}\right)\left(\begin{array}{c}{n}_{1}\\ {n}_{2}\end{array}\right)\frac{\beta {\alpha}_{\mathit{sx}}^{2}{\alpha}_{\mathit{sy}}^{2}}{{\left({\gamma}_{n{n}_{1}{n}_{2}x}{\gamma}_{n{n}_{1}{n}_{2}y}\right)}^{1\u20442}}\mathrm{exp}\left(\frac{{\beta}_{n{n}_{1}}{\alpha}_{\mathit{sx}}^{2}{p}_{x}^{2}}{2{\gamma}_{n{n}_{1}{n}_{2}y}}+\frac{{\beta}_{n{n}_{1}}{\alpha}_{\mathit{sy}}^{2}{p}_{y}^{2}}{2{\gamma}_{n{n}_{1}{n}_{2}y}}\right),$$

where,

$\beta ={k}^{2}{\rho}_{0}^{2},\phantom{\rule{.5em}{0ex}}{\beta}_{n}=n\beta ,\phantom{\rule{.5em}{0ex}}{\beta}_{n{n}_{1}}=\left({n}_{1}-2n\right)\beta ,\phantom{\rule{.5em}{0ex}}{\zeta}_{\mathit{nx}}=\left(\beta {\alpha}_{\mathit{sx}}^{2}+2{L}^{2}n\right)+i\left(k{\rho}_{0}^{2}\mathit{Ln}\right),$

${\gamma}_{n{n}_{1}{n}_{2}x}=\beta {\alpha}_{\mathit{sx}}^{4}+2{\alpha}_{\mathit{sx}}^{2}{L}^{2}\left(2n-{n}_{1}\right)+{\rho}_{0}^{2}{L}^{2}\left(n-{n}_{2}\right)\left(n-{n}_{1}+{n}_{2}\right)+\mathit{ik}{\rho}_{0}^{2}{\alpha}_{\mathit{sx}}^{2}L\left(2{n}_{2}-{n}_{1}\right),$

and ${\zeta}_{\mathit{ny}},{\gamma}_{n{n}_{1}{n}_{2}y}$ are the *y* counterparts of ${\zeta}_{\mathit{nx}}{,\gamma}_{n{n}_{1}{n}_{2}x}$ in which *x* is replaced by *y*.

To examine the beam size variation along the propagation axis, Eq. (4) is employed, with the modification that in the numerator *I*_{s}
(*s*) is replaced with *I*_{r}
(**p**) from Eq. (14). Thus for the *p*_{x}
direction, the beam size during propagation will be established as

Owing to the nature of turbulent atmosphere propagation, there exists no attenuation or absorption in the medium. Thus, we are able to retain the integrand in the numerator of Eq. (15) to be that of the source plane. Looking at Eq. (14), it is clear that due to the complexity of the expression for *I*_{r}
(**p**), we cannot proceed analytically; therefore numerical evaluations have to be sought for finding the beam size during propagation.

The kurtosis parameter of the propagating beam is defined as the ratio of the fourth and second moments of the intensity in the following way [18]

*K*_{y}
is found by switching the receiver plane coordinate parameter to *p*_{y}
. Again, the complexity of the *I*_{r}
(**p**) expression does not allow us to acquire an analytic result for the kurtosis parameter.

Finally in this section, we state the power that will fall into a receiver aperture of radius *α*_{r}
, which is also known as power in bucket, that is

In dividing by the source power, *P*_{sN}
, i.e. the term in the denominator of Eq. (17), it is envisaged that each individual FT beam of order *N*, is normalized with respect to its own order of flatness.

## 3. Results and discussions

#### 3.1. Source plane

Taking Eq. (4), in Fig. 1, we initially plot the 3D picture of source intensity at fixed Gaussian source sizes of *α*_{sx}
=*α*_{sy}
=3 cm, against the flatness order N. As seen in this figure, the source intensity of the FT beam will become flatter with escalating N values, eventually assuming the shape of a cylinder if equal Gaussian source sizes are used. Figure 1 also makes it evident that higher *N* values will give rise to larger beam sizes and source powers.

Figure 2 contains the overlaid contour plots of the beams belonging to Fig. 1, cut along the slanted axis of the transverse plane. From Fig. 2, the increases in beam size and source power of the FT beam at higher values of N are detected once more.

The observations in Figs. 1 and 2 about the variations of beam size and source power against flatness order N are in conformity with the associated mathematical formulations found in Eqs. (5) and (6). To assess these variations quantitatively, using Eqs. (5) and (6), Figs. 3 and 4 are plotted. Comparing Figs. 3 and 4, we conclude that the proportional increases in the source power against N, are more rapid than those of the beam size. It is to be understood from Figs. 3 and 4 that, to arrive at meaningful results for the receiver plane parameters, each FT beam must be normalized with respect to its own source plane parameters, rather than the normalization with respect to the Gaussian case for which N is unity.

According to Fig. 5, ${M}_{x}^{2}$ factor will become larger as the flatness order increases. This finding is consistent with the results in Ref. 10.

#### 3.2. Receiver plane

Throughout the receiver plane illustrations, we have utilized the wavelength and structure constant of *λ*=1.55 µm, ${C}_{n}^{2}$=10^{-15} m^{-2/3}. Therefore, these parameters are not quoted in our figures. By taking Eq. (14) and setting the source sizes and the flatness order as stated in the inset of the figure, Fig. 6 provides the views of the FT beam at selected propagation distances. Figure 6 shows that, the FT beam will initially develop a circular ring in the center. As the propagation advances, the circumference of this ring will become narrower, while a downward peak will gradually emerge from the center of the beam. Eventually this peak will smooth out the initial ring, with the profile turning into a pure Gaussian shape. This phenomenon may be considered to be akin to acquiring the spot of Arago as described in [19]. This type behavior of the FT beam is also commonly observed in other studies of free space propagation [4, 7, 8, 9, and 11]. For the FT beam of Fig. 6, the overlaid views of the beam progress cut along the slanted axis are displayed in Fig. 7. Introducing the Fresnel parameter *F*=${\alpha}_{s}^{2}$/(*λL*) with *α*_{s}
=*α*_{sx}
=*α*_{sy}
being the source size for a symmetrical source, in the plots of Figs. 6 and 7, the propagation distances of *L*=0.7 km, 1.4 km and 5 km correspond to *F* of 0.83 (i.e., close to Fresnel diffraction region meaning near field), 0.42 and 0.12 (i.e., in the Fraunhoffer diffraction region meaning far field), respectively.

Upon choosing asymmetric source sizes, an elliptical source beam is accomplished [8]. Further we can create an FT beam with sharp edges by exaggerating this asymmetry feature. The receiver plane intensity of such a beam is offered in Fig. 8, along with the associated source plane intensity. Here, the effect of the propagation channel implementing a Fourier transform action on the source beam field, thus converting it into an Airy function like shape in intensity during propagation, is slightly visible from the side lobes present in the lower picture of Fig. 8. Here we note that the most general case of elliptical flat-topped beam whose profile variables x and y cannot be separated is characterized by *α*_{sx}
, *α*_{sy}
and *α*_{sxy}
as shown in Ref. [9]. The propagation of such astigmatic flat-topped beam in the turbulent atmosphere is recently reported by Cai and He [20, 21]. The elliptical flat-topped beam of our derivation is obtained by setting *α*_{sx}
≠*α*_{sy}
.

Next we explore the beam size and power in bucket variations along the propagation axis. The beam size variation is shown in Fig. 9 against the flatness order, where the Gaussian source sizes and propagation length are adjusted as given in the same figure Here for better assessment, based on Eqs. (5) and (15), we have resorted to the relative beam size definition given by

Figure 9 reveals that for the selected range of source and propagation parameters, the relative beam size drops sharply at lower values of *N*, while this trend later continues at much slower pace. Turning to power in bucket variation and choosing a receiving aperture size of *α*_{r}
=0.5(${\alpha}_{\mathit{\text{sx}}}^{2}$+${\alpha}_{\mathit{\text{sy}}}^{2}$)1/2 in Eq. (17), we come across a similar picture, as illustrated in Fig. 10, where the falling trend at small *N* values is slightly retarded. Of course we should be careful here in the interpretation of receiving less power with increasing *N*. Less power occurs because, as explained underneath Eq. (17), *P*_{αN}
is determined by dividing the received power of each individual FT beam by its own source power. If this normalization were to be carried out with respect the Gaussian case, i.e., *N*=1, then we would come across rising power in bucket values with increasing *N*, instead of the presently falling trend.

Finally in this section, we analyze the kurtosis parameter which is known to be linked to the flatness or sharpness of the beam [18]. Figure 11 exhibits the variation of kurtosis parameter versus the propagation length at selected *N* values and fixed Gaussian source sizes. It is to be noted from Figs. (1) and (2) that on the source plane, the beam profile becomes sharper with larger *N*. This fact is reflected in Fig. 11 such that at short propagation distances, kurtosis parameter becomes smaller as *N* gets bigger. On the other hand, kurtosis parameter of all beams, except for *N*=1, begins to rise at medium propagation ranges, later reaching a maximum. This is a predictable behavior, where the rise corresponds to the formation of the ring, the subsequent narrowing action and then the appearance of the peak in Figs 6 and 7. When these curly surfaces on the beam start to erode with the beam evolving towards a Gaussian profile, the kurtosis parameter declines, finally approximating to the *N*=1 (Gaussian) case for all *N* values, as seen to the right of Fig. 11.

## 4. Conclusion

The source and receiver plane characteristics of flat topped (FT) beam propagating in turbulent atmosphere are studied. Source beam intensity, source beam size and power, *M*
^{2} factor are formulated and plotted against flatness order. From these plots, it is seen that beam intensity profile will become flatter, while source size, beam power and *M*
^{2} factor will increase against growing value of flatness order. Via the use of Huygens Fresnel integral, the average intensity on the receiver plane for a turbulent propagation environment is derived. Based on this intensity expression, the kurtosis parameter, beam size variation along the propagation axis and power in bucket formulations are stated. Plotting the receiver average intensity at selected source and propagation parameters, it is observed that, initially a circular ring will be formed in the center. As we go along the propagation axis however, the circumference of this ring will become narrower, while a downward peak will gradually emerge from the center of the beam. Eventually this peak will smooth out the initial ring, with the profile converging towards a pure Gaussian shape. These beam profile changes are also reflected in the plots of kurtosis parameter against propagation length. Furthermore, when the source beam is made highly asymmetric, the conversion of receiver intensity profile into an Airy function becomes more visible. The relative beam size becomes smaller, as the order parameter is raised, indicating that flatter beams will be subjected to less spreading during propagation in turbulence. For a fixed receiver aperture radius, less power is captured against the rising flatness order.

## References and links

**1. **F. M. Dickey and S. C. Holswade, *Laser beam shaping: theory and techniques* (Marcel Dekker, New York, 2000). [CrossRef]

**2. **D. L. Shealy and J. A. HoffnagleF. M. Dickey and D. L. Shealy, “Beam shaping profiles and propagation,” in *Laser Beam Shaping VI*, eds., Proc. SPIE5876, 1–11 (2005).

**3. **Y. Li, “Light beams with flat-topped profiles,” Opt. Lett. **27**, 1007–1009 (2002). [CrossRef]

**4. **Y. Li, “New expressions for flat-topped light beams,” Opt. Commun. **206**, 225–234 (2002). [CrossRef]

**5. **F. Gori, “Flattened Gaussian beams,” Opt. Commun. **107**, 335–341 (1994). [CrossRef]

**6. **Y. Baykal and H. T. Eyyuboğlu, “Scintillation index of flat-topped-Gaussian beams,” Appl. Opt.45 (2006). [CrossRef] [PubMed]

**7. **X. Ji and B. Lü, “Propagation of a flattened Gaussian beam through multi-apertured optical *ABCD* systems,” Optik **114**, 394–400 (2003). [CrossRef]

**8. **Y. Cai and Q. Lin, “Light beams with elliptical flat-topped profiles,” J. Opt. A: Pure Appl. Opt. **6**, 390–395 (2004). [CrossRef]

**9. **Y. Cai and Q. Lin, “A partially coherent elliptical flattened Gaussian beam and its propagation,” J. Opt. A: Pure Appl. Opt. **6**, 1061–1066 (2004). [CrossRef]

**10. **V. Bagini, R. Borghi, F. Gori, M. Santarsiero, D. Ambrosini, and G. S. Spagnolo, “Propagation of axially symmetric flattened Gaussian beams,” J. Opt. Soc. Am. A **13**, 1385–1394 (1996). [CrossRef]

**11. **N. Zhou, G. Zeng, and L. Hu, “Algorithms for flattened Gaussian beams passing through apertured and unapertured paraxial *ABCD* optical systems,” Opt. Commun. **240**, 299–306 (2004). [CrossRef]

**12. **B. Lü and S. Lou, “General propagation equation of flattened Gaussian beams,” J. Opt. Soc. Am. A. **17**, 2001–2004 (2000). [CrossRef]

**13. **B. Lü and H. Ma, “Coherent and incoherent off-axial Hermite-Gaussian beam combinations,” Appl. Opt. **39**, 1279–1289 (2000). [CrossRef]

**14. **D. Ge, Y. Cai, and Q. Lin, “Partially coherent flat-topped beam and its propagation,” Appl. Opt. **43**, 4732–4738 (2004). [CrossRef] [PubMed]

**15. **J. Zhang and Y. LiD. Lu and G. G. Matvienko, “Atmospherically turbulent effects on partially coherent flat-topped Gaussian beam,” in *Optical Technologies for Atmospheric, and Environmental Studies*, eds., Proc. SPIE5832, 48–55 (2005). [CrossRef]

**16. **W.H. Carter, “Spot size and divergence for Hermite Gaussian beams of any order,”Appl. Opt. **19**, 1027–1029 (1980). [CrossRef] [PubMed]

**17. **H. T. Eyyuboğlu and Y. Baykal, “Hermite-sine-Gaussian and Hermite-sinh-Gaussian laser beams in turbulent atmosphere,” J. Opt. Soc. Am. A **22**, 2709–2718 (2005). [CrossRef]

**18. **H. Mao, D. Zhao, F. Jing, and H. Liu, “Propagation characteristics of the kurtosis parameters of flat-topped beams passing through fractional Fourier transformation systems with a spherically aberrated lens,” J. Opt. A: Pure Appl. Opt. **6**, 640–650 (2004). [CrossRef]

**19. **J. E. Harvey and J. L. Forgham, “The spot of Arago: New relevance for an old phenomenon,” Am. J. Phys. **52**, 243–247 (1984). [CrossRef]

**20. **Y. Cai and S. He “Average intensity and spreading of an elliptical Gaussian beam in a turbulent atmosphere,” Opt. Lett. **31**, 568–570 (2006). [CrossRef] [PubMed]

**21. **Y. Cai and S. He “Propagation of various dark hollow beams in a turbulent atmosphere,” Opt. Express **14**, 1353–1367 (2006). [CrossRef] [PubMed]