# The quality monitoring of paracetamol medicament using a noninvasive microwave sensor

Youness Zaarour, Fatimazahrae El Arroud, Hafid Griguer, Rafiq El Alami, Mohammed El Kohen, Wiam Salhi, Abdessamad Faik, M'Hamed Drissi

► To cite this version:

Youness Zaarour, Fatimazahrae El Arroud, Hafid Griguer, Rafiq El Alami, Mohammed El Kohen, et al.. The quality monitoring of paracetamol medicament using a noninvasive microwave sensor. Scientific Reports, 2023, 13 (1), pp.17443. 10.1038/s41598-023-43409-y. hal-04282515

HAL Id: hal-04282515

https://hal.science/hal-04282515v1

Submitted on 13 Nov 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
<figure><img src="image_2.png" /></figure>

OPEN

# The quality monitoring of paracetamol medicament using a noninvasive microwave sensor

Youness Zaarour<fnref n="1" />✉, Fatimazahrae EL Arroud<fnref n="1" />, Hafid Griguer<fnref n="1" />, Rafiq El Alami<fnref n="1" />, Mohammed El Kohen<fnref n="1" />, Wiam Salhi<fnref n="1" />, Abdessamad Faik² & M'hamed Drissi<fnref n="1" />,³

Environmental conditions, including temperature, humidity, and light, can impact the quality of drugs. Microwave-based approaches offer a fast and cost-effective way to detect quality variations, providing an alternative to traditional techniques in the pharmaceutical and cosmetic industries. This article proposes the use of a microwave sensor for monitoring the quality of pharmaceutical drugs at distinct temperature levels. A small planar sensor based on three hexagonal split ring resonators (TH-SRR) is fabricated. The design is manufactured on an FR-4 dielectric substrate. The sensor is tested on a 1000 mg paracetamol tablet, at temperatures ranging from 40 to 80 °C. The Variation in the permittivity that characterizes product degradation is translated into a shift in the frequency of the scattering matrix elements. To validate the microwave approach, drug quality is examined with the laser-induced breakdown spectroscopy (LIBS) technique, an optical emission laser used for both qualitative and quantitative investigations of elements contained in a sample. The existing elements are classified using the National Institute of Standards and Technology (NIST) database and categorized according to their spectral line wavelengths. The experiments show the presence of optimal wavelength values for carbon (C), hydrogen (H), nitrogen (N), and oxygen (O) at 247.92 nm, 656.49 nm, 244.23 nm, and 777.48 nm, respectively. The microwave experimental results show a shift frequency of approximately 1 MHz on average when the tablet is heated at 80 °C for 15 min. Meanwhile, the LIBS measurement shows a remarkable shift in terms of intensity of approximately 8884 and 812 for carbon and hydrogen, respectively. Understanding how paracetamol dries under high temperatures and improving the process settings of the microwave sensor are investigated and assessed in this work.

Maintaining pharmaceutical quality during manufacturing has always been a priority for the World Health Organization (WHO)<fnref n="1" />. Every factor, whether individually or collectively, that affects product quality falls under the broad definition of quality control. Therefore, health outcomes can be plainly at risk if these products are not assessed to meet minimum criteria for quality, safety, efficacy, and significant medical needs. Good manufacturing practices (GMP) and other features, such as product design and development, are thus incorporated into quality assurance²,³. According to Diener et al., paracetamol is an analgesic used to relieve pain from a headache or other discomfort, but if it is not stored properly under predetermined conditions, its use could lead to kidney failure, liver damage, or the complete collapse of other essential organs⁴.

Temperature and humidity are the most crucial factors affecting the effectiveness of medicinal formulations. The aim of stability testing is to show the quality variation of a drug product over time under the effects of various environmental conditions⁵, which can accelerate oxidation, reduction, and hydrolysis reactions that lead to drug degradation. The findings consistently revealed that increasing the temperature from 25 to 45 °C across all storage conditions resulted in an extended disintegration time for the paracetamol tablets. These results underscore the significance of temperature as a crucial factor affecting the effectiveness and stability of medicinal formulations, emphasizing the importance of rigorous stability testing to ensure product quality and efficacy over time⁶. In addition to losing some of their effectiveness and decreasing their concentration, improperly stored drugs might even induce new or dangerous impacts. For this reason, Big Pharma (a collective term for major multinational pharmaceutical companies) is looking for effective solutions to control environmental conditions throughout the

<footnote n="1">¹Microwave Energy Sensing (MSE), DICE-University of Mohammed VI Polytechnic, 43152 Benguerir, Morocco. ²Laboratory for Inorganic Materials for Sustainable Energy Technologies (LIMSET), University of Mohammed VI Polytechnic, 43152 Benguerir, Morocco. ³Univ Rennes, INSA Rennes, IETR, UMR CNRS 6164, 35000 Rennes, France. ✉email: youness.zaarour@um6p.ma</footnote>
supply chain (production, storage, transport, etc.) in general but especially in tropical areas where humidity and temperature are high⁷⁻⁹. During processing, drugs are generally sampled at specific periods to measure their quality in relation to their moisture content. This method is time-consuming and only supplies an average result for the analyzed tablets¹⁰. Pharmaceutical quality analysis protocols are essentially based on offline and cumbersome laboratory methods. Although these methods are relatively effective and offer high-precision qualitative analysis, these so-called destructive methods remain unsuitable for online and portable analysis. In fact, Big Pharma is looking for online analysis solutions capable of providing drug quality data throughout the supply chain.

