OSA's Digital Library

Optics Express

Optics Express

  • Editor: C. Martijn de Sterke
  • Vol. 19, Iss. 19 — Sep. 12, 2011
  • pp: 18216–18228
« Show journal navigation

Mathematic model analysis of Gaussian beam propagation through an arbitrary thickness random phase screen

Yuzhen Tian, Jin Guo, Rui Wang, and Tingfeng Wang  »View Author Affiliations


Optics Express, Vol. 19, Issue 19, pp. 18216-18228 (2011)
http://dx.doi.org/10.1364/OE.19.018216


View Full Text Article

Acrobat PDF (1037 KB)





Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

In order to research the statistical properties of Gaussian beam propagation through an arbitrary thickness random phase screen for adaptive optics and laser communication application in the laboratory, we establish mathematic models of statistical quantities, which are based on the Rytov method and the thin phase screen model, involved in the propagation process. And the analytic results are developed for an arbitrary thickness phase screen based on the Kolmogorov power spectrum. The comparison between the arbitrary thickness phase screen and the thin phase screen shows that it is more suitable for our results to describe the generalized case, especially the scintillation index.

© 2011 OSA

1. Introduction

In this paper, we will, based on the Rytov weak fluctuation method and the thin phase screen model of Andrews et al. [5

5. L. C. Andrews, R. L. Phillips, and A. R. Weeks, “Propagation of a Gaussian-beam wave through a random phase screen,” Waves Random Media 7(2), 229–244 (1997). [CrossRef]

], research the statistical property of Gaussian beam propagation through an arbitrary thickness phase screen, establish the mathematic models of involving statistical quantities, and develop the analytical results. The application scope of our results will be discussed based on the comparison between the thin phase screen models and arbitrary thickness models.

2. Gaussian beam parameters

We assume the transmitted wave is a lowest-order Gaussian beam wave that, also called TEM00 wave, in the plane of the emitting aperture of the transmitter at z = 0, can be expressed as [4

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

,12

12. M. Gao and Z. Wu, “Experiments of effect of beam spreading of far-field on aiming deviation,” Opt. Precision Eng. 18(3), 602–608 (2010).

,13

13. Y.-h. Gao, Z.-y. An, N.-n. Li, W.-x Zhao, and J.-s. Wang, “Optical design of Gaussian beam shaping,” Opt. Precision Eng. 19(7), 1464–1471 (2011). [CrossRef]

]

z=0:U0(r,0)=a0exp(rW02ikr22F0)=a0exp(12α0kr2),
(1)

where r = (x 2 + y 2)1/2, α 0 = 2/(kW 0 2) + i/F 0; a 0 is the optical amplitude, F 0 and W 0 denote the phase front radius of curvature and beam radius at which the field amplitude falls to 1/e of that on the beam axis, respectively, and k is the optical wave number. The particular cases F 0 = ∞, F 0>0, and F 0<0 correspond to collimated, convergent, and divergent beam formats, respectively. Based on paraxial theory, the free space optical wave in the receiver plane at z = L is another Gaussian beam described by [4

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

,12

12. M. Gao and Z. Wu, “Experiments of effect of beam spreading of far-field on aiming deviation,” Opt. Precision Eng. 18(3), 602–608 (2010).

]

z=L:U0(r,L)=a0p(L)exp[ikL+(rW2ikr22F)]=a0p(L)exp[ikL12(α0kr2p(L))],
(2)

where W and F denote the receiver plane phase front radius of curvature and beam radius, respectively, and p(L) = 1 + 0 L is the propagation parameter.

Figure 1
Fig. 1 The model of Gaussian beam propagation through a random phase screen
shows the model of Gaussian beam propagation through an arbitrary thickness random phase screen, which the input and output plane are fixed and move the phase screen with thickness d = L 2 to modulate the phase of optical field. In the paper, all the three beam case, F 0 = ∞,F 0>0 and F 0<0,will be contained, and the beam waist appears in the middle of the optical path. Based on this figure, the nondimensional input plane beam parameters [4

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

7

7. L. C. Andrews and W. B. Miller, “Single- and double-pass propagation through complex paraxial optical systems,” J. Opt. Soc. Am. A 12(1), 137–150 (1995). [CrossRef]

] are

Θ0=1LF0Λ0=2LkW02.
(3)

For mathematical simplification, express W and F as

W=W0Θ02+Λ02F=F0(Θ02+Λ02)(Θ1)Θ02+Λ02Θ0.
(4)

