OSA's Digital Library

Optics Express

Optics Express

  • Editor: Martijn de Sterke
  • Vol. 16, Iss. 20 — Sep. 29, 2008
  • pp: 15563–15575
« Show journal navigation

Beam propagation factor of partially coherent flat-topped beams in a turbulent atmosphere

Youquan Dan and Bin Zhang  »View Author Affiliations


Optics Express, Vol. 16, Issue 20, pp. 15563-15575 (2008)
http://dx.doi.org/10.1364/OE.16.015563


View Full Text Article

Acrobat PDF (378 KB)





Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

The Wigner distribution function (WDF) has been used to study the beam propagation factor (M2-factor) for partially coherent flat-topped (PCFT) beams with circular symmetry in a turbulent atmosphere. Based on the extended Huygens–Fresnel principle and the definition of the WDF, an expression for the WDF of PCFT beams in turbulence has been given. By use of the second-order moments of the WDF, the analytical formulas for the root-mean-square (rms) spatial width, the rms angular width, and the M2-factor of PCFT beams in turbulence have been derived, which can be applied to cases of different spatial power spectra of the refractive index fluctuations. The rms angular width and the M2-factor of PCFT beams in turbulence have been discussed with numerical examples. It can be shown that the M2-factor of PCFT beams in turbulence depends on the beam order, degree of global coherence of the source, waist width, wavelength, spatial power spectrum of the refractive index fluctuations, and propagation distance.

© 2008 Optical Society of America

1. Introduction

The propagation of light beams in atmospheric turbulence has been extensively studied for a long time [1

A. Ishimaru, “Theory and application of wave propagation and scattering in random media,” Proc. IEEE 65, 1030–1061 (1977). [CrossRef]

-3

L. C. Andrews and R. L. Phillips, Laser Beam Propagation through Random Media , 2nd ed. (SPIE, Bellingham, 2005). [CrossRef]

]. In recent years, because of many practical applications such as free-space optical communications, it is especially interesting that what types of light beams and under what conditions the effects of the turbulent atmosphere on the beams are the smallest [4

J. Wu, “Propagation of a Gaussian-Schell beam through turbulent media,” J. Mod. Opt. 37, 671–684 (1990). [CrossRef]

-8

T. Shirai, A. Dogariu, and E. Wolf, “Mode analysis of spreading of partially coherent beams propagating through atmospheric turbulence,” J. Opt. Soc. Am. A 20, 1094–1102 (2003). [CrossRef]

]; i.e., variations of the parameters (such as beam width etc.) characterizing beams due to turbulence are the smallest. Preliminary works [4

J. Wu, “Propagation of a Gaussian-Schell beam through turbulent media,” J. Mod. Opt. 37, 671–684 (1990). [CrossRef]

, 5

J. Wu and A. D. Boardman, “Coherence length of a Gaussian-Schell beam and atmospheric turbulence,” J. Mod. Opt. 38, 1355–1363 (1991). [CrossRef]

] have shown that the Gaussian-Schell model (GSM) beam with the worse coherence is less affected by turbulence. Consequently, lots of theoretical studies on the spatial and angular spreading of partially coherent beams in atmospheric turbulence have been carried out, indicating that partially coherent beams are less influenced by turbulence than fully coherent beams [6

G. Gbur and E. Wolf, “Spreading of partially coherent beams in random media,” J. Opt. Soc. Am. A 19, 1592–1598 (2002). [CrossRef]

-8

T. Shirai, A. Dogariu, and E. Wolf, “Mode analysis of spreading of partially coherent beams propagating through atmospheric turbulence,” J. Opt. Soc. Am. A 20, 1094–1102 (2003). [CrossRef]

]. Subsequently, experiments supporting these findings have been performed [9

A. Dogariu and S. Amarande, “Propagation of partially coherent beams: turbulence-induced degradation,” Opt. Lett. 28, 10–12 (2003). [CrossRef] [PubMed]

]. It is well known that the Wigner distribution function (WDF) can characterize partially coherent beams in space and in spatial frequency domain simultaneously [10

A. Walther, “Radiometry and coherence,” J. Opt. Soc. Am. 58, 1256–1259 (1968). [CrossRef]

-13

R. Martínez-Herrero, G. Piquero, and P. M. Mejías, “On the propagation of the kurtosis parameter of general beams,” Opt. Commun. 115, 225–232 (1995). [CrossRef]

]. Consequently, the beam propagation factor (M 2-factor) for characterizing laser beams can be studied. The M 2-factor of partially coherent beams propagating in free space and in gain or absorbing media has been investigated [14

A. E. Siegman, “New developments in laser resonators,” Proc. SPIE. 1224, 2–14 (1990). [CrossRef]

-18

B. zhang, Q. Wen, and X. Guo, “Beam propagation factor of partially coherent beams in gain or absorbing media,” Optik 117, 123–127 (2006). [CrossRef]

]. Very recently, the M 2-factor of fully coherent Gaussian beams in atmospheric turbulence has been studied by numerical approach [19

M. H. Mahdieh, “Numerical approach to laser beam propagation through turbulent atmosphere and evaluation of beam quality factor,” Opt. Commun. 281, 3395–3402 (2008). [CrossRef]

]. However, to our knowledge, there have been no reports on the beam propagation factor of partially coherent beams propagating in atmospheric turbulence by analytical approach.

The main purpose of this paper is to generalize the definition of M 2-factor in terms of second-order moments of the WDF to the case in atmospheric turbulence and to study the propagation law of M 2-factor for partially coherent flat-topped (PCFT) beams with circular symmetry in a turbulent atmosphere. First, by use of the extended Huygens–Fresnel principle and the definition of the WDF, an expression for the WDF of PCFT beams in a turbulent atmosphere has been given. Then, according to method of second-order moments of the WDF, the analytical formulas for the rms spatial width, the rms angular width and the M 2-factor of PCFT beams in turbulence have been derived, which can be applied to cases of different types of spatial power spectra of the refractive index fluctuations. Finally, the rms angular width and the M 2-factor of PCFT beams in turbulence have been discussed with numerical examples.

