Next Article in Journal
Multi-Class Classification and Multi-Output Regression of Three-Dimensional Objects Using Artificial Intelligence Applied to Digital Holographic Information
Next Article in Special Issue
The Luminescence Hypothesis of Olfaction
Previous Article in Journal
Selective Voltammetric Sensor for the Simultaneous Quantification of Tartrazine and Brilliant Blue FCF
Previous Article in Special Issue
Development of a Low-Cost Electronic Nose with an Open Sensor Chamber: Application to Detection of Ciboria batschiana
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Gas Adsorption Response of Piezoelectrically Driven Microcantilever Beam Gas Sensors: Analytical, Numerical, and Experimental Characterizations

1
SDU NanoSYD, Mads Clausen Institute, University of Southern Denmark, Alsion 2, 6400 Sønderborg, Denmark
2
SDU Mechatronics, Department of Mechanical and Electrical Engineering, University of Southern Denmark, Alsion 2, 6400 Sønderborg, Denmark
3
AmiNIC ApS, Jernbanegade 75, 5500 Middlefart, Denmark
4
Fraunhofer Institute for Silicon Technology, Fraunhoferstraße 1, 25524 Itzehoe, Germany
*
Author to whom correspondence should be addressed.
Sensors 2023, 23(3), 1093; https://doi.org/10.3390/s23031093
Submission received: 25 November 2022 / Revised: 11 January 2023 / Accepted: 12 January 2023 / Published: 17 January 2023
(This article belongs to the Special Issue Recent Advancements in Olfaction and Electronic Nose)

Abstract

:
This work presents an approach for the estimation of the adsorbed mass of 1,5-diaminopentane (cadaverine) on a functionalized piezoelectrically driven microcantilever (PD-MC) sensor, using a polynomial developed from the characterization of the resonance frequency response to the known added mass. This work supplements the previous studies we carried out on the development of an electronic nose for the measurement of cadaverine in meat and fish, as a determinant of its freshness. An analytical transverse vibration analysis of a chosen microcantilever beam with given dimensions and desired resonance frequency (>10 kHz) was conducted. Since the beam is considered stepped with both geometrical and material non-uniformity, a modal solution for stepped beams, extendable to clamped-free beams of any shape and structure, is derived and used for free and forced vibration analyses of the beam. The forced vibration analysis is then used for transformation to an equivalent electrical model, to address the fact that the microcantilever is both electronically actuated and read. An analytical resonance frequency response to the mass added is obtained by adding simulated masses to the free end of the beam. Experimental verification of the resonance frequency response is carried out, by applying known masses to the microcantilever while measuring the resonance frequency response using an impedance analyzer. The obtained response is then transformed into a resonance frequency to the added mass response polynomial using a polynomial fit. The resulting polynomial is then verified for performance using different masses of cantilever functionalization solution. The functionalized cantilever is then exposed to different concentrations of cadaverine while measuring the resonance frequency and mass of cadaverine adsorbed estimated using the previously obtained polynomial. The result is that there is the possibility of using this approach to estimate the mass of cadaverine gas adsorbed on a functionalized microcantilever, but the effectiveness of this approach is highly dependent on the known masses used for the development of the response polynomial model.

1. Introduction

Piezoelectrically driven microcantilevers (PD-MCs) have been widely used in gas sensing applications through the adsorption of target gas molecules on functionalized surfaces [1,2,3]. The adsorption of target gas molecules causes a change in the surface stress of the functionalized microcantilever, resulting in a change in resonance frequency and deflection [4,5,6,7]. We previously studied the application of PD-MCs in the electronic nose for the detection of 1,5-diaminopentane(Cadaverine), as a measure of the extent of meat spoilage [8,9]. Cadaverine is a biogenic amine formed by the decarboxylation of amino acids by enzymes of bacterial origin [10], whose concentration has been observed to increase during the spoilage of meat [11]. Therefore, the presence of few target gas molecules is reflected by changes in the resonance frequencies being measured and impedance output actuated through piezoelectric material [12].
Sensing and detection are typically performed using two approaches, the dynamic approach and the static approach. The dynamic approach relies on the analysis of change in resonance frequency while the static approach uses the change in deflection of the microcantilever. In both cases, the change in property measured is proportional to the mass of adsorbed [13,14,15]. The dynamic approach has recently been favored over the static approach due to varying complexities in methods used in the measurement of the deflection of microcantilever sensors. Optical deflection measurement is the most employed method for static approach and has been and has been reported to come with complications in readout and complexity in set-up and alignment, affecting the sensitivity [16,17]. The dynamic approach, on the other hand, offers other readout parameters besides changes in resonance frequency. These include the piezoresistive readout, where the change in resistance of the microcantilever is measured using a Wheatstone bridge with three reference resistors, of which one is adjustable [18], the piezoelectric readout where the piezoelectric microcantilever is driven with the inverse piezoelectric effect (self-excitation) by applying an electric AC voltage to the piezoelectric material, and sensing microcantilever bending by recording the piezoelectric current change [4,19]. Previous studies on gas-sensing applications using PD-MCs have employed the dynamic approach, where the microcantilever is actuated by applying a sinusoidal voltage at a frequency near its natural frequency, to a piezoelectric layer located at the clamped end [20,21,22].
While the PD-MC-based sensors are considered to perform generally well, they have challenges that affect the quality of the measurements. The piezoelectric material layer provides non-linear responses originating from non-linear elastic vibrations [23], and forced vibrations [24]. Other notable sources of nonlinearity and noise in the response signal include flexural rigidity of the beam [25] and adsorption surface-induced stresses [4,26], changes in temperature causing a thermal drift in the output signal, noise in the driving voltage, and surface stress changes originating from the functionalization surface-coating modification due to target molecular binding and unbinding [27,28,29,30,31]. It has been vastly documented that the resonance frequency of a microcantilever beam varies with mass added by using a relationship that can be summarized as [32,33,34]:
f m = 1 2 π k m + δ m
where f m is the resonance frequency, k is the lumped stiffness of the microcantilever system, m is the lumped mass, and δ m is the added mass. The major drawback of this relationship is that it assumes constant stiffness k, neglecting the effects of changes to surface stress and subsequently k.
Various analytical studies on mass sensing using dynamically actuated PD-MCs have been conducted based on the relationship in Equation (1). Bouchala et al. (2017) analytically investigated the linear dynamic response of beams to added mass modeled as discrete masses. The study found that the nonlinearity in the response signal was introduced by the electrical excitation when the microcantilever was dynamically excited [35]. Lam et al. (2022) further investigated the effect of driving a MEMS resonator subjected to an external dynamic load, with a signal buried in a multi-frequency noise floor. Despite the findings of improved signal identification when the system is actively driven near its natural frequency, there was a result showing noise in the output signal, requiring post-processing for improved signal-to-noise ratios [36]. Qiao et al. (2022) analyzed the response to thermal noise in electrostatic microcantilever MEMS sensors and found that mechanical–thermal and electrical–thermal noise dominated the output signal, especially when the stability of co-existing system motions were significantly smaller than the noise-driven motions [37].
Various strategies to address the nonlinearity and noise in dynamically actuated PD-MCs have been attempted, with particular emphasis on the development of mathematical response process models for nonlinearity compensation [38,39,40,41] and microcantilever non-linear vibration control [42,43]. However, they all rely on the determination of corrective system parameters during calibration, which are kept constant over time. This makes them unsuitable for systems whose mechanical characteristics—such as rigidity and Young’s modulus of the functionalization layer—evolve with usage. Alternative approaches (that have been advanced by some studies) are focused on directly addressing nonlinearity in PD-MCs response signals from specific sources, including the development of mathematical response models to address the fluctuations in PD-MC response signals originating from molecular binding and unbinding interactions with functionalized surfaces [44]. The rigidity of adsorbates on the functionalized surface has also been studied, resulting in a method to accurately determine the adsorbed mass, inhomogeneously distributed over the cantilever surface [45]. A simulation model for adsorption-induced stress due to molecular adsorption on the microcantilever surface has also been developed, taking into account the inter-molecular interactions between adsorbed atoms and functionalized surfaces [46]. The predicted deflection from the stress affects resonance frequency and can be used to accurately estimate the adsorbed mass [47]. Research on the response of dynamically actuated PD-MCs has also revealed that fluctuations in the response signal are caused by the noise present in the actuating and readout equipment used [48]. Although the studies enumerated have provided models for correcting the nonlinearity in PD-MC-based sensor response signals and improving on adsorbed mass estimation, they are all modeled around structural parameters pertinent to the PD-MC system and neglect contributions from actuating and readout equipment.
This study aims to demonstrate that by parameterizing the resonance frequency response to the mass added of a complete piezoelectrically driven microcantilever-based sensor system into a polynomial of the form:
F M = α M 3 + β M 2 + γ M + λ
where M is the total mass of the PD-MC, F m is the resonance frequency response from the system, and α , β , γ , and λ are polynomial coefficients determined experimentally, it is possible to estimate the mass of adsorbed target gas. It is our hypothesis the resulting polynomial is a true response of the system with all source response signal disturbances accounted for. Therefore, instead of modeling all of the possible sources of noise and nonlinearity, the response of the PD-MC sensor system to known masses can be used to develop a mathematical response model on which the estimation of the mass of the target gas molecules adsorbed could be based. We limit our study to the adsorption of 1,5-diaminopentane(Cadaverine) on a functionalized PD-MC, based on our previous studies where a microcantilever-based sensor for meat and fish freshness was developed [8,9].
Subsequent parts of this article are organized as follows: In Section 2, we carry out a free and forced numerical vibration analysis for a chosen cantilever design with desired resonance frequency. In this analysis, we derive general PDEs for a stepped beam, which can be extended to beams or arbitrary shapes and sizes. We then transform the resultant model into an equivalent electrical model and simulate resonance frequency to the added mass response. The section is concluded with experiments on the adsorption of cadaverine using a fabricated microcantilever. In Section 3, we present the results obtained from the numerical and experimental studies, and we discuss them in Section 4. In Section 5, we present our conclusions and recommend ideas for future research.

2. Materials and Methods

In order to parameterize the non-linear cantilever response to adsorbed mass, we used an analytical approach to determine the free and forced vibration resonance frequencies and corresponding mode shapes and verified the model using the finite element analysis in ANSYS. The forced microcantilever vibrations due to piezo-actuation were then analytically evaluated, and the consequent model was transformed into an equivalent electrical circuit based on the Butterworth–Van-Dyke (BVD) model [49], by evaluating the impedance and resonance frequency response for different masses added to the system. The experimental verification of the cantilever resonance frequency to adsorbed gas mass was then conducted. The obtained response profile was then transformed into a polynomial, to form a basis for estimating unknown masses.

2.1. Analytical Approach

Figure 1 shows the chosen cantilever baseline. The cantilever beam material is Poly-silicon with a piezoelectric material (AlN) deposited and tightly bonded on the fixed end of the microcantilever. This creates jumped non-uniformities in terms of geometrical structure and materials, making the uniform beam model inadequate for the analysis. The beam is therefore divided into N segments. Appropriate governing conditions are assumed with constant parameters for respective sections, with continuity conditions at stepped points.

2.1.1. Transverse Vibration Analysis

The analytical model used assumes that the piezoelectric-driven cantilever obeys the Euler–Bernoulli (EB) beam theory, where the thickness of the beam is considered very small compared to its length L. The microcantilever beam is considered to have a total length L, width w b , and thickness t b . The piezoelectric patch is of the length l 1 , width w p , and thickness t p . E and I denote Young’s modulus and variable moment of inertia, respectively. From this point, henceforth, subscripts p and b indicate the designation of the property to the piezoelectric layer or elastic beam layer, respectively, and n denotes the numbering of the section under consideration.
The transverse vibration equation of the beam is derived using Hamilton’s principle. Small transverse deflections w ( x , t ) are considered at any point x along the beam, to allow for the assumption of linear properties on the system [50,51]. Assuming that a plane section of the beam remains plane after deformation by axial displacement z from the neutral axis, the axial strain can be expressed as
ε x = z d 2 w d x 2
and the axial stress as
σ x = E ε x
Considering the section of the beam with the piezoelectric patch, the neutral axis is shifted. The strain in that section is modified and can be expressed as:
ε x = ( y y n ) d 2 w ( x , t ) d x 2
where y n is the term connoting the shift in the neutral axis and can be defined as:
y n = 1 2 E p t p ( t p + t b ) E p t p + E b t b
The strain energy of the beam π is expressed as:
π = 0 L d 2 d x 2 [ ( E I ) ( x ) d 2 w ( x , t ) d x 2 ] d x
Since the beam is characterized by material and structural discontinuities along its length, the variable stiffness and the moment of inertia are expressed as:
E I ( x ) = E b I b ( x ) + E p I p ( x ) I b ( x ) = 1 12 b t p 3 + H ( x ) b t b y n 2 I p ( x ) = 1 12 b t p 3 + b t p y n 2 ( 1 2 ( t b + t p ) y n ) 2 b H ( x )
The work of the driving force F p 0 is given by
δ W ( t ) = F p 0 V ( t ) 0 l 1 d 2 H ( x ) d x 2 d x
where F p 0 defines the strain developed from the application of an electric field along the electrodes of the piezoelectric element patch and is defined as
F p 0 = b E p d 31 1 2 t b + t p y n
where H ( x ) is the heaviside function.
The kinetic energy K E is written as
K E = 1 2 0 L ρ ( x ) A ( x ) d w ( x , t ) d x 2 d x
where ρ and A are the varying density and area of the cantilever beam, respectively. Using the Hamiltonian principle, we consider the total energy experienced by the vibrating system, which is expressed as:
d 2 d x 2 E ( x ) I ( x ) d 2 d x 2 w ( x , t ) + m ( x ) d 2 d t 2 w ( x , t ) = F p 0 V ( t ) 0 l 1 d 2 H ( x ) d x 2 d x
Considering free vibrations, Equation (12) is evaluated to
d 2 d x 2 E ( x ) I ( x ) d 2 d x 2 w ( x , t ) + m ( x ) d 2 d t 2 w ( x , t ) = 0
Equation (13) can be separated into its spatial and temporal components resulting in
d 2 d x 2 E ( x ) I ( x ) d 2 d x 2 w ( x , t ) / ( m ( x ) ϕ ( x ) ) = q ¨ / q ( t ) = ω 2
d 2 d x 2 E ( x ) I ( x ) d 2 d x 2 w ( x , t ) = ω 2 m ( x ) ϕ ( x )
where ω is the natural frequency in radians and can explicitly be expressed using β n as a function of resonance frequency as
ω 2 = ( β n ) 4 ( E I ) n m n
Equation (15) can further be simplified into
d 4 d x 4 ϕ n ( x ) β n 4 ϕ n ( x ) = 0
With the general solution:
ϕ n ( x ) = A n s i n β n ( x ) + B n c o s β n ( x ) + C n s i n h β n ( x ) + D n c o s h β n ( x )
where ϕ n ( x ) is the solution to the mode shape function at a given point x. By applying appropriate boundary conditions for a clamped-free beam on the first section (fixed end of the beam) (Equation (19)), solutions to the mode shape equation coefficients are obtained using Equation (20).
ϕ 1 ( 0 ) = d ϕ 1 ( 0 ) d x = 0
B 1 + D 1 = 0 , A 1 + C 1 = 0
By applying continuity equations at discontinuity points, we obtain mode shape coefficients for each n t h section using Equations (21) and (22).
ϕ n ( l n ) = ϕ n + 1 ( l n ) d ϕ n ( l n ) d x = d ϕ n + 1 ( l n ) d x ( E I ) n d 2 ϕ n ( l n ) d x 2 = ( E I ) n + 1 d 2 ϕ n + 1 ( l n ) d x 2 ( E I ) n d 3 ϕ n ( l n ) d x 3 = ( E I ) n + 1 d 3 ϕ n + 1 ( l n ) d x 3
A n s i n β n l n + B n c o s β n l n + C n s i n h β n l n + D n c o s h β n l n = A k s i n β k l n + B k c o s β k l n + C k s i n h β k l n + D k c o s h β n + 1 l n A n c o s β n l n B n s i n β n l n + C n c o s h β n l n + D n s i n h β n l n = A k c o s β k l n B k s i n β k l n + C k c o s h β k l n + D k s i n h β n + 1 l n ( E I ) n β n 2 ( A n s i n β n l n B n c o s β n l n + C n s i n h β n l n + D n c o s h β n l n ) = ( E I ) k β k 2 ( A k s i n β k l n B k c o s β k l n + C k s i n h β k l n + D k c o s h β k l n ) ( E I ) n β n 3 ( A n c o s β n l n + B n s i n β n l n + C n c o s h β n l n + D n s i n h β n l n ) = ( E I ) k β k 3 ( A k c o s β k l n + B k s i n β k l n + C k c o s h β k l n + D k s i n h β k l n )
where k = n + 1 . Finally at the free end ( x = L ), the boundary conditions are:
d 2 ϕ N ( l N ) d x 2 = d 3 ϕ N ( l N ) d x 3 = 0
Resulting in
β N 2 ( A N s i n β N l N B N c o s β N l N + C N s i n h β N l N + D N c o s h β N l N ) = 0 β N 3 ( A N s i n β N l N B N c o s β N l N + C N s i n h β N l N + D N c o s h β N l N ) = 0
A characteristic matrix G of size ( 4 N x 4 N ) is derived from mode shape solutions equations with applied boundary conditions. Solutions to β n are obtained by multiplying G with a vector of mode shape coefficients P and equating it to zero, at resonance frequencies (Equation (25)).
G 4 N x 4 N P 1 x 4 n = 0
where
P = [ A 1 , B 1 , C 1 , D 1 , A 2 , B 2 , C 2 , D 2 , . , A N , B N , C N , D N ]
To obtain non-trivial solutions to Equation (25), varying values of β in small steps spanning a desired range are applied to matrix G and by finding zeros for its determinant, the natural frequencies can be found (see Equation (27)).
d e t [ G ( β ) ] = 0
The free vibration natural frequencies are then obtained from:
ω r 2 = ( β r ) 4 ( E I ) 1 m 1 = ( β n r ) 4 ( E I ) n m n
where β r are solution and ω r is the r t h natural frequency.
The system of ODE developed from Equations (15)–(28) is solved in MATLAB using the values given in Table 1 and the first three natural frequencies of the free vibration of the system obtained.

