A. Mathur is a graduate student at the Sensors, Energy, and Automation Laboratory, Department of Electrical Engineering, University of Washington.

A. Mathur

process analytical technology (PAT) initiative

 reflects the need for improved quality control in the pharmaceutical industry. As the complexity of pharmaceutical products increases, quality control becomes more costly and difficult. New analytical instrumentation should be developed to address the quality control problems associated with more complex pharmaceutical products. Dielectric spectroscopy, one of many process analytical technologies, has been used for several decades. As any other technology, it is continuously evolving. New instrumentation, algorithms, and materials broaden the applicability and capabilities of this method. Selecting the optimal architecture for a sensing system intended for each specific application requires a good understanding of main principles and engineering trade-offs common for this field. This article is written primarily for field practitioners who need to make a judicious choice of the best measurement technique for their application. Main principles and design features of dielectric sensors are discussed in the framework of pharmaceutical applications.  

The range of potential applications of dielectric spectroscopy is quite broad. Virtually any physical process change leads to changes in dielectric properties of samples. Process variability is a primary concern for the pharmaceutical industry (1). Exposure to mechanical and thermal stress can cause a change in the physical properties of pharmaceuticals. Such variations are important to control because physical properties generally determine the efficacy of the drug. For instance, tablet coatings control the rate of drug delivery within the body of a patient (2–5) and influence the bioavailability of the drug. In addition, coatings protect the active ingredient from chemically harsh environments in the body (6). Similarly, API content in a given sample determines the potency of the drug. The drying process of pharmaceuticals is critical because 70% of global granulated pharmaceutical product is made using wet granulation (7). In wet granulation, liquid binding gels are used to facilitate bonding between active ingredients. At this stage, it is important to measure moisture because specific moisture levels are required for the formation of the correct-size granules. Further, in high-shear wet granulation methods, incorrect moisture levels can lead to a process-induced transformation (PIT) (8). In these transformations, the properties of the active ingredients can change, resulting in reduced efficacy.  

 A number of sensing technologies can be used to detect physical properties of pharmaceuticals. Most commonly used technologies include near infrared spectroscopy (NIR) and Raman spectroscopy. Digital imaging methods, optical methods, and dielectric spectroscopy also are used to measure various physical properties.  

 NIR spectroscopy is widely used to measure API content, distribution of contaminants, moisture content, and polymorphism determination (9). Similarly, Raman spectroscopy has demonstrated success in measuring pharmaceutical reaction times, polymorphism, and differences between solid-state forms (8). 

 Digital imaging methods are used to monitor particle size in pharmaceutical powders (10). Optical methods have been used for inspection of tablet coatings (6), measurement of constituent concentration of pharmaceutical powder mixtures (11), and quantitative analyses of ascorbic acid in pharmaceuticals (11, 12). For example, laser induced breakdown spectroscopy (LIBS) is used to measure coating thickness of tablets (13). 

Previous work using dielectric spectroscopy.

 Dielectric spectroscopy is a promising method to study physical properties of pharmaceutical solutions, colloids, microcapsules, gels, and emulsions. Craig provides a thorough review of these applications (14), which are summarized in this section. 

Pharmaceutical solutions are important because the solubility of a drug in the solvent directly influences its rate of disassociation in the body, which in turn impacts bioavailability. Dielectric properties such as the static dielectric constant of solutions has been used to study the solubility of cosolvent systems, reaction rates of pharmaceutical solutions, and drug stability (14).   

Colloids are two-phase systems that consist of particles suspended in a medium. Each particle in a colloidal solution has an electrical layer surrounding it that determines the stability of the system by affecting the aggregation or disassociation of colloidal particles. As a result of these forces, dielectric response of colloidal solutions provides information about particle size in the colloidal solutions. For example, a study conducted by Paul and Vogey suggested an inverse-square dependency between particle radius and dielectric relaxation frequency of the medium (15). This relationship is illustrated in the following equation:  

 is the diffusion coefficient of counterions with a particle radius of 

 is the loss peak frequency. The inverse-square relationship between relaxation frequency and particle radius was confirmed in a subsequent study (16). 

