Temperature-dependent complex dielectric permittivity: a simple measurement strategy for liquid-phase samples

Microwaves (MWs) are an emerging technology for intensified and electrified chemical manufacturing. MW heating is intimately linked to a material’s dielectric permittivity. These properties are highly dependent on temperature and pressure, but such datasets are not readily available due to the limited accessibility of the current methodologies to process-oriented laboratories. We introduce a simple, benchtop approach for producing these datasets near the 2.45 GHz industrial, medical, and scientific (ISM) frequency for liquid samples. By building upon a previously-demonstrated bireentrant microwave measurement cavity, we introduce larger pressure- and temperature-capable vials to deduce temperature-dependent permittivity quickly and accurately for vapor pressures up to 7 bar. Our methodology is validated using literature data, demonstrating broad applicability for materials with dielectric constant ε' ranging from 1 to 100. We provide new permittivity data for water, organic solvents, and hydrochloric acid solutions. Finally, we provide simple fits to our data for easy use.


Methods
The experimental setup employed in this study is a modified system 39 , an autonomous system capable of acquiring dielectric properties of materials near the 2.45 GHz ISM band.This design allows convenient measurement of materials inside tubes, particularly suitable for liquids, semisolids, powders, and granular materials.The MW structure is based on a single-post re-entrant cavity 40 , CNC-machined from aluminum (± 5 µm precision), and featuring an open port on the top cover to accommodate the tubes along the central axis of the cell.The measurement procedure relies on the precise measurement of the MW cavity's reflection or impedance to determine the resonant frequency and quality factor and the subsequent use of the resonance parameters to calculate the dielectric properties of the material from a rigorous numerical method 39 .Permittivity characterization with MW cavities is typically restricted to low-loss materials 41 .However, high-coupling networks enable measurements of both low-loss and high-loss materials using the same experimental setup.To streamline and simplify the measurement process, the device includes a tailor-made MW reflectometer, which assesses the MW cavity's reflection coefficient, eliminating the need for a full-featured VNA.
The developed electromagnetic (EM) numerical method of the MW cavity to calculate the dielectric properties of the materials made use of circuit analysis and mode-matching techniques 42 , as described in detail in previous work 39 .This technique involves decomposing the complex MW cavity into simpler canonical circuits, which are analyzed and computed independently and then connected to assemble the complete structure.The relationship between the complex resonance frequency (resonant frequency and quality factor) and complex permittivity (dielectric constant and loss factor) is solved by applying the resonance condition to the numerical response of the MW cavity.Unlike other dielectric properties measurement methods, such as the cavity perturbation methods (CPM), the approach described in reference 39 does not require calibration with reference materials.We refer readers to our previous demonstration of such a system for a full description 39 and that we build upon here.
In the previous setup, commercial Pyrex tubes were utilized to house the materials.However, the permittivity of Pyrex is highly temperature-dependent 43 .As the temperature rises, the resonance measurement becomes influenced by both the changes in material permittivity and the variations in the tube's permittivity.This interplay complicates the development of an accurate EM model and its practical implementation.On the other hand, while the small size of the Pyrex tubes facilitates the use of minimal sample quantities, it limits the incorporation of external elements for high-pressure containment or atmospheric control.To overcome these limitations, an upgraded MW cavity has been developed to accommodate more convenient larger tubes, with the dimensions of the tubes and cavity described in Fig. 1.These tubes are crafted from materials like quartz and alumina, which exhibit more consistent permittivity across different temperatures 44,45 .
An EM numerical model was developed in MATLAB and employed to ascertain the optimal dimensions of the bireentrant resonator, aiming to achieve a sensitivity level comparable to that of the Pyrex tube version, attaining a single resonant mode encompassing the operational range of the MW reflectometer, spanning from 1.5 to 2.6 GHz. Figure 1 depicts a cross-sectional view of the resonant geometry and the dimensions.The relationship between the permittivity and the resonant frequency of the resulting bireentrant MW structure is shown in Fig. 2. Figure 2a illustrates the frequency shift of the cavity loaded with quartz tubes, and Fig. 2b corresponds to the frequency variation using alumina tubes.The resonance frequencies of the cavity using alumina tubes for the same dielectric material are lower than those using quartz tubes, attributed to the higher dielectric constant of alumina.Nevertheless, in both cases, a maximum frequency shift of approximately 600 MHz is achieved for materials with a dielectric constant ranging from 1 to 100, like the frequency shift obtained with the Pyrex tubes in the original bireentrant setup 39 .
The described MW setup allows for flexible sample handling.A vessel was then developed to facilitate highpressure containment of liquid samples at elevated temperatures while still allowing for accurate measurement of dielectric permittivity.To maintain permittivity uncertainties within a 2% margin for the dielectric constant and a 5% margin for the loss factor, the uncertainty associated with the internal diameter of the tube, with respect to the nominal value used in the model, was determined to be ± 20 µm 39 .The pressure-temperaturedielectric (PTD) vessel is illustrated in Fig. 3 and fulfills the requirements as follows: Alumina vials (Fig. 3c(iii)) are designed to fit the access hole of the bireentrant cavity (iv) and are equipped with a flanged top for sealing via an o-ring fixture.A threaded assembly (ii) ensures pressure containment within the vessel.Additionally, a temperature-sensing optical fiber (i) with a silica cladding is inserted and sealed within the vial using soft ETFE ferrules, as described in prior work 31,46 .Specific dimensions of the alumina vial are offered in Fig. S1.The PTD vessel is heated to a predetermined temperature using an external heating block, after which it is transferred to the cavity.As the vessel cools, permittivity-temperature data points are collected.Throughout this process, the pressure changes to reflect the vapor pressure of the solvent at each temperature.We do not expect this to affect the measurements for sub-critical liquids as demonstrated by Dimitrakis and coworkers 35 who showed no difference in the permittivity for alcohols at pressures ranging from 1 to 60 bar.

Validation of permittivity measurements
Although the MW setup can operate with both quartz or alumina tubes which result in nearly identical measurements (Fig. S2), quartz tubes are limited to near atmospheric pressures.This study focuses on alumina vials since they have superior pressure resistance and can reliably be used up to 7 bar, providing access to wider temperature ranges for liquid permittivity measurements.
To validate the temperature measurement scheme, the dielectric permittivity of ultrapure water was measured from 20 to 160 °C.While some references in the literature report water permittivity up to 100 °C33,34 , the frequency at which they provide the permittivity values is significantly different.This discrepancy hinders the proper comparison of permittivity values due to the strong dependency of permittivity with frequency for polar materials 47 .Then the measured permittivity was compared to the values given by the Debye-type relaxation spectral function provided by Kaatze 48 , allowing the calculation of the permittivity of water as a function of frequency and temperature up to 60 °C.Figure 4 compares our measurements taken at 1.85 GHz to the permittivity values obtained from Kaatze's Debye model at the same frequency.Our experimental results align remarkably well with those values throughout the entire temperature range of the reference 48 .This agreement in data demonstrates the successful efficacy of our temperature measurement scheme in accurately capturing the temperature-dependent dielectric permittivity.
To demonstrate the efficacy of alumina vials for measuring various solvents across a wide range of ε′ values, we conducted room temperature dielectric constant measurements as depicted in Fig. 5.We compare our experimental results to literature values obtained from different sources (refer to Table S1 for details).Notably, certain solvents exhibit a significant ε' variation depending on the measurement frequency according to literature data (namely alcohols).After accounting for this, the agreement between our experimental measurements and the reported literature values is excellent (an average of 5.8% relative error) across the entire ε' range of approximately 1 to 100.Unfortunately, literature data for solvents at elevated temperatures and the same frequency as our measurements are lacking, so additional comparisons cannot be made.
In comparison to the most contemporary approach 35 , our experimental setup does not require a full-featured VNA, an HPLC pump, a dedicated temperature controller, or a customized high-pressure adapter for the coaxial line.Rather, temperature is continuously monitored internally in the alumina vial which is placed on a hot plate for heating and then removed and placed into the measurement cell at the desired time.This process can be repeated quickly for various liquid samples, and the experimental data presented in this work was collected in a single day.However, our approach can only reach moderate pressures (~ 7 bar) with the current design which can limit the accessible temperature range.Furthermore, our method does not allow for traceable measurements which may be important given the sensitivity to frequency exhibited by certain compounds.The low experimental variance in our measurements also represents an improvement, possibly owing to our seal placement scheme which avoids the introduction of geometric changes to the measurement cell upon pressurization 35 .

Dielectric permittivity of non-alcohol organic solvents
We tested a variety of organic solvents to investigate their temperature-dependent permittivity behavior and establish a small repository of such data.Most solvents generally exhibit a consistent, monotonic decrease in both ε' and ε" as temperature increases.Figure 6a and Fig. 6b showcase this trend with representative examples of MIBK and ethyl acetate permittivity.For such cases, a straightforward three-parameter exponential fit, as presented in Eq. (1), provides an accurate, empirical model of permittivity across the entire range of measured temperatures.
Here, A 1 , A 2 , and A 3 are the regressed parameters, T is the measured temperature [°C], and ε is either the real or imaginary component of the dielectric permittivity.The temperature-dependent permittivity data for each non-alcohol organic solvent is listed in Figs.S3 and S4, while the regressed parameters A 1 , A 2 , and A 3 from Eq. (1) are available for each solvent in Table 1.The standard error of the fitted Eq. (1) for each solvent is less www.nature.com/scientificreports/than 3% in all cases.The suitable temperature range for each solvent was chosen to limit the vapor pressure to 7 bar for safe operation.