As nondimensional input plane beam parameters, introduce the nondimensional parameters of the front plane and posterior plane of the phase screen and the output plan [4

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

7

7. L. C. Andrews and W. B. Miller, “Single- and double-pass propagation through complex paraxial optical systems,” J. Opt. Soc. Am. A 12(1), 137–150 (1995). [CrossRef]

]

z=L1Θ1=1+L1F1Λ1=2L1kW02,
(5)
z=L1+L2Θ2=1+L2F2Λ1=2L2kW02,
(6)
z=L1+L2+L3Θ3=1+L3F3Λ1=2L3kW02.
(7)

We can calculate the receiver plane propagation parameter p(L), according to Eqs. (5) to (7),

1p(L)=ΘiΛ=i=13(ΘiiΛi).
(8)

3. Rytov method

For the model of Gaussian beam propagation through a random phase screen, showed in Fig. 1, a single phase screen usually can represent weak fluctuation conditions, therefore, the statistical property will be discussed based on Rytov weak fluctuation method [4

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

,14

14. V. I. Tatarskii, Wave Propagation in a Turbulent Medium (New York: McGraw-Hill, 1961).

]. In the weak fluctuation condition, propagation path L, the optical field can be expressed as:

U(r,L)=U0(r,L)exp[Ψ(r,L)]=U0(r,L)exp[ψ1(r,L)+ψ2(r,L)+],
(9)

where U 0(r,L) is the free space Gaussian beam, described in Eqs. (2), Ψ(r,L) is a complex phase perturbation due to turbulence, and ψ 1(r,L) and ψ 2(r,L) are the first-order and second-order complex phase perturbations, respectively. It is known that the first two order complex phase perturbations can nearly determinate all the long time statistical property, so we just consider ψ 1(r,L) and ψ 2(r,L) in this paper. For the case of homogeneous and isotropic extended turbulence, it has been shown by Andrews and Miller [6

6. L. C. Andrews, W. B. Miller, and J. C. Ricklin, “Propagation through complex optical system: a phase screen analysis,” SPIE 2312, 122–129 (1994). [CrossRef]

,7

7. L. C. Andrews and W. B. Miller, “Single- and double-pass propagation through complex paraxial optical systems,” J. Opt. Soc. Am. A 12(1), 137–150 (1995). [CrossRef]

] that all statistical quantities of interest for an optical wave in the weak fluctuation regime are based on linear combinations of the following three second-order moments:

E1(0,0)=ψ2(r,L)+12ψ12(r,L)=2π2k20Ldz0κΦn(κ,z)dk,
(10)
E2(r1,r2)=ψ1(r1,L)ψ1*(r2,L)=4π2k20Ldz0dκκΦn(κ,z)exp[ΛLκ2(1ZL)2k]×J0{κ|[1Θ¯(1ZL)]p2iΛ(1ZL)r|},
(11)
E3(r1,r2)=ψ1(r1,L)ψ1(r2,L)=4π2k20Ldz0dκκΦn(κ,z)exp{iκ2γLk(1ZL)[1Θ¯(1ZL)]}×exp[ΛLκ2(1ZL)2k])J0{κρ[1(Θ¯+iΛ)(1ZL)]},
(12)

where r = (r 1 + r 2)/2, p = r 1-r 2, r = |r|, ρ = |p|, κ is the space wave number, Θ¯=1Θ, J 0(x) is the Bessel function of the first kind, the brackets < > denote an ensemble average, and the asterisk * denotes the complex conjugate field.

3.1 Arbitrary thickness phase screen

As shown in Fig. 1, the Gaussian beam propagates through an arbitrary thickness phase screen, with thickness d = L 2, thus we can assume that random media exist only in [L 1, L 1 + L 2] and change the integral variable:

1zL=ξ(0zL,L3Lξ(L2+L3)L).
(13)

In this paper, our analysis is based on Kolmogorov spectrum,

Φn(κ)=0.033Cn2κ11/3,
(14)

where the refractive-index structure parameter C n 2 is assumed constant, because the random media only exit over L 1zL 1 + L 2, so we can set Φn(κ) as:

Φn(κ)={0.033C^n2κ11/3L1zL1+L20others,
(15)

where C^n2 denotes the structure constant of the random phase screen, ^ is the corresponding phase screen statistical quantities. Based on these hypothesis and Eqs. (10), (11) and (12), we can obtain the second-order moments of the arbitrary thickness phase screen:

E^1(0,0)=ψ2(r,L)+12ψ12(r,L)=2π2k2LL3LL2+L3Ldξ0κΦn(κ)dk,
(16)
E^2(r1,r2)=ψ1(r1,L)ψ1*(r2,L)=4π2k2LL3LL2+L3Ldξ0dκκΦn(κ)exp[ΛLκ2ξ2k]×J0{κ|[1Θ¯ξ]p2iΛξr|},
(17)
E^3(r1,r2)=ψ1(r1,L)ψ1*(r2,L)=4π2k2LL3LL2+L3Ldξ0dκκΦn(κ)exp{iκ2γLkξ[1Θ¯ξ]}×exp[ΛLκ2ξ2k])J0{κρ[1(Θ¯+iΛ)ξ]},
(18)

The three second-order moments are very important for the statistical property of the arbitrary thickness model, for example the mutual coherent function and scintillation index are both the liner combination of the three moments.

3.2 Rytov Variance of Phase Screen

For comparison, introduce Rytov variation of phase screen as Rytov variance in atmosphere. The Rytov variance [15

15. D. L. Fried and J. B. Seidman, “Laser-beam scintillation in the atmosphere,” J. Opt. Soc. Am. 57(2), 181–185 (1967). [CrossRef]

] defines as, for extend random media on whole optical path,

σR2=8π2k20Ldz0dκκΦn(κ)[1cos(κ2(Lz)k)]=1.23Cn2k76L116.
(19)

Changing the integral variable and interval, we can get the Rytov variance for the arbitrary thickness model:

σ^R2=8π2k2LL3LL2+L3Ldξ0dκκΦn(κ)[1cos(κ2Lξk)].
(20)

It is clear that the integral in Eq. (20) does not converge for Kolmogorov spectrum, but we can utilize von Kármán spectrum to deduce approximate result, because note that von Kármán spectrum,

Φn(κ)=0.033Cn2exp(κ2/κm2)(κ2+κ02)11/6,
(21)

approximates to Kolmogorov spectrum, when κ 0→0+,κm→∞,and κ 0/κm→0+. Based on von Kármán spectrum and using the integral definition of confluent hypergeometric function and asymptotic formula [4

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

],

U(a;c;x)=1Γ(a)0exp(xt)ta1(1+t)ca1dt(a>0,Re(x)>0),
(22)
U(a;c;x)Γ(1c)Γ(1+ac)+Γ(c1)Γ(a)x1c(|x|1),
(23)

and set κ 0→0+,κm→∞,and κ 0/κm→0+, then we can calculate the integral (20):

σ^R24.7466Cn2k76L116Re{i56L3LL2+L3Lξ56dξ}=1.23Cn2k76L116[(L2+L3L)116(L3L)116],
(24)

where ≅ denotes the approximate method. By enforcing the equivalence [5

5. L. C. Andrews, R. L. Phillips, and A. R. Weeks, “Propagation of a Gaussian-beam wave through a random phase screen,” Waves Random Media 7(2), 229–244 (1997). [CrossRef]

] of Eqs. (19) and (24), one is led to the relation between structure constants

Cn2=C^n2[(L2+L3L)116(L3L)116].
(25)

Equations (25) shows that using a phase screen to simulate turbulence depends on the thickness of the phase screen and distance from the phase screen to receiver. It is not practical for doing this comparison, however, the phase screen always used indoors, which the simulation propagation path is shorter than atmosphere propagation path. For the model in Fig. 1, the comparison should be implemented as

Cn2LA116=C^n2L116[(L2+L3L)116(L3L)116],
(26)

where LA is the propagation path in atmosphere, and L is that in the phase screen model. The left side of Eq. (26) represents the turbulence effects, which are turbulence strength and propagation path, on optical field, and the right side of Eq. (26) represents simulation effects of the phase screen.

4. Mutual coherent function

The second order moment of optical filed on receiver plane, known as mutual coherent function (MCF), is defined as in a weak fluctuation condition,

Γ2(r1,r2,L)=U(r1,L)U*(r2,L)=U0(r1,L)U0*(r2,L)exp[2E1(0,0)+E2(r1,r2)]=Γ20(r1,r2,L)exp[2E1(0,0)+E2(r1,r2)],
(27)

where Γ2 0(r 1,r 2,L) is the MCF in free space

Γ20(r1,r2,L)=U0(r1,L)U0*(r2,L)=W02Wexp(2r2W2ρ22W2ikFpr),
(28)

Normally we don’t calculate 2E 1(0,0) + E 2(r 1 + r 2) directly, but make some simplification

Γ2(r1,r2,L)=Γ20(r1,r2,L)exp[2E1(0,0)+E2(r1,r2)]=Γ20(r1,r2,L)exp[σr2(r1,L)+σr2(r1,L)T12Δ(r1,r2,L)].
(29)

The definitions of variables in right side of Eq. (29) are discussed in [4

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

], and, for a phase screen with thickness d = L 2, they are defined as

σ^r2(r,L)=12[E2(r,r)E2(0,0)]=2π2k2LL3LL2+L3Ldξ0dκκΦn(κ)exp[ΛLκ2ξ2k][I0(2Λrξκ)1],
(30)
T^=2E1(0,0)E2(0,0)=4π2k2LL3LL2+L3Ldξ0dκκΦn(κ)[1exp(ΛLξ2κ2k)],
(31)
Δ^(r1,r2,L)=E2(r1,r1)+E2(r2,r2)2E2(r1,r2)=4π2k2LL3LL2+L3Ldξ0dκκΦn(κ)exp(ΛLκ2ξ2k)×{I0(2Λr1ξκ)+I0(2Λr2ξκ)2J0[|(1Θ¯ξ)p-2iΛξr|κ]},
(32)

where I 0(x) = J 0(ix) is the modified Bessel function of the first kind.

4.1 Mean irradiance

For identical observation points, r 1 = r 2 = r, the MCF determines the mean irradiance

I^(r,L)=Γ2(r,r,L)=W02Wexp(2r2W2)exp[2σ^r2(r,L)T^].
(33)

Using integral Eqs. (30) and (31) can calculate the mean irradiation. Based on the definition of Bessel function, I 0(x), and the property of Gamma function, Γ(x),

I0(x)=n=0(x/2)2nn!Γ(n+1)0esttx1dt=Γ(x)sx[Re(x)>0,Re(s)>0],
(34)

we can obtainσ^r2(r,L)

σ^r2(r,L)=0.815C^n2k76L116Λ56[(L3+L2L)83(L3L)83][1F11(56;1;2r2W2)],
(35)

where F11(a;c;x) is the confluent hypergeometric function.

For the integral in Eq. (31), it is not converge for Kolmogorov spectrum, but we can utilize von Kármán spectrum to approximate as the calculation of σ^R2 in Eq. (20),

T^1.63C^n2k76L116Λ56[(L3+L2L)83(L3L)83].
(36)

Thus the mean irradiance for the arbitrary thickness model is

I(r,L)=W02Wexp(2r2W2)exp[1.63Cn2k76L116Λ56[(L3+L2L)83(L3L)83]F11(56;1;2r2W2)].
(37)

For the thin phase screen model, the mean irradiance is given in [5

5. L. C. Andrews, R. L. Phillips, and A. R. Weeks, “Propagation of a Gaussian-beam wave through a random phase screen,” Waves Random Media 7(2), 229–244 (1997). [CrossRef]

]

I(r,L)=W02Wexp(2r2W2)exp[1.93σ^R2(Λd3)56F11(56;1;2r2W2)],
(38)

where

σ^R2=2.25C^n2k76L2L356
(39)

is the Rytov variance for thin phase screen, and d 3 = L 3/L.

In order to compare the two model, as shown in Fig. 1, the Gaussian beam propagates through a random phase screen with different thickness and propagation path L = 1m, and move the phase screen between transmitter and receiver to change the receive plane optical field. Figure 2
Fig. 2 Comparison between mean irradiance of thin phase screen and that of arbitrary thickness (T: thin phase screen; A: arbitrary thickness phase screen; d 3 = L 3/L)
gives the comparison of the mean irradiance for two models, as a function of ( 2/L)1/2. We can find out that there are few differences between Eqs. (37) and (38), even for thick phases (d/L = 0.1,d = 100mm), so the thin phase screen model for the mean irradiance can approximate to the arbitrary model well.

4.2 Mutual coherent function

Because the MCF of Gaussian beam propagation through an arbitrary thickness phase screen depends on the position of the two observation points, it is statistically inhomogeneous. However, it is a function of only the distant of the two points ρ, when they are central symmetry about optical axis, r 1 = -r 2,

Γ^2(ρ,L)=W02W2exp(14Λkρ2L)exp[T^12d^(ρ,L)],
(40)

where

d^(ρ,L)=8π2k2LL3LL2+L3Ldξ0dκκΦn(κ)exp(ΛLκ2ξ2k){1J0[(1Θ¯ξ)κρ]}.
(41)

It is hard to develop the analytic result for integral (41), but we can simplify it for numerical calculation. Use the series representation of Bessel function and confluent hypergeometric function

J0(x)=n=0(1)n(x/2)2nn!Γ(n+1)F11(a;c;x)=n=0(a)nxn(c)nn!
(42)

to deduce Eq. (41),

d^(ρ,L)=8.70C^n2k76L116Λ56H,
(43)

where

Η=L3LL2+L3Lξ53{F11[56;1;(1Θ¯ξ)2kρ24ΛLξ2]1}dξ.
(44)

According to these results, it is easy to get the normalized MCF

Γ^2(ρ,L)Γ^2(0,L)=exp[14Λkρ2L12d^(ρ,L)],
(45)

where Γ^2(0,L)is on axis value of the MCF. [5

5. L. C. Andrews, R. L. Phillips, and A. R. Weeks, “Propagation of a Gaussian-beam wave through a random phase screen,” Waves Random Media 7(2), 229–244 (1997). [CrossRef]

] gives the normalized MCF for thin phase screen model:

Γ^2(ρ,L)Γ^2(0,L)=exp[14Λkρ2L1.179σ^R2|1Θ¯d3|53(kρ2L)56],
(46)

Figure 3
Fig. 3 Comparison between normalized MCF of thin phase screen and that of arbitrary thickness (T: thin phase screen; A: arbitrary thickness phase screen; d 3 = L 3/L)
plots the comparison of MCF, which is shown as a function of ( 2/L)1/2, for arbitrary thickness model and thin phase screen model. From the comparison in Fig. 3, the difference of MCF for two models, Eqs. (45) and (46), is not apparent when the phase screen relative thickness d/L is less than 0.01 (d = 10mm), but it increase as the thickness increasing. The MCF for thin phase screen approximation is more accurate when the phase screen is near to the receiver plane than it is in other position. Therefore, it is better to use Eq. (45), MCF for arbitrary thickness model, to describe the normalized MCF when the relative thickness d/L is greater than 0.01. As a whole tendency, the MCF decreases as the separation distance ρ of two observation points increasing, which means that the further the two observation points separate, the worse the relativity of the two observation points is. And the weakening of the relativity of the two points is more severe when the phase screen is close to input plane.

4.3 Modulus of the complex degree of coherent

In addition to obtain the mean irradiance, the MCF can also be used to calculate the loss of spatial coherence of an initially coherent beam, which can be deduced form the modulus of the complex degree of coherence (DOC) [4

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

]

DOC(r1,r2,L)=|Γ2(r1,r2,L)|Γ2(r1,r1,L)Γ2(r2,r2,L)=exp[12D(r1,r2,L)],
(47)

where D(r 1,r 2,L) = Re[∆(r 1,r 2,L)] is the wave structure function (WSF). Written as a function of separation distance ρ, DOC(ρ,r,L) and D(ρ,r,L), the spatial coherence radius ρ 0 is defined by DOC(ρ 0,r,L) = 1/e or D(ρ 0,r,L) = 2.

For central symmetry observation points, r 1 = -r 2, the WSF for Gaussian beam propagation through an arbitrary thickness phase screen can be expressed as

D^(ρ,L)=8π2k2LL3LL2+L3Ldξ0dκκΦn(κ)exp(ΛLκ2ξ2k){I0(Λρξκ)J0[(1Θ¯ξ)κρ]}=d^(ρ,L)+4σ^r2(ρ2,L).
(48)

For thin phase screen approximation, the WSF [5

5. L. C. Andrews, R. L. Phillips, and A. R. Weeks, “Propagation of a Gaussian-beam wave through a random phase screen,” Waves Random Media 7(2), 229–244 (1997). [CrossRef]

] is

D^(ρ,L)=7.074σ^R2(Λd32)56{F11[56;1;(1Θ¯d3)2kρ24ΛLd32]F11(56;1;kΛρ24L)}.
(49)

Based on Eqs. (48) and (49), we can calculate the DOC by the two methods. Figure 4
Fig. 4 Comparison the DOC of thin and arbitrary thickness phase screen (T: thin phase screen; A: arbitrary thickness phase screen)
plots the

