#### Abstract

We present a methodology and a circuit to extract liquid resistance and capacitance simultaneously from the same output signal using interdigitated sensing electrodes. The principle consists in the generation of a current square wave and its application to the sensor to create a triangular output voltage which contains both the conductivity and permittivity parameters in a single periodic segment. This concept extends the Triangular Waveform Voltage (TWV) signal generation technique and is implemented by a system which consists in a closed-loop current-controlled oscillator and only requires DC power to operate. The system interface is portable and only a small number of electrical components are used to generate the expected signal. Conductivities of saline NaCl and KCl solutions, being first calibrated by commercial equipment, are characterized by a system prototype. The results show excellent linearity and prove the repeatability of the measurements. Experiments on water-glycerol mixtures validate the proposed sensing approach to measure the permittivity and the conductivity simultaneously. We discussed and identified the sources of measurement errors as circuit parasitic capacitances, switching clock feedthrough, charge injection, bandwidth, and control-current quality.

#### 1. Introduction

The quantification and qualification of liquid properties such as ionic concentration, permittivity, or pH are of major interest for the fluid identification in several domains like medicine, biology, or chemistry. In order to face this complex challenge, instrumentation solutions are constantly developed. This paper focuses on electrical measurements of liquid mixtures in order to extract the impedance modeled as an parallel dipole [1–3].

Classic sensing device used in most of these achievements is a solid-state interdigitated electrodes (IDE) sensor which is stimulated by a sinusoidal wave. The results in terms of fluid impedance are extracted out of the modulus and the phase of the response. A more effective way to stimulate the sensor is the use of a triangular signal (Triangular Waveform Voltage (TWV)) as described in [4] and applied to DNA hybridization detection [5]. In this paper, the proposed sensor interface implements the principle of a triangular wave voltage generation by a square current signal. This technique has been used successfully to sense constant phase element impedances but the extraction of the sensed data required the counting of all the elapsed periods of the triangular signal after its generation [6]. The present work uses this technique for measuring the conductivity and permittivity of liquids modeled as dipole. The first amelioration consists in isolating each voltage segment into a compact equation independent of the previous states of the oscillating wave. Therefore, the signal analysis is highly simplified and only requires sampling the data of a signal period to get all the expected parameters of a tested liquid. Secondly, our method implements a current-controlled oscillator (CCO) circuit of simple architecture so that the resulting solution requires a single DC source to operate. Autonomous systems are indeed highly advantageous to improve the compactness and the portability of sensing systems [7].

This paper first describes the IDE liquid sensor structure and defines its equivalent electric impedance; a finite element model (FEM) is used to compute the sensor electrical potential profile. The device stimulation theory is then demonstrated, resulting in a formal equation of a wave segment. The architecture of the CCO that implements the technique is chosen and fabricated using printed circuit board (PCB) technology. Solutions of NaCl and KCl are first measured to test the system performances for conductivity sensing. Mixtures of water and glycerol are finally prepared to show the ability of the system to sense both the conductivity and permittivity measurement. The last part of the paper highlights the error sources in terms of circuit parasitic capacitances, bandwidth, and quality of the current source.

#### 2. Sensitive Element Model

This first section defines the sensor equivalent model in terms of electrical components and highlights the sensed parameters. The device which is in contact with the fluid under test consists of an interdigitated electrodes (IDE) structure whose metal patterns are deposited on an insulating substrate. In order to establish the sensitive element model, we will consider the fluidic mixture as a material with a finite conductivity and a given relative permittivity.

##### 2.1. Interdigitated Sensor Structure

The structure view is illustrated in Figure 1 which depicts on the left the device top and cross-views and on the right the periodic cell structure cross-section; the symbols and represent the device electric contacts. A cell is defined as the periodic element of an IDE structure. The calculation domain of the impedance is therefore bounded on a single cell structure while supposing the number of fingers is high enough to neglect the edge effects. All the geometrical parameters are reported in Table 1.

The electric model of the structure is given in Figure 2 and its impedance by (1). It consists of a global parallel dipole, with the components and being, respectively, linked to the fluid conductivity and permittivity by the geometrical factor defined in (2) as the cell constant expressed in m^{−1} [1]. The component is the parasitic substrate capacitance. In the present context, we suppose the use of techniques and materials that allow for neglecting the interface solid-liquid effects. The electrodes are made of gold, a highly cathodic material. Due to this property, the ionic adsorption on gold is low, thus lowering the faradaic current transfer [2] and the electrical double layer (EDL) formation [1]. The fingers are also considered wide enough in comparison with the typical thickness of the EDL to allow for neglecting the impact of double layer capacitance (i.e., >1 *μ*m) [8]:

##### 2.2. Finite Element Modeling

The impedance values are computed by finite elements modeling (FEM) with the use of *COMSOL Multiphysics* software. The time-harmonic Laplace equation (3) is solved over a meshed domain bounded by a cell structure periodic element as shown in Figures 3(a) and 3(b). The global electric potential profile (Figure 3(c)) obtained after Laplace equation solution computation gives access to the electric field present between the metal fingers. The integration of the electric field over the straight boundaries and shown in Figure 3(a) gives access to both cell constants in the bulk and in the sensitive area, and (4). The bulk and fluid capacitances and the fluid resistance are given by (5). The height of the resolution domain is set to be twice the cell width to enclose 99% of the electric energy [9]. FEM has been chosen because the IDE classical analytical formulas given in [9] do not consider the finger height and make the cell constant estimation less accurate:

#### 3. Sensor Stimulation Theory

The theory exposed in this section explains how the sensitive element is stimulated. The sensing principle is based on the Triangular Waveform Voltage technique (TWV) as described in [4]. The goal of the proposed demonstration is to provide a compact mathematical expression of the output signal that contains the fluid parameters.

##### 3.1. Impedance Transient Analysis

The IDE structure model developed in the previous section is included in an electrical circuit with an additional serial capacitance as shown in Figure 4. The insertion of this capacitance is necessary to generate the triangular shape of the output signal; this cuts the DC path from the analog voltage node to the ground. The circuit loads the device under test with a square current signal and the response voltage is then probed. The total impedance loaded by the square current is given by Applying the TWV theory in the case of the impedance , a square-shaped current wave (7) of amplitude and period stimulates the impedance, with being the unit hard-step function. As shown in [6], the two first steps of the current wave have, respectively, and of amplitude to keep mean value of the voltage response at zero: The Laplace transform of the current wave is expressed by The application of the inverse Laplace transform on the product gives the time domain voltage response wave expression given by (9), introducing three shaping parameters named in Table 2 as , , and . The plots sketched in Figure 5 illustrate the time progression of the current wave and the voltage output wave for several values of the model electrical components:

**(a)**

**(b)**

**(c)**

##### 3.2. Response Wave Segmentation

The transient development gave the output voltage wave as a series expression. Each rising or falling wave segment contains all the information we are seeking about the tested impedance. A compact equation of wave segments thus needs to be calculated to extract the curve parameters , , and . The demonstration reported in the appendix derives the mathematical expression (10) of wave segments. The curve plotted in Figure 6 summarizes the parametric representation of a wave segment with the definition of a variable (11). The segment wave comprises two distinct parts, an exponential one and a linear one. Their separation time point depends on the impedance parameters values:

##### 3.3. Model Interpretation

Equation (10) is the compact mathematical form of each wave segment. Three particular cases can be observed from this expression.

###### 3.3.1. The Exponential Part of Is Predominant

This first case appears when the time constant is greater than the current wave oscillation half period (). Figure 7(a) illustrates how this first case is reached by increasing the liquid resistance : the exponential part is dominant and makes the linear part disappear. The derivative of (10) at the time is computed to get the equation of the approximating straight line, resulting in the initial slope expression of the wave segment function According to the comparison between and , the quantity can be approximated to one. The approximating straight voltage equation is thus given by (13). Replacing , , and by their respective electrical component expression, the total sensor capacitance can be extracted following the input values of and . The parasitic capacitance is calculated according to formula (5); the parasitic permittivity is fixed by the IDE structure board material:

**(a)**

**(b)**

###### 3.3.2. The Linear Part of Is Predominant