2.1.2. Forced Vibration Analysis

Forced vibration formulation of the system was developed by assuming the application of a sinusoidal driving voltage to the electrodes of the piezoelectric patch, with a frequency close to the previously obtained natural frequencies ( ω n ) of the beam. Using results from the modal analysis, the solution of Equation (15) is assumed to be a linear combination of the normal modes of the beam, i.e.,
w ( x , t ) = r = 1 ϕ ( r ) ( x ) q ( r ) ( t )
where ϕ r is the r t h mode shape function found by solving Equation (15) using the boundary conditions (see Equations (19), (21), (23)), and q r ( t ) are generalized coordinates or modal participation coefficients. Considering Equation (12) with the damping coefficient, it is written as:
r = 1 d 2 d x 2 E ( x ) I ( x ) d 2 ϕ ( r ) ( x ) d x 2 q ( r ) ( t ) + c ( x ) ϕ ( r ) ( x ) q ˙ ( r ) ( t ) + m ( x ) ϕ ( r ) ( x ) q ¨ ( r ) ( t ) = F P ( x , t )
By multiplying the s t h mode shape on each side of the equation, the term with stiffness ( E I ) can be isolated from the equation.
l 0 l N r = 1 d 2 d x 2 E ( x ) I ( x ) d 2 ϕ ( r ) ( x ) d x 2 ϕ ( s ) ( x ) q ( r ) ( t ) + c ( x ) ϕ ( s ) ( x ) ϕ ( r ) ( x ) q ˙ ( r ) ( t ) + m ( x ) ϕ ( r ) ( x ) ϕ ( s ) ( x ) q ¨ ( r ) ( t ) = l 0 l N F P ( x , t ) ϕ ( s ) ( x ) d x = 0
By isolating and evaluating the term with stiffness alone, Equation (31) is simplified into:
l 0 l N d 2 d x 2 E ( x ) I ( x ) d 2 ϕ ( r ) ( x ) d x 2 ϕ ( s ) ( x ) q ( r ) ( t ) d x = n = 1 N l 0 l N ( E I ) n d 4 ϕ n ( r ) ( x ) d x 4 ϕ n ( s ) ( x ) q ( r ) ( t ) d x = n = 1 N l n = 1 l n ω r 2 m n ϕ n ( r ) ϕ n ( s ) ( x ) q ( r ) ( t ) d x
The term with stiffness alone is, therefore, simplified and written as:
ω r 2 q ( r ) ( t ) l 0 l N m ( x ) ϕ ( r ) ( x ) ϕ ( s ) ( x ) d x
By inserting Equation (33) into (31) and re-arranging, the forced vibration equation becomes
q ¨ ( r t ) ( t ) + s = 1 c ( x ) ϕ ( r ) ( x ) ϕ ( s ) ( x ) d x + ω r 2 q ( r ) ( t ) = l 0 l N F p ( x , t ) ϕ ( r ) ( x ) d x
Equation (34) is therefore simplified into
q ¨ ( r t ) ( t ) + s = 1 c r s q ˙ ( t ) + ω r 2 q ( r ) = F ( r ) ( t )
where
c r s = l 0 l N c ( x ) ϕ ( r ) ( x ) ϕ ( s ) ( x ) d x , F ( r ) ( t ) = F p ( x , t ) ϕ ( r ) ( x ) d x
The damping coefficient (c) of the cantilever was derived from the combination of both air damping along the length of the cantilever, and material damping, with AlN considered the only contributor due to material damping of the polysilicon, used as the substrate, being considered negligible for these simulations, owing to its low intrinsic damping [52]. Air damping is calculated using kokubun et al’s model [53], based on stake’s law to derive the air damping coefficient. Experimental data points carried out by Hartono in [54] show that the air damping coefficient is given by:
C A i r = 0 L b ( x ) 8 π m ( x ) k b T P
where b ( x ) , m ( x ) , T , k b , and P are the variable beam width, variable beam sectional mass, temperature, Boltzmann’s constant, and atmospheric pressure in Pascals, respectively, at which the damping is evaluated. Aluminum nitride (piezoelectric material) is considered the sole contributor to material damping. The material damping is determined as shown in Equation (38), where B is a material damping constant and has been determined to be the ratio of the material viscosity to the material Young’s modulus. For aluminum nitride, B has been determined to be 3.84 × 10 14 [55].
C m a t = B A l N π ω n
The total damping c of the system is the total of the material and the air damping as shown in Equation (26).
c = C m a t + C A i r
The ODE resulting from Equation (35) was transformed into the state space representation and solved in MATLAB for the first three natural frequencies.

2.2. FEM Analysis in ANSYS

A geometric model of the PD-MC with dimensions used in the numerical analysis was made using Ansys parametric design language (APDL). Two bodies were designed (the elastic beam body and the piezoelectric body) with the respective materials defined. The resultant volumes were glued together with the piezoelectric volume at the clamped end of the beam. Boundary conditions were applied to the resulting geometric model and imported to the ANSYS Workbench 19.2 modal module, where a mesh was applied and corresponding modal resonance frequencies were obtained. The result from the modal module was then connected to the harmonic response module and resonance frequency was obtained using the largest displacement of the node at the midpoint of the free end edge. Since the Euler–Bernoulli beam theory was applied, resonance frequencies due to transverse vibration without shear deformations were the only ones considered.

2.3. Equivalent Circuit Analysis

The piezoelectrically driven microcantilever was considered to be driven electrically by applying a known sinusoidal voltage at a frequency spanning its resonance frequency. The response obtained is in the form of complex impedance and the phase output, from which the resonance frequency can be deduced. The electrical equivalent model of the piezoelectric activated cantilever was, therefore, derived to allow for the analysis of the complete electromechanical system interactions for all modes of interest. The equivalent circuit model is then solved numerically, consequent to the determination of the lumped parameters of the equivalent circuit model.
The used model for the electrical response is the Butterworth–Van Dyke model [49], with lumped element parameters, as shown in Figure 2, where the lumped element parameters of the equivalent circuit are the inductance ( L m ), resistance ( R m ), and capacitance ( C m ), which form the motional arm of the circuit. Parallel capacitance C 0 forms the parallel arm of the circuit and is equivalent to the static capacitance of the piezoelectric patch between the electrodes.
The equivalent circuit modeling was evaluated using the outputs of the analytical model. The static capacitance C 0 can be defined according to the capacitances constitutive Equation (40) [56,57] where A P , t p , ϵ 33 T , d 31 and S 11 E are the areas covered by the piezoelectric patch, thickness of the piezoelectric patch, dielectric constant under constant stress, piezoelectric strain coefficient, and the elastic compliance under the constant electric field, respectively.
C 0 = A p t p ϵ 33 T d 31 2 S 11 E
R m is dubbed the motional resistance. According to circuit theory, R m is reduced to its minimum value at the resonant frequency of ω n and is inversely proportional to the quality factor ( Q ) [57], which can be evaluated as in Equation (41). Using the solutions to the mode shape equations ϕ n ( x ) resulting from analytical vibration analysis, the R m is defined as in Equation (43).
Q = 1 2 ζ
R m = ( s 11 E w p d 31 ) 2 ( 2 t p m n e ( E b t b + E p t p ) ( E p t p + E p t b ) E b t b t p ϕ ´ n ( l 1 ) Q ( t p + t b ) 2 0 l 1 2 ϕ n ( x ) x 2 d x )
where ζ is the damping coefficient, m e is the effective mass of the microcantilever, and ϕ ˙ n ( l 1 ) is the spatial derivative of the deflection of the cantilever at position l 1 (the endpoint of the piezoelectric element patch). The motional inductance L m and motional capacitance C m can be defined in terms of R m [49] as follows:
L m = R m Q w 0
C m = 1 ( w 0 R m Q )
The equivalent resistance and reactance are expressed as
R e = R m / ( ω n 2 C 0 2 ) R m 2 + ω n L m 1 ω n C m 1 ω n C 0 2
X e = 1 ω n C 0 R m 2 + ( ω n L m 1 ω n C m ) + ( ω n L m 1 ω n C m 1 ω n C 0 ) R m 2 + ( ω n L m 1 ω n C m 1 ω n C 0 ) 2
Moreover, the complex impedance is evaluated as:
Z = R e + j X e
The electronic response in the form of impedance and phase of the cantilever was evaluated using MATLAB, with the basis on the simulated cantilever deflections from the forced vibration analysis. A range of masses from 0–200 μ g was added to the free end of the cantilever, to simulate the change in resonance frequency, and the changes in the impedance and phase response curves with the mass are recorded.

2.4. Experimental Verification

A microcantilever provided by Fraunhofer ISIT with the same dimensions as those used in the analytical model formulation and with a measured resonance frequency of 10.02 kHz was used. The impedance and phase response were measured using an Agilent 42941A precision impedance analyzer [58] with a measurement tolerance of ±3.0 Hz. The impedance analyzer was set to use the Butterworth–Van Dyke equivalent circuit model, with a driving source voltage of 1 V. The impedance analyzer was set to use the Butterworth–Van Dyke equivalent circuit model, with a driving source voltage of 1 V. For ease of reading, we use the term “calibration” to mean the application of known masses to obtain a resonance frequency response that is later used to develop the response polynomial.
Two sets of calibrations were done, one with calibration generic calibration masses and the other with measured masses of the binder solution. The provided binder solution had strong adhesive properties with glue-like characteristics to enable it to stick to the microcantilever surface. The response profiles from the two sets of functionalization masses were, therefore, expected to be different owing to the difference in surface stress each category of mass induces on the microcantilever beam.
For the first set of calibration, henceforth referred to as “calibration set-1”, the microcantilever was incrementally dozed with 0.01 μ L up to a homogeneous solution of density 2.1 g/mL using a Hamilton Gastight Syringe, mounted on the OCA 15 EC, optical contact angle meter from data physics [59]. The maximum mass used for calibration was 271.96 μ g. The corresponding mass to each dozing was calculated, and the resulting resonance frequency was extracted from the impedance analyzer (see Figure 3 ). The performance of the obtained polynomial was tested by adding three different volumes of 98 % purity 1-dodecanethiol obtained from Sigma Aldrich, equating to 169 μ g, 338 μ g, and 507 μ g, while recording the resulting resonance frequency.
For the second set of calibration, henceforth referred to as “calibration set-2”, the microcantilever was cleaned in an ultrasonic bath of isopropanol and left for 1 h to dry, prior to each functionalization. The clean microcantilever was then functionalized by drop-casting with the binder solution of 1,5-pentanediamine (cadaverine), provided by AmiNIC ApS. The same setup for the micro-dosing set-up as that in Figure 4. For each functionalization, the mass of the binder on the microcantilever was measured by thermogravimetric analysis, using a TGA-550, from Waters/Thermogravimetric analyzers Instruments [60], with a measurement tolerance of ±0.1 μ g. The measured masses were 2.1 μ g, 7.3 μ g, 13.3 μ g, and 18.6 μ g. The resonance frequency response of the complete setup was measured and recorded for each mass of the binder solution. The resulting data points of resonance frequency response to the added mass from both sets of calibration were then fit to respective polynomials using P o l y f i t and P o l y v a l functions in MATLAB, and the values to the α , β , γ , and λ coefficients to Equation (2) were obtained. The norm of residuals was used to determine the goodness of fit of the polynomial, with a lower norm signifying a better fit. The performance of the resulting polynomial was tested by adding 3.6 μ g, 3.1   μ g, and 3.1 μ g of the binder solution to the cleaned microcantilever. The resonance frequency for each mass was recorded and estimated using the polynomial.
For the estimation of the mass of adsorbed cadaverine, a functionalized microcantilever with 3.4 μ g of the binder solution was exposed to incremental volumes of 0.2 mL, 0.3 mL, and 0.4 mL of 95 % pure cadaverine, purchased from Sigma Aldrich, in a closed 750 mL container. The concentrations of the cadaverine for each volume, in parts per million-volume ( p p m v ), were calculated to be 0.0551 , 0.0827 , 0.1103 p p m v using Equation (48) as described in [61].
C o ( p p m v ) = 22.4 ρ T V 273.15 M V × 1000
where C o is the sample concentration, ρ is the density of the sample, V is the volume of sample put in the container, T is the temperature in Kelvins, M is the molecular weight of the sample, and V is the volume of the chamber in Litres.

3. Results

3.1. Analytical Vibration Analysis

Table 2 shows the resonance frequencies obtained from the numerical analysis by finding solutions to Equation (25), in comparison to those obtained from the simulation in ANSYS. The error between the simulation results and the analytical modal analysis results is also shown.

3.2. Simulation of the Resonance Frequency Response of Microcantilever to the Added Mass

Consequent to the forced vibration analysis, an analysis of the change in resonance frequency with added mass, and the numerical analysis of the corresponding electrical impedance and phase response was performed by adding masses between 0 and 200 μ g to the free end of the microcantilever model. For ease of the simulation, the homogeneous distribution of the masses on the cantilever surface was assumed.
Figure 5 shows the profile of the change of resonance frequency with mass added to the cantilever. There is a linear trend for very small masses, i.e., between 0 and 10 μ g. The response, however, becomes non-linear with an increase in mass, in a polynomial-like trend. Based on our hypothesis, we see a possibility of parameterizing this response trend into a polynomial, to provide a basis on which mass adsorbed on the cantilever can be determined.

3.3. Experimental Results