2. WDF of PCFT beams in a turbulent atmosphere

According to the beam model proposed by Li for a flat-topped (FT) beam [20

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

, 21

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

] and assuming that FT beams of any order contain equal power and under the assumption of the quasi-monochromatic scalar field, the field of an FT beam with circular symmetry can be expressed as [22

Y. Dan, B. Zhang, and P. Pan, “Propagation of partially coherent flat-topped beams through a turbulent atmosphere,” J. Opt. Soc. Am. A 25, 2223–2231 (2008). [CrossRef]

]

EN (ρ)= m=1N αmexp [ (m pN ρ2 w02)]
(1)

where σ is the radial coordinate, w 0 is the waist width of the Gaussian beam, N is the order of the FT beam,

αm= ( 1) m+1 N! m! ( Nm)!
(2)

and

pN=2 m=1N m=1N αm α m m+ m.
(3)

For the sake of convenience, a constant is omitted in Eq. (1). For the case of N=1, Eqs. (1) and (3) are reduced to the field of a Gaussian beam and p 1=1, respectively.

The fully coherent beam can be extended to the partially coherent beam by introducing a Gaussian coherence function, and this type of partially coherent beam can be produced by the fully coherent beam passing through a random phase plate or a liquid crystal [23

X. Ji, X. Chen, S. Chen, X. Li, and B. Lü, “Influence of atmospheric turbulence on the spatial correlation properties of partially coherent flat-topped beams,” J. Opt. Soc. Am. A 24, 3554–3563 (2007). [CrossRef]

, 24

M. Zahid and M. S. Zubairy, “Directionality of partially coherent Bessel-Gauss beams,” Opt. Commun. 70, 361–364 (1989). [CrossRef]

]. Consequently, based on the beam model described by Eqs. (1)-(3) for FT beams and the Schell model used for describing the partially coherent light [25

A. Starikov and E. Wolf, “Coherent-mode representation of Gaussian Schell-model sources and of their radiation fields,” J. Opt. Soc. Am. 72, 923–928 (1982). [CrossRef]

, 26

L. Mandel and E. Wolf, Optical Coherence and Quantum Optics (Cambridge University Press, Cambridge, 1995).

], the cross-spectral density function of PCFT beams at the plane of z=0 can be written as [22

Y. Dan, B. Zhang, and P. Pan, “Propagation of partially coherent flat-topped beams through a turbulent atmosphere,” J. Opt. Soc. Am. A 25, 2223–2231 (2008). [CrossRef]

, 23

X. Ji, X. Chen, S. Chen, X. Li, and B. Lü, “Influence of atmospheric turbulence on the spatial correlation properties of partially coherent flat-topped beams,” J. Opt. Soc. Am. A 24, 3554–3563 (2007). [CrossRef]

]

Γ ( ρ1, ρ2,0)=<E ( ρ1,0) E* ( ρ2,0) >m

= m=1N m=1N αm α mexp [ (m pN ρ12 w02+ m pN ρ22 w02+ ρ1 ρ22 2 σ02)],
(4)

where the angle brackets with subscript m denote averaging over the field ensemble, ρ 1 and ρ 2 are two vectors defining points on the transverse plane of z=0, w 0 represents the waist width of the Gaussian beam, σ 0 is the correlation length of the source. For the case of N=1, Eq. (4) is reduced to the cross-spectral density function of GSM beams.

Let us consider a PCFT beam propagating in a turbulent atmosphere from the plane z=0 into the half-space z>0. By using the paraxial form of the extended Huygens–Fresnel principle [3

L. C. Andrews and R. L. Phillips, Laser Beam Propagation through Random Media , 2nd ed. (SPIE, Bellingham, 2005). [CrossRef]

], the cross-spectral density function of PCFT beams through the turbulence can be expressed as [2

A. Ishimaru, Wave Propagation and Scattering in Random Media (Academic, New York, 1978).

, 3

L. C. Andrews and R. L. Phillips, Laser Beam Propagation through Random Media , 2nd ed. (SPIE, Bellingham, 2005). [CrossRef]

]

W (ρ, ρd,z)= ( k 2πz)2W (ρ, ρd,0)

×exp { ikz [ ( ρρ)· ( ρd ρd)]H ( ρd, ρd,z)} d2ρ d2 ρd,
(5)

where z is the propagation distance, k=2π/λ is the wave number, λ is the wavelength, and

W (ρ, ρd,0)= Γ ( ρ1, ρ2,0)=Γ ( ρ+ ρd2, ρ ρd2,0)
(6)

where the central abscissa coordinate systems are chosen, that is,

ρ= ( ρ1+ ρ2)2,  ρd= ρ1 ρ2,
(7)
ρ= ( ρ1+ ρ2)2,  ρd= ρ1 ρ2.
(8)

In Eq. (5) the term exp[-H(ρ d , ρ d , z)] represents the effect of the turbulence, and H can be written as

H ( ρd, ρd,z)=4 π2 k2z 01 0 [ 1 J0 (κ ρdξ+ ( 1ξ) ρd)] Φn (κ)κdκ,
(9)

where J0 is the Bessel function of the first kind and zero order, Φ n is the spatial power spectrum of the refractive index fluctuations of the turbulent atmosphere.

It is well known that the WDF is especially suitable for the treatment of partially coherent beams. The WDF can be expressed in terms of the cross-spectral density W(ρ, ρ d ,z) as [11

M. J. Bastiaans, “Application of the Wigner distribution function to partially coherent light,” J. Opt. Soc. Am. A 3, 1227–1238 (1986). [CrossRef]

, 12

J. Serna, R. Martínez-Herrero, and P. M. Mejías, “Parametric characterization of general partially coherent beams propagating through ABCD optical systems,” J. Opt. Soc. Am. A 8, 1094–1098 (1991). [CrossRef]

]

