A stochastic beam generated by a recently introduced source of Schell type with cosine-Gaussian spectral degree of coherence is shown to possess interesting novel features on propagation in isotropic and homogeneous atmospheric turbulence with general non-Kolmogorov power spectrum. It is shown that while at small distances from the source the beam’s intensity exhibits annular profile with adjustable area of the dark region, the center disappears at sufficiently large distances and the beam’s intensity tends to Gaussian form. Hence the 3D bottle beam is produced by the cumulative effect of the random source and the atmosphere. The distances at which the on-axis beam intensity has local minima and maxima are shown to have analytic dependence on the source and the atmospheric parameters. And the influence of the fractal constant of the atmospheric power spectrum and refractive-index structure constant on beam characteristics is analyzed in depth. The novel double-cycle qualitative change in the degree of coherence is shown to occur on atmospheric propagation which was not previously known for any other beams.
© 2013 Optical Society of America
Several models for scalar random sources have been recently introduced in addition to well-explored Gaussian-Schell model source  differing by the form of their spectral degree of coherence. Among such classes are the Bessel-correlated sources [2, 3], the non-uniformly correlated sources  and the Multi-Gaussian Schell-model sources [5, 6]. Apart from these classes a family of Schell-model sources stands out in which the classic Gaussian degree of coherence is modulated by a well-behaved function, leading to Bessel-Gaussian-correlated, Laguerre-Gaussian-correlated  and Cosine-Gaussian-correlated Schell-model sources . The sources in  and  all generate random beams which can have arbitrary intensity distribution in the source plane (say, Gaussian) but acquire propagation-invariant ring-shaped intensity profile in the far-field (on propagation in free space). Such behavior distinguishes these random beams from the well-known deterministic and random dark-hollow beams (DHB) which possess ring-shaped intensity profiles already at the source plane but loose them on free-space propagation (see  and references wherein). It is shown that the DHB fills in the center and takes on the Gaussian form at shorter distances if the atmosphere is present . One should not confuse the Cos-Gaussian-correlated Schell-model beams with the beams of Ref. [11, 12] in which the degree of coherence is just Gaussian but the intensity is Gaussian function modulated by cosine function.
The purpose of this paper is to tackle the clear-air atmospheric propagation of Cos-Gaussian-correlated Schell-model beams which we will briefly refer to as CGC beams and to show whether and how the optical turbulence modifies the rings. Compared to other sources producing the ring-shaped beams  the CGC beam  is the most attractive because its modulation by cosine function leads to analytic results for beams propagating in free-space and in the atmosphere with quite arbitrary power spectra. Further, the ability of the CGC sources will be demonstrated to produce beams which build up and sustain robust intensity rings at intermediate distances in the turbulent medium. Moreover, we will show that the filling-in-the-center effect occurring at large distances imply that a 3D bottle-like beam  can be obtained which circumvents objects on propagation in random medium. This feature makes the CGC beams attractive in situations when an optical signal must not be intercepted by a detector located on the optical axis of the beam. The other goal of this work is to explore the evolution in coherence properties of the CGC beam on its atmospheric propagation. In order to cover a wide scope of atmospheric conditions we employ the atmospheric spectrum in which in addition to the structure constant, the inner and the outer scales  the slope of the spectrum in the inertial range can be tuned .
2. Analytic solutions for the CGSM propagating in turbulent atmosphere
The cross-spectral density function of a wide-sense statistically stationary field generated by the CGC source located in the plane z = 0, at points specified by two-dimensional position vectors ρ′1 = (x′1, y′1) and ρ′2 = (x′2, y′2) and oscillating at angular frequency ω has form 
The paraxial form of the extended Huygens-Fresnel principle which describes the interaction of waves with isotropic and homogeneous random medium implies that the elements of the cross-spectral density matrix at two positions r1 = (ρ1, z) and r2 = (ρ2, z) in the same transverse plane of the half-space z > 0 are related to those in the source plane as Eq. (2) can be given by the expression  Eq. (5) is represented by model in [15, 16], in which the slope 11/3 of the conventional van Karman spectrum is generalized to an arbitrary parameter α, i.e. Eq. (7) is a generalized refractive-index structure parameter with units m3−α. With the power spectrum in Eq. (6) the integral in Eq. (5) becomes 
3. Numerical examples
Now we will analyze the behavior of the second-order statistical properties of the CGC beams on propagation in free space and in turbulent atmosphere by a set of numerical examples based on the analytical formulas derived in the previous section. Without loss of generality the parameters of the source and of atmosphere are chosen to be σ = 1 cm, δ = 1 mm, λ = 632.8 nm, , L0 = 1 m, l0 = 1 mm, α = 3.67 unless other values are specified in figure captions.
In Fig. 1 the evolution of the spectral density S calculated from Eq. (9) of several CGC beams with different values of index n, along the optical axis (ρ = 0) is shown. Free-space Fig. 1(a) and atmospheric Fig. 1(b) clearly show the qualitative difference in the beam behavior. While on propagation in free space the on-axis intensity decreases monotonically for all values of n, on propagation in turbulence at the intermediate distances from the source the dark region is formed, but it is gradually filled in at sufficiently large distances, on the order of tens of kilometers. The higher the value of n the more efficiently the dark center is formed: while for n = 0 (classic Gaussian Schell-model beam) and n = 1 the on-axis intensity drop does not occur, for n = 4 it becomes significant. The distances at which local minima and maxima are reached for each order n can be found analytically. Indeed on finding the derivative of S(ρ = 0, z) with respect to z, setting it to zero and making some simplifications one obtains the fourth-order equation of the formEq. (11) has two positive real roots, the smaller/larger root corresponding to the local intensity minimum/maximum, respectively. For example, for the parameters used for the curves in Fig. 1 one finds that the local minima/maxima of on-axis spectral density S occur at distances zmin(n = 2) = 3891.37 m and zmax(n = 2) = 28210.98 m; zmin(n = 3) = 1304.26 m and zmax(n = 3) = 97622.85 m; zmin(n = 4) = 912.36 m and zmax(n = 4) = 191510.70 m. There are no maxima and minima for orders n < 2.
In Fig. 2 the evolution of the whole transverse cross-section of the CGC beams with different values of parameter n is shown at several selected distances z from the source plane. Unlike for mode n = 0, which remains Gaussian for all distances, the beams with n > 0 first form the dark centers, but at larger distances regain the Gaussian shape. The higher the order n the wider and the sooner the ring-shaped profile is formed and remains more robust to optical turbulence. As is seen from Fig. 2(d) plotted for distance z = 1000 km the CGC beams’ intensity still remain a ring-like for n > 1, while the beam with n = 1 reconstructs its original Gaussian profile. We note that the horizontal and the vertical scales of parts (a), (b), (c) and (d) are different addressing the substantial expansion of the beams on propagation.
Figure 3 illustrates the behavior of the transverse cross-sections of the beam with n = 2 as a function of α (a) and (b) and implies that the influence of both parameters on the intensity distribution is crucial. In Fig. 4 the contour plots of the spectral density as a function of z and ρ are shown in two regions: close to the source (a) and sufficiently far from the source (b). The formation and disappearance of the ring structure is evident.
We now turn our attention to the changes in the spectral degree of coherence calculated with the help of Eq. (10). Figure 5 illustrates the evolution of the transverse degree of coherence evaluated at separation distance ρd between two points ρ1 = ρ/2 and ρ2 = −ρ/2 located symmetrically with respect to the optical axis: ρd = |ρ1 − ρ2|. The degree of coherence is shown for the same five CGC modes as the spectral density in Fig. 1, at several selected distances. In the source plane [Fig. 5(a)] the degree of coherence of the beam differs with parameter n, having Gaussian form for n = 0 and starting to resemble a sinc-function as n increases. For relatively short ranges from the source [Fig. 5(b)] the degree of coherence changes its form to Gaussian for all n. At a distance of approximately 600 m (for the given source and atmospheric parameters) they converge to a single Gaussian distribution [Fig. 5(c)]. This modification can be attributed to the source correlations. Further, the degree of coherence is being qualitatively modified for the second time by the turbulence: as is seen from Figs. 5(d)– 5(f) it goes through another change depending on the mode n, constituting a novel effect in atmospheric propagation. At distances of tens of kilometers the sinc-like profile is formed [Fig. 5(e)], which later converts to Gaussian profile [Fig. 5(f)]. In contrast, it can be readily seen from Fig. 5 that for the GSM beams (n = 0, brown solid thick curve) propagating in the atmosphere the degree of coherence goes through a single cycle: while the Gaussian profile is always preserved its transverse r.m.s. correlation width first grows due to source correlations and then reduces due to turbulence.
4. Concluding remarks
In this article we have investigated the effect of atmospheric turbulence on propagation characteristics of a recently introduced class of Schell-model beams with Cos-Gaussian degree of coherence (CGC beams). The analytical formula for the cross-spectral density function of such a beam for propagation in the atmosphere with general, non-classic power spectrum of the refractive index fluctuations is derived and used to explore the evolution of the spectral density and the spectral degree of coherence in various atmospheric conditions.
We have found that the CGSM source can produce an annular beam in weak turbulence, robust at least for several kilometers in typical Kolmogorov atmosphere. As the propagation distance grows and/or atmospheric fluctuations become more severe (α reduces and/or grows) the hollow beam profile gradually disappears being again converted into a wide Gaussian distribution. Thus on propagation in optical turbulence the CGC beam has the 3D bottle-beam-like structure, convenient for circumventing objects at intermediate distances. The dimensions and the shape of such a 3D optical bottle can be adjusted by the choice of source parameters σ, δ and n.
Another discovered effect that is not pertinent to other deterministic and random beams is the double-cycle qualitative change of the transverse spectral degree of coherence. Namely, for all modes with n > 0, after changing the profile from cos-Gaussian to Gaussian at some intermediate distance, the degree of coherence suffers another round of transformation: it losses the Gaussian profile at distances on the order of several kilometers (to a function resembling a sinc-Gaussian function with low orders) and regains it only at very large distances. The appearance of the double cycle in the evolution of the degree of coherence can be explained by two mechanisms, source correlations and atmospheric turbulence, subsequently affecting the beam. At small propagation distances the beam is mostly modified by the source correlations, leading to the ring-shaped intensity and the Gaussian degree of coherence. In the absence of atmospheric fluctuations this would be the single cycle everywhere in the far field. However, as the beam passes to larger distances through turbulence, the latter modifies the beam intensity distribution from ring back to Gaussian while narrowing the degree of coherence. The combination of these two atmospheric effects produces the second cycle. Note that for classic Gaussian Schell-model beams, when neither intensity nor degree of coherence are modified in shape there is still double modification in the width of the transverse degree of coherence: it first widens due to source correlations and then gradually becomes narrower on propagation at sufficiently long distances in the atmosphere . Thus, our new effect can be regarded as the generalization of the one discussed in , relating not to scale change but to shape transformation.
Thus, the analytic and numerical results that we obtained relating to the atmospheric propagation of recently introduced CGC beams  revealed that the turbulence can suppress transformation in the intensity profile (e.g. Gaussian to annular), initially formed by source correlations, converting it back (from annular to Gaussian) at large distances. Such complex intensity transformation was found to be accompanied by the double cycle change in the shape of the degree of coherence. Our finding about the reconstruction mechanism of the random beam intensity on propagation in atmospheric turbulence supports the similar results reported for other statistics, such as spectral composition  and polarization .
Z. Mei’s research is supported by the National Natural Science Foundation of China (NSFC) ( 11247004) and Zhejiang Provincial Natural Science Foundation of China ( Y6100605). E. Shchepakina is supported in part by the Russian Foundation for Basic Research (grants 12-08-00069a, 13-01-97002-p) and AFOSR ( FA9550-12-1-0449). O. Korotkova’s research is supported by US AFOSR (FA9550-12-1-0449).
References and links
1. L. Mandel and E. Wolf, Optical Coherence and Quantum Optics(Cambridge University, Cambridge, 1995) [CrossRef] .
2. F. Gori, G. Guattari, and C. Padovani, “Modal expansion for J0-correlated Schell-model sources,” Opt. Commun. 64, 311–316 (1987) [CrossRef] .
6. O. Korotkova, S. Sahin, and E. Shchepakina, “Multi-Gaussian Schell-model beams,” J. Opt. Soc. Am. A 29, 2159–2164 (2012) [CrossRef] .
8. Z. Mei and O. Korotkova, “Cosine-Gaussian Schell-model sources,” (Opt. Lett ., in press).
9. Y. Cai and F. Wang, “Partially coherent anomalous hollow beam and its paraxial propagation,” Phys. Lett. A 372, 4654–4660 (2008) [CrossRef] .
10. M. Alavinejad, G. Taherabadi, N. Hadilou, and B. Ghafary, “Changes in the coherence properties of partially coherent dark hollow beam propagating through atmospheric turbulence,” Opt. Commun. 288, 1–6 (2013) [CrossRef] .
11. H. T. Eyyuboğlu and Y. Baykal, “Cosine-Gaussian laser beam intensity in turbulent atmosphere,”Proc. SPIE 5743, 131–141 (2004) [CrossRef] .
12. H. T. Eyyuboğlu and Y. Baykal, “Transmittance of partially coherent cosh-Gaussian, cos-Gaussian and annular beams in turbulence,” Opt. Commun. 278(1), 17–22 (2007) [CrossRef]
13. J. Pu, X. Liu, and S. Nemoto, “Partially coherent bottle beams,” Opt. Commun. 252, 7–11 (2005) [CrossRef] .
14. L.C. Andrews and R. L. Phillips, Laser Beam Propagation in Random Media(SPIE Press, Bellington, 2005) [CrossRef] .
15. I. Toselli, L. C. Andrews, R. L. Phillips, and V. Ferrero, “Free-space optical system performance for laser beam propagation through non-Kolmogorov turbulence,” Opt. Eng. 47, 026003 (2008) [CrossRef] .
17. W. Lu, L. Liu, J. Sun, Q. Yang, and Y. Zhu, “Change in degree of coherence of partially coherent electromagnetic beams propagating through atmospheric turbulence,” Opt. Commun. 271(1), 1–8 (2007) [CrossRef] .
19. E. Shchepakina and O. Korotkova, “Second-order statistics of stochastic electromagnetic beams propagating through non-Kolmogorov turbulence,” Opt. Express 18, 10650–10658 (2010) [CrossRef] [PubMed] .