A plot of the actuated microcantilever resonance frequency response from calibration set-1 (where generic masses are used) shows a reduction in the resonance frequency with the increase in the micro-doses of known masses (see Figure 6). A linear connection of the resonance frequency data points shows varying gradients between each point (dotted line). By fitting the resonance frequency response data points into a polynomial using P o l y f i t and P o l y v a l functions in MATLAB, a third-order polynomial was obtained. The resulting norm of residuals for the polynomial fit was 0.16309 (see Figure 6 solid line). Table 3 shows the resonance frequencies read from the equipment consequent test masses added, in comparison with estimated masses based on the previously obtained resonance frequency to the added mass response polynomial. The table also shows the estimation of added masses outside the calibration range.
Using the resulting polynomial to estimate the mass when a test mass of 69 μ g of the test mass is used, shows an estimated mass of 170 μ g, giving a difference of 1.0 μ g (Figure 7). Moreover, 169 μ g is in the range of calibration masses, ranging from 0 to 271.96 μ g. Estimating masses outside of the calibration range shows a significant difference in the masses estimated and known masses added. The difference increases with the increase in mass added.
Figure 8 shows the response polynomial resulting from the use of known masses of functionalization binder solution. The resulting norm of residuals for the polynomial fit was 0.056417 . The results from testing the resultant response polynomial with measured masses of 3.1 μ g, 3.4 μ g, and 4.0 μ g of the binder solution show estimated masses of 3.3 μ g, 3.6 μ g, and 4.07 μ g, respectively. Table 4 and Figure 9 show the estimations of the added masses using the resulting polynomial.
Figure 10 shows the comparison between the parameterized resonance frequency to the added mass response resultant polynomial from calibration set-1 using known masses, and calibration set-2 using the binder solution for the first 20 μ g. This is because calibration set-2 was done with a maximum mass of 20 μ g. The graph also shows that the analytically obtained resonance frequency response to the added mass response obtained analytically is lower than those obtained from both parameterized polynomial responses from calibration set-1 and calibration set-2. The responses from the two parameterized polynomials also show significant differences in the rate of the resonance frequency drop-off with mass added. The response polynomial from calibration set-1 has a lower drop-off rate compared to the response polynomial from calibration set-2.
Figure 11 shows the resonance frequencies obtained when a functionalized cantilever with 3.4 μ g of the binder was exposed to 0.0551 , 0.0827 , 0.1103 m v calculated concentrations of the cadaverine were 8.9618 kHz, 8.8833 kHz, and 8.8520 kHz, respectively. These resonance frequencies translated to 4.7 μ g, 5.7 μ g, and 6.2 μ g of cadaverine adsorbed (see Table 5).

4. Discussion

The analytical simulation of the micro-cantilever resonance frequency to the added mass response showed a non-linear polynomial-like trend with an increase in the mass. The results confirmed the possibility of applying a polynomial fit to the actuated microcantilever resonance frequency to the added mass response. The rate of change of the resonance frequency response with the known added mass in the polynomial realized from calibration set-1 is lower than that realized from calibration set-2 and the simulated response. This is attributed to the differences in damping accruing from each setup. The analytical model used estimates the damping from predefined models of air and material damping, excluding non-classical sources that were not accounted for. There are other contributions to damping especially those originating from non-classical damping contributions, such as the coupling of the piezoelectric layer dynamics to the mechanical dynamics, for instance, the counter strain due to self-induced voltage [62], adsorption-induced damping [45], and shear deformations, especially for beams that are widely compared to the length [63]. There is, therefore, a risk of always over-estimating or under-estimating the damping coefficient. This phenomenon has been widely investigated and different methods of approximations were developed (although they are lacking in accuracy) [64].
Most of the piezoelectrically driven microcantilever sensor system damping is challenging to model accurately; therefore, it is not fully accounted for in most analytical and simulation models. By obtaining the response of a piezoelectrically driven microcantilever in a complete setup with actuating and readout equipment, it was possible to deduce the actual response of the system combined with all sources of noise, damping, and electrical losses. By parameterizing the resonance of the microcantilever resonance frequency to the added mass, a third-order polynomial for use in the estimation of the mass of the target gas was adsorbed.
The differences in parameterized response profiles from both calibration sets are attributed to the differences in adhesive properties of the two materials used for calibration. The binder solution used in calibration set-2 has glue-like sticky properties to allow for immobilization when applied to the microcantilever surface. This results in the introduction of surface stresses on the microcantilever, thereby damping the actuated vibrations. The damped actuated vibrations in turn result in damped resonance frequencies. This observation led to the deduction that the material used for calibration should have microcantilever surface adhesive properties similar to those of the material used in the sensing application.
On testing the performances of the two microcantilever response polynomials to added mass, it was found that each polynomial could only be used for the estimation of masses with similar adhesive properties to those of materials of known masses. Looking at the response polynomial developed from calibration set-1, a known mass of 169 μ g of 1-dodecanethiol was estimated as 170 μ g, with an error of 1 μ g. Without regard for experimental errors, this can be considered accurate for masses in the order of tens of micrograms. The highest mass used in the calibration was 271 μ g, which lies within the calibration range. An attempt to use the polynomial to estimate masses outside the calibration range was futile, with estimated masses and resonance frequencies significantly different from known test masses used in the experiment and resonance frequencies read from the equipment. This same result was observed when known masses of the binder solutions were used. The estimations of three differently known functionalization masses with sub-microgram differences were more accurate with an average error of 0.16 μ g. This improvement in accuracy was attributed to the use of the same solutions as both calibration masses and test masses. For accurate use of the resulting response polynomials, the masses to be estimated should be within the range of masses used for calibration, and the material of mass under the test should have the same adhesive properties as the microcantilever surface, such as the masses used in calibration.
Consequent to the performance of the response polynomial from calibration set-2, it was then used to estimate the mass of the adsorbed cadaverine on a functionalized microcantilever in different cadaverine concentrations. There was a reduction in the microcantilever response resonance frequency, with an increase in concentration. The estimated masses were 4.7 μ g, 5.7 μ g, and 6.2 μ g for the 0.0551, 0.0827, and 0.1103 m v calculated concentrations of cadaverine, with resonance response frequencies of 8.9618 kHz, 8.8833 kHz, and 8.8520 kHz, respectively. These results are withing the assumption that there is no desorption and that the interaction between the binder and adsorbed cadaverine molecules has a negligible effect on the viscosity of the binder. The change in the binder viscosity influences the overall surface stress, subsequently affecting the resulting resonance frequency. The estimated mass of the adsorbed cadaverine is in micrograms, which can be deemed high for a gas. However, it is important to note that synthetic cadaverine was used. The cadaverine has a molecular weight of 102.18 g/mol and a settling density of 870 kg/mol [65]. Despite the fact that the cadaverine is a volatile gas, it crystallizes when exposed to air and forms heavy particles. This made it possible to estimate the mass of the adsorbed cadaverine since the polynomial was developed using masses in the order of micrograms. For an accurate response in nanograms or picograms, it is imperative that known masses (in the order of target masses for the estimation) be used for developing the resonance frequency to the added mass response polynomial.

5. Conclusions

An analytical solution for the transverse vibration analysis of a stepped beam extendable to beams with desired structures and shapes was derived. A free and forced vibration analysis for a provided cantilever with known structures and dimensions was conducted. An Ansys simulation was carried out to verify the developed analytical model and they were proximate, especially for the first mode, which is of our interest. The resonance frequency for the first mode from the analytical solution was 10.553 kHz while that from the simulation in ANSYS was 10.553 kHz. An analytical resonance frequency to the mass added to the model was also conducted and a polynomial-like response profile of the microcantilever was obtained.
The system set-up of the piezoelectrically driven microcantilever, together with the driving and read-out equipment, was calibrated with known masses to obtain a resonance frequency to the added response. The response data points were then used as the bases (i.e., parameterizing the microcantilever response into a polynomial of the order of three), based on which unknown masses could be estimated. The approach was exclusively applicable to the estimation of masses within the range of masses used for calibration. It was also found that the accuracy of the resulting polynomial used for the estimation of adsorbed masses is highly dependent on the order of masses and adhesive properties of masses used for calibration.
It can be concluded that there is a possibility of using this approach to estimate the cadaverine adsorbed on a microcantilever. The performance of this approach is, however, limited by the known masses used for calibration.
The findings from this study are the bases on which the calibration of microcantilevers used for biogenic gas sensing can be done. This will allow for accurate estimation of the biogenic gas concentration in a sample, as derived from the mass of the gas adsorbed onto the functionalized cantilever.

Author Contributions