h (ρ,θ,z)= ( k 2π)2 W (ρ, ρd,z)exp ( ikθ· ρd) d2 ρd
(10)

where vector θ=(θ x ,θ y ), k θ x and y is the wave vector component along the x-axis and y-axis, respectively. Hence θ=(θ 2 x +θ 2 y )1/2 represents an angle of propagation (without taking the evanescent waves into account). Integration of function h over the angular variables θ x and θ y gives the beam intensity, and its integral over the spatial variables x and y is proportional to the radiant intensity of the field.

On the basis of inverse Fourier transform of the Dirac delta function and its property of even function (Chap. 5 of Ref. [27

R. N. Bracewell, The Fourier Transform and Its Applications , 2nd ed. (McGraw-Hill, New York, 1986).

]), we obtain

δ ( ρρ)= 1 (2π)2 exp [ ±i κd· ( ρρ)] d2 κd.
(11)

Then, the cross-spectral density of the beams in the source plane (z=0) can be rewritten as

W (ρ, ρd,0) 1 (2π)2 W (ρ, ρd,0)exp [ i κd· ( ρρ)] d2 κd d2ρ.
(12)

On substituting from Eq. (12) into Eq. (5) and using Eq. (11) and the property of the Dirac delta function (p. 75 of Ref. [27

R. N. Bracewell, The Fourier Transform and Its Applications , 2nd ed. (McGraw-Hill, New York, 1986).

])

f(x)δ (x)d x=f (0),
(13)

we obtain

W (ρ, ρd,z)= 1 (2π)2W (ρ, ρd + zk κd,0)

×exp [ i ρ· κd+ i κd· ρH( ρd, ρd + zk κd,z)] d2 κd d2ρ
(14)

where the term W(ρ, ρd + zk κd,0) is given by the expression from Eqs. (4) and (6), i.e.,

W(ρ, ρd + zk κd,0) = m=1N m=1N αm α mexp ( { pN w02 (m+m) ρ2

+ pN w02( m m)ρ· ( ρd+ zk κd)

+[ pN 4 w02 ( m+ m)+ 1 2 σ02] ρd+ zk κd2 }).
(15)

On substituting from Eq. (14) into Eq. (10) and calculating the integral with respect to ρ by using the formula (3.3232) in Ref. [31

I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series, and Products , 5th ed. (Academic, New York, 1980).

]

exp ( s2 x2±qx) dx= πsexp ( q2 4 s2),  ( s>0)
(16)

the WDF of PCFT beams through a turbulent atmosphere can be expressed as

h (ρ,θ,z)= w02 k2 16 π3 pN m=1N m=1N [ αm α m m+ m

× exp ( a ρd2 b κd2+c ρd · κd iρ· κd ik θ· ρdH) d2 κd d2 ρd],
(17)

where κ d =|κ d |=(κ 2 dx +κ 2 dy )1/2, and

H=4 π2 k2z 01dξ 0 [ 1 J0( κ zk κdξ+ ρd)] Φn (κ)κdκ,
(18)
a= 1 w02 ( mm pN m+ m+ 1 2 α2),
(19)
b= w02 4 ( m+ m) pN+ z2 k2 w02 ( mm pN m+ m+ 1 2 α2) iz 2k ( m m m+ m),
(20)
c= 2z k w02 ( mm pN m+ m+ 1 2 α2) + i2 ( m m m+ m),
(21)

where

α= σ0 w0
(22)

denotes the degree of global coherence of the source.

3. The rms spatial and angular spreads of PCFT beams in a turbulent atmosphere

The moments of the order n 1+n 2+m 1+m 2 of WDF h for three-dimensional beams are given by the expression [12

J. Serna, R. Martínez-Herrero, and P. M. Mejías, “Parametric characterization of general partially coherent beams propagating through ABCD optical systems,” J. Opt. Soc. Am. A 8, 1094–1098 (1991). [CrossRef]

, 13

R. Martínez-Herrero, G. Piquero, and P. M. Mejías, “On the propagation of the kurtosis parameter of general beams,” Opt. Commun. 115, 225–232 (1995). [CrossRef]

]

< x n1 y n2 θx m1 θy m2>= 1P x n1 y n2 θx m1 θy m2h (ρ,θ,z) d2ρ d2θ
(23)

where P=∬h(ρ,θ,z)d2 ρd2 θ is the total irradiance of the beams, and P=πw 2 0/2 for the case of PCFT beams.

As far as the second-order field correlations are concerned, a partially coherent beam can be characterized by the WDF. It can be seen from Eq. (23) that the moments at the plane of z=const are determined by the WDF at the same plane. Consequently, using Eqs. (17) and (23), we can obtain the moments of PCFT beams at the plane of z=const in turbulence.

On substituting from Eq. (17) into Eq. (23), and recalling the formulas (p. 82, p. 95, and p. 141 of Ref. [27

R. N. Bracewell, The Fourier Transform and Its Applications , 2nd ed. (McGraw-Hill, New York, 1986).

])

δ (s)= 1 2π exp ( isx)dx,
(24)
δ (n) (s)= 1 2π ( ix)nexp ( isx)dx, ( n=1,2)
(25)
f(x)δ (n) (x)dx= ( 1)n f (n) ( 0), ( n=1,2)
(26)

(where δ (n) is nth derivative of the Dirac delta function, f(x) is an arbitrary function and f (n)(x) is its nth derivative), some second-order moments of PCFT beams in turbulence turn out to be

< ρ2>= < x2>+< y2>

= 2 w02 pN2 [ m=1N m=1N αm α m ( m+ m)2] + 2 k2 w02[ 1 α2

+4 m=1N m=1N mm αm α m (m + m)2] z2 + 43 π2T z3,
(27)
< θ2>= < θx2>+< θy2>

= 2 k2 w02 [ 1 α2+4 m=1N m=1N mm αm α m ( m+ m)2]+4 π2Tz,
(28)
< ρ·θ>= < x θx>+< y θy>

= 2 k2 w02 [ 1 α2+4 m=1N m=1N mm αm α m ( m+ m)2]z+2 π2 Tz2,
(29)

whereas the other second-order moments (i.e., <xy>, <θ x θ y >, < y > and < x >) and all the first-order moments turn out to be zero. In Eqs. (27)-(29), the quantity T is given by the expression, i.e.,

T= 0 Φn (κ) κ3dκ
(30)

which denotes the intensity of the turbulent atmosphere.

It should be noted that the quantity T, which can only be determined by the spatial power spectrum of the refractive index fluctuations, is independent of the parameter of the beams.

According to the definition of rms spatial and angular width of light beams [12

J. Serna, R. Martínez-Herrero, and P. M. Mejías, “Parametric characterization of general partially coherent beams propagating through ABCD optical systems,” J. Opt. Soc. Am. A 8, 1094–1098 (1991). [CrossRef]

], the rms spatial and angular width of PCFT beams in turbulence can be expressed as

wN (z) ( < ρ <ρ>2>) 12 = ( < ρ2>) 12

={ 2 w02 pN2 [ m=1N m=1N αm α m ( m+ m)2] + 2 k2 w02[ 1 α2

+4 m=1N m=1N mm αm α m (m + m)2] z2+ 43 π2T z3 } 12,
(31)
θN (z) ( < θ <θ>2>) 12= ( < θ2>) 12

={ 2 k2 w02 [ 1 α2+4 m=1N m=1N mm αm α m ( m+ m)2]+4 π2Tz } 12.
(32)

Equation (31), which is same as Eq. (44) of Ref. [22

Y. Dan, B. Zhang, and P. Pan, “Propagation of partially coherent flat-topped beams through a turbulent atmosphere,” J. Opt. Soc. Am. A 25, 2223–2231 (2008). [CrossRef]

], represents the spatial spreading of PCFT beams in turbulence; and it is also consistent with the formula (30) of Ref. [6

G. Gbur and E. Wolf, “Spreading of partially coherent beams in random media,” J. Opt. Soc. Am. A 19, 1592–1598 (2002). [CrossRef]

] derived by method of the “equivalent source”. Equation (32) represents the angular spreading of PCFT beams in turbulence, and is consistent with the formula (25a) of Ref. [7

S. A. Ponomarenko, J. -J. Greffet, and E. Wolf, “The diffusion of partially coherent beams in turbulent media,” Opt. Commun. 208, 1–8 (2002). [CrossRef]

] derived by method of the “phase-space dynamics of beam”. The angular width given by Eq. (32) can be regarded as the width of the beams in spatial frequency domain. In other words, the angular width θ N (z ) given by Eq. (32) is equal to the divergence of the “source” at the plane of z=z′ in the free space, but the “source” is described by the WDF h(ρ,θ,z ) in turbulence, i.e., Eq. (17).

It can be readily seen from Eq. (32) that the rms angular spread of PCFT beams propagating in turbulence increases with the propagation distance z, whereas it maintains invariance in free space (i.e., for the case of T=0).

For the sake of comparison, the normalized rms angular width, which represents the relative angular spreading of the beams, is defined as

θ rN (z) θN (z) θN (0)

={ 1+2 π2 k2 w02[ 1 α2+4 m=1N m=1N mm αm α m ( m+ m)2 ] 1Tz } 12.
(33)

From Eq. (33), it is clear that the normalized rms angular spreading of PCFT beams in turbulence decreases with decreasing α; i.e., the beams with the worse coherence are less affected by turbulence. For a given propagation distance z, Eq. (33) indicates that the relative angular spreading increases with waist width w 0 and decreasing wavelength λ, whereas based on Eq. (31) the relative spatial spreading w N (z)/w N increases with decreasing w 0 and increasing λ. That is, the influence of atmospheric turbulence on the relative angular spreading is larger for PCFT beams with the larger w 0 and the smaller λ, but the influence on the relative spatial spreading is smaller. Consequently, in these cases neither the normalized rms spatial nor angular spreading can represents comprehensively the effect of turbulence on PCFT beams.

For the case of N=1, Eqs. (31)-(33) are easily simplified to the results of GSM beams, i.e.,

w1 (z)= [ w022+ 2 k2 w02 ( 1+ 1 α2) z2 + 43 π2T z3] 12,
(34)
θ1 (z)= [ 2 k2 w02 ( 1+ 1 α2) +4π2Tz] 12,
(35)
θ r1(z)= [ 1+2 π2 k2 w02 ( 1 α2+1) 1Tz] 12.
(36)

Taking into account the modified von Karman spectrum [3

L. C. Andrews and R. L. Phillips, Laser Beam Propagation through Random Media , 2nd ed. (SPIE, Bellingham, 2005). [CrossRef]

], Eq. (34) is equal to Eq. (13) in [28

X. Ji and B. Lü, “Turbulence-induced quality degradation of partially coherent beams,” Opt. Commun. 251, 231–236 (2005). [CrossRef]

]. For the case of α→∞, Eqs. (34)-(36) are further reduced to the expressions of fully coherent Gaussian beams in turbulence.

4. Beam propagation M2-factor of PCFT beams in a turbulent atmosphere

For three-dimensional beams, the beam propagation M 2-factor can be defined as [12

J. Serna, R. Martínez-Herrero, and P. M. Mejías, “Parametric characterization of general partially coherent beams propagating through ABCD optical systems,” J. Opt. Soc. Am. A 8, 1094–1098 (1991). [CrossRef]

]

M2 (z)= k( < ρ2> < θ2> < ρ·θ>2) 12

=k [ ( < x2>+ < y2>) ( < θx2>+ < θy2>) ( <x θx>+ <y θy>)2] 12.
(37)

It should be noted that in Eq. (37) the M 2-factor of a beam can be defined in terms of second-order moments of the WDF in the same plane, and it can also be calculated from the WDF h at any plane of z=const. Consequently, the definition of M 2-factor, i.e. Eq. (37), is valid for the case of the presence of the turbulence. In general, the value of M 2(z) in turbulence is equal to the one that is measured at the distance z.

On substituting from Eqs. (27)-(29) into Eq. (37), the beam propagation M 2-factor of PCFT beams through a turbulent atmosphere can be expressed as

MN2 (z)= [ MN4 (0)+ k2 w02ATz+ B w01 Tz3+ 43 π4 k2 T2 z4] 12
(38)

where T is given by Eq. (30), and

MN4 (0)= 3AB 16 π4,
(39)
A= 8 π2 pN2 m=1N m=1N αm α m ( m+ m)2,
(40)

and

B= 8 π23[ 1 α2 +4 m=1N m=1N mm αm α m ( m+ m)2].
(41)

Note that the value of M 2(z ) given by Eq. (38) is equal to the M 2-factor of the “source” at the plane of z=z′ in the free space, but the “source” is described by the WDF h(ρ,θ,z ) described by Eq. (17) in turbulence. Equation (38), which is the main result of this paper, indicates that the M 2-factor of PCFT beams in turbulence increases with the propagation distance z. For the case of T=0 (in free space), Eq. (38) shows that the M 2-factor, which is only determined by the beam order and coherent parameter of the source, maintains constant (cf. Refs. [29

G. Wu, H. Guo, and D. Deng, “Paraxial propagation of partially coherent flat-topped beam,” Opt. Commun. 260, 687–690 (2006). [CrossRef]

, 30

H. T. Eyyuboğlu, Ç. Arpali, and Y. Baykal, “Flat topped beams and their characteristics in turbulent media,” Opt. Express 14, 4196–4207 (2006). [CrossRef] [PubMed]

]).

For the sake of convenience, the normalized M 2-factor, which represents the relative M 2-factor of the beams, is introduced by the expression

M rN2 (z) MN2 (z) MN2 (0)

= ( 1+ 16 π4 k2 w02 3BTz+ 16 π4 3 w02A Tz3+ 64 π8 k2 9AB T2 z4) 12.
(42)

It can be shown from Eqs. (40) and (41) that in Eq. (42), the coefficient A only depends on N, the coefficient B only depends on N and α, and the value of B increases with decreasing α. Further, it can be readily seen from Eq. (42) that the normalized M 2-factor of PCFT beams decreases with decreasing α and increasing λ. Consequently, the effect of atmospheric turbulence on the normalized M 2-factor of PCFT beams for the worse coherence and the longer wavelength is smaller. It can be also seen from Eq. (42) that for a given propagation distance z, there exists a waist width w 0=w 0min corresponding to the minimal value of the normalized M 2-factor for the beams with different waist widths.

For the cases mentioned in section 3, neither the spatial nor angular spreading can denote comprehensively the effect of atmospheric turbulence on PCFT beams at the plane of z=const for different values of w 0 and of λ, but now we can use the normalized M 2-factor for these cases. Namely, in these cases the smaller value of the normalized M 2-factor of a beam indicates that the beam is less affected by turbulence. Therefore, the normalized M2-factor can show more comprehensively the effect of the turbulence on the optical beams than the normalized spatial and angular width.

For the case of N=1, Eqs. (38) and (42) are simplified to the M 2-factor and the normalized M 2-factor of GSM beams, respectively, namely

M12 (z)= [ ( 1+ 1 α2) +2 π2 k2 w02Tz

+ 8 π2 3 w02 ( 1+ 1 α2)T z3 + 4 π4 k23 T2 z4 ] 12,
(43)
M r12(z)=[ 1 +2 π2 k2 w02 ( 1+ 1 α2) 1Tz

+ 8 π2 3 w02T z3 + 4 π4 k23 ( 1+ 1 α2) 1 T2z4 ] 12.
(44)

The formula (43) is a generalization of known formula for the M 2-factor of GSM beams in free space [16

M. Santarsiero, F. Gori, R. Borghi, G. Cincotti, and P. Vahimaa, “Spreading properties of beams radiated by partially coherent Schell-model sources,” J. Opt. Soc. Am. A 16, 106–112 (1999). [CrossRef]

]. Further, for the case of α→∞, Eqs. (43) and (44) are reduced to the M 2-factor and the normalized M 2-factor of fully coherent Gaussian beams in turbulence respectively.

5. Numerical calculation results and analysis

In this section, the Tatarskii spectrum is used as a model for atmospheric turbulence [3

L. C. Andrews and R. L. Phillips, Laser Beam Propagation through Random Media , 2nd ed. (SPIE, Bellingham, 2005). [CrossRef]

], i.e.,

Φn (κ)= 0.033Cn2 κ 113exp ( κ2 κm2)
(45)

where C 2 n is the structure constant of the refractive index fluctuations of the turbulence and κ m =5.92/l 0 with l 0 being the inner scale turbulence.

On substituting from Eq. (45) into Eq. (30) and by using the formula 3.4621 in Ref. [31

I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series, and Products , 5th ed. (Academic, New York, 1980).

], the quantity T can be easily simplified to

T=0.1661 Cn2 l0 13
(46)

Numerical calculations for the rms angular width and normalized rms angular width of PCFT beams were performed on the basis of Eqs. (32), (33) and (46). Some illustrative examples are compiled in Figs. 1 and 2. The calculation parameters are α=2, λ=850 nm, w 0=0.05 m, C 2 n=10-15 m-2/3.

