Focal plane analysis of optical power-related metrics for atmospheric turbulence-affected laser beams

Abstract. The analysis of budget link and free space optical system performances requires the calculation of several metrics of the atmospheric turbulence-affected collected light. In this work, assuming the collected light is focused into an optical fiber or over a sensor positioned in the focal plane, we use the ABCD ray-matrix representation to calculate the impact of atmospheric turbulence on the power in the fiber or power over the sensor. Calculation of such metrics requires the knowledge of the transmitted average power that enters the receiver aperture (power in the bucket) and the long-term beam spread in the focal plane, from which the Strehl ratio can be obtained.


Introduction
2][3] This parameter holds significance in various applications, particularly in directed energy scenarios.In the realm of free-space optical communications, it becomes essential to extend the theoretical modeling of the long-term beam spread to the focal plane, where an optical detector or an optical fiber is commonly positioned. 4his paper extends the analysis of the long-term beam spread and Strehl ratio (SR) in the focal plane, as explained in Ref. 4, to conduct an examination of other optical power-related important metrics for budget link calculations.In particular, we calculate the average power in the bucket (PIB) captured by the receiver aperture and, depending on the presence of a fiber or a sensor in the focal plane, the average power into the fiber (PIF) or the average power over the sensor (POS).
We perform this analysis for Gaussian beams with general geometries, encompassing focused, collimated, and divergent configurations.Our methodology relies on the ABCD matrix representation 1 and accommodates scenarios where turbulence can deviate from the assumptions of the Kolmogorov model. 2

Kolmogorov and Non-Kolmogorov Turbulence
Let us focus on the following power law-dependent non-Kolmogorov power spectrum: 2 E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 1 ; 1 1 4 ; 5 4 3 where AðαÞ ¼ Γðα−1Þ 4π 2 cosðα π 2 Þ, κ ≡ ðκ x ; κ y ; κ z Þ is the spatial wavenumber vector, α is the power law, C2 n ¼ βðαÞ • C 2 n is the generalized structure parameter with units ½m 3−α , βðαÞ is a constant depending on α and has units ½m 11∕3−α , and symbol ΓðxÞ denotes the Gamma function.When power law assumes value α ¼ 11∕3, the generalized structure parameter reduces to the structure parameter C 2 n with units ½m −2∕3 and Eq. ( 1) reduces to the Kolmogorov power spectrum, Φ n ðκÞ ¼ 0.033 The schematic of a free space optical communication system is shown in Fig. 1.A Gaussian beam with spot radius W 0 exits the transmitter aperture with geometry defined by a couple of beam parameters Here L is the path length, F 0 is the radius of curvature of the wavefront at the transmitter, and k is the wavenumber.After propagation, the beam is captured at the receiver aperture by a collecting lens of diameter D G .The beam on the collecting lens has geometry described by and it is focused at distance L f behind the lens, where a detector or an optical fiber is located.
Note that here we consider a hard aperture with diameter D G related to the soft aperture Gaussian lens of radius W G by the known expression ).The beam in the focal plane has the geometry ðΛ is the focal length of the collecting lens.The diffraction-limited beam spot radius, Note that all the geometric parameters, Θ 0 ; Λ 0 ; Θ 1 ; Λ 1 and Θ 2 ; Λ 2 are free space parameters, however they enter the equation describing the change in the on-axis mean irradiance in the focal plane after propagation in turbulence, Eq. (8).

