## Abstract

Recently, a new and generalized statistical model, called ℳ or Málaga distribution, was proposed to model the irradiance fluctuations of an unbounded optical wavefront (plane and spherical waves) propagating through a turbulent medium under all irradiance fluctuation conditions in homogeneous, isotropic turbulence. Málaga distribution was demonstrated to have the advantage of unifying most of the proposed statistical models derived until now in the bibliography in a closed-form expression providing, in addition, an excellent agreement with published plane wave and spherical wave simulation data over a wide range of turbulence conditions (weak to strong). Now, such a model is completed by including the adverse effect of pointing error losses due to misalignment. In this respect, the well-known effects of aperture size, beam width and jitter variance are taken into account. Accordingly, after presenting the analytical expressions for the combined distribution of scintillation and pointing errors, we derive its centered moments of the overall probability distribution. Finally, we obtain the analytical expressions for the average bit error rate performance for the ℳ distribution affected by pointing errors. Numerical results show the impact of misalignment on link performance.

© 2012 OSA

## 1. Introduction

Atmospheric optical communication (AOC) is receiving considerable attention recently for use in high data rate wireless links. However, even in clear sky conditions, wireless optical links experience fluctuations in both the intensity and the phase of an optical wave propagating through this medium [1] due to time varying inhomogeneities in the refractive index of the atmosphere. The reliability of an optical system operating in such an environment can be deduced from a mathematical model for the probability density function (pdf) of the randomly fading irradiance signal (scintillation). For that reason, one of the goals in studying optical wave propagation through turbulence is the identification of a tractable pdf of the irradiance under all intensity fluctuation regimes.

Over the years, many irradiance pdf models have been proposed with different degrees of success. Perhaps the most successful models are the lognormal [2] and the gamma-gamma [3] ones. The scope of the lognormal model is restricted under weak irradiance fluctuations. Conversely, the gamma-gamma pdf was suggested by Andrews et al. as a reasonable alternative to Beckmann’s pdf [4] because of its much more tractable mathematical model. Nevertheless, both lognormal and gamma-gamma models are particular cases of the new and recently proposed ℳ (Málaga) distribution model [5]. ℳ distribution was validated by comparing its pdf with published simulation data, showing that it unifies in an analytical expression most of the irradiance statistical models proposed in literature by the scientific community for more than four decades. So, for instance: Rice-Nakagami [1], gamma [3], shadowed-Rician [6], K [3, 7], homodyned-K [7], exponential [3] or Gamma-Rician [5] pdfs are considered as particular cases of this Málaga distribution.

In addition to the scintillation effect, misalignment between transmitter and receiver due to building sway causes vibrations of the transmitted beam and, thus, pointing errors that limit the performance of atmospheric free-space optical links. The impact of pointing error (jitter) was widely investigated for intersatellite space-based free-space optical links [8, 9] although without providing closed-form expressions. In a further step, Kiasaleh [10] deduced approximated expressions for the pdf of optical signal intensity in an optical communication channel impaired by motion-induced beam jitter and turbulence, assuming that the beam scintillation is governed by a lognormal pdf and a K distribution. Recently, a new procedure to derive the statistical distribution of pointing errors was proposed in [11]. Results derived in this last work were employed by Sandalidis et al. [12] to study the average bit error rate (ABER) expression for an intensity-modulation/direct detection (IM/DD) atmospheric optical system with on-off keying (OOK), assuming pointing errors, but restricted to strong turbulence.

Now, through this paper we widely generalize the results obtained in [12] for all turbulence regimes, thus completing our works reported in [5, 13] with ℳ distributed turbulence, by including the impact of misalignment fading on the performance of AOC systems. In this sense, we present a statistical AOC channel model considering the joint effects of ℳ distributed turbulence and pointing errors. As detailed in [5], ℳ distribution is valid under any turbulence regime, unifying in a closed-form expression the link performance of most of the irradiance statistical models proposed in literature. In addition, we derive the centered moments of the overall probability distribution and, finally, we provide closed-form expressions for the ABER performance of an IM/DD system over the combined effect of Málaga turbulence and misalignment fading, showing the impairment in performance of pointing errors in atmospheric optical links.

## 2. System model