In the context of drugs and equivalent substances, many interesting technologies can be found in the literature for qualitative online analysis, such as optical and infrared spectroscopy¹¹, chromatographic analysis¹², and microwave sensing spectroscopy¹³. Microwave sensing technology has long been utilized as a fast and accurate method in research and industrial equipment, such as medical devices for glucose measurement¹⁴⁻¹⁶, gas sensing¹⁷,¹⁸, humidity sensing¹⁹, or breast tumors²⁰,²¹. This technique may offer an effective way to monitor drug quality in real-time and consequently reduce the time required for this task. Overall, this approach involves a direct interaction between an electromagnetic (EM) field created by a sensor and the material under test (MUT) for better extraction of the required information. It operates by monitoring resonance parameters (resonant frequency, amplitude, or quality factor) for permittivity variations caused by variations in the characteristics of the substance. Various studies have described this technique, with microwave sensors used for different applications. In²², a miniature nondestructive sensor using a microwave microstrip ring resonator was developed to assess the relative humidity of a single grain of wheat. This method was based on a microstrip line's electrical characteristics, which depend on the dielectric constant and thickness of the material that is being layered over it. Additionally, a microstrip-based sensor without a resonator that could be inserted directly into a container was suggested in²³ for grain quality evaluation and process control²⁴. Similarly, a test was carried out to examine the impact of microwave drying on the degradation of acetylsalicylic acid, a moisture-sensitive model medication²⁵. Additionally, the impact of various microwave powers (300, 600, and 900 W) on the degradation of the medication was examined. To the best of our knowledge, we did not identify in the open literature any other drugs with online qualitative analysis solutions based on a planar microwave sensing approach, especially for solid paracetamol. The present paper proposes a design, fabrication, and experimental study of a new low-cost miniature and planar microwave sensor deployed for solid paracetamol quality analysis under specific heating conditions. It mainly focuses on developing an online real-time portable device that can track drug quality by the calculated transmitted signal of the sensor while it is propagating through the drug. The functioning of the proposed microwave sensor, which operates at 3.48 GHz, is examined using a three-hexagonal split ring resonator (TH-SRR). The primary parameter studied here is the temperature condition, which directly affects the quality of drugs transported from Big Pharma to local shops. Temperatures can be increased by 22–27 °C in an hour and reach a maximum of 89 °C given the right circumstances, such as no ventilation and complete sun exposure. For this reason, we are estimating 80 °C as the highest temperature that can be achieved within a drug transport vehicle. In addition, it is necessary to avoid approaching the melting point of paracetamol, which is near 100 °C²⁶,²⁷ With the objective of confirming the heating impact on paracetamol drug quality, we chose a quantitative experimental analysis setup based on laser-induced breakdown spectroscopy (LIBS) that has been used in the pharmaceutical industry for qualitative drug analysis²⁸. This technique is capable of providing us with a spectrum that depicts several intensities. These intensities are situated at standardized wavelengths corresponding to the atomic level of each chemical element of the material under analysis.

In this work, the design rules and the main optimized structure of the sensor are presented. In addition, parametric simulation and the principal work are highlighted. Finally, the potential use of the sensor and the experimental results of both microwave and LIBS measurements are analyzed and assessed using the same heated sample of solid paracetamol.

# Material and methods

## Sample preparation and heating procedure

Paracetamol, also known as acetaminophen, is a chemical compound found in many drugs, including the most commonly prescribed by doctors. Paracetamol has long been the most widely used medication for relieving pain due to its excellent safety and tolerability record. Surprisingly, subsequent investigations have raised concerns about its therapeutic efficacy, especially under extreme environmental storage conditions. It is safe and effective when used and stored in an appropriate place but could possibly suffer from a loss of potency or even a reduction in efficacy, which could be harmful when stored at high temperatures. In our experiment, we employed solid effervescent paracetamol tablets (Doliprane)²⁹ that were purchased from a local pharmacy in Benguerir (Morocco) and kept under standard storage conditions. These 1000 mg tablets were obtained from the same production batch. To avoid any contamination, sample preparation was carried out in a controlled chamber with a set temperature (27 °C) and 30% humidity, and sample manipulation was performed while wearing protective gloves. For the analyses of temperature impacts on drug quality, several samples of heated paracetamol tablets were prepared. These drug samples were uniformly heated at 40 °C, 60 °C, and 80 °C with varying lengths of time from 15 to 60 min with a step of 15 min. As the objective was to simulate the temperature variations that might occur during the transport of the drugs, the chosen temperatures did not exceed 80 °C. The samples underwent the heating procedure described in Table 1:

* First, at a fixed temperature (80 °C), a set of four samples S1, S2, S3, and S4 was prepared to analyze the impact of the heating period on the paracetamol tablets.

* Furthermore, at a fixed period, three additional samples (S5, S6, S7) were heated at 80 °C to confirm the accuracy, readability, and reproducibility of the paracetamol drug behavior toward the sensor's response.
<table><caption>**Table 1.** The 10 prepared samples for studying the effect of the heating period and various temperatures.</caption><thead><tr><td>Samples</td><td>Period (min)</td><td>Temperature (°C)</td></tr></thead><tbody><tr><td>S1</td><td>15</td><td>80</td></tr><tr><td>S2</td><td>30</td><td>80</td></tr><tr><td>S3</td><td>45</td><td>80</td></tr><tr><td>S4</td><td>60</td><td>80</td></tr><tr><td>S5</td><td>15</td><td>80</td></tr><tr><td>S6</td><td>15</td><td>80</td></tr><tr><td>S7</td><td>15</td><td>80</td></tr><tr><td>S8</td><td>15</td><td>80</td></tr><tr><td>S9</td><td>15</td><td>40</td></tr><tr><td>S10</td><td>15</td><td>60</td></tr></tbody></table>

* Additionally, sample S8 was heated to 80 °C to track the effect of the temperature on the paracetamol tablet.

* Finally, two further samples S9 and S10 were heated at 40 °C and 60 °C, respectively, at a fixed heating period.

## Microwave setup

