Abstract
Expressions for the correlation coefficient between light-flux fluctuations of two waves counter-propagating along a common path in weak turbulence are developed. Only the aperture and inner-scale Fresnel parameters are needed for evaluation of the correlation coefficient if the turbulence spectrum has no path dependence, and of the path weighting functions for the cross-covariance and variances of normalized light-flux fluctuations if the turbulence spectrum is dependent on path locations. Under the condition that atmospheric turbulence is statistically homogeneous over a path, although good correlation between light-flux fluctuations of two counter-propagating spherical waves may be achieved for a relatively small aperture Fresnel parameter or relatively large inner-scale Fresnel parameter, the correlation coefficient between light-flux fluctuations of two counter-propagating plane waves is always lower than 1 obviously. When the aperture Fresnel parameter becomes larger than the inner-scale Fresnel parameter, the inner scale of turbulence tends to play an unimportant role in determining the correlation coefficient.
© 2017 Optical Society of America
1. Introduction
Over the years, the scientific community has paid much attention to optical wave propagation in atmospheric turbulence. Considerable research effort has been devoted to the cases of single-direction propagation. However, the situations where two optical waves counter-propagate along a common path can also be encountered in many practical applications [1–6]; among them are bidirectional free-space optical (FSO) links, in which two waves emanating from two separated transmitters propagate in opposite directions. For a bidirectional FSO link, the correlation coefficient between turbulence-induced light-flux fluctuations of the two counter-propagating waves is actually a statistical quantity of practical interest. On the one hand, the correlation between received light-flux fluctuations at two link ends may be utilized to obtain the instantaneous channel state directly at the transmitter, hence eliminating the need for a dedicated channel-state feedback link; the instantaneous channel-state information is useful for adaptive transmission in FSO communication systems impaired by atmospheric turbulence [5,7]. On the other hand, the correlation between received light-flux fluctuations at two link ends reveals that a turbulent FSO channel has the potential for generating common randomness, which can be used to extract random secret keys shared by the two link ends [8]; shared random secret keys are valuable in cryptographic applications for secure communications. To better exploit the correlation between received light-flux fluctuations at two link ends, an understanding of its behavior is imperative.
For several decades, the point-source point-receiver (PSPR) reciprocity has proven to remain valid even when atmospheric turbulence exists over the propagation path [9–11], meaning that the received light signals at the two ends of a bidirectional FSO link with two point transceivers are perfectly correlated, irrespective of atmospheric turbulence. Nevertheless, experimental measurement has shown that the correlation coefficient between received light-flux fluctuations at two link ends may obviously decrease below 1 when the two transceivers have a finite-size aperture [12]. An important exception reported in [1,4,5] is that the registered light signals of two transceivers consisting of single-mode fiber collimators (SMFC) in a bidirectional FSO link can always preserve prefect correlations, regardless of atmospheric turbulence; the experimental observations presented in [1,13] provide evidence to support this theoretical prediction. Coupling an optical wave into a single-mode fiber in an FSO receiver may result in a significant power loss [14–18], implying that the light signal received by an SMFC transceiver is not equivalent to the light flux entering its aperture; to overcome this problem, power-in-the-bucket (PIB) receivers [1], which can measure all the light flux collected by a receiving aperture, are commonly used in FSO communication systems. Accordingly, it is interesting to explore the behavior of the statistical correlation between light-flux fluctuations received at the two ends of a bidirectional FSO link through atmospheric turbulence.
By taking into account finite-size receiving apertures, Perlot and Giggenbach [3] addressed the correlation between light-flux fluctuations of two counter-propagating spherical waves in weak atmospheric turbulence; they formulated the correlation coefficient based on the spatial covariance function of irradiance fluctuations and obtained a relevant expression whose numerator includes a five-dimensional integral (see Eq. (22) in [3] where bF (·) and HFR (·) are essentially expressed in the form of a one- and two-dimensional integral, respectively). To·derive simpler and more insightful mathematical models, unlike the traditional treatment of light-flux fluctuations, we will first develop the relationship between the normalized light flux and aperture-averaged first-order log-amplitude fluctuations, and then formulate the correlation coefficient between light-flux fluctuations of two counter-propagating waves in terms of the statistics of aperture-averaged first-order log-amplitude fluctuations. By doing so, the spatial covariance function of irradiance fluctuations is not necessary for our theoretical development, hence simplifying the derivation. Furthermore, although a zero inner scale of turbulence has actually been implicitly assumed in [3], the inner scale may have an impact on optical scintillations. In this work, we will incorporate the inner scale of turbulence into our theoretical derivation and elucidate its role on the correlation coefficient between light-flux fluctuations of two counter-propagating waves. Compared with the existing effort, we intend to deal with the correlation between light-flux fluctuations of two counter-propagating waves in weak turbulence from different points of view. Specifically, we attempt to understand the physical aspects of this issue in both wavenumber and path-space domains.
2. Theoretical analysis
Figure 1 illustrates the counter-propagation geometry under consideration here for two waves with the same wavelength in atmospheric turbulence. The propagation from the plane at z = 0 to the plane at z = L is referred to as “forward propagation” and that from the plane at z = L to the plane at z = 0 is called “inverse propagation”. Plane, spherical and Gaussian-beam waves are the common model waves considered in the existing literature [19–21] related to optical propagation through atmospheric turbulence. Among them, the Gaussian-beam model may be a tool most appropriate for describing the propagated waves from the perspective of many practical applications. However, for the sake of tractability, here we restrict our analysis to the first two kinds of model waves, i.e., the plane and spherical waves. Charnotskii [22] performed a comprehensive asymptotic analysis of beam-wave scintillations under both weak-and strong-turbulence conditions and presented a complete set of asymptotes. As shown in [23], under the near- and far-field conditions, i.e., the initial-beam-size Fresnel parameter is much larger and smaller, respectively, than 1, scintillation of collimated-Gaussian-beam waves in weak-turbulence cases, viz., when the coherence-radius Fresnel parameter is much greater than 1, approaches that of plane and spherical waves, respectively; additionally, the light-flux fluctuation is closely related to the scintillation. With these in mind, our later treatment of the correlation between light-flux fluctuations of two counter-propagating plane and spherical waves in weak turbulence would be expected to provide some somewhat reasonable estimates for the cases of two identical counter-propagating collimated Gaussian beams in the near- and far-field regions, respectively.
2.1. Two counter-propagating plane waves
By making use of the Rytov approximation, the first-order complex phase perturbation of a plane wave propagating in atmospheric turbulence from z = 0 to z = L can be expressed by [19,20]
with where r denotes a point in the plane at z = L, k is the optical wavenumber of the propagating wave, κ = |κ| is the magnitude of κ, dv (κ, z) is related to the refractive-index fluctuation n1 (r′, z) by Performing the substitution z → (L − z) in Eq. (2) leads to the expression for the first-order complex phase perturbation of a plane wave propagating from z = L to z = 0 with where ρ is a point in the plane at z = 0. Notice that, the subscripts “F” and “I” in the above expressions imply “forward propagation” and “inverse propagation”, respectively.Note that, ψ1,pl,F (r, z = L) = χ1,pl,F (r, z = L) + iS1,pl,F (r, z = L), where χ1,pl,F (·) and S1,pl,F (·) are the first-order log-amplitude and phase fluctuations of the forward-propagating plane wave, respectively. By recalling the line of reasoning presented in Sec. 17–6 of [20], after some manipulations, the first-order log-amplitude fluctuation for the forward-propagating wave can be written by
where the asterisk stands for the complex conjugate. In a similar manner, one can formulate the first-order log-amplitude fluctuation for the inverse-propagating wave byIn what follows, we assume that two circular apertures with the same size, centered on the z-axis and located at z = 0 and z = L, respectively, are used to collect the counter-propagated waves; the aperture planes are perpendicular to the z-axis. It is well known that a finite-size receiving aperture may cause a decrease in optical scintillations, which is called aperture averaging. Conceptually, to determine the statistics of light-flux fluctuations received by an aperture, one needs to integrate the irradiance instead of the log-amplitude fluctuation within the aperture. Thus, it would seem that the statistics of light-flux fluctuations cannot be directly derived from the log-amplitude fluctuation within the aperture. To proceed further, below we first deduce a relationship between the normalized light flux and aperture-averaged first-order log-amplitude fluctuation under weak-turbulence conditions. For an incident plane wave distorted by atmospheric turbulence, the normalized light flux received by a circular aperture can be expressed as
where I (r) denotes the turbulence-induced irradiance at position r, is the plane-wave irradiance at the aperture plane in the absence of atmospheric turbulence, a0 is the amplitude of the plane wave in vacuum which is a constant within the aperture, D represents the region of the aperture, and SA is the area of the aperture. In the Rytov approximation, it follows that , where χ1,pl (r) and χ2,pl (r) denote the first- and second-order log-amplitude fluctuations of a plane wave, respectively; 〈χ1,pl (r)〉 = 0, 〈χ2,pl (r)〉 is of the same order of magnitude as , and the angle brackets represent an ensemble average. The second-order Rytov perturbation is essential to evaluation of various statistical moments of optical fields; however, it can be omitted in calculation of the logamplitude variance and irradiance correlation function. In this paper, we only focus on statistics related to the log-amplitude variance and irradiance correlation function. Hence, thereinafter, we leave the terms , χ2,pl (r), ⋯ out of I (r)/I0. With these in mind, one finds where represents the aperture-averaged value of χ1,pl (r). Equation (9) manifests the relationship between the instantaneous normalized light flux and aperture-averaged first-order log-amplitude fluctuation. In accordance with the statistical theory, the correlation coefficient between light-flux fluctuations of two counter-propagating plane waves can be formulated by where , and ; and represent the aperture-averaged values of χ1,pl,F (r, z = L) and χ1,pl,I (ρ, z = 0), respectively. Notice that, in fact, the cross-covariance and variances of normalized light-flux fluctuations of two counter-propagating plane waves are 4 times Bpl,FI, and , respectively; similar relation holds true for the case of two counter-propagating spherical waves that will be addressed later. It should be emphasized here that both and have a mean value of zero and are real, implying and which can be readily deduced from the definition of the first-order log-amplitude fluctuation. However, for convenience of later mathematical treatments, here we formally use the complex-conjugate operation. and can be expressed by with where Ra denotes the radius of the aperture used to collect the wave. By introducing Eqs. (6) and (7) into Eqs. (11) and (12), respectively, and using the Fourier transform relation where the jinc function is defined by jinc (x) = 2J1 (x) /x with J1 (·) being a Bessel function of the first kind, one finds With Eqs. (15) and (16) in hand, one can develop the expressions for Bpl,FI, and . After some tedious mathematical manipulations, we find where Φn (κ, z) represents the three-dimensional turbulence spectrum. As in the mathematical treatment done by Tatarski [19] and Ishimaru [20], in arriving at Eqs. (17) – (19), we have used the following: the first is the relationship with Fn (·) representing the two-dimensional turbulence spectrum and δ () being the delta function; the second is the fact that Fn () takes a negligible value when |z − z′| is beyond the correlation distance of refractive-index fluctuations; the third is the assumption that atmospheric turbulence is statistically homogeneous and isotropic in each plane perpendicular to the z-axis. If we let Ra = 0, and become the variances of log-amplitude fluctuations for the forward- and inverse-propagating waves, respectively. It is straightforward to find that our result for with Ra = 0 is consistent with the formula for the variance of plane-wave log-amplitude fluctuations presented in [19].According to the existing literature [19–21], the expression for the turbulence spectrum can generally be separated into two terms, i.e., , where the first term, referred to as the refractive-index structure constant, characterizes the strength of refractive-index fluctuations at the path location z, and the second term depicts the shape of the spectrum. With this observation in mind, by interchanging the order of integrations in Eqs. (17) – (19), we can elucidate the underlying physics associated with light-flux fluctuations of two counter-propagating plane waves from different perspectives.
2.1.1. Wavenumber domain perspective
When is independent of the path location z, we can first evaluate the integrations on z in Eqs. (17) – (19) and hence eliminate the dependence of Bpl,FI, and on z. Performing this operation permits us to express Bpl,FI, and in the forms of an integral over the wavenumber domain, which are given by
with where is the scaled wavenumber, is called the aperture Fresnel parameter, and the sinc function is defined by sinc (x) = sin (x)/x. It is straightforward to find that fpl(·) and f′pl(·) actually act as two spectral weighting functions, which may emphasize or deemphasize the role of refractive-index fluctuations associated with certain spatial wavenumber κ, implying that they can filter a specific portion of appearing in Eqs. (17) – Eqs. (21) and Eqs. (17) – (22). In this sense, they can be regarded as spectral filter functions.We note that fpl(·) and f′pl(·) factorize into a product of two different terms. The first term emphasizes the role of refractive-index fluctuations associated with the scaled wavenumber and basically neglects those associated with the scaled wavenumber . This statement becomes apparent if we recognize the approximation jinc2(x) ≈ exp(− β2x2) with β = 0.5216 [24]. With this in mind, one can deduce that the portion of with κ > 1/(βRa) actually contributes little to both Bpl,FI and . Indeed, the term formulates the role of a receiving aperture in the treatment of the scintillation problem. It is obvious that a greater Ra leads to a smaller portion of that can produce scintillations effectively. This is the essential reason for the reduction in scintillations induced by aperture averaging. Strictly speaking, the use of the aforementioned approximation of jinc2 (·) means that the “hard aperture” described by Eqs. (13) is replaced by a Gaussian aperture, which has been previously employed by Kon [25]; doing so will facilitate our theoretical development. It is noted that fpl() is distinguished from f′pl(·) only by the two terms hpl(·) and h′pl(·). To obtain some specific idea about the difference between Bpl,FI and , below we compare the behaviors of hpl(·) and h′pl(·) in terms of . Figure 2 shows the shapes of the functions hpl(·) and h′pl(). It is found from Fig. 2 that, with an increasing value of , both hpl(·) and h′pl () first grow monotonically from 0 to a peak in the vicinity of ; after that hpl(·) begins to oscillate around 0 with almost constant amplitude and h′pl(), however, begins to oscillate around 1 with diminishing amplitude. Hence, refractive-index·fluctuations associated with all spatial wavenumbers make positive contribution to , whereas there may exist refractive-index fluctuations related to some spatial wavenumbers that make negative contribution to Bpl,FI. It is apparent that , implying μpl < 1, because the curve corresponding to h′pl(·) basically always lies above that corresponding to hpl(·). As a result, perfect correlation between light-flux fluctuations of two counter-propagating plane waves in atmospheric turbulence is impossible in general cases. Moreover, by recalling the aforesaid approximation of jinc2(·) and the curves shown in Fig. 2, it is straightforward to infer that, if , the specific behaviors of hpl(·) and h′pl(·) beyond only play a trivial role in determining Bpl,FI, and ; it is also noted·that and for . Hence, we can reason from the above analysis that and thus μpl ≃ 1/2 when , irrespective of what mathematical form takes.
To gain more details about μpl, we need to consider the specifics of the turbulence spectrum. In the literature, the von Kármán spectrum is usually used to describe the refractive-index fluctuations, and its spectral shape is characterized by , where κ0 = 2π/L0 and κm = 5.92/l0 with l0 and L0 being the inner and outer scales of turbulence, respectively. However, examination of the integrands in Eqs. (17) – (19) reveals that L0 can be reasonably specified by ∞ in our case because both Ra and are much smaller than L0 for typical optical-wave propagation scenarios, implying that indeed we can use the Tatarskii spectrum, which is equivalent to the von Kármán spectrum with L0 → ∞, in the following theoretical development. To facilitate the mathematical treatment, we divide hpl(·) into two terms, i.e., with . Accordingly, Bpl,FI is expressed as Bpl,FI = b1 − b2 with b1 given by
where qp = qm + β2qR, is called the inner-scale Fresnel parameter, , and . We have used the aforesaid approximation of jinc2 (·) in the second step of Eqs. (27), meaning that the Gaussian-aperture approximation has been employed; indeed, this is the reason for the fact that qm and qR can be combined together to form qp. The term b2 is expressed by This result is obtained by recognizing that, with the use of the said approximation of jinc2 (·), the integral in the first step of Eqs. (28) is functionally of a form akin to that dealt with by Ishimaru ([20], Eqs. (17–97) and (17–98)). In a similar way, one can work out the integral associated with to yield With the above results in hand, we find where, on the left-hand side, we explicitly show the dependence of μpl on the nondimensional parameter qp.Equation (30) shows our closed-form expression for the correlation coefficient between light-flux fluctuations of two counter-propagating plane waves in weak turbulence. It is worth noting that μpl is actually completely determined by the nondimensional parameter qp. By examining the mathematical form of qp, we can deduce that the inner scale l0 plays a role similar to that the aperture radius Ra does on the correlation coefficient μpl. It is easy to find that μpl ≃ 0.468 when qp → 0, i.e., both qm → 0 and qR → 0. Further, it should be pointed out that the use of the aforementioned approximation of jinc2 (·) has actually no impact on the value of μpl with qp → 0, thus meaning qR → 0, because with qR = 0. Figure 3 shows the correlation coefficient μpl in terms of under the condition that the turbulence spectrum has no path dependence. It is found from Fig. 3 that μpl is peaked at with its maximum value equal to 0.574, and μpl achieves its minimum value at . Furthermore, μpl tends to 0.5 when , which is consistent with the previous deduction (i.e., μpl ≃ ½ when ). The above analyses and calculations provide us a clear understanding of the effects of various propagation parameters on μpl when atmospheric turbulence is statistically homogeneous over the path.
2.1.2. Path-space domain perspective
Unlike the treatment above, at this point, we first evaluate the integrations on κ in Eqs. (17) – (19) and then express Bpl,FI, and in the forms of an integral over the normalized path-space domain, which are given by
with ξ = z/L being the normalized path location, , , , and In arriving at Eqs. (31) – (34), the said approximation of jinc2 (·) has been utilized.It is worth noting that gpl(ξ), g′pl(ξ) and g″pl(ξ) play the role of a path weighting function which emphasizes or deemphasizes the contributions of turbulent eddies at different locations to Bpl,FI, and , respectively. The path weighting functions appearing in Eqs. (31) – (33) in terms of ξ are plotted in Fig. 4 with various qp. It is seen from Fig. 4 that gpl(ξ) is symmetric about ξ = 0.5, and g′pl(ξ) is the reflection of g″pl(ξ) in the vertical line ξ = 0.5. Examination of Fig. 4 reveals that, on the one hand, turbulent eddies located in the vicinity of the path midpoint are most effective in increasing Bpl,FI and those located in the vicinities of the two path endpoints nearly do not contribute to Bpl,FI ; on the other hand, turbulent eddies near the transmitting plane are weighted most strongly for and , and those near the receiving plane contribute little to and . By recalling Eqs. (10), a statement can be made that turbulent eddies near the path midpoint and endpoints tend to make positive and negative contributions to μpl, respectively. Note that, the integrals in Eqs. (31) – (33) can be regarded as an inner product of the two terms contained in the integrands. It is obvious that the shape of in terms of ξ has an important impact on the results of the inner products, implying that μpl will depend on the path variations of . If the strength of refractive-index fluctuations over the path featured a single large sharp peak at the midpoint, μpl would be likely to take a value close to 1. Indeed, it is observed from Fig. 4 that, with the same qp, the curves corresponding to gpl(ξ), g′pl(ξ) and g″pl(ξ), respectively, come together at ξ = 0.5. This phenomenon is consistent with what could be expected by intuition; i.e., if only one random phase screen with a negligible thickness is positioned at the midpoint of a path in vacuum, one can find μpl ≡ 1. In the same context, if the random phase screen is moved away from the midpoint, we can infer from Fig. 4 that μpl < 1 and . By the way, there are researches on experimental simulations of laser beam propagation through atmospheric turbulence in the laboratory by using a liquid-crystal spatial light modulator acting as a random phase screen [26]; according to the above analysis, it becomes apparent that the position where the random phase screen is located plays an important role in the simulated optical scintillations, meaning that the position of the random phase screen needs to be chosen carefully in experiments in order to properly emulate the target propagation scenarios. Finally, one can see from Fig. 4 that a larger qp leads to smaller magnitude of the path weighting functions. This means that an increase in qp generally lowers the absolute magnitude of contributions from turbulent eddies over the path to Bpl,FI, and .
2.2. Two counter-propagating spherical waves
Now we turn to the situation of two counter-propagating spherical waves. Following an approach similar to that used in Sec. 2.1, the aperture-averaged values of the first-order log-amplitude fluctuations of the forward- and inverse-propagating spherical waves can be formulated by
respectively, where Based on Eqs. (35) and (36), analogous to the plane-wave case, we can further develop the expressions for , and , which take the forms The correlation coefficient between light-flux fluctuations of two counter-propagating spherical waves is formulated by2.2.1. Wavenumber domain perspective
If the turbulence spectrum Φn (·) has no path dependence, we can reformulate Bsp,FI, and as follows:
with Equations (45) and (46) show the spectral filter functions for Bsp,FI, and , respectively. When Ra = 0, i.e., qR = 0, it follows that , implying that and thus μsp ≡ 1, regardless of the mathematical form of . This result is indeed what we will expect from the reciprocity principle [9–11] by recognizing that Ra = 0 means a point receiver. Under general conditions, it is formidable to analytically evaluate the integral of Eq. (45) in a simple closed form. An analytical expression for Eq. (46) can be developed if the said approximation of jinc2 (·) is employed; however, here we do not show it in detail. Thereinafter, this approximation will not be used. In addition, for the same reason stated in Sec. 2.1, we utilize the Tatarskii spectrum in the following theoretical analysis.The profiles of and are illustrated in Fig. 5 with various qR. It is found from Fig. 5 that the curves corresponding to and , respectively, basically coincide with each other when qR = 0.04 and obviously deviate from each other within a range of the scaled wavenumber when qR = 1 or 16. In fact, it is the deviation of from that causes μsp to decrease below 1. Further, a larger qR moves the said deviation deeper into the low-wavenumber range. Notice that, has a greater value with a smaller κ in the inertial range. With these observations in mind, it is straightforward to deduce that a larger qR generally leads to a smaller μsp. On the other hand, as for the two curves associated with qR = 1 in Fig. 5, it is seen that they begin to obviously deviate from each other only when . Hence, if κm = 5.92/l0 ≲ 2/(L/k)1/2, the deviation portion of these two curves only has a negligible impact on μsp because rapidly drops close to 0 when κ increases beyond κm. This leads us to an inference that a nonzero inner scale l0 may make the decrease in μsp caused by a nonzero aperture radius Ra negligible. As a result, it is necessary to consider the role of a nonzero inner scale l0 besides the receiver aperture radius Ra in determining μsp under some conditions. Moreover, it is seen from Fig. 5 that always takes a positive value, whereas with qR > 0 may take a negative value at some scaled wavenumbers. Thus, like the plane-wave case, refractive-index fluctuations related to some spatial wavenumbers may make negative contributions to Bsp,FI.
According to Eqs. (43) – (46), one can examine the correlation coefficient in the spherical-wave case by making use of numerical calculations. By considering the Tatarskii spectrum and making the change of variable κ = t/(L/k)1/2 in Eqs. (43) and (44), it is easy to find that μsp is completely determined by the two nondimensional parameters qm and qR. Figure 6 depicts the contours of μ as a function of and . It is seen from Fig. 6 that, with the same qm, μsp decreases when qR increases; in contrast, with the same qR, a larger qm results in a greater μsp. Hence, an increase in qm may cancel out the reduction in μsp caused by enlarging the value of qR. The underlying physical reason for this behavior of μsp has been elaborated in the analysis of Fig. 5. Comparison of these results with those in the plane-wave case shows that the inner scale l0 plays a different role in determining the correlation coefficient. Notice that, Fig. 6 indeed illustrates the behavior of μsp in most parameter regions of practical interest to us. The contour lines in Fig. 6 are nearly straight lines except for the parts corresponding to relatively small and .
At this point, a question that needs to be addressed further is what value μsp approaches when qR → ∞. By introducing Eqs. (45) and (46), respectively, into Eqs. (43) and (44), taking the Tatarskii spectrum into account and making the change of variable κ = t′/Ra = t′/(qR L/k)1/2, for a large enough qR, it follows that exp(−qmt′2/qR) ≈ 1 and sinc[ξ (1 − ξ) t′2/(2qR)] ≈ 1 within most of the integration range of practical interest. With this observation in mind, we can deduce that and when qR ≫ 1, where denotes the spherical-wave log-amplitude variance. Consequently, μsp approaches 0.122 as qR tends to infinity. For two counter-propagating spherical waves in weak atmospheric turbulence, Perlot and Giggenbach inferred the same conclusion from their numerical results (see page 2891 in [3]). In addition, it is easy to find from our above results that the aperture-averaging factor for spherical-wave scintillations with qR ≫ 1 actually approaches ; a similar result has been presented in [19] previously.
2.2.2. Path-space domain perspective
Like the plane-wave case, Bsp,FI, and can also be expressed in the forms of an integral over the normalized path-space domain, i.e.,
where , , and withIt can be deduced from Eq. (50) that both qm and qR play a role in determining the shapes of gsp(ξ), g′sp(ξ) and g″sp(ξ) as a function of ξ. Different from the plane-wave case, qm and qR in Eq. (50) cannot be combined together to define a single nondimensional parameter. It is obvious that gsp(ξ) ≡ g′sp(ξ) ≡ g″sp(ξ) when qR = 0. This indeed corresponds to the case of the PSPR reciprocity. Figure 7 demonstrates the path weighting functions gsp(ξ), g′sp(ξ) and g″sp(ξ) in terms of the normalized path location ξ with various qm and qR. Similar to the plane-wave case, in general, gsp(ξ) always features a symmetric shape with its peak located at ξ = 0.5 and g′sp(ξ) is the reflection of g″sp(ξ) in the vertical line ξ = 0.5; further, g′sp(ξ) and g″sp(ξ) are peaked at a path location dependent on qm and qR. It is seen from Fig. 7 that gsp(ξ), g′sp(ξ) and g″sp(ξ) basically take on the same shape with qR much smaller than 1, implying that μsp is close to 1. One can observe that, as qR grows larger with qm fixed, although gsp(ξ) remains peaked at ξ = 0.5, the peaks of g′sp(ξ) and g″sp(ξ) gradually move from the midpoint to the left and right sides of the abscissa, respectively; this behavior means that turbulent eddies at the path midpoint are always weighted most heavily for Bsp,FI, whereas the location of turbulent eddies that are weighted most strongly for or gradually moves toward the transmitting plane with an increasing value of qR. In addition, it can be seen from Fig. 7 that qm produces effects opposite to qR, which has been found in the preceding analysis. This reveals that, from comprehensive points of view, the inner scale of turbulence should be considered in the treatment of the correlation between light-flux fluctuations of two counter-propagating waves in atmospheric turbulence. Finally, like the plane-wave case, the shape of in terms of ξ obviously plays an important role in determining μsp; if featured a single large sharp peak at the midpoint of a path, μsp would be close to 1 even with a somewhat great qR.
3. Discussions
The optical signal collected by the aperture of a PIB receiver can be viewed as the light flux entering it [19,27–29], which is equal to the integrated irradiance over the aperture. The problem as to aperture-collected optical scintillations has been commonly treated based on the spatial covariance function of irradiance fluctuations. In the previous section, we employed the relation given by Eq. (9) to relate the normalized light flux with the aperture-averaged first-order log-amplitude fluctuation, and subsequently developed the expressions for the correlation coefficient between light-flux fluctuations of two counter-propagating plane or spherical waves in weak turbulence. In our theoretical treatment, by utilizing the Fourier transform, the role of a finite receiving aperture is formulated as the spectral filtering factors presented by Eq. (14) and (38). This permits us to use a straightforward procedure to deal with the light-flux fluctuations.
In Sec. 2, we have occasionally analyzed the asymptotic behavior of the correlation coefficient between light-flux fluctuations of two counter-propagating waves. Charnotskii [28] presented a method for performing asymptotic analysis of the variance of light-flux fluctuations. Thereinafter, we follow the procedure used by this method to find a full set of asymptotes of the correlation coefficient. By considering the Tatarskii spectrum with a refractive-index structure constant independent of the path location, Eqs. (17) – (19) and (39) – (41) can be rewritten in the form
Specifically, Bpl,FI = G (1, 1, ξ, 1 − ξ), , , Bsp,FI = G (ξ, 1 − ξ, ξ (1 − ξ), ξ (1 − ξ)), , and . There are three spatial scales included in Eq. (51), i.e., Ra, and . Three asymptotic cases are related to three different situations in which one of these spatial scales determines the effective integration domain of the inside integral in Eq. (51). They are given as follows. When qm ≪ 1 and qR ≪ 1, one finds When qm ≪ qR and qR ≫ 1, one finds When qm ≫ qR and qm ≫ 1, one findsA 2D map of asymptotes, previously introduced in [28], for our case is plotted in Fig. 8(a) by combining the above three asymptotes together. From Fig. 8(a), we can verify that the three asymptotes completely cover the parameter space. Asymptotes of μpl and μsp found immediately from Eqs. (52) – (54) are exhibited by Figs. 8(b) and 8(c). Indeed, Fig. 6 reflects the structure of Fig. 8(c). For many practical optical-wave propagation scenarios, qR is larger than qm; it can be deduced from Figs. 8(b) and 8(c) that the inner scale tends to become unimportant if qR > qm. Additionally, from Figs. 8(b) and 8(c), one can find that in many situations the correlation between light-flux fluctuations of two counter-propagating waves is lower than 1 obviously. However, as shown in [1,30], good correlation between two counter-propagating waves measured by SMFC transceivers can be preserved for large receiving apertures even under strong-turbulence conditions; as far as the preservation of good correlation is concerned, bidirectional FSO systems based on SMFC transceivers appear more advantageous than those based on PIB receivers.
As an important comment, in arriving at Eq. (8), the irradiance in vacuum is indeed thought of as being independent of r. This is reasonable for the plane- and spherical-wave cases under the paraxial approximation. Because the irradiance of a collimated-Gaussian-beam wave in vacuum does depend on r, strictly speaking, the result shown by Eq. (8) is not completely applicable to the collimated-Gaussian-beam cases. However, if the beam radius at the receiving plane is much greater than the aperture radius, the variation of the irradiance in the absence of atmospheric turbulence within the aperture will be small, meaning that the use of Eq. (8) to estimate the normalized light flux of a collimated-Gaussian-beam wave will not incur a significant error. It has been stated previously that the results obtained in Sec. 2 may provide some somewhat reasonable estimates for the collimated-Gaussian-beam cases under the near-and far-field conditions; with the above discussion in mind, it is apparent that an additional condition necessary for this statement is that the aperture radius is much smaller than the beam radius at the receiving plane.
4. Conclusions
We have formulated the correlation coefficient between light-flux fluctuations of two counter-propagating plane or spherical waves in weak atmospheric turbulence. Based on the obtained formulae, by considering the Tatarskii turbulence spectrum, we have elucidated the behavior of the correlation coefficient from both the wavenumber domain and path-space domain perspectives. If the turbulence spectrum has no path dependence, the correlation coefficient is completely determined by the aperture and inner-scale Fresnel parameters. If the turbulence spectrum is dependent on path locations, the aperture and inner-scale Fresnel parameters completely determine the shape of the path weighting functions for the cross-covariance and variances of normalized light-flux fluctuations. For two counter-propagating plane waves, turbulent eddies located in the vicinities of the path midpoint and endpoints are most effective in making positive and negative contributions to the correlation coefficient, respectively. On the other hand, for two counter-propagating spherical waves, turbulent eddies located in the vicinity of the path midpoint are still weighted most strongly in producing the correlation, whereas the role of those near the path endpoints on the correlation coefficient depends on the value of the aperture Fresnel parameter.
Funding
National Natural Science Foundation of China (61007046, 61275080 and 61475025); Natural Science Foundation of Jilin Province of China (20150101016JC); Specialized Research Fund for the Doctoral Program of Higher Education of China (20132216110002).
Acknowledgments
The authors are very grateful to the reviewers for valuable comments.
References and links
1. J. Minet, M. A. Vorontsov, E. Polnau, and D. Dolfi, “Enhanced correlation of received power-signal fluctuations in bidirectional optical links,” J. Opt. 15(2), 022401 (2013). [CrossRef]
2. M. A. Vorontsov, “Conservation laws for counter-propagating optical waves in atmospheric turbulence with application to directed energy and laser communications,” in Proceedings of Imaging and Applied Optics, OSA Technical Digest (Optical Society of America, 2013), paper PW3F.1.
3. N. Perlot and D. Giggenbach, “Scintillation correlation between forward and return spherical waves,” Appl. Opt. 51(15), 2888–2893 (2012). [CrossRef] [PubMed]
4. J. H. Shapiro and A. L. Puryear, “Reciprocity-enhanced optical communication through atmospheric turbulence -part I: reciprocity proofs and far-field power transfer optimization,” J. Opt. Commun. Netw. 4(12), 947–954 (2012). [CrossRef]
5. L. Puryear, J. H. Shapiro, and R. R. Parenti, “Reciprocity-enhanced optical communication through atmospheric turbulence - part II: communication architectures and performance,” J. Opt. Commun. Netw. 5(8), 888–900 (2013). [CrossRef]
6. C. Chen, H. Yang, S. Tong, and Y. Lou, “Mean-square angle-of-arrival difference between two counter-propagating spherical waves in the presence of atmospheric turbulence,” Opt. Express 23(19), 24657–24668 (2015). [CrossRef] [PubMed]
7. A. Khatoon, W. G. Cowley, and N. Letzepis, “Capacity of adaptive free-space optical channel using bi-directional links,” Proc. SPIE 8517, 85170X (2012). [CrossRef]
8. N. Wang, X. Song, J. Cheng, and V. C. M. Leung, “Enhancing the security of free-space optical communications with secret sharing and key agreement,” J. Opt. Commun. Netw. 6(12), 1072–1081 (2014). [CrossRef]
9. J. H. Shapiro, “Reciprocity of the turbulent atmosphere,” J. Opt. Soc. Am. 61(4), 492–495 (1971). [CrossRef]
10. R. F. Lutomirski and H. T. Yura, “Propagation of a finite optical beam in an inhomogeneous medium,” Appl. Opt. 10(7), 1652–1658 (1971). [CrossRef] [PubMed]
11. V. P. Lukin and M. I. Charnotskii, “Reciprocity principle and adaptive control of optical radiation parameters,” Sov. J. Quantum Electron. 12(5), 602–605 (1982). [CrossRef]
12. D. Giggenbach, W. Cowley, K. Grant, and N. Perlot, “Experimental verification of the limits of optical channel intensity reciprocity,” Appl. Opt. 51(16), 3145–3152 (2012). [CrossRef] [PubMed]
13. R. R. Parenti, J. M. Roth, J. H. Shapiro, F. G. Walther, and J. A. Greco, “Experimental observations of channel reciprocity in single-mode free-space optical links,” Opt. Express 20(19), 21635–21644 (2012). [CrossRef] [PubMed]
14. P. J. Winzer and W. R. Leeb, “Fiber coupling efficiency for random light and its applications to lidar,” Opt. Lett. 23(13), 986–988 (1998). [CrossRef]
15. Y. Dikmelik and F. M. Davidson, “Fiber-coupling efficiency for free-space optical communication through atmospheric turbulence,” Appl. Opt. 44(23), 4946–4952 (2005). [CrossRef] [PubMed]
16. M. Charnotskii, “Coupling turbulence-distorted wave front to fiber: wave propagation theory perspective,” Proc. SPIE 7814, 78140I (2010). [CrossRef]
17. C. Chen, H. Yang, H. Wang, S. Tong, and Y. Lou, “Coupling plane wave received by an annular aperture into a single-mode fiber in the presence of atmospheric turbulence,” Appl. Opt. 50(3), 307–312 (2011). [CrossRef] [PubMed]
18. C. Chen and H. Yang, “Characterizing the radial content of orbital-angular-momentum photonic states impaired by weak-to-strong atmospheric turbulence,” Opt. Express 24(17), 19713–19727 (2016). [CrossRef] [PubMed]
19. V. I. Tatarskii, Wave Propagation in a Turbulent Medium (McGraw-Hill, 1961).
20. A. Ishimaru, Wave Propagation and Scattering in Random Media, Vol. 2 (Academic, 1978).
21. L. C. Andrews and R. L. Phillips, Laser Beam Propagation through Random Media, 2nd ed. (SPIE, 2005). [CrossRef]
22. M. I. Charnotskii, “Asymptotic analysis of finite-beam scintillations in a turbulent medium,” Waves Random Media 4(3), 243–273 (1994). [CrossRef]
23. M. I. Charnotskii and G. J. Baker, “Practical calculation of the beam scintillation index based on the rigorous asymptotic propagation theory,” Proc. SPIE 8038, 803804 (2011). [CrossRef]
24. Y. Cheon and A. Muschinski, “Closed-form approximations for the angle-of-arrival variance of plane and spherical waves propagating through homogeneous and isotropic turbulence,” J. Opt. Soc. Am. A. 24(2), 415–422 (2007). [CrossRef]
25. A. I. Kon, “Focusing of light in a turbulent medium,” Radiophys. Quantum. Electron. 13(1), 43–50 (1970). [CrossRef]
26. J. D. Schmidt, M. E. Goda, and B. D. Duncan, “Emulating bulk turbulence with a liquid-crystal spatial light modulator,” Proc. SPIE 6306, 63060O (2006). [CrossRef]
27. L. C. Andrews, “Aperture-averaging factor for optical scintillations of plane and spherical waves in the atmosphere,” J. Opt. Soc. Am. A 9(4), 597–600(1992). [CrossRef]
28. M. I. Charnotskii, “Asymptotic analysis of the flux fluctuations averaging and finite-size source scintillations in random media,” Waves Random Media 1(4), 223–243 (1991). [CrossRef]
29. C. Chen, H. Yang, Z. Zhou, W. Zhang, M. Kavehrad, S. Tong, and T. Wang, “Effects of source spatial partial coherence on temporal fade statistics of irradiance flux in free-space optical links through atmospheric turbulence,” Opt. Express 21(24), 29731–29743 (2013). [CrossRef]
30. M. A. Vorontsov, S. L. Lachinova, and A. K. Majumdar, “Target-in-the-loop remote sensing of laser beam and atmospheric turbulence characteristics,” Appl. Opt. 55(19), 5172–5179 (2016). [CrossRef] [PubMed]