In this paper, we consider IM/DD channels using on-off keying (OOK) modulation, which is widely employed in practical systems. Laser beams propagate along a horizontal path through a turbulence channel governed by an ℳ distribution in the presence of pointing errors. An additive white Gaussian noise (AWGN) is assumed that includes any shot noise caused by ambient light that may be much stronger than the desired signal as well as front-end thermal noise in the electronics following the photodetector. The receiver integrates the photocurrent signal which is related to the incident optical power by the detector responsivity, *R*, for each bit period. As explained above, the received optical power, *P _{R}*(

*t*), suffers from a fluctuation in signal intensity due to atmospheric turbulence and misalignment, as well as additive noise, and can be modeled, from [11], as:

*P*(

_{T}*t*) being the transmitted optical power, with

*R*being the responsivity, whereas

*N*(

*t*) is the AWGN with variance ${\sigma}_{N}^{2}$, and

*h*(

*t*) is the normalized channel fading coefficient representing the effect of the intensity fluctuations on the transmitted signal due to the combined actions of scintillation and misalignment, in the way explained in [11]. This channel state is considered to be a product of two random components with a third and deterministic factor. The random components that take part in

*h*are geometric spread and pointing errors, denoted by

*h*and atmospheric turbulence, represented by

_{p}*h*, following the same notation employed in [11]. The deterministic component that must be included in

_{a}*h*arises due to atmospheric path loss,

*h*. This atmospheric attenuation,

_{l}*h*, is described by the exponential Beers-Lambert Law and can be expressed in terms of the visibility [11]. However, although weather-induced attenuation can also degrade the performance of AOC systems in the way shown in [14], due to its deterministic nature

_{l}*h*acts as a scaling factor as shown in [15] and for this reason, it is not considered in this manuscript. Thus, through this paper,

_{l}*h*(

*t*) =

*h*(

_{a}*t*) ·

*h*(

_{p}*t*), and each of these two factors are explained in detail through Sections 3 and 4, respectively.

Finally, and as in [11], it is considered that the pointing error is due to building sway, with sub-Hz bandwidth (correlation times of a few seconds), so that the scintillation process must have a smaller correlation time (10–100 ms) as compared with building sways. Hence, the two above impairments are stated to be independent.

## 3. Generalized atmospheric distribution model

In this section, we detail the analytical development for the factor *h _{a}* introduced in Section 2. The starting point is the new statistical model for the irradiance fluctuations of an unbounded optical wavefront propagating through a turbulent medium under all irradiance fluctuation conditions in homogeneous, isotropic turbulence that was presented in [5]. For this reason, a brief description of such a model is achieved through this section in order to improve the overall understanding of this paper. The proposed model, labeled Málaga or simply ℳ, is based on a physical model of the scattering processes that is seen as an extension of the previous work developed by Churnside and Clifford [4]. The main advantage of the Málaga distribution associated to the proposed physical model is that it does reduce to a simple closed-form analytical formula valid from every turbulence regime.

Thus, assume an electromagnetic wave propagating through a turbulent atmosphere with a random refractive index. As the wave passes through this medium, part of the energy is scattered and the form of the irradiance probability distribution is determined by the type of scattering involved. The basis of the ℳ distribution is the conception of an alternative physical model in the generation of small-scale fluctuations. The propagation scheme is illustrated in Fig. 1. As detailed in [5] the observed field at the receiver is supposed to consist of three terms: the first one is the line-of-sight (LOS) contribution, *U _{L}*, the second one is the component which is quasi-forward scattered by the eddies on the propagation axis,
${U}_{S}^{C}$, and coupled to the LOS contribution; whereas the third term,
${U}_{S}^{G}$, is due to energy which is scattered to the receiver by off-axis eddies, this latter contribution being statistically independent from the previous two other terms. The inclusion of the coupled to the LOS scattering component,
${U}_{S}^{C}$, is the main novelty of the model and it can be justified by the high directivity and the narrow beamwidths of laser beams in atmospheric optical communications.

Accordingly, the total observed field when atmospheric turbulence is the only adverse effect considered in *h* can be written as:

*U*and ${U}_{S}^{G}$ are also independent random processes. In Eq. (2),

_{L}*G*is a real variable following a gamma distribution with

*E*[

*G*]=1. It represents the slow fluctuation of the LOS component. Following the notation of [6], the parameter Ω =