Long Term Beam Spread Analysis in the Focal Plane
The long-term beam spread is the parameter that physically describes the mean radius of the Gaussian beam after its propagation in optical turbulence.The long-term spot radius in the focal plane can be defined as in Refs. 1, 2, and 4 E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 2 ; 1 1 4 ; 1 4 6 where T focal describes the change in the on-axis mean irradiance in the focal plane and characterizes the beam spread due to turbulence on a long-term average.Indeed, such a term includes the spread due to beam wander and small-scale turbulence cells diffraction.Toselli: Focal plane analysis of optical power-related metrics for atmospheric. . .
As previously introduced, to calculate T focal we use the ABCD ray-matrix representation shown in Ref. 1, and we neglect the effect of turbulence on the beam behind the lens (beam path from the lens to the detector).Under this condition, we use two statistical terms 1 E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 3 ; 1 1 7 ; 7 0 0 (3) and E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 4 ; 1 1 7 ; 6 5 3 where k ¼ 2π∕λ is wave number, λ is the wave length, ξ ¼ 1 − z∕L is the normalized z-distance.
For a beam focused on the photodetector/optical fiber we impose 1 E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 5 ; 1 1 7 ; 5 9 1 and γ 1 and B 1 ðξÞ reduce to [see Eq. ( 9), page 398 in Ref. 1 for details] E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 6 ; 1 1 7 ; 5 4 3 where Therefore, the exponential term can be expressed as ; t e m p : i n t r a l i n k -; e 0 0 7 ; 1 1 7 ; 4 8 9 and we obtain E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 8 ; 1 1 7 ; 4 3 9 which means collimated beam incidents the collecting lens.For the general case of a beam not arriving collimated on the lens ðΘ 1 ≠ 1Þ, to insure the focusing on sensor/optical fiber, the distance behind the lens has to be We remark that Eq. (3) diverges using the spectrum of Eq. ( 1) because outer scale is supposed to be infinite.However, it can be shown that such a singularity is removed in Eq. ( 8) and a finite result is found.

Kolmogorov Turbulence Case
Introducing Eq. ( 1) in Eq. ( 8) with α ¼ 11∕3, we obtain the Kolmogorov turbulence case result for Gaussian beams E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 9 ; 1 1 7 ; 2 3 8 where σ 2 R is the Rytov variance.For a beam arriving collimated on the collecting lens, Θ 1 ¼ 1, the integral in Eq. ( 9) can be expressed in closed form using the hypergeometric function 2 F 1 ða; b; c; xÞ as E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 1 0 ; 1 1 7 ; 1 6 3 For a plane wave incident on the collecting lens, ðΛ ; t e m p : i n t r a l i n k -; e 0 1 1 ; 1 1 7 ; 1 1 3 Finally, for a spherical wave incident on the collecting lens, ðΛ 1 ¼ 0; Θ 1 ¼ 0Þ we obtain the same result as in Ref. 1 E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 1 2 ; 1 1 4 ; 7 1 2 is the spherical Fried parameter.For this specific spherical wave beam geometry, Andrews 1,3 modified Eq. ( 12) to include also the case when tip and tilt Zernike modes are both removed from the wavefront.
He introduced the coefficient 0.28 in the equation obtaining

Non-Kolmogorov Turbulence Case
Using the non-Kolmogorov turbulence power spectrum, Eq. ( 1) we obtain ; t e m p : i n t r a l i n k -; e 0 1 3 ; 1 1 4 ; 5 7 0 is the non-Kolmogorov Rytov variance defined in Ref. 2.
For a beam arriving collimated on the collecting lens, ; t e m p : i n t r a l i n k -; e 0 1 4 ; 1 1 4 ; 4 9 4 For a plane wave incident on the collecting lens, ðΛ 1 ¼ 0; Θ 1 ¼ 1Þ and we obtain E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 1 5 ; 1 1 4 ; 4 4 3 For a spherical wave incident on the collecting lens, ðΛ 1 ¼ 0; Θ 1 ¼ 0Þ and we obtain E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 1 6 ; 1 1 4 ; 3 8 6 Although we found expressions also for the non-Kolmogorov case, we show only results for Kolmogorov turbulence, α ¼ 11∕3.
For a specific set of parameters, we plot in Fig. 2 the term T focal as a function of the diffraction parameter, Λ 0 ¼ 2L kW 2 0 .Note that in our plots we change the spot size radius at the transmitter, W 0 from 0.001 m to 0.5 m covering the whole range from the near field Λ 0 ≪ 1 to the far , which is the radius of the first Fresnel zone.Note also that for the case of a collimated beam, the value Λ 0 ¼ 1 corresponds to the Rayleigh distance, z R ¼ 1 2 kW 2 0 , which forms a separation line between the near and far field.All the following plots will show only results for collimated and focused beams because of their relevance to free space communications.In terms of the formalism used in this paper, it should be clear that the collimated beam at the transmitter has a wavefront radius of curvature F 0 ¼ ∞; the focused beam has F 0 ¼ L. Note furthermore that F 0 ¼ ∞; Λ 0 ¼ 0 is equivalent to the plane wave case at the transmitter and F 0 ¼ L; Λ 0 → ∞ is equivalent to the spherical wave case.
We deduce from Fig. 2 that for a collimated beam and for this set of parameters, the spread term T focal is maximum in the near field.Note that in the region Λ 0 ≪ 1 a collimated beam (Θ 0 ¼ 1) is essentially a plane wave.This explains why the spread is higher: the radial dimension of the beam geometrically intercepts a higher number of turbulence cells along the path.The minimum value of T focal is reached at about 2.5 times the Fresnel distance (Λ 0 ¼ 1) and, for Λ 0 ≫ 1, T focal increases up to the saturation value double than unity (it approaches the spherical wave case, Λ 0 → ∞; Θ 1 → 0).On the other hand, a focused beam shows a lower spread than a collimated one all over the range from the near field to the far field, where it approaches again the spherical wave case Λ 0 → ∞; Θ 0 ¼ 0. As for the collimated beam case, this is physically explained by the radial dimension of a focused beam, which intercepts a lower number of turbulence cells along the path.
In Fig. 3, we plot, for the collimated beam case, two spot radii in the focal plane, respectively, the long-term beam spread W LT;f and the diffraction limited spot radius W 2 , both of them as a function of the diffraction parameter, Λ 0 ¼ 2L kW 2 0 .We deduce from Fig. 3 that in the near field Λ 0 ≪ 1 (where plane wave model holds well) the long-term beam spread in the focal plane is higher than in the far field Λ 0 ≫ 1, where the spherical wave model holds.This explains the saturation of both plots as well.Note also that the maximum diffraction-limited (no turbulence) spot radius W 2 can be obtained when the receiver is located at Fresnel distance, Λ 0 ¼ 1.However, in turbulence it is more beneficial to locate the receiver in the far field to reduce the long-term spot size of the beam.
This conclusion also applies to a focused beam as shown in Fig. 4.Here we plot the same parameters as in Fig. 3 but now for a beam focused on the receiver.Note how the value of the long-term beam spread in the focal plane is essentially driven by the diffraction parameter or, in other words, it is anchored on the diffraction-limited spot size W 2 .
To highlight the differences between the two beam geometries, we plot again in Fig. 5 the long-term beam spread in the focal plane W LT;f as a function of the diffraction parameter Λ 0 ¼ 2L kW 2 0 for both a collimated and a focused beam (these two plots are the blue curves shown also in Figs. 3 and 4).We deduce from Fig. 5 that in the near field a collimated beam shows a lower long-Fig.4 Long-term beam spread term W LT;f (blue) and the diffraction limited spot size W 2 (red) as a function of the diffraction parameter Λ 0 for a beam focused on the receiver (Kolmogorov turbulence case).Distance is fixed to the value L ¼ 2 km.term beam spread than a focused beam.However, they approach the same value when Λ 0 ¼ 1 (Fresnel distance) or higher.In general, Fig. 5 suggests that independently from using a collimated or focused beam, to obtain the smallest long-term spot size in the focal plane the receiver should be located in the far field, Λ 0 ≫ 1 (spherical wave model).We remark that in the far field Eq. ( 12) for spherical wave model holds well.

Strehl Ratio Analysis
The SR as one of the most used metrics at the focal plane is defined as 1 E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 1 7 ; 1 1 7 ; 6 3 1 where Ið0Þ is the on-axis intensity in the focal plane of the beam after propagation in free space (no turbulence) and I turb ð0Þ is the on-axis intensity of the beam in the focal plane after propagation in turbulence (we remark here that turbulence acts only from the transmitter to the receiver aperture, we neglect its effect behind the collecting lens).Using Eq. ( 9), the SR can be expressed as 1,3 E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 1 8 ; 1 1 7 ; 5 3 5 For the spherical wave case ðΛ 1 ¼ 0; Θ 1 ¼ 0Þ and Kolmogorov turbulence, α ¼ 11∕3, Eq. ( 18) reduces to the same result shown in Ref. 1

