Preprint
Article

This version is not peer-reviewed.

Ultralinear Glucose Sensing Device with Membrane Supported Complementary Split-Ring Resonator on SIW

Submitted:

22 May 2026

Posted:

26 May 2026

You are already at the latest version

Abstract
The linearity of a glucose sensor is one of the most important indicators for ensuring measurement accuracy, system stability and post-measurement data processing efficiency. As yet, this issue has not been seriously investigated in the field of glucose sensing. In this work, the use of a substrate-integrated-waveguide (SIW) together with a specially tailored membrane-supported resonator (CSRR) is proposed for realizing an ultra-linear glucose sensor. Method: An SIW waveguide with two coplanar ports were realized on an FR4 substrate. On one side of the SIW waveguide, a CSRR was fabricated at the central region, which also served as the sensing region. On the other side of the SIW waveguide, the substrate at the central sensing region was thinned down to a negligible thickness using a nail grinder, thereby forming a circular recess with a membrane-thick bottom for holding a glucose-loaded test solution. The glucose concentration was determined by obtaining the resonant frequency shift in the reflection coefficient (i.e. S11). Results: The proposed glucose sensor has exhibited a highly linear correlation between the glucose concentration and the resonant frequency, with a Pearson correlation coefficient (r) reaching 0.99, even though the magnitudes of S11 minima were subjected to ambient electromagnetic interference. The sensitivity was found to be 1.6 MHz/(mg/dL). The resonant frequency was very weakly dependent on the volume of the test solution. Conclusion: Overall, the proposed glucose sensor has exhibited not only a high sensitivity but also a highly linear characteristics even in the presence of external electromagnetic interference.
Keywords: 
;  ;  ;  ;  
Subject: 
Engineering  -   Bioengineering

Introduction

Glucose sensing is currently one of the most popular topics of research. The public perception towards COVID has not only generated renewed research interests in the field of glucose sensing, but also a lot of demands in terms of the quality of a glucose sensing device. In the old days, the research community around the globe focused mainly on the accuracy of a glucose sensor, and any sensing device that satisfies the Clarke Error Grid criteria was usually considered acceptable. Nowadays, due to the practical requirements from the public, and in the wake of the global coronavirus outbreak in 2019, the quality requirements of a glucose sensor are much broader and stricter. Unlike before, not only should we pay attention to the in-vivo blood glucose measurement, but we should be mindful of the sugar concentration in the food or drink.
Of all the quality requirements under consideration by the research community, the most researched one is perhaps the sensitivity. In recent years, almost all the studies on experimental glucose sensing devices have devoted almost exclusively to sensitivity. What is relatively under-explored is the linearity of a glucose sensor. The linearity of a sensor is one of the important indicators for ensuring measurement accuracy, system stability, and data processing efficiency. Poor linearity in a sensor is an indication that more work will be required to get X if we know Y. With a highly linear glucose sensor, we can (with in reasonable limits) determine the glucose concentration from the measured samples by simple interpolation. If a glucose sensor has a poor linearity, on the other hand, we have to go to greater lengths to determine the glucose concentration when you know just the measured parameters. In general, the more linear the glucose sensor’s output, the easier it is to calibrate and to minimize uncertainty in its output scaling.
Lack of linearity frequently occurs in measurement processes. A poorly linear glucose sensor may simply have non-linear elements inside the sensor but this problem is difficult to correct. On the other hand, another more likely cause of non-linearity in a sensor is existence of measurement errors. In a contact-based glucose sensing process, one of the most common causes of errors that leads to non-linearity is volume discrepancy. Most of the experimental contact based glucose sensing device is usually designed to conduct measurements of glucose concentrations for samples of equal volume. When the volume of a sample is slightly different, the measured glucose concentration can be very different from the expected value.
Glucose sensing is undoubtedly a highly important field of science. Although this field has been established for more than a decade, there has not been any serious attempt to resolve the issues of volume discrepancies among different sample. Motivated by this concern, this study attempted to resolve the volume discrepancies among different sample solutions using the concept of spoof surface plasmon polaritons.
Reference Sensitivity (MHz/(mg/dL)) Lowest Detection Quantity Methodology Freq (GHz) Linearity
[1] 1.218 0.774 mg/dL SIW 4.65 Not known
[9] 0.42% 100 mg/dL SIW 3.43 Not known
[10] 0.310 20 mg/dL SIW 1.7 Not known
[11] ~2 at lowest concentration
~0.4 for concentration range (50 mg/dL - 100 mg/dL)
0 mg/dL CSIW 4.4 Not linear
[12] 0.25 for concentration range (106 mg/dL - 162 mg/dL)
0.235 for concentration range (162 mg/dL - 196 mg/dL)
0.257 for concentration range ( 196 mg/dL - 231 mg/dL)
106 mg/dL SIW 5.36 Not linear
This work 1.6 for all glucose concentrations 0 SIW with membrane-supported CSRR 15.4 Highly linear

Methodology

This section discusses the methodology of this work. To start with, the basic construction of the proposed sensor is first explained. An analysis that justifies the proposed design will be given. Finally, the measurement process focusing on the sensitivity and the linearity will be described in details.
A)
Basic Construction of the Proposed Glucose Sensor
The dimensions of the proposed glucose sensor is 26 mm x 8 mm x 1mm. Figure 1a and Figure 1b show the bottom side of the proposed glucose sensor facing upward. The proposed glucose sensor has been realized in FR4. The bottom of the sensing region was a complementary split ring resonator (CSRR) at the centre of a full-wave resonator realized in substrate-integrated waveguide (SIW). The CSRR was sitting on an ultra-thin FR4 membrane, of which the thickness was no more than 10 microns. The SIW full-wave resonator has a length of exactly 1 guided wavelength of 15.2 GHz. Connected to the input/output ports of the SIW full-wave resonator was coplanar waveguide ports. When used, port 1 is used as a sensing port to be connected to a vector network analyzer (VNA), whilst port 2 is to be terminated with a 50 ohm shunt resistor.
Figure 1c,d show the top side of the proposed glucose sensor. The top of the sensing region was a circular recess designed as a container to hold the test solution. In Figure 1d, the FR4 substrate is made transparent so that the vias are highly visible. During the fabrication process, the circular recess was ground using a finger nail grinder until the bottom of the sensing area was partially transparent (see Figure 1e for illustration). Under this condition, the depth of the circular recess was slightly less than the thickness of the FR4 substrate, and the bottom of the sensing area was literally an ultra-thin water-proof membrane which not only prevented water leakage but also protected the CSRR from any form of chemical attacks. According to our measurement, the thickness of this membrane was slightly less than 100 microns.
During the process of a measurement, the circular recess was filled with a glucose-water solution using a mini-pipette until it was fully filled. Port 2 will be terminated with a 50 ohm shunt resistor, whilst the reflection coefficient will be measured at port 1 using a VNA.
B)
Formation of Surface Waves along the FR4/CSRR Interface
At the interface, where z = 0 , the electric fields are continuous, i.e., E dz ( z = 0 ) = E mz ( z = 0 ) and E dy ( z = 0 ) = E my ( z = 0 ) . Using these boundary conditions, the following well known existence condition of spoof surface polaritons can be obtained:
k d k m = - ε d ε m
Substituting k d 2 = k y 2 - k o 2 ε d and k m 2 = k y 2 - k o 2 ε m into (2.1), we obtain the following dispersion relationship:
s k y k o = ε m ε d ε m + ε d
Equation (2.2) dictates that either ε m or ε d should be negative in order to make k y real and positive. ε d is positive and cannot be changed because it is the permittivity of the FR4 membrane suspending the CSRR. The only permittivity that is possibly negative is ε m , which is the permittivity of the CSRR. Metal is considered as a plasma with a permittivity given by:
ε m = 1 ω p 2 ω 2
where ω p is the plasma frequency of the metal. On the other hand, the CSRR is essentially a diluted plasma, of which the formula for the CSRR’s permittivity is of the same form but the plasma frequency ω p is reduced and obtained as a dual for the split-ring resonator (SRR). By substituting equation (2.3) into (2.2), we obtain:
k y k o = ( ω 2 - ω p 2 ) ε d ( ω 2 - ω p 2 ) + ω 2 ε d
For the purpose of glucose sensing, the magnitude of the transverse electric field is preferred to be maximized. The transverse electric field is given by:
E z = ± E 0 k y k z e x p j ( k y y + k z z ) ω t
In other words, k y is preferably maximum, whilst k z is preferably minimum. By setting k x in equation (2.1.16b), which means that ( ω 2 - ω p 2 ) + ω 2 ε d   should be set to zero, we have:
ω = ω p 1 + ε d
In equation (2.6), ω is the frequency limit where the velocity of the spoof surface plasmon polariton becomes minimum in the direction of propagation. This frequency limit is commonly known as surface plasma frequency. ω is also the resonant frequency where the combined permittivity along the interface between the CSRR and the FR4 membrane, ( ε d ε m ) / ( ε d + ε m ) , changes from very positive to very negative. The condition where ( ε d ε m ) / ( ε d + ε m ) undergoes such a positive-to-negative transition is commonly known as anomalous dispersion.
It is difficult to use equation (2.6) to obtain the surface plasma frequency without any knowledge of the CSRR’s plasma frequency. However, we can still obtain the resonant frequency of the CSRR using an analytical formula for the conventional split-ring resonator (SRR). According to Babinet’s principle of duality [3], the resonance frequencies of a split-ring resonator (SRR) and its complementary counterpart are identical if their corresponding geometries are the same. The major difference between SRR and its counterpart (i.e., CSRR) lies in the fact that SRR gives rise to a bandpass characteristic whilst CSRR gives rise to a bandstop characteristic [4].
Let’s apply the modified formula given in [6,7] to Figure 2a, in which the basic geometry of a dummy SRR is shown. The complementary counterpart of the SRR in Figure 2a is shown in Figure 2b. In Figure 2a, r , w   and g respectively stand for the inner-most radius of the circular metal track, the width of the metal track and the gap width. Another parameter that is not shown in Figs 2a and 2b is the metal thickness, h .
The resonant frequency is given by:
f o = 1 2 π L t ( C g + C s u r )
where L t is the total inductance of the metallic part given by:
L t = μ 0 ( r + 0.5 w ) l n r + 0.5 w h + w 0.5
C g is the capacitance due to the fringing electric in the gap, g, given as:
C g = ε 0 w h g + ( h + w + g )
C s u r is the capacitance due to the surface, which is given as:
C s u r = 2 ε 0 h π l n 4 r g
In this study, r , g and w were respectively 1mm, 0.2mm and 0.5mm. The metal thickness, h, is approximately 34 microns. With these data substituted to the formula in equation (2.7), By substituting the these data into equation (2,7), we obtain the resonant frequency to be 70.72 GHz. But this frequency ignores the dielectric substance attached to the CSRR, which is the FR4 membrane separating the CSRR and the water column.
The relative permittivity of water at around 15GHz is approximately 38. The thickness of the FR4 membrane (denoted as d f ) was 0.1 mm. The thickness of the water column (dentoted as d w ) was 1mm. Using Maxwell-Garnett approximation, the aggregate permittivity due to the combination of glucose and the FR4 membrane becomes:
ε d = d d + d w d f ε f + d w ε w 1 20.9
By applying substituting ε d into equation (2.6), the resonant frequency (or surface plasma frequency) becomes exactly 15.1 GHz. According to the Babinet’s principle of duality, this calculated resonant frequency is applicable to not only the SRR in Figure 3a but also to the CSRR in Figure 3b.
The calculated resonant frequency has been further verified by simulating the structure as shown in Figure 4a, where two CSRR’s are connected back to back on a microstrip line. The simulated performance is shown in Figure 4b, which has clearly shown that the simulated resonant frequency was 15.9 GHz, as opposed to the calculated value of 15.1 GHz.
The simulated S11 for the proposed glucose sensor, including the SIW together with a test solution sitting on the membrane-supported CSRR, is shown in Figure 5a. Consistent with our calculated value, Figure 5a suggests that the simulated resonant frequency is also 15.4 GHz when the sensing area of the proposed glucose sensor is fully filled with a test solution.
Figure 5b shows the 3-dimensional view of the simulated electric field distribution right at the interface between the CSRR and the FR4 membrane. Because the thickness of the FR4 membrane is much less than the three-dimensional guided wavelength, which should be around 0.9mm, the electric field in the vertical direction is largely evanescent in nature. Most of the electromagnetic power is trapped in the FR4 membrane in the form of a spoof surface plasmon polariton along the interface between the glucose-water and the FR4 membrane. Once the recess in the sensing region is fully filled with the test solution, the discrepancy in the volume of the test solution filling the sensing region is expected to have non-zero but relatively less influential effects on the resonant frequency.
C)
Sensitivity Analysis
The circular recess was filled with a glucose-water solution, of which the permittivity is much larger than that of FR4. In line with most of the dielectric waveguides, the electromagnetic energy should be preferentially stored in the recess or the sensing area. The liquid in the sensing region can be thought as a collection of different types of atoms, each of which is associated with its own resonance frequency (ωj), damping factor (γj) and the oscillator strength ( S j ). Here, j is just the index of the resonance corresponding to a particular type of atoms. The relative permittivity of the liquid can be readily modeled using Sellmeier Equation [8] as:
ε r = 1 + ω pm 2 j = 0 M S j / S 0 ω j 2 - ω 2 - j ω γ j
where ω pm is the plasma frequency. Since the liquid in the sensing region is directly attached to the FR4 membrane and indirectly attached to the complementary split-ring resonator, the plasma frequency should be the plasma frequency derived using equation (2.7). Accordingly, ω pm becomes:
ω pm =   2 π f o = 1 L t ( C g + C s u r )
As emphasized previously, the surface plasma frequency f o has already been determined to be around 15.1 GHz. At the surface plasma frequency, the electric field in the normal direction should be minimum, and the majority of the electromagnetic power is expected to be highly confined along the interface along the FR4 and the complementary split-ring resonator.
The surface plasma frequency ( ω pm or f o ) gives rise to a condition known as anomalous dispersion, under which the group velocity becomes greater than the speed of light, or even negative. As a result of this anomalous dispersive condition, there will be a perturbation in refractivity.
Suppose the operating frequency approaches a particular resonant frequency, ω j . We can assume that ω = ω j . Equation (3.1) can be reduced to:
ε r = 1 + ω pm 2 j S j / S 0 ω j γ j
Differentiating Equation (3.3) with respect to ω j , we get:
d ε r d ω j = ω pm ω j 2 - j S j / S 0 γ j
Using the well-known chain rule in calculus, equation (3.4) can be rewritten as:
d ε r d ω j = d ( ε + j ε ) dC dC d ω j = ω p , m ω j 2 - j S j / S 0 γ j
where C is the concentration of the glucose solution in the recess under a anomalous dispersion. The sensitivity can be thought as a change of resonant frequency per unit concentration, i.e., d ω j dC . The sensitivity should be a real number. It should not be imaginary. With a bit of algebraic re-arrangement, we can isolate the d ω j / d C from equation (3.5) in the following manner:
d ω j dC = S 0 S γ j ω j ω pm 2 d dC ( ε )
where γ j ,   S ,   S 0 ,   ω pm ,   C respectively stand for the damping factor, the oscillator strength at resonant frequency ( ω j ), the oscillator strength at the highest resonance, the plasma frequency and the glucose concentration.
Equation (3.6) suggests that the sensitivity of the glucose sensor is a function of the damping factor γ j , the surface plasma frequency ω pm , the resonant frequency ω j .
The sensitivity is directly proportional to the damping factor ( γ j ) caused by the conductivity of the glucose solution and the radiation resistance due to the antenna effects. Higher glucose concentration is in general indicative of a higher conductivity, which means a higher damping factor. The damping factor is inversely proportional to Q-factor. This Q-factor can be directly estimated by ratio of ( ω j ω ) over the frequencies in the neighborhood of the resonant frequency in the S11 plot.
In equation (3.6), one of most obvious factor influencing the sensitivity was obviously the ratio ω j / ω pm . For achieving a higher sensitivity, the plasma frequency, ω pm ,should be minimized, whilst the resonant frequency, ω j , should be maximized.
The term ω j / ω pm in equation (3.6) also suggests that not all the resonances are equal. The resonant frequency is always lower than or equal to the plasma frequency. Assuming that there exists more than one resonant frequency in the system, equation (3.6) suggests that the resonances at lower frequencies tend to yield a lower sensing sensitivity.
D)
Measurement Setup and Procedure
The performance of the proposed sensing device has been verified using the measurement setup as illustrated in Figure 6, where port 1 of the proposed sensor was connected to port 1 of a Rohde Schwarz vector network analyzer (ZVA 40), and port 2 of the same was terminated with a shunt resistor at 50 ohm.
Before conducting any measurement, a range of sample solutions of known glucose concentrations were prepared. During the measurement process, the sensing area (or the recess) of the proposed sensor was loaded with each of the samples using a pipette until the recess was fully filled (See Figure 6). The resulting S-parameters as displayed by ZVA 40 was then saved, whilst the resonant frequency was noted and recorded.
The concentrations of these samples include 0, 60mg/dL, 80mg/dL, 106mg/dL, 111mg/dL and 127mg/dL. The measurements have been replicated about 5 times but they were not done in one-go.

Measured Results and Discussion

The main finding of this work is depicted in Figure 7. The measured reflection coefficient (i.e., S11 parameter) for each of the sample glucose water has been plotted against frequency in Figure 7a, whilst Figure 7b shows a plot of resonant frequency as a function of glucose concentration. The resonant frequencies along the Y-axis of Figure 7b were originally extracted from Figure 7a.
At the first glance, the minima of the S11 in Figure 7a look highly erratic, with virtually no observable linearity in the correlation between the magnitude of S11 and the glucose concentration. The irregular magnitudes of S11 among different samples have been observed because the measurements were done at different times. On the other hand, the roles CSRR resonator are both a radiator and an absorbent of electromagnetic radiations. The ambient electromagnetic interference has changed from time to time. Under such a condition, the magnitudes of the minima of S11 were irregular and expected to be unstable.
After a closer look on the Figure 7b, however, it was found that the correlation between the resonant frequency and the glucose concentration was almost 100% linear, with a Pearson correlation coefficient (r) reaching 0.99. The high linearity as observed in Figure 7b suggests that the permittivity of the sensing region was relatively unaffected by the ambient electromagnetic interference.
The sensitivity was found to be 1.6 MHz/(mg/dL). Compared to counterparts as mentioned in Table 1, the sensitivity of the proposed was relatively high. The high sensitivity obtained in this work was attributed to the use of complementary split-ring resonator, which directly contributed a negative permittivity at a resonance.
Linearity is one of the most desired characteristics in sensing. With a highly linear sensing characteristic, a sensor usually exhibits a high measurement accuracy, high system stability and high post-measurement data processing efficiency. In this work, a high linearity was achieved even though the magnitudes of the reflective coefficient (S11) was erratic, unstable and subjected to the influences of unpredictable ambient electromagnetic interference.
The linearity of the proposed glucose sensor has been tested from 0 to a maximum glucose concentration of 127 mg/dL. Higher glucose concentrations have not been tested in this work. More research is needed to confirm if the linearity can be extended beyond 127 mg/dL.

Conclusion

In this work, a substrate-integrated-waveguide (SIW) has been used in conjunction with a specially tailored membrane-supported resonator to form an ultra-linear glucose sensor. The proposed glucose sensor has exhibited an ultra-high linearity in the correlation between resonant frequency and glucose concentration, with a Pearson correlation coefficient (r) reaching 0.99, even though the S11 minima have unstable magnitudes. The sensitivity was found to be 1.6 MHz/(mg/dL).

Acknowledgments

This research is funded by the Vietnamese-German University under grant number DTCS2025-008.

References

  1. Hamid Allah A, Ayissi Eyebe G, Domingue F. Improved Fully 3D-Printed SIW-Based Sensor for Non-Invasive Glucose Measurement. Sensors (Basel). 2025 Apr 9;25(8):2382. [CrossRef] [PubMed] [PubMed Central]
  2. Ghosh, Prabir & Nayak, Avinash & Gorai, Abhik. (2025). Slit Capacitor-Based SIW Sensor with Dual Sensitivity for Non-Invasive Glucose Detection. 1-5. [CrossRef]
  3. Julián D. Ortiz, Juan D. Baena, Ricardo Marqués, Amarachukwu N. Enemuo, Jonah Gollub, Roman, Akhmechet, Boyan Penkov, Chris Sarantos, and David T. Crouse, Babinet’s principle and saturation of the resonance frequency of scaled-down complementary metasurfaces, Appl. Phys. Lett. 118, 221901 (2021). [CrossRef]
  4. S. Moitra, S. Nayak, R. Regar, S. Kumari, K. Kumari, F. Parween and F. Naseem, “Circular Complementary Split Ring Resonators (CSRR) based SIW BPF”, 2019 Second International Conference on Advanced Computational and Communication Paradigms (ICACCP), 2019.
  5. Sydoruk, Oleksiy & Tatartschuk, E. & Shamonina, E. & Solymar, Laszlo. (2009). Analytical formulation for the resonant frequency of split rings. Journal of Applied Physics. 105. 014903-014903. [CrossRef]
  6. Cadence Design System, “Single Split-Ring Resonator Design”, Retrieved on 29 April, 2026, Retrieve from https://resources.system-analysis.cadence.com/blog/msa2021-single-split-ring-resonator-design.
  7. Smolyaninov, I., Balzano, Q. & Kozyrev, A. B. Surface electromagnetic waves at seawater-air and seawater-seafloor interfaces. IEEE Open. J. Antennas Propag. 4, 51–59 (2023). [CrossRef]
  8. Baird, C. Lecture 2 Notes, Electromagnetic Theory II (University of Massachusetts, 2024).
  9. H. M. Marzouk, A. S. A. El-Hameed, A. Allam, R. K. Pokharel and A. B. Abdel-Rahman, “Substrate Integrated Waveguide Sensor for Noninvasive Glucose Measurement,” 2023 11th International Japan-Africa Conference on Electronics, Communications, and Computations (JAC-ECC), Alexandria, Egypt, 2023, pp. 46-49. [CrossRef]
  10. Ameer B. Alsultani, Hussam Al-Saedi, Omer S. Alkhafaf, Katalin Kovács, J. Geoffrey Chase, Balazs Benyo, Development and experimental evaluation of single-port substrate integrated waveguide resonator with dual-parameter sensitivity for non-invasive blood glucose monitoring, Measurement, Volume 278, 2026, 121635, ISSN 0263-2241. [CrossRef]
  11. Mohd Bahar, A.A., Zakaria, Z., Md. Arshad, M.K. et al. Real Time Microwave Biochemical Sensor Based on Circular SIW Approach for Aqueous Dielectric Detection. Sci Rep 9, 5467 (2019). [CrossRef]
Figure 1. Construction of the proposed glucose sensor. (a) the bottom view of the glucose sensor; (b) the isometric view of the glucose sensor with the bottom metal facing upward; (c) the top view of the glucose sensor; (d) the isometric view of the glucose sensor with the top metal plate facing upward; (e) the photo showing the top metal plate of the glucose sensor; (f) the photo showing the bottom view of the glucose sensor.
Figure 1. Construction of the proposed glucose sensor. (a) the bottom view of the glucose sensor; (b) the isometric view of the glucose sensor with the bottom metal facing upward; (c) the top view of the glucose sensor; (d) the isometric view of the glucose sensor with the top metal plate facing upward; (e) the photo showing the top metal plate of the glucose sensor; (f) the photo showing the bottom view of the glucose sensor.
Preprints 214875 g001
Figure 3. Resonators: (a) split-ring resonator; (b) complementary resonator of the same geometry.
Figure 3. Resonators: (a) split-ring resonator; (b) complementary resonator of the same geometry.
Preprints 214875 g003
Figure 4. Simulation: (a) simulated circuit containing two CSRR’s connected back to back on a microstrip line; (b) a plot S-parameters as a function of frequency.
Figure 4. Simulation: (a) simulated circuit containing two CSRR’s connected back to back on a microstrip line; (b) a plot S-parameters as a function of frequency.
Preprints 214875 g004
Figure 5. Simulation results: (a) Simulated performance of the glucose sensor when the sensing area is loaded with a tap water. (b) 3-dimensional view of the simulated electric field distribution right at the interface between the CSRR and the FR4 membrane.
Figure 5. Simulation results: (a) Simulated performance of the glucose sensor when the sensing area is loaded with a tap water. (b) 3-dimensional view of the simulated electric field distribution right at the interface between the CSRR and the FR4 membrane.
Preprints 214875 g005
Figure 6. Measurement setup.
Figure 6. Measurement setup.
Preprints 214875 g006
Figure 7. Measured results: (a) S11 as a function of frequency for different glucose concentrations; (b) A plot of resonant frequency as a function of glucose concentration.
Figure 7. Measured results: (a) S11 as a function of frequency for different glucose concentrations; (b) A plot of resonant frequency as a function of glucose concentration.
Preprints 214875 g007aPreprints 214875 g007b
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.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

Accessibility

Disclaimer

Terms of Use

Privacy Policy

Privacy Settings

© 2026 MDPI (Basel, Switzerland) unless otherwise stated