Dielectric permittivity of alcohols
The permittivity of alcohols displays a nonmonotonic dependence on temperature 35,60 .As temperature increases, so do ε′ and εʺ until the trend is reversed at some point and ε′ and εʺ decrease with further increments of temperature.This effect is highlighted in Fig. 7a and Fig. 7b where ethanol and isopropanol permittivity are displayed.In these cases, a modified lognormal fit is used with a total of 6 regressed parameters as shown in Eq. ( 2) where B 1 , B 2 , B 3 , B 4 , σ, and µ are the regressed parameters.In this case, σ and µ are the traditional descriptors of a lognormal fit while B 1 , B 2 , B 3 , and B 4 allow x-axis and y-axis stretching and offset.The temperature-dependent permittivity data for each alcohol are listed in Fig. S5, while the regressed parameters B 1 , B 2 , B 3 , B 4 , σ, and µ from Eq. ( 2) are available for each solvent in Table 2.The standard error of the fitted Eq. ( 2) for each solvent is less than 3.5% in all cases.

Dielectric permittivity of aqueous hydrochloric acid solutions
Next, we measured permittivity of water and hydrochloric acid solutions of various strengths.Figure 8a shows the permittivity for a continuum of HCl solutions ranging from pH 7 to 1. ε′ for pure water descends monotonically with temperature.Up to pH ≈ 1.6, HCl addition has little effect on εʺ.However, for stronger and stronger HCl solutions, smaller and smaller ε' drops are observed with increasing temperature.Finally, at pH 1, ε′ increases with increasing temperature.ε′ at room temperature decreases slightly by lowering pH.In Fig. 8b, εʺ is reported for the same solutions.HCl addition monotonically increases εʺ-even at room temperature.As HCl is added, the decrease in εʺ with temperature is less prominent until pH 2, where increases in temperature cause increases Table 1.Regressed parameters for non-alcohol organic solvents fitted to Eq. (1).The standard error relative to the mean is reported for each parameter set along with the maximum temperature measured for each solvent.www.nature.com/scientificreports/ in εʺ.This behavior in which ionic concentration changes the temperature effect on permittivity has not been previously reported to our knowledge.The permittivity data for HCl solutions is fitted using Eq. ( 2).The regressed parameters for each HCl solution are presented in Table 3.In each case, the standard error of the fitted Eq. ( 2) is less than 4%.

Conclusions
We demonstrate a simple benchtop approach for evaluating temperature-dependent dielectric permittivity near the 2.45 GHZ ISM frequency, constituting a substantial simplification of existing equipment and ease of use.By building upon the existing benchtop standalone dielectric measurement kit with larger, easier-to-use vials, a simple pressurization scheme was implemented.Improvements in vial materials were also made by tuning the electromagnetic model for quartz and alumina tubes, which exhibit less sensitivity of dielectric permittivity to temperature.Ultimately, liquids can be heated via a hotplate to temperatures significantly above their boiling point and up to 7 bar, and the pressurized and heated vial can be transferred to the measurement cavity for permittivity evaluation.The methodology was validated using the known temperature-dependent permittivity of water and the room temperature permittivity of various solvents.Excellent agreement with the literature data was demonstrated.
We then report the permittivity of water, alcohols, organic solvents, and hydrochloric acid solutions across temperatures up to 7 bar of vapor pressure.While water and non-alcohol organic solvents display monotonic Table 2. Regressed parameters for alcohols fitted to Eq. ( 2).The standard error relative to the mean is reported for each parameter set along with the maximum temperature measured for each solvent.decreases in ε′ and εʺ with increasing temperatures, more surprising behavior emerges for other compounds.Alcohols seemingly share a separate behavior in which ε′ and εʺ increase to a maximum and then decrease with increasing temperature.Hydrochloric acid solutions also exhibit non-monotonic behavior which changes across pH values.Based on the permittivity behavior for the various categories of liquids, exponential and lognormal fits were provided.

Solvent
Table 3. Regressed parameters for HCl solutions fitted to Eq. ( 2).The standard error relative to the mean is reported for each parameter set along with the maximum temperature measured for each solvent.

Figure 2 .
Figure 2. Frequency shift as a function of permittivity for the designed bireentrant cavity.(a) Dielectric materials enclosed in quartz tubes and (b) dielectric materials inside alumina tubes.Contours are placed at 0.1 GHz frequency increments.

Figure 3 .Figure 4 .
Figure 3. Depiction of the PTD vessel.The external (a) and internal cross section (b) views of the vessel are shown as designed for machining, illustrating the simplicity of the approach.A photograph of the vessel (c) with the temperature-sensing optical fibers (i), screw-topped sealing fixture (ii), alumina vials (iii), and the cavity access hole (iv).

Figure 5 .Figure 6 .
Figure 5. Parity of measured and literature ε′ values for various solvents at room temperature utilizing the alumina PTD vessel.Including data from references 49-59 .

Figure 7 .
Figure 7. Sample permittivity profiles for ethanol (a) and isopropanol (b).Grey squares represent ε′ (left axis) while purple diamonds represent εʺ (right axis).Solid and dashed lines represent the lognormal fits for ε′ and εʺ data, respectively.

Figure 8 .
Figure 8. Temperature-dependent dielectric permittivity of hydrochloric acid solutions including dielectric constant (a) and loss factor (b).