(see page 623)
E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 1 9 ; 1 1 7 ; 4 6 5 We plot in Fig. 6 the SR as a function of the diffraction parameter Λ 0 ¼ 2L

Power in the Bucket
The transmitted mean power of a Gaussian beam with peak intensity A 2 0 and spot radius W 2 0 is The fraction of the mean power entering the receiver aperture after propagation (distance L) is the PIB. 3  (20) where we ignored the atmospheric transmission loss caused by aerosols, etc., and the receiver optical element transmission loss.For a large beam (close to a plane wave) incident on the lens, an approximation is 3 where is the long term beam spread at the receiver aperture, W is the diffraction limited spot radius at distance L; and T pupil ðαÞ ¼ is the analogous of T focal ðαÞ at the receiver (pupil) aperture (see Refs. 1 and 2).Note that for Kolmogorov turbulence . We plot in Figs.7 and 8 the percentage of PIB scaled by the transmitted power, P Tx as a function of the diffraction parameter Λ 0 for several distances of the receiver for a collimated beam at the transmitter or focused on the receiver (Kolmogorov turbulence case only).We deduce from Figs. 7 and 8 that the maximum PIB is reachable when the beam is focused on the receiver (see Fig. 8) and the collecting lens is located in the near field.Also, when the receiver is positioned in the far field the spherical wave model [please see Eqs. ( 12) and ( 19)] holds well and there is essentially no difference in using either beam geometry.6 Power over the Sensor Similar to the PIB but now on the focal (detector) plane, the fraction of transmitted mean power illuminating the sensor/detector (POS) can be expressed as E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 2 2 ; 1 1 7 ; 4 3 5 where W S is the radius of the sensor (supposing it has a circular geometry) and, as previously mentioned, we ignored any additional loss, such as circulator loss and optical loss. 3upposing the beam illuminating the sensor is consistently larger than the sensor itself (the Gaussian beam across the sensor is almost a plane wave), an approximated expression of Eq. ( 22) is E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 2 3 ; 1 1 7 ; 3 3 9 We plot in Figs. 9 and 10, respectively, the percentage of POS scaled by the transmitted power, P Tx as a function of the diffraction parameter Λ 0 for several distances of the receiver for a collimated beam at the transmitter and for a beam focused on the receiver (Kolmogorov turbulence case only).
We deduce from Figs. 9 and 10 that, for both geometries and for this specific set of parameters, the maximum POS is reachable when Λ 0 assumes values in the interval from two to five times the Fresnel distance ðΛ 0 ¼ 1Þ, depending on path distances (note that the peak shifts slightly to the left at shorter distances).Also, there is essentially no difference of using the two different beam geometries (collimated or focused beam) in the far field (spherical wave model holds well in such a case).

Power into the Fiber
Supposing an optical fiber positioned in the focal plane, the mean power coupled into the fiber can be expressed as 3 where η turb ¼ η • SR focal , η turb is the fiber coupling efficiency including the turbulence-induced beam spread (it includes the reduction of power due to turbulence effects in focal plane, SR focal ) and η is the free space (defined without keeping into account the turbulence) fiber coupling efficiency.In our analysis, as already mentioned, we ignored any additional loss such as circulator loss and optical loss. 3e plot in Figs.11 and 12 the PIF as a function of the diffraction parameter Λ 0 for several distances of the receiver for a collimated beam at the transmitter or for a beam focused on the receiver (Kolmogorov turbulence case only).Similar to the POS in Fig. 9, we deduce from Fig. 11 (collimated beam) that, for this specific set of parameters, the maximum PIF is reachable  when Λ 0 assumes values in the interval from two to five times the Fresnel distance ðΛ 0 ¼ 1Þ, depending on path distances (note that the peak shifts slightly to the left at shorter distances).
Also, when the beam is focused on the receiver, we deduce form Fig. 12 that the maximum PIF (for this specific set of parameters) is reached mostly in the near field.Finally, similarly to previous metrics, we note that when the receiver is located in the far field the spherical wave model holds well and there is essentially no difference in using either beam geometry (collimated or focused beam case).

Conclusion
In this paper, we used the ABCD matrix formulation to obtain the analytical expressions of the main optical power-related metrics for a Gaussian beam after propagation through atmospheric  turbulence.Specifically, these optical power-related metrics are: the long-term beam spread in the focal plane, SR, PIB, POS, and PIF.We investigated those metrics as a function of the diffraction parameter, Λ 0 ¼ 2L kW 2 0 .For a specific set of parameters, we found that both the maximum SR and PIB are reachable when a beam focused on the receiver is used and the collecting lens is located in the near field.In addition, we found that for a collimated beam the maximum POS and PIF are reachable when the diffraction parameter Λ 0 assumes values in the interval two to five times the Fresnel distance ðΛ 0 ¼ 1Þ depending on path distances (the POS and PIF peaks shifts slightly to the left at shorter distances).However, when the beam is focused on the receiver, the maximum PIF is mostly reached in the near field.
Finally, in the far field and for all metrics analyzed in this paper, the spherical wave model holds well and there is essentially no difference of using the two different beam geometries (collimated or focused beam).Our results can be useful for the budget link analysis for free space optical communications.

Fig. 3
Fig.3Long-term beam spread term W LT;f (blue) and the diffraction limited spot radius W 2 (red) as a function of the diffraction parameter Λ 0 for a beam collimated at the transmitter (Kolmogorov turbulence case only).Distance is fixed to the value L ¼ 2 km.

Fig. 2
Fig. 2 Spread term T focal as a function of the diffraction parameter Λ 0 for a beam collimated at the transmitter (blue) or focused on the receiver (red) (Kolmogorov turbulence case).Distance is fixed to the value L ¼ 2 km.

Fig. 5
Fig. 5 Long-term beam spread term W LT;f as a function of the diffraction parameter Λ 0 for both beam geometries: collimated and focused (Kolmogorov turbulence case).Distance is fixed to the value L ¼ 2 km.

kW 2 0
Fig. 6 SR as a function of the diffraction parameter Λ 0 for a beam collimated at the transmitter (blue) or focused on the receiver (red) (Kolmogorov turbulence case only).Distance is fixed to the value L ¼ 2 km.

2 p
Considering a Gaussian lens (soft aperture) with radius W G related to the hard aperture diameter, D G by W G ¼ D G 2 ffiffi and supposing the turbulence affected beam still Gaussian, 1 the average hPIBi can be expressed as E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0

Fig. 7
Fig. 7 Percentage of PIB scaled by the transmitted power as a function of the diffraction parameter Λ 0 for several distances of the receiver for a collimated beam at the transmitter (Kolmogorov turbulence case only).

Fig. 8
Fig. 8 Percentage of PIB scaled by the transmitted power as a function of the diffraction parameter Λ 0 for several distances of the receiver for a beam focused on the receiver (Kolmogorov turbulence case only).

Fig. 10
Fig. 10 Percentage of POS scaled by the transmitted power as a function of the diffraction parameter Λ 0 for several distances of the receiver for a beam focused on the receiver (Kolmogorov turbulence case only).

Fig. 9
Fig. 9 Percentage of POS scaled by the transmitted power as a function of the diffraction parameter Λ 0 for several distances of the receiver for a collimated beam at the transmitter (Kolmogorov turbulence case only).

Fig. 12
Fig. 12 Percentage of PIF scaled by the transmitted power as a function of the diffraction parameter Λ 0 for several distances of the receiver for a beam focused on the receiver (Kolmogorov turbulence case only).

Fig. 11
Fig. 11 Percentage of PIF scaled by the transmitted power as a function of the diffraction parameter Λ 0 for several distances of the receiver for a collimated beam at the transmitter (Kolmogorov turbulence case only).