*E*[|

*U*|

_{L}^{2}] represents the average power of the LOS term whereas the average power of the total scatter components is denoted by $2{b}_{0}=E\left[{|{U}_{S}^{C}|}^{2}+{|{U}_{S}^{G}|}^{2}\right]$.

*ϕ*and

_{A}*ϕ*are the deterministic phases of the LOS and the coupled-to-LOS scatter terms, respectively. On another note, 0≤

_{B}*ρ*≤1 is the factor expressing the amount of scattering power coupled to the LOS component. Finally,

*U*′

*is a circular Gaussian complex random variable, and*

_{S}*χ*and

*S*are real random variables representing the log-amplitude and phase perturbation of the field induced by the atmospheric turbulence, respectively. A plausible justification for the coupled-to-LOS scattering component, ${U}_{S}^{C}$, is given in [16]. There, it is said that if the turbulent medium is so thin that multiple scattering can be ignored, the multipath delays of the scattered radiation collected by a diffraction-limited receiver will usually be small relative to the signal bandwidth. Then the scattered field will combine coherently with the unscattered field and there will be no-“interfering” signal component of the field, in a similar way as ${U}_{S}^{C}$ combines with

*U*in our proposed model. Of course, when the turbulent medium becomes so thick, then the unscattered component of the field can be neglected.

_{L}Now, and from Eq. (2), the observed irradiance of the proposed propagation model (*h _{a}*, as we have adopted the notation given in [11]) can be expressed as:

*U*′

*) and a Nakagami distribution, $\sqrt{G}$. Then, we can apply the same procedure exposed in [6] consisting in calculating the expectation of the Rayleigh component with respect to the Nakagami distribution and then deriving the pdf of the instantaneous power. Hence, the pdf of*

_{S}*Y*is given by:

*β*≜ (

*E*[

*G*])

^{2}/Var[

*G*] is the amount of fading parameter with Var[·] being the variance operator, where | · | denotes the absolute value. For simplicity, we have denoted $\gamma =E\left[{|{U}_{S}^{G}|}^{2}\right]=2{b}_{0}(1-\rho )$ whereas ${\mathrm{\Omega}}^{\prime}=\mathrm{\Omega}+\rho 2{b}_{0}+2\sqrt{2{b}_{0}\mathrm{\Omega}\rho}\text{cos}({\varphi}_{A}-{\varphi}_{B})$represents the average power from the coherent contributions, as detailed in [5]. Finally,

_{1}

*F*

_{1}(

*a; c; x*) is the Kummer confluent hypergeometric function of the first kind.

Otherwise, the large-scale fluctuations, *X* ≜ exp (2*χ*), is widely accepted to be a lognormal amplitude [4] but, as in [3, 6], this distribution is approximated by a gamma one, this latter with a more favorable analytical structure. For further details about the fact that the gamma distribution closely approximates the lognormal distribution, we suggest to consult references [17, 18]. Then:

*α*is a positive parameter related to the effective number of large-scale cells of the scattering process, as in [3].

Thus to obtain the statistical characterization of the intensity, *h _{a}*, the mixture of

*f*(

_{X}*x*) and

*f*(

_{Y}*y*) must be accomplished. After the corresponding and non-obvious algebraic manipulation, extensively detailed in [5], the pdf of the irradiance is represented by:

*β*being a natural number. In Eq. (7),

*K*(·) is the modified Bessel function of the second kind and order

_{ν}*ν*.

A generalized expression of Eq. (7) was also given in [5] with *β* being a real number:

In Eq. (10), (*β*)* _{k}* represents the Pochhammer symbol.

## 4. Misalignment fading model

As explained in Section 2, the channel fading coefficient, *h*(*t*) is considered through this paper to be a product of two random components: the scintillation, *h _{a}*, detailed in Section 3, and geometric spread and pointing errors, denoted by

*h*. This last factor is now analytically detailed through this section. Thus, in atmospheric optical links, pointing accuracy is an important issue in determining link performance. For this reason, Farid and Hranilovic [11] derived its statistical nature,

_{p}*f*

_{hp}(

*h*), assuming a Gaussian spatial intensity profile of beam waist,

_{p}*w*, at distance

_{z}*z*from the transmitter:

**is the radial vector from the beam center. In addition, a circular detection aperture of radius**

