Propagation Characteristics of Hermite–Gaussian Beam under Pointing Error in Free Space

: Hermite–Gaussian (HG) beams have signiﬁcant potential to improve the capacity of free-space optical communication (FSOC). The inﬂuence of pointing error on the propagation characteristics of an HG beam cannot be ignored in the FSOC system. Although the average irradiance of the HG beam under a small pointing error from the FSOC tracking mechanism has been investigated through Taylor series approximation, here, we propose that the average irradiance of the HG beam under an arbitrary magnitude pointing error can be deduced through a statistical averaging method. We ﬁrstly found that the average irradiance proﬁle of an HG beam ﬁnally changes into an approximately Gaussian shape with the increase in pointing error and propagation distance and a larger beam waist at the transmitter could mitigate the proﬁle change. The correlation coefﬁcient between deduced theoretical expression and Monte Carlo simulation reaches 0.999. Additionally, the effective spot size, location of the local extreme value, average received power and signal-to-noise ratio (SNR) loss for an HG beam under pointing error were theoretically deduced and analyzed for the ﬁrst time. We found that the effective spot size of the higher-order HG beam experiences less broadening under the pointing error than that of the lower-order HG beam. The fundamental theoretical expressions of average irradiance for an HG beam under pointing error have provided effective guidance for analyzing the propagation characteristics and link performance.


Introduction
Free-space optical communication (FSOC) has the advantages of high bandwidth, high data rate, low energy consumption and high confidentiality and is considered a promising technology for next-generation communication [1][2][3][4].
A Hermite-Gaussian (HG) beam, as a type of higher-order TEM mn mode Gaussian beam, can be generated using higher-order solutions of the paraxial equation with Hermite polynomials (CO 2 laser) in rectangular coordinates [5]. A higher-order HG beam can form a multiple-spot pattern for the irradiance rather than a single spot as generated by lowest-order Gaussian beam, and furthermore, the HG beam has unique orthogonal spatial modes and well-preserved irradiance distribution in free-space propagation [5][6][7][8][9][10]. The definition of the HG beam's spot size is different from that of the Gaussian beam which was proposed by Carter [11]. For FSOC system application, the HG beam has garnered attention due to its potential to improve the capacity using mode-division multiplexing (MDM) [12][13][14][15]. Meanwhile, the HG beam can be also used in FSOC with a single-input multiple-output (SIMO) system owing to its propagation characteristics. Therefore, research on the propagation characteristics and expressions of HG beams, which include average irradiance, spot size and location of the local extreme value in the average irradiance for the HG beam, is critical for HG beam applications in FSOC.
Notably, the mobile FSOC link must be maintained by the beam-tracking mechanism, and the pointing error comes from its photoelectric detector noise and platform vibration residual noise which cause the propagation characteristics and irradiance distribution to be different from previous research [16][17][18]. Hence, the pointing error of the beam-tracking mechanism inevitably influences the irradiance distribution, which is one of the most critical characteristics of an HG beam. Kiasaleh first studied the influence of small pointing error (approximately 2.5 µrad with 3000 m) on the average irradiance distribution in free space through a Taylor series approximation [8]. However, the expression of the average irradiance under a small pointing error might become inaccurate as the pointing error increases. Due to the different application requirements of FSOC link systems, different pieces of equipment are selected as the coarse pointing and fine pointing systems. The magnitudes of pointing error are equal to several microradians (µrad) to several tens of microradians (µrad) [4], and a general expression for arbitrary magnitude pointing error is needed.
In this paper, we adopted the statistical averaging method instead of Taylor series approximation to deduce the expression of average irradiance for an HG beam under arbitrary pointing error. Based on average irradiance, the propagation characteristics, average received power and signal-to-noise ratio (SNR) loss of an HG beam have been investigated.
The paper is organized as follows: In Section 2, the statistical averaging method is adopted to establish the modeling of the average irradiance for the HG beam under pointing error. In Section 3, Monte Carlo simulation is used to verify the results of the deduced theoretical expression and prove its accuracy. In Section 4, based on average irradiance, the effective spot size, location of local extreme values in the average irradiance, average received power and signal-to-noise ratio (SNR) loss are discussed. Finally, the conclusions are presented in Section 5.

Modeling of the Average Irradiance under Pointing Error
The irradiance of a higher-order TEM mn HG beam can be treated as a generalized eigenfunction of the optical field equation in free space, which is obtained as follows [5,8]: where (x, y) denotes the Cartesian coordinates at distance L; x and y represent the horizontal and vertical directions, respectively; and H m (x) is a Hermite polynomial of the order m. In particular, for m = n = 0, the HG beam is reduced to the lowest-order TEM 00 Gaussian beam. Additionally, W 0 denotes the TEM 00 beam waist at the transmitter, and W denotes the TEM 00 spot size at the receiver, which is W = W 0 1 + Λ 2 0 , where Λ 0 = 2L/(kW 2 0 ), k = 2π/λ denotes the optical wave number and λ is the wavelength. Owing to the orthogonal spatial modes of the HG beam, its irradiance in Equation (1) allows for splitting into the horizontal (x-coordinate) and vertical (y-coordinate) direction components [5][6][7][8][9][10], and we observe the following: Figure 1, the pointing error causes the center of the HG beam to move randomly at the receiving plane. For convenience, a TEM 11 HG beam was selected to illustrate this process. The term S denotes the center of the HG beam at the transmitter, and O denotes the center of the receiver plane. The term SO denotes the propagation axis without pointing error (Figure 1a). The center of the HG beam moves to position O j under the pointing error angle θ j at time ∆t j , and the cases of j = 1, 2 and 3 are shown in Figure 1b. The term SO j (j = 1, 2, 3, . . . ) denotes the propagation direction under the pointing error, and θ j denotes the pointing error angle between SO j and SO. The average irradiance under the pointing error can be obtained by statistically averaging multiple HG beams at different positions ( Figure 1c). Because the probability of appearance at a specific position is determined by the probability density function (PDF) of the pointing error angle, the statistical averaging results can be formulated through position probability weighting for the HG beam irradiance at different positions (see Figure 1c).
illustrate this process. The term S denotes the center of the HG beam at the transmitter, and O denotes the center of the receiver plane. The term SO denotes the propagation axis without pointing error (Figure 1a). The center of the HG beam moves to position Oj under the pointing error angle θj at time Δtj, and the cases of j = 1, 2 and 3 are shown in Figure  1b. The term SOj (j = 1, 2, 3, …) denotes the propagation direction under the pointing error, and θj denotes the pointing error angle between SOj and SO. The average irradiance under the pointing error can be obtained by statistically averaging multiple HG beams at different positions (Figure 1c). Because the probability of appearance at a specific position is determined by the probability density function (PDF) of the pointing error angle, the statistical averaging results can be formulated through position probability weighting for the HG beam irradiance at different positions (see Figure 1c). The pointing error angle can be generally composed of two independent random variables along the horizontal and vertical axes, and the PDF of the pointing error angle for both the horizontal and vertical axes follows a nonzero-mean Gaussian distribution [5,15]. The pointing error-induced displacement is the product of the pointing error angle and propagation distance, which also follows a zero-mean Gaussian distribution for both the horizontal and vertical axes as follows [5,16,17]: where Then, the statistical averaging results for multiple HG beams at different positions can be expressed as follows: By substituting Equation (1) into Equation (4), the average irradiance of the HG beam under the pointing error can be rewritten as follows: The pointing error angle can be generally composed of two independent random variables along the horizontal and vertical axes, and the PDF of the pointing error angle for both the horizontal and vertical axes follows a nonzero-mean Gaussian distribution [5,15]. The pointing error-induced displacement is the product of the pointing error angle and propagation distance, which also follows a zero-mean Gaussian distribution for both the horizontal and vertical axes as follows [5,16,17]: where r x = θ x L and r y = θ y L denote displacements along the horizontal and vertical axes, respectively, and θ x and θ y represent the pointing error angles along the horizontal and vertical directions, respectively. Correspondingly, σ r = σ θ L denotes the standard variance of r x and r y , and σ θ represents the standard variance of θ x and θ y . Then, the statistical averaging results for multiple HG beams at different positions can be expressed as follows: By substituting Equation (1) into Equation (4), the average irradiance of the HG beam under the pointing error can be rewritten as follows: where I m (x, L) PE = ∞ −∞ I m (x − r x , L) × f (r x )dr x and I n (y, L) PE = ∞ −∞ I n (y − r y , L) × f (r y )dr y . Equation (5) is composed of two separate components determined by m and n, respectively. However, Equation (5) cannot be integrated when m and n are uncertain. For an FSOC system, the values of m and n are commonly equal to 1, 2 and 3 [8][9][10]. When the values of m and n are certain, the Hermite polynomials have a certain expression, and Equation (5) can be deduced as a closed-form expression. Table 1 summarizes the integration results of I p (s, L) PE for p = 1, 2 and 3, where p represents either m or n, and s represents either x or y. This article only provides the expression for the orders equal to 1, 2 and 3, but the expression for higher orders can also be derived using the statistical averaging method of this paper. Table 1. Expression of I p (s, L) PE for p = 1, 2 and 3.
Thus, the closed-form expression for the average irradiance of the HG beam under a pointing error can be formulated by multiplying I m (x, L) PE and I n (y, L) PE . For instance, the average irradiance of the TEM 21 HG beam under the pointing error is equal to the product of I 2 (x, L) PE and I 1 (y, L) PE as follows:

Numerical Simulation Results and Verification
The Monte Carlo simulation is widely used in FSOC system link simulation [19,20]. In this section, we adopted Monte Carlo simulation to design the propagation process of an HG beam under pointing error.
The random variation at the receiver caused by the pointing error can be generated using the Monte Carlo method. The detailed simulation process is shown as follows: Combined with Equations (1)-(3), the instantaneous HG beam irradiance deviating under the influence of pointing error can be calculated as follows: where r x and r y are generated by Equations (2) and (3), respectively. Then, the simulated average irradiance under pointing error can be obtained statistically, which is shown as follows: where j denotes a serial number from 1 to M and I p mn (x, y, L) denotes the j-th simulated speckle irradiance moving under the j-th simulated pointing error.
The numerical simulation parameters were set as follows: beam waist at transmitter W 0 = 0.05 m, wavelength λ = 850 nm and propagation distance L changes from 1000 to Photonics 2022, 9, 478 5 of 14 3000 m. According to Equations (2) and (3), 10,000 groups of random displacements along the horizontal and vertical axes were generated, and the standard variance of the pointing error angle was set as 0, 5, 10 and 20 µrad. The simulation parameters are set in Table 2; they are often used in practical FSOC systems. For a satellite-satellite FSOC link, owing to the difference in the selection of tracking mechanism, the magnitudes of pointing error are equal to several microradians (µrad) to several tens of microradians (µrad) [4]. Especially for small satellite FSOC links, the fine pointing of small satellite FSOC links is primarily dependent on the fast-steering mirror (FSM) of microelectromechanical systems (MEMSs) with lower precision and smaller size, and the pointing error commonly reaches several tens of microradians (µrad) [21]. In this condition, the deduced average irradiance for the HG beam under arbitrary pointing error can provide guidance for theoretical analysis.
In Figures 2-4, the lines represent the theoretical equations in Table 1, and the dots represent the results of the Monte Carlo simulation. Red, blue, green and yellow represent the results of the standard variance of the pointing error angle equal to 0, 5, 10, and 20 µrad, respectively. The correlation coefficients between the theoretical and simulation results were all greater than 0.999, which proves that the simulation results coincide with the deduced closed-form expressions. Notably, the numerical results of the proposed model are the same as Kiasaleh's results under a small pointing error, and we provide the results of a large pointing error. , x FOR PEER REVIEW 6 of 16 increase in pointing error and propagation distance was first discovered. Additionally, the profile variation of a higher-order HG beam is slower than that of a lower-order HG beam under the same propagation conditions.  In the absence of the pointing error, the HG beam irradiance distribution is well preserved with an increase in the propagation distance when p = 1, 2,and 3. With the increase in pointing error and propagation distance, the valleys of the HG beam profile rise and converge toward the central origin, the profile changes into concave shape and then from concave shape to an approximately flat-top shape and, finally, the profile changes into an approximately Gaussian shape. This variation of the average irradiance profile with the increase in pointing error and propagation distance was first discovered. Additionally, the profile variation of a higher-order HG beam is slower than that of a lower-order HG beam under the same propagation conditions.    The average irradiance of the HG beam under the pointing error with different beam waists at the transmitter is illustrated in Figure 5. The beam waist at the transmitter was set as 0.05, 0.1 and 0.2 m. The standard variance of the pointing error was 20 μrad, and the propagation distance L = 2000 m. The other parameters were the same as those described initially. The lines represent the theoretical equations listed in Table 1, and the dots represent the results of the Monte Carlo simulation. Red, blue and green represent the beam waists at transmitter W0 = 0.05, 0.1, and 0.2 m, respectively. Figure 5a-c shows the cases of p = 1, 2 and 3, respectively. We found that a larger beam waist at the transmitter could mitigate the profile change of the average irradiance of the HG beam under an increasing pointing error. Therefore, the original profile of the HG beam can be maintained by  Figure 5a-c shows the cases of p = 1, 2 and 3, respectively. We found that a larger beam waist at the transmitter could mitigate the profile change of the average irradiance of the HG beam under an increasing pointing error. Therefore, the original profile of the HG beam can be maintained by properly increasing the beam waist and controlling the standard variance of pointing error.
The average irradiance of the HG beam under the pointing error with different beam waists at the transmitter is illustrated in Figure 5. The beam waist at the transmitter was set as 0.05, 0.1 and 0.2 m. The standard variance of the pointing error was 20 μrad, and the propagation distance L = 2000 m. The other parameters were the same as those described initially. The lines represent the theoretical equations listed in Table 1, and the dots represent the results of the Monte Carlo simulation. Red, blue and green represent the beam waists at transmitter W0 = 0.05, 0.1, and 0.2 m, respectively. Figure 5a-c shows the cases of p = 1, 2 and 3, respectively. We found that a larger beam waist at the transmitter could mitigate the profile change of the average irradiance of the HG beam under an increasing pointing error. Therefore, the original profile of the HG beam can be maintained by properly increasing the beam waist and controlling the standard variance of pointing error.

Effective Spot Size
Different from the lowest-order TEM00 Gaussian beam, the HG beam forms a pattern of spots rather than a single spot of light; the conventional spot size W of a Gaussian beam

Effective Spot Size
Different from the lowest-order TEM 00 Gaussian beam, the HG beam forms a pattern of spots rather than a single spot of light; the conventional spot size W of a Gaussian beam is not suitable for an HG beam, and a new definition of the spot size of an HG beam is needed.
Carter et al. defined the "effective spot size" of an HG beam without pointing error [5,7,11]: The terms ρ x,m (L) and ρ y,n (L) represent the spot size of the p-th mode [5,7,11]: where p represents either m or n, and s represents either x or y. Based on Equations (1) and (10), ρ x,m (L) and ρ y,n (L) can be rewritten as follows [5,7]: In the presence of a pointing error, substituting Equation (5) into Equation (10), the spot size under the pointing error can be deduced as follows: Equations (13) and (14) show that the pointing error increases the spot size of the HG beam. The spot size ratio is defined as ρ 2 s,p (L)/ρ 2 s,p (L) PE , where p represents either m or n, s represents either x or y, and m can be replaced by n with the same result. In Figure 6, red, blue and green represent the beam waist at transmitter W 0 equal to 0.05, 0.1 and 0.2 m, respectively. We found that the effective spot size of the higher-order HG beam experiences less broadening under the pointing error than that of the lower-order HG beam. Meanwhile, the pointing error-induced HG beam spot size broadening decreased with an increase in the beam waist at the transmitter.
Equations (13) and (14) show that the pointing error increases the spot size of the HG beam. The spot size ratio is defined as , where p represents either m or n, s represents either x or y, and m can be replaced by n with the same result. In Figure 6, red, blue and green represent the beam waist at transmitter W0 equal to 0.05, 0.1 and 0.2 m, respectively. We found that the effective spot size of the higher-order HG beam experiences less broadening under the pointing error than that of the lower-order HG beam. Meanwhile, the pointing error-induced HG beam spot size broadening decreased with an increase in the beam waist at the transmitter.

Location of the Local Extreme Value in the Average Irradiance under Pointing Error
The higher-order HG beam forms a pattern of spots rather than a single spot of an optical beam. Therefore, the location of each spot is necessary for the formulation of the receiving scheme.
According to the expressions in Table 1, the location of the local extreme value in the average irradiance under the pointing error can be obtained by solving the equation d I p (s, L) PE /ds = 0 and determining by the positive and negative of equation d 2 I p (s, L) PE /ds 2 , where p represents either m or n, and s represents either x or y.
The locations of the local maxima and minima in the average irradiance under the pointing error are listed in Table 3. Table 3. Locations of local maxima and minima value in the average irradiance under pointing error. The term q in Table 3 represents an expression associated with W and σ r and is obtained as follows:

Locations of the local maxima
Because m and n represent the order in the x and y axes, the coordinates of local extreme value locations can be determined by the location along the x and y axes. Taking I 21 (x, y, L) PE for instance, the locations of local maxima and local minima of irradiance are shown in Table 4. Table 4. Locations of local maxima and minima for TEM 21 Hermite-Gaussian beam.
In the presence of pointing error Coordinates of the local maximum irradiance , y .

Coordinates of the local minimum irradiance
In the absence of pointing error Coordinates of the local maximum irradiance Coordinates of the local minimum irradiance

Average Received Power and SNR Loss
Based on the aforementioned research on the average irradiance of the HG beam under the pointing error and location of the local maxima, the average received power at the location of the local maxima can be further deduced for the HG-based FSOC with (SIMO).
As shown in Figure 7, the three common modes of the HG beam, TEM 11 , TEM 22 and TEM 33 , are selected for the average received power investigation. Each array receiving the aperture center is located at the location of the local maxima. The TEM 11 HG beam irradiance without pointing error had four local maximum points, and their peak irradiances were equal. Therefore, the receiving aperture D 0 is arranged at four locations of the local maxima. The TEM 22 HG beam without pointing error has nine local maximum points, which can be divided into three groups according to their peak irradiance. Therefore, the receiving apertures D 1 , D 2 and D 3 were arranged at nine locations of the local maxima (see Figure 7b). Similarly, the TEM 33 HG beam without pointing error has 16 local maxima points, and it can also be divided into three groups according to their peak irradiance: the receiving apertures D 1 ', D 2 ' and D 3 ' are arranged at 16 locations of local maxima (see Figure 7c). Taking the TEM 33 HG beam as an example, the receiving apertures are the same but are divided into three groups according to the peak irradiance (red, yellow and orange), and the average received power values of similar color receiving apertures are equal. Figure 7c). Taking the TEM33 HG beam as an example, the receiving apertures are the same but are divided into three groups according to the peak irradiance (red, yellow and orange), and the average received power values of similar color receiving apertures are equal.
where D denotes the diameter of the circular receiving aperture. However, Equation (16) cannot be integrated over a circular area when the receiving aperture is not located at the center of the main propagation axis. Inspired by [21,22], we use a square within the same area to substitute the circle as the integration area. The square length l is equal to 0.5 D π . Assuming that the coordinates of the center of the receiving aperture are (xf, yf), Equation (16) can be approximated as follows: where D denotes the diameter of the circular receiving aperture. However, Equation (16) cannot be integrated over a circular area when the receiving aperture is not located at the center of the main propagation axis. Inspired by [21,22], we use a square within the same area to substitute the circle as the integration area. The square length l is equal to 0.5 √ πD. Assuming that the coordinates of the center of the receiving aperture are (x f , y f ), Equation (16) can be approximated as follows: When taking the values of m and n, the expressions for the average received power corresponding to TEM 11 , TEM 22 and TEM 33 are as follows: where and er f () is an error function.
As shown in Figure 8a, for the TEM 11 HG beam, the average received power on the receiving aperture D 0 decreases with an increase in the pointing error. As shown in Figure 8b, for the TEM 22 HG beam, the average received power on D 1 and D 2 decreased with an increase in the pointing error, and the average received power on D 1 decreased faster than the average received power on D 2 . The average received power on D 3 fluctuated with an increase in the pointing error. Notably, the average received powers on D 1 , D 2 and D 3 are equal when σ r = 12.67 µrad; we define it as the "equal power point". As shown in Figure 8c, for the TEM 33 HG beam, the variation trend of the average received power on D 1 ', D 2 ' and D 3 ' with an increase in the pointing error was the same as that of the TEM 22 HG beam. The "equal power point" appears at σ r = 14.05 µrad. As shown in Figure 9, the average irradiance of the TEM22 HG beam at the "equal power point" r σ = 12.67 μrad exhibits a flat-topped profile. As shown in Figure 10, the average irradiance of the TEM33 HG beam at the "equal power point" r σ = 14.05 μrad exhibits a slightly hollow shape. The average irradiance of the HG beam appears to have a nearly flat-topped shape at the "equal power point". As shown in Figure 9, the average irradiance of the TEM 22 HG beam at the "equal power point" σ r = 12.67 µrad exhibits a flat-topped profile. As shown in Figure 10, the average irradiance of the TEM 33 HG beam at the "equal power point" σ r = 14.05 µrad exhibits a slightly hollow shape. The average irradiance of the HG beam appears to have a nearly flat-topped shape at the "equal power point". As shown in Figure 9, the average irradiance of the TEM22 HG beam at the "equal power point" r σ = 12.67 μrad exhibits a flat-topped profile. As shown in Figure 10, the average irradiance of the TEM33 HG beam at the "equal power point" r σ = 14.05 μrad exhibits a slightly hollow shape. The average irradiance of the HG beam appears to have a nearly flat-topped shape at the "equal power point".  The TEM33 HG beam is taken as an example to analyze the received power loss ratio (also SNR loss ratio) under the influence of pointing error in the corresponding receiving aperture. The numerical simulation parameters are the same as above, and the transmission distance is 3000 m.
As shown in Figure 11, with the increase in pointing error, the loss of average SNR is 3.00 dB when r σ = 7.89 μrad at receiving aperture D1', the loss of average SNR is 2.91 dB when r σ = 20 μrad at receiving aperture D2' and the average SNR decaying speed of receiving aperture D1' is faster than that of receiving aperture D2'. For receiving aperture D3', the average SNR decreases and then increases with the increase in pointing error, and the increase in average SNR is 0.87 dB when r σ = 16.97 μrad. This shows that the transmitted power is homogenized at the receiver under the influence of pointing error. In order to achieve better reception efficiency for the HG beam, we should control the transmission parameters. It can be seen that when the propagation parameter is set as shown in Table 2, the standard variance of pointing error should be less than 7.5 μrad for the TEM33 HG beam so that the SNR loss at each receiving aperture is less than 3.00 dB. The TEM 33 HG beam is taken as an example to analyze the received power loss ratio (also SNR loss ratio) under the influence of pointing error in the corresponding receiving aperture. The numerical simulation parameters are the same as above, and the transmission distance is 3000 m.
As shown in Figure 11, with the increase in pointing error, the loss of average SNR is 3.00 dB when σ r = 7.89 µrad at receiving aperture D1', the loss of average SNR is 2.91 dB when σ r = 20 µrad at receiving aperture D 2 ' and the average SNR decaying speed of receiving aperture D 1 ' is faster than that of receiving aperture D 2 '. For receiving aperture D 3 ', the average SNR decreases and then increases with the increase in pointing error, and the increase in average SNR is 0.87 dB when σ r = 16.97 µrad. This shows that the transmitted power is homogenized at the receiver under the influence of pointing error. In order to achieve better reception efficiency for the HG beam, we should control the transmission parameters. It can be seen that when the propagation parameter is set as shown in Table 2, the standard variance of pointing error should be less than 7.5 µrad for the TEM 33 HG beam so that the SNR loss at each receiving aperture is less than 3.00 dB. r mitted power is homogenized at the receiver under the influence of pointing error. In order to achieve better reception efficiency for the HG beam, we should control the transmission parameters. It can be seen that when the propagation parameter is set as shown in Table 2, the standard variance of pointing error should be less than 7.5 μrad for the TEM33 HG beam so that the SNR loss at each receiving aperture is less than 3.00 dB.

Conclusions
The closed-form expression of the average irradiance of the HG beam under the pointing error was deduced through a statistical averaging method and verified by Monte Carlo simulation. The correlation coefficient between deduced theoretical expression and Monte Carlo simulation reaches 0.999. We firstly found that the valleys of the HG beam profile rise and converge toward the central origin, the profile changes into a concave shape and then from a concave shape to an approximately flat-top shape and, finally, the profile changes into an approximately Gaussian shape, and the larger beam waist at the transmitter could mitigate the profile change of the average irradiance of the HG beam under increasing pointing error. The effective spot size, location of the local extreme value and average received power for the average irradiance under the pointing error were theoretically deduced for the first time. The spot size will be broadened under the influence of the pointing error, the location of the local extreme value will change with the variation of the pointing error and the average received power at the local extreme value fluctuates under different pointing errors. Moreover, the pointing error causes SNR loss in the receiving aperture that has been analyzed. For the TEM33 HG beam, the standard variance of pointing error should be within 7.5 μrad so that the SNR loss at each receiving aperture is less than 3.00 dB with the transmission distance of 3000 m. The fundamental theoretical expressions of average irradiance for an HG beam under pointing error have provided effective guidance for analyzing the propagation characteristics and link performance.
Further research on the average irradiance for an HG beam under pointing error in atmospheric turbulence should be conducted. Under the influence of atmospheric turbulence, the optical field is disturbed and more random processes are introduced. Further research may prove theoretical guidance for HG beams used in atmospheric turbulence.

Conclusions
The closed-form expression of the average irradiance of the HG beam under the pointing error was deduced through a statistical averaging method and verified by Monte Carlo simulation. The correlation coefficient between deduced theoretical expression and Monte Carlo simulation reaches 0.999. We firstly found that the valleys of the HG beam profile rise and converge toward the central origin, the profile changes into a concave shape and then from a concave shape to an approximately flat-top shape and, finally, the profile changes into an approximately Gaussian shape, and the larger beam waist at the transmitter could mitigate the profile change of the average irradiance of the HG beam under increasing pointing error. The effective spot size, location of the local extreme value and average received power for the average irradiance under the pointing error were theoretically deduced for the first time. The spot size will be broadened under the influence of the pointing error, the location of the local extreme value will change with the variation of the pointing error and the average received power at the local extreme value fluctuates under different pointing errors. Moreover, the pointing error causes SNR loss in the receiving aperture that has been analyzed. For the TEM 33 HG beam, the standard variance of pointing error should be within 7.5 µrad so that the SNR loss at each receiving aperture is less than 3.00 dB with the transmission distance of 3000 m. The fundamental theoretical expressions of average irradiance for an HG beam under pointing error have provided effective guidance for analyzing the propagation characteristics and link performance.
Further research on the average irradiance for an HG beam under pointing error in atmospheric turbulence should be conducted. Under the influence of atmospheric turbulence, the optical field is disturbed and more random processes are introduced. Further research may prove theoretical guidance for HG beams used in atmospheric turbulence. The average bit error rate and fade probability for an HG beam-based FSOC link should also be studied in following work.