Free-Space Materials Characterization by Reflection and Transmission Measurements using Frequency- by-Frequency and Multi-Frequency Algorithms

The knowledge of the electromagnetic constitutive properties of materials is crucial in many applications. Free-space methods are widely used for this purpose, despite their inherent practical difficulties. This paper describes an affordable free-space experimental setup for the characterization of flat samples in 1–6 GHz in a non-anechoic environment. The extracted properties are obtained from the calibrated Scattering Parameters, using a frequency-by-frequency solution or a multi-frequency reconstruction. For the first, we describe how the Time-Domain Gating can be implemented and used for filtering the signals. For the latter, a weighting factor is introduced to balance the reflection and transmission data, allowing one to have a more favorable configuration. The different role of transmission and reflection measurements on the achievable results is analyzed with regard to experimental uncertainties and different noise scenarios. Results from the two strategies are analyzed and compared. Good agreement between simulation, measurement and literature is obtained. According to the reported results for dielectric materials, there is no need of filtering the data by a Time-Domain Gating in case of the multi-frequency approach. Experimental results for Polymethylmethacrylate (PMMA) and Polytetrafluorethylene (PTFE) samples validate both the setup and the processing.


Introduction
Characterization of the electric permittivity ε, the magnetic permeability µ and the electric conductivity σ is of fundamental importance in physics and engineering, since these constitutive properties determine the response of the Material Under Test (MUT) to the electromagnetic fields [1,2].Usually, such properties are extracted from the Scattering (S-) Parameters [1][2][3][4], which can be measured by a Vector Network Analyzer (VNA) [3].Several measurement setups can be conceived for this purpose, including the Free-Space Method (FSM) [1].
The FSM consists in connecting antennas to the ports of a VNA and using them to sense the propagation path of interest; in this way, different materials can be characterized by means of proper reconstruction algorithms [1,[3][4][5][6].The FSM is contactless, requires little sample preparation and is suitable for heterogeneous samples.It can be used for measuring large flat solids, materials under high temperature [7], liquids [8], gases [1,6], radomes in aerospace [9], multi-layered dielectrics [10], composite materials and metamaterials [11][12][13][14].One should note that the effective response, the non-invasive nature and the wideband operation are essential in many actual applications.
The FSM is commonly used to characterize building materials in order to estimate their influence on telecommunications [15][16][17] or in through-wall imaging [18].Moreover, the knowledge of the electromagnetic properties can be used to evaluate the physical conditions of materials in different areas, such as masonry structures in civil engineering [17,19,20], the characteristics of snow [21] or even agricultural materials [22].There are also several works in millimeter wavelengths, for instance, for characterizing dielectrics [23], electromagnetic absorbers [24] and nanocomposites [25].
The concept of transmission and reflection FSM is shown in Figure 1.Antennas and the space therein can be considered as a two port network whose behavior is characterized by its four S-Parameters: S 11 , S 12 , S 21 and S 22 [3,26].After the measurement, the S-Parameters are used for extracting the electromagnetic properties of the sample.However, the measured S-Parameters are influenced by multipath propagation, unwanted reflections and other interferences, especially in a non-anechoic environment.As a consequence, additional signal processing is often required.Two techniques are important to deal with the above mentioned free-space configuration: (a) the free-space calibration [1,14,27], intended to normalize the measurements and to dislocate the calibration planes to the surfaces of the sample; and (b) the Time-Domain Gating (TDG) [1,3,[26][27][28][29][30][31], intended to filter the measured signals.After such preprocessing, the resulting signal can be used to recover the material properties by different extraction techniques.
In this paper two recovery approaches are presented and analyzed, using simulated and experimental data.One approach assumes to deal with non dispersive materials and its multifrequency solution is based on the minimization of a properly defined functional.The other is a frequency-by-frequency approach, in which a possible dispersion in the measurement range is assumed, without a previous knowledge of the dispersion model of the material.Both the approaches rely on the same one dimensional model of scattering, where the measured S-Parameters are interpreted as the reflection and transmission coefficient of an equivalent transmission line.
In particular, some contributions of this work are: (a) a weighted and normalized cost function is introduced for the multi-frequency extraction from S 11 and S 21 ; it is studied in terms of signal to noise levels, different weights for transmission/reflection and it is applied to free-space measurements; (b) an affordable experimental setup is presented and its calibration is discussed; (c) an analysis of uncertainty is presented; (d) the implemented TDG is explained and its applicability to the extraction methods is discussed.
The paper is organized as follows.In Section 2 we present the mathematical model, the extraction algorithms, the experimental setup and its calibration and, at last, the signal processing.The effect of the measurement uncertainties is analyzed in Section 3. Afterwards, numerical and experimental results are presented and discussed in Section 4. The conclusions are in Section 5.

Slab Scattering Model
The relationship between the scattered field and the electromagnetic properties of the material under test is provided by assuming an one dimensional geometry, where an incident plane wave is reflected (by) and transmitted (through) a planar homogeneous slab of relative electric permittivity ε r , relative magnetic permeability µ r and electric conductivity σ.This allows to apply the equivalent transmission line equations for time harmonic fields at the frequency f [2].
In this way, the reflection coefficient, defined as the ratio between the reflected and the incident field at the first slab interface, can be expressed as the following function of the slab properties: where the auxiliary variables y = µ r /(ε r − jσ/ωε 0 ), τ = tan(ωxd/c) and x = µ r (ε r − jσ/ωε 0 ) are introduced; d is the thickness of the sample, j is the imaginary unit, ω = 2π f is the angular frequency; c and ε 0 are the free space electromagnetic wave speed and dielectric permittivity, respectively.Note that y and 1/x are the impedance and wave speed of the material, normalized to free space, respectively.In its turn, the transmission coefficient, defined as the ratio between the transmitted field at the second slab interface and the incident field at the first slab interface, can be expressed as:

Multi-Frequency Extraction
When the MUT is non dispersive in the investigated frequency band, a possible strategy in order to estimate the unknown material properties from the reflection and transmission coefficients is to minimize the "distance" between theoretical and measured quantities.Optimization-based reduction is interesting because correlations between frequency measurements are allowed [32].It can be done by searching for the global minimum of a normalized cost function, here defined as: where || • || 2 is the quadratic norm.The scattering parameters S 11 and S 21 play the role of data of the problem, that can be either measured by the VNA or numerically simulated.The value of the weighting coefficient w can be fixed in the interval 0 ≤ w ≤ 1, in order to enhance the role of either the reflection or the transmission.
The non linear nature of the relationships (1) and ( 2) makes the functional (3) non-quadratic, so that local minima can occur.Moreover, ill conditioning may affect the minimization procedure when reflection and transmission weakly depend on the material properties.In this case, a proper tuning of w can help in improving the functional behavior.
If the MUT is dispersive but its properties obey a dispersion model, the Debye or Lorentz for example, it can be included in the procedure and its characteristic parameters must be calculated [1,4].Otherwise, a frequency-by-frequency approach can be used, but this one does not exploit the multi-frequency nature of the measurement, as explained next.

Frequency-by-Frequency Extraction
Here, the complex scattering parameters S 11 and S 21 are input data to three extraction algorithms, as illustrated in Figure 2 and explained afterwards: (a) the Nicolson-Ross-Weir (NRW) Algorithm, (b) The Reflection-Only (RO) Algorithm, and (c) The Transmission-Only (TO) Algorithm.These are frequency-by-frequency methods, since the properties are extracted for each particular frequency [1], in a point-to-point sweep.In the following the relative permittivity and permeability are assumed to be complex, that is, ε r = ε r − jε r and µ r = µ r − jµ r , respectively.The NRW [33,34] is the most traditional algorithm for retrieving ε r and µ r from S 11 and S 21 .In a transmission line, these parameters are defined as in ( 4) and (5), where Γ and T are the interfacial reflection and the propagation coefficients, respectively.Because it is a well-established algorithm, its development is left to Appendix A.
However, when using NRW the results diverge for low loss materials in multiple integers of half wavelength in the sample.At these frequencies |S 11 | is too small and the equations become algebraic unstable as |S 11 | → 0.Moreover, the uncertainty in the measured phase is high for small |S 11 | [1,6,35,36].In contrast, when the sample is thick and lossy, S 21 is small and its measurement can present a low Signal-to-Noise Ratio (SNR).
Alternative algorithms were also developed, exploiting only reflection or only transmission measurements [1,6,35,36].However, in this case only one between ε r and µ r can be determined, due to the reduced available information.Because ordinary materials are usually non-magnetic in microwave frequencies, µ r = 1 is often assumed and ε r is calculated by finding the root of a function F(ε r ), as described next.
The Reflection-Only (RO) method is particularly interesting when S 21 is unavailable or unreliable.In this case only S 11 is used.In order to have the function F RO in (6), one should consider µ r = 1 in (1) and calculate the difference between the datum S 11 and the function F R .The solution is found by forcing it to zero: where ν = d(ω/c) √ ε r .The Transmission-Only (TO) method uses only S 21 and it is more interesting when S 11 is more susceptible to errors.The function F TO in (7) is obtained by considering µ r = 1 in (2) and comparing the datum S 21 and the function F T .Again, after some rearrangements the solution of the problem is the value of ε r that satisfies the equation: One could use also the noniterative method in [36] to overcome the instabilities of NRW; it is restricted, however, for low-loss non-magnetic materials.

Experimental Setup
The experimental setup is at the Laboratory of Applied Electromagnetics of the Federal University of Minas Gerais.In order to reduce the reflections on the structure, metallic materials are avoided.So, as can be seen in Figure 3, the structure is made with wood and plywood.Each antenna is attached to a custom-made plastic telescoping stand.Aiming simplicity and low cost, the setup does not include electromagnetic absorbers.The measurement setup has 2.1 × 0.85 × 0.85 m (height, width, depth).The sample holder is a squared frame with 0.5 m edges and the sample is supported by a 0.5 cm border around the squared hole in the middle of the setup.Because of that, the sample should be quite rigid.Alternatively, a rigid block of styrofoam can be used to support smaller samples or samples that bend.
We use two broadband horns with dielectric lenses (model Q-Par WBH1-18S14L).The lenses help to concentrate more radiation on the sample.The test range is 1-6 GHz, limited by the antennas (1-18 GHz) and by the VNA (9 kHz-6 GHz).
It is crucial that most of the electromagnetic energy passes through the sample.Thus, the dimensions of the setup are defined in a way that the half-power beamwidth (HPBW) [37] is entirely on the sample for the frequencies in our range.
On the actual configuration, the distance r between the antennas and the sample can vary between 0.3 m and 0.48 m.Furthermore, the design also permits that the main part of the radiation reaches the sample with similar phases, that is, the far-field condition [37].Thus, r satisfies ( 8), where D is the largest aperture of the antenna and λ is the wavelength [1,37].Satisfying (8), the free-space can be approximated as an uniform transmission line [1].In our case, the last condition in ( 8) is barely satisfied with r = 0.3 m when D is large compared to λ (above 2 GHz, approximately).

Calibration
Two physical calibrations are performed in order to obtain the scattering coefficients of the MUT: (a) the VNA calibration, and (b) the free-space calibration.The first one is the Full 2-Port Calibration, in which the standards Through, Open, Short and Match (TOSM) of the Calibration Kit are used to compensate the effects of the cables [3,26].Thus, the reference planes are shifted to the antennas connectors.The second one is the free-space calibration, where two reference measurements are performed, one in absence of the sample and the other in presence of a metallic plate.The aim is to remove undesired contributions and to move the calibration plane to the faces of the sample [27,38].
Each calibrated scattering parameter is obtained by three measurements: (a) measuring a metal plate with the same cross section of the sample; (b) leaving the setup empty; and, finally, (c) measuring the material under test (sample) [27,38].Figure 4 illustrates possible contributions to the measured signals: reflections at the antennas due to mismatches (arrow 1); diffracted signal (arrow 2); signal transmitted through the sample and its secondary components (arrow 3); multiple reflections on the sample (arrow 4); signal reflected somewhere and passed through the sample area (arrow 5); signal reflected somewhere but did not cross the sample area (arrow 6); signal reflected somewhere and returned to antenna 1 (arrow 7); external influences that are not systematic and must be treated as noise (arrow 8); signal reflected on the metallic plane (arrow 9); and the transmitted signal at the empty setup (arrow 10).Observe that the arrows in Figure 4 can be interpreted as the error terms of the signal flow graphs (errors-adapters) in microwave networks [3,17,26].
It has to be stressed that the reversed path, with S 12 and S 22 , is not used in this work in order to reduce the measurement complexity.
In particular, measurement with empty setup (air) S (air) 11 allows to evaluate the contributions indicated by the arrows 1 and 7.In this way, S (air) 11 should be subtracted to S 11 for calibration, as in (9).Then, MUT measurement is normalized with that in the presence of the metal plate (metal), whose ideal reflection coefficient is −1.In a dual way, the contributions indicated by the arrows 2 and 6 are evaluated by measurement with metal S (metal) 21 .As in (10), it is subtracted for calibrating S 21 and we normalize it with the measurement in air, whose ideal transmission coefficient is 1 [17,38].
The exponential terms in ( 9) and ( 10) adjust the phases, since we have to consider that the sample and the metal plate actually substitute a slab of air [15].As illustrated in Figure 5, the correction is made using the thicknesses L 1 of the metal plate and d of the sample.In our experimental configuration, L 1 has a minor effect, since we use a very thin aluminum foil (0.01 mm).However, the effect could be not negligible in case one use a electrically thicker metal plate.In other words, it is important that the calibration planes have a considerable coincidence with the faces of the sample.This is challenging, since a different bending or deformity on the metallic plate with respect to the sample (or vice-versa), affect the phase, especially for higher frequencies.Although there are different calibrations applicable to FSM, such as the Thru-Reflect-Line (TRL), Thru-Reflect-Match (TRM) [1,3] and the Gated-Reflect-Line (GRL) [39,40], the described free-space calibration is a post-processing procedure that does not involve movements of the antennas.Again, it is important to prevent movements on the setup, since small mispositionings can generate great errors when using free-space methods [1,27].The classical TRL and TRM calibrations demand three reference standards.In particular, TRL demands an accurate positioning system to perform the Line standard.By its turn, TRM demands a well-matched absorber to perform the Match standard, which is a challenge for broadband measurements in free-space [40].
Similarly to the GRL, our calibration demands only the Thru (air) and Reflect (metal) standards.It is a two-tier procedure, since we firstly perform the TOSM to remove the effect of the cables.Note that the new GRL, that is, the GRL-1 (one-tier) described in [40] does not need the TOSM at the end of the cables.One should remember, however, that the measurement bandwidth and the number of points are limited, so a long network could result in a low Time-Domain (TD) resolution.
Because of the practical difficulties, thickness and position invariant calibrations has been investigated both in FD [41][42][43][44] and in TD [45].Observe that the TD analysis has been used in the GRL calibration [39,40], for extracting a desired part of a network [46] and also in TD extraction from measured FD data [45].In this way, it is interesting to have a flexible procedure for the signal processing.
Note that our extraction described in Section 2.2 differs from those introduced elsewhere [32,[41][42][43][44] in a way we do not use the reverse path (S 12 and S 22 ), since it would demand the metal plate to be positioned also on the bottom face of the sample in order to perform the calibration.It would require an additional repositioning of the sample, which is, as explained, a source of errors.
We can assume that the systematic errors are removed by TOSM and free-space calibrations.The random errors, in their turn, are usually unpredictable and cannot be removed by the calibration [1].However, the filtering property of the TDG can minimize part of them.Here we do not consider drift errors, since we assume that the working conditions do not change substantially from the calibration ones.In other words, measurements are performed no more than few minutes after the calibration and after the warming-up time of the VNA.

Filtering the S-Parameters
After the calibration, the S-Parameters mainly represent the sample itself.Despite of that, they are quite noisy, so that a filtering procedure is introduced here.The TDG is applicable for different transmission line setups.It is, however, an optional tool with a specific license in a typical VNA.Other authors have implemented gating codes, for instance [27,30], but yet some details are not much discussed in literature.Our gating algorithm is described in Appendix B, in which a flow chart is used to illustrate the overall process.
In Figures 6 and 7 both the noisy and the filtered measurements referring to S 11 and S 21 of a Polymethylmethacrylate (PMMA) sample with d = 10.2 mm are shown.For comparison, reference signals simulated by ( 4) and ( 5) with ε r = 2.61 and µ r = 1 are also reported.The VNA provides the S-Parameters in frequency-domain (FD), but in many cases it is worthy to analyze them also in time-domain (TD).In this way, one can remove unwanted contributions, criteriously.Due to the multiple reflections and transmissions, we have multiple peaks in TD, with different magnitudes and delays.This separation in TD is useful because it can facilitate the identification of the unwanted contributions.The spreaded noise in Figures 6 and 7, for instance, corresponds to side lobes in TD (see Figures 8 and 9), which can be reduced by a filtering process (gating).Note that in practice this identification is impaired by frequency bandwidth (BW) limitations, since the TD resolution is inversely proportional to BW [1,3].
One should note that it is not necessary to calculate the transform in the whole default range, since we are only interested in the details in the vicinity of the sample.The inverse Chirp Z-Transform (CZT) is then used [47].It is more general and flexible than the Discrete Fourier Transform (DFT) because it controls the mapping between the domains by an initial point t 1 and an increment δt [48,49].Note by Figures 6 and 7 that the filtered data match very well the frequency behavior of the simulated scattering parameters.However, a spurious low-level oscillation is still present due to: (a) the inherent difficulty in sharply separating the noise from the genuine signal when they are too close in TD.In particular, this gets worse when the TD resolution is not so high to avoid overlapping, that is, unresolved peaks; and (b) the gated signal presents errors concentrated on the frequency band edges and also a small oscillation over the whole band [50].Therefore, a residual distortion is a typical characteristic of gated signals and, in our context, it leads to some distortion on the extracted properties.

Analysis of the Uncertainties in the Model
For the estimation of the uncertainty in the extracted ε r due to the uncertainties in the measured S-Parameters and in the thickness d, we assume µ r = 1 and use a differential uncertainty analysis [41].Then, the total uncertainty in ε r is where ∆ means the uncertainty in the respective quantity, θ m is the phase of S m , m = 11 or 21 and p stands for the real or imaginary part of ε r .Because we do not have an explicit function for ε r , the partial derivatives in (11) are calculated using the implicit function theorem [51], after calculating ∂F RO /∂|S 11 |, ∂F RO /∂θ 11 , ∂F RO /∂ε r , ∂F RO /∂d, ∂F TO /∂|S 21 |, ∂F TO /∂θ 21 , ∂F TO /∂ε r and ∂F TO /∂d.Due to the non linear nature of the involved functions, the total uncertainty depends on the working point.Therefore, in order to quantify (11), we use the VNA specifications [52] and assume that ∆d is the standard deviation of d, which we measured in different points of the sample.

Individual Contributions on the Total Uncertainty
For the analysis of the behavior of ∆ε r as a function of frequency, we fix the evaluation point at ε r In Figures 10 and 11 the individual contribution of ∆|S m |, ∆θ m and ∆d to the uncertainty on the real and imaginary part of ε r are shown (lines with markers), together with their total effect (dashed line).In general, we observe that: (a) the major contribution to ∆ε r is due to |S 11 | when using reflection and is due to θ 21 when using transmission (see Figure 10); thus, S 11 data is very important when θ 21 is not reliable; (b) the major contribution to ∆ε r is due to uncertainties in the magnitude of both the S-Parameters, but their effect as a function of frequency is much different in reflection and transmission configurations, i.e. one increases with f , the other decreases (see Figure 11).

Uncertainty for Different Materials and Thicknesses
Here we study the uncertainty when the evaluation point in (11) changes.Specifically, we keep fixed ε r = 0.005 and make ε r and thickness d vary within a predefined range.In order to simply quantify the uncertainty taking into account the whole frequency range, we defined it as

Results and Discussion
In this section two materials are used for the measurements: (a) Polytetrafluorethylene (PTFE), and (b) Polymethylmethacrylate (PMMA).The samples are 50 × 50 cm with d = 10.5 mm (PTFE) and d = 10.2 mm (PMMA).The measured standard deviations of d are 1.59% and 0.97%, respectively.Two distances r between the antennas and the sample are considered: 30 cm and 48 cm.
When TDG is used, it is assumed the suitable set [31]: 5 GHz, g t 1 = −0.5 ns, g t 2 = +0.5 ns, 68 dB of sidelobe attenuation and a 0.5 ns cut-off time [28].Details about these parameters are in Appendix B.

Frequency-by-Frequency with Simulated Data
Here we use simulated data to validate the extraction with TDG in different SNR scenarios.Arbitrary white gaussian noise has been added to the theoretical S-Parameters with the following SNR: 5, 10 and 20 dB.In order to simplify the study, we set ε r = 2.05 (PTFE), ε r = 2.61 (PMMA) and µ r = 1.00 for both materials.
The average value and the standard deviation of the recovered properties, over the frequency range, using the three methods under consideration are in Table 1.We can observe that, in general, the standard deviation decreases as the SNR increases.It does not imply, however, that the average value tends to the nominal one, since the TDG may remove not only the noise, but also some genuine information of the signal.
The worst permittivity value obtained for PMMA, that is, the most different with respect to the nominal one, is ε r = 2.58 (5 dB SNR), which is 1.15% lower than the expected value.The worst case for PTFE is ε r = 2.03, that is 1% below the adopted reference.The worst difference in µ r is 2% for PMMA and 1% for PTFE.
Note that as the method does not restrict the extracted values and the imaginary parts are very small, a variation of the order of 10 −2 easily introduces negative values, without physical meaning.
Table 1.Over the frequency average and standard deviation of the extracted properties for the simulated materials, using the frequency-by-frequency approaches with different SNR.The extracted ε r and µ r of PMMA using the NRW algorithm are in Figure 13, whilst the results for the iterative methods (transmission or reflection) are in Figure 14.It is possible to observe some oscillation on the curves, which is mainly the consequence of using the TDG.The average value over the frequency range and its standard deviation are in Table 2, for different configurations and for both tested materials.For brevity, the results for the PTFE are only in the table.Because the tested materials are low-loss, the relative uncertainty in the imaginary part is pronounced.Actually, it is also coherent with the fact that free-space methods have a limited capability of measuring imaginary parts lower than 0.01 [5,13,27].

PTFE (NRW)
Regarding the two distances, we observe that the results are more consistent when r = 48 cm, especially for PMMA, confirming that the plane-wave premise is quite important in our setup.We should also note that the standard deviations of the extracted properties using experimental data are comparable with those obtained in the simulations, particularly for SNR = 10 or 20 dB.
Table 2.The average and the standard deviation of the extracted parameters of the PTFE and PMMA samples, using the frequency-by-frequency approaches in two distances r.Regarding the extraction methods, the results for PMMA are more uniform compared with those for PTFE.This can be justified by the higher ∆d on the PTFE sample and by its possible bending during the calibration (the PTFE sample is not so rigid as the PMMA one, so it bends slightly when supported only by its borders).

PTFE (NRW)
It must be noted that the results obtained with the experimental data are better than that expected by the uncertainty analysis.This is due to the fact the errors in the actual measurements are lower than the uncertainties declared by the VNA specifications.In this way, the influence of ∆d becomes predominant.Additionally, when performing the free-space calibration more uncertainties may be introduced in the data.Note, for instance, that the shift in the phase depends on knowing d.Considering the discussed aspects, in general, the experimental results are in good agreement with simulations and literature.

Multi-Frequency Reconstruction with Simulated Data
Here we simulate the same scenarios of Section 4.1, but with different weights w in (3).Observe that when using w = 0 or 1, either transmission or reflection is used, respectively.The dependence of the scattering parameters on the searched for ε r , µ r and σ varies with frequency, providing more or less sensitivity for either S 11 or S 21 .In order to avoid ill conditioning we limit 0.1 ≤ w ≤ 0.9.Thus, we always use both transmission and reflection, but with different weights.
The MATLAB function fmincon is used to minimize Ψ, with the bounds 1 ≤ ε r ≤ 10, 1 ≤ µ r ≤ 10 and 0 ≤ σ ≤ 0.1 S/m; the initial point is ε r = µ r = 1 and σ = 0. Observe that local minimum can occur, so the choice of the initial point is crucial as we are using a deterministic algorithm.Alternatively, one can use a global search algorithm, which typically increase the computational cost.
The extracted properties of the theoretical samples are in Figure 15.Intuitively, more reliable results are obtained for the higher SNR.In this case, the weight w has less influence, after all, the signals are good.In contrast, in a lower SNR scenario w is, in general, more influent.

Multi-Frequency Reconstruction with Experimental Data
After the use of theoretical data, we apply the minimization strategy to extract the properties from the experimental S-Parameters, measured with r = 48 cm.In order to compare the approaches and because the experimental noise in data is not necessarily gaussian, here we test the extraction using the raw signals (without gating) and also using the gated ones.
Table 3 summarizes the extracted properties for both samples, for different w, with and without gating.In some cases the results for different w vary in the order of 10 −3 or less, so the rounded values at the table are the same.Such a weak dependence on w indicates that our measured data perform similar to simulated high SNR data.Table 3. Extracted parameters of the samples using the multi-frequency reconstruction with the measurements performed at distance r = 48 cm, using different weights w in (3).We observe that w = 0.9 can give wrong solutions due to local minima presence.A simple inspection of the value of Ψ can help in discarding abnormal solutions.One should, therefore, test different weights w in order to detect the occurrence of a local minimum.

PTFE (without
The extracted properties are quite good and we see that, indeed, the TDG seems to be not necessary when using this approach.This is because, in a certain sense, the minimization of the functional Ψ acts like a filter.Actually, not using the gate in PMMA achieved better results.
The better results achieved by the multi-frequency approach are due to the a priori assumption that the material is non dispersive, so that its properties are constant with frequency and all the measured data, as a whole, contributes to the final results.Conversely, the frequency-by-frequency approach exploits each data individually, so that it has "less information" at disposal.The counterpart, of course, is that it can treat dispersive materials without the need of knowing the model.

Conclusions
In this paper we explored two complementary approaches for extracting the properties of materials using a free-space experimental setup.The presented system avoids the cost of an anechoic environment and the costs of commercial solutions for post-processing and extraction.
The implemented Time-Domain Gating has proven to be an useful tool when using frequencyby-frequency approaches, filtering out much of the interferences that affect the free-space measurements.Some advantages of an TDG implementation are: (a) an algorithm more flexible than proprietary solutions, (b) equipment independent, (c) no costs, and (d) useful for research and academic purposes.
A differential uncertainty analysis has been presented to quantify the influence of measured data and slab thickness uncertainties in the final results.The extracted properties are in good agreement with the theoretical ones.Variations in the experimental results are, in general, consistent with those obtained with simulation, which indicates that the measurement setup can indeed provide good data, despite not having electromagnetic absorbers or a sophisticated positioning system.
Results are consistent within the hundredth, but the relative error in the imaginary part is significant when dealing with low loss materials.This is also consistent with the expectation that such small imaginary parts are difficult to be measured by free-space methods.
The multi-frequency approach has proven to be more robust with respect to noisy data.We validate its applicability to FSM, with non gated experimental data.In general, it is preferable using the two S-Parameters, since it can avoid worse results and under-determination when one of them is low.On the other hand, the frequency-by-frequency approach does not make use of restrictions, so that it can be applied to dispersive materials straightforwardly.
In order to improve the measurement system, some topics for further investigations are: (a) using the additional data S 12 and S 22 to increase the robustness of the measurements; (b) testing dispersive materials in the measurement band; (c) optimizing the gate parameters and the variables d and L 1 .
Since T is complex, the logarithm in (A8) have periodic solutions, that is, ln(T −1 ) = ln |T −1 | + j(φ T + 2πn), where φ T is the angle of T −1 and n should be incremented when the phase wraps [1].

Appendix B. Filtering the S-Parameters with the Time-Domain Gating
Our gating is implemented in MATLAB and the flow chart in Figure A1 illustrates its main steps.Note the circled numbers that guide the reader throughout the process.The process starts in 1 with the input S-Parameter, that are measured at N f points between the frequencies f 1 and f 2 .It is known, however, that the TDG introduces errors concentrated approximately in the first and in the last 10% of the gated signal in FD, that is, in the lower and in the higher edges of the frequency band [28,53].
Therefore, aiming to preserve the original signal from excessive distortions, a wider frequency range is considered for now.To this end, an extrapolation is done in 2 , extending BW as f 1e = f 1 − ∆ f e and f 2e = f 2 + ∆ f e .
The extrapolated signal in 3 has N f e points.The padded values are based on the edge parts of S(ω) with a smooth attenuation.A constant or zero padding does not yield good results [27], since they can abruptly change the tendency of the signal.Moreover, the discontinuity between the edges of the signal causes strong side-lobes in TD.This is known as leakage or Gibbs' phenomena and it is related to the non-periodicity of the signal in the observation interval [54].Because of that, prior to the transform, it is important to use a window to smooth the edges.Fraudulent echoes are then avoided in TD, that is, there is a reduction in the side-lobe levels [1,54].For this reason, in 5 the signal 3 is multiplied by the window 4 , resulting in the weighted signal S W (ω) in 6 .
The transformation to TD is performed in 7 by Inverse CZT (ICZT) [47], with N t points between t 1 and t 2 .It results in 8 , that is s(t).Observation of the signal in time domain allows to specify the limits of the gate time interval g t 1 and g t 2 , in 9 .For instance in Figures 8 and 9 this interval could be chosen as g t 1 = −0.5 ns and g t 2 = 0.5 ns, so that the main lobe of the signal is preserved.The gate limits define the rectangular signal g(t) in 10 , which is 1 between g t 1 and g t 2 , and 0 outside it.An analytical transform is done in 11 , obtaining the sinc G(ω) in 12  [48,49].
Depending on the specifications of the filter, a Kaiser-Bessel window [48,54] is created in 13 and used in 14 to weight G(ω).Because G K (ω) in 15 is the filter in FD, it should be convolved with the signal to be filtered, as in 16 , resulting in the gated one in 17 .Since both convolved signals have the length N f e , the result has 2N f e − 1 samples [48], where the range of interest is the middle portion and it is easily extracted.
In order to minimize the effects introduced by the TDG, an equalization is done, based on the rectangular signal in 18 .The idea is to pass it through the same gate-the convolution in 19 -and using the output 20 to get the effects caused by the gating [30,50].Thus, the gated equalization signal in 20 is compared to the gated S-Parameter by the division in 21 , resulting in the equalized signal in 22 .The process in 23 discards the extrapolated edges, achieving the final result in 24 .At this point, the filtered measured data can be used as input of the extraction algorithms.

Figure 1 .
Figure 1.Representation of the transmission and reflection Free-Space Method.The sample is placed between two antennas, which are connected to the ports of the VNA.

Figure 2 .
Figure 2. Extraction of the permittivity ε r and permeability µ r from the S-Parameters S 11 and S 21 .

Figure 3 .
Figure 3.The sketch (left) and a photograph (right) of the experimental setup.

Figure 4 .
Figure 4.The three configurations used in the free-space calibration: (a) metal plate, (b) air, and (c) material under test.The arrows indicate possible contributions on the measurements.

Figure 5 .
Figure 5. Representation of the ticknesses d and L 1 , which are used to adjust the phase.

Figure 10 .
Figure 10.The uncertainty in the real part of permittivity (∆ε r ) due to the uncertainties in the magnitude and phase of the S-Parameters and in the thickness d.The reflection configuration is on the top, whilst the transmission is on the bottom.

Figure 11 .
Figure 11.The uncertainty in the imaginary part of permittivity (∆ε r ) due to the uncertainties in the magnitude and phase of the S-Parameters and in the thickness d.The reflection configuration is on the top; the transmission is on the bottom.

|| f 2 /Figure 12 .
Figure 12.Total uncertainty over frequency ∆ε r f F RO (solid lines) and ∆ε r f F TO (dashed lines), for different thicknesses d and ε r .The different markers distinguish the thickness d.

Figure 13 .
Figure 13.The properties of PMMA extracted by the NRW using experimental data with r = 48 cm.

Figure 14 .
Figure 14.The permittivity of PMMA extracted using Reflection Only (RO) and Transmission Only (TO) methods, using experimental data with r = 48 cm.

Figure 15 .
Figure 15.The extracted properties of the theoretical samples, for different SNR and weights w.The simulation for PMMA are on the (left), whereas for PTFE are on the (right).