1. Introduction
The electric complex permittivity is one of the most important parameters of material characterization. It is utilized in a large range of applications such as: Material description [
1,
2], tests of organic tissue [
3,
4] microfluidics [
5,
6,
7,
8], bio sensing [
9,
10,
11], ecological operators [
12,
13], and quality control in the food industry [
14,
15]. Accurate determination of the permittivity is an important task and a great challenge for microwave engineering, in general, and therefore many solutions have been investigated lately.
A relatively new option for implementing sensitive planar sensors is to use metamaterials. Metamaterial structures present a major advantage over other conventional options: They can be artificially tailored to achieve better resolution and accuracy. In the last few years, an increased interest for studying the sensors based on resonant metamaterial structures such as Split Ring Resonator (SRR) and Complementary Split Ring Resonator (CSRR) has been noticed due to minimal efforts in sample preparation, nondestructive effect, ability to characterize both low and high losses materials, and higher sensitivity [
2,
16].
Considering these important advantages, we present a modified SRR sensor for accurate complex permittivity of solid dielectrics. It is implemented in planar technology and it exhibits two resonant frequencies, which are used to overcome numerical limitations that may appear in real-life measurements. The results obtained for our sensor are compared with similar approaches existing in the literature.
There are many SRRs sensors depending on a single resonant frequency, which is produced by the related resonator circuit, and the main focus of the authors is only on the detection of the real part of permittivity.
For example, in reference [
5], an SRR sensor working at 2 GHz is proposed. It is implemented in planar technology, on Rogers AD1000 substrate, and it is used to measure the thickness of thin films, as well as the electric permittivity for both dielectrics and liquids, being fully submersible. The extraction of the imaginary part of the permittivity is not rigorous as the authors conclude [
5]. They also argue that a more thorough and systematic study investigating many more combinations of real and imaginary parts of the permittivity [
5] should be considered for further improvements. Nevertheless, the authors suggest that another version of the sensor should be developed to obtain two resonance frequencies, in order to determine the complex permittivity and the thickness of solids at the same time [
5]. The sensor we propose has two resonant frequencies able to measure complex permittivity with errors smaller than 4.5% and for samples with thicknesses from 1 mm to 10 mm.
In ref. [
16], an interdigital capacitor based SRR (IDC-SRR) sensor for dielectric testing is investigated. The authors also propose and analyze a meandered line based split ring resonator (ML-SRR) RF sensor for magnetic testing. Both sensors work at 2.45 GHz and are implemented on RT/Duroid 6006 substrate. The accuracy of the real part of permittivity measurements is more than 94% [
16], but still, these sensors are not able to measure the complex permittivity as our proposed sensor. Our sensor can measure the imaginary part of permittivity with comparable accuracy.
Another SRR based sensor for magnetodielectric substrates characterization is the one presented in [
17]. The device is fabricated using the microstrip technology on a 1.27mm-thick RT/Duroid 6006 substrate and working at 2.5 GHz [
17]. The SRRs are magnetically coupled to the microstrip line allowing both electric permittivity and magnetic permeability measurements. Still, only the real parts are measured. Our sensor can measure both real and imaginary parts for the complex permittivity and the errors for real part measurements are comparable for the same samples as in [
17].
A different approach for the sensor design is presented in [
18]. It uses a two-layer and three-layer magnetic coupled SRR for higher sensitivity, better stability, and stronger anti-jamming ability from the external interface. These sensors have dimensions of 0.052λ × 0.052λ allowing miniaturization, but no imaginary parts of permittivity measurements are carried out [
18]. The influence of the thickness of the sample is not considered in this study. On the other hand, the errors for real part of permittivity measurements are similar to ours.
Another modified SRR sensor is presented in [
19], but it is used for thin-film detection, not for thicker dielectrics as ours. Furthermore, for this thin-film sensor, only the frequency shifts are investigated [
19], without determining the complex permittivity of the MUT as in our case.
In [
6], an SRR-based sensor is presented for measurement of complex permittivity of liquids. The sensitivity of the sensor is improved by overlapping the middle part of the outer ring of the SRR and part of the feeding line. It is best suited for measuring mixed liquids and determining the complex permittivity for each component, but no study regarding its application for solid materials characterization is done.
Regarding the characterization of material under test (MUT), several techniques have been proposed and employed for the permittivity. The most important ones can be categorized as free-space methods, transmission-line methods, and resonant cavity methods [
20].
The free-space method commonly employs using extremely directive lens and horn antennas placed on both sides of the MUT. The vector network analyzer (VNA) is connected to the antennas to measure the scattering parameters and phase constant to characterize the sample [
21]. This technique has the advantage of being contactless and not wasteful, but it demands the usage of expensive lenses and horn antennas, as well as the need for a large sample.
Another technique for measuring the material’s electric permittivity is the transmission-line method. In this method, the MUT is used as a loading material for transmission lines, such as a slice of material that can be incorporated to a waveguide [
22] or the deposing materials of a coaxial line that can be replaced by the MUT [
23]. The scattering parameters from the MUT-filled region provide the data necessary to extract the material’s properties. This method is comparatively lower cost than the free-space method. However, the sensitivity of the scattering parameters approach is not very efficient for low loss samples, and the sample elaboration is also very often a challenging task [
24]. The structures of microstrip-line and stripline are also used for this technique [
25].
One quite accurate technique is the resonant cavity method [
26]. In this method, a cavity resonator is loaded with the MUT, and the shift in the resonance frequency and the variation in the quality factor are determined. Circular resonators and microstrip-line resonators have been also employed for this purpose [
27] other than a conventional box resonator. This technique also needs accurate sample elaboration.
In order to overcome the limitations described above and to obtain all the information required to accurately characterize the complex electrical permittivity of a solid material, we propose a modified SRRs planar sensor for noninvasive complex permittivity measurements of solid materials. The resonant structures studied in this paper are based on the well-known SRR structures, with a modified topology, improving the selectivity and assuring two resonant frequencies. The measurement technique adopted in our research is the resonant cavity method.
The sensor is implemented on FR-4 substrate, offering portability, low-cost manufacturing, and easy-to-use and easy-to-interpret results. The structure is designed to be easily manufactured on a single metal layer, while allowing the easiness of integration of resonator elements at the same level.
The propagation phenomena occurring in these modified resonant structures can also be used to create RF absorbing materials, which can lead to designing efficient microwave absorbers for different applications, such as 5G antennas and automotives.
3. Numerical Analysis
In our further numerical analysis, we investigate both resonant frequencies as argued above. For each of them, we consider the data obtained after full wave electromagnetic simulation, and using a curve fitting tool, we determine analytical expressions for both the real and imaginary parts of the permittivity. These are expressed as a function of the resonant frequencies, the MUT’s thickness, and the quality factor of the loaded sensor. A curve fitting tool is often used in the literature [
1,
16,
17] for successfully estimating numerical expressions based on collected data from simulation or measurements.
In the unloaded situation, the simulated resonant frequencies (fr,1 and f’r,1) of the VS-SRR sensor are 1.24 GHz and 2.08 GHz.
Knowing that the quality factor for general resonators, Q can be written [
30]:
where
fr is the resonant frequency and
∆f represents the relative 3dB bandwidth of the resonator’s frequency response; we determine the quality factors corresponding to the two resonances as being equal to 35.4 and to 65, respectively.
After loading the sensor with the material under test, a shift in the resonant frequency as well as a change in the magnitude of S
21 (dB) are recorded as mentioned earlier in
Figure 5. The values for the resonant frequencies (
fr,1,
fʹr,1) and the quality factors are calculated from the response of the transmission coefficients and are then used to achieve a numerical expression with the aid of a curve fitting tool. A commercially available software OriginPro 8 [
16] is used as a curve fitting tool for the data obtained after full wave electromagnetic simulation in Ansys HFSS. The equations are generated using the sets of obtained data. A particular profile curve is chosen based on the least error between the chosen profile and the sets of numerically obtained data [
30], as it will be presented in the next sections.
3.1. Deduction of the Real Part of the Permittivity
To determine the type of dependency between the resonant frequency and the real part of the permittivity, we consider the expression for the resonant frequency [
31]:
The capacitance introduced by the load,
Cload depends directly proportional to the real part of the electrical permittivity [
31], so considering relation (8) as well, we obtain
f−2r ∝
εʹ
r.
So, for the proposed sensor, the affected transmission coefficient due to sample loading can be observed in
Figure 7. The inverse squares of the resonant frequencies (
fr,1 and
f ʹ
r,1) are extracted from the simulated transmission coefficient data and the results with the corresponding real permittivity (
εʹ
r) of the MUT are plotted and showed in
Figure 7a,b.
In
Figure 7a,b, it can be observed that the slope of the plotted curves relies on the thickness of MUT (
th) and increases as the thickness of the MUT increases. However, in
Figure 7a, it can be observed that the slope of the curve residues roughly constant for the sample thickness (
th) larger than 9 mm. This specific behavior may be noticed from the two lines corresponding to the MUT thickness of 9 mm and 12 mm, where both curves overlap.
Moreover, in
Figure 7a, it may be observed that all plotted curves corresponding to the MUT permittivity variation (
εʹ
r = 2 to
εʹ
r = 12) have a linear dependency, while in
Figure 7b, all plotted curves corresponding to the MUT permittivity (
εʹ
r > 10) have a roughly exponent dependency while the MUT thickness is increasing. In order to combine all the above effect, the dielectric constant of the specimen is mathematically expressed in terms of the family of straight lines as well as a family of exponential curves, where the freelance parameters are the resonant frequencies (
fr,1 and
fʹ
r,1 expressed directly in GHz) and the sample thickness (
th expressed directly in mm). By taking this aspect into consideration when using the fitting tool practically, the accuracy of the numerical model increases. So, based on the plotted data and using the curve fitting tool, we obtain the expressions for the real permittivity as a function of the MUT’s thickness and the two resonant frequencies:
These expressions will be used in our real-life measurements to determine which resonant frequency provides more accurate results. The mathematical limitation of relation (9) is for MUTs with thickness of 1 mm. In this case, relation (9) cannot be used, but we can use relation (10) instead. This case proves the limitation of the sensor, but also the importance of having an alternative analytical formula, based on the second resonant frequency.
3.2. Deduction of the Imaginary Part of the Permittivity
After finding the numerical relations (9) and (10) for determining the real part of the permittivity of the material under test, an identical analysis is completed to find a numerical relation for computing the loss tangent of the tested sample, which will give us information for determining the imaginary part of the permittivity.
The relation between the loss tangent, tan
δ, the quality of the proposed sensor after loading the MUT,
QMUT, the real part of the permittivity,
εʹ
r, and the imaginary part of the permittivity,
εʹʹ
r is given by [
32]:
where
QMUT states the quality factor of the proposed sensor after loading the MUT, which may be determined applying relation (7). The imaginary part of the complex permittivity is therefore obtained using (7) and (11).
At first, the real part of the permittivity takes values in the range of 2 to 12. For each value, the loss tangent is changed in the range from 0 to 0.12. For each change, the resonant frequencies (
fr,1 and
f ʹ
r,1) are recorded from the S
21 parameter’s simulation. Then, the quality factor is extracted from the simulation result of S
21 as depicted earlier in
Figure 5, for each resonance frequency. After that, the inverse of the quality factors for each corresponding loss tangent are plotted in
Figure 8a,b.
In
Figure 8a, it may be observed that all plotted curves corresponding to the MUT loss tangent variation (tan
δ = 0 to tan
δ = 0.12) have a semi-linear low slope component.
In
Figure 8b, all plotted curves remain on a semi-linear high slope component. Therefore, to deduce the tan
δ of the MUT, which relies on the loaded quality factor as well as the
εʹ
r of the MUT, a curve fitting tool is utilized. As in the previous case of the real part of permittivity, a commercially available software OriginPro 8 is used to determine the numerical model for both extracted results in
Figure 8a and b as presented in expressions (12) and (13), respectively:
After determining the εʹr from (9), (10) and tan δ from (12), (13) the imaginary part of the complex permittivity can be obtained using (11).
The mathematical limitation of relations (12) and (13) appear indirectly through the value of the real part of permittivity, εʹr. If this quantity cannot be determined using the first resonant frequency, as in the case of 1 mm thick MUTs, then automatically neither the loss tangent using relation (12) can be determined. As in the case of real part of permittivity, a second option, one using the second resonant frequency is very useful in practical applications.
4. Results
The sensor proposed in
Figure 4 is now implemented and measured. The substrate used is FR-4 (relative permittivity ε
r = 4.4 and the dissipation factor, tan δ, is approximately 0.02), with a thickness of 1.6 mm and cooper metallization electrodeposited on both sides of the substrate, with a thickness of 18 µm.
The technological development and manufacturing of the PCB sensor structure was made using Press-n-Peel Blue transfer foil from Techniks, with the etching process being done in turbulent and warm (approximately 50 °C) ferric chloride (FeCl
3), with the concentration of 38%. Press-n-Peel method uses a coated Mylar (Polyester) foil base material in which several layers of release agents and resist coatings are applied, as shown in
Figure 9a.
The width of the copper traces has been set to 1.52 mm, with the spacing between two concentric squares also being 1.52 mm, while the spacing between strips is 0.2 mm. The openings in the copper squares are equal to 1.22 mm. The ground plane size is 77 × 18 mm, being developed as a full copper plane onto the bottom side of the PCB. The overall dimensions of the sensor, access transmission lines included, is 76 mm × 26 mm, while the resonant structure itself, where the MUT is placed has a dimension of 44 mm × 18 mm.
The SMA (SubMiniature version A) connecters, which are classical semi-precision coaxial RF connectors used as interface for coaxial cables with screw-type coupling mechanism are mounted on the structure using mechanical welding. The SMA has a 50 Ω characteristic impedance and is designed to work in the range 0–18 GHz, fully matched with the necessities of the current sensing structure. The manufactured sensor is presented in
Figure 9b.
The measurement setup consists of the sensor connected to the Agilent E5071C (9 kHz to 6.5 GHz) network analyzer through 50 Ω cables, as shown in
Figure 9c. Before starting the measurements, a short-open-load-through (SOLT) calibration was performed using the Agilent calibration Kit. The number of sweep points is chosen 1601.
A set of materials under test: Transparent glass [
32], acrylic glass [
33], Teflon [
32], and Plexiglas [
34], with different thicknesses (
th) of 1 mm, 2 mm, 5 mm, and 10 mm is selected and used for measurements. For each measurement, the sensor is placed on a rough, stable surface and the MUT is carefully placed to cover the whole sensing area. Then, using the Agilent E5071C network analyzer, the magnitude of S
21 parameter is measured. Further, it is inspected and, using a marker, the resonant frequencies and the relative 3 dB bandwidth of the resonator’s frequency response are read. Then, using relation (7), the quality factor of the loaded sensor is determined. The quality factor for the loaded sensor with real MUTs is determined based only on measurements. Both resonant frequencies obtained after measurement for different types of MUT are considered and using relations (9) and (10), two possible values for the real part of the permittivity are obtained. They are compared with reference values [
32,
33,
34] and the results, including errors, are synthesized in
Table 2.
Analyzing the data obtained after measurements, it can be observed that both resonant frequencies can be used to compute the real part of the permittivity, except for the case of 1 mm thickness Plexiglas MUT. In this case, relation (9) cannot be used, because of numerical limitations, but relation (10) gives a value, measured with an error less than 3.5% than the reference value, proving the importance of an extra resonant frequency. The best results are obtained for samples with thicknesses of 2 mm and 5 mm. Another important observation is that when using the second resonant frequency, the errors are slightly larger than those corresponding to using the first resonant frequency. This can be explained because of technological imperfections and placing the probe in direct contact with the sensor. Still, taken into consideration that the errors are quite small for both frequencies, smaller than 4%, it can be concluded that the sensor is suitable for accurate real part of permittivity measurements. The small errors show that the gap between the sensor and the MUT can be ignored.
For the measurement of the imaginary part of permittivity, first, the measured quality factor Q
MUT is replaced in relations (12) and (13) and the value of the loss tangent, tan
δ is obtained. It is compared to the reference values [
32,
33,
34] and the results are given in
Table 3.
Analyzing the results in
Table 3, we notice that for acrylic glass of 5 mm thickness, the error when using the second resonant frequency is smaller than for the first one. Overall, the measurements were done with less than 4% errors. Again, because we could not determine the real part of the permittivity for the first resonant frequency, we could not determine the loss tangent either. A good solution for such cases is to use the alternative, given by the second frequency.
Next, using relation (11), the imaginary part of the permittivity is determined. The results of the measurements are given in
Table 4.
The errors in
Table 4 are smaller when using the second resonant frequency, except for the transparent glass case, proving the importance of the second resonant frequency. Additionally, the errors determined in
Table 4 show both the impact of approximations due to computing and the impact of measuring two parameters with different errors: The real part of permittivity and the loss tangent. So, we find cases when the errors are smaller than 1%, even if the corresponding errors for the loss tangent measurements alone are not that small. The observations regarding the technological imperfections and the placing procedure of the probe remain valid. Further, the impact of the air gap over the measurements was not considered and, still, the results are very good, much better than the ones in literature [
1]. For example, for Teflon, we have obtained measurement errors of 1.19% and 3.91% for the real permittivity and 3.88% and 1.231% for the imaginary part, proving the accuracy of the results. In reference [
1], the errors are 1.9% for the real part of permittivity and 8.6 2% for the imaginary part of permittivity. Nevertheless, it is worth observing that the thickness of the MUT has an impact on the overall response of the sensor. If the thickness of the MUT is increased, the interaction of the electromagnetic field is enhanced, so a change in the sensor’s frequency response is more obvious.
The results in the two tables show that the sensor can be used successfully to accurately characterize the dielectric parameters (dielectric constant and loss tangent) for both low-losses and lossy dielectrics, as well as for high dielectric constants dielectrics and small dielectric constants dielectrics.
5. Conclusions
In this paper, we present a modified SRRs planar sensor for noninvasive, accurate complex permittivity measurements of solid dielectrics. Starting from the classical SRR, a modified structure, using vertical strips added at a close distance of 0.2 mm to the SRR is investigated both from the enhanced selectivity perspective and from the overall dimensions. The result is a sensor made of two modified SRRs with lateral vertical strips, exhibiting high sensitivity for two resonant frequencies, at 1.24 GHz and 2.08 GHz.
A simplified equivalent circuit model is used to explain the microwave sensor’s design, and a very good agreement between the circuit model and the full electromagnetic simulation results is achieved. After a careful investigation, the two VS-SRRs sensor is selected to be further investigated. For each resonant frequency, we consider the data obtained after full wave electromagnetic simulation and using a curve fitting tool, we determine analytical expressions for both the real and imaginary parts of the permittivity. These are expressed as a function of the resonant frequencies, the MUT’s thickness, and the quality factor of the loaded sensor.
The sensor is implemented on an affordable, commercial substrate, FR-4 substrate, with a thickness of 1.6 mm, with reduced dimensions and being able to measure the real and imaginary parts of the permittivity for different solid dielectric samples, with errors less than 4.5% for both resonant frequencies in all analyzed cases. In our work, we have considered a large range of samples, with different thicknesses, different loss tangents, and dielectric constants to better investigate the sensor’s capabilities in real-life scenarios. The diversity of the samples helped us to observe the limitations of the numerical model developed in
Section 3 and find solutions to overcome them, such as successfully using the second resonant frequency.
Also, we have measured the quality factor both for the unloaded and loaded sensor using the resonant frequency and the relative 3dB bandwidth of the resonator’s frequency response. This approach added more practical consistency to our investigation. Still, some improvements can be done with respect to further miniaturization and the possibility to use this sensor for liquid dielectric characterization.
In future, we will investigate if this sensor can be used to measure the permeability for magnetic samples and if a lower losses substrate will improve the results. Last, but not least, the modified structure will be investigated if it is suitable for other resonant applications, which require the usage of similar configurations.