*ρ**a*was considered and a Gaussian beam profile at the receiver,

*I*

_{beam}. Then, the fraction of the collected power at a receiver of radius

*a*in the transverse plane of the incident wave when a pointing error,

**r**, is present can be approximated as the Gaussian form:

*r*= ||

**r**|| is the radial distance. In particular, due to the symmetry of the beam shape and the detector area,

*h*depends only on this magnitude,

_{p}*r*. On another note,

*A*

_{0}= [erf(

*v*)]

^{2}is the fraction of the collected power at

*r*= 0, where $v=\left(\sqrt{\pi}a\right)/\left(\sqrt{2}{w}_{z}\right)$, with erf(·) denoting the error function, whereas the square of the equivalent beam width is given by:

By considering independent identical Gaussian distributions for the elevation and the horizontal displacement (sway), the radial displacement, *r*, at the receiver follows a Rayleigh distribution.

*h*is given by:

_{p}*g*=

*w*/(2

_{zeq}*σ*) is the ratio between the equivalent beam radius at the receiver,

_{s}*w*, and the pointing error displacement standard deviation at the receiver,

_{zeq}*σ*. See [11] for more details.

_{s}## 5. Combined channel statistical model

In this section, we derive a complete statistical model of a stochastic atmospheric optical channel taking into account both turbulence-induced scintillation and misalignment-induced fading. Therefore, the unconditional pdf, *f _{h}*(

*h*), for the channel state,

*h*, is obtained by calculating the mixture of the two distributions presented above in Eqs. (7) or (9) and (15):

*f*

_{h|ha}(

*h*|

*h*) is the conditional probability given a turbulence state,

_{a}*h*, and it is expressed as:

_{a}*f*(

_{ha}*h*) was written in Eqs. (7) and (9) depending on whether parameter

_{a}*β*is a natural number or a real one, respectively. For the first case, we substitute Eq. (7) into the integral expression given by Eq. (16) to finally obtain

*f*(

_{h}*h*):

*A*and

*a*were already defined in Eq. (8). From [19, Eq. (07.34.03.0605.01)], the modified Bessel function of the second kind,

_{k}*K*(·) can be expressed as a special case of the Meijer G function, given by the following relationship: where a definition of the Meijer G function can be found at [19, Eq. (07.34.02.0001.01)]. Taking

_{v}*a*= −

*b*= (

*α*−

*k*)/2 in Eq. (19) and substituting such an equation into Eq. (18), then:

Finally, to solve the last integrate, we can employ Eq. (07.34.21.0085.01) in [19] and, thus, the unconditional pdf, *f _{h}*(

*h*), for the channel state,

*h*, can be written in a closed-form expression:

For the generic case of Eq. (9) with *β* being a real number, the obtention of *f _{h}*(

*h*) is not so evident because it is firstly mandatory to verify if the infinite summation of Eq. (9) and the integrate operator that appears when building

*f*(

_{h}*h*) can be interchanged. To prove that, by normalization, it is possible to express Eq. (9) as a discrete mixture [20] in the form of $f({h}_{a})={\sum}_{k=1}^{\infty}{\omega}_{k}{f}_{k}({h}_{a})$ where ${\omega}_{k}\ge 0\forall k,{\sum}_{k=1}^{\infty}{\omega}_{k}=1$ and with

*f*(·) being a Gamma-Gamma pdf. Since the cumulative distribution function (CDF) can be written as $F({h}_{a})={\sum}_{k=1}^{\infty}{\omega}_{k}{F}_{k}({h}_{a})$, where

_{k}*F*(

_{k}*h*) is the CDF of each

_{a}*f*(·), then we can build a new CDF in the form: ${\varphi}_{n}({h}_{a})={\sum}_{k=1}^{n}{\omega}_{k}{F}_{k}({h}_{a})$. By applying the ratio criterion [21] it is easy to test that

_{k}*F*(

_{k}*h*) converges absolutely so

_{a}*ϕ*(

_{n}*h*) is said to converge weakly to

_{a}*F*(

*h*). This condition and the fact that Eq. (17) is bounded and continuous allow us to apply the Helly-Bray theorem [22]. By proceeding in an analogous manner as in [13], we can conclude that all the conditions needed to apply the Helly-Bray theroem [22] are satisfied. Thus, operating in a similar way to obtain Eq. (21), then

_{a}*f*(

_{h}*h*) is written as an infinite series of Meijer-G functions as

#### 5.1. Moments of the combined ℳ probability distribution

In this subsection, the *k ^{th}* moment of the channel state,

*h*, following a combined ℳ probability distribution,

*f*(

_{h}*h*), as given in Eq. (23) is derived. As we indicated in Section 2, the channel state is considered to be a product of

*h*and

_{p}*h*, where

_{a}*h*represents the random effect of pointing errors whereas

_{p}*h*denotes the turbulence-induced scintillation suffered by the observed irradiance. Furthermore,

_{a}*h*can be decomposed into two other components, as indicated in Eq. (3):

_{a}*h*=

_{a}*YX*. The first component, denoted by

*Y*and following a shadowed-Rician distribution [6], is associated to small-scale fluctuations of the irradiance, i.e., the small-scale contributions to scintillation associated with turbulent cells smaller than either the first Fresnel zone or the transverse spatial coherence radius, whichever is smallest; whereas the second component of

*h*, denoted by

_{a}*X*and following a gamma distribution, is related to to large-scale fluctuations of the irradiance that are generated by turbulent cells larger than that of either the Fresnel zone or the so-called “scattering disk”, whichever is largest.

The *k ^{th}* centered moments of

*h*, denoted by

_{a}*m*(

_{k}*h*), following an ℳ distribution was derived in [5]. Based on assumptions of statistical independence for the underlying random processes,

_{a}*X*and

*Y*, then:

*E*[·] stands for ensemble average. From [23, Eq. (2.23)], the moment of a Nakagami-m pdf is given by:

In addition, the moment of a shadowed-Rician distribution with the restriction of *β* being a natural number was derived in [5]. Then:

Hence, when performing the product of Eq. (25) by Eq. (26), we certainly obtain:

And for the generic case of Eq. (9) with the parameter *β* being a real number, the moment of a shadowed-Rician distribution was given in [23, Eq. (2.69)]:

On the other hand, as discussed in [11] and indicated in Section 2, it is considered that pointing errors and turbulence-induced scintillation are independent processes. As in the case of *h _{a}*, we can obtain the

*k*centered moments of the channel state,

^{th}*h*, as

To calculate *m _{k}* (

*h*) we must solve:

_{p}Inserting Eq. (15) into Eq. (31), then

Hence, for the case of *β* being a natural number, we can obtain the *k ^{th}* centered moments of

*h*by performing the product of Eqs. (27) and (32):

Finally, for the generic case of *β* being a real number, the *k ^{th}* centered moments of

*h*is obtained by the product of Eqs. (29) and (32):

## 6. Average BER

The study of the ABER of the ℳ probability distribution in the absence of pointing errors was derived in [13]. Now, in this section, we complete such a previous work by detailing the procedure to derive the analytical expressions for the ABER of the ℳ probability distribution in the presence of misalignment fading. In this sense, an IM/DD link using OOK signalling technique is considered again. For such links, each bit symbol is transmitted by pulsing the light source either on or off during each bit time. Since in most practical systems, the receiver signal-to-noise ratio (SNR) is limited by shot noise caused by ambient light that may be much stronger than the desired signal and/or by thermal noise in the electronics following the photodetector, then noise can usually be modeled to high accuracy as AWGN that is statistically independent of the desired signal [2]. Because of AWGN, errors may be made at the receiver in determining the actual symbols transmitted. The overall probability of error, *P _{b}*(

*e*) is expressed as the following weighted sum:

*p*

_{0}is the transmission probability of a binary “0” and

*p*

_{1}is the transmission probability of a binary “1”, as indicated in [24]. On another note,

*p*(

*e*|0) and

*p*(

*e*|1) denote the conditional bit error probabilities when the transmitted bit is “0” and “1”, respectively. We assume that either symbol (binary “0” or binary “1”) are equally likely to be transmitted, i.e.,

*p*

_{0}=

*p*

_{1}= 1/2. Considering that

*p*(

*e*|0) =

*p*(

*e*|1), then it is straightforward to show that the conditional BER is expressed, from [13, 24], as where

*P*is the average of optical power, with

*R*being the responsivity and ${\sigma}_{N}^{2}$ being the AWGN variance while

*Q*(·) is the Gaussian Q function which is related to the complementary error function, erfc(·), by $\text{erfc}(X)=2Q\left(\sqrt{2}x\right)$. The ABER,

*P*(

_{b}*e*), can be obtained by averaging

*P*(

_{b}*e*|

*h*) over

*f*(

_{h}*h*): as shown in [12, 13].

Thus, for the particular case of *β* being a natural number, the pdf of the channel state, *f _{h}*(

*h*), is taken from Eq. (21). Now we can express the Gaussian Q function appeared in Eq. (36) in terms of the complementary error function, erfc(

*x*). Next, we can write both the erfc(·) and the modified Bessel integrands as Meijer G functions by employing [19, Eqs. (07.34.03.0619.01), (07.34.03.0605.01) ], respectively, as in [12, 13]. Such relationships are written here for added convenience: in fact, we have already written in Eq. (19) the relationship between a modified Bessel function and a Meijer G function. In the same way, the erfc(·) and the Meijer G function are related by:

Finally, using [19, Eq. (07.34.21.0013.01)], the ABER is derived in a closed-form expression as:

For the generic case of the irradiance pdf represented by Eq. (9) with *β* being a real number, we verify again that the infinite summation and the integrate operator can be interchanged. Operating in a similar way as to obtain Eq. (39), the ABER is written as:

## 7. Numerical results and discussions

Some numerical results are shown in Figs. 2–3, where the IM/DD AOC system is evaluated for its error performance capabilities against average SNR, and in terms of the Rytov variance, ${\sigma}_{R}^{2}$, as in [12, 13]. Effectively, as indicated in [24], atmospheric turbulence is well characterized by such a Rytov variance when using the optical wave models of an infinite plane wave or a spherical wave (point source). The Rytov variance is defined, for a plane wave, by [24, 11]:

where ${C}_{n}^{2}({\text{m}}^{-2/3})$ is the index of refraction structure parameter of the atmosphere,*k*= 2

*π*/

*λ*is the optical wave number,

*λ*(m) is the wavelength, and

*L*(m) is the propagation path length between transmitter and receiver. To obtain typical values of ${\sigma}_{R}^{2}$, we have employed some experimental data of ${C}_{n}^{2}$ provided by University of Waseda, in Japan, on October 15th, 2009, obtained from atmospheric parameters. Thus, near midday, ${C}_{n}^{2}$ parameter oscillated from 2.8 to 11 · 10

^{−14}m

^{−2/3}, and the temperature oscillated from 22 °C to 26.4 °C, and the received power was recorded from −22.52 dBm to −20.48 dBm. At night (near 1 a.m.), ${C}_{n}^{2}$ parameter registered its minimum value (∼ 7 ·10

^{−15}m

^{−2/3}), for a temperature of 15.1 °C and a received power of −22.91 dBm. At sunrise (6.45 a.m.), ${C}_{n}^{2}=1.2\cdot {10}^{-14}\hspace{0.17em}{\text{m}}^{-2/3}$, for a temperature of 15.9 °C and with a received power of −23.47 dBm. The wavelength employed in all cases was 785 nm and the receiver aperture had a diameter of 100 mm. The propagation path length is

*L*= 1 km at a height of 25 m. The optical power transmitted is 11.5 dBm with a responsivity of 0.8 A/W. With such a set of parameters, if we consider ${C}_{n}^{2}$ to 7.2 ·10

^{−15}m

^{−2/3}(night), then ${\sigma}_{R}^{2}\approx 0.32$; secondly, for the case of sunrise when ${C}_{n}^{2}=1.2\hspace{0.17em}{\text{m}}^{-2/3}$, then ${\sigma}_{R}^{2}\approx 0.52$; finally, for a ${C}_{n}^{2}=2.8\cdot {10}^{-14}{\text{m}}^{-2/3}$ (near midday), then ${\sigma}_{R}^{2}\approx 1.2$. This will be the values of Rytov variance employed through this section. As a final remark, all cases presented were obtained by employing the same normalized average optical power (Ω + 2

*b*

_{0}= 1).

Figure 2 shows the ABER for various values of the normalized beamwidth: *w _{z}*/

*a*=10, 20, and 25. Parameters of ℳ turbulence model were fixed to:

*α*= 10 and

*β*= 5 in all cases. Moreover, for each

*w*/

_{z}*a*value, three different behaviors were obtained:

*ρ*= 1, which corresponds to a Gamma-Gamma distribution;

*ρ*= 0.75, and

*ρ*= 0.25. For such values, the associated Rytov variance is ${\sigma}_{R}^{2}=0.32,0.52$ and 1.2, respectively, coinciding with the values selected above. As a reference, the behavior of the system in the absence of misalignment fading is also provided (dotted curves) for the three intensities of turbulence mentioned before [13]. Hence, it is straightforward to check the impairment in performance when pointing errors are considered: for instance, the penalties in SNR with respect to the ideal case where no pointing errors are considered are 17.03, 23.02 and 24.95 optical dB, for

*w*/

_{z}*a*= 10, 20 and 25, respectively, for a bit error rate of 10

^{−6}.

In Fig. 3, different ABERs are displayed for *w _{z}*/

*a*= 20 and a same intensity of turbulence ( ${\sigma}_{R}^{2}=0.36$, ${C}_{n}^{2}=0.83\cdot {10}^{-14}{\text{m}}^{-2/3}$, measured on October 15th at 23:10), a typical value in terrestrial horizontal links. It turns out that the behavior of each curve is different depending on

*ρ*. Thus, when

*ρ*= 1 the overall scattering power travels through ${U}_{S}^{C}$, so it is entirely coupled to the LOS component, minimizing its adverse effect, as displayed in Fig. 3. When

*ρ*is getting lower then the scattering power is distributed both into ${U}_{S}^{C}$ and into ${U}_{S}^{G}$. The presence of a higher power from the ${U}_{S}^{G}$ component has a more harmful effect in the performance of an AOC system because it is not correlated to the LOS contribution, as corroborated in Fig. 3. Obviously, obtained results are consistent with the ones achieved in [13], being the results shown in this paper a 23.03 optical dB displaced version with reference to the ones displayed in Fig. 3 of [13]. For this reason, same conclusions derived there are applicable when pointing errors are considered.

On the other hand, all BER curves tend to a same diversity order, where the diversity order is defined in [25].

Finally, we must recall that Eq. (7) is a particularization of Eq. (9) valid when *β* is restricted to a natural number. Nevertheless, due to the high degree of freedom inherent to the ℳ distribution, demonstrated in the set of results shown in [5], this particularization allows us to reproduce every turbulent scenario, avoiding the necessity of employing the infinite summation included in Eqs. (9) or (40).

## 8. Concluding remarks

In summary, in this paper, two analytical closed-form representations for the ABER performance of an AOC system operating over a generalized turbulence in the presence of pointing errors are derived. Due to the high degree of freedom inherent to the ℳ distribution [5], the particularization given in Eq. (39) let us reproduce every turbulent scenario, avoiding the infinite summation included in Eq. (40).

As discussed in Section 7, the effect of pointing errors on the final performance of an AOC system can be seen as a penalty in SNR with respect to the ideal case where no pointing errors are considered, but with the special feature that such a penalty in SNR is seen as a constant value for every different turbulence regime when the ratio *w _{z}*/

*a*is fixed, and for a given value of bit error rate.

Finally, as can be deduced, the consideration of both factors, turbulence and misalignment, affecting the AOC systems performance, is an essential AOC link design criterion. Through this paper, we have derived a generalized analytical expression to represent the ABER of an AOC system including the presence of pointing errors, so that a designer can obtain in a straightforward manner the ABER related to a particular location affected for both turbulence-induced scintillation and pointing errors due to misalignment between the transmitter and receiver.

## Acknowledgment

This work was supported by the Spanish Ministerio de Ciencia e Innovación, Project TEC2008-06598. The authors thank Prof. Mitsuji Matsumoto from University of Waseda, in Japan, for his kind support and for providing some experimental measurements of the ${C}_{n}^{2}$ parameter needed to calculate realistic values of Rytov variance.

## References and links

**1. **L. C. Andrews and R. L. Phillips, *Laser Beam Propagation Through Random Media* (SPIE, 1998).

**2. **X. Zhu and J. M. Kahn, “Free-space optical communication through atmospheric turbulence channels,” IEEE Trans. Commun. **50**, 1293–1300 (2002). [CrossRef]

**3. **M. A. Al-Habash, L. C. Andrews, and R. L. Phillips, “Mathematical model for the irradiance probability density function of a laser beam propagating through turbulent media,” Opt. Eng. **40**, 1554–1562 (2001). [CrossRef]

**4. **J. H. Churnside and S. F. Clifford, “Log-normal Rician probability-density function of optical scintillations in the turbulent atmosphere,” J. Opt. Soc. Am. A **4**, 1923–1930 (1987). [CrossRef]

**5. **A. Jurado-Navas, J. M. Garrido-Balsells, J. F. Paris, and A. Puerta-Notario, “A unifying statistical model for atmospheric optical scintillation,” in: *Numerical Simulations of Physical and Engineering Processes* (Intech, 2011). [CrossRef]

**6. **A. Abdi, W. C. Lau, M. S. Alouini, and M. A. Kaveh, “A new simple model for land mobile satellite channels: first- and second-order statistics,” IEEE Trans. Wireless Commun. **2**, 519–528 (2003). [CrossRef]

**7. **E. Jakerman, “On the statistics of K-distributed noise,” J. Phys. A **13**, 31–48 (1980). [CrossRef]

**8. **S. Arnon and N. S. Kopeika, “Laser satellite communication network-vibration effect and possible solutions,” Proc. IEEE **85**, 1646–1661 (1997). [CrossRef]

**9. **S. Arnon, “Optimization of urban optical wireless communication systems,” IEEE Trans. Wireless Commun. **2**, 626–629 (2003). [CrossRef]

**10. **K. Kiasaleh, “On the probability density function of signal intensity in free-space optical communications systems impaired by pointing jitter and turbulence,” Opt. Eng. **33**, 3748–3757 (1994). [CrossRef]

**11. **A. A. Farid and S. Hranilovic, “Outage capacity optimization for free-space optical links with pointing errors,” J. Lightwave Technol. **25**, 1702–1710 (2007). [CrossRef]

**12. **H. G. Sandalidis, T. A. Tsiftsis, G. K. Karagiannidis, and M. Uysal, “BER performance of FSO links over strong atmospheric turbulence channels with pointing errors,” IEEE Commun. Lett. **12**, 44–46 (2008). [CrossRef]

**13. **A. Jurado-Navas, J. M. Garrido-Balsells, J. F. Paris, M. Castillo-Vázquez, and A. Puerta-Notario, “General analytical expressions for the bit error rate of atmospheric optical communication systems,” Opt. Lett. **36**, 4095–4097 (2011). [CrossRef] [PubMed]

**14. **M. Al Naboulsi and H. Sizun, “Fog attenuation prediction for optical and infrared waves,” Opt. Eng. **43**, 319–329 (2004). [CrossRef]

**15. **ITU-R Report F.2106-1 “Fixed service applications using free-space optical links,” Nov. 2010.

**16. **R. S. Kennedy, “Communication through optical scattering channels: an introduction,” Proc. IEEE **58**, 1651–1665 (1970). [CrossRef]

**17. **T. Rappaport, *Wireless Communications: Principles and Practice*, 2nd ed. (Prentice Hall, 2001.).

**18. **J. R. Clark and S. Karp, “Approximations for lognormally fading optical signals,” Proc. IEEE **58**, 1964–1965 (1970). [CrossRef]

**19. **Wolfram, http://functions.wolfram.com/

**20. **C. A. Charalambides, *Combinatorial Methods in Discrete Distributions* (John Wiley & Sons, 2005) [CrossRef]

**21. **L. C. Andrews, *Special Functions of Mathematics for Engineers*, 2nd ed. (SPIE, 1998).

**22. **P. Billingsley, *Convergence of Probability Measures*, 2nd ed. (John Wiley & Sons, 2005)

**23. **M. K. Simon and M. S. Alouini, *Digital Communications over Fading Channels*, 2nd ed. (Wiley-Interscience, 2005).

**24. **L. C. Andrews, R. L. Phillips, and C. Y. Hopen, *Laser Beam Scintillation with Applications* (SPIE, 2001). [CrossRef]

**25. **L. Zheng and D. N. C. Tse, “Diversity and multiplexing: A fundamental tradeoff in multiple antenna channels,” IEEE Trans. Inf. Theory **49**, 1073–1096 (2003). [CrossRef]