A two-port vector network analyzer (PICO VNA 108) was employed for microwave analysis. The VNA provides 0.03 MHz of frequency resolution and 0.01 dB of magnitude resolution. A SOLT (Short, Open, Load, Through) kit was used to calibrate the VNA in a frequency range from 2 to 4 GHz with a step of 0.5 MHz. This frequency band was chosen to provide a simple method for dielectric dispersion analysis. Our choice of operating frequency bands was guided by key considerations, particularly the extensive study of dielectric properties of organic materials over 2 GHz. This range aligns with molecular processes in compounds like Paracetamol, where observable dielectric changes are induced by molecular interactions and the realignment of polar functional groups. Paracetamol's possession of polar functional groups, such as amide and hydroxyl groups, makes it plausible to anticipate a dielectric response in this frequency band³⁰,³¹. Our study aims to connect these dielectric property changes to molecular structure and interactions, offering valuable insights into its behavior and applications. Additionally, the chosen frequency range of 2–4 GHz aligns practically with cost-effective commercial microwave components, making sensor implementation more accessible and reinforcing our research goals. To avoid any Fairfield coupling or external noises within the electromagnetic environment surrounding the sensor, the VNA microwave power was set at – 10 dBm to cover the selected frequency band. The VNA was able to provide real and imaginary parts of the Sij (i,j: 1,2) parameter matrix. In our study, we chose the S21 parameter, which represents the magnitude in dB of the transmission coefficient between the two VNA ports. For microwave transductance, we chose a two-port planar microwave sensor based on a TH-SRR (three hexagonal split ring resonator). The detailed design and results of the TH-SRR will be presented in the next sections. Our MUT was a 1000 mg paracetamol tablet. To ensure the proper location of the tablet above the sensor, a holder made of plexiglass with a relative permittivity of εr = 3.4 was installed in this area for MUT sensing. The holder had a suitable receptacle to secure the paracetamol tablet, which was in direct contact with the upper part of the sensor above the sensing area. The diameter of the holder receptacle was 23 mm to match the tablet diameter. The two ports of the sensor were linked to the VNA with a coaxial cable matched with 50Ω impedance. Real-time processing software was launched on a laptop to receive data from the VNA via serial communication (Fig. 1).

The following parameter is defined to qualitatively analyze the performance of the proposed microwave sensor in different heating procedures. The microwave quality shifting $Q_{μW}$ is defined as:

$$Q_{μW} = f_{ref} - f_{heated} \quad (MHz) \qquad (1)$$

Where $f_{ref}$ is the resonance frequency of the sensor loaded with a paracetamol drug at room temperature. This frequency is considered as the reference. $f_{heated}$ is the resonance frequency of the sensor for a drug heated at a specific value of temperatures (Fig. 2).

## Laser induced breakdown spectroscopy (LIBS) technique

To validate the microwave approach used for tracking the drug microwave quality shift introduced in the previous section, the LIBS technique was used to provide a reference laboratory analysis. LIBS is a sensitive optical method that can perform fast multielement analysis. It is based on establishing a qualitative and quantitative examination of the substance under test. This method of elemental analysis for pharmaceutical drug production can ultimately improve human health. In the past few years, LIBS has attracted interest from the pharmaceutical industry. This method has been utilized only in the examination of solid dosage forms for the purposes of quantifying active pharmacological ingredients in pharmaceuticals up to this point³². In the present research, we used the LIBS setup to identify the elements present in pharmaceutical 1000 mg paracetamol tablets. A pulsed laser was used as an excitation source in this investigation, with a wavelength of 1,064 nm and an energy expenditure of 25–50 mJ per pulse. A spectrometer is used with a wide spectral range, high sensitivity, and fast response rate
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Setup used for microwave measurement.</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>**Figure 2.** Definition of microwave quality shifting $Q_uW$.</figcaption></figure>

with a time-gated detector. It was connected to a computer for rapid processing and interpretation of the acquired data. LIBS was employed to examine and provide an analytic spectrum with quantitative intensities for the main chemical elements forming paracetamol (C, H, N, and O) at different temperature levels. These elements were extensively analyzed to ensure that they could be identified with high accuracy and little uncertainty. The data were gathered using a laser source and spectrometer to examine the spectral lines associated with different elements according to their wavelength and energy.

## Sensor design

The proposed symmetrical sensor is designed to operate at approximately 3.48 GHz. The establishment of the sensor geometry followed a systematic approach to attain an optimal configuration. Initially, we opted for a circular patch identical in size to the Paracetamol tablet. This choice not only facilitated direct contact between the sensor and the drug, but also ensured surface compatibility, enabling accurate dielectric property sensing. In conjunction with the circular patch, our design incorporated three hexagonal split ring resonators (TH-SRR) excited by a double transmission line (Fig. 3a), enhancing the sensor's overall performance. The sensor is fabricated in an FR-4 dielectric substrate with $ε_r = 4.4$ and tan$δ = 0.02$, as shown in Fig. 3b. This array of hexagonal resonators is adopted to create a concentrated electromagnetic near field, which is the result of the coupling between the adjacent hexagonal edges. This EM near field is focused and oriented perpendicular to the resonator plane, thus penetrating the tablets. To achieve impedance matching, the transmission line is constructed with an 18 mm strip length and a 3 mm strip width coupled to a 50Ω coaxial SMA connection along both sides. This hexagonal resonator array remains among the best shapes to acquire a strong coupling between individual resonator edges and produce a high-sensitivity behavior to detect slight variations in the microwave properties of MUTs. A simulation was performed with a high-performance HFSS 20.1 full-wave electromagnetic EM field simulator.
<figure><img src="image_5.png" /><figcaption>**Figure 3.** (a) Top view and (b) Cross-sectional view of the sensor. (Length L = 60 mm, width W = 60 mm, thickness h = 1.6 mm, mutual coupling distance a = 9 mm, resonator's side s = 4.6 mm, diagonal length d = , gap g = and path width p = 0.4 mm).</figcaption></figure>

The arrangement of three resonators offers various possibilities to control the sensor's response. Two of the critical factors that significantly influence the sensor's response are the gap of the resonator and the distance separating them, also known as mutual coupling distance. These factors significantly shaped response and overall sensitivity. Through meticulous optimization, we fine-tuned these parameters, leading to the final sensor geometry. In our case, as shown in Fig. 3a, the chosen array of three resonators is crafted to form a maximum interaction between six resonator edges. All resonators are separated by a mutual coupling distance of a = 9 mm between their geometric centers. Each resonator has a side of s = 4.6 mm with a gap of g = 0.6 mm and a path width of p = 0.4 mm. The diagonal length is d = 8 mm, and the resonators are coupled to each other with c = 1 mm. As shown in Fig. 4, in comparison to a single hexagonal resonator, the TH-SRR design exhibits a greater level of electric field confinement in the 20 × 20 mm² sensing area, with an intensity reaching 10⁴ V/m around the dielectric gaps, as well as the pathways connecting them. This high intensity can be explained by the significant interaction between the highly localized fields resulting from the greater electric field in the coupling region of the TH-SRR. The small width of the slots makes the resonators extremely sensitive to variations in the characteristics of the drug under examination and improves global sensitivity.

<figure><img src="image_6.png" /><figcaption>**Figure 4.** Full-wave simulation of the electric field distribution $\vec{E}$ in the upper surface plane for individual and array of resonators setups at 3.48 GHz with and without drug samples.</figcaption></figure>
## Results and discussion

To describe the intrinsic behavior of the unloaded sensor in depth, a parametric study using a geometric tuning process was carried out on a variety of parameters. It is known that the geometrical shape and resonator parameters define the resonance frequency of the proposed sensor when facing a specific material under test³³,³⁴. Additionally, determining the sensor resonance behavior depends critically on the design characteristics of each resonator and the coupling topology, such as the split gaps and the coupling distance between the resonators. The objective of this parametric study through simulation analysis is to obtain the optimal geometrical configuration that exhibits a higher sensitivity defined by a narrow frequency bandwidth and attenuated S₂₁ response. The bandwidth is defined here by the difference in the cutting frequencies at – 3 dB of the S₂₁ magnitude peak. In general, the reached array of the split ring resonators can be modeled as a shunt LC³⁵, where L is the equivalent inductance formed by the mutual inductances obtained by every two adjacent resonators and C is the equivalent capacitance formed by all the split gaps of the resonator's array.

In our case, as shown in Fig. 5, the full-wave simulation shows that the unloaded sensor exhibits band gap resonance behavior at approximately 3 GHz. Furthermore, the resonance peak moves toward higher frequencies while the split gap g is changed from 0.4 to 1 mm with a step of 0.2 mm. This is because the equivalent capacitance decreases as the gap size increases, hence increasing the dependency of the transmission response (Fig. 5a). The optimal sensitivity is obtained with the g = 0.8 mm configuration, which provides a narrower bandwidth at a 3.48 GHz resonance frequency. The impact of mutual coupling distance c is examined in Fig. 5b, revealing heightened sensitivity at c = 1 mm. Notably, it is important to emphasize that beyond this point, specifically after c = 2 mm, the resonance frequency shows no substantial influence on its magnitude.

After a thorough geometrical study, we proceeded to simulate and analyze the sensor's response to changes in the dielectric properties of the drug. As part of this analysis, we modeled the drug as a material under test, varying its permittivity within the range of 2 to 4. By examining the sensor's sensitivity across different permittivity values, we gained valuable insights into its performance characteristics and its ability to detect and measure variations in the dielectric properties of the drug. Figure 6 illustrates the sensitivity of the TH-SSR sensor based on its resonance frequency for different permittivity values ranging from 2 to 4 with a step of 0.1. The sensitivity (S) is expressed as the change in resonance frequency divided by the change in permittivity $S=\Delta f_r/\Delta\varepsilon_r$. The results clearly demonstrate a non-linear relationship between the sensor's sensitivity and the permittivity value. The sensitivity is highest around a permittivity value of 2, reaching approximately 150 MHz. However, as the permittivity value increases, the sensitivity gradually decreases, reaching around 98 MHz for higher permittivity values.

Following the simulations and analysis, we proceeded to experimentally evaluate the fabricated sensor using the test configuration depicted in Fig. 7a. As illustrated in Fig. 7b, the measured transmission coefficients S₂₁, as a function of frequency (dashed line) for the unloaded sensor and the numerical findings (continuous line), are in good agreement. The unloaded sensor's resonance is approximately 3.48 GHz, and its resonance peak is approximately – 27.6 dB within a BW of 2 MHz. This optimized sensor design will be used in future investigations since it provides the highest sensitivity performance.

The initial test is performed under a loaded sensor with a 1000 mg paracetamol tablet at room temperature (27 °C). The tablet and the active zone of the sensor are not separated by any air gap. Figure 8 presents the experimental results of the loaded sensor. The measured transmission coefficient response is – 26.22 dB at a resonance frequency of $f_{ref} = 3.0355$ GHz. The test clearly shows that there is a sensitivity to the presence of the paracetamol tablet with a frequency shift of approximately 444 MHz compared to the unloaded sensor response. This frequency shift is naturally due to the interaction of the radiated near field within the MUT, which is considered an additive dielectric substrate to the sensor. This obtained reference frequency $f_{ref}$ will be used to evaluate any frequency shifting when the loaded tablets are heated to different temperature levels.

<figure><img src="image_7.png" /><figcaption>**Figure 5.** Full-wave simulation of unloaded TH-SSR sensor for (a) split gap variation from g = 0.4 mm to 0.8 mm, with c = 1 mm; (b) coupling distance variation from c = 1 mm to 2 mm, with g = 0.8 mm.</figcaption></figure>
<figure><img src="image_8.png" /><figcaption>**Figure 6.** Sensitivity of the TH-SSR based sensor versus the different value of the real part of permittivity $ε'$.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Figure 7.** (a) Experimental setup, (b) simulated (solid line) and measured (dashed line) result of the unloaded sensor.</figcaption></figure>

In what follows, we study the correlation between different heating levels of the paracetamol tablets and their impact on the measured responses in terms of frequency shifting of the microwave sensor. After heating the S8 paracetamol sample at 80 °C over 15 min, its temperature decreased back to room temperature (27 °C). Every state's resonance frequency as the temperature drops are measured versus its particular temperature value. The surface temperature of the drug was measured by a Medisana TM A77 thermometer. Figure 9a provides the resonance frequency response of the microwave sensor at different temperatures (80, 67, 54, 46, 38, 34, and 27 °C, which is room temperature). Higher temperature values cause a decrease in the resonance frequency.
<figure><img src="image_10.png" /><figcaption>**Figure 8.** Measured result (solid line) of the loaded sensor with holder. The dashed line presents the measured results of the unloaded sensor with the holder.</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Figure 9.** (a) Resonance frequency of the sensor versus magnitude of $S_{21}$, as temperature is decreased from 80 to 27 °C. (b) Comparison between resonance frequency versus magnitude of $S_{21}$ before and after heating at 27 °C.</figcaption></figure>

Table 2 shows the resonance frequencies for each distinct temperature level. Figure 9b shows that after the tablet heating operation and stabilization of temperature at 27 °C, which is considered the reference temperature, we note that the microwave response obtained by the sensor does not correspond to that taken before the heating operation at the same reference temperature. In fact, the resonance frequency drops to 3.0345 GHz, and then the corresponding microwave quality shifting is approximately 1 MHz. A positive $Q_{μW}$ value (in this case 1 MHz) distinctly serves as an indicative marker of potential temperature-driven variations within the pharmaceutical

<table><caption>**Table 2.** Heating temperature versus resonance frequency (after heating).</caption><thead><tr><td>Temperature °C</td><td>Frequency of resonance = f<sub>heated</sub> (GHz)</td><td>Q<sub>μW</sub>=f<sub>ref</sub>-f<sub>heated</sub> (MHz)</td></tr></thead><tbody><tr><td>80</td><td>3.0212</td><td>14.30</td></tr><tr><td>67</td><td>3.0235</td><td>12.00</td></tr><tr><td>54</td><td>3.0249</td><td>10.60</td></tr><tr><td>46</td><td>3.0260</td><td>9.50</td></tr><tr><td>38</td><td>3.0293</td><td>3.30</td></tr><tr><td>34</td><td>3.0322</td><td>3.30</td></tr><tr><td>27</td><td>3.0345</td><td>1.00</td></tr></tbody></table>
<table><caption>**Table 3.** RMSE value for each proposed model.</caption><thead><tr><td>Model</td><td>Linear regression</td><td>Ridge regression</td><td>Random forest</td><td>Huber regression</td></tr></thead><tbody><tr><td>RMSE</td><td>2.1085</td><td>3.0529</td><td>1.9223</td><td>3.0585</td></tr></tbody></table>

compound, consequently impacting its overall quality. This observed variation arises from the dynamic alteration of the compound's dielectric properties in response to temperature changes, which in turn influences the dependence of its resonance frequency.

This phenomenon can be explained by the known heat effect caused by evaporating the tablet's internal humid solvent. Higher temperatures could cause the creation of a differential pressure, which encourages the quick evacuation of any humid substance related to the internal water amount inside the tablet. In our context, this evaporation could decrease the effective dielectric permittivity of the heated tablet and consequently lower the resonance frequency of the sensor.

A set of four 1000 mg drug samples S1, S2, S3, and S4 were prepared and heated at 15 min, 30 min, 45 min, and 60 min, respectively, at a fixed temperature of 80 °C. As shown in Fig. 10a, the $Q_{μW}$ becomes saturated after a 15 min heating period. As explained before, the heating procedure can cause evaporation of the internal moisture of the drug sample, but beyond 15 min, no additional amount of moisture evaporates. Thus, 15 min is the heating period used in future studies.

To evaluate the reproducibility of the sensor behavior, three identical 1000 mg paracetamol tablets were prepared and analyzed under the following procedure. S5, S6, and S7 were heated to 80 °C for a 15 min period. Directly after reaching this heating period, the heated samples were carefully placed on the same sensing zone of the sensor. The $Q_{μW}$ was continuously measured during the tablet's natural cooling, which itself is dependent on a room temperature of approximately 27 °C. At each $Q_{μW}$ measured value, the corresponding temperature was recorded. Figure 10b shows the microwave quality shifting $Q_{μW}$ recorded from the three iterations using S5, S6, and S7 versus different heating temperatures. Each tablet's temperature is lowered during the cooling process from 80 °C to the room temperature of 27 °C. The fitting curve (2) presents a linear regression model that better describes the variation of the temperature as a function of microwave quality shifting $Q_{μW}$ for each used sample. The exposure of drugs to high temperatures while being stored or transported could reduce their efficacy, and most licenses specify storage at 27 °C or less. Knowing the instantaneous recorded value of $Q_{μW}$ versus temperature would allow us to identify all of the temperatures to which the drug was subjected while being transported.

In practice, it is critical to optimize the parameters of a predictive model using a variety of performance metrics. The root-mean-squared error (RMSE) is then determined from the linear regression model and the collected data. Additionally, the RMSE calculation was employed to examine the effect of changing the drug position effect in the measurement. 2.1 °C of RMSE appeared and it could be due to the cooling ratio that can differ from each drug as well as the microwave measurement errors. However, the three iterations demonstrate that an identical $Q_{μW}$ value of roughly 14.3 MHz was obtained at the maximum heating temperature of 80 °C.

$$T = 26.108 * Q_{μW}^2 + 1.438 * Q_{μW} + 0.154 \quad (°C) \quad (2)$$

<figure><img src="image_12.png" /><figcaption>**Figure 10.** (a) Microwave quality shifting $Q_{μW}$ versus heating period at 80 °C, (b) Average microwave quality shifting $Q_{μW}$ versus temperature for 3 different paracetamol tablets.</figcaption></figure>
$$RMSE = \sqrt{\frac{1}{4} * \sum_{i=1}^{k} (T_i - T)} \qquad (3)$$

Along with linear regression, several other models of machine learning were used in this study to calculate and compare various RMSE values. Ridge regression³⁶ is a linear regression model with a quadratic constraint on the coefficients. A random forest regressor³⁷ is an algorithm built from a collection of decision trees. It is simple to grasp, quick to learn, and yields outcomes that are transferable. A Huber regression³⁸ is a type of linear regression used in robust regression that is less sensitive to data outliers than squared error loss. Table 3 provides an overview of the RMSEs of each method, demonstrating that the random forest regression and linear regression perform best by approximately 2 °C of RMSE and are the appropriate models to correct such errors.

## Heating temperature

In this experiment, we compared different heating temperatures (80 °C, 60 °C, and 40 °C) for 15 min, as shown in Fig. 11. Heating a drug at a high temperature leads to a higher shift frequency. At 80 °C, 14.3 MHz of microwave quality shifting $Q_{\mu W}$ is achieved. Therefore, after having cooled the sample, the resonance frequency of the sensor arrives at 1 MHz shifted to the reference. In contrast, 60 °C and 40 °C reach a shift frequency compared to the reference of 300 kHz and 50 kHz, respectively.

As a result, the change in water amount brought on by temperature vaporization shifts the resonance frequency of the loaded sensor, thus changing the dielectric characteristics of the sample. As a consequence, the rate of paracetamol variation (concentration) increased as the temperature increased. Moreover, for each of the three-sample situations, the suggested microwave sensor readings are successful in discriminating between different temperature levels.

To confirm the results obtained by the developed sensor, the LIBS technique was employed to examine the presence of C, H, N and O under different temperature conditions. The ingredients were extensively analyzed to ensure that they could be identified with high accuracy. To examine the spectral lines associated with different elements according to their wavelength and energy, the data were gathered using a laser source and spectrometer. The initial sample under test was a standard 1000 mg paracetamol tablet kept under the same conditions regarding the temperature (room temperature) and humidity. The spectra generated by LIBS were used to identify the substances contained in this medication. Regardless of the laser's direction, the emission spectra were captured using a spectrometer positioned at a 90° angle to the sample's surface.

Figure 12a shows the spectra produced by LIBS to identify the elements present in a sample containing active agents and additional elements as inactive agents. This is the first investigation of such a sample. The identified elements are crucial in understanding the composition of the drugs. Furthermore, as shown in Fig. 12b, the sample was located in the LIBS instrument during the investigation. Fig. 13 shows a typical spectrum for each element present in the samples, displayed individually for the sake of exposition. Peaks belonging to nitrogen, oxygen, hydrogen, and carbon could be seen in every spectrum. These lines can be linked to the main constituents (active medicinal substances) of the drug, which are organic compounds. We believe that the additional metals might be impurities or sodium-containing excipients. The findings of the required constituents of the paracetamol samples are properly matched to the results from the NIST database in Table 4, compared to the observed peaks with LIBS in Fig. 11a. The chemical formula of paracetamol is C8H9NO2, which means that the molecule is composed of 8 carbon atoms associated with 9 hydrogen atoms, 1 nitrogen atom, and 2 oxygen atoms. The same procedure as the microwave measurement was applied to characterize the drug at different temperature levels. We heated the drug at 80 °C for 15 min.

Figure 13 shows different peak levels for the three different temperature conditions. The black line indicates the peak corresponding to the drug at a room temperature of 27 °C (before heating). The red and blue lines refer

<figure><img src="image_13.png" /><figcaption>Figure 11. Comparison of different heating temperature (80 °C, 60 °C and 40 °C).</figcaption></figure>
<table><caption>**Table 4.** Peaks and their corresponding atomic elements (NIST DATABASE).</caption><thead><tr><td>Element</td><td>Wavelength (nm)</td></tr></thead><tbody><tr><td>Carbon</td><td>247.85</td></tr><tr><td>Hydrogen</td><td>656.27</td></tr><tr><td>Nitrogen (N1, N2)</td><td>742.36, 744.23</td></tr><tr><td>Oxygen</td><td>777.53</td></tr></tbody></table>

<figure><img src="image_14.png" /><figcaption>**Figure 12.** (a) LIBS spectra of the samples used in the study. (b) Drug location in LIBS.</figcaption></figure>

to a heated drug at 80 °C and the drug after heating while it returns to 27 °C, respectively. Figure 13 illustrates the changes in the atomic emission line intensities for carbon, hydrogen, nitrogen (N1, N2), and oxygen upon varying the heating temperature.

The parameters presented below are defined to qualitatively analyze the performance of the LIBS instrument. The intensity quality shifting $Q_i$ is expressed as:

$$Q_i = I_{\text{ref}} - I_{\text{heated}} \quad (\text{MHz}) \quad (4)$$

where $I_{\text{ref}}$ is the measured intensity of paracetamol at room temperature before heating. This intensity is considered as the reference. $I_{\text{heated}}$ is the measured intensity of the drug heated at 80 °C.

It is notable from the values shown in Table 5 taken from peaks that carbon and hydrogen have large $Q_i$ values of 8884 and 812, respectively, between the intensities of the normal drug (before heating) and after heating at 80 °C, unlike nitrogen (N1, N2) and oxygen, which provide small $Q_i$ values of approximately 16/140 and 198, respectively. Additionally, paracetamol consists of 8 atoms of carbon and 9 atoms of hydrogen, which are the most common elements in the structure. Thus, we consider carbon and hydrogen as the representative elements to characterize the phenomena in this work. The measured values shown in Fig. 12 represent the average of three repetitions; therefore, the risk of random fluctuations is small.

It is clear in Fig. 14a,b that changes in temperature heating, which modifies the structure and dipole moment of water, definitely have an impact on the intensity of atomic emission lines from carbon and hydrogen. Increased heat gradients may create significant electrostatic fields that change the energy radiated by the sample, leading to an alteration in intensity level. This is due to the polarization of water induced by temperature fluctuations, which is an out-of-equilibrium phenomenon. The maximum intensity of both elements is observed during the measurement of the drug under heating at 80 °C (red line), whereas the electric field induced by the temperature gradient should be larger. Otherwise, the intensity of the drug after heating (blue line) decreased because of the drop in temperature.
<figure><img src="image_15.png" /><figcaption>Figure 13. Different peaks for each detected element.</figcaption></figure>

<table><caption>Table 5. Intensities for peaks of each elements in different condition).</caption><thead><tr><td>ElementsConditions</td><td>Before heating 27°C</td><td>Heated 80°C</td><td>After heating 27°C</td><td>Q<sub>i</sub></td></tr></thead><tbody><tr><td>Carbon</td><td>33,396</td><td>50,673</td><td>42,280</td><td>8884</td></tr><tr><td>Hydrogen</td><td>1952</td><td>3144</td><td>2764</td><td>812</td></tr><tr><td>Nitrogen (N1/N2)</td><td>284/405</td><td>222/526</td><td>300/545</td><td>16/140</td></tr><tr><td>Oxygen</td><td>1400</td><td>1866</td><td>1202</td><td>198</td></tr></tbody></table>

## Conclusion

In this paper, an SRR-loaded microwave planar sensor is proposed for tracking drug quality under different heat treatment conditions. A triple hexagonal split ring resonator (TH-SRR) is loaded in a circular patch, which is the sensing area in our system. The evaluation of sensor measurements centers around the microwave quality shifting ($Q_{μW}$), a parameter indicating how microwaves change behavior with shifting temperatures. This shift serves as an indicator of potential variations in the drug's properties. By analyzing these shifts, we indirectly assess if the drug has been affected by temperature-induced changes, providing insights into its stability and potential damage. The TH-SRR-based sensor has the benefit of a very strong and confined intensity between the resonator edges, which allows for a good interpretation of the smallest variations in drug properties caused by different heating conditions from 40 to 80 °C. By locating the drug in this region and evaluating the change in the frequency response of the sensor, the microwave quality shifting $Q_{μW}$ is determined. In particular, heating a drug at 80 °C has shown a noticeable $Q_{μW}$ of approximately 1 MHz. The sensor exhibits $Q_{μW}$s of 300 kHz and 50 kHz at 60 °C and 40 °C, respectively. Additionally, a comparison study is conducted using an analytical technique for determining the elemental composition of the drug. It is shown that LIBS presents an outstanding intensity quality shift of 8884 and 812 for carbon and hydrogen, respectively. Moreover, it demonstrates good
<figure><img src="image_16.png" /><figcaption>**Figure 14.** Comparison between $Q_{\mu W}$ and $Q_i$ for (a) Carbon (b) hydrogen.</figcaption></figure>