L.N.: Conceptualization, formal analysis, experimentation, writing, review, and editing. L.D.: methodology, writing, review, and editing. T.L.M.: investigation, visualization. S.H.: resources, validation. J.M.: resources, validation. F.L.: resources, validation, H.-G.R.: project administration, resources. R.d.O.H.: supervision, methodology, validation, writing, editing, funding acquisition. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by Danmarks Frie Forskningsfond (DFF) in relation to the project “High-accuracy Prediction of Meat Expiration Dates by Overcoming Nonlinearity Barriers for Microcantilever Biosensors” (0217-00023A-DFF-Tematisk forskning i Grøn omstilling).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors acknowledge AmiNIC ApS (https://www.aminic.eu, accessed on 1 January 2020) for providing their cadaverine-specific binder.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chaudhary, M.; Gupta, A. Microcantilever-based sensors. Def. Sci. J. 2009, 59, 634–641. [Google Scholar] [CrossRef] [Green Version]
  2. Yoo, K.; Na, K.; Joung, S.; Nahm, B.; Kang, C.; Kim, Y. Microcantilever-based biosensor for detection of various biomolecules. Jpn. J. Appl. Phys. 2006, 45, 515. [Google Scholar] [CrossRef]
  3. Vashist, S. A review of microcantilevers for sensing applications. J. Nanotechnol. 2007, 3, 1–18. [Google Scholar]
  4. Afshari, M.; Jalili, N. Towards non-linear modeling of molecular interactions arising from adsorbed biological species on the microcantilever surface. Int. J. -Non-Linear Mech. 2007, 42, 588–595. [Google Scholar] [CrossRef]
  5. Thundat, T.; Oden, P.; Warmack, R. Microcantilever sensors. Microscale Thermophys. Eng. 1997, 1, 185–199. [Google Scholar] [CrossRef]
  6. Sone, H.; Fujinuma, Y.; Hosaka, S. Picogram mass sensor using resonance frequency shift of cantilever. Jpn. J. Appl. Phys. 2004, 43, 3648. [Google Scholar] [CrossRef]
  7. Yang, J.; Ono, T.; Esashi, M. Mechanical behavior of ultrathin microcantilever. Sens. Actuators A Phys. 2000, 82, 102–107. [Google Scholar] [CrossRef]
  8. Alexi, N.; Hvam, J.; Lund, B.; Nsubuga, L.; de Oliveira Hansen, R.M.; Thamsborg, K.; Lofink, F.; Byrne, D.; Leisner, J. Potential of novel cadaverine biosensor technology to predict shelf life of chilled yellowfin tuna (Thunnus albacares). Food Control 2021, 119, 107458. [Google Scholar] [CrossRef]
  9. Alexi, N.; Thamsborg, K.; Hvam, J.; Lund, B.; Nsubuga, L.; de Oliveira Hansen, R.M.; Byrne, D.; Leisner, J. Novel cadaverine non-invasive biosensor technology on the prediction of shelf life of modified atmosphere packed pork cutlets. Meat Sci. 2022, 192, 108876. [Google Scholar] [CrossRef]
  10. Maijala, R.; Eerola, S. Contaminant lactic acid bacteria of dry sausages produce histamine and tyramine. Meat Sci. 1993, 35, 387–395. [Google Scholar] [CrossRef]
  11. Edwards, R.; Dainty, R.; Hibbard, C. The relationship of bacterial numbers and types to diamine concentration in fresh and aerobically stored beef, pork and lamb. Int. J. Food Sci. Technol. 1983, 18, 777–788. [Google Scholar] [CrossRef]
  12. Lakshmanan, R.; Xu, S.; Mutharasan, R. Impedance change as an alternate measure of resonant frequency shift of piezoelectric-excited millimeter-sized cantilever (PEMC) sensors. Sens. Actuators B Chem. 2010, 145, 601–604. [Google Scholar] [CrossRef]
  13. Zhang, X.; Yang, M.; Vafai, K.; Ozkan, C. Design and analysis of microcantilevers for biosensing applications. JALA 2003, 8, 90–93. [Google Scholar]
  14. Zhu, Q. Microcantilever sensors in biological and chemical detections. Sens. Transducers 2011, 125, 1. [Google Scholar]
  15. Narducci, M.; Figueras, E.; Lopez, M.; Gràcia, I.; Fonseca, L.; Sant, J.; Cané, C. A high sensitivity silicon microcantilever-based mass sensor. In Proceedings of the SENSORS, 2008 IEEE, Lecce, Italy, 26–29 October 2008; pp. 1127–1130. [Google Scholar]
  16. Lang, H.; Gerber, C. Microcantilever sensors. In STM And AFM Studies on (bio) Molecular Systems: Unravelling the Nanoworld; Springer: Berlin/Heidelberg, Germany, 2008; pp. 1–27. [Google Scholar]
  17. Zinoviev, K.; Plaza, J.; Cadarso, V.; Dominguez, C.; Lechuga, L. Optical biosensor based on arrays of waveguided microcantilevers. Silicon Photonics II 2007, 6477, 348–360. [Google Scholar]
  18. Xu, J.; Bertke, M.; Wasisto, H.; Peiner, E. Piezoresistive microcantilevers for humidity sensing. J. Micromechanics Microengineering 2019, 29, 053003. [Google Scholar] [CrossRef]
  19. Raiteri, R.; Grattarola, M.; Butt, H.; Skládal, P. Micromechanical cantilever-based biosensors. Sens. Actuators B Chem. 2001, 79, 115–126. [Google Scholar] [CrossRef]
  20. Debéda, H.; Dufour, I. Chapter 9-Resonant microcantilever devices for gas sensing. In Collection Advanced Nanomaterials for Inexpensive Gas Microsensors; Elsevier: Amsterdam, The Netherlands, 2020; pp. 161–188. Available online: https://www.sciencedirect.com/science/article/pii/B9780128148273000098 (accessed on 19 December 2022).
  21. Dong, Y.; Gao, W.; Zhou, Q.; Zheng, Y.; You, Z. Characterization of the gas sensors based on polymer-coated resonant microcantilevers for the detection of volatile organic compounds. Anal. Chim. Acta 2010, 671, 85–91. [Google Scholar] [CrossRef]
  22. Xu, J.; Bertke, M.; Setiono, A.; Schmidt, A.; Peiner, E. ZNO Nanostructures Functionalized Piezoresistive Silicon Microcantilever Platform for Portable Gas Sensing. In Proceedings of the 2019 IEEE 32nd International Conference On Micro Electro Mechanical Systems (MEMS), Seoul, Republic of Korea, 27–31 January 2019; pp. 456–459. [Google Scholar]
  23. Hosseini, S.; Kalhori, H.; Shooshtari, A.; Mahmoodi, S. Analytical solution for nonlinear forced response of a viscoelastic piezoelectric cantilever beam resting on a nonlinear elastic foundation to an external harmonic excitation. Compos. Part B Eng. 2014, 67, 464–471. [Google Scholar] [CrossRef]
  24. Hosseini, S.; Shooshtari, A.; Kalhori, H.; Mahmoodi, S. Nonlinear-forced vibrations of piezoelectrically actuated viscoelastic cantilevers. Nonlinear Dyn. 2014, 78, 571–583. [Google Scholar] [CrossRef]
  25. Mahmoodi, S.; Daqaq, M.; Jalili, N. On the nonlinear-flexural response of piezoelectrically driven microcantilever sensors. Sens. Actuators A Phys. 2009, 153, 171–179. [Google Scholar] [CrossRef]
  26. Mahmoodi, S.; Jalili, N.; Daqaq, M. Modeling, nonlinear dynamics, and identification of a piezoelectrically actuated microcantilever sensor. IEEE/ASME Trans. Mechatron. 2008, 13, 58–65. [Google Scholar] [CrossRef]
  27. Ke, L.; Wang, Y.; Wang, Z. Nonlinear vibration of the piezoelectric nanobeams based on the nonlocal theory. Compos. Struct. 2012, 94, 2038–2047. [Google Scholar] [CrossRef]
  28. Ribeiro, P.; Thomas, O. Nonlinear modes of vibration and internal resonances in nonlocal beams. J. Comput. Nonlinear Dyn. 2017, 12, 3. [Google Scholar] [CrossRef]
  29. Yong, Y.; Moheimani, S.; Kenton, K. Invited review article: High-speed flexure-guided nanopositioning: Mechanical design and control issues. Rev. Sci. Instrum. 2012, 83, 121101. [Google Scholar] [CrossRef] [PubMed]
  30. Villanueva, L.; Karabalin, R.; Matheny, M.; Chi, D.; Sader, J.; Roukes, M. Nonlinearity in nanomechanical cantilevers. Phys. Rev. B 2013, 87, 024304. [Google Scholar] [CrossRef]
  31. Venstra, W.; Capener, M.; Elliott, S. Nanomechanical gas sensing with nonlinear resonant cantilevers. Nanotechnology 2014, 25, 425501. [Google Scholar] [CrossRef] [PubMed]
  32. Burra, R.; Vankara, J.; Reddy, D. Effect of added mass using resonant peak shifting technique. J. Micro/Nanolithography Mems. Moems 2012, 11, 021203. [Google Scholar] [CrossRef]
  33. Rabih, A.; Dennis, J.; Khir, M.; Abdullah, M. Mass detection using a macro-scale piezoelectric bimorph cantilever. In Proceedings of the 2013 IEEE International Conference on Smart Instrumentation, Measurement and Applications (ICSIMA), Kuala Lumpur, Malaysia, 25–27 November 2013; pp. 1–6. [Google Scholar]
  34. Subhashini, S.; Juliet, A. CO2 Gas Sensor Using Resonant Frequency Changes in Micro-Cantilever. In Computer Networks & Communications (NetCom); Springer: New York, NY, USA, 2013; pp. 75–80. [Google Scholar]
  35. Bouchaala, A.; Nayfeh, A.; Younis, M. Analytical study of the frequency shifts of micro and nano clamped–clamped beam resonators due to an added mass. Meccanica 2017, 52, 333–348. [Google Scholar] [CrossRef]
  36. Lam, N.; Hayashi, S.; Gutschmidt, S. A novel MEMS sensor concept to improve signal-to-noise ratios. Int. J. -Non-Linear Mech. 2022, 139, 103863. [Google Scholar] [CrossRef]
  37. Qiao, Y.; Wei, W.; Arabi, M.; Xu, W.; Abdel-Rahman, E. Analysis of response to thermal noise in electrostatic MEMS bifurcation sensors. Nonlinear Dyn. 2022, 107, 33–49. [Google Scholar] [CrossRef]
  38. Karparvarfard, S.; Asghari, M.; Vatankhah, R. A geometrically nonlinear beam model based on the second strain gradient theory. Int. J. Eng. Sci. 2015, 91, 63–75. [Google Scholar] [CrossRef]
  39. Kong, S.; Zhou, S.; Nie, Z.; Wang, K. The size-dependent natural frequency of Bernoulli–Euler micro-beams. Int. J. Eng. Sci. 2008, 46, 427–437. [Google Scholar] [CrossRef]
  40. Park, S.; Gao, X. Bernoulli–Euler beam model based on a modified couple stress theory. J. Micromechanics Microengineering 2006, 16, 2355. [Google Scholar] [CrossRef]
  41. Awrejcewicz, J.; Krysko, V.; Pavlov, S.; Zhigalov, M.; Krysko, A. Mathematical model of a three-layer micro-and nano-beams based on the hypotheses of the Grigolyuk–Chulkov and the modified couple stress theory. Int. J. Solids Struct. 2017, 117, 39–50. [Google Scholar] [CrossRef]
  42. Sarrafan, A.; Zareh, S.; Zabihollah, A.; Khayyat, A. Intelligent vibration control of micro-cantilever beam in MEMS. In Proceedings of the 2011 IEEE International Conference on Mechatronics, Istanbul, Turkey, 13–15 April 2011; pp. 336–341. [Google Scholar]
  43. Vatankhah, R.; Karami, F.; Salarieh, H.; Alasty, A. Stabilization of a vibrating non-classical micro-cantilever using electrostatic actuation. Sci. Iran. 2013, 20, 1824–1831. [Google Scholar]
  44. Jokić, I.; Djurić, Z.; Frantlović, M.; Radulović, K.; Krstajić, P. Fluctuations of the mass adsorbed on microcantilever sensor surface in liquid-phase chemical and biochemical detection. Microelectron. Eng. 2012, 97, 396–399. [Google Scholar] [CrossRef]
  45. Chen, G.; Thundat, T.; Wachter, E.; Warmack, R. Adsorption-induced surface stress and its effects on resonance frequency of microcantilevers. J. Appl. Phys. 1995, 77, 3618–3622. [Google Scholar] [CrossRef]
  46. Dareing, D.; Thundat, T. Simulation of adsorption-induced stress of a microcantilever sensor. J. Appl. Phys. 2005, 97, 043525–043526. [Google Scholar] [CrossRef]
  47. Ramos, D.; Calleja, M.; Mertens, J.; Zaballos, Á.; Tamayo, J. Measurement of the mass and rigidity of adsorbates on a microcantilever sensor. Sensors 2007, 7, 1834–1845. [Google Scholar] [CrossRef] [Green Version]
  48. Mouro, J.; Pinto, R.; Paoletti, P.; Tiribilli, B. Microcantilever: Dynamical Response for Mass Sensing and Fluid Characterization. Sensors 2021, 21, 115. [Google Scholar] [CrossRef]
  49. Sherrit, S.; Wiederick, H.; Mukherjee, B. Accurate equivalent circuits for unloaded piezoelectric resonators. In Proceedings of the 1997 IEEE Ultrasonics Symposium Proceedings. An International Symposium (Cat. No. 97CH36118), Toronto, ON, Canada, 5–8 October 1997; Volume 2, pp. 931–935. [Google Scholar]
  50. Jalili, N. Piezoelectric-Based Vibration Control: From Macro to Micro/Nano Scale Systems; Springer Science & Business Media: New York, NY, USA, 2009. [Google Scholar]
  51. Rao, S. Vibration of Continuous Systems; John Wiley & Sons: Hoboken, NJ, USA, 2019. [Google Scholar]
  52. Sosale, G.; Das, K.; Fréchette, L.; Vengallatore, S. Controlling damping and quality factors of silicon microcantilevers by selective metallization. J. Micromechanics Microengineering 2011, 21, 105010. [Google Scholar] [CrossRef]
  53. Kokubun, K.; Hirata, M.; Murakami, H.; Toda, Y.; Ono, M. A bending and stretching mode crystal oscillator as a friction vacuum gauge. Vacuum 1984, 34, 731–735. [Google Scholar] [CrossRef]
  54. Nor, N.; Khalid, N.; Osman, R.; Sauli, Z. Estimation of Material Damping Coefficients of AlN for Film Bulk Acoustic Wave Resonator. Mater. Sci. Forum 2015, 819, 209–214. [Google Scholar] [CrossRef]
  55. Uzunov, I.; Terzieva, M.; Nikolova, B.; Gaydazhiev, D. Extraction of modified butterworth—Van Dyke model of FBAR based on FEM analysis. In Proceedings of the 2017 XXVI International Scientific Conference Electronics (ET), Sozopol, Bulgaria, 13–15 September 2017; pp. 1–4. [Google Scholar]
  56. Zhao, J.; Liu, S.; Huang, Y.; Gao, R. A new impedance based sensitivity model of piezoelectric resonant cantilever sensor. In Proceedings of the World Congress on Advances In Structural Engineering And Mechanics, Ilsan(Seoul), Republic of Korea, 28 August–1 September 2017. [Google Scholar]
  57. Burda, I. Quartz Crystal Microbalance with Impedance Analysis Based on Virtual Instruments: Experimental Study. Sensors 2022, 22, 1506. [Google Scholar] [CrossRef]
  58. Keysight 4294a Precision Impedance Analyzer, 40 Hz to 110 mhz [obsolete]. Keysight 2021, 2. Available online: https://www.keysight.com/us/en/product/4294A/precision-impedance-analyzer-40-hz-to-110-mhz.html (accessed on 19 December 2022).
  59. Neurtek dat_oca_15ec_en. Dataphysics Products. Available online: https://www.neurtek.com/descargas/dat_oca_15ec_en.pdf (accessed on 19 December 2022).
  60. TA-Instruments tga-55. TA Instruments. Available online: https://www.tainstruments.com/tga-55/ (accessed on 19 December 2022).
  61. Dipak, P.; Tiwari, D.; Samadhiya, A.; Kumar, N.; Biswajit, T.; Singh, P.; Tiwari, R. Synthesis of polyaniline (printable nanoink) gas sensor for the detection of ammonia gas. J. Mater. Sci. Mater. Electron. 2020, 31, 22512–22521. Available online: https://link.springer.com/article/10.1007/s10854-020-04760-2 (accessed on 22 November 2022).
  62. Hagood, N.; Chung, W.; Von Flotow, A. Modelling of Piezoelectric Actuator Dynamics for Active Structural Control. J. Intell. Mater. Syst. Struct. 1990, 1, 327–354. [Google Scholar] [CrossRef]
  63. Reddy, R.; Panda, S.; Gupta, A. Nonlinear dynamics and active control of smart beams using shear/extensional mode piezoelectric actuators. Int. J. Mech. Sci. 2021, 204, 106495. [Google Scholar] [CrossRef]
  64. Mevada, H.; Patel, D. Experimental Determination of Structural Damping of Different Materials. Procedia Eng. 2016, 144, 110–115. [Google Scholar] [CrossRef]
  65. Pubchem Cadaverine. National Center For Biotechnology Information. PubChem Compound Database. Available online: https://pubchem.ncbi.nlm.nih.gov/compound/Cadaverine (accessed on 21 December 2022).
Figure 1. Chosen microcantilever beam design used for the baseline simulations.
Figure 1. Chosen microcantilever beam design used for the baseline simulations.
Sensors 23 01093 g001
Figure 2. Lumped element for the Butterworth–Van Dyke model.
Figure 2. Lumped element for the Butterworth–Van Dyke model.
Sensors 23 01093 g002
Figure 3. The set-up used for the calibrating microcantilever response to mass added.
Figure 3. The set-up used for the calibrating microcantilever response to mass added.
Sensors 23 01093 g003
Figure 4. (a) Piezoelectrically driven microcantilever beam mounted on the PCB board. (b) Functionalization solution on the microcantilever before exposure to the cadaverine. (c) After exposure to the cadaverine.
Figure 4. (a) Piezoelectrically driven microcantilever beam mounted on the PCB board. (b) Functionalization solution on the microcantilever before exposure to the cadaverine. (c) After exposure to the cadaverine.
Sensors 23 01093 g004
Figure 5. Simulated polynomial-like microcantilever resonance frequency response to the mass added, for the modeled piezoelectrically driven microcantilever beam.
Figure 5. Simulated polynomial-like microcantilever resonance frequency response to the mass added, for the modeled piezoelectrically driven microcantilever beam.
Sensors 23 01093 g005
Figure 6. Parameterized resonance frequency to the added mass response calibration set-1 using generic masses.
Figure 6. Parameterized resonance frequency to the added mass response calibration set-1 using generic masses.
Sensors 23 01093 g006
Figure 7. Estimation of 169 μ g of added mass using the polynomial from calibration set-1.
Figure 7. Estimation of 169 μ g of added mass using the polynomial from calibration set-1.
Sensors 23 01093 g007
Figure 8. Calibration set-2 resonance frequency response to the added mass response of the microcantilever characterized into a polynomial (solid line).
Figure 8. Calibration set-2 resonance frequency response to the added mass response of the microcantilever characterized into a polynomial (solid line).
Sensors 23 01093 g008
Figure 9. Estimation of the known test masses of the binder solution using the polynomial from calibration set-2.
Figure 9. Estimation of the known test masses of the binder solution using the polynomial from calibration set-2.
Sensors 23 01093 g009
Figure 10. Comparison of simulated resonance frequency-mass response to parameterised response from calibration set-1 and calibration set-2.
Figure 10. Comparison of simulated resonance frequency-mass response to parameterised response from calibration set-1 and calibration set-2.
Sensors 23 01093 g010
Figure 11. Estimation of the adsorbed mass on the microcantilever, using the polynomial from calibration set-2, after exposure to different concentrations of cadaverine.
Figure 11. Estimation of the adsorbed mass on the microcantilever, using the polynomial from calibration set-2, after exposure to different concentrations of cadaverine.
Sensors 23 01093 g011
Table 1. Values used to solve ODEs in MATLAB.
Table 1. Values used to solve ODEs in MATLAB.
PropertyValueunit
l10.5mm
l21.0mm
w b 1.0mm
w p 1.0mm
tb20.0 μ m
tp2.0 μ m
E b 70.0GPa
E P 302.0Gpa
d31−1.9159 × 10−12C/N
ɛ P 11.5 × 10−12-
ρ P 3260Kg/m3
ρ b 2270Kg/m3
Vt5V
Table 2. Comparison of resonance frequencies of the first three cantilever transverse vibration modes.
Table 2. Comparison of resonance frequencies of the first three cantilever transverse vibration modes.
ModeResonance Frequency (kHz) Analytical ApproachResonance Frequency (kHz) Numerical ApproachDifference in (kHz)Percentage Error
Mode 110.55310.5560.0030.028
Mode 252.98752.7070.280.005
Mode 2135.37132.922.450.018
Table 3. Resonance frequency response to added mass test masses with expected resonance frequency from the parameterized response. The (*) shows values due to masses outside the calibration range.
Table 3. Resonance frequency response to added mass test masses with expected resonance frequency from the parameterized response. The (*) shows values due to masses outside the calibration range.
Known Mass ( μ g)Measured Resonance Frequency (kHz)Deduced Adsorbed Mass Based on Polynomial ( μ g)Difference in Known Mass and Deduced MassExpected Resonance Frequency to Fit Polynomial (kHz)Mass Estimation Tolerance (%)
1698.577417018.56710.59
338 *8.2645222 *116 *7.746452.25 *
507 *7.9992277 *230 *7.117483.03 *
Table 4. Response of the second set calibration polynomials based on the binder solution with test masses of the binder solution.
Table 4. Response of the second set calibration polynomials based on the binder solution with test masses of the binder solution.
Known Mass ( μ g)Measured Resonance Frequency (kHz)Deduced Adsorbed Mass Based on Polynomial ( μ g)Difference in Known Mass and Deduced MassExpected Resonance Frequency to Fit Polynomial (kHz)Mass Estimation Tolerance (%)
3.19.47423.30.29.50276.06
3.49.43723.60.29.46205.56
4.09.37464.070.079.38411.72
Table 5. Estimation of the mass of cadaverine adsorbed when the functionalized microcantilever with 3.4 μ g of the binder is exposed to different concentrations of cadaverine.
Table 5. Estimation of the mass of cadaverine adsorbed when the functionalized microcantilever with 3.4 μ g of the binder is exposed to different concentrations of cadaverine.
Concentration of the Cadaverine ( ppm v g )Measured Resonance Frequency (kHz)Deduced Total on the Microcantilever Based on the Polynomial ( μ g)Deduced Adsorbed Mass of Cadaverine ( μ g)
0.05518.96188.14.7
0.08278.88339.15.7
0.11038.85209.56.2
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Nsubuga, L.; Duggen, L.; Marcondes, T.L.; Høegh, S.; Lofink, F.; Meyer, J.; Rubahn, H.-G.; de Oliveira Hansen, R. Gas Adsorption Response of Piezoelectrically Driven Microcantilever Beam Gas Sensors: Analytical, Numerical, and Experimental Characterizations. Sensors 2023, 23, 1093. https://doi.org/10.3390/s23031093

AMA Style

Nsubuga L, Duggen L, Marcondes TL, Høegh S, Lofink F, Meyer J, Rubahn H-G, de Oliveira Hansen R. Gas Adsorption Response of Piezoelectrically Driven Microcantilever Beam Gas Sensors: Analytical, Numerical, and Experimental Characterizations. Sensors. 2023; 23(3):1093. https://doi.org/10.3390/s23031093

Chicago/Turabian Style

Nsubuga, Lawrence, Lars Duggen, Tatiana Lisboa Marcondes, Simon Høegh, Fabian Lofink, Jana Meyer, Horst-Günter Rubahn, and Roana de Oliveira Hansen. 2023. "Gas Adsorption Response of Piezoelectrically Driven Microcantilever Beam Gas Sensors: Analytical, Numerical, and Experimental Characterizations" Sensors 23, no. 3: 1093. https://doi.org/10.3390/s23031093

APA Style

Nsubuga, L., Duggen, L., Marcondes, T. L., Høegh, S., Lofink, F., Meyer, J., Rubahn, H. -G., & de Oliveira Hansen, R. (2023). Gas Adsorption Response of Piezoelectrically Driven Microcantilever Beam Gas Sensors: Analytical, Numerical, and Experimental Characterizations. Sensors, 23(3), 1093. https://doi.org/10.3390/s23031093

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop