Centroid Error Analysis of Beacon Tracking under Atmospheric Turbulence for Optical Communication Links

: Optical satellite communication has received considerable attention as a promising al-ternative to radio frequency communication because of its potential advantages including higher data rates and license free spectrum. Many studies have conducted performance analyses of optical communication channels, but few have investigated beacon tracking channels under atmospheric turbulence. The centroid accuracy of beacon tracking channels is limited by not only noise sources, but also a ﬁnite delay time, which also ﬂuctuates due to atmospheric turbulence. Consequently, the centroid error is an important ﬁgure of merit when evaluating the performance of a beacon tracking system. In this study, the closed-form expressions were derived for average centroid error and fade probability, based on received photoelectron counts depending on exposure time, taking into account the log-normal tracking channels. We analyzed the angular positioning performance of beacon tracking detectors onboard small satellites in the presence of atmospheric turbulence, in terms of centroid error and fade probability. We found that an optimal exposure time exists, which minimizes the centroid error, and that fade probability is inversely proportional to the exposure time. These are signiﬁcant properties to consider in the design of beacon tracking systems.


Introduction
Optical satellite communication has received significant attention as a promising advance in satellite communication, with the potential advantages of higher data rate, license free spectrum, better security, smaller size, lower mass, and lower power consumption [1][2][3][4][5]. Because of the extremely narrow optical beamwidth, a very precise beam steering system is required to ensure a reliable and stable optical communication link, since this is usually beyond the satellite's body pointing capability. For this purpose, pointing, acquisition, and tracking (PAT) technology, considered one of the key technologies in optical satellite communication, must be employed to maintain an optical line-of-sight (LOS) between the transmitter and receiver, even in the face of pointing jitter caused by satellite platform vibrations and atmospheric turbulence. The PAT system is normally composed of three components that deal with different requirements: a coarse pointing assembly (CPA), a fine pointing assembly (FPA), and a point-ahead assembly (PAA).
The CPA provides wide angle steering, but is a low bandwidth mechanism for initial pointing and spatial acquisition. It employs wide dynamic range gimbals and wide fieldof-view (FoV) image sensors, based on charge coupled device (CCD) or complementary metal-oxide-semiconductor (CMOS). The CCD and CMOS sensors allow a low frame rate, on the order of a dozen frames per second, even though they are capable of covering the uncertainty region of the receiver. With the use of a windowed readout mode for the array containing the beacon laser spots, the image sensors can also provide a sufficiently high frame rate and good positioning accuracy required for fine tracking [6][7][8].
array containing the beacon laser spots, the image sensors can also provide a sufficiently high frame rate and good positioning accuracy required for fine tracking [6][7][8].
The FPA detector provides not only more precise and faster LOS alignment, but also a narrower steering angle. To accomplish this, it employs a fast steering mirror and a position-sensitive detector such as quadrant PIN photodiode (QPD) or quadrant avalanche photodiode (QAPD). The position-sensitive detectors have several advantages over classical image sensors including better positioning accuracy, a faster response time, and very simple signal conditioning circuits [9]. The CPA have to be able to control several-centimeter-class telescopes within a pointing accuracy of sub-micro radians for planetary missions, or several micro radians for Earth orbiting missions, and its bandwidth should also be larger than 1 kHz to accommodate the pointing jitters caused by the satellite platform vibration and atmospheric turbulence [4,8,[10][11][12].
Because of the relative velocity of the two terminals and the finite speed of light, an angular offset exists between the transmitting and received beams, the so-called pointahead angle. The magnitude of angular offset is achievable with the help of known ephemeris data, up to 70 micro radians for inter-satellite links and as much as 500 micro radians for interplanetary links [13][14][15]. The PAA compensates the point-ahead angle by controlling the fast steering mirror if the angle is comparable to or larger than the transmitting beamwidth. Figure 1 shows the block diagram of the optical communication terminal including the PAT system with two types of actuators: gimbal and fast steering mirror (FSM). When a beacon beam propagates through the atmosphere, it can be adversely affected by atmospheric turbulence including atmospheric scintillation such as intensity and phase fluctuations, angle of arrival fluctuations, beam spot wandering, and spatial beam spreading [16][17][18]. The beacon image follows the point spread function of the imaging optics since the uplink beacon can be regarded as a point source rather than an extended object [19], which can be distorted and blurred by atmospheric turbulence. Therefore, the centroid error of beacon images on a beacon detector depends on atmospheric turbulence as well as the exposure time, which plays a great role in the signal-to-noise ratio (SNR).
The limited accuracy of beacon detectors has been analyzed in terms of the noise equivalent angle, known as the random centroid error [20][21][22], particularly in the field of adaptive optics to evaluate the Shack-Hartmann wavefront sensors [23][24][25][26][27]. However, the intensity fluctuations resulting from atmospheric turbulence were not taken into consideration in these reports to precisely assess the centroid error. Over the years, many mathematical models have been proposed to describe the turbulence-induced scintillation in When a beacon beam propagates through the atmosphere, it can be adversely affected by atmospheric turbulence including atmospheric scintillation such as intensity and phase fluctuations, angle of arrival fluctuations, beam spot wandering, and spatial beam spreading [16][17][18]. The beacon image follows the point spread function of the imaging optics since the uplink beacon can be regarded as a point source rather than an extended object [19], which can be distorted and blurred by atmospheric turbulence. Therefore, the centroid error of beacon images on a beacon detector depends on atmospheric turbulence as well as the exposure time, which plays a great role in the signal-to-noise ratio (SNR).
The limited accuracy of beacon detectors has been analyzed in terms of the noise equivalent angle, known as the random centroid error [20][21][22], particularly in the field of adaptive optics to evaluate the Shack-Hartmann wavefront sensors [23][24][25][26][27]. However, the intensity fluctuations resulting from atmospheric turbulence were not taken into consideration in these reports to precisely assess the centroid error. Over the years, many mathematical models have been proposed to describe the turbulence-induced scintillation in terms of probability density functions including the log-normal, K-channel, negative exponential, and gamma-gamma models [28][29][30][31][32]. Specifically, the log-normal and gammagamma models are widely used due to their simplicity and good agreement between the theoretical and experimental data. Consequently, they have been found to be suitable for weak-to-moderate and moderate-to-strong turbulence conditions, respectively [33].
An accurate centroid estimation is critical to the beacon tracking system to point the optical beam toward the partner terminal, so that sufficient optical power is received to meet the data rate and bit error rate performance requirements. This makes it essential, when evaluating the pointing capability of a PAT system, to analyze the centroid error intrinsic to noise sources and the finite delay time under atmospheric turbulence.
Korea has a plan to conduct an optical downlink experiment employing a CubeSat that requires not only the use available commercial off-the-shelf (COTS) components for a reduction in development costs, but also accurate pointing and tracking capabilities under the size, weight, and power constraints for a reliable communication link. Thus, the beacon tracking system onboard a CubeSat should be designed through a trade-off between two requirements for COTS components as well as pointing and tracking capabilities over fade channels.
In this study, the SNR formula of position-sensitive detectors was given in terms of the received photoelectron count depending heavily on the exposure time, which also includes shot noise, background noise, dark current, and thermal noise. Taking into account the log-normal tracking channels, the closed-form models were derived for two metrics, average centroid error and fade probability in the function of the received photoelectron count. For practical applications, the received photon count was estimated by considering turbulence-induced beam spreading and atmospheric attenuation due to atmospheric constituents and cirrus clouds.
We analyzed the angular positioning performance of an uplink beacon system for two metrics in terms of exposure time and atmospheric turbulence strength using CMOS and QAPD detectors. We found from numerical simulations that an optimal exposure time exists, which minimizes the centroid error resulting from noise sources and the exposure time latency. However, fade probabilities decreased rapidly as the exposure time was extended and atmospheric turbulence became strong. These characteristics should be taken into account in a trade-off of the design of the beacon tracking system requiring high bandwidth and high positioning accuracy.

Link Equation
The link equation is widely used in communications to predict the power received at an optical detector, because it enables the analysis of the trade-offs between system parameters and channel performance. The intensity of a beacon beam is non-uniformly distributed in the spatial domain along the direction of beam propagation. Taking into account the non-uniform irradiance distribution and beam spreading due to atmospheric turbulence, the link equation is given by [18,34] where P T and P R are the transmitting and received optical powers, respectively; W e f f is the effective spot size including the beam divergence and spreading; z is the axial distance from the beam waist; η T is the transmitter optical efficiency; η R is the receiver optical efficiency including the optical filter transmittance; A R is the receiving area; d is the radial distance from the axis of beam propagation (i.e., d = zθ e for a small pointing error of θ e ); and T A and T C are the atmospheric attenuations produced by atmospheric constituents and cirrus clouds, respectively. The last term indicates the pointing loss due to LOS misalignment between the transmitter and the receiver.
The beacon beam suffers from power losses resulting from absorption and scattering by molecules and aerosols suspended in the atmosphere as well as cirrus clouds. The atmospheric attenuations resulting from the constituents and cirrus clouds are expressed as [35] where σ(λ, V, 0) and V are the atmospheric attenuation coefficient and visibility at sea level, respectively, and λ is the beacon wavelength. h SH = 1.2 km is the scale height; θ Z is the zenith angle of the satellite; h s is the height of the optical ground station; and C T is the cirrus cloud thickness. It is worth noting that in the low elevation angle regimes, atmospheric attenuation by cirrus clouds is dominant compared to attenuation by atmospheric constituents.

Turbulence-Induced Beam Spreading
When the beacon beam propagates through the atmosphere, it undergoes beam spreading due to small-scale turbulent eddies, in addition to the far-field angular spreading coming from the divergence angle. This beam spreading makes the spot size larger, resulting in a reduction in beam irradiance at the detector. Assuming a Gaussian beam, the effective spot size is approximated in the uplink as follows [32]: where T = 4.35k 7/6 (ΛH) 5/6 sec 11/6 (θ z ) Here, T is the additional spreading caused by atmospheric turbulence; k is the optical wavenumber; W 0 is the beam waist (i.e., beam size at z = 0); θ D is the far-field divergence half-angle of the beacon beam; Λ is the Fresnel ratio defined as Λ = 2z/kW 2 (z); H is the difference in altitude between the receiver and the optical ground station at h s ; and C 2 n (h) is the refractive-index structure parameter, which characterizes the strength of atmospheric turbulence.
Many models have been proposed to predict the behavior of the refractive-index structure parameter, which is a quantitative measure of atmospheric turbulence. The Hufnagel-Valley model is widely used since it not only allows for easy variation of the day-time and night-time profile by adjusting various site parameters including wind speed, isoplanatic angle, and altitude, but is also best suited for ground-to-satellite uplink [4]. The model is given as where v is the mean square value of the wind speed in m/s, and A 0 is the nominal value of the refractive-index structure parameter at the ground in m −2/3 . It is noteworthy that the turbulence strength is dominated by A 0 near the ground, but by v at high altitude (i.e., h > 5 km) [36]. The value of A 0 is typically in the range from 10 −17 m −2/3 for weak turbulence to 10 −12 m −2/3 for strong turbulence, with a typical average value being 10 −15 m −2/3 [37].

Log-Normal Tracking Channel
The centroid accuracy of detectors including image sensors or position-sensitive detectors can be significantly degraded by beacon intensity fluctuations due to atmospheric turbulence over the tracking channel. Over the years, a number of statistical models have been introduced to describe intensity fluctuations of a propagated beam in terms of probability density function. The log-normal model is restricted to weak turbulence conditions with scintillation indexes of less than 1.0. Taking into account the exposure time, the average number of received beacon photoelectrons can be calculated as where η Q is the quantum efficiency of the detector; T s is the exposure time in seconds; and and c are the Plank's constant and light speed, respectively. Without loss of generality, K s can also be considered a log-normal distributed random variable instead of the received intensity. Therefore, the probability density function of the log-normal model is defined as [31] f (K s ) = 1 where Here, m k and σ 2 k denote the mean and variance of ln(K s ), respectively, and σ 2 si is the scintillation index to describe the intensity fluctuation. The scintillation index is expressed by various formulas depending on turbulence strength, wave types, and wave propagation path such as uplink, downlink, and horizontal path. In the case of both weak turbulence condition and uplink, σ 2 si can be written for a spherical wave as [12]:

Signal-to-Noise Ratio
The centroid error can be estimated based on the image size as well as SNR, which is defined as the ratio of the signal photoelectron counts to the noise fluctuation intensity during the integration time [27]. The centroid accuracy is proportional to SNR, but various noise sources during the integration time may lead to significant degradation of the centroid accuracy. In the quantitative signal analysis of beacon images on the beacon tracking detector, the centroid algorithm is used to compute the beacon location within the centroid window, which should be carefully selected so as not to deteriorate the centroid accuracy. Assuming the background noise is uniformly distributed over each pixel, the CMOS sensor's SNR corresponding to the exposure time is given by where N p is the number of pixels involved in the centroid window; N d is the dark current in electrons per second per pixel; N r is the read noise in root-mean-square electrons per pixel; and N b is the background noise in photoelectrons per second per pixel and can be calculated from η Q λP b / cN FoV , where N FoV is the number of pixels corresponding to the FPA detector's FoV, and P b is the received power from background radiance. Note that Remote Sens. 2021, 13,1931 6 of 14 the SNR equation for the CCD detector is exactly the same as for the CMOS detector. The background power is given by [4] where L λ is the background spectral radiance; Ω FoV is the detector FoV in steradians; and ∆λ f ilter is the optical bandwidth. A quadrant position-sensitive detector consists of four individual yet identical photodiodes, and each quadrant generates an output photocurrent that is proportional to the amount of incident optical power, depending on the beam spot position and shape on the photodiode surface. Because of these properties, it provides continuous positioning data, higher positioning accuracy, and high-speed response. Unlike QPD detectors, QAPD detectors are used to detect the position of a faint beacon spot due to its high sensitivity. Assuming the responsivity and dark current are the same in each quadrant, the QAPD detector's SNR corresponding to the exposure time is derived as where F is the excess noise factor, and G is the gain. N d is the APD dark current consisting of two components (i.e., surface dark current I ds and bulk dark current I db ), which can be calculated from the equation FG 2 I db + I ds /q, because the surface dark current does not trigger an avalanche process within the substrate but the bulk dark current does [38]. σ 2 n is the thermal noise, defined as 2k B T/R L , where k B is the Boltzmann constant; T is the absolute temperature; and R L is the load resistance of the detector circuit. Note that we can also get the SNR equation of the QPD detector from Equation (15) by taking F = 1 and G = 1.

Centroid Error under Atmospheric Turbulence
Regardless of sensor type (i.e., image sensors or position-sensitive detectors), the centroid error is dominated by not only the angular size of the beacon image, but also SNR. Assuming there is no intensity fluctuation due to atmospheric turbulence, the instantaneous centroid error of a CMOS detector is expressed as [21,22] where iFoV is the FoV of a single pixel; N is the truncated half width of a centroid window with the relation N p = (2N + 1) 2 ; and K nC is defined as N p N b T s + N d T s + N 2 r . Assuming diffraction-limited optics, the instantaneous centroid error of a QAPD detector is, in the absence of intensity fluctuation, given by [20,39] where D is the diameter of the receiving aperture, and K nQ is defined as 4 FG 2 N b T s + N d T s + σ 2 n T s /q 2 /G 2 . If the optics is not diffraction-limited, Equation (17) should be modified by replacing the factor (2.44λ/D) with the actual spot size in angular diameter [40].
However, the beacon beam does experience intensity fluctuation during atmospheric propagation, so the centroid error also fluctuates due to the varying SNR. Thus, the Remote Sens. 2021, 13,1931 7 of 14 probability density function should take into account the intensity fluctuation. With some mathematical effort, the average centroid error of the CMOS detector can be written as where w i and x i denote the weight factors and the zero points of the nth-order Hermite polynomial, respectively, whose values are well tabulated in [37]. The order of Hermite polynomial should be chosen depending on the desired accuracy, which was taken as 20 in the numerical simulation. Following a procedure similar to Equation (18), the average centroid error of the QAPD detector can also be written as There is a finite delay time between the beacon sensing and the pointing of the optical beam toward the partner terminal. Due to the temporal evolution of atmospheric turbulence, the finite delay time also produces a centroid error, which involves the detector exposure time, the detector readout time, and the control-loop response time [41]. Taking only the exposure time latency into consideration in this study, the centroid error is given by [23,24] where r 0 and t 0 are the atmospheric coherence length called the Fried parameter, and the coherence time, respectively. The coherence time is reported to be on the order of milliseconds. For a spherical wave, they are written by [12,42] where v is the average wind speed over different atmospheric layers. So, the total centroid error, due to noise sources and exposure time latency, is written as

Fade Probability
Fade probability is a significant metric when evaluating the reliability of a tracking channel. It is defined as the probability when the SNR falls below a certain threshold over the tracking channels. It is related to the detector's sensitivity limit, and can also be described in terms of received photon counts. Thus, the fade probability of a CMOS detector, over the log-normal tracking channels, is expressed as where γ thCMOS and K thCMOS are the thresholds in terms of SNR and received photoelectron counts, respectively, and er f c(·) is the complementary error function.
Applying the same conditions as Equation (24), the fade probability of a QAPD detector is also expressed as where K thQAPD is the threshold of received photoelectron counts. Given that the SNR threshold is 1 (i.e., γ thCMOS = 1 and γ thQAPD = 1), we can compute two thresholds for the CMOS and QAPD detectors from Equations (13) and (15):

Numerical Results
In this section, we present the numerical results of an analysis of the angular positioning performance of two beacon detectors (CMOS and QAPD) over log-normal tracking channels in terms of centroid error and fade probability. The numerical simulation was performed by considering two detectors onboard a small satellite placed at an altitude (H + h s ) of 500 km and a zenith angle θ z of 70 • , leading to a distance z of 1193 km from an optical ground station located at a height h s of 100 m. The station transmits a beacon laser to the satellite with a wavelength λ of 850 nm, the transmitting power P T of 10 W [19,43], the transmitting optical efficiency η T of 0.95, the divergence θ D of 2.5 mrad in half-angle, and the pointing error θ e of 1.5 mrad, leading to a received power loss of 3 dB. Taking a clear sky into account, the atmospheric attenuation coefficient and cirrus cloud thickness were σ(λ, V, 0) = 0.25 and C T = 1100 m, respectively. The spectral radiance of the background is assumed to be L λ = 71 W/m 2 · ster · µm at 850 nm.
To compare the angular positioning performance of the two detectors in the presence of atmospheric turbulence, the two detectors were considered to have the same optical configurations including receiving aperture and optical efficiency, so that the number of received photons on the detectors were of the same order. We selected the MT9P031 model as the CMOS detector [43], and the S4402-SPL model as the QAPD detector [44]. Considering the optical communication terminal onboard a CubeSat, the common parameters related to the receiver are summarized in Table 1. Taking the beam spreading due to atmospheric turbulence into account, the average count of received photoelectrons can be calculated from the link equation, which is dominated by the detector exposure time. It is remarkable that the turbulence-induced beam spreading angle is so small (i.e., on the order of micro radians) that it is negligible compared to the beam divergence of 5 mrad in this simulation, but it has a great effect on the link equation for beam divergence with the order of micro radians. Figure 2 shows the received photoelectron counts and SNRs under atmospheric turbulence with A 0 = 1.7 × 10 −14 m −2/3 and v = 21 m/s in terms of the exposure time. The pixel number N p used for the centroid window was assumed to be 9. As expected from Equation (8), K s is proportional to the exposure time, and the K s of the QAPD detector Remote Sens. 2021, 13, 1931 9 of 14 is larger than that of the CMOS detector due to the high quantum efficiency. However, the SNR of the QAPD detector was higher than that of the CMOS detector in the short exposure time regime, and vice versa in the long exposure time regime. The results come from the fact that the noise counts (i.e., background, dark current, and thermal noises) are significantly dominated by the exposure time for the QAPD detector, while the readout noise is independent of the exposure time for the CMOS detector excluding background and dark current noises.
spreading angle is so small (i.e., on the order of micro radians) that it is negligible compared to the beam divergence of 5 mrad in this simulation, but it has a great effect on the link equation for beam divergence with the order of micro radians. Figure 2 shows the received photoelectron counts and SNRs under atmospheric turbulence with is larger than that of the CMOS detector due to the high quantum efficiency. However, the SNR of the QAPD detector was higher than that of the CMOS detector in the short exposure time regime, and vice versa in the long exposure time regime. The results come from the fact that the noise counts (i.e., background, dark current, and thermal noises) are significantly dominated by the exposure time for the QAPD detector, while the readout noise is independent of the exposure time for the CMOS detector excluding background and dark current noises. Given the optical configuration, the centroid error is influenced by not only noise sources, but also the finite delay time causing the temporal evolution of atmospheric turbulence. Note that the exposure time latency that causes the finite delay time is taken into consideration in this simulation. Figure 3 shows the average centroid errors due to noise Given the optical configuration, the centroid error is influenced by not only noise sources, but also the finite delay time causing the temporal evolution of atmospheric turbulence. Note that the exposure time latency that causes the finite delay time is taken into consideration in this simulation. Figure 3 shows the average centroid errors due to noise sources over log-normal tracking channels, and Figure 4 shows the centroid errors due to the exposure time latency, for two atmospheric turbulence conditions (i.e., A 0 = 1.0 × 10 −16 m −2/3 and v = 5 m/s for the weak turbulence condition, A 0 = 1.7 × 10 −14 m −2/3 and v = 21 m/s for the moderate turbulence condition). Applying the average wind speed v of 15 m/s, these two turbulence conditions provide values of the scintillation index, coherence length, and time (i.e., σ 2 si = 0.28, r 0 = 14.1 cm and t 0 = 2.9 ms for the weak turbulence condition, and σ 2 si = 0.79, r 0 = 7.4 cm and t 0 = 1.6 ms for the moderate turbulence condition).
The average centroid error over log-normal tracking channels was obviously smaller under the weak turbulence condition than the moderate condition. Compared to the CMOS detector, in terms of the average centroid error, the QAPD detector showed better performance regardless of the exposure time, even though its SNR was lower than the CMOS detector in the long exposure time regime, as shown in Figure 2.
As shown in Equation (20), the centroid error resulting from the exposure time latency is totally unrelated to the detector types, but depends on exposure time, coherence length, and diffraction-limited optics (λ/D). The centroid error resulting from the exposure time latency increases significantly as the exposure time is extended, and atmospheric turbulence becomes stronger because the coherence length is inversely proportional to the strength of atmospheric turbulence.
The total centroid errors are shown in Figures 5 and 6 for weak and moderate atmospheric turbulence conditions, respectively. The average centroid error due to noise sources was inversely proportional to exposure time, but the centroid error resulting from exposure time latency increased exponentially in accordance with exposure time. Consequently, an optimal exposure time exists, which produces the total smallest centroid error under the designed optical configuration and atmospheric turbulence. This optimal exposure time is one of the critical design factors considered to meet the pointing requirements of the beacon tracking system. ms 6 . 1 0 = t for the moderate turbulence condition).
The average centroid error over log-normal tracking channels was obviously smaller under the weak turbulence condition than the moderate condition. Compared to the CMOS detector, in terms of the average centroid error, the QAPD detector showed better performance regardless of the exposure time, even though its SNR was lower than the CMOS detector in the long exposure time regime, as shown in Figure 2.
As shown in Equation (20), the centroid error resulting from the exposure time latency is totally unrelated to the detector types, but depends on exposure time, coherence length, and diffraction-limited optics ( D λ ). The centroid error resulting from the exposure time latency increases significantly as the exposure time is extended, and atmospheric turbulence becomes stronger because the coherence length is inversely proportional to the strength of atmospheric turbulence.   The total centroid errors are shown in Figures 5 and 6 for weak and moderate atmospheric turbulence conditions, respectively. The average centroid error due to noise sources was inversely proportional to exposure time, but the centroid error resulting from exposure time latency increased exponentially in accordance with exposure time. Consequently, an optimal exposure time exists, which produces the total smallest centroid error under the designed optical configuration and atmospheric turbulence. This optimal exposure time is one of the critical design factors considered to meet the pointing requirements of the beacon tracking system.
Under weak turbulence conditions, the total centroid accuracy was 1.8 μrad at an exposure time of 0.5 ms for the CMOS detector, and 1.2 μrad at the exposure time of 0.3 ms for the QAPD detector. In moderate turbulence conditions, the total centroid accuracy was 2.7 μrad at the exposure time of 0.4 ms for the CMOS detector, and 1.8 μrad at the Under weak turbulence conditions, the total centroid accuracy was 1.8 µrad at an exposure time of 0.5 ms for the CMOS detector, and 1.2 µrad at the exposure time of 0.3 ms for the QAPD detector. In moderate turbulence conditions, the total centroid accuracy was 2.7 µrad at the exposure time of 0.4 ms for the CMOS detector, and 1.8 µrad at the exposure time of 0.2 ms for the QAPD detector. The QAPD detector provided better performance than the CMOS detector under atmospheric turbulence in terms of the total centroid accuracy, especially in the short exposure time domain. This means that the QAPD detector is more suitable for the PAT system, which requires accurate pointing capabilities and higher bandwidth for closed-loop beacon tracking.
of the beacon tracking system.
Under weak turbulence conditions, the total centroid accuracy was 1.8 μrad at an exposure time of 0.5 ms for the CMOS detector, and 1.2 μrad at the exposure time of 0.3 ms for the QAPD detector. In moderate turbulence conditions, the total centroid accuracy was 2.7 μrad at the exposure time of 0.4 ms for the CMOS detector, and 1.8 μrad at the exposure time of 0.2 ms for the QAPD detector. The QAPD detector provided better performance than the CMOS detector under atmospheric turbulence in terms of the total centroid accuracy, especially in the short exposure time domain. This means that the QAPD detector is more suitable for the PAT system, which requires accurate pointing capabilities and higher bandwidth for closed-loop beacon tracking.  The fade probabilities for the two turbulence conditions are shown in Figure 7, defined as the probability when the SNR falls below 1 over log-normal tracking channels due to the intensity fluctuations produced by atmospheric scintillation. In other words, they mean probabilities when the received photoelectron counts are less than the threshold values of thCMOS K and thQAPD K for the CMOS and QAPD detectors, respectively.
These threshold values increase as the exposure time becomes longer, since noise count is proportional to exposure time. at the exposure time of 10 ms, respectively. The big difference between the threshold values of two detectors at the exposure time of 10 ms occurred because the noise sources of the QAPD detector were dominated more by the exposure time than the CMOS detector. This means that the CMOS detector will provide a better performance than the QAPD detector at long exposure times, in terms of fade probability. As the exposure time increases, the fade probability decreases, since the increase in received photoelectron counts is higher than the threshold value increment, which has to be taken The fade probabilities for the two turbulence conditions are shown in Figure 7, defined as the probability when the SNR falls below 1 over log-normal tracking channels due to the intensity fluctuations produced by atmospheric scintillation. In other words, they mean probabilities when the received photoelectron counts are less than the threshold values of K thCMOS and K thQAPD for the CMOS and QAPD detectors, respectively. These threshold values increase as the exposure time becomes longer, since noise count is proportional to exposure time. The threshold values for fade probabilities were K thCMOS = 65.6 and K thQAPD = 226.4 at the exposure time of 0.1 ms, and K thCMOS = 141.8 and K thQAPD = 2262.2 at the exposure time of 10 ms, respectively. The big difference between the threshold values of two detectors at the exposure time of 10 ms occurred because the noise sources of the QAPD detector were dominated more by the exposure time than the CMOS detector. This means that the CMOS detector will provide a better performance than the QAPD detector at long exposure times, in terms of fade probability. As the exposure time increases, the fade probability decreases, since the increase in received photoelectron counts is higher than the threshold value increment, which has to be taken into consideration in the design of the beacon tracking system.

Conclusions
An accurate centroid estimation of beacon images is required to ensure a beacon tracking system will precisely point an optical beam toward a partner terminal, and provide a reliable optical communication link even in the presence of atmospheric turbulence. However, centroid accuracy is limited by finite delay time as well as noise sources including detector and background noises. Centroid error is also influenced by intensity fluctuations resulting from atmospheric turbulence, which can terminate the beacon tracking link.
In this study, we analyzed the performance of two beacon tracking detectors onboard small satellites under atmospheric turbulence in terms of two metrics: centroid error and fade probability. Taking into account the log-normal tracking channels, the closed-form expressions were derived for not only the average centroid error, but also fade probability. We investigated two metrics for CMOS and QAPD detectors through numerical simulations by considering the windowed readout mode of the CMOS detector for high bandwidth, and confining finite delay time to the exposure time latency. As the exposure time became longer, the average centroid error due to noise sources decreased, but the centroid error caused by exposure time latency increased exponentially.
Atmospheric turbulence degrades the performance of the beacon tracking system, which is a crucial factor in causing the large centroid error in the regimes of long exposure time rather than a short one (i.e., > l ms) due to the angle of arrival fluctuations. In addition, it makes a big difference in fade probabilities between weak and moderate turbulence conditions, on 3 or 4 orders at the exposure time of 0.1 ms, and the longer the exposure time, the greater the difference.
An optimal exposure time exists that produces the smallest total centroid error in the presence of atmospheric turbulence. It ranges from 0.2 ms to 0.5 ms, depending on the detector type and the strength of atmospheric turbulence, and where the total centroid error is to the order of micro radians. The QAPD detector showed better performance for total centroid error than the CMOS detector, while the CMOS detector provided superior performance in terms of the fade probability. Thus, two metrics have to be taken into account in the PAT system design to improve beacon tracking accuracy and to ensure reliable beacon tracking link.

Conclusions
An accurate centroid estimation of beacon images is required to ensure a beacon tracking system will precisely point an optical beam toward a partner terminal, and provide a reliable optical communication link even in the presence of atmospheric turbulence. However, centroid accuracy is limited by finite delay time as well as noise sources including detector and background noises. Centroid error is also influenced by intensity fluctuations resulting from atmospheric turbulence, which can terminate the beacon tracking link.
In this study, we analyzed the performance of two beacon tracking detectors onboard small satellites under atmospheric turbulence in terms of two metrics: centroid error and fade probability. Taking into account the log-normal tracking channels, the closed-form expressions were derived for not only the average centroid error, but also fade probability. We investigated two metrics for CMOS and QAPD detectors through numerical simulations by considering the windowed readout mode of the CMOS detector for high bandwidth, and confining finite delay time to the exposure time latency. As the exposure time became longer, the average centroid error due to noise sources decreased, but the centroid error caused by exposure time latency increased exponentially.
Atmospheric turbulence degrades the performance of the beacon tracking system, which is a crucial factor in causing the large centroid error in the regimes of long exposure time rather than a short one (i.e., >1 ms) due to the angle of arrival fluctuations. In addition, it makes a big difference in fade probabilities between weak and moderate turbulence conditions, on 3 or 4 orders at the exposure time of 0.1 ms, and the longer the exposure time, the greater the difference.
An optimal exposure time exists that produces the smallest total centroid error in the presence of atmospheric turbulence. It ranges from 0.2 ms to 0.5 ms, depending on the detector type and the strength of atmospheric turbulence, and where the total centroid error is to the order of micro radians. The QAPD detector showed better performance for total centroid error than the CMOS detector, while the CMOS detector provided superior performance in terms of the fade probability. Thus, two metrics have to be taken into account in the PAT system design to improve beacon tracking accuracy and to ensure reliable beacon tracking link.