Microcapsules have attracted attention as an improved drug delivery method. Coating thickness and uniformity of microspheres directly influences the rate of dissolution in the body, making it very important to study its structural properties. Previous studies show methods to relate the relative dielectric permittivity and thickness of the capsule wall (17). Additional information provided by dielectric measurements includes volume fraction of microspheres in solution and the dielectric properties of the interior core of microcapsules (14).  

Dielectric spectroscopy has shown promise for the measurement of structural and interfacial properties of gels. Gels are complex systems with a number of dielectric relaxation processes that can be related to unique physical characteristics. For example, previous studies showed the response of aqueous gels containing cetosteryl alchohol and cetrimide show four distinct bulk processes (18), These responses were attributed to the relaxation of water, ions in bilayers formed resulting from the presence of cetosteryl alchohol, and migration of ions in channels between these bilayers. Therefore, dielectric spectroscopy can be used to monitor multiple physical properties of complex structures. Dielectric spectroscopy was used to study the likelihood of high surface contact between nasal mucosa and pharmaceutical gels by measuring conductivity as pharmaceutical drug ions transported from the gel into the nasal mucous membrane (19). 

Emulsions and microemulsions have attracted interest as an effective drug delivery method for drugs that are insoluble in water. Dielectric spectroscopy has been used to classify the types of emulsions. For example, dielectric properties are used to distinguish between water-in-oil and oil-in-water emulsions (14).  

The remainder of this article discusses key issues in dielectric sensor selection. Measurements of tablet coating thickness are used to illustrate methods of dielectric data processing.  

There are two major properties of dielectric spectroscopy that are typically varied to suit the desired application: the spatial distribution of the interrogation field and the interrogation frequency range. Changes to the field's spatial distribution of the field used in dielectric spectroscopy enables measurements in the following modes: bulk, surface/low-penetration, multiple penetration depth, imaging, and a combination of these. The frequency range can be optimized for a particular material of interest and can vary from 10

  Hz, although a frequency range of only a few orders of magnitude is typically used.  

Based on the measurement constraints for the application, the spatial distribution of the field and the frequency range are ascertained by answering a series of questions, which are discussed in the following sections. 

 In dielectric spectroscopy, the attenuation and phase shift of an electric field at specific frequencies are measured as it passes through a material of interest. Attempts are typically made to either linearize the field or to limit the penetration of the field to a known depth into the material of interest. In general, through-field implies a "pitch and catch" configuration in which one electrode delivers an excitation signal and another electrode measures the resulting signal after passing in an approximately linear fashion through the material of interest such as in a parallel plate capacitor configuration (see Figure 1a). In the case of fringing-field sensing, one electrode delivers an excitation signal, and a second electrode is placed such that the original signal passes through the material of interest in a nonlinear fashion. The most common fringing field configuration is one in which the excitation and sensing electrodes are mounted in a coplanar fashion on the same substrate (see Figure 1b). 

Figure 1 (ALL FIGURES ARE COURTESY OF THE AUTHORS.) 

Through-field or fringing-field sensing techniques are typically chosen on the basis of the following: only one sided access to sample is available, significant variations in thickness exist, and the sensor must operate in an electrically noisy environment. To determine whether through-field or fringing-field is desired for an application, access to the material of interest must first be determined. If two-sided access is possible and the sample thickness does not vary significantly, through-field is recommended. If only one-sided access is available or if the sample thickness varies significantly, then a fringing-field measurement is recommended. To obtain an effective measurement with either sensor type, a basic knowledge of the system is important. The nominal thickness of the sample should be known as well as the primary sample constituents and the likelihood of metallic impurities. 

A common electrode arrangement for through field imaging is seen in electrical tomography (ET). In this method, the electrodes are distributed around the edges of a cylindrical surface. Individual electrodes are electrically excited and the current or voltage induced on other electrodes is measured. The measured values are used to reconstruct distribution of physical properties in the material by solving the inverse problem. ET is used for biomedical, geophysics, and industrial process control applications (20–22). Compared with other imaging methods such as X-ray tomography and magnetic resonance imaging, ET is fast, noninvasive, and inexpensive in instrumentation. Although ET is useful for measurements in cylindrical surfaces, it is not viable for scenarios in which one-sided measurements are required.  

In the case of fringing-field measurements, a common example is referred to as a fringing electric field (FEF) sensor. An FEF sensor can perform three distinct measurement functions depending on how it is configured: defined penetration depth measurements, multiple penetration depth measurements, and one-sided measurements.  

In the case of the FEF sensor, the field is nonlinear and often has smaller signal strength than that of the through-field configuration. However, the field is concentrated much closer to the electrodes and therefore has a defined penetration depth. Penetration depth depends primarily upon how close the electrodes are to each other.  

Dielectric spectroscopy can be performed with FEF sensors such that only the layer closest to the electrodes is examined. By altering the geometry of the electrodes or by selecting the electrodes that are being excited and those being measured, the penetration of the field can be progressively moved further into the material. By using multiple channel devices, various channels can be used to measure dielectric properties of the material at different penetration depths (see Figure 2).  

Figure 2 (ALL FIGURES ARE COURTESY OF THE AUTHORS.) 

FEF sensors measure the material of interest from a single side. They are formed by a coplanar distribution of electrodes. Figure 2 shows the cross-sectional view of an FEF sensor. The spatial wavelength of the sensor, which is defined as the distance between centers of electrodes with the same polarity, can be varied depending on the excitation of the electrodes. The penetration depth of the sensor is mainly dependant on the spatial wavelength and is approximated as one-third the spatial wavelength. Conceptually, the penetration depth is an effective measure of the rate at which the electric field decays away from the sensor surface. Variations in penetration depth also can be obtained from the various electrode excitation patterns to provide the sensor with access to different layers of the test specimen. In Figure 2, the parameters 

 represent three different ways of connecting the electrodes. 

 represents the driving electrodes; an AC voltage signal is applied to these. 

 represents the sensing electrodes. At these electrodes, current–voltage measurements are made. 

 represents the guard electrodes. In the table below the figure, the cross-sectional view of the sensor in Figure 2, each row corresponds to one of the three connection schemes; each column corresponds to the types of connections for the electrode directly above the column used in the different connection schemes. For 

, every other electrode is driven and the undriven electrodes are treated as sense electrodes. For 

 several undriven fingers are chosen as guard electrodes and only the current–voltage at the sensing electrodes is measured. The guard electrodes are either connected to ground or kept at the same voltage potential as the sensing electrodes by using unity-gain buffer amplifiers. The spatial wavelength of the sensor is increased (λ

) by using various connection schemes. A ground plane below the electrode surface eliminates the effects of underlying material and stray noise in the medium. Major advantages of this technology include one-sided access to materials under test, profiling and imaging capabilities, and model-based signal analysis. 

 There are many cases in which gas or liquid analytes must be sensed in small concentrations close to, or below, the detection threshold of a dielectric spectroscopy measurement device. Direct sensing is simpler than preconcentration; however, direct sensing is not always possible. Preconcentration should be performed in the following cases:   

When the analyte is at a concentration at or below the measurement threshold

To select an analyte of interest from a mixture of gases or liquids.  

The advantage of direct sensing is that it is does not require the modification of the measurement setup for a particular analyte. However, because the preconcentration process acts as an amplifier, the signal-to-noise ratio may better than the direct-sensing case.  

To give a feel for typical magnitude of the effect, consider polyimide used for humidity measurements. Polyimide is hydrophilic and can absorb water for as much as 3% of its dry weight (23, 24). Because of the high dielectric permittivity of water (ε

≈ 80) in comparison with the dielectric permittivity of dry polyimide (ε

≈3), even relatively small changes in humidity may increase the dielectric permittivity of polyimide by as much as 30% (25, 26).  

Dielectric spectroscopy has been used with gaseous analytes as well. With gases, the preconcentrator material will typically form a compound with the gas and may require heating of the substrate to increase the analyte absorption. Metal compounds are used as preconcentrators for gaseous analytes (27, 28). For example, use of Ni-PC (tetra-t-butyl phthalocyaninatonickelCII) coating as a preconcentrator improves the sensitivity and signal-to-noise ratio of interdigital sensors toward ethanol fumes (27). Preconcentrators improve the selectivity of the sensor for a particular analyte because it is possible to select a preconcentrator with an affinity for only the analyte of interest. Hagen et al., provided examples of highly specific preconcentrators in which the addition of a chromium layer between interdigitated gold electrodes and a zeolite Pt-Na-Y membrane resulted in highly selective sensors for hydrocarbon vapors such as butane (29). 

 Benefits associated with preconcentrators include improved sensitivity and versatility with various analytes. However, with the use of preconcentrators, sensor response is subject to hysterisis because the sensitive layer cannot reach immediate equilibrium with the environment.  

Relaxation processes in dielectric spectroscopy are very similar to relaxation processes in the optical regime. However, the interrogation frequencies used in dielectric spectroscopy are lower than optical frequencies, so this technique studies molecular interactions such as polymer reconfiguration within a matrix, percolation processes, and moisture diffusion. A major advantage of dielectric spectroscopy is that it can be performed over a wide band of measurement frequencies. The lowest boundary for frequency in existing dielectric spectroscopy is around 1 μHz, and the highest is in THz range. It is rarely practical to go to such extremes; most practical industrial measurements are accomplished in the range from 1 Hz to 100 MHz. In the PAT context, time available for inline measurements can be a limitation for a spectroscopic approach. This can be overcome by selecting frequencies at which the response to properties of interest provides adequate sample resolution. In this section, experimental data is analyzed using a spectroscopic approach, and, after some postprocessing, a single frequency optimum for measurements is selected. 

Figure 3 (ALL FIGURES ARE COURTESY OF THE AUTHORS.) 

A spectroscopic approach involves a detailed study of various processes that occur at different frequencies. Application of an oscillating electromagnetic field to samples results in either ionic charge transport, at low frequencies, or dielectric loss caused by energy dissipation at higher frequencies. These phenomena are dependant on physical properties of the dielectrics such as viscosity, density, and moisture content. Therefore, the response of a sample over a range of excitation frequencies provides information pertaining to the physical structure and composition. Figure 3, Figure 4, and Figure 5 provide an example of spectroscopic measurement. The impedance spectra of tablets with 25, 50, 75, and 100% coatings were measured. Figure 3 shows a monotonic trend in capacitance with coating levels. In Figure 3, the smooth trend in the frequency spectrum is a characteristic of dielectric spectroscopy measurements. Unlike other technologies such as NIR, measurement peaks resulting from relaxation times are visible with wide frequency bands spanning several decades. Moreover, the separation between samples at lower frequencies is greater because at lower frequencies the system is not as dissipative. It should be noted that the discontinuity in measurement from 1 kHz to 2 kHz is a result of the switching of internal circuitry in the measurement instrumentation and does not reflect dielectric behavior. Figure 4 shows the conductance spectrum for the four samples. At lower frequencies, the conductance measurements are too low for instrumentation. At low frequencies, conductance is a result of charge transport processes. Figure 5 shows the variation in phase spectrum with various coating thicknesses.  

Figure 4 (ALL FIGURES ARE COURTESY OF THE AUTHORS.) 

All chemicals that constitute a pharmaceutical sample influence the output of the sensor. Hence it is important to isolate the effect of the analyte under study on the sensor output from those from other chemicals. This is a complex problem that requires extraction and isolation of the contributions of each constituent based on a weighted convolution of their known relaxation processes and that observed in the mixture.  

Figure 5 (ALL FIGURES ARE COURTESY OF THE AUTHORS.) 

A simple solution is to create a calibration-based estimation model that uses the spectroscopic data from a narrow frequency range and makes approximations on the type of relaxations, transport mechanism, and chemical processes being observed (30–32). This method is most commonly used in situations in which the variations in the composition are in a narrow range and speed is of more importance than accuracy. 

In a rigorous and exhaustive approach to model multianalyte concentration variations, the authors consider the material under test to be a homogeneous mixture of the analyte, other constituent chemicals, and air. By varying the proportion of this mixture, we can vary the relative concentrations of the materials. The effect of variation in proportion of the mixture on its dielectric property is nonlinear and requires considerable effort to model. The static (frequency-independent) relative dielectric permittivity for an isotropic ensemble of 

 equivalent molecules each having a dipole moment, μ, and contained in a spherical volume, 

(0) is the instantaneous dipole moment of the macroscopic sphere at the arbitrary time 

=0. The equilibrium value of the scalar product ‹

(0)› taken over all the complexions of the ensemble will be independent of time for a stationary thermodynamic system. Considering a system of 

 (0) is the instantaneous liquid dipole moment of the 

th constituent of the mixture. Detailed analysis of mixtures and their dielectric properties can be found in References 33–37. 

Obtaining wideband spectroscopic measurements can be lengthy process because the time available for inline measurements is limited. Alternatives for such measurements are single-frequency or selected-frequency excitations. In this case, the response of a sample at a single or selected frequency can be used to monitor changes in physical properties. It should be noted that the observation of dielectric spectra at just selected frequencies could result in measurement artifacts being misinterpreted as dielectric phenomena; therefore, such an approach is possible if there is some prior knowledge of the sample response at the chosen frequency. Off-line studies can establish such characteristics of the sample. To illustrate a single- and selected-frequency approach, data presented in Figure 3 is analyzed further to establish a dependency between capacitance and coating thickness. Ideally the mapping between measured capacitance, 

The first term of this equation represents the exponential decrease in the effect of the core layers as the thickness increases. The second term is a function that accounts for the escaping electronic fields. Because these are mostly owing to the zero-order field lines and their dispersion depends on the material property of the coating, this loss also can be written as a function of the thickness, χ

; α and β  are constants. To obtain thickness as a function of capacitance, the previous equation can be rearranged, expanded as a series in 

, and the higher order terms neglected. This gives a polynomial fit. Because there are only four data points and at least one degree of freedom is needed to validate the fit, a quadratic expansion is chosen. Thus, the thickness of the coating can be obtained from the measured capacitance using the equation 

 is the capacitance measured at frequency 

, respectively, are the quadratic, linear, and offset coefficients. 

Table I: Results of the ranking algorithm showing the top three frequencies ranked by mean error and standard deviation. Error values are normalized. 

To determine these coefficients and the frequency to be used, a simple algorithm is performed. First, the measured capacitances are averaged over the four samples for each of the coating levels. For each frequency, the minimum capacitance is subtracted to remove offset errors. Then the values are divided by the maximum capacitance at each frequency to remove the linear errors. The resulting capacitance ratios at each frequency are individually fitted to the quadratic equation above, and the resulting error in estimation is computed. The mean of the estimation errors and the standard deviation is tabulated for each frequency. The frequencies are ranked on the basis of least mean error and the least standard deviation. The product of these rankings is used as a figure of merit for selection of the frequency 

. Table I shows the results obtained from the algorithm detailed above. It shows that the measurements at 100 kHz are the most suitable for the estimation process. The coefficients 

 were determined by least square methods. The resulting equation for estimation of coating thickness is 

Articles

This tutorial paper is meant to aid in dielectric-sensor selection