The exponential part of the wave segment becomes negligible in comparison with the linear one (). Figure 7(b) shows how the linear part is becoming predominant when the exponential part is going down. The wave segment equation is written as (14). This case occurs when the IDE sensor equivalent impedance magnitude is negligible in (6). The system response is perfectly triangular with a slope equal to the parameter and therefore inversely proportional to the additional capacitance :

###### 3.3.3. Both Parts of Are Appearing Distinctly

This case corresponds to the general form of the segmented voltage and is used to extract the fluid resistance . Focusing on the parameter, the exponential term can be approximated to zero with an error of less than 1% when the condition is satisfied. The parameter is therefore substituted by (15), with the parameter being directly linked to the fluid resistance :

#### 4. System Design

Working in the time domain is highly advantageous because of the immediate method to generate a square signal current and thus constantly charge and discharge the sensitive device. The system that we propose implements the modified TWV technique within a current controlled oscillator (CCO) so that no variable signal source is needed to stimulate the circuit, with the oscillating wave being generated by the system itself.

##### 4.1. Block Diagram

The system interface block diagram is given in Figure 8.

Two opposite current sources are alternatively enabled and disabled by a digital command which is the system feedback control signal. The swing of the analog wave is limited by the hysteresis of a Schmitt trigger. The system polarity changes each time reaches the upper or the lower limit of the hysteresis voltage.

##### 4.2. Circuit Implementation

The CCO architecture of the closed-loop system generates an oscillating signal and has been designed on CMOS advanced technologies in [10]. In the present work, as our goal is to prove an instrumentation concept and not yet to design a CMOS integrated sensor interface, we propose a similar architecture implemented with bipolar transistors to mirror the load current. The MOSFET devices act as switches to invert the circuit polarity. The analog voltage is monitored by a Schmitt NAND gate which yields the feedback digital control to change the current load polarity. A high input impedance analog buffer has been added to avoid deforming the triangular wave due to the relatively low input impedance of the probing equipment.

#### 5. Measurements Results

The experiments presented in this section are conducted to validate the detection methodology applied by the closed-loop system proposed in the previous section. The system circuit architecture is implemented on printed circuit board (PCB) with discrete components. Such a system can also be integrated on advanced CMOS technologies. As the main goal of this research work is to prove the concept of the sensing methodology, the discrete PCB implementation is the less costly choice to validate the sensing principle. The CMOS integration is considered as a future work that will improve the system performances in terms of detection resolution, power consumption, and size reduction. This section first describes the characterization setup and the PCB platform. Preliminary tests on lumped capacitance and resistance elements are performed to illustrate the evolution of the analog output voltage. Two types of experiments are presented to extract the conductivity and permittivity of liquids. The first type of experiment characterizes the conductivity variation of NaCl and KCl solutions while the second test highlights the ability of the system to detect both the permittivity and the conductivity of glycerol-water mixtures.

##### 5.1. Characterization Setup

The circuit architecture of Figure 9 has been implemented on a standard PCB platform (Multi Circuit Boards Ltd.) using discrete components listed in Table 3. The board includes two output connectors used to probe the analog voltage and the digital control signal. The analog buffer and the Schmitt trigger are supplied by the voltage values listed in Table 4. The trigger hysteresis amount is fixed by the supply voltage amplitude. The serial capacitance value of the system is set to 90 nF to fix an output frequency range starting from 100 Hz to 10 kHz respecting the load current amplitude range of 10 *μ*A to 1 mA. The interdigitated structures are directly printed on the circuit board.

##### 5.2. Lumped Elements Values Extraction

The analog voltage measurements have been performed using lumped components of high precision (<0.1%) that represent the IDE sensor electrical model. The tested components are listed in Table 5; the analog voltage is probed with the *Agilent Infiniium MSO8104A* oscilloscope and plotted in Figure 10. These results validate the analog output wave shape variations as predicted by the analytical model of Section 3. As the system works in closed loop, the analog wave oscillates between the two defined hysteresis levels. The values of the lumped elements have been extracted according to the compact expression of the analog voltage (15). The measured values given in Table 5 are compared to the calibrated values of the lumped components. We see that the largest error is observed for the capacitance extraction. The next section will explain in more detail the origin of this error.

**(a)**

**(b)**

##### 5.3. Liquid Conductivity Measurements

The objective of the electrolyte conductivity tests is to prove the efficiency and robustness of the present sensing system. Several NaCl and KCl solutions have been prepared and calibrated using an industrial electrical conductivity meter (*EUTECH CON 270*, 1% full scale accuracy). The IDE sensor dimensions reported in Table 6 are chosen according to the theoretical prediction of the solutions conductivities using formula (16) in order to set a total range of liquid resistance measurable by the system, with being the saline concentration (mol^{−1}) and the molar conductivities being given in Table 7. The conductivities are then extracted using expression (5); the cell constant links the conductivity to the measured liquid resistance. The value of the cell constant gives also the IDE capacitance which is about 50 pF for water with relative permittivity of 80. Since this last value is out of the PCB system capacitance detection range, an additional capacitance of 5 nF has been added in parallel with the IDE sensor:
The results plotted in Figures 11 and 12 show the progression of the liquid conductance extracted from the system triangular output voltage. The measurements have been repeated ten times. The -axis scale is the electrolyte salt concentration measured by the industrial electrical conductivity meter and expressed in mol·L^{−1} according to (16). We observe, at first sight, an excellent linearity and repeatability of the liquid measurements. The standard deviation of the measurements increases with the electrolyte concentration; a high electrolyte concentration means that the liquid resistance is small and therefore more difficult to extract. In addition, the results are compared to the expected values (continuous curves) predicted by the theoretical estimation of the cell constant; (5). The theoretical predictions and the measurements are in excellent agreement.

##### 5.4. Liquid Permittivity Measurements

The characterization of liquid permittivity is performed by mixing two liquids of different permittivities. The two liquids have to be miscible in each other to be homogeneous and thus to create the most linear permittivity range of mixtures. We start from pure glycerol in which we add water to prepare several mixtures of different water-glycerol ratio. Knowing that the relative permittivities are about 41 for glycerol and 80 for water, the IDE sensor cell constant has to be chosen to have a significant capacitance to detect. The IDE parameters used to achieve permittivity measurements are listed in Table 8; those dimensions are chosen to increase the nominal capacitance of the sensing element; a too small liquid capacitance involves a bad output voltage parameter extraction as explained further in Section 6. The measurements depicted in Figure 13 show the evolution of the dielectric constant versus the water-glycerol ratio of the mixture. As these two liquids have a different natural conductivity, the evolution of conductivity is also shown for the mixture.

**(a)**

**(b)**

The permittivity curve versus the water-glycerol ratio was supposed to be linear due to the homogeneous prepared mixtures. The measurements shown in Figure 13 confirm this linearity. The error observed on the last point of the measured conductivity is due to the liquid resistance extraction error; the smaller the liquid resistance, the higher the impact of the extraction error.

#### 6. Discussion

In every system implementation, the circuit elements are nonideal and limit the system performances. Our major concern is to insure the good quality of the analog output signal to detect capacitances and resistances on the largest measurement range; it means that the probed wave has to fit as most as possible the modeling equation (9). Therefore all parasitic impacts on the analog output voltage have to be discussed. As the system is implemented on a PCB platform by discrete components, it will be easier to highlight the system limitations. Having a good identification of these problems is important for the system implementation on any technology.

##### 6.1. Parasitic Capacitance Effects

The most significant limitation on circuit implementations is the parasitic capacitances of the discrete components. These capacitances considerably affect the analog measurement and make the sensor capacitance more difficult to quantify, lowering the sensitivity. The parasitic effects are represented in Figure 14 with all the capacitances connected to the analog voltage node. The Schmitt trigger and the analog buffer have their own input capacitances noted ; the capacitances , , and are the MOSFET switches internal capacitances. Depending on where these capacitances are connected, their impact is different on the analog voltage.

The capacitances are connected to DC nodes in parallel with the sensing element; they consequently impact only the system sensitivity. The drain to source capacitances and is very small and presents a negligible impact on the measurements. The capacitance is connected between the digital command terminal and the analog voltage. Its impact deforms the analog voltage by clock feedthrough and charge injection each time the system switches [11]; the analog output response presents a voltage drop and small peaking effects compensated by oscillations which disappear after a few microseconds. The effects of the clock feedthrough and charge injection are depicted in Figure 15. The major consequence of these phenomena is bad quantification of the capacitance during the curve parameter estimation. We observed a voltage drop of 90 mV for a tested capacitance of 200 pF. Considering the analog oscillation amplitude of 1 V (Figure 10), the liquid capacitance is extracted with error of up to 9%. The best manner to attenuate these two parasitic effects is to integrate the circuit interface on CMOS technology. The most impacting devices that increase this effect are the CMOS switches; their surface can be drastically downscaled on integrated technology. This will reduce the coupling capacitance . These effects can also be lowered by increasing the sensor nominal capacitance (i.e., increasing the finger number) to make the element small enough in comparison with the capacitance.

##### 6.2. Circuit Bandwidth

Concerning the system frequency bandwidth, the limiting factors are the Schmitt trigger switching times and the second-order oscillations caused by the parasitic capacitances of the switches. The PCB implementation of the interface includes a 4093 NAND gate component as Schmitt trigger. Its transition time depends on the supply voltage. For the tests conditions presented in Section 5, both the rising and falling times of the trigger are of 600 ns. In comparison with the duration of the second-order oscillations which exceeds 1 *μ*s, the bandwidth limitation is therefore due to the switches parasitic capacitances. For the proposed PCB implementation, we set the working frequency limit under 50 kHz; the quality of the analog output voltage is considered to be too much deteriorated beyond this limit.

##### 6.3. Current Mirrors

The last impact on the measurements is due to the current copy performed by the bipolar transistors. Regarding the datasheet of the transistors (*BCV6X*) the minimum value of the forward current gain is 100; this means that the base current error is lower than 1%. The quality of the current copy is also altered by the Early effect and the mismatch between the transistors. We compensated for this error by performing a calibration which consists in measuring the mirrored amount of current injected in a reference test capacitor. In a future CMOS design, the quality of the current copy will be improved since the gate current of MOS devices is negligible in comparison with the base current of bipolar devices. In addition, the silicon integration of the circuit brings a much better matching of the mirror than discrete implementation. Nevertheless, it is also important to note that the Early effect is worse using CMOS transistors. Fortunately, this effect can be simply corrected by cascading the current mirror.

#### 7. Conclusions

A modified TWV detection technique has been demonstrated to sense the impedance model of liquids. The theoretical result of wave segmentation includes in one compact and simple expression all the information that has to be extracted from an parallel equivalent sensor like IDE structures. The implementation of a closed-loop self-oscillating circuit interface makes the system autonomous and compact and requires only a few elements to be manufactured. Measurements on saline concentration have been measured with high accuracy in the 10^{−6} to 10^{−5} mol^{−1} input range and water-glycerol mixtures validated the sensing technique in respect to the correct range of permittivity (i.e., from 40 for pure glycerol to 80 for pure water). The results discussions highlighted the error sources in terms of circuit bandwidth, switching charge injection due to the parasitic capacitances, and quality of the current copy in the sensing device. The CMOS integration of the whole system would be the best solution to considerably attenuate the parasitic effects of the devices, to extend the circuit bandwidth, and also to reduce the system power consumption.

#### Appendix

#### Calculation of the Response Wave Segments

*Response Wave Segmentation*. Assuming that the output wave (9) is probed on a rising segment after half periods, all the hard-step functions for are elapsed and equal to 1:
A time translation is then applied following the variable change (A.2):
The expression (A.1) can be rewritten as
If we consider the square current signal as an infinite sum of hard-step functions, the expression (A.3) can be computed for . The two exponential terms of respective coefficients and thus disappear; a last serial progression has to be computed in the resulting equation
The last term of (A.4) corresponds to the well-known geometrical progression (A.5) where is the sum value, is a scale factor, and is the common ratio. The progression converges only if the common ratio modulus is smaller than unity. This is always verified in the present case because the parameters and are always greater than zero:
The final segmented expression of the output voltage wave is thus given by expression (A.6); the equation has been vertically shifted to set the starting voltage at zero:

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

The authors would like to thank the Wallonia Region of Belgium for funding this research work via the *Skywin HM+* project. They also thank the staff of the WELCOME characterization platform that allowed them to perform the measurements.