1. Introduction
The QCM sensor has been used as a label-free biosensor in recent decades, after its ability to operate in a liquid medium [
1,
2,
3] was demonstrated. The QCM sensor has been adopted as a biosensor because it is versatile and durable, involving both a simple and very stable structure. The QCM sensor has numerous nanoscale applications, including cells [
4], antibody interactions [
5], detection of DNA hybridizations [
6], and viruses [
7]. However, the application of piezoelectric materials as biosensors is a challenge due to losses induced by contact with the liquid medium. The detection limit for the QCM sensor is lower in comparison with surface acoustic wave (SAW) sensors [
8]. Furthermore, with regard to label-free biosensing in liquids, surface plasmon resonance (SPR) spectroscopy [
9] has a limit of detection lower than the QCM sensor. On the other hand, the QCM sensor provides access to the physical parameters of the sample by measuring the dissipation factor or another equivalent electrical parameter, ensuring a detailed analysis of the surface and interactions on the surface without being limited to mass measurement per unit area. Based on the dissipation parameter, the viscoelastic and conformational properties of the sample [
10] are monitored. The most powerful methods based on passive interrogation are impedance analysis (QCMI) [
11,
12] and the ringdown method (QCM-D) [
13,
14,
15].
The methods regularly found in the literature focus exclusively on improving the accuracy of resonant frequency and the dissipation factor measurement. This is because resonance frequency and the dissipation factor are very valuable parameters for evaluating the behavior of the QCM sensor in a liquid medium. It is important to note at this time that the measurement of all parameters of the Butterworth–van Dyke (BVD) model ensures a clear picture of the interactions at the surface of the QCM sensor, which is always harshly dependent on experimental circumstances. The use of the QCM sensor in experimental circumstances characterized by temperature and flow gradients is subject to uncertainty.
Eliminating uncertainty implies an ideal configuration of the experimental setup for the QCM sensor. An ideal configuration, from which both internal and external factors that disturb the stability of the QCM sensor are eliminated, cannot be achieved experimentally because it is impossible to make portable equipment for large-scale applications of the QCM sensor. Temperature changes [
16] mostly affect the behavior of the QCM sensor, and active temperature control can in turn become a noise source. If the liquid medium is biological, then the working temperature is imposed by it, and its stabilization usually involves significant experimental resources. The measuring cell, regardless of its topology, exerts mechanical stress on the QCM sensor, having the effect of reducing its quality factor (
) and finally resulting in the degradation of the frequency stability. Despite all the advances made in the field of electronic devices, their thermal noise negatively influences the stability of the QCM sensor, especially in active measurement configurations. The jitter of the periodic signal used as a clock reference in modern passive measurement instruments is also a disturbing parameter that has the effect of reducing the resolution of the QCM sensor. If the QCM sensor is used in a liquid environment, then the pressure [
17] or flow [
18,
19] differences induced by the pumping system will also have a negative influence on the stability.
The use of external sensors to monitor environmental influences is possible, but the evaluation of real-time corrections is difficult due to the propagation times of these influences in the experimental setup being simultaneous with the changes that take place on the surface of the QCM sensor. An instrument based on impedance analysis that can monitor in real time all the parameters of the BVD model would eliminate a misinterpretation of the external events that occurred during the measurement. In this way, even if the external influences are not eliminated, their identification would always be possible by measuring the impedance of the QCM sensor. With each decade, the cost of impedance analyzers has decreased [
20], so the major disadvantage related to their cost has now been eliminated.
Based on the recent concept in instrumentation, a fast and accurate virtual impedance analyzer (VIA) was developed [
21]. Today we can consider that a modern VIA must ensure the real-time calculation of the equivalent electrical parameters for the QCM sensor. The concept of a virtual instrument implies the existence of a field programmable gate array circuit (FPGA) in which configurable hardware modules are implemented and analog support is provided by an analog-to-digital converter (ADC) and digital-to-analog converters (DACs).
This paper aims to demonstrate the benefits of advanced impedance spectroscopy relative to the Butterworth–van Dyke (BVD) model for QCM sensors immersed in a liquid medium. The VIA defined by software with real-time computation capabilities for QCM sensor parameters covers a wide range of measurements for the impedance. Advanced software methods, such as self-compensation, post-compensation, and simultaneous computation by using a multimethod based on the BVD model, are the source of experimental and theoretical results presented in this paper. In the literature, few advanced electrical models were proposed [
22,
23,
24,
25,
26] to solve the limitations of the BVD model in a liquid medium. Some electrical models [
27], inspired by the field of impedance spectroscopy, are not related to the BVD model or based on mechanical or acoustic analogy.
This paper makes the following contributions: (i) advanced QCM sensor characterization based on BVD model in a liquid medium; (ii) valuable consideration of the BVD model in accord with the energetic transfer model; (iii) an innovative post-compensation strategy based on the Nyquist plot; and (iv) a theoretical and experimental study of the capabilities of the BVD model to describe the QCM sensor immersed in a liquid medium. Experimental validation of proposed new post-compensation strategies based on the Nyquist plot was performed in a liquid medium with different viscosities. In this way, the advantages of compensation software techniques for the realization of high-performance applications in the area of QCM sensors are confirmed without resorting to dedicated electronics analog circuits or specialized hardware resources.
2. Materials and Methods
2.1. Virtual Impedance Analyzer
Traditionally, a simplified “active” method, in which the QCM sensor is part of an oscillator circuit, is frequently used. The amplifier contained in this circuit has a certain influence on the oscillation frequency. For this reason, “passive” methods are more accurate, and even if they are more sophisticated, they are preferred. Of these, the impedance analysis method based on passive interrogation is commonly used in applications [
20] to investigate the behavior of the QCM sensor.
Figure 1a shows the basic configuration for passively measuring the impedance of a QCM sensor in air. This configuration can provide a wide impedance measurement range. The configuration shown is of the voltage divider or half-bridge type with a passive interrogation in the resonant frequencies range of the QCM sensor. The high impedance analog inputs ensure the measurement, in digital synchronism conditions, of the voltages
and
. The impedance of the QCM sensor
is calculated, taking into account a reference resistance
. In this case, the impedance of the QCM sensor, in terms of voltage, is given by the relation
The configuration shown in
Figure 1a is recommended for measuring the impedance of the QCM sensor with a high-quality factor (
). This situation is encountered in the air or gaseous medium when the QCM sensor is usually functionalized by using a thin layer. The possibility of being able to amplify (or by using another voltage scale) the voltage on the reference resistance against zero background transforms this configuration into the basic method in impedance measurements.
A QCM sensor in the air has a high parallel resonance impedance and must be wired in the configuration shown in
Figure 1a. It is recommended that the electrode in contact with the liquid medium be connected to the ground to provide electrochemistry applications. If grounding of the liquid-exposed electrode is important, a balun transformer and appropriate compensations can be used as shown in
Figure 1b. This configuration is mandatory for electrochemistry and recommended in the case of a liquid medium in which its dielectric properties can disturb the resonance of the QCM sensor.
The most important criteria for a virtual instrument, used in monitoring a QCM sensor’s impedance, are the lowest possible power consumption and the smallest possible physical size. The extensive presentation of VIA used in the experimental setup and the standard compensation concepts [
21,
28,
29,
30] are available in the literature.
Analog Discovery 2 (AD2) from Digilent Inc. (Pullman, WA, USA) [
31,
32,
33,
34] is the selected virtual instrument for performing a VIA. The choice of AD2 was determined by the performance of the analog input–output channels. The analog channels have two ADC converters and two DACs, each with a resolution of 14 bits at 100 MSPS. The hardware modules required to make the virtual instrument have been implemented in a Xilinx Spartan 6 FPGA (XC6SLX16-1L). The user has two high-impedance input channels (1 MΩ with 24 pF in parallel) with voltage ranges up to ±25 V. For a voltage range ≤ ±5 V, the absolute resolution is 0.32 mV. The reading of the analog input channels is digitally synchronized, which allows the accurate determination of the phase difference between the input signals. The two analog outputs are typically used by the arbitrary waveform generator (AWG) which provides output signals in the range of ±5 V.
The hardware interface with the QCM sensor is provided by the shield for the impedance analyzer application [
28,
29], transforming AD2 into two versions of VIA, as shown in
Figure 2. In both versions, an analog-output channel provides the sinusoidal signal of passive interrogation. The AWG interrogation sine wave is applied to the QCM sensor, and its response is measured in digital synchronization conditions through the two analog inputs. The measurement takes place directly on the QCM sensor or with a balun transformer.
The selection of the reference resistance is done with relays, controlled by the digital outputs of AD2. Both topologies shown in
Figure 1 have been implemented, as can be seen in
Figure 2. The experimental support thus created allows a comparative study about the performance of the advanced software-compensation methods.
For the topology shown in
Figure 2b, the compensation is implicit and without a standard compensation. It is practically impossible to measure the electrical parameters of the BVD model. The experimental setup is completed by the QCM flow cell kit (011121, ALS Co., Ltd., Tokyo, Japan), mounted in its static measurement mode [
35].
The PC-type host computer provides the configuration commands, data transfer, and power supply through the USB2 interface of the AD2 virtual instrument. Data acquisition and processing are performed by a Python module that exploits the software development kit (SDK) functions of AD2. The Python module provides real-time processing of experimental raw data and calculates the raw parameters of the BVD model simultaneously by several methods. The acquired raw experimental data are graphically represented together with the parameters of the BVD model. To ensure advanced post-processing in the MATLAB® environment, the raw experimental data are recorded in * .csv files.
2.2. The QCM Sensor and BVD Model
Based on the energy transfer model, the QCM sensor generates and stores acoustic energy. For this reason, we can consider an analogy with an electrical model consisting of an
series circuit to describe its behavior. If the QCM sensor is transformed into a functional sensor, it is covered with material, e.g., a sensitive film. A small part of the acoustic energy is transferred to this sensitive film, which has partly elastic and partly viscoelastic properties. This energy is stored in its elastic part and dissipated in its viscoelastic part. In an ideal approximation of an infinite-quality Q factor of the QCM sensor, the applied electrical power is equal to the sound power at the resonant frequency. Taking this into account, the acoustical impedance is analogous to the electrical impedance. Furthermore, in this process the liquid medium, in which the QCM sensor is used, intervenes. It interacts with the respective materials that are attached to the sensitive film. For the attached materials we must perform an identical consideration about their elastic and viscoelastic properties. Moreover, the dielectric properties of these materials can induce unexpected behavior in the QCM sensor that is reflected in the modification of the electrical parameters. The BVD-lumped electrical model [
11] describes the modes of interest of the QCM sensor by using 4 parameters, as shown in
Figure 3.
The arm consisting of the series combination of a resistance
, inductance
, and a capacitance
is responsible for modeling the storage and dissipation of acoustic energy. Parallel to this arm is a shunt capacitance determined by the existence of electrodes on the two sides of the quartz crystal. The
Z(
ω) impedance of the QCM sensor based on the BVD model is determined by the series
and parallel
impedance combination of the two arms:
The reactance of the shunt capacitance
is given by
and the impedance of the series arm is given by the following equation
where
,
, and
are the parameter of the series arm of the BVD model.
The BVD model for the QCM sensor with a single electrode in contact with the liquid medium is composed of the impedance of the load
in the series arm and a stray capacitance
in the parallel arm. The experimental measured raw data for the QCM sensor in the liquid medium is expected to fit with BVD model parameters calculated considering the following equation:
,
, and
, where
,
, and
are the parameter calculated from experimental measurements
The constraint for the motional capacitance to be constant and equal to its calculated value for the undisturbed state is established by the mechanical or acoustic analogy. The equivalent BVD model derived from theoretical approximations fits the experimental data for the QCM sensor used in the liquid medium. Equation (5) is the key objective of all active or passive methods of measuring the parameters of the QCM sensor. To achieve this objective, many analog topologies [
3] and many calibrations and compensation methods have been proposed [
36].
In the literature, more models are presented, not relatable by a mechanical or acoustic analogy, inspired by the field of impedance spectroscopy [
27]. The analogy of the electrical model to the acoustic load impedance generated by the interaction of the QCM sensor with the liquid medium is mandatory. Any proposed electrical model must be derived from a mechanical or acoustic analogy to be considered as a realistic approach. Experimental investigations have shown a complex behavior for the QCM sensor, such as increasing the resonance frequency [
37] under a silicon oil load. In general, for piezoelectric material coupled with other materials, we observe a severe modification in equivalent electrical parameters.
To validate the software advanced impedance analysis based on the BVD model, very precise experimental data obtained by scanning with a small step is used. Another key parameter strongly affected for QCM sensor in liquid medium is the quality factor,
. The quality factor is related to the half-power spectrum at the resonance frequency. This is also known in terms of the network analyzer as the −3 dB points because the response is down from the maximum by three decibels. In the ringdown measured method [
13,
14,
15], the dissipation factor (
) is preferred to describe the natural acoustic energy losses at the resonance, and it is very simply related to the quality factor by
.
2.3. Computation of the QCM Sensor Parameters from Raw Experimental Data
2.3.1. Based on the Experimental Quality Factor
The value of the
series resistance and the
series frequency is determined after considering the sequence of the raw experimental data of the impedance measurements, in the frequency range of the resonant frequencies of the QCM sensor with the help of the minimum search function. These are key parameters of the BVD model, and many instruments are limited to measuring them. The bandwidth
is calculated from the intersection of the raw experimental data with a line parallel to the abscissa at
. The anti-resonance frequency
is determined by a maximum search function in the sequence of raw experimental data. Based on the quality factor equation
QCM sensor parameters that do not result directly from the raw experimental data can be calculated based on the following equations:
The accuracy of the BVD parameters calculated from the measured raw data depends on the frequency-swiping step as well as on the calibration and compensation procedures. This computation method shown above cannot be applied only after the end of the acquisition process. In other words, the method based on quality factor cannot be used for real-time compensation of the shunt capacitance .
2.3.2. Based on the Experimental Resonance Frequencies
A very efficient method is based on the direct measurement from the impedance raw experimental data of the series resonance frequency
and the anti-resonance frequency
with the help of minimum and maximum search functions. The determination of the series resistance
is implicit in this case. A measurement of the impedance of the QCM sensor at the frequency
, less than or higher away from the resonant frequencies where its reactance
is purely capacitive, allows the calculation of the shunt and stray capacitance
based on the equation
The calculation of the parameters that cannot be determined directly from the raw experimental data is done based on the following equations:
and
Because the shunt and stray capacitance is known at the beginning of the process of measuring the impedance of the QCM sensor, in the area of resonance frequencies, a real-time compensation of the sensor can be obtained.
2.4. Standard Method for Impedance Compensation
Modern advanced instruments with data-processing capabilities ensure automatic calibration or self-calibration procedures. This is undoubtedly the main advantage of digital technology for the production of high-precision instruments with extensive data-communication facilities through high-speed interfaces. An advanced impedance analyzer measures the electrical impedance of the QCM sensor around the resonant frequencies. The measurement process is followed by a complete characterization of its response by calculating the electrical parameters of the BVD model. Moreover, if the impedance analyzer is used to identify the parameters of the BVD model for the QCM sensor in a liquid medium, the compensation procedures are of the utmost importance. A minimum and necessary standard compensation is described below.
The standard compensation procedure [
30] is shown in
Figure 4, which involves the following processes: (i) open-circuit parasitic impedance compensation,
(
Figure 4a) followed by (ii) short-circuit parasitic impedance compensation,
(
Figure 4b). After
and
have been measured, the actual impedance of the QCM sensor is calculated by using the following equation:
where
is the QCM sensor raw measured impedance.
2.5. Real-Time Shunt and Stray Capacitance Compensation
Real-time analog compensation of shunt and stray capacitance is a common topic in the literature [
3], and active or passive circuits are used for this purpose. The digital method proposed here is based on Equation (2), and the method of calculating the parameters of the BVD model is presented in
Section 2.3.2. Equation (2) can be rewritten in the following form:
where the
is the measured raw impedance. From Equation (14) we can calculate the raw experimental data of the series arm (motional arm):
where
computed in every measurement iteration based on previous measurements of the
as was presented in Equation (10).
The most important advancement of the VIA is related to the capabilities of real-time data manipulation, thus ensuring the functions of compensation through software virtualization and the component of front-end electronics is kept to a minimum. The advantages of software-based instruments in order to customize the hardware of virtual instruments of general use are not only reduced to an economy of electronic devices through the dematerialization of applications, but also through advanced mathematical methods that generate a new level of performance. It is extremely beneficial to replace electronic devices or complex analog circuits with a few program lines. In addition, thermal noise or other parasitic effects specific to analog compensation methods are eliminated.
2.6. Advanced Residual Capacitance Compensation in VIA
A natural question about modern virtual instruments based on software and intensive data processing is how far they can go with compensation and self-calibration methods. The only true answer is that technologically we can go up to the physical limit of the investigated phenomena. This situation becomes practically evident when virtual instruments capable of performing complex calculations of a large volume of acquired experimental data are involved. If these calculations are performed in real time, then the measurement and interpretation of the results of measurement becomes a single operation with huge benefits in practice.
Commonly in the literature [
3] related to the QCM sensor, QCM sensor frequency response is shown as a Bode plot. This representation of the frequency response of the QCM sensor is very useful in practice but at the same time limited in being able to depict the complexity of the response. The Nyquist Plot or complex-impedance plane representation, in which the data from each frequency point is plotted by the imaginary part on the ordinate and the real part on the abscissa, is very sensitive to changes. It is a common convention in the electrochemistry community to plot
on the y-axis. After the standard compensation procedure described in
Section 2.5, the residual capacitance is generated by the imprecision of shunt and stray capacitance measurement. The simulated effect of the residual shunt and stray capacitance is shown in
Figure 5.
The residual shunt and stray capacitance
are defined by the following equation:
The simulation from
Figure 5 is based on the BVD model considering the values of the parameters typical for a QCM sensor with 10 MHz serial resonant frequency and one electrode immersed in a liquid medium. This residual shunt and stray capacitance are inherently due to the uncertainty of the
measurement and perturbations from the experimental environment during raw data acquisition. The Nyquist plot based on the BVD model is made to identify the contribution of the residual capacitance. For
, the diameter of the circle in the Nyquist plot is infinite for measurements in an infinite frequency range. In practice, the impedance of the QCM sensor is measured in a narrow range around resonant frequencies. In this situation, the condition for perfect compensation for the shunt and stray capacitance of the QCM sensor is
or in other words, we must have a straight line parallel to the imaginary axis in the Nyquist plot. To minimize residual capacity, an algorithm of successive approximations is used. Because the experimental raw data have a statistical dispersion and are subject to uncertainty, the residual capacity cannot be reduced to zero. In an extremely short post-processing time, the real-time compensated experimental raw data for shunt and stray capacitance is refined to reduce the residual capacitance to less than
.
In this new approach specific to virtual instruments, we can make very precise measurements based both on the physical qualities of the electronic devices used and on the computing performance involved in the process. This beneficial change from the perspective of the costs of achieving the experimental setup ensures a new level of performance inaccessible to traditional methods. The self-calibration procedures followed by compensation or adaptive compensation are applicable in different stages: (i) before the data acquisition, (ii) during the data acquisition, or (iii) after the data acquisition based on algorithms usually implemented in the host computer. The ultimate goal of these calibration, compensation, and post-compensation procedures is to obtain a set of experimental data to satisfy Equation (5), the dream of many analog methods. This approach ensures better efficiency of the measurement process than the involvement of methods based on nonlinear fits. The nonlinear fit involves a long calculation time and the final result depends on the initial conditions which in turn are extracted from the raw experimental data.
In conclusion, the Nyquist plot is valuable for identifying how many characteristic features are exhibited by the QCM sensor, and all frequency information is inherently lost. Even if the Nyquist plot is a representation of the same experimental data, or the result of a simulation based on the BVD model, the depiction of the imaginary and real component using two linear axes changes the perspective in impedance analyses.
4. Discussion
Electrical systems, mechanical systems, and acoustical systems are three physically different systems. Frequently, it is sensible to describe some physical systems in terms of electrical analogies. This analogic approach can be considered the birth certificate for many innovative scientific fields such as mechatronics or piezotronics, to name just two. The different physical systems are analogous to each other if the differential equation modeling of these systems is identical. The direct analogy between resistance in the BVD model and acoustic energy dissipation (acoustic resistance), for the QCM sensor, is very useful in the right interpretation of the experimental results. The frequency shift of the resonant frequency
is correlated with a mass attached to the sensor surface. Both acoustic circuit elements and electrical equivalent elements in a complex conjugate form are responsible for storing energy and are accountable for the resonance frequency. The real source of frequency shift is a modification of the acoustic load impedance, and in this situation, the best fitting is ensured by the best equivalent electrical impedance composed by resistance, inductance, and capacitance. The frequency shift is related to a modification in the elements of the system in charge of energy storage. Acoustical resistance, or in analogy, electrical resistance, reflects all dissipation processes from the QCM sensor, functional layers, and its contact with other materials in a liquid medium. The BVD model is capable of correctly describing the modification in electrical impedance for the QCM sensor around resonant frequencies in a high-viscosity liquid medium after a full compensation of shunt and stray capacitance. The theoretical approach based on the Nyquist plot was experimentally validated in the very harsh medium for the QCM sensor, pure glycerin.