DOC, as a function of 2L 1/kW 0 2, for different thickness phase screen. The DOC, deduced by thin phase screen approximation, can only approximate the spatial coherent loss when the phase screen is located near the beam waist, and the further away from beam waist the phase screen is, the greater the differences between the two models, Eqs. (48) and (49), are. For the arbitrary thickness model, the DOC is near unit when the phase screen is close to beam waist, and it gets down as the phase screen is far away from beam waist. So the DOC [5

5. L. C. Andrews, R. L. Phillips, and A. R. Weeks, “Propagation of a Gaussian-beam wave through a random phase screen,” Waves Random Media 7(2), 229–244 (1997). [CrossRef]

] for thin phase screen is suitable for describing the beam waist coherence, but for describing the coherence of other position, we should use the DOC for arbitrary thickness model.

5. Scintillation index

Scintillation index σI 2(r,L), one of the most important components of fourth order field moment, describes the irradiance fluctuation on receiver plane. When the log-amplitude variance is sufficiently small σχ 2 = 1, which is usually valid, the scintillation index [4

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

] is defined as

σI2(r,L)=exp[4σχ2(r,L)]14σχ2(r,L)=2Re[E2(r,r)+E3(r,r)].
(50)

For mathematic simplification, it is convenient to express the scintillation index as a sum of radial and longitudinal components, σI,r 2(r,L) and σI,l 2(L), in the form [4

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

]

σI2(r,L)=σI,r2(r,L)+σI,l2(L)=4σr2(r,L)+σI,l2(L).
(51)

For phase screen case, σ^r2(r,L) has been given by Eq. (35). σ^I,l2(L), the longitudinal component, is defined as

σ^I,l2(L)=8π2k2LL3LL2+L3Ldξ0dκκΦn(κ)exp[ΛLξ2kκ2]{1cos[Lκ2kξ(1Θ¯ξ)]}.
(52)

Based on the property of the Gamma function Γ(x) [15

15. D. L. Fried and J. B. Seidman, “Laser-beam scintillation in the atmosphere,” J. Opt. Soc. Am. 57(2), 181–185 (1967). [CrossRef]

]

0xλ1(eμxeνx)dx=Γ(x)(μλνλ)Re(λ)>2,Re(μ)>0,Re(ν)>0,
(53)

and the integral equation [4

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

]

0xtμ1(1+βt)νdt=xμμF21(ν,μ;1+μ;βx)Re(μ)>0,
(54)

we can deduce the longitudinal component σ^I,l2(L)

σ^I,l2(L)=4.74C^n2k76L116Re{i56(L3+L2L)11/6F21[56,116;176;(Θ¯+iΛ)L3+L2L]i56(L3L)11/6F21[56,116;176;(Θ¯+iΛ)L3L]1116Λ56[(L3+L2L)8/3(L3L)8/3]}.
(55)

So the scintillation index for Gaussian beam propagation through an arbitrary thickness phase screen is the sum of Eqs. (35) and (55)

σ^I2(r,L)=4σ^r2(r,L)+σ^I,l2(L)=3.26C^n2k76L116Λ56[(L3+L2L)83(L3L)83][1F11(56;1;2r2W2)]+4.74C^n2k76L116Re{i56(L3+L2L)11/6F21[56,116;176;(Θ¯+iΛ)L3+L2L]i56(L3L)11/6F21[56,116;176;(Θ¯+iΛ)L3L]1116Λ56[(L3+L2L)8/3(L3L)8/3]}.(rW)
(56)

Taking the thin phase screen approximation, the scintillation index is given in [5

5. L. C. Andrews, R. L. Phillips, and A. R. Weeks, “Propagation of a Gaussian-beam wave through a random phase screen,” Waves Random Media 7(2), 229–244 (1997). [CrossRef]

]

σ^I2(r,L)=6.45σ^R2(Λd3)56(r2W2)+3.87σ^R2{[(Λd3)2+(1Θ¯d3)2]512cos[56arctan(1Θ¯d3Λd3)](Λd3)56}.(rW)
(57)

Figure 4 plots the scintillation index, for Gaussian beam propagation through different thickness phase screens, which relativity thickness d/L are 10−4 and 0.01, as a function of the phase screen position 2L 1/kW 0 2. From Fig. 4, it is clear that there are apparent differences between the two models, Eqs. (56) and (57), even for the thin phase screen, d/L = 10−4. For the relationship described by Eq. (56), the scintillation index is get down firstly and then increasing, as the phase screen moving from input plane to output plane, and it obtain the minimum value when the phase screen is located in beam waist. The scintillation index is in direct proportion to the thickness of phase screens when they are located in the same position. The results predicted by scintillation index model, Eq. (56), are according with those shown by DOC model, Eq. (47), because the weakening of coherence results in the fluctuation of irradiance as shown in Fig. 4 and Fig. 5
Fig. 5 Comparison of scintillation index with different thickness (T: thin phase screen; A: arbitrary thickness phase screen)
.

6. Summary

In this paper, we, based on Rytov method and Andrews et al’s [4

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

7

7. L. C. Andrews and W. B. Miller, “Single- and double-pass propagation through complex paraxial optical systems,” J. Opt. Soc. Am. A 12(1), 137–150 (1995). [CrossRef]

] works, establish the mathematic models of statistical quantities for Gaussian beam propagation through an arbitrary thickness phase screen and developed the analytical results. The comparison between Andrews’ thin phase screen approximation results and our results shows that Andrews’ results are concision but has its own limitation, and our results are suitable for more generalized case, though complicate, especially for the scintillation index. For mean irradiance on receiver plane and the MCF of thin phase screen (d/L < 0.01), the two models have few differences, so both of them can be used to describe these quantities accurately. But for MCF of thick phase screen, DOC, and scintillation index, we’d better to utilize the arbitrary thickness phase screen model to describe them, because the thin phase screen model has intrinsic error due to the approximation. The paper offers a theory fundamental, for AO and laser communication application research in laboratory.

Acknowledgments

The authors are very thankful to the reviewers for valuable comments. This work was supported by the Innovation Foundation of the Chinese Academy of Sciences III under grant O98Y32C100.

References and links

1.

C. Chao, L. Hu, and Q. Mu, “Bandwidth requires of adaptive optical system for horizontal turbulence correction,” Opt. Precision Eng. 18(10), 2137–3142 (2010).

2.

B. Wang, Z.-y. Wang, J.-l. Wang, J.-Y. Zhao, Y.-H. Wu, S.-X Zhang, L. Dong, and M. Wen, “Phase-diverse speckle imaging with two cameras,” Opt. Precision Eng. 19(6), 1384–1390 (2011). [CrossRef]

3.

H. G. Booker, J. A. Ferguson, and H. O. Vats, “Comparison between the extended-medium and the phase-screen scintillation theories,” J. Atmos. Terr. Phys. 47(38), 1–399 (1985). [CrossRef]

4.

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

5.

L. C. Andrews, R. L. Phillips, and A. R. Weeks, “Propagation of a Gaussian-beam wave through a random phase screen,” Waves Random Media 7(2), 229–244 (1997). [CrossRef]

6.

L. C. Andrews, W. B. Miller, and J. C. Ricklin, “Propagation through complex optical system: a phase screen analysis,” SPIE 2312, 122–129 (1994). [CrossRef]

7.

L. C. Andrews and W. B. Miller, “Single- and double-pass propagation through complex paraxial optical systems,” J. Opt. Soc. Am. A 12(1), 137–150 (1995). [CrossRef]

8.

S. V. Mantravadi, T. A. Rhoadarmer, and R. S. Glas, “Simple laboratory system for generating well-controlled atmospheric-like turbulence,” Proc. SPIE 5553, 290–300 (2004). [CrossRef]

9.

X. J. Gan, J. Guo, and Y. Y. Fu, “The simulating turbulence method of laser propagation in the inner field,” J. Phys. Conf. Ser. 48, 907–910 (2006). [CrossRef]

10.

B. D. Zhang, S. Qin, and X. S. Wang, “Accurate and fast simulation of Kolmogorov phase screen by combining spectral method with Zernike polynomials method,” Chin. Opt. Lett. 8(10), 969–971 (2010).

11.

M. X. Qian, W. Y. Zhu, and R. Rao, “Phase screen distribution for simulating laser propagation along an inhomogeneous atmospheric path,” Acta Phys. Sinica 58(9), 6633–6639 (2009).

12.

M. Gao and Z. Wu, “Experiments of effect of beam spreading of far-field on aiming deviation,” Opt. Precision Eng. 18(3), 602–608 (2010).

13.

Y.-h. Gao, Z.-y. An, N.-n. Li, W.-x Zhao, and J.-s. Wang, “Optical design of Gaussian beam shaping,” Opt. Precision Eng. 19(7), 1464–1471 (2011). [CrossRef]

14.

V. I. Tatarskii, Wave Propagation in a Turbulent Medium (New York: McGraw-Hill, 1961).

15.

D. L. Fried and J. B. Seidman, “Laser-beam scintillation in the atmosphere,” J. Opt. Soc. Am. 57(2), 181–185 (1967). [CrossRef]

OCIS Codes
(010.1080) Atmospheric and oceanic optics : Active or adaptive optics
(010.1290) Atmospheric and oceanic optics : Atmospheric optics
(010.1300) Atmospheric and oceanic optics : Atmospheric propagation
(010.1330) Atmospheric and oceanic optics : Atmospheric turbulence
(060.4510) Fiber optics and optical communications : Optical communications

ToC Category:
Atmospheric and Oceanic Optics

History
Original Manuscript: July 7, 2011
Revised Manuscript: August 15, 2011
Manuscript Accepted: August 22, 2011
Published: September 1, 2011

Citation
Yuzhen Tian, Jin Guo, Rui Wang, and Tingfeng Wang, "Mathematic model analysis of Gaussian beam propagation through an arbitrary thickness random phase screen," Opt. Express 19, 18216-18228 (2011)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-19-19-18216


Sort:  Author  |  Year  |  Journal  |  Reset  

References

  1. C. Chao, L. Hu, and Q. Mu, “Bandwidth requires of adaptive optical system for horizontal turbulence correction,” Opt. Precision Eng.18(10), 2137–3142 (2010).
  2. B. Wang, Z.-y. Wang, J.-l. Wang, J.-Y. Zhao, Y.-H. Wu, S.-X Zhang, L. Dong, and M. Wen, “Phase-diverse speckle imaging with two cameras,” Opt. Precision Eng.19(6), 1384–1390 (2011). [CrossRef]
  3. H. G. Booker, J. A. Ferguson, and H. O. Vats, “Comparison between the extended-medium and the phase-screen scintillation theories,” J. Atmos. Terr. Phys.47(38), 1–399 (1985). [CrossRef]
  4. L. C. Andrews and R. L. Phillips, Laser Beam Propagation through Random Media, 2nd ed. (SPIE, 2005).
  5. L. C. Andrews, R. L. Phillips, and A. R. Weeks, “Propagation of a Gaussian-beam wave through a random phase screen,” Waves Random Media7(2), 229–244 (1997). [CrossRef]
  6. L. C. Andrews, W. B. Miller, and J. C. Ricklin, “Propagation through complex optical system: a phase screen analysis,” SPIE2312, 122–129 (1994). [CrossRef]
  7. L. C. Andrews and W. B. Miller, “Single- and double-pass propagation through complex paraxial optical systems,” J. Opt. Soc. Am. A12(1), 137–150 (1995). [CrossRef]
  8. S. V. Mantravadi, T. A. Rhoadarmer, and R. S. Glas, “Simple laboratory system for generating well-controlled atmospheric-like turbulence,” Proc. SPIE5553, 290–300 (2004). [CrossRef]
  9. X. J. Gan, J. Guo, and Y. Y. Fu, “The simulating turbulence method of laser propagation in the inner field,” J. Phys. Conf. Ser.48, 907–910 (2006). [CrossRef]
  10. B. D. Zhang, S. Qin, and X. S. Wang, “Accurate and fast simulation of Kolmogorov phase screen by combining spectral method with Zernike polynomials method,” Chin. Opt. Lett.8(10), 969–971 (2010).
  11. M. X. Qian, W. Y. Zhu, and R. Rao, “Phase screen distribution for simulating laser propagation along an inhomogeneous atmospheric path,” Acta Phys. Sinica58(9), 6633–6639 (2009).
  12. M. Gao and Z. Wu, “Experiments of effect of beam spreading of far-field on aiming deviation,” Opt. Precision Eng.18(3), 602–608 (2010).
  13. Y.-h. Gao, Z.-y. An, N.-n. Li, W.-x Zhao, and J.-s. Wang, “Optical design of Gaussian beam shaping,” Opt. Precision Eng.19(7), 1464–1471 (2011). [CrossRef]
  14. V. I. Tatarskii, Wave Propagation in a Turbulent Medium (New York: McGraw-Hill, 1961).
  15. D. L. Fried and J. B. Seidman, “Laser-beam scintillation in the atmosphere,” J. Opt. Soc. Am.57(2), 181–185 (1967). [CrossRef]

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.

Figures

Fig. 1 Fig. 2 Fig. 3
 
Fig. 4 Fig. 5
 

« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited