Effect of Surface Error Distribution and Aberration on Electromagnetic Performance of a Reflector Antenna

In this paper, a method based on the Zernike distribution and the optical aberration is proposed to investigate the effects of the distribution characteristics of surface distortions of a reflector antenna on its electromagnetic performance (EMP). For large-scale errors, an analytical model is introduced to describe the arbitrary distortions, based on the orthogonal Zernike polynomials. The effects of error distribution described by the Zernike series on typical EMP are analyzed. The numerical results indicate that the distortions in the distribution feature of defocus or spherical aberration have a greater impact on gain, and the distribution feature of tilt or coma mainly influences boresight offset; the distribution feature of defocus, astigmatism, and spherical aberration has a greater impact on sidelobe levels. The results indicate that the beam contour patterns are related to the distribution forms of distortions and are similar with the same aberration feature. On the basis of the Seidel aberration, the relationships between typical EMP and the aberration coefficient are presented. Based on these, the error profile of the primary influence and the approximate feature of EMP can be determined, and the antenna performance can be predicted in a simple manner.


Introduction
With the development of radio astronomy and deep space exploration, the requirements for the resolution and sensitivity of the reflector antenna are raised.This requires antennas to have a larger aperture size and better electromagnetic performance (EMP) to facilitate obtaining high-quality scientific data [1,2].It is known that the surface quality is one of the key parameters used to evaluate an antenna's mechanical performance [3,4].Reflector surface errors are inevitable, and they may be introduced by manufacturing tolerances or structural deformations caused by environmental loads (e.g., gravity, temperature, wind, and accretions of snow and ice) [5,6].Reflector surface errors or aperture phase errors directly influence the antenna gain, boresight pointing, sidelobe level, and beam shape [7].Many EMP degradations depend on the magnitude, shape, distribution of errors, aperture illumination taper, etc.Hence, understanding the effect of errors on an antenna's EMP is very important.
In general, the mathematical description of the surface error model and the radiation integral can provide a direct connection between structural design methods such as finite element analysis and the resulting EMP of the antenna.However, for convenience, the root mean square (RMS) of surface errors is still currently considered as a significant metric to evaluate the mechanical behavior of antennas [4].Many studies focus on the determination of the loss in antenna gain from the RMS value.Ruze's formula [8] can be used effectively as a simple estimate and is very critical.However, the gain loss can only be predicted from the knowledge of the RMS of surface errors; other EMP features, such as the pointing offset and sidelobe level, cannot be obtained from the RMS of errors.At present, antenna design demands that the estimates of boresight pointing offset, sidelobe level, and beam pattern deformation must be obtained rapidly and accurately.These estimates are very important to evaluate an antenna's overall sensitivity.
Besides Ruze, many other researchers have carried out studies on relevant topics.Rahmat-Samii [9] modified Ruze's formula considering the focal ratio (F/D) and the edge taper: ΔG = 10 log e − 4πk g ε/λ 2 , 1 where ε is the RMS of half of the optical pathlength difference (OPD) and k g is a correction factor that depends on the F/D and the edge taper.In that paper, the effects of random errors on the average power pattern of reflectors were analyzed by using geometrical optics.Corkish [10] presented a survey of the literature relating to the effects of reflector surface distortions on sidelobe levels.The errors are classified on the basis of whether they are random or systematic.Systematic errors are further considered as periodic or nonperiodic errors.Bahadori and Rahmat-Samii [11] studied the influence of periodic and aperiodic surface distortions of cylindrical reflector antennas on the performance of membrane reflector antennas described by a sinusoidal function, such as Δz = Δ⋅ sin 2πNx/D .Bahadori proposed a closed-form expression to predict some of the unique features of these patterns including the location and the level of the observed grating lobes.Duan and Wang [12] studied the influence of surface random errors and systematic errors on the EMP and established an optimization model for the integrated mechanical-electromagnetic performance.Lian et al. [13] developed a multiobjective optimization model to analyze the effects of surface errors' nonuniformity along the radius on the reflector's radiation characteristic and discussed that for the same RMS of errors for the whole reflector, different distributions along the radius result in different gain losses and sidelobe levels.However, these works consider only one type of error distribution, such as the variations in radial or azimuthal directions, and hence, they cannot represent all situations that influence of the arbitrary error distributions.For large ground-based reflector antennas, surface distortions are very complicated.With the exception of manufacturing tolerances and misalignment of panels, the backup structure and panels will be deformed due to gravity, temperature, wind, and accretions of snow and ice.Their distribution forms and effects on the EMP are different.The purpose of this paper is to study the effect of different error distribution characteristics from the point of view of the scope of the influence and the aberration.The errors are divided into three types based on different scales of the correlation interval.Zernike polynomials are used to describe wavefront error distributions.These polynomials are completely orthogonal, and their combination can describe arbitrary distortions.The effects of error distributions with the Zernike series on gain loss, boresight pointing offset, and sidelobe level increase are analyzed by geometrical optics.These effects of the Zernike distribution errors are classified and further discussed based on the Seidel aberration.Based on many numerical simulations, the relationships between a typical EMP and the aberration coefficients are presented.We think that the aberration coefficients can be considered as a metric to evaluate the wavefront error behavior.It means that in the presence of wavefront distortions, one can determine the error profile and the approximate feature of the EMP based on the decomposition of the Zernike series and analysis of the aberration.We think that this method is useful for structural design and performance evaluation of a reflector antenna.

Classification in the Correlation Scale for Typical Errors
Reflector surface distortions can be often defined on the basis of the characteristic correlation interval, C, for the phase errors or OPDs (also called wavefront errors) in the aperture plane.It is assumed that the OPD values are completely correlated within a radius of C and completely uncorrelated otherwise [8,14].It is an existent and realistic assumption for segmented surfaces of large ground-based reflector antennas and is a practical way of defining surface irregularity.
From the point of view of the correlation interval scale, the distortions can be categorized into three main subdivisions: small-scale errors, intermediate-scale errors, and large-scale errors, shown in Figure 1.For the segmented surface of reflector antennas, the distinction of different scales is related to the aperture diameter, D, and length of the reflector panels, L. It is assumed as follows: If C ≤ L, the surface irregularity can be distributed in the high spatial frequency form in the aperture surface, as shown in Figure 1(a), and the errors are correlated over areas of the order of the individual panel or less.They may be introduced by manufacturing tolerances, fabrication errors in the reflector panels, and local deformation of the individual panels due to gravity or temperature as well as the errors in the alignment of the panels to the prescribed reflector contour.The small-scale errors are generally considered randomly distributed as panel fabrication and setting errors.The magnitude of errors is usually much less than the wavelength.
If L < C ≤ D/2, the surface errors can be distributed in the medium spatial frequency form in an aperture surface, as shown in Figure 1(b), and the errors are correlated over areas of the order of the several panels to half of the circle aperture.They may be caused by local deformation of the backup structure under the influence of local nonuniform temperature force distribution.The time of occurrence and the influence region of the temperature force distribution are random and nonrepeatable as the change of the antenna working case (azimuth and elevation) and environmental conditions.Consequently, the medium-scale errors are generally considered randomly distributed.
If D/2 ≤ C, the surface errors can be distributed in the low spatial frequency form in the aperture surface, as shown in Figure 1(c), and the errors are correlated over areas of the 2 International Journal of Antennas and Propagation order of almost the entire circle aperture.They may be caused by deformation of the backup structure due to gravity or wind.The appearance of wind is usually random, and its effect may be rapid variations.As gravitational loading is constant and its effects are repeatable, these errors are systematic.In this study, large-scale errors only refer to the gravitational deformation of the backup structure and will affect almost the entire reflector surface.In this paper, the errors will be considered on the basis of whether they are random or systematic, since that is what determines the method of treatment.Unless the random errors of a particular reflector have been characterized by measurement, they must be treated in a probabilistic way to find the most likely effects.Systematic errors can be treated in a deterministic way.

Radiation Formulation
According to the Fourier transform relationship between the aperture field and the far field, the radiation integral of a standard reflector can be expressed as where F r, θ is the aperture field distribution function on the aperture surface A; Δφ is the phase error distribution function in the aperture plane caused by surface errors; ψ, γ are the coordinates of the observation direction in far field; r, θ are the coordinates of the point in the aperture plane; k is the free space wave constant, k = 2π/λ; and λ is the working wavelength.The schematic diagram of the reflector antenna is shown in Figure 2.
The phase error is proportional to the OPD and is given by where δ is the OPD caused by surface deformations.The OPD at a point on a reflector surface is greatly exaggerated in Figure 3.It can be seen in Figure 3 that a point p 0 on the ideal surface is changed to p 0 ′ , and the ray pathlength is increased by the sum of the distance of L1 (along the incident ray) and L2 (along the reflected ray) [15]; moreover, the OPD where d 0 is the surface deformation vector, d 0 = Δdx, Δdy, Δdz T ; N 0 is a unit normal vector; and Δdn is the normal component of the deformation vector.
Equation 4 can be used for a primary reflector to compute the OPD for any deformed surface in the normal vector and the axial or normal component of the deformation vector.In this paper, we mainly focus on the descriptions of different scale surface errors and the effects on EM preformation, and the surface error will be expressed in the OPD in general.
In this study, the patterns will be obtained by equation ( 2) with numerical calculation.In order to verify the validity of EMP results, we analyze the far-field patterns with the numerical calculation and commercial electromagnetic software named GRASP, respectively.The comparison of farfield patterns with radiation integral and GRASP software is shown in Figure 4. We can see that the results from numerical calculation and GRASP software are almost consistent; there are only little differences in null depth and side lobes in the far region.So it indicates that the numerical results are effective and accurate.

Modeling and Analysis of Large-Scale Errors
The large-scale systematic errors can be described by combination of trigonometric and radial polynomials.It can describe the errors distributed in the radial and azimuthal directions, respectively.However, these polynomial combinations may not be completely orthogonal and cannot accurately describe a complicated distortion, which may be obtained by finite element analysis or measured data.
In general, a set of orthogonal polynomials is needed to model the complicated distortion, and the polynomials have attractive mathematical properties in representing complex surfaces in the form of linear expansion of their basis functions.Recently, almost all researchers in the field of optic telescope have analyzed an arbitrary surface deformation or wavefront/OPD described by a linear combination of a set of orthogonal polynomial series named Zernike polynomials.

Zernike Polynomial.
The Zernike polynomial is a set of complete orthogonal polynomials defined on a unit circle, and the low-level terms correspond with the optical primary aberrations [16][17][18].The polynomials are usually expressed in polar coordinates (r, θ) as products of angular functions and radial polynomials, where r is the normalization radial coordinate ranging from 0 to 1 and θ is the azimuthal 4 International Journal of Antennas and Propagation component ranging from 0 to 2π.In general, the Zernike polynomials are defined as where the radial function R m ′ n r is described by  1, and the corresponding shapes of the first 21 polynomials (excluding the first term) are shown in Figure 5.A significant property of the Zernike polynomials is the orthogonality [17], which is expressed by where if j = j ′ , δ jj′ = 1, and if j ≠ j ′ , δ jj′ = 0.The normalization property of the Zernike polynomials is given by From Zernike circle polynomials, the associated RMS can be calculated, and the Zernike RMS polynomials where A is the area of a unit circle.The Zernike RMS k j is listed in Table 1, and 1/k j is usually called the normalization factor.The reflector surface distortion can be expressed in terms of a sequence of Zernike polynomials in polar coordinates: where a i is the coefficient of Zernike polynomials.It is useful to decompose any distortion function (wavefront errors) into components represented by the primary aberrations (tilt, defocus, astigmatism, coma, and spherical aberration), and each polynomial function of the Zernike expansion can provide a good description of each aberration component and extract a component of the wavefront errors.This provides an insight into the mechanical or structural causes of the error.
The RMS of distortion is given by where A is the area of the antenna's aperture.The RMS value can also be calculated from a vector of Zernike coefficients and the Zernike RMS values as follows: The Zernike polynomials can be normalized to ensure that the RMS value of the phase error caused by any polynomial is one.The normalized Zernike polynomials can be described as follows: The attractive features of the normalized Zernike polynomial representation are that for each polynomial or Zernike mode (apart from the piston term), the mean value of the aberration across the aperture is zero, and the value of the coefficient gives the RMS wavefront error for that particular Zernike mode.

Effects of Error Distributions
As mentioned above, the RMS of wavefront errors distributed in a normalized Zernike polynomial is unity; therefore, any wavefront errors distributed in the normalized Zernike polynomials would have the same RMS.Then, it is helpful to study how these errors influence the EMP.
The geometry of a parabolic reflector antenna with a large-scale distortion is depicted in Figure 5.It is assumed that the OPD is described by different normalized Zernike polynomials in the circle aperture.Furthermore, the reflector has a diameter D = 26 m and a focal length F = 7 8 m.The operation wavelength is λ = 0 06 m, and the aperture field is a Gaussian distribution with -12 dB edge taper.To obtain the obvious effect, the normalized Zernike polynomials are multiplied by a weight factor ke or in other words the actual distortion RMS.It is assumed ke = λ/20 = 0 3 cm, and the RMS of half of the distortion is λ/40.The estimated gain loss with RMS/λ = 1/40 can be calculated by using modified Ruze's formula (equation ( 1)), and the result is -0.32 dB.Then, the OPD is modeled by The far-field patterns with the OPD distributed in the first 21 normalized Zernike polynomials are computed by using the radiation formulation (equation 2).The results of the EMP-gain loss, boresight offset, and sidelobe level degradation-are listed in Table 2, and the far-field patterns in the 0 °and 90 °cut planes, shown in Figure 6, are normalized to the peak of the ideal reflector, without distortions.The beam pattern contours are shown in Figure 7.
The loss item in Table 2 shows that the gain losses are different.Note that the gain losses in z5 and z13 are more     The offset item in Table 2 shows that z2, z3, z8, and z9 have a greater impact on the boresight pointing offsets, and other distributions have almost no effect on the pointing offset.This means that the boresight offset is mainly influenced by tilt and coma.
As shown in Figure 6, almost all the sidelobe levels are increased to different extents.The results of sidelobe level degradation (ΔSLD) in Table 2 indicate that the sidelobe levels of z5 and z13 are increased to some extent, and these are symmetrical.The sidelobe levels of z6 and z14 are increased to some extent, and these are asymmetrical; it means that their sidelobe levels in 0 °and 90 °cut planes are different, but the sidelobe levels in the same cut plane are equal.The peak levels of sidelobe degradation of z8, z9, z18, and z19 are fairly high, and their levels of left and right sidelobes in the 0 °or 90 °cut plane are different; it means that one side is high and the other one is low.From the distribution characteristic, except the characteristics of defocus and spherical aberration (excluding piston), the other characteristics will lead to sidelobe level asymmetry.From the overall results of Table 2, we can see that the greatest influences on sidelobe levels are the coma (the peak of ΔSLD is about 10.7 dB) and 2nd-coma (the peak of ΔSLD is about 11.6 dB), followed by 2nd-astigmatism (the peak of ΔSLD is about 9.3 dB) and spherical aberration (the peak of ΔSLD is about 8 dB); next are defocus (the peak of ΔSLD is about 6.8 dB) and astigmatism (the peak of ΔSLD is about 6.1 dB), and the tilt has the least influence (the peak of ΔSLD is about 3.4 dB).
The beam contour patterns of the far-field radiation are shown in Figure 7. Patterns are plotted in ranges of magnitudes from "peak" to "peak-30 dB," and the magnitude difference of the adjacent contours is 3 dB.We can see in Figure 7 that the error distribution forms have a significant effect on pattern shapes and cause different distortions of pattern shapes.The shapes of patterns with the same aberration type are similar, such as the patterns caused by z2 and that caused by z3.
The results are shown in Figures 5 and 6, and we can know that the gain losses, sidelobe levels, and pattern shapes caused by wavefront errors with the same RMS value are different and are distributed with different distribution characteristics.It shows that the beam patterns are correlated with the distribution forms of distortions.Thus, the coefficient of Zernike polynomials can be applied in modeling the wavefront shapes to decompose the wavefront aberrations.

Effects of Aberrations
To further investigate the effect of the Zernike coefficient on EMP, the distributions of wavefront errors are analyzed from the optical aberration point of view.The Zernike decomposition provides a good connection between EMP and the wavefront error distribution, which represents the aberration component [19].Each polynomial function of the Zernike expansion extracts a component of the wavefront error; these components can be grouped into classical aberrations called Seidel aberration.
According to the Seidel aberration theory, the first-order wavefront properties and third-order wavefront aberration coefficients can be obtained from the Zernike polynomial coefficients [17].The nine Zernike polynomials are adopted to describe the Seidel aberration.The aberrations and properties corresponding to these Zernike terms are shown in   where a i represents the new Zernike polynomial coefficients.
Rewriting the new Zernike expansion of equation ( 15), the firstand third-order field-independent wavefront aberration terms can be obtained.This is achieved by grouping the similar terms and equating them with the wavefront aberration coefficients: W r, θ = a 0 − a 3 + a 8 piston + a 1 − 2a 6 r cos θ + a 2 − 2a 7 r sin θ tilt + 2a 3 − 6a 8 r 2 focus + a 4 r 2 cos 2θ + a 5 r 2 sin 2θ astigmatism + 3a 6 r 3 cos θ + 3a 7 r 3 sin θ coma + 6a 8 r 4 spherical 16 This equation can be rearranged corresponding to the field-independent wavefront aberration coefficients as follows: From equation (17), the magnitude, sign, and angle of these field-independent aberration terms are listed in Table 4, and the corresponding shapes of the Seidel wavefront aberration are shown in Figure 8.
As mentioned above, the effects of the normalized Seidel wavefront aberrations on EMP-which are multiplied by a weight factor (1/RMS)-are computed.In the analysis, the reference angles of tilt, astigmatism, and coma are set to zero.The results of EMP are listed in Table 5, and the far-field patterns in the 0 °and 90 °cut planes are shown in Figure 9 and are normalized to the peak of the ideal reflector.The beam contour patterns are shown in Figure 10.
The result of Table 5 indicates that the terms s2, s3, and s5 have a greater impact on gain loss; the terms s1 and s4 mainly have effect on the pointing offset and, moreover, cause the asymmetric sidelobes in the cut planes.
Based on many numerical simulations, the relationships between gain loss ΔG and the different aberration coefficients a have been determined and are shown in Figure 11.The relationship can be given by where h is the gain loss influence factor and a′ is the aberration coefficient normalized to the wavelength, a/λ.Through simulations, the gain loss influence factors for the focus, astigmatism, and spherical aberration have be obtained, and the values of h are -8.777,-3.657, and -7.132, respectively.The relationships between the aberration coefficients and the boresight pointing offsets have been determined and are shown in Figure 12 and can be expressed as follows:

19
where ΔP t and ΔP c are the boresight pointing offsets caused by tilt and coma, respectively.k t , k c1 , and k c2 are the corresponding pointing offset influence factors; k t = -3.526,k c1 = 3 94, and k c2 = -2.2.
Similarly, the relationships between the normalized aberration coefficients and the sidelobe level increase (ΔSLL) have been determined and are shown in Figure 13, and they     So far, the effects of wavefront errors with the aberration forms have been studied.The arbitrary wavefront error can be expressed by the Zernike expansions, and the Seidel aberration coefficient can be computed by using equation (17) and Table 4. Based on the aberration coefficients, the weight of the effects of different types of aberrations can be determined.Thus, the aberration coefficients can serve as a useful metric to evaluate the wavefront distortion behavior and to determine the primary feature of the influence of distortion.

Conclusion
To understand the effects of arbitrary wavefront errors, it is necessary to investigate the effects of different distribution characteristics of surface distortions on EMP.On the basis of the correlation interval of error distribution, the surface  13 International Journal of Antennas and Propagation distortions are divided into small-scale, medium-scale, and large-scale errors.We proposed that the large-scale errors can be described by the Zernike polynomials.The analysis indicated that distortions in the distribution features of the defocus and the spherical aberration have a greater impact on the gain, and the distribution features of the tilt and the coma mainly influence the boresight offset; the distribution features of the defocus, astigmatism, and spherical aberration have a greater impact on sidelobe levels.Besides, the beam patterns show the similarity between different beam shapes with the same distribution forms of distortions.Based on the Seidel aberration, the relationships between typical EMPs and the aberration coefficients are presented.The gain loss is proportional to the square of the aberration coefficients of the focus, the astigmatism, and the spherical aberration.The boresight pointing offset is proportional to the tilt aberration coefficient, and the relationship between the aberration coefficient and the boresight pointing offset of the coma is a quadratic polynomial.The relationships between the normalized aberration coefficients and the sidelobe level increase are also a quadratic polynomial.We can predict the antenna performance with the help of the aberration coefficients.Therefore, the aberration coefficient can be considered as a metric to evaluate the wavefront distortion behavior and to determine the primary influence feature of distortions.

3
International Journal of Antennas and Propagation is equal to twice the axial component of the normal deviation of the surface at that point and is given by

r n−2l , 6 where
the indices n and m are the radial order and the azimuthal frequency, respectively, m ′ = m , and j is a polynomial number.The values of n and m ′ are always integers and satisfy m ′ ≤ n and n − m ′ = even, and the corresponding pairs of n and m are called Zernike modes.Consequently, only polynomials with certain combinations of n and m exist.The first 21 Zernike polynomials, modes, and distribution features are listed in Table

Figure 5 :
Figure 5: Shape characteristics of the first 21 (excluding the first one) Zernike circle polynomials.

Figure 6 :
Figure 6: 0 °and 90 °cut plane far-field patterns of the OPD distributed in the first 21 normalized Zernike polynomials.The patterns are normalized to the peak of the ideal pattern.

Figure 7 :
Figure 7: Beam contour patterns of the far-field radiation with the OPD distributed in the first 21 normalized Zernike polynomials.The magnitudes of contours range from "peak" to "peak-30 dB." sign of the magnitude of the focus term is chosen to minimize the magnitude of the coefficient, and astigmatism uses the sign opposite that.

Figure 8 :
Figure 8: Shape characteristics of the Seidel wavefront aberration.

Figure 9 :
Figure 9: 0 °and 90 °cut plane far-field patterns of the normalized Seidel wavefront aberration.The patterns are normalized to the peak of the ideal pattern.

Figure 10 :
Figure 10: The beam contour patterns of the normalized Seidel wavefront aberration.The magnitudes of contours range from "peak" to "peak-30 dB."

Figure 13 :
Figure 13: Sidelobe level increase ΔSLL as a function of a/λ.
can also be obtained.By integration in a unit circle, the Zernike RMS polynomials are given by

Table 2 :
The result of effect of errors distributed in the normalized Zernike polynomials on EMP.

Table 3 :
The nine Zernike terms and aberrations.
a ⋅ a′ 2 + ks b ⋅ a′, 20 where ks a and ks b are the sidelobe level influence coefficients, which are different for different aberrations.The values of ks a for tilt, focus, astigmatism, coma, and spherical aberration

Table 4 :
The Seidel aberration in terms of Zernike coefficients.

Table 5 :
The result of effect of the normalized Seidel wavefront aberrations on EMP.