Next Article in Journal
Single-Pixel Background Modeling Algorithm for Strong Sky Scenes Based on Local Region Spatial Bases
Next Article in Special Issue
Over the Limits of Traditional Sampling: Advantages and Issues of AICs for Measurement Instrumentation
Previous Article in Journal
Development of Multilayer Transducer and Omnidirectional Reflection Model for Active Reflection Control
Previous Article in Special Issue
Unlocking All-Solid Ion Selective Electrodes: Prospects in Crop Detection
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Design and Test of a New Dielectric-Loaded Resonator for the Accurate Characterization of Conductive and Dielectric Materials

1
Department of Industrial, Electronic and Mechanical Engineering, Roma Tre University, Via Vito Volterra 62, 00146 Roma, Italy
2
Istituto Nazionale di Fisica Nucleare INFN, Sezione Roma Tre, Via della Vasca Navale 84, 00146 Roma, Italy
*
Author to whom correspondence should be addressed.
Sensors 2023, 23(1), 518; https://doi.org/10.3390/s23010518
Submission received: 21 November 2022 / Revised: 28 December 2022 / Accepted: 28 December 2022 / Published: 3 January 2023
(This article belongs to the Collection Advanced Techniques for Acquisition and Sensing)

Abstract

:
The spread of additive manufacturing techniques in the prototyping and realization of high-frequency applications renewed the interest in the characterization of the electromagnetic properties of both dielectric and conductive materials, as well as the design of new versatile measurement techniques. In this framework, a new configuration of a dielectric-loaded resonator is presented. Its optimization, realization, and use are presented. A measurement repeatability of about one order of magnitude lower than the commonly found values ( 10 3 on the Q-factor and 15 × 10 6 on the resonance frequency, given in terms of the relative standard deviations of repeated measurements) was reached thanks to the design of a closed resonator in which the samples can be loaded without disassembling the whole measurement fixture. The uncertainty levels, the ease of use, and the versatility of the realized system make its use of potential interest in numerous scenarios.

1. Introduction

The characterization at microwave frequencies of the electromagnetic (e.m.) properties of materials has always been a field of great interest due to the relevance that these measurements can have in the design and performances of telecommunication systems. For this reason, a large number of different high-frequency characterization techniques have started to appear in the literature since the 1950s. Several techniques were designed to meet the different measurement needs in terms of the operative band, accuracy, and characteristics of the material under investigation [1,2,3,4,5,6,7,8,9,10,11,12].
In recent years, new materials and manufacturing techniques, of interest also for telecommunication and sensing applications, have emerged. For example, the possibility of printing dielectric structures either with complex shapes or with geometrically controlled e.m. properties was investigated for the realization of radar sensors, antennas, graded-index lenses, etc. [13,14,15,16]. Moreover, printed conductive materials were under study for high-frequency flexible sensing applications, guiding structures, antennas, etc. [17,18,19]. For these reasons, in recent years, the interest in the e.m. characterization of conductive and dielectric materials was revived.
In this work, we describe the design and test of a new application of the Hakki–Coleman configuration of a dielectric-loaded resonator useful for the characterization of both dielectrics and conductors. Thus, the physical quantities of interest in this study are the complex relative permittivity ε ˜ = ε i ε [20] for dielectrics and the surface resistance R s [20] for conductors. The here presented measurement fixture was designed with the aim of obtaining a versatile system for laboratories committed to the prototyping of high-frequency systems also through the use of additive manufacturing techniques. For this reason, the three main design constraints were (i) the need to obtain a system able to work both with conductors and dielectrics, (ii) the ease of use, and (iii) a target uncertainty of ∼1 m Ω for R s measurements on normal conductors and <5% on the real and imaginary part of ε ˜ .
Due to the required accuracy level, a resonant system was selected on purpose. In particular, the dielectric-loaded resonator (DR) fixture was chosen for its high sensitivity, which is in fact generally exploited, and also for the characterization of low-loss materials and even superconductors [2,21,22,23,24,25,26,27,28,29,30,31,32]. However, due to their high sensitivity, microwave resonant techniques generally suffer from a low measurement repeatability, being particularly sensitive to small variations in the sample mounting. For this reason, we designed a DR which does not require a complete disassembling for each measurement: the idea was to produce a closed DR on which the samples under investigation can be mounted from the outside thanks to the presence of windows on the cavity and the dedicated sample holders. In this way, in addition to improving the mounting repeatability, it was possible to obtain an easy-to-use fixture even for non-expert users and to fasten the measurement procedure with respect to other common approaches. Finally, for what concerns the need of characterizing both conductors and dielectrics, the resonator geometry and the sample holders were optimized to guarantee the best sensitivity on R s and ε ˜ considering the properties of the materials of interest.
The paper is organized as follows: In Section 2, the measurement method is introduced; then the design and geometry optimization and realization of the measurement fixture are shown in Section 3; and finally, the experimental use of the realized DR on the R s and ε ˜ measurements is shown in Section 4. A brief summary of the work is provided in Section 6.

2. Measurement Methods

The used technique is based on the measurement of the quality factor Q and the resonance frequency f 0 of a dielectric-loaded resonator (see Figure 1 for a sketch of the structure). Because Q and f 0 depend on the e.m. properties and geometries of every component of the resonator, then the insertion of a sample under test into the resonator will cause a change in Q and f 0 . From the change in these quantities, the e.m. properties of the sample can be, at least in principle, obtained [2].
In particular, it can be shown that the Q factor of a resonator is determined by the components of the resonator itself as follows [2]:
1 Q = i R s , i S i | H τ | 2 d S i 4 π f 0 W + j ε 0 ε j V j | E | 2 d V j 2 W tan δ j : = i R s , i G i + j η j tan δ j ,
where W is the energy stored in the DR at the resonance and H τ is the magnetic field that is tangential to the i-th conductive surface S i with surface resistance R s , i , in the DR. Meanwhile, E is the electric field in the j-th dielectric element of volume V j and loss tangent tan δ j : = ε j / ε j .
Thus, from a measurement of Q, the measurements of R s or tan δ s of a sample under study could be obtained once the R s , i , G i , tan δ j and η j of all the other components of the resonator are known. This would require a complete calibration of the resonator, also taking into account the dependencies of these quantities on all the possible variables of influence (e.g., temperature, pressure, humidity). Because this procedure generally provides large measurement uncertainties, a different strategy, based on a perturbation method, is usually preferred.
Because different measurement procedures are followed to obtain the quantities of interest, i.e., ε ˜ or R s , the measurement methods used for the two classes of materials are described separately in the next subsections.

2.1. Surface Resistance R s Measurement Method

The end-wall perturbation method [2] is used to obtain the variation in R s of the material under investigation with respect to a reference sample of surface resistance R s , r e f :
R s R s , r e f = G s 1 Q 1 Q r e f ,
where Q r e f is the quality factor measured once the reference sample is loaded in the DR. Equation (2) is obtained from Equation (1), assuming that during the two measurements (i.e., that with the sample of unknown R s and that with the reference one R s , r e f ) the change in all the other quantities can be neglected.
The resonance frequency can be used to measure the variations in the surface reactance X s of the material as follows:
X s X s , r e f = 2 G s f 0 f 0 , r e f f 0 , r e f .
Actually, in conventional conductive materials, at the frequency of interest, X s = R s [20,33]; thus, once R s is determined from Q, there is no further need for using f 0 measurements. The use of Q instead of f 0 for R s measurements is motivated by the different sensitivities of Q and f 0 to small variations both in the mounting and positioning of the sample and in influence quantities, such as temperature and pressure. The higher sensitivity of f 0 to all these phenomena generally makes Q measurements more reliable.
The identification of a reliable reference sample goes beyond the scope of this work. Only in principle R s = μ 0 π f ρ , with ρ the dc resistivity of the material. Just as a counter example, it is well known from the literature that the surface roughness R g strongly affects R s [34,35,36,37], but only empirical models R s ( R g ) exist [38,39,40]. Thus, the simple derivation of R s from ρ measurements cannot be considered reliable for the realization of a measurement standard. For this reason, in this work, the metrological performances—in terms of measurement precision and accuracy—of the designed DR will be evaluated on Δ R s measurements and not on absolute R s . However, a procedure for the evaluation of the absolute R s , without the need for a reference, will be presented at the end.

2.2. Complex Permittivity ε ˜ Measurement Method

In this case, because two independent quantities (i.e., ε and ε ) must be measured, both Q and f 0 of the DR are exploited.
The volume perturbation method is used [2]: a part of the inner volume of the DR is substituted with the sample under investigation and then, in the very same volume, the reference dielectric material is loaded. From the variations in Q and f 0 measured in the different configurations, ε ˜ can be obtained. The resonance frequency f 0 can be used to measure ε , while Q, as also shown in Equation (1), is used to obtain tan δ and thus ε , once ε is known. In general, measurements of f 0 with respect to the reference are linked to ε by
ε = g 1 f 0 f 0 , r e f f 0 , r e f ,
where g 1 is a calibration curve that depends on the geometry, the position of the sample in the DR, and on ε r e f . The calibration function g 1 can be obtained through e.m. simulations of the DR.
Once ε is obtained, the filling factor can also be determined using the calibration curve η s = g 2 ( ε ) , similarly obtained by simulations. With η s , Equation (1) can be computed to obtain tan δ , which in the small perturbation limit can be reduced to
tan δ s 1 η s 1 Q 1 Q r e f η s tan δ r e f .
When the small perturbation limit is no longer valid, Equation (1) as a whole must be used, taking into account the variation in G i and η j caused by the insertion of a sample with a different ε with respect to the reference.
A very practical reference can be the air itself: within the typical measurement uncertainties, one can assume ε r e f 1 and tan δ r e f 0 without introducing significant errors.

3. Design and Realization of the DR

The DR was designed starting from these requirements:
  • Improvement in the measurement repeatability: a closed configuration in which the sample can be loaded from the outside without the need of disassembling the whole DR for each measurement is preferred;
  • Possibility of hosting two samples at the same time: this can be used to perform multiple-sample comparisons [41] or to increase the sensitivity when needed;
  • Contactless measurements: the sample holder must be designed to support the samples without letting the probed area of the sample touch other surfaces. This is useful to avoid damaging delicate sample surfaces and/or coatings.
To meet these requirements, a closed structure was designed with open windows on the resonator bases from which the samples can be exposed to the resonator inner volume, as shown in Figure 2. In order to keep the dielectric crystal in position, despite the openings on both bases, a polytetrafluoroethylene (PTFE) holder, pressed into the metallic cavity, was used. The PTFE was chosen because of its good rigidity and relatively low ε and tan δ . In order to reduce the detrimental impact of the PTFE holder on the Q-factor of the DR, its geometry was carefully designed to obtain a small filling factor: the volume occupied by the holder is far from peaks of the E field.
The bases of the resonator were designed to fill the dual role of the sample holder and mask: a rectangular housing was performed to host and center the samples from the outside of the DR. In its center, a circular hole allowed the e.m. mode to probe the surface of the sample without disturbing the mode symmetry. This also allowed the contactless mounting of the sample: in fact, the surface of the sample facing the central hole was not in contact with any surface.
Finally, the resonator was excited in the TE011 mode through coaxial cables (see Figure 2) ended with magnetic loops. The TE011 mode was chosen because, due to the field configuration, it was not necessary to ensure a perfect electrical contact between the bases and the lateral wall. In addition to this, the circular symmetry implies that no electrical contact was needed between the sample and the mask. Moreover, the TE011 is usually well separated from other spurious modes, helping its recognition and avoiding disturbances of near modes.

3.1. Dimensions Optimization— R s Measurements

Once the DR structure was chosen, the optimization of its dimensions to maximize the measurement sensitivity of typical conductive samples was performed through e.m. simulations. We simulated a full 3D structure of the resonant cell, exploiting the Finite Element Method (FEM) with the eigenmode Comsol solver. Because the design of the structure required the optimization of a large set of parameters, we developed a custom external script to automatically identify the TE011 mode-related parameters. The automatically adaptive finer mesh was used for all the structures and an extra fine mesh was used for the dielectric and mask parts.
In particular, the crystal dimensions, the PTFE thickness, and air gaps between the crystal and the samples have to be carefully chosen. The radius of the wall of the metallic enclosure of the DR was simply chosen large enough to make negligible the conduction losses on the lateral wall itself while avoiding other cavity modes close to the TE011. Thus, the cavity radius was fixed at R c a v = 15 mm.
The DR measurement sensitivity c R to the sample surface resistance R s is the parameter to maximize: from Equation (1), c R = | Q / R s | = Q 2 / G s . The optimization was performed by varying the heights of the gaps and analyzing c R given by one available sapphire single-crystal cylinder of dimensions = 8.00 ( 1 ) mm and height h = 5.00 ( 1 ) mm. The simulations were performed setting R = 60 m Ω for all the metallic (brass) surfaces as a realistic surface resistance value [42].
The optimization of the thickness h P T F E of the PTFE holder between the bases and the sapphire crystal was performed by simulating the c R of the DR, changing h P T F E in the interval 0 < h P T F E / mm < 2 , as shown in curve 1 in Figure 3. A maximum, at h P T F E 0.75 mm, exists in c R ( h P T F E , h m a s k ) because of the different sensitivity of Q and G s to h P T F E . With this optimum h P T F E value, the effect of h m a s k on c R was also studied (see curve 2 in Figure 3). It was chosen h m a s k = 0.3 mm as the best compromise between the maximum c R and the minimum thickness reachable to guarantee a sufficiently rigid brass surface usable for the sample holder.

3.2. Dimensions Optimization— ε ˜ Measurements

Once the DR was optimized for R s measurements, the best dimensions of the measurable dielectric samples, and thus of the sample holder, were assessed from the point of view of the sensitivity on the tan δ s measurement. The sensitivity coefficient is obtained from Equation (1) as follows:
c tan δ = Q tan δ s = η s ( i R s , i G i + j η j tan δ j ) 2 = η s Q 2 .
Then, to maximize the sensitivity c tan δ , η s must be optimized depending on the tan δ s of the sample. Thus, a study of c tan δ ( tan δ s , η s ) must be performed.
Defining the quantity b r : = Q 1 η s tan δ s , i.e., a measure of the overall losses apart from those due to the sample, it is possible to identify two different asymptotic limits for c tan δ , one in which η s tan δ s b r and the second for η s tan δ s b r . In the first case, one obtains | c tan δ | η s 1 tan δ s 2 , while in the second | c tan δ | η s b r 2 which, as expected, is no longer dependent on tan δ s . In other words, when the losses in the dielectric sample are large with respect to the other losses into the resonator, then the measurement sensitivity decreases due to the lowering of the Q factor: in this case, a smaller sample is useful to reduce η s . Whereas, when the losses in the sample are smaller than the others, then the sensitivity of the resonator is limited by its intrinsic Q: in this case, a bigger sample is useful to increase the amount of losses in the sample. Hence, the best η s , o p t for a fixed tan δ s is obtained at the crossover of these two opposite scenarios, that is, when the losses in the sample are the same as those in the cavity, η s , o p t tan δ s / b r = 1 , which can be analytically obtained from Equation (6). The corresponding maximum sensitivity is | c tan δ | m a x = 1 / ( 4 b r tan δ s ) .
Thus, the best dimension of the sample can be evaluated starting from b r and a rough estimation of tan δ s . From these, the optimum filling factor is evaluated η s , o p t = b r / tan δ s . Finally, through e.m. simulators, the dimensions of the sample can be assessed in order to obtain η s , o p t using an expected value for ε . Hence, there is not an optimum configuration for every kind of dielectric sample, but the design of the DR should remain versatile enough to accommodate samples of different dimensions. This means that a set of different sample holders can be designed and realized in order to accommodate samples with different e.m. properties, and thus to always exploit the maximum measurement sensitivity of the resonator. In particular, sample holders with holes of different diameters can be useful for this application. In addition to this, it is also possible to act on the thickness of the samples when they are specifically prepared for characterization.

3.3. Realization of the DR

The so-optimized design was realized by also taking into account constraints regarding the cost, hardness, and minimum workable dimensions of the different materials. The final dimensions of the optimized configuration are reported in Figure 4 and pictures of the realized DR are shown in Figure 5.

4. Experimental Tests and Performances Analysis

In this section, the metrological characteristics of the realized resonator are first experimentally assessed, and then the examples of the use of the DR for the characterization of both the conductive and dielectric materials are provided.
The transmission and reflection scattering S-parameters are measured with an Anritsu 37269D Vector Network Analyzer (VNA), setting the center frequency on the peak of the resonance, a span of ∼5 times the full-width half-maximum (FWHM) (as shown in [43,44,45] for an improved measurement accuracy) and 1601 evenly spaced points (the maximum allowed number of points for a frequency sweep). The emitted power is set at −10 dBm and the Intermediate Filter (IF) bandwidth at 10 kHz. The DR is linked to the VNA with 60 cm long phase-stable semi-rigid K coaxial cables, Anritsu 3670KF50-2. The transmission line is calibrated with the standard Short-Open-Load-Through (SOLT) calibration procedure before performing the measurement. At this calibration plane, the DR is connected through 10 cm semi-rigid K coaxial cables ended with the couplers into the DR cavity.
From the transmission scattering parameters, the loaded Q l and f 0 are evaluated, fitting the resonance curves with the following modified Lorentzian model [46,47,48]:
A 1 + i 2 Q l f f 0 f 0 + B e i ( α + f β ) ,
where B is a complex constant representing the cross-coupling between the two ports, A the peak of the unperturbed Lorentzian curve, and α + f β the phase delay given by the uncalibrated final part of the transmission line.
The coupling is set to be so low to allow the approximation Q Q l . In the case of non-negligible coupling, the unloaded Q is evaluated through the TMQF algorithm [11,49].
An example of the measurement of the transmission S21 parameter around the resonance frequency is shown in Figure 6. The position of the out-of-resonance points shows the small cross-coupling of the DR, B 0.00026 + i 0.00002 .

4.1. Measurement Repeatability

Due to the typical high sensitivity of the DR on the mounting, the measurement repeatability performances of the DR were tested by repeating measurements disassembling and re-mounting a conductive sample each time. The repeatability test was performed in the single-sample configuration (closing the lower base with the brass cap) and using a copper sample of nominal dimensions 15 × 15 × 3 mm 3 . A mass of ∼200 g was placed on the sample to improve the repeatability.
In Figure 7, the Q and f 0 measurement repetitions are shown. With this new DR design, it was possible to reduce the Q mounting repeatability (evaluated as the standard deviation of the experimental points) from a typical 5% [25] to 0.1%. Meanwhile, the resonance frequency repeatability is here assessed to be ∼ 15 × 10 6 . These are the largest contributions on u ( Q ) and u ( f 0 ) . In fact, the uncertainties u ( Q ) and u ( f 0 ) evaluated by the fitting algorithm from the residuals are more than an order of magnitude lower than those provided by the mounting repeatability. In addition to this, u ( Q ) and u ( f 0 ) obtained by the fitting procedure can be easily reduced by narrowing the IF bandwidth or performing an average on the acquisition of the single points in the frequency, whereas the contributions of the mounting are not improved by this. Thus, in other words, the mounting repeatability is an intrinsic characteristic of the designed measurement fixture, whereas the other contributions on u ( Q ) and u ( f 0 ) depend on the measurement procedure and instrumentation used.

4.2. R s Measurements and Uncertainty Evaluation

The performances of the designed resonator on the measurements of R s are here reported. However, because the measurement of the absolute R s requires a complete calibration of the resonator, an accurate measurement of all R s , i and tan δ j (see Equation (1)) is necessary, in addition to the accurate determination of the geometrical and filling factors. For this reason, to highlight the proper metrological characteristics of the DR to R s measurements, we first show the application to differential Δ R s measurements and then to absolute R s measurements. The uncertainty affecting this last type of measurement is in fact largely determined by the uncertainties on the e.m. properties of the components of the resonator: these could improve with the development of novel measurement techniques, and they are not determined by the design of the here presented measurement fixture.
In the next sections, the measurement performances of the realized DR will be experimentally tested both on the Δ R and R measurements. The samples used for this evaluation are reported in Table 1.

4.2.1. Differential Δ R s Measurement

In several cases, one is not interested in the absolute R s value but instead in the differences Δ R s between the samples in order to evaluate, for example, which surface treatment gives the best results. In this case, no calibration of the resonator is needed so that the uncertainties u ( Δ R s ) can be very small.
To lighten the mathematical notation, we define the sum of the whole dielectric losses as l d : = j η j tan δ j , the inverse of the unloaded quality factor as l : = Q 1 , and that of the geometrical factors as A i : = G i 1 .
Because with differential Δ R s measurements l d and the geometrical factors can be assumed constant among the different measurements (if the difference Δ R s in the samples is not so large to change the field configuration and negligible variations in the measurement influence quantities are present), one obtains:
Δ R s = R s , i R s , i i = l i l i i A s ,
where l i and l i i are the inverses of the quality factors measured when the samples with surface resistance R s , i and R s , i i are loaded.
The measurements are performed with a brass mask with 1 = 13.00 ( 1 ) mm with a geometrical factor G 1 = 10.2 ( 1 ) k Ω . From this, the geometrical factor of the sample is evaluated, obtaining G s 1 = 1.31 ( 1 ) k Ω . The number in parentheses is the numerical value of the standard uncertainty u ( G ) referred to the corresponding last digits of the quoted parameter. u ( G ) are evaluated by the e.m. simulations, taking into account the uncertainties on the physical dimensions of the resonator u ( h c a v ) / h c a v = 0.2 % , where h c a v is the height of the cavity, the uncertainties on the dielectric crystal relative permittivity ( 9 < ε < 10 ), and the variation in the effective length of the cavity given by the field penetration length in Cu at 13 GHz ( δ C u , 13 GHz 0.6  μm). The contribution to the geometrical factors uncertainties given by these three effects are, respectively, u ( G ) d i m / G 0.5 % , u ( G ) ε / G 1 % , and u ( G ) δ C u , 16 GHz / G 0.3 % . Thus, the main contribution is given by u ( ε ) of the sapphire crystal. The uncertainties on the physical dimensions of the cavity and of the mask are the mechanical tolerances of the numerical control tool used for manufacturing.
Taking into account the type-A evaluation of the measurement uncertainty reported in Section 4.1, the obtained Δ R s measurements are shown in Table 2:
As expected, the uncertainties are “small” in this case because those are determined only by u ( Q ) and u ( G i ) . Once again, this highlights the importance an R s measurement standard, currently missing, would have in microwave measurements. Thanks to its high measurement repeatability, the DR here presented allows discriminating samples with Δ R 2 m Ω which is a sufficiently accurate measurement for standard microwave applications.

4.2.2. Absolute R s Measurement

Because it is well known that the actual R s of metallic samples is often far from the nominal R s value evaluated from the dc resistivity ρ [50,51], an accurate calibration cannot rely on an analytically computed R s . Moreover, for l d , the use of the literature values is not recommended due to the large dispersion of the tan δ values associated with sapphire and different kinds of PTFE. For these reasons, an in situ calibration procedure was developed, exploiting the properties of the new DR.
The aim is to find a combination of independent measurements that allows the determination of l d and the surface resistance R 0 of the brass bases of the resonator (the contribution of the aluminum lateral wall on the Q is negligible in this configuration), once the geometrical factors are estimated by numerical simulations (indeed, there is no way to measure geometrical factors without a calibrated standard [25]). The minimum set of independent measurements is obtained by performing the first two measurements using two pairs of samples of different materials (one of them of the same material as the cavity R 0 ) and performing the third measurement by changing the DR geometrical factors. In this way, there is no need to use a third unknown material, and l d can be determined. The geometrical factors of the resonator can be changed using masks with different diameters of the central hole, through which the samples are exposed to the DR.
Mask M i ( i = 1 , 2 ) has a hole diameter i , a geometrical factor G i = A i 1 , and exposes a sample area of the geometrical factor G s i = A s i . Due to the structure symmetry, the bases of the resonator are considered equivalent. The set of measurements is represented by the following system:
2 A s 1 2 A 1 1 0 2 A s 1 + 2 A 1 1 A s 2 A s 1 + A 1 + A 2 1 · R x R 0 l d = l 1 l 2 l 3 ,
where R x is the surface resistance of an unknown sample and the conduction losses of the lateral wall of the resonator are neglected. The system (9) allows to find R 0 , l d , and R x .
The following calibration procedure was used with two pairs of brass and copper samples. The masks are entirely made of brass and the DR is used in the dual sample configuration (i.e., with samples mounted on both bases). One of these is that used in the previous measurement (see Section 4.2.2 for the dimensions), while the second one has the central hole of diameter 2 = 9.00 ( 1 ) mm and a geometrical factor G 2 = 2.84 ( 3 ) k Ω . The samples’ geometrical factors are then G s 1 = 1.31 ( 1 ) k Ω (with the first mask) and G s 2 = 1.96 ( 2 ) k Ω (with the second mask).
By solving system (9), the measured volume losses are l d = 4 ( 2 ) × 10 5 and R 0 = 92 ( 12 ) m Ω . Assuming that the sapphire crystal filling factor is ∼1 (in dielectric-loaded resonators, thanks to the high dielectric permittivity ε of the crystal used, the electromagnetic field is confined in the volume of the crystal itself; thus, the energy stored in the resonator is in the first approximation that is contained in the dielectric crystal [52]), then l d tan δ s a p . The obtained tan δ s a p of the sapphire crystal is well in agreement with the room temperature value reported in [53]. Moreover, R 0 is compatible with that of brass, considering that the surface roughness can double [38] the nominal R evaluated from the dc resistivity of the material (see Table 1).
The DR is then tested with the samples in Table 1. The results obtained with the single-sample configuration are reported in Table 3.
The large uncertainties are given in this case by the large uncertainty on l d . It must be noticed that the measured R seems to be particularly far from the ideal values shown in Table 1. This evident discrepancy can be interpreted according to the empirical models [54] that describe the dependence of R on the surface roughness R g of the sample. R g is estimated to be in the range ( 0.9 ÷ 3.0 )   μ m through these same models which is a roughness level well compatible with what is expected on the used samples.

4.3. ε ˜ Measurements and Uncertainty Evaluation

The versatility of the designed measurement fixture allows its use both for the R s measurements of conductive samples and for the ε ˜ of dielectric materials. In this section, we show an example of the use of the designed DR for the characterization of dielectric materials.
In order to experimentally test the DR in a wide ε ˜ space, we exploited the possibility of geometrically controlling ε ˜ by just printing dielectric samples with different amounts of vacuum inside. The samples were printed with a high-temperature photopolymer material, using the PolyJet deposition technique. The high spatial resolution of the printer allowed the realization of artificially porous samples. The effective permittivity ε ˜ e f f of the samples was controlled by printing empty columns across the whole sample thickness. The columns are arranged on square or hexagonal lattices: varying the lattice parameter l p and the column diameter h , samples with different filling percentages are obtained. The samples are electromagnetically homogeneous because l p λ with λ the wavelength in the medium. At ∼13 GHz and with ε 3 (anticipating the results shown in Figure 8), λ 13 mm. A sketch of the samples is shown in Figure 9, and the characteristics of the samples are reported in Table 4.
The following Section 3.2 samples of nominal thickness t = 1.75 mm were printed to match the maximum measurement sensitivity, considering the mounting of the sample on the mask with the 13 mm central hole. The variations in the resonance frequency Δ f 0 = f 0 , s f 0 , r e f measured with or without a loaded sample, and considering air as a reference, are simulated for different ε values and shown in Figure 10. Hence, the calibration curve ε ( Δ f 0 ) is obtained.
ε ( Δ f 0 ) = 151.7 ( 1.6 ) Δ f 0 2 67.50 ( 8 ) Δ f 0 + 1.0020 ( 10 ) ,
where the numbers in parentheses are the numerical values of the standard uncertainties referred to the corresponding last digits of the quoted parameters. For each of these points, the sample filling factor η s is computed (see Figure 10) and the corresponding calibration curve η s ( Δ f 0 ) obtained:
η s ( Δ f 0 ) = 0.4620 ( 5 ) Δ f 0 2 0.17390 ( 5 ) Δ f 0 + 0.0023610 ( 5 ) ,
The simulated points used to obtain the 2nd-order polynomial calibration curves are shown in Figure 10.
So, once Δ f 0 are measured, with Equations (10) and (11), ε and η s of the sample can be obtained.
In order to measure the resonance frequency variation Δ f 0 with respect to an air-reference sample, a ring of the same thickness as the samples is printed in order to compare f 0 and Q, keeping the upper metallic closing cap at the same height and thus not changing the geometry of the DR. This ring has a central circular hole with an inner diameter = ( 14.00 ± 0.01 ) mm, and thus larger than the mask hole, in order not to interfere with the e.m. magnetic field. Δ ( Q 1 ) and Δ f 0 are then determined by measuring Q and f 0 of the DR loaded with the full square samples and with the ring.
The measurement procedure and uncertainty evaluation are now described: The VNA IF bandwidth was set to 10 kHz and every point was averaged over five acquisitions. The frequency span was chosen as ∼5 FWHMs [43,44,45] of the resonance curve. For each sample, the mounting was repeated 10 times. All the measurements were performed after a 12-terms VNA calibration performed with SOLT standards. Each sample (including the reference rings) was loaded into the DR in place of the metallic sample and then covered with a brass cap. The thickness t of every sample was measured with a micrometer, and the standard deviation of 10 repeated measurements was used as the uncertainty u ( t ) .
The f 0 standard deviation of the 10 different mountings is always u m ( f 0 ) / f 0 < 10 6 . Because in real cases the height of the DR changes when the reference or the sample is mounted, the uncertainty u ( t ) was taken into account and propagated on the overall u ( f 0 ) through the sensitivity function f 0 / t evaluated with numerical simulations. Thus, u 2 ( f 0 ) = u m 2 ( f 0 ) + u t 2 ( f 0 ) where u t ( f 0 ) = u ( t ) f 0 / t . Then, Δ f 0 = f 0 , R f 0 , S between the samples (subscript s) and rings (subscript R) of the same nominal thickness is evaluated and the uncertainty is obtained with the standard propagation procedure u 2 ( Δ f 0 ) = u 2 ( f 0 , R ) + u 2 ( f 0 , s ) .
Once f 0 and u ( f 0 ) are correctly evaluated, these can be used to evaluate ε and u ( ε ) with Equation (10) and η , u ( η ) with Equation (11). Finally, using Equation (5), tan δ s , and hence ε , and their uncertainties, can be evaluated. The results are reported in Figure 8.
The uncertainties on the filling percentage of the samples are obtained starting from the 3D printer spatial accuracy and then analytically propagated using the geometrical properties of the lattices shown in Table 4.
Both ε and ε show a linear dependence on the filling percentage γ of the sample in agreement with the effective medium theory and the upper limit of the Wiener model [55,56] and other literature works [57,58,59].
The results are also in agreement with those shown in other works. In [60], acrylonitrile butadiene styrene (ABS) samples doped with different quantities of BaTiO3 microparticles were measured with a split-post resonator at 15 GHz, obtaining 2.6 < ε < 8.7 and 0.005 < tan δ < 0.027 , whereas a broadband characterization (from 1 MHz to 11 GHz) was performed in [61]. The high-frequency range ( 8.2 ÷ 11 ) GHz was analyzed with the Nicholson–Ross–Weir reflection method [61]. Even if no measurement uncertainties are indicated, the results show 2.50 < ε < 3.29 and 0.005 < tan δ < 0.037 .
Finally, an experimental validation of the measurement accuracy of the new DR configuration on ε was shown in [62], comparing the ε measurements obtained with this resonant set-up and with a transmission/reflection standard technique based on the so-called “NIST-precision” method [63], which is an improved version of the Nicholson–Ross–Weir, and using a WR90 waveguide. With this set-up, we obtained ε 3.1 [62], which is fairly in agreement with the here shown value measured on the full sample ε = 2.9 ( 1 ) .

5. Comparison with the State of the Art

The measurement fixture and measurement methods shown in this paper are the results of careful work, carried out in the last years, oriented to the development of microwave measurement techniques and systems for the characterization of materials (both conductors and dielectrics). In particular, and differently from many other works, we put a strong focus on the analysis and improvement in the metrological characteristics. In fact, it is clear how materials characterization methods at microwave frequencies were an evolving field in the last years: these techniques had great success in the experimental investigation of the e.m. properties of matter [64,65,66,67,68], and nowadays, microwave measurements are instrumental both for the design and testing of devices and complex systems whose performances can be particularly sensitive to the so-measured quantities [24,69,70,71,72,73,74,75,76] and in metrology [77,78,79]. This is the general framework into which this work fits: it shows a novel configuration of a dielectric-loaded resonator, and its possible applications, with a thorough evaluation of the measurement uncertainty.
Because one of the most detrimental sources of uncertainty in microwave resonant systems is the generally poor mounting repeatability, the presented fixture was specifically thought to overcome this point. Thus, a good figure of merit for this comparison can be the uncertainty u ( Q ) on the quality factor measurement. From the repetitions shown in Figure 4, we assessed that the standard relative uncertainty on the Q measurements for this fixture is u ( Q ) / Q = 0.1 % , including also the mounting repeatability. This can be compared with the “few” percent declared in [25], the 1% in [80], or the 4% in the international standard [81], even if in all these works no information about the repeatability is provided. We can conclude that with this new design, we improved the uncertainty on the Q-factor measurement of about one order of magnitude with respect to other literature works.
The complete description of the resonator and of its metrological characteristics here reported were not presented elsewhere, although some potential of the fixture was anticipated. Some preliminary results, based on differential surface resistance measurements, were previously reported in [82]. We recently proposed the use of this DR as a contactless surface roughness measurement fixture [83], taking advantage of the high measurement repeatability. For what concerns its use for the characterization of dielectric materials, we previously showed a preliminary optimization of the fixture in [62], which was then adapted to the configuration here presented, that was more extensively tested in this paper.
A discussion of the general metrological performances on the measurements of the physical quantities reported in this paper is in order. The relative measurement uncertainties u ( ε ) / ε 3 % and u ( tan δ ) / tan δ 10 % obtained with this resonator are aligned with other works: u ( ε ) / ε 1 % and u ( tan δ ) / tan δ 20 % are declared in [84] on ABS samples at 10 GHz, u ( ε ) / ε 1 % and u ( tan δ ) / tan δ 0.2 % (but simply and incompletely evaluated as the standard deviation of six measurements) on doped ABS samples in [60], and u ( ε ) / ε 5.8 % and 10 % < u ( tan δ ) / tan δ 0.2 % < 200 % on 3D printed ABS with the Nicholson–Ross–Weir waveguide method in [61]. Looking at measurements of ε ˜ only, the method here presented does not significantly improve the uncertainties with respect to other methods. However, the novelty of this method is still of interest in several applications due to the ease of use and to the possibility of characterizing with the same fixture both dielectric materials (even in the presence of back metal plates) and conductors. A planned application of this system is the measurement of the surface impedance/complex permittivity of bad conductors, such as conductive paints or a semiconductor.
Finally, it is worth noting that the uncertainties shown in the previous sections also take advantage of the optimization on the Q and f 0 measurement procedure and algorithm, based on the complex Lorentzian fit of the transmission S-parameters, that we extensively discussed in [43].

6. Summary

Due to the rapid spread of additive manufacturing techniques, and also, for perspective, high-frequencies applications, the interest in microwave measurement systems able to characterize the e.m. properties of these materials is increasing. For this reason, we designed a new measurement fixture, based on a dielectric-loaded resonator DR, able to measure either the surface resistance R s of conductive samples or the complex permittivity ε ˜ of dielectric samples. The DR was designed to obtain a measurement system that is easy to use and versatile enough to quickly test the properties of different materials used in prototyping laboratories: in the proposed configuration, the samples can be loaded from the outside of the resonant cell without the need to disassemble the whole resonator for each measurement. This feature also allowed to enhance the typical low measurement repeatability provided by this kind of fixture.
In this paper, the optimization procedure followed to maximize the DR measurement sensitivity is discussed and the final geometry is reported. The designed DR was then realized and its metrological characteristics, in terms of the measurement repeatability and uncertainties, were carefully evaluated. The relative standard deviations of the Q and f 0 measurements obtained in 20 repeated mountings are 10 3 and 15 × 10 6 , respectively. These are about one order of magnitude smaller than the typical values obtainable with this kind of fixture. The realized DR was then experimentally tested on the measurement of R s of different conductive samples and on the ε ˜ measurement of dielectric samples. In particular, to highlight the performances of the designed fixture in terms of the effect of the excellent measurement repeatability, the DR was first used for the measurement of the R s variation among different samples, using one of these as the reference. This approach allowed us to evaluate the uncertainty u ( Δ R s ) = 1.1 m Ω which is given only by the uncertainties u ( Q ) and u ( G ) . In the text, a procedure for an in situ calibration of the resonator is shown to obtain absolute R s measurements. With this procedure, a relative standard uncertainty of 20% was obtained on copper. This highlights the importance that an R s standard would have in microwave measurements.
For what concerns the characterization of dielectric materials, the use of the designed DR for ε ˜ measurements of 3D printing materials was shown. Due to the large variability that the loss tangent tan δ s can have among different materials, here an optimization of the dimensions of the sample to be measured must be performed, starting from an a priori estimation of ε ˜ (if necessary, this optimization can be refined in successive steps). In this paper, the complex permittivity ε ˜ of 3D printed samples with different porosity was measured. The linear dependence of both the real and imaginary parts of ε ˜ found on the filling amount of the samples was shown to be in good agreement with the theoretical models. Thus, in addition to the good agreement between the obtained measurements and the literature values, that confirms the accuracy of the developed system also in the ε ˜ measurements.
In summary, we developed a versatile microwave measurement fixture for the characterization of conductive and dielectric materials. The measurement uncertainties provided, the simple conceptual approach, and the ease of use make this fixture of potential interest for test laboratories involved in the design and prototype of microwave components and systems.

Author Contributions

Conceptualization, A.A., K.T., N.P. and E.S.; methodology, A.A. and K.T.; validation, A.A. and K.T.; formal analysis, A.A.; investigation, A.A. and K.T.; resources, A.A.; data curation, A.A., K.T. and P.V.G.; writing—original draft preparation, A.A.; writing—review and editing, A.A., K.T., P.V.G., N.P. and E.S.; visualization, A.A., K.T. and P.V.G.; supervision, N.P. and E.S.; project administration, E.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data are available on request to the corresponding author.

Acknowledgments

The authors warmly thank the “C.M. srl—Costruzioni Meccaniche Latina” for the support in the realization of the dielectric-loaded resonator.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kaatze, U. Techniques for measuring the microwave dielectric properties of materials. Metrologia 2010, 47, S91. [Google Scholar] [CrossRef]
  2. Chen, L.F.; Ong, C.; Neo, C.; Varadan, V.; Varadan, V.K. Microwave Electronics: Measurement and Materials Characterization; John Wiley & Sons: Hoboken, NJ, USA, 2004. [Google Scholar]
  3. Courtney, W.E. Analysis and evaluation of a method of measuring the complex permittivity and permeability microwave insulators. IEEE Trans. Microw. Theory Tech. 1970, 18, 476–485. [Google Scholar] [CrossRef]
  4. Krupka, J. Contactless methods of conductivity and sheet resistance measurement for semiconductors, conductors and superconductors. Meas. Sci. Technol. 2013, 24, 062001. [Google Scholar] [CrossRef]
  5. Adam, S.F. A new precision automatic microwave measurement system. IEEE Trans. Instrum. Meas. 1968, 17, 308–313. [Google Scholar] [CrossRef]
  6. Karbowiak, A. An instrument for the measurement of surface impedance at microwave frequencies. Proc. IEE-Part B Radio Electron. Eng. 1958, 105, 195–203. [Google Scholar] [CrossRef]
  7. Roussy, G.; Felden, M. A sensitive method for measuring complex permittivity with a microwave resonator. IEEE Trans. Microw. Theory Tech. 1966, 14, 171–175. [Google Scholar] [CrossRef]
  8. Harris, F.E.; O’Konski, C.T. Measurement of High Permittivity Dielectrics at Microwave Frequencies. Rev. Sci. Instrum. 1955, 26, 482–485. [Google Scholar] [CrossRef]
  9. Krupka, J. An accurate method for permittivity and loss tangent measurements of low loss dielectric using TE/sub 01 delta/dielectric resonators. In Proceedings of the 1988 Fifth International Conference on Dielectric Materials, Measurements and Applications, Canterbury, UK, 27–30 June 1988; pp. 322–325. [Google Scholar]
  10. Afsar, M.N.; Birch, J.R.; Clarke, R.; Chantry, G. The measurement of the properties of materials. Proc. IEEE 1986, 74, 183–199. [Google Scholar] [CrossRef]
  11. Leong, K.; Mazierska, J. Precise measurements of the Q factor of dielectric resonators in the transmission mode-accounting for noise, crosstalk, delay of uncalibrated lines, coupling loss, and coupling reactance. IEEE Trans. Microw. Theory Tech. 2002, 50, 2115–2127. [Google Scholar] [CrossRef]
  12. D’Alvia, L.; Pittella, E.; Rizzuto, E.; Piuzzi, E.; Del Prete, Z. A portable low-cost reflectometric setup for moisture measurement in cultural heritage masonry unit. Measurement 2022, 189, 110438. [Google Scholar] [CrossRef]
  13. Sene, B.; Reiter, D.; Knapp, H.; Pohl, N. Design of a cost-efficient monostatic radar sensor with antenna on chip and lens in package. IEEE Trans. Microw. Theory Tech. 2021, 70, 502–512. [Google Scholar] [CrossRef]
  14. Nassar, I.T.; Weller, T.M.; Tsang, H. 3-D printed antenna arrays for harmonic radar applications. In Proceedings of the WAMICON 2014, Tampa, FL, USA, 6 June 2014; pp. 1–4. [Google Scholar]
  15. Wirth, S.G.; Morrow, I.L. Near-field microwave sensor composed of 3D printed antennas and lenses. In Proceedings of the 2018 IEEE MTT-S International Conference on Numerical Electromagnetic and Multiphysics Modeling and Optimization (NEMO), Reykjavik, Iceland, 8–10 August 2018; pp. 1–4. [Google Scholar]
  16. Harter, M.; Hildebrandt, J.; Ziroff, A.; Zwick, T. Self-calibration of a 3-D-digital beamforming radar system for automotive applications with installation behind automotive covers. IEEE Trans. Microw. Theory Tech. 2016, 64, 2994–3000. [Google Scholar] [CrossRef]
  17. Zhang, B.; Zirath, H. Metallic 3-D printed rectangular waveguides for millimeter-wave applications. IEEE Trans. Compon. Packag. Manuf. Technol. 2016, 6, 796–804. [Google Scholar] [CrossRef]
  18. Patel, A.M.; Grbic, A. A printed leaky-wave antenna based on a sinusoidally-modulated reactance surface. IEEE Trans. Antennas Propag. 2011, 59, 2087–2096. [Google Scholar] [CrossRef]
  19. Castelli Gattinara Di Zubiena, F.; D’Alvia, L.; Del Prete, Z.; Palermo, E. A static characterization of stretchable 3D-printed strain sensor for restoring proprioception in amputees. In Proceedings of the 2022 IEEE International Conference on Flexible and Printable Sensors and Systems (FLEPS), Vienna, Austria, 10–13 July 2022; pp. 1–4. [Google Scholar]
  20. Jackson, J.D. Classical Electrodynamics; John Wiley & Sons: Hoboken, NJ, USA, 2007. [Google Scholar]
  21. Baker-Jarvis, J.; Geyer, R.G.; Grosvenor, J.; Janezic, M.D.; Jones, C.A.; Riddle, B.; Weil, C.M.; Krupka, J. Dielectric characterization of low-loss materials a comparison of techniques. IEEE Trans. Dielectr. Electr. Insul. 1998, 5, 571–577. [Google Scholar] [CrossRef]
  22. Krupka, J.; Klinger, M.; Kuhn, M.; Baryanyak, A.; Stiller, M.; Hinken, J.; Modelski, J. Surface resistance measurements of HTS films by means of sapphire dielectric resonators. IEEE Trans. Appl. Supercond. 1993, 3, 3043–3048. [Google Scholar] [CrossRef]
  23. Mazierska, J. Dielectric resonator as a possible standard for characterization of high temperature superconducting films for microwave applications. J. Supercond. 1997, 10, 73–84. [Google Scholar] [CrossRef]
  24. Romanov, A.; Krkotić, P.; Telles, G.; O’Callaghan, J.; Pont, M.; Perez, F.; Granados, X.; Calatroni, S.; Puig, T.; Gutierrez, J. High frequency response of thick REBCO coated conductors in the framework of the FCC study. Sci. Rep. 2020, 10, 1–12. [Google Scholar] [CrossRef]
  25. Mazierska, J.; Wilker, C. Accuracy issues in surface resistance measurements of high temperature superconductors using dielectric resonators (corrected). IEEE Trans. Appl. Supercond. 2001, 11, 4140–4147. [Google Scholar] [CrossRef]
  26. Lee, S.Y.; Soh, B.; Ahn, J.; Cho, J.; Park, B.; Jung, C.; Fedorov, V.; Denisov, A.; Kim, Y.; Hahn, T.; et al. Use of a dielectric-loaded cylindrical cavity in measurements of the microwave surface resistances of high-T/sub c/superconducting thin films. IEEE Trans. Appl. Supercond. 1997, 7, 2013–2017. [Google Scholar]
  27. Kobayashi, Y.; Kayano, H. An improved dielectric resonator method for surface impedance measurement of high-T/sub c/superconductors. In Proceedings of the 1992 IEEE MTT-S Microwave Symposium Digest, Albuquerque, NM, USA, 1–5 June 1992; pp. 1031–1034. [Google Scholar]
  28. Pompeo, N.; Torokhtii, K.; Alimenti, A.; Sylva, G.; Braccini, V.; Silva, E. Pinning properties of FeSeTe thin film through multifrequency measurements of the surface impedance. Supercond. Sci. Technol. 2020, 33, 114006. [Google Scholar] [CrossRef]
  29. Pompeo, N.; Alimenti, A.; Torokhtii, K.; Bartolomé, E.; Palau, A.; Puig, T.; Augieri, A.; Galluzzi, V.; Mancini, A.; Celentano, G.; et al. Intrinsic anisotropy and pinning anisotropy in nanostructured YBa2Cu3O7−δ from microwave measurements. Supercond. Sci. Technol. 2020, 33, 044017. [Google Scholar] [CrossRef]
  30. Pompeo, N.; Alimenti, A.; Torokhtii, K.; Sylva, G.; Braccini, V.; Silva, E. Microwave properties of Fe(Se,Te) thin films in a magnetic field: Pinning and flux flow. J. Phys. Conf. Ser. 2020, 1559, 012055. [Google Scholar] [CrossRef]
  31. Alimenti, A.; Pompeo, N.; Torokhtii, K.; Spina, T.; Flükiger, R.; Muzzi, L.; Silva, E. Microwave measurements of the high magnetic field vortex motion pinning parameters in Nb3Sn. Supercond. Sci. Technol. 2020, 34, 014003. [Google Scholar] [CrossRef]
  32. Pompeo, N.; Alimenti, A.; Torokhtii, K.; Silva, E. Physics of vortex motion by means of microwave surface impedance measurements. J. Low Temp. Phys. 2020, 46, 343–347. [Google Scholar] [CrossRef]
  33. Hagen, E.; Rubens, H. Über Beziehungen des Reflexions-und Emissionsvermögens der Metalle zu ihrem elektrischen Leitvermögen. Ann. Phys. 1903, 316, 873–901. [Google Scholar] [CrossRef] [Green Version]
  34. Tischer, F. Effect of surface roughness on surface resistance of plane copper surfaces at millimetre waves. In Proceedings of the Institution of Electrical Engineers; IET: London, UK, 1974; Volume 121, pp. 333–336. [Google Scholar]
  35. Morgan, S.P., Jr. Effect of surface roughness on eddy current losses at microwave frequencies. J. Appl. Phys. 1949, 20, 352–362. [Google Scholar] [CrossRef] [Green Version]
  36. Ansuinelli, P.; Schuchinsky, A.G.; Frezza, F.; Steer, M.B. Passive intermodulation due to conductor surface roughness. IEEE Trans. Microw. Theory Tech. 2018, 66, 688–699. [Google Scholar] [CrossRef]
  37. Garcia, C.; Rumpf, R.; Tsang, H.; Barton, J. Effects of extreme surface roughness on 3D printed horn antenna. Electron. Lett. 2013, 49, 734–736. [Google Scholar] [CrossRef]
  38. Saad, T.S. Microwave Engineers’ Handbook; Artech House Publishers: London, UK, 1971; Volume 2. [Google Scholar]
  39. Biot, M.A. Some new aspects of the reflection of electromagnetic waves on a rough surface. J. Appl. Phys. 1957, 28, 1455–1463. [Google Scholar] [CrossRef]
  40. Gu, X.; Tsang, L.; Braunisch, H. Modeling effects of random rough interface on power absorption between dielectric and conductive medium in 3-D problem. IEEE Trans. Microw. Theory Tech. 2007, 55, 511–517. [Google Scholar] [CrossRef]
  41. Alimenti, A.; Torokhtii, K.; Pompeo, N.; Silva, E. Sensitivity limits comparison of surface resistance measurements based on dielectric loaded resonators. J. Phys. Conf. Ser. 2018, 1065, 052029. [Google Scholar] [CrossRef]
  42. Rossiter, P.L. The Electrical Resistivity of Metals and Alloys; Cambridge University Press: Cambridge, UK, 1991. [Google Scholar]
  43. Torokhtii, K.; Pompeo, N.; Silva, E.; Alimenti, A. Optimization of Q-factor and resonance frequency measurements in partially calibrated resonant systems. Meas. Sens. 2021, 18, 100314. [Google Scholar] [CrossRef]
  44. Torokhtii, K.; Alimenti, A.; Pompeo, N.; Leccese, F.; Orsini, F.; Scorza, A.; Sciuto, S.A.; Silva, E. Q-factor of microwave resonators: Calibrated vs. uncalibrated measurements. J. Phys. Conf. Ser. 2018, 1065, 052027. [Google Scholar] [CrossRef] [Green Version]
  45. Torokhtii, K.; Alimenti, A.; Pompeo, N.; Silva, E. Frequency span optimization for asymmetric resonance curve fitting. In Proceedings of the 2021 IEEE International Instrumentation and Measurement Technology Conference (I2MTC), Glasgow, UK, 17–20 May 2021; pp. 1–5. [Google Scholar] [CrossRef]
  46. Petersan, P.J.; Anlage, S.M. Measurement of resonant frequency and quality factor of microwave resonators: Comparison of methods. J. Appl. Phys. 1998, 84, 3392–3402. [Google Scholar] [CrossRef] [Green Version]
  47. Kajfez, D. Q Factor Measurements Using MATLAB; Artech House: Norwood, MA, USA, 2011. [Google Scholar]
  48. Alimenti, A.; Torokhtii, K.; Silva, E.; Pompeo, N. Challenging microwave resonant measurement techniques for conducting material characterization. Meas. Sci. Technol. 2019, 30, 065601. [Google Scholar] [CrossRef]
  49. Leong, K.T.M. Precise Measurements od Surface Resistance of High Temperature Superconducting Thin Films Using a Novel Method of Q-Factor Computations for Sapphire Dielectric Resonators in the Transmission Mode. Ph.D. Thesis, Electrical and Computer Engineering, School of Engineering, James Cook University, Townsville, Australia, 2000. [Google Scholar]
  50. Maxwell, E. Conductivity of metallic surfaces at microwave frequencies. J. Appl. Phys. 1947, 18, 629–638. [Google Scholar] [CrossRef] [Green Version]
  51. Tischer, F.J. Measurement of the surface resistance of single-crystal copper in the millimeter-wave region at room temperature. Appl. Phys. 1974, 5, 285–286. [Google Scholar] [CrossRef]
  52. Hakki, B.W.; Coleman, P.D. A Dielectric Resonator Method of Measuring Inductive Capacities in the Millimeter Range. IEEE Trans. Microw. Theory Tech. 1960, 8, 402–410. [Google Scholar] [CrossRef]
  53. Krupka, J.; Derzakowski, K.; Tobar, M.; Hartnett, J.; Geyer, R.G. Complex permittivity of some ultralow loss dielectric crystals at cryogenic temperatures. Meas. Sci. Technol. 1999, 10, 387. [Google Scholar] [CrossRef]
  54. Hammerstad, E.O. Microstrip Handbook; Norvegian Istitute of Technology-Trondheim University: Trondheim, Norway, 1975; Volume 74169. [Google Scholar]
  55. Dang, Z.M.; Yuan, J.K.; Zha, J.W.; Zhou, T.; Li, S.T.; Hu, G.H. Fundamentals, processes and applications of high-permittivity polymer–matrix composites. Prog. Mater. Sci. 2012, 57, 660–723. [Google Scholar] [CrossRef]
  56. Isakov, D.; Lei, Q.; Castles, F.; Stevens, C.; Grovenor, C.; Grant, P. 3D printed anisotropic dielectric composite with meta-material features. Mater. Des. 2016, 93, 423–430. [Google Scholar] [CrossRef]
  57. Hu, H.; Sinha, S.; Meisel, N.; Bilén, S.G. Permittivity of 3d-printed nylon substrates with different infill patterns and densities for design of microwave components. Designs 2020, 4, 39. [Google Scholar] [CrossRef]
  58. Arroyave, G.R.; Quijano, J.A. Broadband characterization of 3D printed samples with graded permittivity. In Proceedings of the 2018 International Conference on Electromagnetics in Advanced Applications (ICEAA), Cartagena, Colombia, 10–14 September 2018; pp. 584–588. [Google Scholar]
  59. Goulas, A.; Zhang, S.; McGhee, J.R.; Cadman, D.A.; Whittow, W.G.; Vardaxoglou, J.C.; Engstrøm, D.S. Fused filament fabrication of functionally graded polymer composites with variable relative permittivity for microwave devices. Mater. Des. 2020, 193, 108871. [Google Scholar] [CrossRef]
  60. Castles, F.; Isakov, D.; Lui, A.; Lei, Q.; Dancer, C.; Wang, Y.; Janurudin, J.; Speller, S.; Grovenor, C.; Grant, P.S. Microwave dielectric characterisation of 3D-printed BaTiO3/ABS polymer composites. Sci. Rep. 2016, 6, 22714. [Google Scholar] [CrossRef] [Green Version]
  61. Deffenbaugh, P.I.; Rumpf, R.C.; Church, K.H. Broadband microwave frequency characterization of 3-D printed materials. IEEE Trans. Comp. Pack. Man. 2013, 3, 2147–2155. [Google Scholar] [CrossRef]
  62. Alimenti, A.; Torokhtii, K.; Pompeo, N.; Piuzzi, E.; Silva, E. Characterisation of dielectric 3D-printing materials at microwave frequencies. Acta IMEKO 2020, 9, 26–32. [Google Scholar] [CrossRef]
  63. Baker-Jarvis, J.; Vanzura, E.J.; Kissick, W.A. Improved technique for determining complex permittivity with the transmission/reflection method. IEEE Trans. Microw. Theory Tech. 1990, 38, 1096–1103. [Google Scholar] [CrossRef] [Green Version]
  64. Bagguley, D.M.S.; Owen, J. Microwave properties of solids. Rep. Prog. Phys. 1957, 20, 304. [Google Scholar] [CrossRef]
  65. Krishna, K.G. A method of determining the dipole moment and relaxation time from microwave measurements. Trans. Faraday Soc. 1957, 53, 767–770. [Google Scholar] [CrossRef]
  66. Dakin, T.; Works, C. Microwave dielectric measurements. J. Appl. Phys. 1947, 18, 789–796. [Google Scholar] [CrossRef]
  67. Anlage, S.M. Microwave superconductivity. IEEE J. Microw. 2021, 1, 389–402. [Google Scholar] [CrossRef]
  68. Krupka, J. Microwave Measurements of Electromagnetic Properties of Materials. Materials 2021, 14, 5097. [Google Scholar] [CrossRef] [PubMed]
  69. Krkotić, P.; Tagdulang, N.; Calatroni, S.; O’Callaghan, J.M.; Pont, M. Potential impedance reduction by REBCO coated conductors as beam screen coating for the future circular hadron collider. Europhys. Lett. 2022, 140, 64001. [Google Scholar] [CrossRef]
  70. Hempel, H.; Savenjie, T.J.; Stolterfoht, M.; Neu, J.; Failla, M.; Paingad, V.C.; Kužel, P.; Heilweil, E.J.; Spies, J.A.; Schleuning, M.; et al. Predicting solar cell performance from terahertz and microwave spectroscopy. Adv. Energy Mater. 2022, 12, 2102776. [Google Scholar] [CrossRef]
  71. Alimenti, A.; Torokhtii, K.; Di Gioacchino, D.; Gatti, C.; Silva, E.; Pompeo, N. Impact of Superconductors’ Properties on the Measurement Sensitivity of Resonant-Based Axion Detectors. Instruments 2021, 6, 1. [Google Scholar] [CrossRef]
  72. Braine, T.; Rybka, G.; Baker, A.; Brodsky, J.; Carosi, G.; Du, N.; Woollett, N.; Knirck, S.; Jones, M. Multi-mode Analysis of Surface Losses in a Superconducting Microwave Resonator in High Magnetic Fields. arXiv 2022, arXiv:2208.11799. [Google Scholar]
  73. Abdolrazzaghi, M.; Nayyeri, V.; Martin, F. Techniques to improve the performance of planar microwave sensors: A review and recent developments. Sensors 2022, 22, 6946. [Google Scholar] [CrossRef]
  74. Dargent, T.; Haddadi, K.; Lasri, T.; Clément, N.; Ducatteau, D.; Legrand, B.; Tanbakuchi, H.; Theron, D. An interferometric scanning microwave microscope and calibration method for sub-fF microwave measurements. Rev. Sci. Instrum. 2013, 84, 123705. [Google Scholar] [CrossRef] [Green Version]
  75. Happy, H.; Haddadi, K.; Theron, D.; Lasri, T.; Dambrine, G. Measurement techniques for RF nanoelectronic devices: New equipment to overcome the problems of impedance and scale mismatch. IEEE Microw. Mag. 2014, 15, 30–39. [Google Scholar] [CrossRef]
  76. Kharkovsky, S.; Zoughi, R. Microwave and millimeter wave nondestructive testing and evaluation-Overview and recent advances. IEEE Instrum. Meas. Mag. 2007, 10, 26–38. [Google Scholar] [CrossRef]
  77. Teppati, V.; Ferrero, A.; Sayed, M. Modern RF and Microwave Measurement Techniques; Cambridge University Press: Cambridge, UK, 2013. [Google Scholar]
  78. Rüfenacht, A.; Flowers-Jacobs, N.E.; Benz, S.P. Impact of the latest generation of Josephson voltage standards in ac and dc electric metrology. Metrologia 2018, 55, S152. [Google Scholar] [CrossRef]
  79. Gugliandolo, G.; Tabandeh, S.; Rosso, L.; Smorgon, D.; Fernicola, V. Whispering gallery mode resonators for precision temperature metrology applications. Sensors 2021, 21, 2844. [Google Scholar] [CrossRef] [PubMed]
  80. Krupka, J.; Nguyen, D.; Mazierska, J. Microwave and RF methods of contactless mapping of the sheet resistance and the complex permittivity of conductive materials and semiconductors. Meas. Sci. Technol. 2011, 22, 085703. [Google Scholar] [CrossRef]
  81. International Electrotechnical Commission (IEC). International Standard 61788—15 Superconductivity—Part 15: Electronic 308 Characteristic Measurements—Intrinsic Surface Impedance of Superconductor Films at Microwave Frequencies; Technical Report; IEC: Geneva, Switzerland, 2011. [Google Scholar]
  82. Alimenti, A.; Torokhtii, K.; Pompeo, N.; Silva, E. High precision and contactless dielectric loaded resonator for room temperature surface resistance measurements at microwave frequencies. In Proceedings of the 2022 IEEE International Instrumentation and Measurement Technology Conference (I2MTC), Ottawa, ON, Canada, 16–19 May 2022; pp. 1–6. [Google Scholar]
  83. Torokhtii, K.; Alimenti, A.; Vidal García, P.; Pompeo, N.; Silva, E. Proposal: Apparatus for Sensing the Effect of Surface Roughness on the Surface Resistance of Metals. Sensors 2023, 23, 139. [Google Scholar] [CrossRef]
  84. Riddle, B.; Baker-Jarvis, J.; Krupka, J. Complex permittivity measurements of common plastics over variable temperatures. IEEE Trans. Microw. Theory Tech. 2003, 51, 727–733. [Google Scholar] [CrossRef]
Figure 1. Sketch of a cylindrical dielectric-loaded resonator (DR). Between the metallic bases, a low losses dielectric crystal is loaded to increase the quality factor Q, and thus the sensitivity, of the DR. When used as a fixture for the measurement of the surface resistance R s of conductors, a part of the surface S i of the metallic enclosure is substituted with the sample under study while, when used for the characterization of dielectric materials, the sample of volume V j is loaded in the DR.
Figure 1. Sketch of a cylindrical dielectric-loaded resonator (DR). Between the metallic bases, a low losses dielectric crystal is loaded to increase the quality factor Q, and thus the sensitivity, of the DR. When used as a fixture for the measurement of the surface resistance R s of conductors, a part of the surface S i of the metallic enclosure is substituted with the sample under study while, when used for the characterization of dielectric materials, the sample of volume V j is loaded in the DR.
Sensors 23 00518 g001
Figure 2. Three-dimensional exploded view of the DR design.
Figure 2. Three-dimensional exploded view of the DR design.
Sensors 23 00518 g002
Figure 3. The measurement sensitivity is maximized through the optimization of the geometry of the resonator. Once the dimensions of the dielectric crystal are fixed, the thickness h P T F E of the PTFE layer between the crystal and the sample holder, and that of the metallic mask h m a s k , can be optimized. The optimization is performed with e.m. simulations by exploring the effects on the sensitivity coefficient c R = Q 2 / G s when h P T F E and h m a s k are varied within mechanically reasonable limits. Curve 1 (blue circles): sensitivity c R optimization by varying h P T F E (lower horizontal axis). Curve 2 (orange squares): sensitivity c R optimization by varying h m a s k (upper horizontal axis) once fixed h P T F E = 0.75 mm, which corresponds to the maximum c value shown in curve 1.
Figure 3. The measurement sensitivity is maximized through the optimization of the geometry of the resonator. Once the dimensions of the dielectric crystal are fixed, the thickness h P T F E of the PTFE layer between the crystal and the sample holder, and that of the metallic mask h m a s k , can be optimized. The optimization is performed with e.m. simulations by exploring the effects on the sensitivity coefficient c R = Q 2 / G s when h P T F E and h m a s k are varied within mechanically reasonable limits. Curve 1 (blue circles): sensitivity c R optimization by varying h P T F E (lower horizontal axis). Curve 2 (orange squares): sensitivity c R optimization by varying h m a s k (upper horizontal axis) once fixed h P T F E = 0.75 mm, which corresponds to the maximum c value shown in curve 1.
Sensors 23 00518 g003
Figure 4. Vertical straight section of the designed DR. All the dimensions are in mm.
Figure 4. Vertical straight section of the designed DR. All the dimensions are in mm.
Sensors 23 00518 g004
Figure 5. Pictures of the realized DR. On the left, the open cavity with the sapphire crystal. On the right, the closed DR with the external supporting structure.
Figure 5. Pictures of the realized DR. On the left, the open cavity with the sapphire crystal. On the right, the closed DR with the external supporting structure.
Sensors 23 00518 g005
Figure 6. Measurement of the S12 parameter on the complex plane.
Figure 6. Measurement of the S12 parameter on the complex plane.
Sensors 23 00518 g006
Figure 7. Q (red circles) and f 0 (blue squares) measurement repeatability evaluation.
Figure 7. Q (red circles) and f 0 (blue squares) measurement repeatability evaluation.
Sensors 23 00518 g007
Figure 8. Real part of the relative complex permittivity (blue circles, left axis) and imaginary part (red triangles, right axis) measured on samples prepared with different porosity. The horizontal axis shows the filling percentage γ of the samples.
Figure 8. Real part of the relative complex permittivity (blue circles, left axis) and imaginary part (red triangles, right axis) measured on samples prepared with different porosity. The horizontal axis shows the filling percentage γ of the samples.
Sensors 23 00518 g008
Figure 9. Three-dimensional drawing of the printed sample at 68 % filling.
Figure 9. Three-dimensional drawing of the printed sample at 68 % filling.
Sensors 23 00518 g009
Figure 10. Calibration curves ε ( Δ f 0 ) (red circles) and η s ( Δ f 0 ) (blue squares). The points are obtained with e.m. simulations and fitted with a 2nd-order polynomial. The obtained calibration curves are reported in Equations (10) and (11).
Figure 10. Calibration curves ε ( Δ f 0 ) (red circles) and η s ( Δ f 0 ) (blue squares). The points are obtained with e.m. simulations and fitted with a 2nd-order polynomial. The obtained calibration curves are reported in Equations (10) and (11).
Sensors 23 00518 g010
Table 1. List of the measured sample. The nominal R is evaluated at ∼ 12.9 GHz from the dc resistivity ρ of the material, assuming R s = μ 0 π f ρ . The samples are square in shape and of dimensions 15 × 15 × 3 mm 3 .
Table 1. List of the measured sample. The nominal R is evaluated at ∼ 12.9 GHz from the dc resistivity ρ of the material, assuming R s = μ 0 π f ρ . The samples are square in shape and of dimensions 15 × 15 × 3 mm 3 .
Ref.MaterialNominal R (m Ω )
R0Brass 55 ÷ 68
R1Copper29
R2Aluminum38
R3Zinc55
Table 2. Δ R s measurement results.
Table 2. Δ R s measurement results.
Sample Δ R (m Ω ) u ( Δ R ) (m Ω )
R035.41.1
R1ref-
R211.01.1
R335.01.1
Table 3. Measurement results.
Table 3. Measurement results.
Sample R SM (m Ω ) u ( R ) (m Ω )
R09212
R15812
R26812
R39212
Table 4. Lattice characteristics of the samples.
Table 4. Lattice characteristics of the samples.
Lattice Type l p (mm) h (mm)Filling %
---100
square1.600.4081
hexagonal1.000.3068
hexagonal1.500.5650
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Alimenti, A.; Torokhtii, K.; Vidal García, P.; Pompeo, N.; Silva, E. Design and Test of a New Dielectric-Loaded Resonator for the Accurate Characterization of Conductive and Dielectric Materials. Sensors 2023, 23, 518. https://doi.org/10.3390/s23010518

AMA Style

Alimenti A, Torokhtii K, Vidal García P, Pompeo N, Silva E. Design and Test of a New Dielectric-Loaded Resonator for the Accurate Characterization of Conductive and Dielectric Materials. Sensors. 2023; 23(1):518. https://doi.org/10.3390/s23010518

Chicago/Turabian Style

Alimenti, Andrea, Kostiantyn Torokhtii, Pablo Vidal García, Nicola Pompeo, and Enrico Silva. 2023. "Design and Test of a New Dielectric-Loaded Resonator for the Accurate Characterization of Conductive and Dielectric Materials" Sensors 23, no. 1: 518. https://doi.org/10.3390/s23010518

APA Style

Alimenti, A., Torokhtii, K., Vidal García, P., Pompeo, N., & Silva, E. (2023). Design and Test of a New Dielectric-Loaded Resonator for the Accurate Characterization of Conductive and Dielectric Materials. Sensors, 23(1), 518. https://doi.org/10.3390/s23010518

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop