Theoretical and Experimental Comparison Results of Dual-Channel 3 D Quasi-Optical Network System between Frequency Selective Surface and Wire Grid Polarizer

BUPT-QMUL EM Theory and Application International Research Lab, Beijing University of Posts and Telecommunications, Beijing 100876, China College of Physics and Electronic Information, Anhui Normal University, Wuhu, Anhui Province 241002, China China Academy of Electronics and Information Technology, Beijing 100041, China School of Electronic Engineering and Computer Science, Queen Mary University of London, London E1 4NS, UK


Introduction
In recent years, along with the rapid development of millimeter wave (MMW) and submillimeter wave technology, the QON system with the features of multifrequency and multipolarization is widely applied in the radio astronomy, meteorological satellites, remote sensing, and imaging system [1].For example, the QO demultiplexer in MASTER has three FSS and one polarizer to perform the separation of five bands for measuring the atmospheric chemistry objectives [2].The MWHS in FY-3A meteorological satellite provides vertical information about water vapor, which has five frequency and dual polarizations at 150 and 183.31 GHz to monitor natural disasters [3,4].A wire grid in the radioastronomy interferometer of ALMA Band 10 has been used to separate the two linear polarizations fields after two elliptical mirrors into corrugated horns, where the cross-polarization is lower than −28 dB at all measured frequencies [5].The channels of demultiplexing and filtering are also performed using FSS and polarizing grid in the MARSCHALS with three center bands at 300, 325, and 345 GHz when the incident angle has 45 ° [6] and so on.On the whole, their development mainly benefits from the progress of optic components and the application of Gaussian beam (GB) model analysis method [7,8].
In the MMW and sub-MMW band, the atmosphere has the absorption characteristic when the electromagnetic wave propagates forward.The QON system often works at the peak absorption frequencies to measure those compositions on the microwave radiometer, such as the water vapor absorption lines have some center frequencies at 183.310, 323.826, and 325.152GHz [9][10][11][12].In order to integrate more channels in a very compact space, the design of QON system often selects double-layer structure that involves the separation, guidance, and bunching of beams with the help of a serious of feed horns, different cascaded mirrors, and filters, which makes it possible that all channel signals with different frequencies and polarizations can mechanically share a single main scanning reflector to satisfy satellite payload constraints on cost [13].
Typically, the feed in QON systems often selects the circular corrugated horn.It can finish a transformation of circular waveguide mode from fundamental TE 11 mode to balanced hybrid HE 11 mode in the horn aperture, where the aperture field has the characteristics of very low crosspolarization level and circular symmetrical radiation pattern with Gaussian distribution characteristic [14][15][16][17].In the higher frequency band, limited by high machining accuracy in the grooves and junction, the machining quality of corrugated horn can easily affect its far-field pattern that leads to nonuniformity symmetrical Gaussian distribution due to the decreasing of main mode HE 11 .This problem often appears on the manufacturing during the process of spark erosion, classical milling, and electroforming [18].It can bring greater distortion and deformation of beam after passing cascade mirrors.
On the other hand, the FSS comprises periodic arrays of metallic patches (reflection) or apertures (transmission) that are fabricated in a conductor surface directly or printed on a dielectric-slab indirectly, where the different periodic structures contain center connected, loop, solid interior, or combination types.Although they can finish a frequency selection perfectly with low-pass or high-pass features, obtaining a narrow band in the MMW range and very rapid cut-off characteristic outside the working frequency band is difficult [19,20], so it is necessary to introduce the polarizer due to its frequency irrelevant attribute.
The polarizing grid is mainly used to separate double signals with horizontal and vertical polarizations [21].The most optimal structure is the freestanding grid with round wire for higher transmissivity, such as copper or tungsten wire.They can be wounded across a supporting metal frame and then glued in place just like being sandwiched between two frame halves [22], but it should accord with the requirement that ensures the accuracy in the mounting and positioning attachments.However, in the MMW and sub-MMW range, it is limited by the processing and manufacturing condition because the wire spacing g < λ/2 and round radius a < λ/10 for a good polarization discrimination, which leads to spending more time in the wire-drawing, stretched tight, and vulnerability when it has periodic array distribution in a compact space.Another selectable component consists of many parallel metallic strips and plays the same vital role in polarization processing.Those parallel metallic strips can be directly manufactured by chemistry corrosion technology or printed on a dielectric substrate following a photolithography process [23][24][25] and become easier to ensure the precision of surface.Although the introduction of dielectric substrate makes the transmission coefficient reduced and become related to frequency in a certain extent, the phenomenon of multiple reflections in the substrate can be solved analogously by an antireflection coating with smaller dielectric constant, where the thickness of dielectric substrate can choose the multiple of quarter of a wavelength.What is more, the cross isolation between both directions should be further discussed, which is related to the ratio of wire period and width.
The FSS and the polarizer should be tested in the QO test bench, where the transmitted and received feeds are simultaneously placed in the double focus of elliptical mirrors as a test system if the device under test is also put in the public focus place [26,27].During the FSS and grid measurement, the multipath stray signal and mutual coupling effect often occur from some reflections or transmission influencing the measurement results, which can be solved by a timedomain analysis technology [28].Meanwhile, the measurement of QON system is essential for assessing the overall system working performance.However, there is rarely information to describe such system test process both in the nearfield and far-field measurement, and the difference between the double methods should also be more discussed.Based on the advanced MMW and sub-MMW instruments in the BUPT, the far-field radiation patterns of the corrugated horns can be directly implemented in the tri-reflectors CATR [29][30][31].The near-field of dual-channel QON system at the output position has been obtained by a planar scanning in the two dimensions, where the sampling space can be chosen as dx = dy < λ/2 according to the Nyquist sampling theory; meanwhile, their far-field pattern of dual channels is also directly measured in the tri-CATR compared with nearto-far-field transformation by Fourier transform (FT) [32].Because inevitable stray signals existed in the surroundings, another method using aperture field integral with ideal plane wave is adopted to evaluate the far-field performance [33].
This paper mainly focuses on the design principle, measurement, and contrast method towards the dual-channel 3D QON system between FSS and polarizer.This dual-layer mirror structure can transmit 324 GHz channel on the top layer and reflect the 183 GHz to the bottom layer.The double channels are finally synthesized into a light path and then tested in the output position by near-and far-field.The mirror system based on the tracing of beam radius is firstly presented.The mode contents of feeds in the horn aperture are analyzed by the spherical wave expansion (SWE) method for verifying the asymmetrical pattern characteristic and energy loss during the near-field measurement.The principle of FSS and polarizer and their measurement results are also discussed in Section 2. The direct near-field planar measurement and three kinds of far-field results are compared in Sections 3 and 4, respectively.The corresponding experimental 2 International Journal of Antennas and Propagation results verify the WGP has better work performance in the good pattern symmetry, high cross-polarization isolation, low energy loss, and side-lobe level.At last, we share relevant design and testing experience about the QON system.

The Design Principle of QON System
The QON system mainly controls the change of beam radius and propagation direction.The conical corrugated horn feed can radiate the circular symmetrical GB.When the GB propagates forward from the waist ω 0 near the phase center of horn aperture, the beam radius ω changes with the distance of z and has a relationship of w = w 0 1 + λz/πw 2 0 2 For controlling the size of beam radius during the whole propagation path, the reflector mirror is mainly treated as a focusing unit to converge the beam radius thereby forming a new waist through it.When the lengths of ellipsoidal mirror from its double focal points to the center of surface are expressed as RI and RE, its equivalent focal length f can be expressed as 1/f = 1/RI + 1/RE.According to the paraxial approximation and the geometrical optic limit, ellipsoidal mirror can be treated as a thin lens, where the beam radius changes with input distance d in in the free space, followed by the lens with f , and finally forms a new waist in the output distance d out .After passing the lens through the distance d, the system transfer matrix using ABCD law can be expressed as [7,8] where 0 ≤ d ≤ d out .Due to the beam radiating from the waist position, when the complex beam parameter for the incident direction has q in = jz c = jπw 2 0 /λ, the parameter for output beam at any position d has Then the radius of output beam can be obtained by the expression of w out = sqrt λ/π/Im −1/q out .This ABCD matrix technique provides a good estimate for the characterization of beam size at different positions throughout the whole QON system.
2.1.The Dual-Channel 3D QON System.The dual-channel 3D QON system is modeled and simulated in the commercial GRASP software.It is shown in Figure 1, where the double corrugated horn feeds, three elliptical mirrors, and one filter using FSS or polarizer are installed.The double corrugated feed horns of T-183G and T-324G are similarly placed on the one focus point of elliptical mirrors M1 and M2, respectively, with distance D1 and D2.The lengths of long and short axis in mirrors have 218.6 mm and 180.9 mm, and the half rim of mirrors equals to Rim1 = Rim2 = Rim3 = 35 mm = 21.3λ183G = 37.8λ 324G .
When beams spread to the filter using FSS or polarizer with distance D3 = 172.1 mm and D4 = 162.86mm, the beam propagation direction of T-183G and T-324G is reflected and transmitted dividedly, where the reflection angle has θ = 30 °and then the double beams are synthesized into an emergent field by sharing the same reflection mirror M3 with a distance D5 = 92.12mm.
However, the double beam waists are not placed in one of the focus points of the next mirror M3, which brings a phenomenon of beam distortion due to existing higher-order Gauss-Hermite modes GB(0,1) and GB(0,3) [7].The distortion parameter U can be defined by U = w M tan θ i / 2 2f , and the reflected power of beam is given by K f = 1 − U 2 .A criterion of K f ≥ 0 99 requires that the diameter D of the mirrors should be four times the beam radius of ω M on the mirror surface, and it also needs to satisfy with a relationship of f /D ≥ 0 5 at the case of an incident angle θ i ≤ 30 °.In this dual-channel 3D QON system, the double beams passing the M3 have the distortion with U 183G = 6.69% and U 324G = 6.19%,where the corresponding reflection power is K f,183G = 99.5514% and K f,324G = 99.6162%,respectively.

2.2.
The Corrugated Feed Horn.The corrugated feed horn plays a very important role in the development of QON systems with high characteristics that its radiation pattern is Gaussian shaped and offers the potential with high gain, axial symmetry, low cross-polarization, and wide bandwidth.Generally, the corrugated horn is mainly divided into three sections that contain mode converter, corrugated profile, and phase stability section, respectively [14].The input waveguide only needs to be excited by quite pure TE 11 mode, which has dissymmetrical copolarization field in the E-plane and Hplane.The mode converter is required to finish the TE 11 -to-HE 11 mode conversion over a specified number of slots with a choice of variable-depth slot, ring-loaded slot, or variablepitch-to-width slot [15,34].The theoretical knowledge comes from the circular waveguide mode in the horn aperture that can be achieved into the balance condition with hybrid mode that has zero cross-polar field and is marked as HE 11 similar to 85% × TE 11 + 15% × TM 11 approximatively, when the deep of slots changes from λ/2 to λ/4 at the center frequency [35][36][37].The types of internal corrugated profiles can select linear, sinusoid, tangential, exponential, hyperbolic, or polynomial and so on.Nevertheless, the phase of both fields is not coincidence matching thus increases a phase stability period to achieve the field in phase at the aperture, where the phase stability section has a single-slot depth [14,16].

International Journal of Antennas and Propagation
In the 3D QON system, the double ultra-Gaussian corrugated feed horns with the cross of sine-squared and parallel section are designed and shown in Figures 3(a) and 3(b), respectively, where the profile section generates the HE 11 and HE 12 modes, and the parallel section brings them in phase [38].The design principle comes from the circular waveguide mode matching technique, where the waveguide step discontinuity is analyzed by the scattering matrix and optimized in the CHAMP software for low side-lobe and cross-polar levels.The both WR4 and WR3 feeds after machining are directly measured by the tri-reflectors CATR in the BUPT [30,39,40].The test results of co-and cxpolarization radiation field are presented in Figures 3(c) and 3(d) working at the 183 GHz and 324 GHz, where the side-lobe levels have −39 dB and −42.3 dB, and the crosspolarizations are −30 dB and −29 dB, respectively.Furthermore, by the spherical wave expansion (SWE) method, the mode components are obtained at the aperture of horns, where the spherical expansions expressed in the Q-modes with coefficients Q smn are written as [41,42] where the number of θ samples and ϕ samples over 360 °is n theta and n phi .The spherical modes have an azimuthal index m and a polar index n, the maximum polar index N is selected as 1 ≤ N max ≤ n theta /2 , and the maximum azimuthal index M is 0 ≤ M max ≤ min n phi − 1 /2, N max .The k 0 and η 0 are the wavenumber 2π/λ and free-space impedance.The F smn r, θ, ϕ represents the power-normalized spherical vector wave functions by the associated Legendre polynomials P m n cos θ and the spherical Hankel functions of the second kind h 2 n kr [41].When the Q-mode coefficients Q smn are solved, the tabulated SWE coefficients are fed as the tabulated feed.Then the aperture fields of double feed horns are obtained after passing own length of WR4 and WR3.The mode components of aperture field in the vector direction about WR4 at 183 GHz and WR3 at 324 GHz are calculated out and made up of HE 11 : HE 12 : HE 13 = 0.971 : 0.032 : 0.001 and HE 11 : HE 12 : HE 13 = 0.872 : 0.128 : 0.001.The planar aperture field is also shown in Figures 3(e) and 3(f), where the WR4 horn has axis-symmetrical pattern features both in far-field and aperture field.Due to higher difficulty in the processing with the rise of frequency, the aperture field of the WR3 feed has some noncircular symmetric Gaussian decay characteristics in the amplitude, which is caused by the decreasing of the level of main mode HE 11 and an increasing of another HE 12 mode.

The FSS.
The FSS in the dual-channel QON system selects a perforated metal plat plate with periodical circular hole array structure.This type of FSS has a bandpass filter characteristic, and the main design principle is the Floquet theory in the short circuit and the waveguide mode theory in the open circuit, where the minimum and the highest cut-off frequencies during the passband are determined by the cut-off frequency of waveguide unit and the higherorder Floquet modes, respectively [43].For the case of normal incidence, the reflection and the transmission coefficients of the circular openings arranged in equilateral triangular lattice can be reduced to following forms  where h is the thickness of metal plate, and the other parameters of A, B, and β can be found in [43].According to the transmission line theory, considering the circular hole as a circular waveguide, when the FSS needs to transmit f c = 324 GHz, its cut-off lambda λ c is linked to the hole radius a and satisfies the relationship with a = λ c /3.41 = 0.27 mm.From (4), if the spacing d between adjacent holes becomes smaller, the transmission bandwidth of passband is wider.The thickness h of FSS can influence the energy loss.What is more, it can further decrease the ripple and reflection coefficient in the passband, and the single thickness can select h = λ c /5 ∼ λ c /2.The structure of FSS having circular hole and equilateral triangular lattice is shown in Figure 4(a), where it has the radius a = 0.26 mm, the thickness h = 0.4 mm, and the size 4 cm × 4 cm and is adjacent to surrounding borders with 1 cm extension.Meanwhile, the incidence angle can influence its frequency response curve both in transmission and reflection but is not sensitive to incidence polarization wave, which are shown in Figures 4(b) and 4(c) at the angles of 0 °, 15 °, 30 °, and 45 °.The larger incidence angle is, the smaller transmission resonance frequency has.For the requirement of incidence angle 30 °, the FSS needs to play a good role in the beam separation at the 183/324 GHz.The measured transmission coefficient at 30 °has −30.21 dB at 183 GHz and −0.53 dB at 324 GHz with vertical polarization incidence, comparing the simulation values with −29.35 dB at 183 GHz and −0.16 dB at 324 GHz.However, due to the influence of coupling responses on two-feed radiation field and small structure size of holes, the reflection result is only given directly by simulation with −0.011 dB at 183 GHz and −19.74 dB at 324 GHz at 30 °.
2.4.The Polarizer.In this QON system, the broadband WGP printed on the dielectric substrate is designed.The polarizing grid can reflect perfectly the field that electric vector parallels to the direction of grid wires and transmit without loss orthogonally polarized radiation field.The analysis method of grids with dielectric substrate can adopt the transmission line equivalent circuit.The shunt reactance of an inductive grid sandwiched between air and substrate having index of refraction n 1 and n 2 , respectively, is [7] where it has w 0 ′ = w 0 2/ n 2 1 + n 2 2 , g and 2a are the period and width of the polarizing grid.
According to the equivalent transmission line model in Figure 5, the reflection and transmission coefficients can be expressed as R = Z − Z s / Z + Z s and T = 2Z/ Z + Z s , so the polarization isolation between both directions is related to the ratio of g/2a, and the final relation is shown in Figure 6 at different frequencies.Due to the machining becoming more difficult for small radius a, when the grid and the incident angle are fixed with a = 0.04 mm and θ i = 30 °, the polarization isolation decreases with the increase of frequency and has the variation with a period of g/2a = 1.Considering the actual processing precision in the printed circuit board (PCB), the pitch-to-width ratio of g/2a should be selected at the range of [3.2, 3.8], where the theoretical polarization isolation can be less than −42 dB.
Meanwhile, for increasing the transmission coefficient, the phase delay Φ of reflected field on the double interfaces of substrate should satisfy a relationship with Φ = 2βd and equal to odd times of π, so the thickness of dielectric substrate should be selected as a quarter-wave matching layer [44] where K 0 is an odd integral value, the θ i is the incident angle, and the propagation constant in the substrate has β = 2π/λ n 2 − sin 2 θ 0 5 .In the case of vertical incidence, it has min-

Near-Field Planar Testing Results of QON System
The near-field planar measurement of the dual-channel 3D QON system is implemented in the anechoic chamber.The test platform is shown in Figure 8 and is driven by the scanners with two-axis stepper motors.The magnitude and phase information of near-field are continuous sampling at regular intervals over a well-defined output plane with a size of 80 mm × 80 mm.Due to the fact that the maximal sampling interval of △x and △y is λ/2 from the Nyquist sampling rate, the planar testing intervals are selected as being less than λ/2 having dx = dy = 0.8 mm at 183 GHz and dx = dy = 0.46 mm in 324 GHz.The measurement distance of near-field has D6 = 160 mm in both 183 GHz and 324 GHz.Meanwhile, the both copolarization and cross-polarization fields are measured at the same time.When the filter selects the polarizer or FSS, the polarization direction of transmitted feed at 183 GHz channel is rotated as horizontal polarization for 6 International Journal of Antennas and Propagation reflection, and otherwise it is vertical polarization for transmission at 324 GHz.
Of course, the simulation model of QON system is also created in the commercial GRASP software.The single unit cell of FSS with equilateral triangular lattice is the plane with a central circular hole and limited by a tabulated rim that is given by means of a tabulation data.The table values define many points along the rim, in a plane perpendicular to the z-axis of the coordinate system of FSS, and are tabulated in rectangular XOY coordinates.According to the size of beam radius calculated in the filter position, the dimension of FSS structure is selected as the 24.89 mm × 24.89 mm and can satisfy the condition that the length is larger than the four-time beam radius both in double channels.There are a total of 529 units cell, and the single unit cell is defined by a rim in total of 2288 data points.Correspondingly, the polarizer is also established by the rectangular plate with electrical properties of strip grid and dielectric layer based on the simulation parameters.When the beam spreads to filter, the reflection fields of 183 GHz are obtained by the rigorous current analysis using the physical optic (PO) method.In order to compare the consistency among the different planar output fields, the scalar and vector conversion efficiency formulas η s and η v on behalf of beam energy distribution and vector direction are used.Especially, comparing with results of planar measurement f 1 ⋅ exp jφ 1 , 183 GHz channel selects the measured result of aluminum plate and the 324 GHz channel chooses the direct transmission through air as the reference field f 2 ⋅ exp jφ 2 .
The output planar near-field of 183 GHz channel is measured in the four cases using different filters, including FSS, polarizer, and direct reflectors by copper or aluminum plate.Thereinto, the double direct reflectors are adapted to compare the energy loss by the system itself.The testing planar grid has a total of 111 × 111 data number within the output range of 80 cm × 80 cm.The four testing results of two-dimension field are shown in Figure 9, where the amplitude and phase information of co-and cross-polarizer are given together.The more detailed parameters contain the cross-polarization level, the Gaussian beam radiuses, and the consistency of copolarization pattern between simulation and measurement and are summarized in Table 1.
In general, from the measurement results in Figures 9(a)-9(d), the four copolarization fields have circular symmetry Gaussian distribution characteristic in both the amplitude and phase, which is similar to the radiation field of WR4 corrugated horn feed in Figure 3(e).The beam radius can be well controlled after passing through the double reflector mirrors and one filter.Inversely, the amplitude of cross-polarization has double peak levels during the whole planar beam ranges especially at the case of using polarizer.Their phases have the same Gaussian attenuation characteristic with the spiral shape, where the phase converges uniformly to the center position of beam.Furthermore, those radiation patterns of amplitude and phase have the same discontinuity-like ripples in the bottom.This phenomenon can be interpreted by the system diffraction from the structure framework when the received horn is moved to the edge of system baseplate.
For the cross-polarization level, comparing the simulation and measurement results given in Table 1, the actual test levels are generally lower 3 dB∼4 dB than the simulation values.The test using FSS and polarizer can reach up to −18.33 dB and −18.82 dB, respectively, which are close to the direct reflection by plane mirror using copper and aluminum plate with −18.07 dB and −17.76 dB separately.
The accurate calculation of beam radius takes the average position values of both sides at the E-field and H-field.According to the simulation and measurement results, the four beam radiuses have a small difference relative to the theoretical radius value with ω out,183G = 11.07λ183G .This small change mainly comes from the uncertainty of GB waist position in the actual WR4 corrugated horn feed, the assembly location and angle errors of QON components, the shift of measuring planar position being far away from the center of M3 mirror, and the movement precision of scanning motors when it accurately reaches the test plane at the D6.Obviously, this difference of beam radius in D6 can reach the range from 2.89% to 5.03% relative to theoretical value seeing from all measurement at 183 GHz.For the most case of waist position, if the input distance D2 increases +1 mm and +2 mm, the theoretical output radius in the D6 will increase +0.4λ 183G and +0.8λ 183G , which are close to the simulation and measurement values.
In the beam consistency of copolarizer field, the scalar conversion efficiency has η s,polarizer = 98.7954% relative to the reference result of aluminum plate with a 100%.The FSS has also a similar result with η s,FSS = 98.9749%.Although the double energy loss has exceeded 1% with respect to the direct reflection case just only choosing a different filter, the energy loss of reflection with 0.9% should be a common phenomenon in the testing of QON system based on another direct energy calculation using copper plate with η s,copper = 99.0955%.
Meanwhile, in the vector conversion efficiency, the values become lower for obtaining bigger influence from the difference of phase.The three different filters have the highest ratio up to a total of η s,polarizer = 97.9952%,but it has a contrast that the double ideal reflection between copper and aluminum plate is only η s,copper = 97.5476%.It further demonstrates that the filters using FSS with perforated array and polarizer with strip grids can play a perfect reflection performance in the

8
International Journal of Antennas and Propagation 183 GHz channel.Of course, at the same test conditions for beam consistency distribution, the influence from the variation of phase is bigger than that from the change of amplitude when different filters are adopted.

The 324 GHz Channel Near-Field Planar Measurement.
The output planar near-field of 324 GHz channel on the front layer is similarly measured with the three different cases, including FSS, polarizer, and direct transmission through air as a reference.The planar grid has 201 × 201 data number in total within the output range of 80 cm × 80 cm, and those testing results are shown in Figure 10.The comparative values are also summarized in Table 2. Due to the fact that the transmissive effect traversing the filter cannot be simulated in the software, the double output fields in simulation are calculated by the same PO currents when the beam passes through the FSS or polarizer at 324 GHz.
From the planar measurement results shown in Figures 10(b) and 10(c), the double fields have a small noncircular symmetric characteristic at the copolarization amplitude when the level of beam pattern is far away from the center of −20 dB, which is caused by the same nonuniform    3(f).This asymmetry effect mainly comes from the corrugated horn because the consistency energy coefficient η s of polarizer is close to 99.96% relative to the direct transmission through air, which illustrates the correctness of transmission channel.Another interesting fact is that the diffraction from the mirrors' edge and system frames is less than that occurring at the 183 GHz channel obviously, so the amplitude and phase ripple have smaller ranges in the plane.
However, when the filter uses FSS exactly shown in Figure 10(a), the distorted beam appears more serious both in copolarization and cross-polarization with a beam-expanded peculiarity.Unfortunately, both amplitudes have a slower attenuation and larger beam shape when the amplitude levels are far away from the center of −10 dB in the horizontal direction, although both phases have small difference relative to the case of using polarizer and air obviously.This phenomenon results from a worse center flatness of FSS in the horizontal direction due to a small thickness with only 0.4 mm, and meanwhile this perforated array structure has many holes throughout the whole plate surface.On the other hand, with the increasing of insertion loss with −0.53 dB at the 324 GHz, the transmission result of FSS has the worst isolation in its cross-polarization field, which has larger warped profile range both in the amplitude and in the phase comparing the two others.
In more detail, comparing the testing case of direct air, the polarizer can play a good role in the beam transmission.Obviously, the FSS has caused some problems that are energy loss and directional variation both in the double amplitude, where the energy loss and vector change of copolarization amplitude have reached 3.4% and 6.3% relative to the polarizer seen from Table 2. What is more, the three beam radiuses by simulation and measurement have some difference likewise relative to the theoretical value with ω out,324G = 16.35λ324G .The measurement differences reach the range from 1.33% to 4.61% due to some test errors mentioned in the 183 GHz channel.What is more, when the beam passes through the filter, the theoretical calculation of beam radius is direct as air.The change of radius does not consider an influence from internal structure, which leads to bring a larger deviation than the direct reflection during the 183 GHz channel.Of course, for most cases of waist position, if the input distance D1 changes −1 mm and +1 mm, the theoretical output radius in the D6 will change −0.91λ 324G and +1.147λ 324G .Based on such facts, the good flatness and the low insertion loss of filter in the QON system should be considered more to reduce the beam distortion phenomenon especially in the high frequency.
Just like the cross-polarization level, −29 dB, of double transmitted WR3 feeds, the simulation polarization values of output field are also close to it.However, the measured values have a large difference, where the cross-polarization level using the FSS has only −13.45 dB but the polarizer can reach up to −19.93 dB that is better than the direct air with −18.44 dB.The planar cross-polarization level of polarizer also has faster attenuation characteristics than that in the FSS.This phenomenon happened in the case of transmission channel, and it can be understood adequately that the polarizer can isolate another polarization field in the horizontal direction when the polarization of input field has a vertical direction.Furthermore, the cross-polarization using polarizer has even more amplitude plane due to the improvement of polarization purity.

Far-Field Radiation Pattern Results of QON System
Traditionally, the method to indirectly obtain the 3D far-field radiation pattern is to make full use of the results of near-field planar measurement, where it finishes the near-to-far-field transformation from time domain to frequency spectrum domain by the Fourier transform (FT).In this method, the vector spherical wave spectrum function f t k x , k y is firstly obtained by this transformation based on double tangential E-field complex components E tx x, y, z = D6 and E ty x, y, z = D6 in the x and y direction, respectively, over the bounds of output plane S a , and then the far-field radiation pattern E r, θ, ϕ is further obtained in terms of spherical spectrum functions f x θ, ϕ and f y θ, ϕ shown as follows [32,37] Ê r, θ, ϕ = jk 0 exp −jk 0 r 2πr θ f x cos ϕ + f y sin ϕ where k 0 is the wavenumber and k x = k 0 sin θ cos ϕ and k y = k 0 sin θ sin ϕ are propagation constant.The direct far-field pattern measurement in the CATR using an approximate plane wave has become a trend in  International Journal of Antennas and Propagation evaluating the far-field radiation performance.The corresponding principle and test setup for the far-field CATR measurement can be found in [45,46].In this method, the antenna under test (AUT) is fixed on a rotary table in the QZ range.Then it rotates around the emergent QZ field in the different θ orientations, the radiation characteristic of   AUT can be obtained at any direction by recording the received fields in a very short period of time that is shown in Figure 11.The theoretical knowledge comes from the integral between the aperture field distribution A x a of AUT and the QZ field E qz that can be represented as [33] F θ = A x a E qz x, z dx a , 10 where each section position of θ direction needs to have an integral over the whole output aperture of AUT in Figure 11.The emergent direction of QZ field comes from negative z-axis, and the aperture of AUT is rotated in the XOZ plane.E qz x, z = E qz x a cos θ, −x a sin θ = a 0 x a cos θ, −x a sin θ exp jγ 0 exp −jk 0 z is the incident QZ field on the AUT with a θ orientation.a 0 x, z = 1 and γ 0 = 0 are the amplitude and phase of the ideal planar wavefront.
In the actual measurement environment, there are still some stray signals from the environments, such as the diffraction phenomenon from the testing system structures, the influence of mirror surface error from larger reflectors in creating the QZ field, and the random noise in measurement equipment, which can bring the ripple of amplitude and phase in the QZ field, where this typically predicted amplitude ripple is within ±0.5 dB, and phase ripple is within ±5 ° [30,39].Those influences from the quality of   Due to the fact that far-field pattern testing in the CATR needs double same corrugated feed horns in the transmitter and receiver each other, the far-field of 183 GHz channel is merely indirectly obtained based on the near-field measurement and the near-to-far-field transformation for only using one WR4 corrugated horn, however, the 324 GHz channel is directly measured in the tri-CATR by double WR3 corrugated horns, which is implemented from the range −90 deg to 90 deg within a sampling interval of 0.1 deg, and then those results are compared by the same aperture integral method using ideal plane wave.The simulation results of far-field are also similarly obtained by the PO current method in GRASP.The different far-field results in such double channels are also discussed and compared below.
4.1.The 183 GHz Far-Field Results.The four planar far-field amplitudes by near-to-far transformation are shown in Figure 12(a) based on each near-field test results.The cut section results of the FSS and the polarizer are also compared with amplitude and phase in Figures 12(b) and 12(c), where the simulation, near-to-far-field transformation based on the near-flied measurement, and the aperture integral using ideal QZ field are compared and summed up in Table 3.
Benefiting from the near-field planar test with a feature of circular symmetry Gaussian distribution, the far-field copolarization planar amplitude by near-to-far-field transformation is also consistent both in the θ and in the ϕ direction, so the four radiation patterns have small difference in the copolarization amplitude shown in Figure 12(a), where the three scalar and vector coefficients of planar far-field are also calculated relative to the results of aluminum plate that have η s,FSS = 99.5457% and η v,FSS = 98.7896% for FSS, η s,polari- zer = 99.9157% and η v,polarizer = 99.2113% for polarizer, and η s,copper = 99.5442% and η v,copper = 98.8688% for copper plate.Those transformation efficiencies are obviously superior to that in the near-field measurement shown in Table 1, which mainly benefits from the rapid attenuation in far-field.
Comparing with the first side-lobe level of copolarization amplitude both in phi = 0 deg and in phi = 90 deg from Figures 12(b) and 12(c) and summarized in Table 3, the calculation values by the near-to-far-field transformation have −40 dB to −32 dB, which is generally higher 10 dB∼20 dB than the simulation values but close to the aperture integrating values.Both differences with respect to the simulation come from the stray signals in testing environment and are limited by the side-lobe level of transmitted WR4 horn with only −39 dB.It is also obvious that the first side-lobe level cannot be significantly decreased by the ideal QZ field integral for relying on the test condition of near-field.The cut section of copolarization phase has the same shape among the three far-field results, but the phase of simulation has a difference with 2π comparing with the two others.
The cross-polarization level in the simulation is close to the case of near-to-far-field transformation.The isolation effect of the aluminum plate is slightly superior to the copper plate, but the polarizer has the best isolation effect than three other results both in phi = 0 deg and in phi = 90 deg, where the cross-polarization level reaches −19.74 dB and −30.03 dB, respectively, by the near-to-far-field transformation.When calculated by the aperture integrating, the far-field in the double directions is achieved using the same cut section field based on near-field measurement.Due to the QZ field being an ideal plane wave in the aperture integrating, the cross-polarization level can be changed both in the phi = 0deg and in phi = 90 deg.The biggest change happens when the filter uses polarizer, so the quality of QZ field can influence the cross-polarization level of radiation far-field; however, it is difficult to forecast whether it is good or not because there are some randomness amplitude and phase ripples in the actual QZ field when the AUT is measured at the different positions.

4.2.
The 324 GHz Far-Field Results.The three far-field planar amplitudes of copolarization are similarly obtained by a near-to-far-field transformation and shown in Figure 13(a).Just like an influence in the near-field test, the far-field pattern using the FSS has obvious pattern distortion and asymmetry especially in the section direction near phi = 0 °, where the coefficients of far-field planar pattern have η s,FSS = 98.7915% and η v,FSS = 96.9797%relative to air.However, the far-field planar pattern of polarizer has a very good consistent circular symmetry feature both in the θ and in the ϕ direction relative to the air, where it has the coefficients of η s,polarizer = 99.9887% and η v,polarizer = 99.5524% that can verify the polarizer has a better transmission performance at the 324 GHz.Of course, the planar far-field has bigger noncircular symmetry Gaussian distribution in the phi = 90 °due to the asymmetric aperture field of the WR3 horn itself.The cut section results of far-field patterns using the FSS and the polarizer are also presented in Figures 13(b) and 13(c), where the simulation, near-to-far-field transformation, and the aperture integral values using ideal QZ are given.What is more, the horizontal cut section of far-field pattern has been obtained at 324 GHz by the direct far-field measurement by the tri-reflector CATR system in the BUPT [29].Their far-field results of amplitude and phase have been shown in Figure 13(d).Those relevant results are also compared and summed up in Table 4.
Comparing with the side-lobe level of copolarization amplitude, the near-to-far-field transformation can reach −45 dB∼−39 dB, which are general higher 10 dB∼20 dB than     When the 324 GHz channel is directly measured in the tri-CATR, their side-lobe levels have an obvious difference.This variation has been fully shown in Figure 13(d).It is the case of without any filter that has the best side-lobe level of −46.57dB, which has thoroughly closed to the ideal simulation situation.When the filter adopts the polarizer, the effect in the side-lobe becomes excellent as low as up to −45.22 dB, so the polarizer has the same transformation characteristics just like through the air.Of course, the FSS has also the worst side-lobe level for only −22.67 dB.Meanwhile, the tested copolarization phase has the same shape with periodic-serrated oscillation function, where it has a periodic 2π change with the angle θ that are different from the simulation with a Gaussian distribution characteristic in the main-lobe range.
Although the far-field calculation by near-to-far and aperture integral are similar both in cross-polarization levels, their test in the CATR still has a large variation.From the cut section results, the terrible cross-polarization level still happens in the FSS for only −11.64 dB that is lower near −19 dB relative to the case of direct air.The good phenomenon similarly appears on the polarizer, where the crosspolarization level can reach up to −34.93 dB that has a similar and better isolation than that through the air with −30.3 dB.Based on those facts, it can be derived that the direct test of far-field radiation pattern in the CATR still has an important and indispensable role in evaluating the QON system performance especially in the cross-polarization and side-lobe levels just like the simulation.

Discussion and Conclusion
A dual-channel 3D QON system has been analyzed and compared using polarizer and FSS that transmits the 324 GHz signal on the top layer while diverting the 183 GHz signal to the bottom layer.The changes of beam radiuses are firstly traced in such system by the complex beam parameter and system transfer matrix for deciding the size of mirrors and calculating the beam distortion.Of course, the quality of corrugated feed horn can influence the radiation feature of output field, so the mode components of double corrugated feed horns are calculated over the horn aperture by the SWE method, which is used to explain the reason of measurement with asymmetry pattern in the near-field and far-field results in 324 GHz channel.What is more, the design principle of FSS with perforated hexagonal array and WGP printed on the dielectric substrate is also discussed and verified by testing transmission coefficients.
In the planar near-field measurement, the 183 GHz channel is reflected at four different cases with FSS, polarizer, copper, or aluminum plate.Meanwhile, the 324 GHz channel is transmitted through FSS, polarizer, or direct air.Generally, the beam radiuses in both channels can agree very well in the all measurement results, but the difference with the theoretical values is small and depends on the actual measurement environment conditions.In the 183 GHz channel, comparing the planar copolarization and cross-polarization results, the four different reflectors can achieve very good reflection results both in the amplitude and phase.Their symmetry of planar patterns, cross-polarization level, and side-lobe levels is very similar.However, in the 324 GHz channel, the cross-polarization level and side-lobe level using the polarizer have the best performance.What is more, the symmetry of planar near-field using FSS further decreases when the beam is transmitted.
The far-field patterns of 183 GHz are indirectly obtained based on different methods.The near-to-far-field transformation and aperture integral using ideal plane wave depend on the measurement condition of near-field.Their planar and cut section of far-fields are solved out and discussed both in the phi = 0 deg and in the phi = 90 deg, which have similar results just like the near-field.Meanwhile, the direct far-field measurement of 324 GHz is also implemented by the trireflectors CATR with an anechoic chamber in the BUPT.In this case, the cross-polarization level and the side-lobe level are more outstanding.Although the near-field measurement in the QON system is time-consuming for satisfying the Nyquist sampling ratio, it can foresee the working performance within whole output plane especially on the symmetry of planar pattern.However, it is also often limited by many test conditions, such as the stray signals, the position drift from scanning frame, and the diffraction phenomena.In the far-field radiation pattern, the aperture integral using ideal QZ field can obtain equivalent cross-polarization level relative to the near-to-far-field transformation, but it has no very good improvement in the side-lobe level for relying on the near-field results.Comparing with such double far-field

6 Figure 2 :
Figure 2: The changes of beam radius in the dual-channel 3D QON system.
Aperture field of WR4 feed at 183 GHz.From left to right are aperture field, HE 11 , HE 12 , and HE 13 Aperture field of WR3 feed at 324 GHz.From left to right are aperture field, HE 11 , HE 12 , and HE 13

Figure 3 :
Figure 3: The radiation far-field and aperture patterns of corrugated feed horns.
substrate thickness with d ε,min = λ/ 4 ε .The dielectric material selects PTFE with low dielectric constant for good transmission.The structure model of inductive strip grid is shown in Figure7(a), where the size has 6 cm × 6 cm.The simulation of TM wave reflection and TE wave transmission and their polarization isolation at the different incidence angles are also shown in Figures7(b) and 7(c), where the incidence angle has larger influence in the transmission coefficients than that in the reflection.The simulation of TM wave reflection coefficients is directly given that there are 0.022 dB at 183 GHz and −0.076 dB at 324 GHz, and the TE wave transmission coefficients have −0.008dB at 183 GHz and −0.36 dB at 324 GHz for the incidence of 30 °.The polarization isolation of TE wave transmission has been measured from 110 GHz to 325 GHz at the vertical incidence, where its isolation is processed by the time domain analysis technology for reducing the influence of stray wave signals and the multipath coupling effect, which are plotted in right axis of the Figure7(c) and have −25.33 dB at 183 GHz (simulation is −28.41 dB) and −19 dB at 324 GHz (simulation is −18.87 dB) at 30 °.
Single-layer FSS has circular hole with equilateral triangular lattice 150 deg (b) The simulation and measurement results of transmission coefficients Simulation, 45 deg (c) The simulation results of reflection coefficients

Figure 4 :
Figure 4: The structure of FSS and its frequency response curve at the different incidence angles of 0 °, 15 °, 30 °, and 45 °.

Figure 5 :
Figure 5: The structure of grid printed on the dielectric substrate and its equivalent transmission line model.
Sim, 15 deg Sim, 30 deg Sim, 45 deg (b) The results of TM wave reflection and its polarization isolation TE wave transmission coefficient (dB) The results of TE wave transmission and its polarization isolation

Figure 7 :
Figure 7: The structure model of WGP and its frequency response curve with different incidence angles of 0 °, 15 °, 30 °, and 45 °at TE and TM wave.

9
International Journal of Antennas and Propagation symmetrical aperture field distribution of the WR3 corrugated feed horn at 324 GHz shown in Figure

Figure 8 :
Figure 8: The planar near-field measurement of dual-channel 3D QON system.

Figure 9 :
Figure 9: The near-field measurement results of 183 GHz channel in the dual-channel 3D QON system.From left to right are the measured plane of copolarization amplitude, the copolarization phase, the cross-polarization amplitude, and the cross-polarization phase.

Figure 10 :
Figure 10: The near-field measurement results of 324 GHz channel in the dual-channel 3D QON system.From left to right are the measured plane of copolarization amplitude, the copolarization phase, the cross-polarization amplitude, and the cross-polarization phase.

Figure 11 :
Figure 11: The far-field pattern measurement in the CATR.

Figure 12 :
Figure 12: The far-field radiation pattern results of 183 GHz channel.

Figure 13 :
Figure 13: The far-field radiation pattern results of 324 GHz channel.

Table 1 :
The near-field measurement data of 183 GHz channel.

Table 2 :
The near-field measurement data of 324 GHz channel.

Table 3 :
The far-field results of 183 GHz channel.