Figure 1 gives the rms angular and normalized rms angular widths for different values of beam order N versus the propagating distance in atmospheric turbulence, respectively. It can be seen from Fig. 1(a) that differing from the case of free space (absence of the turbulence), the rms angular width in turbulence obviously increases with the propagation distance, and it is larger for the beams with a higher beam order. Figure 1(b) implies that the normalized rms angular width in turbulence is smaller for the beams with a higher beam order; i.e., the relative angular spreading is less affected by atmospheric turbulence for the beams with a higher beam order.

Fig. 1. (a). rms angular width and (b) normalized rms angular width versus propagation distance for different values of beam order N, respectively. The calculation parameters are α=2, λ=850 nm, w 0=0.05 m, C 2 n=10-15 m-2/3.

Figure 2 shows the variation of the normalized rms angular width for different values of coherent parameter α on the propagating distance in turbulence. It follows from Fig. 2 that the relative angular spreading of PCFT beams in turbulence decreases with decreasing α, meaning that the relative angular spreading for the beams with the worse spatial coherence is less affected by turbulence.

Fig. 2. Normalized rms angular width versus propagation distance for different values of α. The calculation parameters are N=10, λ=850 nm, w 0=0.05 m, C 2 n=10-15 m-2/3.

Moreover, based on Eqs. (38), (42) and (46), numerical calculations for the M 2-factor and normalized M 2-factor of PCFT beams in turbulence were carried out. Some typical examples are compiled in Figs. 3-6.

Figure 3 shows the M 2-factor and normalized M 2-factor for different values of beam order N versus the propagating distance in atmospheric turbulence, respectively. The calculation parameters are α=2, λ=850 nm, w 0=0.05 m, C 2 n=10-15 m-2/3. It can be shown in Fig. 3(a) that the M 2-factor in turbulence apparently increases with the propagation distance, indicating that the beam quality obviously becomes worse with increasing propagation distance in turbulence. Figure 3(b) implies that the normalized M 2-factor in turbulence is smaller for PCFT beams with a higher beam order, showing the less influence of the turbulence on the relative M 2-factor of the beams with a higher beam order.

Fig. 3. (a). M 2-factor and (b) normalized M 2-factor versus propagation distance for different values of beam order N, respectively. The calculation parameters are α=2, λ=850 nm, w 0=0.05 m, C 2 n=10-15 m-2/3.

Figure 4 gives the normalized M 2-factor for different values of coherent parameter α versus the propagating distance in atmospheric turbulence. The calculation parameters are N=10, λ=850 nm, w 0=0.05 m, C 2 n=10-15 m-2/3. It can be easily seen that the normalized M 2-factor of PCFT beams in turbulence decreases with decreasing α, showing the less influence of the turbulence on the relative M 2-factor of the beams with the worse spatial coherence.

Fig. 4. Normalized M 2-factor versus propagation distance for different values of α. The calculation parameters are N=10, λ=850 nm, w 0=0.05 m, C 2 n=10-15 m-2/3.

Figure 5 gives the normalized M 2-factor for different values of beam order N at the plane of z=1 km versus the waist width w 0 of the GSM beam. The calculation parameters are α=2, z=1 km, C 2 n=10-15 m-2/3, λ=850 nm. It can also be seen from Fig. 5 that there exists a minimum of the normalized M 2-factor as the waist width w 0 increases. For example, for the case of N=1 (i.e. GSM beams), the minimal value of the normalized M 2-factor is of about 1.11, whereas for the case of N=10, the minimal value of the normalized M 2-factor is of about 1.10.

Fig. 5. Normalized M 2-factor versus w 0 for different values of beam order N. The calculation parameters are α=2, z=1 km, C 2 n=10-15 m-2/3, λ=850 nm.

Figure 6 represents the normalized M 2-factor for different values of beam order N at the plane of z=1 km versus the wavelength of PCFT beams. The calculation parameters are α=2, z=1 km, C 2 n=10-15 m-2/3, w 0=0.05 m. It can be clearly seen from Fig. 6 that the normalized M 2-factor decreases with increasing wavelength λ, indicating that the relative M 2-factor of the beams with the longer wavelength is less affected by turbulence.

Fig. 6. Normalized M 2-factor versus wave length λ for different values of beam order N. The calculation parameters are α=2, z=1 km, C 2 n=10-15 m-2/3, w 0=0.05 m.

6. Conclusions

In the present paper, by use of the extended Huygens–Fresnel principle and definition of the WDF, an expression for the WDF of PCFT beams with circular symmetry in a turbulent atmosphere has been given. The analytical formulas for the rms spatial width, the rms angular width, and the beam propagation M 2-factor of PCFT beams in turbulence have been derived based on the method of second-order moments of the WDF. For the case of beam order N=1, all results obtained in this paper can be simplified to those of GSM beams. Furthermore, the rms angular width and M 2-factor in turbulence have been discussed with numerical examples.

It can be shown that the beam propagation M 2-factor of PCFT beams depends on the beam order, the degree of global coherence of the source, the waist width, the wavelength, and the spatial power spectrum of the refractive index fluctuations of atmospheric turbulence. Moreover, the beam propagation M 2-factor of PCFT beams increases obviously with the propagation distance in turbulence. From the viewpoint of relative variation of the M 2-factor due to the turbulent atmosphere, the effect of the turbulence on PCFT beams is smaller for beams with a higher order, a lower degree of global coherence of the source, and a larger wavelength.

Acknowledgment

This research was supported by the Program for New Century Excellent Talents in University (NCET-05-0784)

References and links

1.

A. Ishimaru, “Theory and application of wave propagation and scattering in random media,” Proc. IEEE 65, 1030–1061 (1977). [CrossRef]

2.

A. Ishimaru, Wave Propagation and Scattering in Random Media (Academic, New York, 1978).

3.

L. C. Andrews and R. L. Phillips, Laser Beam Propagation through Random Media , 2nd ed. (SPIE, Bellingham, 2005). [CrossRef]

4.

J. Wu, “Propagation of a Gaussian-Schell beam through turbulent media,” J. Mod. Opt. 37, 671–684 (1990). [CrossRef]

5.

J. Wu and A. D. Boardman, “Coherence length of a Gaussian-Schell beam and atmospheric turbulence,” J. Mod. Opt. 38, 1355–1363 (1991). [CrossRef]

6.

G. Gbur and E. Wolf, “Spreading of partially coherent beams in random media,” J. Opt. Soc. Am. A 19, 1592–1598 (2002). [CrossRef]

7.

S. A. Ponomarenko, J. -J. Greffet, and E. Wolf, “The diffusion of partially coherent beams in turbulent media,” Opt. Commun. 208, 1–8 (2002). [CrossRef]

8.

T. Shirai, A. Dogariu, and E. Wolf, “Mode analysis of spreading of partially coherent beams propagating through atmospheric turbulence,” J. Opt. Soc. Am. A 20, 1094–1102 (2003). [CrossRef]

9.

A. Dogariu and S. Amarande, “Propagation of partially coherent beams: turbulence-induced degradation,” Opt. Lett. 28, 10–12 (2003). [CrossRef] [PubMed]

10.

A. Walther, “Radiometry and coherence,” J. Opt. Soc. Am. 58, 1256–1259 (1968). [CrossRef]

11.

M. J. Bastiaans, “Application of the Wigner distribution function to partially coherent light,” J. Opt. Soc. Am. A 3, 1227–1238 (1986). [CrossRef]

12.

J. Serna, R. Martínez-Herrero, and P. M. Mejías, “Parametric characterization of general partially coherent beams propagating through ABCD optical systems,” J. Opt. Soc. Am. A 8, 1094–1098 (1991). [CrossRef]

13.

R. Martínez-Herrero, G. Piquero, and P. M. Mejías, “On the propagation of the kurtosis parameter of general beams,” Opt. Commun. 115, 225–232 (1995). [CrossRef]

14.

A. E. Siegman, “New developments in laser resonators,” Proc. SPIE. 1224, 2–14 (1990). [CrossRef]

15.

F. Gori, M. Santarsiero, and A. Sona, “The change of width for a partially coherent beam on paraxial propagation,” Opt. Commun. 82, 197–203 (1991). [CrossRef]

16.

M. Santarsiero, F. Gori, R. Borghi, G. Cincotti, and P. Vahimaa, “Spreading properties of beams radiated by partially coherent Schell-model sources,” J. Opt. Soc. Am. A 16, 106–112 (1999). [CrossRef]

17.

C. Palma, P. De Santis, G. Cincotti, and G. Guattari, “Propagation and coherence evolution of optical beams in gain media,” J. Mod. Opt. 43, 139–153 (1996). [CrossRef]

18.

B. zhang, Q. Wen, and X. Guo, “Beam propagation factor of partially coherent beams in gain or absorbing media,” Optik 117, 123–127 (2006). [CrossRef]

19.

M. H. Mahdieh, “Numerical approach to laser beam propagation through turbulent atmosphere and evaluation of beam quality factor,” Opt. Commun. 281, 3395–3402 (2008). [CrossRef]

20.

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

21.

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

22.

Y. Dan, B. Zhang, and P. Pan, “Propagation of partially coherent flat-topped beams through a turbulent atmosphere,” J. Opt. Soc. Am. A 25, 2223–2231 (2008). [CrossRef]

23.

X. Ji, X. Chen, S. Chen, X. Li, and B. Lü, “Influence of atmospheric turbulence on the spatial correlation properties of partially coherent flat-topped beams,” J. Opt. Soc. Am. A 24, 3554–3563 (2007). [CrossRef]

24.

M. Zahid and M. S. Zubairy, “Directionality of partially coherent Bessel-Gauss beams,” Opt. Commun. 70, 361–364 (1989). [CrossRef]

25.

A. Starikov and E. Wolf, “Coherent-mode representation of Gaussian Schell-model sources and of their radiation fields,” J. Opt. Soc. Am. 72, 923–928 (1982). [CrossRef]

26.

L. Mandel and E. Wolf, Optical Coherence and Quantum Optics (Cambridge University Press, Cambridge, 1995).

27.

R. N. Bracewell, The Fourier Transform and Its Applications , 2nd ed. (McGraw-Hill, New York, 1986).

28.

X. Ji and B. Lü, “Turbulence-induced quality degradation of partially coherent beams,” Opt. Commun. 251, 231–236 (2005). [CrossRef]

29.

G. Wu, H. Guo, and D. Deng, “Paraxial propagation of partially coherent flat-topped beam,” Opt. Commun. 260, 687–690 (2006). [CrossRef]

30.

H. T. Eyyuboğlu, Ç. Arpali, and Y. Baykal, “Flat topped beams and their characteristics in turbulent media,” Opt. Express 14, 4196–4207 (2006). [CrossRef] [PubMed]

31.

I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series, and Products , 5th ed. (Academic, New York, 1980).

OCIS Codes
(010.1300) Atmospheric and oceanic optics : Atmospheric propagation
(010.1330) Atmospheric and oceanic optics : Atmospheric turbulence
(010.3310) Atmospheric and oceanic optics : Laser beam transmission
(030.0030) Coherence and statistical optics : Coherence and statistical optics

ToC Category:
Atmospheric and oceanic optics

History
Original Manuscript: June 26, 2008
Revised Manuscript: August 30, 2008
Manuscript Accepted: September 10, 2008
Published: September 17, 2008

Citation
Youquan Dan and Bin Zhang, "Beam propagation factor of partially coherent flat-topped beams in a turbulent atmosphere," Opt. Express 16, 15563-15575 (2008)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-16-20-15563


Sort:  Author  |  Year  |  Journal  |  Reset  

References

  1. A. Ishimaru, "Theory and application of wave propagation and scattering in random media," Proc. IEEE 65, 1030-1061 (1977). [CrossRef]
  2. A. Ishimaru, Wave Propagation and Scattering in Random Media (Academic, New York, 1978).
  3. L. C. Andrews and R. L. Phillips, Laser Beam Propagation through Random Media, 2nd ed. (SPIE, Bellingham, 2005). [CrossRef]
  4. J. Wu, "Propagation of a Gaussian-Schell beam through turbulent media," J. Mod. Opt. 37, 671-684 (1990). [CrossRef]
  5. J. Wu and A. D. Boardman, "Coherence length of a Gaussian-Schell beam and atmospheric turbulence," J. Mod. Opt. 38, 1355-1363 (1991). [CrossRef]
  6. G. Gbur and E. Wolf, "Spreading of partially coherent beams in random media," J. Opt. Soc. Am. A 19, 1592-1598 (2002). [CrossRef]
  7. S. A. Ponomarenko, J. -J. Greffet, and E. Wolf, "The diffusion of partially coherent beams in turbulent media," Opt. Commun. 208, 1-8 (2002). [CrossRef]
  8. T. Shirai, A. Dogariu, and E. Wolf, "Mode analysis of spreading of partially coherent beams propagating through atmospheric turbulence," J. Opt. Soc. Am. A 20, 1094-1102 (2003). [CrossRef]
  9. A. Dogariu and S. Amarande, "Propagation of partially coherent beams: turbulence-induced degradation," Opt. Lett. 28, 10-12 (2003). [CrossRef] [PubMed]
  10. A. Walther, "Radiometry and coherence," J. Opt. Soc. Am. 58, 1256-1259 (1968). [CrossRef]
  11. M. J. Bastiaans, "Application of the Wigner distribution function to partially coherent light," J. Opt. Soc. Am. A 3, 1227-1238 (1986). [CrossRef]
  12. J. Serna, R. Martínez-Herrero, and P. M. Mejías, "Parametric characterization of general partially coherent beams propagating through ABCD optical systems," J. Opt. Soc. Am. A 8, 1094-1098 (1991). [CrossRef]
  13. R. Martínez-Herrero, G. Piquero, and P. M. Mejías, "On the propagation of the kurtosis parameter of general beams," Opt. Commun. 115, 225-232 (1995). [CrossRef]
  14. A. E. Siegman, "New developments in laser resonators," Proc. SPIE. 1224, 2-14 (1990). [CrossRef]
  15. F. Gori, M. Santarsiero, and A. Sona, "The change of width for a partially coherent beam on paraxial propagation," Opt. Commun. 82, 197-203 (1991). [CrossRef]
  16. M. Santarsiero, F. Gori, R. Borghi, G. Cincotti, and P. Vahimaa, "Spreading properties of beams radiated by partially coherent Schell-model sources," J. Opt. Soc. Am. A 16, 106-112 (1999). [CrossRef]
  17. C. Palma, P. De Santis, G. Cincotti, and G. Guattari, "Propagation and coherence evolution of optical beams in gain media," J. Mod. Opt. 43, 139-153 (1996). [CrossRef]
  18. B. Zhang, Q. Wen, and X. Guo, "Beam propagation factor of partially coherent beams in gain or absorbing media," Optik 117, 123-127 (2006). [CrossRef]
  19. M. H. Mahdieh, "Numerical approach to laser beam propagation through turbulent atmosphere and evaluation of beam quality factor," Opt. Commun. 281, 3395-3402 (2008). [CrossRef]
  20. Y. Li, "Light beams with flat-topped profiles," Opt. Lett. 27, 1007-1009 (2002). [CrossRef]
  21. Y. Li, "New expressions for flat-topped light beams," Opt. Commun. 206, 225-234 (2002). [CrossRef]
  22. Y. Dan, B. Zhang, and P. Pan, "Propagation of partially coherent flat-topped beams through a turbulent atmosphere," J. Opt. Soc. Am. A 25, 2223-2231 (2008). [CrossRef]
  23. X. Ji, X. Chen, S. Chen, X. Li, and B. Lü, "Influence of atmospheric turbulence on the spatial correlation properties of partially coherent flat-topped beams," J. Opt. Soc. Am. A 24, 3554-3563 (2007). [CrossRef]
  24. M. Zahid and M. S. Zubairy, "Directionality of partially coherent Bessel-Gauss beams," Opt. Commun. 70, 361-364 (1989). [CrossRef]
  25. A. Starikov and E. Wolf, "Coherent-mode representation of Gaussian Schell-model sources and of their radiation fields," J. Opt. Soc. Am. 72, 923-928 (1982). [CrossRef]
  26. L. Mandel and E. Wolf, Optical Coherence and Quantum Optics (Cambridge University Press, Cambridge, 1995).
  27. R. N. Bracewell, The Fourier Transform and Its Applications, 2nd ed. (McGraw-Hill, New York, 1986).
  28. X. Ji and B. Lü, "Turbulence-induced quality degradation of partially coherent beams," Opt. Commun. 251, 231-236 (2005). [CrossRef]
  29. G. Wu, H. Guo, and D. Deng, "Paraxial propagation of partially coherent flat-topped beam," Opt. Commun. 260, 687-690 (2006). [CrossRef]
  30. H. T. Eyyubo?lu, C. Arpali, and Y. Baykal, "Flat topped beams and their characteristics in turbulent media," Opt. Express 14, 4196-4207 (2006). [CrossRef] [PubMed]
  31. I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series, and Products, 5th ed. (Academic, New York, 1980).

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.


« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited