# Investigations on Shaped Mirror Systems in Quasi-Optical Mode Converters Based on Irradiance Moments Method

Authors: Hai Wang; Zejian Lu; Xiaoming Liu; Junsheng Yu; Yuan Yao; Xiaodong Chen; Limei Qi; Zhijiao Chen
Journal: International Journal of Antennas and Propagation
 (2016)

---

## Abstract

A method of transforming high-order Gaussian beams (GBs) mode into circular symmetry fundamental Gaussian beam (FGB) mode with arbitrary waist size is presented using irradiance moments method in quasi-optical (QO) mode converters. The double shaped mirrors correcting amplitude and phase simultaneously are generated by a single incidence irradiance sampling data and known ideal output FGB taking advantage of linear moment matching technique and Fresnel diffraction theory, which can be applied to a wide frequency range especially significant for terahertz band. The numerical coding procedure of creating double correcting mirrors and its fast convergence speed are discussed at 325 GHz. Numerical and experimental comparisons reveal the conclusion that enhancing surface precision and increasing moments order can improve main lobe levels.


---

## Body

## 1. Introduction

Gyrotrons are usually used in high power microwave (HPM) devices and excite higher-order cylindrical cavity modeTEmn or TMmn. However, those modes cannot directly transmit in long distance due to nonlinear polarization and divergence of axial radiation [1–4]. What is worse, the existence of high-order spurious modes leads to the asymmetric, offset-focus, and unstable amplitude and phase [5–9]. Quasi-optical (QO) mode converters with correction mirrors are widely utilized in such system to shape the electromagnetic field into desired fundamental Gaussian beam (FGB) TEM00 field distributions.Traditionally, for effectively suppressing beam divergence and focusing adjustment, the QO mode converters usually adopt regular mirrors with quadratic surfaces to transform high-order modes into low-order modes, such as parabolic and elliptic mirrors. However, in some applications, the quality of the output beam from regular mirrors cannot meet the requirements of pure FGB distributions. For example, although the power-conversion efficiency of 220 GHz Denisov-type QO mode convertor is about95% in [9], the output beam still can be corrected. The conversion efficiency of Vlasov-type QO mode convertor for a 0.42 THz in [10] is 88%. So the conversion efficiency needs to be further improved, and the output beam should be effectively converted into the desired FGB profile by introducing shaped correcting mirrors system, which are often irregular structures created by iterative algorithms. There are many theoretical and experimental investigations on shaped correcting mirrors system [11, 12], such as Katsenelenbaum-Semenov algorithm (KSA) method [13–18], the phase gradient method [19], wave-beam shaping technique [20], and the irradiance moments method [21–24]. The goal of KSA method is to obtain the minimum deviation of beam energies between the known incident field and desired output field passing through each mirror. This algorithm needs to be calculated separately for each mirror, and the iterative result of one mirror directly affects the generation of next mirror. Moreover, the time-consuming and optimization results largely depend on the choices of initial perturbed surface.When it knows the incident electric field and the reflected electric field of each mirror, the phase gradient method based on the Fast Fourier Transform (FFT) can be applied to obtaining the express of correction mirror. The perfect electric conductor (PEC) mirror phase corrector has been designed in the sub-THz gyrotron in [19], where it is suitable for smooth phase front for solving derivative directly on the phase difference. The wave-beam shaping technique needs to calculate the forward-backward-forward field by multiple circulatory iterations taking advantage of plane wave expansion for calculation with FFT. The iterative is completed when the wave-beam profiles converge and its typical iteration number is 50–100 times [20].Irradiance moments method, which is originally utilized for phase retrieval based on amplitude measurement data at several planes [21, 22], has also been successfully demonstrated in mirror synthesis in a 110 GHz gyrotron internal mode converter [24]. The correcting mirrors synthesized using the irradiance moments are easy to fabricate because the surface profiles are smooth when compared with surfaces created by the KSA algorithm in [18], making a huge advantage in comparison to the mirrors synthesized using other iteration methods, where the mirror surfaces have to be smoothed out using the phase unwrapping and filtering procedures [25]. However, There are still some problems that remain to be unsolved when applying the irradiance moments method into terahertz (THz) wave band. The increasing frequency posts strict demands on manufacture of the system, which could lead to more stray high-order modes. In this paper, several important parameters are investigated, such as (a) the number of moments orders adapted in the iterative process, (b) the proportion of high-order GBs mode, (c) the mirror manufacturing error, and (d) the available working bandwidth especially in high frequency.In this paper, the characteristics of irradiance moment method are overall investigated, and the numerical example of a double shaped mirrors system analyzing different mode ratio of the first fourth high-order GBs mode is discussed in THz band. The resolution ability of different moments order is also examined according to the quality of reflection field involving double conversion coefficients and main lobe levels. The comparison results among window fields obtain the effect of surface roughness and suitable bandwidth.This paper is organized as follows. The principle of irradiance moments method and coding procedure of numerical optimization creating double shaped mirrors are viewed in Section2. The double correcting mirrors system and related parameters analysis on the irradiance moments method are detailedly described in Section 3. Finally, we discuss and summarize all of the work.


## 2. Theory and Implementation of Irradiance Moments Method

The location relationship between double correcting mirrors system utilized in QO mode convertor and the conversion process of FGB is depicted in Figure1, where Mirror 1 and Mirror 2 are shaped mirrors to be determined. In this system, the field on the incident plane radiates out into Mirror 1 with an angle θ and then is reflected from Mirror 2 into the window plane. Generally, the field on the incident plane can be known by either direct measurement of the amplitude and phase distributions on the plane or, for those cases where only amplitude distribution can be measured, some phase retrieval methods being available to restore the phase information on the plane via sampling two or more arbitrary measurement planes [26, 27]. The field on the window plane is the desired FGB that should have characteristics of the circular symmetry amplitude, low side lobe levels, and stable phase center.Figure 1
The location relationship between double shaped correcting mirrors and the conversion process of the FGB.The correcting mirror system can be regarded as an optimization module in amplitude and phase. According to the difference of input and output field, it needs to design light path by mirrors for obtaining desired FGB. Therefore, the problem stated in here can be simplified equivalently to the wave path model as shown in Figure2, where |U1| and |U2| are the amplitude distributions while Φ1(r)=arg⁡(U1) and Φ2(r)=arg⁡(U2) are the phase distributions on plane 1 and plane 2, respectively.Figure 2
The schematic diagram of equivalent wave path model of double correcting mirrors system.In the equivalent model, the electromagnetic field propagates from Mirror 1 to Mirror 2. In order to synthesize Mirror 1 and Mirror 2, the phase distributions on plane 1 and plane 2 have to be determined. The fields on the input plane and plane 1 are corresponding to incident and reflected fields from Mirror 1 while the fields on plane 2 and output plane are relative to Mirror 2. The field on the input plane can be obtained by propagating the original field on the incident plane forward. Meanwhile the field on the output plane is calculated via propagating the field on the window plane backward. The amplitude distributions on plane 1 and plane 2 are identical to amplitude distributions on the input plane and output plane. The location points of discrete mirrors surfaces are derived, respectively, according to the phase relationship of incident and reflected field in their own mirror coordinate systems.The irradiance momentMpq(z0) over the input plane and its propagation formula Mpq(z) are characterized with two-dimensional moments [22, 24]:(1)Mpqz0=Mpqxp,yq,z0M00z0=∬-∞+∞xpyqUzr2dxdy∬-∞+∞Uzr2dxdy,(2)Mpqz,ϕ=-izkp+q∬-∞+∞UG∂p+q∂xp∂yqUG∗dxdy,where the amplitude of |Uz|2 is the irradiance intensity at input plane and G(r,z) is the Green function in free space. The subscripts of p and q represent different moments order in the x and y direction. Thus the total moments order number N is p+q. The denominator M00(z0) corresponds to the zero-order moment representing total beam energy, and Mpq(z0) is also called the average value of irradiance energy. For a collimated light beam, the wave equation with Gaussian shape at plane 2 can be expressed as(3)U2r=P2rexp⁡iφ2r,φ2r=∑j=1Nmax∑k=0jϕj-k,kxj-kykLj+k,where P2(r) is the Gaussian amplitude distribution and 2L is the dimension of the plane. The phase distribution φ2(r) is an oblique upper triangular matrix except for zero order and should be iterated out by moment propagation theorem [22]. The set {ϕmn} is polynomial coefficient with dimensionless constants and can be solved out by the moments Mpq. The number of Nmax depends on the same total moments order number. The total number of the phase terms is Nt=Nmax(Nmax+3)/2. When Nmax=3, Nmax=4, and Nmax=5, the phase terms, which are to be discussed in detail, involve Nt=9, Nt=14, and Nt=20 polynomial coefficients, respectively.Using (2), the moments in arbitrary location from plane 2 to plane 1 can be represented. According to the known incidence moment information in (1), the next step is to find a convergence way that makes the differences between Mpq(z) and Mpq(z,ϕ) smallest, which is related to variables set {ϕmn}.When it begins the first iteration, the phase in plane 2 can take a guess that the phase polynomial coefficients{ϕmn} are given the same number as initial iteration condition. By evaluating the derivative of function Mpq(z,ϕ) with respect to every set element {ϕmn} and calculating the differences ΔMpq between the Mpq(z,ϕ) and Mpq(z0), the changes of set {Δϕmn} can be found by solving the matrix equations. What is more, the new set {ϕmn+Δϕmn} can be made up in plane 2 for continuing to start next iteration until set convergence. Due to the fact that the set {ϕmn} is real value, the real number of iteration matrix equation can be expressed as follows[24]:(4)ReBNt×Nt·ΔϕmnNt×1=ReΔMpqNt×1,B=∂Mpqz,ϕ∂ϕmn=-izkp+q∬-∞+∞dxdy·ixmynLm+nUG∂p+qUG∗∂xp∂yq+UG∂p+q-ixmyn/Lm+nUG∗∂xp∂yq.The linear matrix can be solved via Newton iteration method and the iterative convergence criteria are that the square sum of amplitudeΔMpq is close to zero. Once the polynomial coefficient set {ϕmn} is solved, the field on plane 2 can be determined and the field in plane 1 can be obtained by the Rayleigh-Sommerfeld diffraction integral formula in [28](5)U1x1,y1,z1=12π∬-∞+∞U2x2,y2,z2·exp⁡jkrRz1-z2R1R-jkdx2dy2,(6)Ewindowx,y=2πw02exp⁡-x2+y2w02,where R is the integral points distance between double planes. According to (5) and input plane filed, the phase distribution variation passing the first mirror can be obtained. In the same principle, we can obtain the back propagation field Eout reflected from the second mirror by the Rayleigh-Sommerfeld diffraction integral only using the window plane field (6), and the phase variation passing the second mirror can also be obtained. Here, the window field is a FGB, the parameters of which, such as the position of window and the size of beam waist w0, are chosen according to practical applications.When the total four plane fields have been solved, the coordinate position of double correcting mirrors profiles(x′,y′,z′) with a certain incidence angle θ can be created by the difference of phase. However, the field calculated directly is wrapped phase, where the range becomes (-π,+π) and leads to discontinuous surfaces, so it needs to manually add the appropriate 2π·g(r) to become the unwrapped phase with new range (-∞,+∞), where g(r) has to be integer values. An example of phase process of exp⁡(j5x2) transformation wrapped phase into unwrapped phase can be found in [25]. The double correcting unwrapped phase now can be generated according to (7) and (8). The three-dimensional correcting mirrors profiles decided by phase shift can be obtained by (9), where the y-axis is system symmetry axis. The unwrapped process of wrapped phase is coded in Algorithm 1. The corresponding numerical coding procedure of irradiance moments method is shown as in Figure 3:(7)ΔΦIr=Φ1r-Φinr+2π·gIr,(8)ΔΦIIr=Φoutr-Φ2r+2π·gIIr,(9)x′y′z′=1cos⁡θ000100012kcos⁡θxyΔΦ.Algorithm 1:The unwrapped process of wrapped phase precoding.
Input: (1) Wrapped phase matrix [A]Nt×Nt; (2) Calculate the center number Ncenter=INTNt2+1. Initialization: Unwrapped phase matrix P=0Nt×Nt; (1) Update matrix P=A; (2) Phase unwrap in x direction: for  (ii=1; ii≤Nt; ii++) for  (jj=Ncenter; jj≥2; jj--) while  abs(P(ii,jj)-P(ii,jj-1))≥180 if  (abs(P(ii,jj))>abs(P(ii,jj-1))) P(ii,jj-1)←P(ii,jj-1) add 360 degree; else P(ii,jj-1)←P(ii,jj-1) minus 360 degree; end if end while end for for  (jj=Ncenter; jj≤Nt-1; jj++) while  abs(P(ii,jj+1)-P(ii,jj))≥180 if  (abs(P(ii,jj+1))<abs(P(ii,jj))) P(ii,jj+1)←P(ii,jj+1) add 360 degree; else P(ii,jj+1)←P(ii,jj+1) minus 360 degree; end if end while end for end for (3) Phase unwrap in y direction by repeat step (2); Output: Unwrapped phase P.Figure 3
Numerical coding procedure of irradiance moments method creating double correcting mirrors.For comparing the consistency between the real reflected fieldsf1·exp⁡(jϕ1) after every mirror and the theoretical design values f2·exp⁡(jϕ2), both the parameters that scalar and vector conversion efficiency formula use (10) on behalf of beam coupling in energy and direction should be selected. They should be referenced at the same time because their center phase of main lobe may be likely to have significant deterioration although it sometimes has very high value of recovered energy:(10)ηs=∬Sf1·f2dr∬Sf12dr·∬Sf22dr,ηv=∬Sf1f2exp⁡iϕ1-ϕ2dr·∬Sf1f2exp⁡iϕ2-ϕ1dr∬Sf12dr·∬Sf22dr.


## 3. Numerical Results and Comparison

A double correcting mirrors system working at the center frequency of 325 GHz is synthesized in this chapter based on the irradiance moments method. The double correcting mirrors have the function of adjusting amplitude and phase simultaneously. The presented irradiance moments method is implemented in a MATLAB code and the system structure is created in commercial simulation software GRASP.


### 3.1. The Geometry Model of the Correcting Mirrors System

The two and three dimensions’ geometry model of the double correcting mirrors system is depicted in Figure4. The distance between feed and the first mirror is D1=6.23 cm, and the double mirrors have the distance of D2=8.31 cm. The surfaces are defined by local coordinate systems (x1,y1,z1) and (x2,y2,z2) in Figure 4 at the center reflection points specified by D1 and D2 in system, where sign⁡(+) and sign⁡(-) express rotation direction of z-axis. The dimensions of both mirrors are 6.71 cm × 5.81 cm in rectangle shape and there are 14641 discrete coordinate points in total for each mirror, where sampling distance between half and one wavelength is advisable. It is worth remembering that the truncated dimension of the mirror should be larger than beam waist for avoiding edge diffraction as far as possible. The window plane is chosen at a distance of D3=9.84 cm away from Mirror 2. The desired beam waist is w0=6λ, which can be adjusted according to the GB propagation theory.Figure 4
The geometry model of the correcting mirrors system: (a) coordinate system and structural parameters; (b) three-dimensional simulation model(θ=45deg) in GRASP.
(a) (b)The correcting shaped mirrors are synthesized via the procedure presented in Figure3, while the whole system is validated using Physical Optics (PO) method in GRASP. In practice, the radiated field from the radiator can be decomposed into a superposition of various GB modes. Therefore, in order to simulate the actual output field of radiators, a mixture of the fundamental and high-order GBs mode by applying different waist radius is utilized as the feed. The selection of GBs mode combination with different waist radius greatly affects the distortion of the feed from FGB. The zero-order and high-order GBs mode in the rectangular coordinate [29] are shown as(11)Emnr,z=1πωxωy2m+n-1m!n!0.5Hm2xωxHn2yωy·exp⁡-x2ωx2-y2ωy2·exp⁡-jkz+πx2λRx+πy2λRy-2m+1ϕ0x2-2n+1ϕ0y2,where ωx and ωy denote the beam radius, Rx and Ry describe the radius of curvature, and ϕ0x and ϕ0y express the phase shift in the x and y direction, respectively. The function of Hm(x) and Hn(y) represents the Hermite polynomials. The field on the incident plane composes zero-order mode GB(0,0) and high-order modes GB(m,n) with different waist radius on the x and y direction and then is inputted into the GRASP software as the feed in tabulated form.Considering different component content of elliptic high-order GBs mode [see (11)], we select GB(0,0) and the first fourth high-order GBs mode under the ratio of 0.7 : 0.06 : 0.03 : 0.02 : 0.002 as inputted feed field, and the total number of GBs mode is 15 (1 : 2 : 3 : 4 : 5). The inputted zero-order mode GB(0,0) of feed field has ωx=8λ and ωy=10λ in the x and y direction, respectively, but the waists of high-order modes vary from 6λ to 15λ for presenting the asymmetry of the feed field. The desired beam waist on the window is ω0=6λ at the center frequency of 325 GHz. The scalar and vector conversion efficiency between the inputted feed field and the desired window field are only 88.84% and 78.79%, respectively.The presented irradiance moments algorithm converges after seven iterations. Using (7), (8), and (9), the double shaped correcting mirrors created in θ=45 degrees by irradiance moments method are shown in Figure 5. The feed field is transformed into a field specified in a plane to equivalent currents in the same plane as tabulated planar source in GRASP. These currents are integrated when the radiated field is to be determined. The system surfaces are defined through interpolation between tabulated z-values given at irregularly spaced (x,y)-points from discrete coordinate points. The actual beam transmission and transformation passing through double correcting mirrors have been shown in Figure 6, where the left is normalized amplitude (dB) and the right is its phase comparing theoretical field with simulation values. In Figures 6(a)-6(b), the feed field is made up of 15 GBs mode with tilted amplitude and asymmetrical phase. Figures 6(c)-6(d) are reflection field passing through distance of D2, where the beam waist is 10λ and it has linear phase in accordance with theoretical change tendency. The window field is obtained through D3 in Figures 6(e)-6(f), where the beam waist is 6λ, and the phase in main lobe range (1 cm × 1 cm) is circular symmetrical like the change of GB amplitude distribution.Figure 5
The double shaped correcting mirrors are created by irradiance moments method at 325 GHz withθ=45 deg and GBs ratio 0.7 : 0.06 : 0.03 : 0.02 : 0.002. (a) is the first correcting mirror. (b) is the second correcting mirror.
(a) (b)Figure 6
The actual beam transmission and transformation process passing through double correcting mirrors. The left is normalized amplitude (dB) and the right is phase with theoretical (red line) and simulation (blue line) values. (a)-(b) are the feed field, (c)-(d) are the reflected field from the first mirror afterD2, and (e)-(f) are the window field (ω0=6λ) reflected from the second mirror after D3.
(a)
Amplitude (b)
Phase (c)
Amplitude (d)
Phase (e)
Amplitude (f)
PhaseIn the process of validation in GRASP, the discrete coordinate points are inputted into the software to construct the mirror surfaces. However, if the double mirrors are mounted according to the original coordinate system in Figure4, the output field on the window is found to be slightly off-center in the simulation, which leads to a loss of conversion efficiency. This phenomenon is caused by discretization mirror surfaces, where they do not completely satisfy the law of reflection for irregular shape. However, the problem can be solved by slightly adjusting the orientation of double mirrors in the opposite direction. According to Figure 4(a), selecting mirror coordinate axes z1 and z2 as the rotation axis is to change their directions meanwhile guaranteeing that the y direction is constant. Obviously, when the zero-order mode GB(0,0) becomes lower, the perturbed surface causes the greater rotation angle. For example, when the ratio of GB(0,0) is 0.7, it needs to adjust the angle (Δθz1,Δθz2)=(+0.2997°,-0.293°), respectively, but selects the angle (Δθz1,Δθz2)=(+0.6426°,-0.5807°) at the case of 0.56.


### 3.2. The Influence of Incident Angles on the Window Field

In (9), the final discrete data points of mirrors surface depend on the incident angle θ. Considering the blockage effect between double mirrors and in order to fully receive the field energy, the incident angle θ can be selected reasonably in the range of 30~50 degrees. Comparison simulation results of window field between scalar and vector conversion efficiency with different input angles and the position of the first side lobe about normalized copolarization amplitude are shown in Table 1. There are minor differences between such results and they could be caused by simulation error, so the incident angles cannot largely change the performance of the window field. In this case, unless the incident angles are otherwise expressed, all the simulation results are obtained in the case of 45 degrees.Table 1
The influence of incident angles on the output window field.
Input angle η s   D 3 η v   D 3 Side lobe levels 30 deg 0.989361 0.941582 −42.12 dB 35 deg 0.990832 0.945805 −42.31 dB 45 deg 0.990953 0.948277 −42.82 dB 50 deg 0.987739 0.938693 −41.03 dB


### 3.3. The Convergence Speed of the Irradiance Moments Method

The convergence speed of the irradiance moments depends on the proportion of different GBs mode and the order of moments chosen. The relationship among the iteration times, the convergence speed, and the computation time with the different component of GBs mode and the order of moments is shown in Figure7. In the picture, the horizontal axis represents the iteration times and the left vertical axis is the convergence criteria while the right vertical axis corresponds to the computation time in hours. In both pictures, the solid lines are mapped based on the number of high-order GBs mode (m+n is less than 2, 3, 4, and 5 where it is marked as 2M, 3M, 4M, and 5M), the ratio of elliptic zero-order mode GB(0,0) in the mixed fields (marked as 0.40, 0.56, and 0.70), and the moments order (marked as 3rd, 4th, and 5th); meanwhile the dotted line having same color and marker symbol is its corresponding iteration time (hours).Figure 7
The relationship among the iteration times, the convergence speed, and the computation time. The horizontal axis is iteration times. In vertical axis, the left is convergence speed(log10(ΔM)), and the right is iteration time (hours): (a) keeps same mode number but has different ratio of GB(0,0) and moments order number; (b) keeps same ratio of GB(0,0) with 0.7 but increases the number of high-order GBs mode and moments order number.
(a) (b)As it is shown in the picture, the iteration times increase when the order of moments is improving. The algorithm can always converge for the first fourth-order GBs mode within 6~14 iteration times. With the same moments order, the lower the mode component of GB(0,0) is, the slower the convergence speed it has. However, when the component of GB(0,0) mode is only 40 percent, the fifth order of moment does not converge because it becomes difficult to restore the accurate phase distribution of fundamental mode when its incident content is very low. From the black line with square maker in Figure 7(a), the value of ΔM rises from the first time and then it begins to float in 104~108 all the time during 12 hours with 9 times. This also shows that the increasing number of moments order does no good to the lower GB(0,0) situation. If the input field ensures the same ratio of GB(0,0) with 0.7 and meanwhile increases the number of high-order modes, the convergence speed becomes faster. Obviously, the third-order moment has the fastest convergence speed, and the convergence speed of the fifth-order moment is 7~9 times that of the fourth order.The required computational time is calculated with121×121 complex data using CPU i7-2600@3.4 GHz and RAM 12 GB in 64-bit operating system.


### 3.4. The Convergence Speed of the Irradiance Moments Method

The comparison results of window field between amplitude and phase under the different moments order are shown in Figure8. The incident field is the same in Section 3.1 (Figures 6(a)-6(b)). The three different mirrors systems are generated by the 3rd-, 4th-, and 5th-order moments method and are required to obtain same window field ω0=6λ, which are plotted at the one picture by shifting 4 cm and 8 cm. The three systems have a good amplitude conversion capability as low as −46 dB and centrically circumferential phase distribution features. In the main lobe range, all phases distribution is consistent with the amplitudes shape. The phase results show that the 5th-order moments have a significant transition effect for possessing better side lobe profile. More copolarization comparison results of window field with different GBs ratio and moment orders are shown in Table 2. In this table, when the content of GB(0,0) is greater than 0.7, selecting the 5th-order moment becomes best relative to the 4th and 3rd in improving the side lobe levels. However, when the GB(0,0) is less than 0.4, the 5th-order moment should consider its nonconvergent character. The 4th is close to 3rd orders on emergent field conversion performance. From the practical conversion efficiency, it is obvious that 4th-order moment method is the most acceptable and efficient.Table 2
The influence of different GBs ratio and moments orders on the window field.
GBs ratio Moment order η s   D 3 η v   D 3 Side lobe levels 0.7 : 0.06 : 0.03 : 0.02 : 0.002 3rd 0.979397 0.927940 −42.50 dB 4th 0.987519 0.938268 −42.86 dB 5th 0.989743 0.940147 −53.01 dB  0.56 : 0.07 : 0.05 : 0.035 : 0.002 3rd 0.977260 0.863493 −33.89 dB 4th 0.977479 0.864714 −34.87 dB 5th 0.974416 0.864449 −34.61 dB  0.4 : 0.095 : 0.06 : 0.055 : 0.002 3rd 0.943856 0.648910 −19.54 dB 4th 0.945071 0.655087 −19.77 dB 5th NonconvergentFigure 8
The copolarization comparison results of window field (ω0=6λ) between amplitude (dB) and phase (degree) under the different moments order number (3rd, 4th, and 5th) but with same incident field.


### 3.5. The Resolution of High-Order GBs Mode

With the mixture of high-order GBs mode, the amplitude and phase of incident field are asymmetric and become obvious tilt. Those phenomena are caused by mirrors diffraction and location mismatching from superior radiators. The mode levels of high-order GBs have different effect on the recovery of window field FGB. Selecting different mode content of the first 15 GBs mode is to consider the solving resolution of 4th moment method. The results of scalar and vector conversion efficiency on the positionP and position W are shown again in Table 3. From the table, it can be found that when the mode level of elliptic zero-order GB(0,0) content is lower, the side lobe levels on the window position will rise. It can be deduced from the results that its energy and direction of main lobe after the second correction are further improved relative to the first mirror results. The improvement is obvious especially for vector conversion, and the highest value can increase by 28% in 0.56 content. However, when the main mode is less than 40 percent, the amplitude and phase correction are not easy to be achieved for having very high side lobe levels.Table 3
The transformation results of the first fourth-order GBs mode with different percentages.
GBs ratio η s   D 2 η v  D2 η s   D 3 η v   D 3 Side lobe levels 1 : 0 : 0 : 0 : 0 0.9945 0.9885 0.9945 0.9978 −83.45 dB 0.78 : 0.055 : 0.02 : 0.01 : 0.002 0.9901 0.9722 0.9905 0.9797 −56.84 dB 0.7 : 0.06 : 0.03 : 0.02 : 0.002 0.9848 0.8953 0.9909 0.9458 −42.86 dB 0.56 : 0.07 : 0.05 : 0.035 : 0.002 0.9544 0.6207 0.9775 0.8647 −34.87 dB 0.4 : 0.095 : 0.06 : 0.055 : 0.002 0.8303 0.3950 0.9451 0.6551 −19.77 dBFurthermore, if it only changes the second correcting mirror, a number of FGB having different waist radius size can be obtained. When the beam waist radius atP location and the window position are known, according to the GB propagation theory, those beam propagation distances can be derived. The beam propagation distances 9.84 cm, 7.52 cm, 5.86 cm, and 4.01 cm correspond to beam waists of 6λ, 4λ, 3λ, and 2λ in the window plane, respectively. Those propagation distances have been changed down to 0.7 times comparing theoretical propagation values for making structure compact. The results converting different waist radius at 325 GHz with regard to the first fourth high-order GBs mode with some ratio are shown in Figure 9. The copolarization results in phi = 0 deg and phi = 90 deg are drawn together. The simulation results show the good consistency in 6λ, 4λ, and 3λ and can be up to −43 dB. The dissymmetry result of 2λ in phi = 90 deg is caused by the very short output distance and the fact that that the focal length of second elliptical concave mirror decreases.Figure 9
Window field is FGB with different waist radius6λ, 4λ, 3λ, and 2λ at 325 GHz (λ=0.0923cm).


### 3.6. The Influence of Surface Roughness

With the increasing of frequency, the smaller the manufacture error of the discrete mirror is, the more accurate the reflection is because mirror surfaces are created by the difference of input and reflection phase distribution. However, the actual manufacturing error is inevitable depending on the equipment accuracy, material thicknesses, or surface flatness although the discretization data points can be selected as much as possible.In order to simulate the actual surface roughness, the surface data is determined by a cubic interpolation procedure at intermediate points. Choosing three groups error parameters with different peak value atλ/100=0.000923 cm, λ/50=0.001846 cm, and λ/25=0.003692 cm, the mirror surface randomly floats between zero and peak value on the basis of theoretical values. Comparing such window field results with the case of zero error value, the results are shown in Figure 10, where the input field is also like in Section 3.1 (Figures 6(a)-6(b)).Figure 10
When surfaces add peak error withλ/100, λ/50, and λ/25, the copolarization and cross polarization results of window field compared with zero error. (a) is the result in phi = 0 deg; (b) is the result in phi = 90 deg.
(a) (b)From the simulation results in Figure10, window field can be found that the beam width and the levels of cross polar are consistent at the center beam range, but the effect on the side lobe levels is obvious and different. Without mirror error, the side lobe levels of copolarization can be lower than −42.86 dB. When the mirror error rises every two times, the side lobe levels of copolarization and cross polarization increase by every 5 dB away from the main lobe at the same time. Due to its phase being chaotic near the side lobe, the rising of side lobe levels makes the width of main lobe further decreased and leads to reducing the central range of stable phase.


### 3.7. The Operating Bandwidth for the Same Correcting Mirrors

For considering suitable operating bandwidth of mirrors in irradiance moments method, a number of feeds working at the different frequencies are provided to radiate the same mirrors system. The same double mirrors are created at the frequency of 325 GHz. The incident feed field also makes up the first fourth high-order GBs mode with same parameters in Section3.1 (Figures 6(a)-6(b)). All the window fields in same lengths are obtained passing through the same distance D3=9.84 cm. The results of the co- and cx-polarization in phi = 0 and phi = 90 deg are shown in Figure 11. More details comparing output window field among different frequencies are listed in Table 4, where all results are compared relative to 325 GHz. According to the simulation results, the change of frequency will not directly influence work performance of mirrors surfaces. Its working frequency can be as low as 200 GHz and up to 700 GHz having some beam size in the window position. When deviating from the center frequency, both scalar and vector conversion efficiencies are reduced. However, the cross polarization isolation will rise when the frequency increases because the decreasing size of feed leads to the more energy focusing on the center of mirror. When it works below 200 GHz, the cross polarization isolation is significantly worse down to −20 dB. Although it has better side lobe levels above 700 GHz, its phase in window plane is distorted seriously for vector conversion efficiency down to 91 percent. Therefore, the corrective mirrors by irradiance moments method can be applied to a wide range of frequency especially significant for THz band.Table 4
The window field performance comparison using same mirrors working at different frequency.
Frequency (GHz) η s   D 3 η v   D 3 Cross polarization isolation (dB) 200 0.990371 0.969840 27.82 250 0.995544 0.985844 29.13 300 0.997027 0.987443 30.33 315 0.997192 0.988120 30.65 320 0.997230 0.988909 30.75 325 1 1 31.06 330 0.997386 0.988345 31.15 350 0.997234 0.986745 31.34 450 0.996939 0.972275 32.98 500 0.996585 0.962383 33.61 700 0.995150 0.912628 35.31Figure 11
A number of feeds working at the different frequencies are provided to radiate the same mirrors system created at 325 GHz. Window field is obtained passing through the same distanceD3. (a) is the copolarization and cx-polarization results in phi = 0 deg; (b) is the results in phi = 90 deg.
(a) (b)


## 3.1. The Geometry Model of the Correcting Mirrors System

The two and three dimensions’ geometry model of the double correcting mirrors system is depicted in Figure4. The distance between feed and the first mirror is D1=6.23 cm, and the double mirrors have the distance of D2=8.31 cm. The surfaces are defined by local coordinate systems (x1,y1,z1) and (x2,y2,z2) in Figure 4 at the center reflection points specified by D1 and D2 in system, where sign⁡(+) and sign⁡(-) express rotation direction of z-axis. The dimensions of both mirrors are 6.71 cm × 5.81 cm in rectangle shape and there are 14641 discrete coordinate points in total for each mirror, where sampling distance between half and one wavelength is advisable. It is worth remembering that the truncated dimension of the mirror should be larger than beam waist for avoiding edge diffraction as far as possible. The window plane is chosen at a distance of D3=9.84 cm away from Mirror 2. The desired beam waist is w0=6λ, which can be adjusted according to the GB propagation theory.Figure 4
The geometry model of the correcting mirrors system: (a) coordinate system and structural parameters; (b) three-dimensional simulation model(θ=45deg) in GRASP.
(a) (b)The correcting shaped mirrors are synthesized via the procedure presented in Figure3, while the whole system is validated using Physical Optics (PO) method in GRASP. In practice, the radiated field from the radiator can be decomposed into a superposition of various GB modes. Therefore, in order to simulate the actual output field of radiators, a mixture of the fundamental and high-order GBs mode by applying different waist radius is utilized as the feed. The selection of GBs mode combination with different waist radius greatly affects the distortion of the feed from FGB. The zero-order and high-order GBs mode in the rectangular coordinate [29] are shown as(11)Emnr,z=1πωxωy2m+n-1m!n!0.5Hm2xωxHn2yωy·exp⁡-x2ωx2-y2ωy2·exp⁡-jkz+πx2λRx+πy2λRy-2m+1ϕ0x2-2n+1ϕ0y2,where ωx and ωy denote the beam radius, Rx and Ry describe the radius of curvature, and ϕ0x and ϕ0y express the phase shift in the x and y direction, respectively. The function of Hm(x) and Hn(y) represents the Hermite polynomials. The field on the incident plane composes zero-order mode GB(0,0) and high-order modes GB(m,n) with different waist radius on the x and y direction and then is inputted into the GRASP software as the feed in tabulated form.Considering different component content of elliptic high-order GBs mode [see (11)], we select GB(0,0) and the first fourth high-order GBs mode under the ratio of 0.7 : 0.06 : 0.03 : 0.02 : 0.002 as inputted feed field, and the total number of GBs mode is 15 (1 : 2 : 3 : 4 : 5). The inputted zero-order mode GB(0,0) of feed field has ωx=8λ and ωy=10λ in the x and y direction, respectively, but the waists of high-order modes vary from 6λ to 15λ for presenting the asymmetry of the feed field. The desired beam waist on the window is ω0=6λ at the center frequency of 325 GHz. The scalar and vector conversion efficiency between the inputted feed field and the desired window field are only 88.84% and 78.79%, respectively.The presented irradiance moments algorithm converges after seven iterations. Using (7), (8), and (9), the double shaped correcting mirrors created in θ=45 degrees by irradiance moments method are shown in Figure 5. The feed field is transformed into a field specified in a plane to equivalent currents in the same plane as tabulated planar source in GRASP. These currents are integrated when the radiated field is to be determined. The system surfaces are defined through interpolation between tabulated z-values given at irregularly spaced (x,y)-points from discrete coordinate points. The actual beam transmission and transformation passing through double correcting mirrors have been shown in Figure 6, where the left is normalized amplitude (dB) and the right is its phase comparing theoretical field with simulation values. In Figures 6(a)-6(b), the feed field is made up of 15 GBs mode with tilted amplitude and asymmetrical phase. Figures 6(c)-6(d) are reflection field passing through distance of D2, where the beam waist is 10λ and it has linear phase in accordance with theoretical change tendency. The window field is obtained through D3 in Figures 6(e)-6(f), where the beam waist is 6λ, and the phase in main lobe range (1 cm × 1 cm) is circular symmetrical like the change of GB amplitude distribution.Figure 5
The double shaped correcting mirrors are created by irradiance moments method at 325 GHz withθ=45 deg and GBs ratio 0.7 : 0.06 : 0.03 : 0.02 : 0.002. (a) is the first correcting mirror. (b) is the second correcting mirror.
(a) (b)Figure 6
The actual beam transmission and transformation process passing through double correcting mirrors. The left is normalized amplitude (dB) and the right is phase with theoretical (red line) and simulation (blue line) values. (a)-(b) are the feed field, (c)-(d) are the reflected field from the first mirror afterD2, and (e)-(f) are the window field (ω0=6λ) reflected from the second mirror after D3.
(a)
Amplitude (b)
Phase (c)
Amplitude (d)
Phase (e)
Amplitude (f)
PhaseIn the process of validation in GRASP, the discrete coordinate points are inputted into the software to construct the mirror surfaces. However, if the double mirrors are mounted according to the original coordinate system in Figure4, the output field on the window is found to be slightly off-center in the simulation, which leads to a loss of conversion efficiency. This phenomenon is caused by discretization mirror surfaces, where they do not completely satisfy the law of reflection for irregular shape. However, the problem can be solved by slightly adjusting the orientation of double mirrors in the opposite direction. According to Figure 4(a), selecting mirror coordinate axes z1 and z2 as the rotation axis is to change their directions meanwhile guaranteeing that the y direction is constant. Obviously, when the zero-order mode GB(0,0) becomes lower, the perturbed surface causes the greater rotation angle. For example, when the ratio of GB(0,0) is 0.7, it needs to adjust the angle (Δθz1,Δθz2)=(+0.2997°,-0.293°), respectively, but selects the angle (Δθz1,Δθz2)=(+0.6426°,-0.5807°) at the case of 0.56.


## 3.2. The Influence of Incident Angles on the Window Field

In (9), the final discrete data points of mirrors surface depend on the incident angle θ. Considering the blockage effect between double mirrors and in order to fully receive the field energy, the incident angle θ can be selected reasonably in the range of 30~50 degrees. Comparison simulation results of window field between scalar and vector conversion efficiency with different input angles and the position of the first side lobe about normalized copolarization amplitude are shown in Table 1. There are minor differences between such results and they could be caused by simulation error, so the incident angles cannot largely change the performance of the window field. In this case, unless the incident angles are otherwise expressed, all the simulation results are obtained in the case of 45 degrees.Table 1
The influence of incident angles on the output window field.
Input angle η s   D 3 η v   D 3 Side lobe levels 30 deg 0.989361 0.941582 −42.12 dB 35 deg 0.990832 0.945805 −42.31 dB 45 deg 0.990953 0.948277 −42.82 dB 50 deg 0.987739 0.938693 −41.03 dB


## 3.3. The Convergence Speed of the Irradiance Moments Method

The convergence speed of the irradiance moments depends on the proportion of different GBs mode and the order of moments chosen. The relationship among the iteration times, the convergence speed, and the computation time with the different component of GBs mode and the order of moments is shown in Figure7. In the picture, the horizontal axis represents the iteration times and the left vertical axis is the convergence criteria while the right vertical axis corresponds to the computation time in hours. In both pictures, the solid lines are mapped based on the number of high-order GBs mode (m+n is less than 2, 3, 4, and 5 where it is marked as 2M, 3M, 4M, and 5M), the ratio of elliptic zero-order mode GB(0,0) in the mixed fields (marked as 0.40, 0.56, and 0.70), and the moments order (marked as 3rd, 4th, and 5th); meanwhile the dotted line having same color and marker symbol is its corresponding iteration time (hours).Figure 7
The relationship among the iteration times, the convergence speed, and the computation time. The horizontal axis is iteration times. In vertical axis, the left is convergence speed(log10(ΔM)), and the right is iteration time (hours): (a) keeps same mode number but has different ratio of GB(0,0) and moments order number; (b) keeps same ratio of GB(0,0) with 0.7 but increases the number of high-order GBs mode and moments order number.
(a) (b)As it is shown in the picture, the iteration times increase when the order of moments is improving. The algorithm can always converge for the first fourth-order GBs mode within 6~14 iteration times. With the same moments order, the lower the mode component of GB(0,0) is, the slower the convergence speed it has. However, when the component of GB(0,0) mode is only 40 percent, the fifth order of moment does not converge because it becomes difficult to restore the accurate phase distribution of fundamental mode when its incident content is very low. From the black line with square maker in Figure 7(a), the value of ΔM rises from the first time and then it begins to float in 104~108 all the time during 12 hours with 9 times. This also shows that the increasing number of moments order does no good to the lower GB(0,0) situation. If the input field ensures the same ratio of GB(0,0) with 0.7 and meanwhile increases the number of high-order modes, the convergence speed becomes faster. Obviously, the third-order moment has the fastest convergence speed, and the convergence speed of the fifth-order moment is 7~9 times that of the fourth order.The required computational time is calculated with121×121 complex data using CPU i7-2600@3.4 GHz and RAM 12 GB in 64-bit operating system.


## 3.4. The Convergence Speed of the Irradiance Moments Method

The comparison results of window field between amplitude and phase under the different moments order are shown in Figure8. The incident field is the same in Section 3.1 (Figures 6(a)-6(b)). The three different mirrors systems are generated by the 3rd-, 4th-, and 5th-order moments method and are required to obtain same window field ω0=6λ, which are plotted at the one picture by shifting 4 cm and 8 cm. The three systems have a good amplitude conversion capability as low as −46 dB and centrically circumferential phase distribution features. In the main lobe range, all phases distribution is consistent with the amplitudes shape. The phase results show that the 5th-order moments have a significant transition effect for possessing better side lobe profile. More copolarization comparison results of window field with different GBs ratio and moment orders are shown in Table 2. In this table, when the content of GB(0,0) is greater than 0.7, selecting the 5th-order moment becomes best relative to the 4th and 3rd in improving the side lobe levels. However, when the GB(0,0) is less than 0.4, the 5th-order moment should consider its nonconvergent character. The 4th is close to 3rd orders on emergent field conversion performance. From the practical conversion efficiency, it is obvious that 4th-order moment method is the most acceptable and efficient.Table 2
The influence of different GBs ratio and moments orders on the window field.
GBs ratio Moment order η s   D 3 η v   D 3 Side lobe levels 0.7 : 0.06 : 0.03 : 0.02 : 0.002 3rd 0.979397 0.927940 −42.50 dB 4th 0.987519 0.938268 −42.86 dB 5th 0.989743 0.940147 −53.01 dB  0.56 : 0.07 : 0.05 : 0.035 : 0.002 3rd 0.977260 0.863493 −33.89 dB 4th 0.977479 0.864714 −34.87 dB 5th 0.974416 0.864449 −34.61 dB  0.4 : 0.095 : 0.06 : 0.055 : 0.002 3rd 0.943856 0.648910 −19.54 dB 4th 0.945071 0.655087 −19.77 dB 5th NonconvergentFigure 8
The copolarization comparison results of window field (ω0=6λ) between amplitude (dB) and phase (degree) under the different moments order number (3rd, 4th, and 5th) but with same incident field.


## 3.5. The Resolution of High-Order GBs Mode

With the mixture of high-order GBs mode, the amplitude and phase of incident field are asymmetric and become obvious tilt. Those phenomena are caused by mirrors diffraction and location mismatching from superior radiators. The mode levels of high-order GBs have different effect on the recovery of window field FGB. Selecting different mode content of the first 15 GBs mode is to consider the solving resolution of 4th moment method. The results of scalar and vector conversion efficiency on the positionP and position W are shown again in Table 3. From the table, it can be found that when the mode level of elliptic zero-order GB(0,0) content is lower, the side lobe levels on the window position will rise. It can be deduced from the results that its energy and direction of main lobe after the second correction are further improved relative to the first mirror results. The improvement is obvious especially for vector conversion, and the highest value can increase by 28% in 0.56 content. However, when the main mode is less than 40 percent, the amplitude and phase correction are not easy to be achieved for having very high side lobe levels.Table 3
The transformation results of the first fourth-order GBs mode with different percentages.
GBs ratio η s   D 2 η v  D2 η s   D 3 η v   D 3 Side lobe levels 1 : 0 : 0 : 0 : 0 0.9945 0.9885 0.9945 0.9978 −83.45 dB 0.78 : 0.055 : 0.02 : 0.01 : 0.002 0.9901 0.9722 0.9905 0.9797 −56.84 dB 0.7 : 0.06 : 0.03 : 0.02 : 0.002 0.9848 0.8953 0.9909 0.9458 −42.86 dB 0.56 : 0.07 : 0.05 : 0.035 : 0.002 0.9544 0.6207 0.9775 0.8647 −34.87 dB 0.4 : 0.095 : 0.06 : 0.055 : 0.002 0.8303 0.3950 0.9451 0.6551 −19.77 dBFurthermore, if it only changes the second correcting mirror, a number of FGB having different waist radius size can be obtained. When the beam waist radius atP location and the window position are known, according to the GB propagation theory, those beam propagation distances can be derived. The beam propagation distances 9.84 cm, 7.52 cm, 5.86 cm, and 4.01 cm correspond to beam waists of 6λ, 4λ, 3λ, and 2λ in the window plane, respectively. Those propagation distances have been changed down to 0.7 times comparing theoretical propagation values for making structure compact. The results converting different waist radius at 325 GHz with regard to the first fourth high-order GBs mode with some ratio are shown in Figure 9. The copolarization results in phi = 0 deg and phi = 90 deg are drawn together. The simulation results show the good consistency in 6λ, 4λ, and 3λ and can be up to −43 dB. The dissymmetry result of 2λ in phi = 90 deg is caused by the very short output distance and the fact that that the focal length of second elliptical concave mirror decreases.Figure 9
Window field is FGB with different waist radius6λ, 4λ, 3λ, and 2λ at 325 GHz (λ=0.0923cm).


## 3.6. The Influence of Surface Roughness

With the increasing of frequency, the smaller the manufacture error of the discrete mirror is, the more accurate the reflection is because mirror surfaces are created by the difference of input and reflection phase distribution. However, the actual manufacturing error is inevitable depending on the equipment accuracy, material thicknesses, or surface flatness although the discretization data points can be selected as much as possible.In order to simulate the actual surface roughness, the surface data is determined by a cubic interpolation procedure at intermediate points. Choosing three groups error parameters with different peak value atλ/100=0.000923 cm, λ/50=0.001846 cm, and λ/25=0.003692 cm, the mirror surface randomly floats between zero and peak value on the basis of theoretical values. Comparing such window field results with the case of zero error value, the results are shown in Figure 10, where the input field is also like in Section 3.1 (Figures 6(a)-6(b)).Figure 10
When surfaces add peak error withλ/100, λ/50, and λ/25, the copolarization and cross polarization results of window field compared with zero error. (a) is the result in phi = 0 deg; (b) is the result in phi = 90 deg.
(a) (b)From the simulation results in Figure10, window field can be found that the beam width and the levels of cross polar are consistent at the center beam range, but the effect on the side lobe levels is obvious and different. Without mirror error, the side lobe levels of copolarization can be lower than −42.86 dB. When the mirror error rises every two times, the side lobe levels of copolarization and cross polarization increase by every 5 dB away from the main lobe at the same time. Due to its phase being chaotic near the side lobe, the rising of side lobe levels makes the width of main lobe further decreased and leads to reducing the central range of stable phase.


## 3.7. The Operating Bandwidth for the Same Correcting Mirrors

For considering suitable operating bandwidth of mirrors in irradiance moments method, a number of feeds working at the different frequencies are provided to radiate the same mirrors system. The same double mirrors are created at the frequency of 325 GHz. The incident feed field also makes up the first fourth high-order GBs mode with same parameters in Section3.1 (Figures 6(a)-6(b)). All the window fields in same lengths are obtained passing through the same distance D3=9.84 cm. The results of the co- and cx-polarization in phi = 0 and phi = 90 deg are shown in Figure 11. More details comparing output window field among different frequencies are listed in Table 4, where all results are compared relative to 325 GHz. According to the simulation results, the change of frequency will not directly influence work performance of mirrors surfaces. Its working frequency can be as low as 200 GHz and up to 700 GHz having some beam size in the window position. When deviating from the center frequency, both scalar and vector conversion efficiencies are reduced. However, the cross polarization isolation will rise when the frequency increases because the decreasing size of feed leads to the more energy focusing on the center of mirror. When it works below 200 GHz, the cross polarization isolation is significantly worse down to −20 dB. Although it has better side lobe levels above 700 GHz, its phase in window plane is distorted seriously for vector conversion efficiency down to 91 percent. Therefore, the corrective mirrors by irradiance moments method can be applied to a wide range of frequency especially significant for THz band.Table 4
The window field performance comparison using same mirrors working at different frequency.
Frequency (GHz) η s   D 3 η v   D 3 Cross polarization isolation (dB) 200 0.990371 0.969840 27.82 250 0.995544 0.985844 29.13 300 0.997027 0.987443 30.33 315 0.997192 0.988120 30.65 320 0.997230 0.988909 30.75 325 1 1 31.06 330 0.997386 0.988345 31.15 350 0.997234 0.986745 31.34 450 0.996939 0.972275 32.98 500 0.996585 0.962383 33.61 700 0.995150 0.912628 35.31Figure 11
A number of feeds working at the different frequencies are provided to radiate the same mirrors system created at 325 GHz. Window field is obtained passing through the same distanceD3. (a) is the copolarization and cx-polarization results in phi = 0 deg; (b) is the results in phi = 90 deg.
(a) (b)


## 4. Discussion and Conclusions

We discuss the correcting mirror system where the incident field is a mixture of high-order GBs mode, and the window field is FGB with different waist radius working at 325 GHz based on the irradiance moments method. Making full use of Rayleigh-Sommerfeld diffraction integral and unwrapped phase, this method is successfully applied to creating double correcting mirrors, where the first mirror produces linear phase and the second can obtain the desired waist radius size.Comparing correcting systems with different incidence angles, it has the best conversion efficiency and main lobe levels when selecting 45 degrees. A double shaped correcting mirrors system dealing with the first fourth high-order GBs and the actual beam transformation process are given in 45 degrees. For solving the problem of center optical axis offsetting the center of window plane at the situation of lowerGB(0,0) ratio, it also needs to adjust the orientation of double mirrors coordinate system in the opposite direction.The convergence speed of algorithm has also been discussed at different moments orders and GBs mode. According to the different simulation results,ΔM can converge to as low as 10-32 after only 6~14 iterations times so the convergence speed is very fast. When it contains more high-order GBs mode and reduces the ratio of GB(0,0), the convergence speed becomes slower. Meanwhile, if the incident field has the same ratio of high-order GBs mode, the simulation time of the fifth-order moment is 7~9 times with respect to the fourth-order moment method.The higher-order moment is significant towards promotion of the performance of the output field. Comparing 5th-order moments with 3rd and 4th, when the content of main modeGB(0,0) is greater than 0.7, selecting the 5th-order moment becomes better in improving both phase distribution and main lobe levels. However, when the content of main mode is very low, it also should consider whether 5th is convergence or not. Although the expression of the field might be more precise using 5th order, the iteration times will increase significantly, which has shown that the convergence becomes more difficult. The 4th is close to 3rd on the conversion results of window field. From practical iterative results and computational efficiency, using 4th-order moment method is more acceptable and its recovery ability can be as low as 40 percent of GB(0,0). What is more, when it needs to acquire different waist radius of window field, it only adjusts the second mirror. The window position and propagation distance can be obtained by GB propagation theory.The surface precision of the practical mirror manufacture needs to be considered importantly. When the mirror error rises every two times, the side lobe levels of copolarization and cross polarization increase by 5 dB away from the main lobe at the same time. Of course, the corrective mirrors can be applied to a wide frequency range from 200 GHz to 700 GHz so it is especially significant for THz band, where low frequency will directly affect main lobe levels and high frequency mainly influences cross polarization. So this QO double correcting mirrors system by irradiance moments method deserves further development.



---

*Source: 8043042-2016-04-07.xml*