agreement with the microwave quality shifting $Q_{\mu W}$. The findings of this study could be regarded as a model shift in microwave sensing for customized pharmaceutical applications such as monitoring the quality of drugs, which opens the door for their commercial exploitation.

## Data availability

The datasets generated and/or analyzed during the current study are available in the OneDrive Cloud repository, you can access all the data through this link https://um6p-my.sharepoint.com/:f:/g/personal/youness_zaarour_um6p_ma/EshQOEITZh9JoT8e9DhsPuMB8mlh8soj82SRejNNC7C4mQ?e=rfa61P.

Received: 5 May 2023; Accepted: 23 September 2023

Published online: 14 October 2023

## References

1. Ağca, A. C. et al. Antioxidant, anti-inflammatory and antidiabetic activity of two *Sternbergia* taxons from Turkey. *South Afr. J. Bot.* **136**, 105–109 (2021).

2. On Specifications for Pharmaceutical Preparations, W. E. C., World Health Organization et al. WHO Expert Committee on Specifications for Pharmaceutical Preparations: Thirty-Sixth Report (World Health Organization, 2002).

3. Qu, J. et al. A nationwide survey exploring physicians' and pharmacists' knowledge, awareness and perceptions regarding generic medicines in China. *BMC Health Serv. Res.* **22**, 1–11 (2022).

4. Diener, H., Pfaffenrath, V., Pageler, L., Peil, H. & Aicher, B. The fixed combination of acetylsalicylic acid, paracetamol and caffeine is more effective than single substances and dual combination for the treatment of headache: A multicentre, randomized, double-blind, single-dose, placebo-controlled parallel group study. *Cephalalgia* **25**, 776–787 (2005).

5. Driad, Y. *Stabilité du paracétamol: Application à un sachet produit en industrie pharmaceutique*. Ph.D. thesis, UHP-Université Henri Poincaré (2009).

6. Ahmad, I. & Shaikh, R. H. Effect of temperature and humidity on hardness and friability of packaged paracetamol tablet formulations. *Pak. J. Pharm. Sci.* **7**, 69–78 (1994).

7. Bate, R., Tren, R., Hess, K. & Attaran, A. Physical and chemical stability of expired fixed dose combination artemether-lumefantrine in uncontrolled tropical conditions. *Malaria J.* **8**, 1–7 (2009).

8. Terzyk, A. P., Rychlicki, G., Biniak, S. & Łukaszewicz, J. P. New correlations between the composition of the surface layer of carbon and its physicochemical properties exposed while paracetamol is adsorbed at different temperatures and pH. *J. Colloid Interface Sci.* **257**, 13–30 (2003).

9. Hodel, E. M., Kaur, H. & Terlouw, D. J. Stability of dihydroartemisinin-piperaquine tablet halves during prolonged storage under tropical conditions. *Am. J. Trop. Med. Hyg.* **96**, 338 (2017).

10. Wei, X.-C. et al. Recent advances of novel technologies for quality consistency assessment of natural herbal medicines and preparations. *Chin. Med.* **15**, 1–24 (2020).

11. Alcalà, M. et al. Qualitative and quantitative pharmaceutical analysis with a novel hand-held miniature near infrared spectrometer. *J. Near Infrared Spectrosc.* **21**, 445–457 (2013).

12. Parys, W., Dołowy, M. & Pyka-Pająk, A. Significance of chromatographic techniques in pharmaceutical analysis. *Processes* **10**, 172 (2022).

13. Gupta, A., Austin, J., Davis, S., Harris, M. & Reklaitis, G. A novel microwave sensor for real-time online monitoring of roll compacts of pharmaceutical powders online-a comparative case study with nir. *J. Pharm. Sci.* **104**, 1787–1794 (2015).

14. Omer, A. E. et al. Low-cost portable microwave sensor for non-invasive monitoring of blood glucose level: Novel design utilizing a four-cell csrr hexagonal configuration. *Sci. Rep.* **10**, 1–20 (2020).

15. García, H., Juan, C. G., Ávila-Navarro, E., Bronchalo, E. & Sabater-Navarro, J. M. Portable device based on microwave resonator for noninvasive blood glucose monitoring. In *2019 41st Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC)*, 1115–1118 (IEEE, 2019).

16. Huang, S. Y. et al. Microstrip line-based glucose sensor for noninvasive continuous monitoring using the main field for sensing and multivariable crosschecking. *IEEE Sens. J.* **19**, 535–547 (2018).

17. Zarifi, M. H. et al. Microbead-assisted high resolution microwave planar ring resonator for organic-vapor sensing. *Appl. Phys. Lett.* **106**, 062903 (2015).

18. Walter, S., Bogner, A., Hagen, G. & Moos, R. Novel radio-frequency-based gas sensor with integrated heater. *J. Sens. Sens. Syst.* **8**, 49–56 (2019).

19. Park, J.-K. et al. Real-time humidity sensor based on microwave resonator coupled with PEDOT: PSS conducting polymer film. *Sci. Rep.* **8**, 439 (2018).
20. Aldhaeebi, M. A., Almoneef, T. S., Attia, H. & Ramahi, O. M. Near-field microwave loop array sensor for breast tumor detection. *IEEE Sens. J.* **19**, 11867–11872 (2019).

21. Wang, L. Microwave sensors for breast cancer detection. *Sensors* **18**, 655 (2018).

22. Abegaonkar, M. P., Karekar, R. & Aiyer, R. A microwave microstrip ring resonator as a moisture sensor for biomaterials: Application to wheat grains. *Meas. Sci. Technol.* **10**, 195 (1999).

23. Jafari, F., Khalid, K., Yusoff, W. D. W. & Hassan, J. The analysis and design of multi-layer microstrip moisture sensor for rice grain. *Biosyst. Eng.* **106**, 324–331 (2010).

24. Sharafadinzadeh, N., Abdolrazzaghi, M. & Daneshmand, M. Highly sensitive microwave split ring resonator sensor using gap extension for glucose sensing. In *2017 IEEE MTT-S International Microwave Workshop Series on Advanced Materials and Processes for RF and THz Applications (IMWS-AMP)*, 1–3 (IEEE, 2017).

25. Chee, S. N., Johansen, A. L., Gu, L., Karlsen, J. & Heng, P. W. S. Microwave drying of granules containing a moisture-sensitive drug: A promising alternative to fluid bed and hot air oven drying. *Chem. Pharm. Bull.* **53**, 770–775 (2005).

26. Maniruzzaman, M. et al. Taste masking of paracetamol by hot-melt extrusion: An in vitro and in vivo evaluation. *Eur. J. Pharm. Biopharm.* **80**, 433–442 (2012).

27. Marty, W., Sigrist, T. & Wyler, D. Temperature variations in automobiles in various weather conditions: An experimental contribution to the determination of time of death. *Am. J. Forensic Med. Pathol.* **22**, 215–219 (2001).

28. Almsellati, N. F., Swesi, O. & Aladuli, S. Quantitative analysis of pharmaceutical tablets using libs-technique. *Int. J. Dev.* **8**, 1–7 (2019).

29. Le Garrec, S., Burnat, P. & Gentes, P. Le paracétamol. *Lyon pharmaceutique* **45**, 227–242 (1994).

30. De los Santos, J., Garcia, D. & Eiras, J. A. Dielectric characterization of materials at microwave frequency range. *Mater. Res.* **6**, 97–101 (2003).

31. Heng, P. W., Loh, Z., Liew, C. V. & Lee, C. Dielectric properties of pharmaceutical materials relevant to microwave processing: Effects of field frequency, material density, and moisture content. *J. Pharm. Sci.* **99**, 941–957 (2010).

32. Nisar, S., Dastgeer, G., Shafiq, M. & Usman, M. Qualitative and semi-quantitative analysis of health-care pharmaceutical products using laser-induced breakdown spectroscopy. *J. Pharm. Anal.* **9**, 20–24 (2019).

33. Wiltshire, B. D., Rafi, M. A. & Zarifi, M. H. Microwave resonator array with liquid metal selection for narrow band material sensing. *Sci. Rep.* **11**, 8598 (2021).

34. Zhou, H. et al. Multi-band sensing for dielectric property of chemicals using metamaterial integrated microfluidic sensor. *Sci. Rep.* **8**, 14801 (2018).

35. Baena, J. D. et al. Equivalent-circuit models for split-ring resonators and complementary split-ring resonators coupled to planar transmission lines. *IEEE Trans. Microw. Theory Tech.* **53**, 1451–1461 (2005).

36. Nouretdinov, I., Melluish, T. & Vovk, V. Ridge regression confidence machine. In *ICML*, 385–392 (Citeseer, 2001).

37. Antoniadis, A., Lambert-Lacroix, S. & Poggi, J.-M. Random forests for global sensitivity analysis: A selective review. *Reliab. Eng. Syst. Saf.* **206**, 107312 (2021).

38. Bickel, P. J. One-step Huber estimates in the linear model. *J. Am. Stat. Assoc.* **70**, 428–434 (1975).

## Acknowledgements

This research was supported by the PhD Student Grant of the University of Mohammed VI Polytechnique (UM6P).

## Author contributions

Y.Z., H.G. and A.F. wrote the main manuscript. Preparation of the samples and processing of laboratory tests were done by F.E.A. and Y.Z. Critical analysis and deep interpolation of data were done by R.E.A. and W.S. LIBS Measurement, data collection and LIBS analysis were performed by M.K., H.G. and Y.Z. All authors reviewed the manuscript. Prof. A.F. and Prof. M.D are the main supervisors and the co-supervisor of the Ph.D. Student Y.Z.

## Competing interests

The authors declare no competing interests.

## Additional information

Correspondence and requests for materials should be addressed to Y.Z.

Reprints and permissions information is available at www.nature.com/reprints.

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_17.png" /></figure>

Open Access This article is licensed under a Creative Commons Attribution 4.0 International

License, which permits use, sharing, adaptation, distribution and reproduction in any medium or

format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the
Creative Commons licence, and indicate if changes were made. The images or other third party material in this
article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the
material. If material is not included in the article's Creative Commons licence and your intended use is not
permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from
the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2023