Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Two-channel Hyperspectral LiDAR with a Supercontinuum Laser Source
Next Article in Special Issue
A Field Programmable Gate Array-Based Reconfigurable Smart-Sensor Network for Wireless Monitoring of New Generation Computer Numerically Controlled Machines
Previous Article in Journal
Fluorescent Chemosensors for Toxic Organophosphorus Pesticides: A Review
Previous Article in Special Issue
Study on a Two-Dimensional Scanning Micro-Mirror and Its Application in a MOEMS Target Detector
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Investigation of the Frequency Shift of a SAD Circuit Loop and the Internal Micro-Cantilever in a Gas Sensor

Department of Precision Instruments and Mechanics, Tsinghua University, Beijing 100084, China
*
Author to whom correspondence should be addressed.
Sensors 2010, 10(7), 7044-7056; https://doi.org/10.3390/s100707044
Submission received: 15 June 2010 / Revised: 30 June 2010 / Accepted: 10 July 2010 / Published: 23 July 2010
(This article belongs to the Special Issue Intelligent Sensors - 2010)

Abstract

:
Micro-cantilever sensors for mass detection using resonance frequency have attracted considerable attention over the last decade in the field of gas sensing. For such a sensing system, an oscillator circuit loop is conventionally used to actuate the micro-cantilever, and trace the frequency shifts. In this paper, gas experiments are introduced to investigate the mechanical resonance frequency shifts of the micro-cantilever within the circuit loop(mechanical resonance frequency, MRF) and resonating frequency shifts of the electric signal in the oscillator circuit (system working frequency, SWF). A silicon beam with a piezoelectric zinc oxide layer is employed in the experiment, and a Self-Actuating-Detecting (SAD) circuit loop is built to drive the micro-cantilever and to follow the frequency shifts. The differences between the two resonating frequencies and their shifts are discussed and analyzed, and a coefficient α related to the two frequency shifts is confirmed.

Graphical Abstract

1. Introduction

Micromachined resonant devices are attracting increasing interest in the field of chemical sensor applications, due to some of their interesting properties such as sensitivity, compactness and low energy consumption [1]. Helped by an actuating and detecting circuit loop, the working frequency of the MEMS resonant devices can be followed automatically, and thus the analyte quantity can be determined [2].
In the study of Zhou et al. of Tsinghua University, a gas sensor based on a MEMS cantilever was introduced [34]. In their work, a piezoelectric micro-cantilever was used, with a sensitive layer deposited on the tip to adsorb gas molecules. The adsorption modifies the mechanical properties of the structure, and therefore the resonance frequency f0 shifts accordingly [5,6]. As the micro-cantilever has no detection structure, a self-actuating and detecting (SAD) circuit loop was designed specifically for this single-port device. The SAD circuit loop drives the micro-cantilever and traces the working frequency fr of the circuit. The details of the SAD circuit loop will be introduced in Section 2.1.
Most of the actual studies follow the same pattern as presented previously [714]: on exposure to an analyte vapor, the additional mass loading of the polymer layer decreases the mechanical resonance frequency (MRF) of the micro-cantilever. This MRF shift, noted Δf0, should be proportional to the mass loading Δm, and thus proportional to the change in the gas concentration Δc [1]. Nevertheless, in real world applications, as MRF usually is not available in operation, the system working frequency (SWF) shift of the actuating and detecting circuit loop, noted as Δfr in this study, is measured instead. The MRF shift and SWF shift are not actually identical. Li et al. first referred to the difference between the MRF shift and SWF shift, in a simulation approach [4].
An investigation is performed in this paper to demonstrate the existence of a coefficient α, which makes Δfr = α·Δf0. The experimental investigation consists of an open-loop test, where the MRF shift Δf0 is measured, and a close-loop test, where the SWF shift Δfr is recorded. The origin of this coefficient is also discussed. The results indicate that the conventional usage of SWF shift as a measurement of gas concentration is not appropriate, as the SWF shift is not identical to MRF shift. Instead, SWF shift measured should be revised by a coefficient to approach the MRF shift. The further discussion in this work also proves that this finding can be applied to all the sensing system using an actuating-detecting oscillator loop.

2. Experimental Details

2.1. Brief Description of the SAD Circuit Loop

With a presumption of small deformation, an electrical model of a piezoelectric micro-cantilever is developed [3].
The mechanical function analogy of the electrical components is presented in Equation (1). Where meffect is the effective mass of the micro-cantilever, k is the spring constant, c0 is the static capacitance of the piezoelectric layer, ξ is the damping coefficient. The k1, k2 are proportionality constants:
{ C 1 = k 1 k 2 L 1 = m effect kk 1 k 2 R 1 = 2 ξ km effect kk 1 k 2 C 0 = c 0 k 1 k 2
Therefore, the electrical components in the model can be determined by the mechanical properties of the beam, among which C1 and C0 can be seen as a constant, R1 and L1 will vary with the frequencies. Furthermore, the impedance of the micro-cantilever can be determined, as shown in Figure 2. As the impedance of C0 is much smaller than that of C1-R1-L1 branch, the impedance of the micro-cantilever is rather determined by that of C0, and has only a local minimum at the MRF f0. This means the micro-cantilever does not have a good capability of selecting the resonance frequency.
Therefore, for ameliorating the frequency selection performance of the beam, a frequency selection network module is designed. The constant capacitor C0 in the electrical model was compensated by an external tunable capacitor. In an ideal situation, the output of the network should perform as a band-pass filter, having maximum amplitude at the mechanical resonance frequency (see Figure 3).
Figure 4 shows the block diagram of the SAD circuit loop. The micro-cantilever is integrated in the module of frequency selection network, which constitutes the frequency determining unit of the oscillator circuit, and serves like a band-pass filter with a central frequency equal to f0. The output signal of the frequency selection network is amplified and then adjusted by the following phase compensator. An amplitude limiting module is used to control the amplitude of the signal. Finally, as in reality the output of frequency selection network cannot be reduced to zero at frequencies far away from f0, a band-pass filter is needed, to reduce the secondary mode oscillation of the beam.
Suppose the gain of the frequency selection network is A, and its phase-frequency curve is φa, the rest of the SAD circuit loop is considered as a feedback module, with a gain of F and phase-frequency curve φb. When the SAD circuit loop working in a close-loop mode, the Barkhausen conditions should be met [15]:
| A ( f r ) F ( f r ) | = 1 φ a ( f r ) + φ b ( f r ) = 2 π
Generally, the amplitude condition in (2) is not a crucial one, as the non-linearity in the feedback module makes the amplitude of the signal always satisfy the first equation. Therefore, the working frequency is rather determined by the second Equation of (2), the phase condition.
Figure 5 illustrates the phase-frequency curves for frequency selection network as well as the feedback module. As seen from the figure, the feedback module shows a good linearity in its phase-frequency curves. And in a neighborhood of the mechanical resonance frequency, at least in the interval of f0 ± 100 Hz, the frequency selection network also shows a good linearity. After mass loading, mechanical resonance frequency decreases, and φa (f) simply has a parallel shift of Δf0, to φa (f).

2.2. Analyte Concentration Controlling

The sensing system is mounted in a glass canister, with a thin glass slide as top window. Drops of liquid ethanol were jetted by a transfer pipette through this window to provide defined vapor concentrations to the canister. Suppose the volume of the canister is Vvessel, and liquid ethanol with a volume Vanalyte is jetted into the canister. When volatilized completely, the concentration of ethanol vapor follows Equations (3) and (4) where M is the molecular weight of the analyte, and P is the ambient atmospheric pressure:
C { mg / m 3 } = ρ analyte V analyte V vessel
C { ppm } = C { mg / m 3 } M 22.4 273 273 + T P 1011325
For the experiments, we calculate the saturated concentration of ethanol vapor, known the ambient temperature and pressure of atmosphere. Then, an appropriate volume Vanalyte is chosen, to make sure that with six liquid injections, the concentration of the ethanol vapor in the canister increases at a constant pace and is finally below the saturated level. The open-loop test and the close-loop one are performed in the same canister, with the same ambient conditions. The transfer pipette has a precision of measurement of 0.3%, which ensures the coherence of the two tests, and the comparability of the results.

2.3. Open-Loop Tests

An Open-loop test records the variation of mechanical resonance frequency f0 of the micro-cantilever. The micro-cantilever used in this work is a probe from a commercial Atomic Force Microscopy (AFM) system supplied by the Veeco® Company. It consists of a micro-machined silicon beam with a piezoelectric zinc oxide layer. The MRF of the AFM cantilever is 50,000 Hz, according to the specificationz and is measured at 66,640 Hz. Polyethyleneoxide (PEO) was chosen as the sensitive material, for its good sensitivity and selectivity to ethanol vapor [16]. An aqueous solution of PEO was jetted on top of the cantilever. After the evaporation of the solvent, the tip of the cantilever was covered with a layer of sensitive material in the form of granule (Figure 6), and the MRF of the cantilever decreases to 65,670 Hz accordingly.
The open-loop test instruments system is illustrated by Figure 7. A signal generator provides an actuating signal, and the laser beam was focused on the free-moving end of the coated cantilever. The laser beam deflection signal was captured by a photodetector, and the amplitude of the oscillation was read out from the screen. By sweeping the actuating signal near the mechanical resonance frequency, f0 can be detected. The Q-factor of the coated probe is near 130 at ambient conditions of standard atmospheric pressure and 293 K temperature.

2.4. Close-Loop Tests

In a close-loop experiment, the system working frequency fr of the SAD circuit loop was detected. Schematic diagram of the test can is illustrated in Figure 4. The resolving power of the tests is approximately 0.5 Hz. To ensure that the SAD circuit loop works in its closed-loop mode, the following two adjustments should be made. Firstly, the tunable capacitor in the frequency selection network should be adjusted to compensate C0, that is to say, at a frequency ff0, the output of the frequency selection network is equal to 0. Secondly, connect the frequency selection network with the feedback module, and adjust the phase compensator, to meet the phase condition in the Barkhausen conditions.

3. Results and Discussion

3.1. Initial Offsets

A pair of experimental datapoints can be observed in Figure 8. The liquid ethanol was jetted every 30 minutes in both of the tests. The set of triangle symbols represents the mechanical resonance frequency in the open-loop test. The step-form points record the working frequency of the SAD circuit loop as a function of time. As the volatilization takes place, the working frequency of the circuit loop decreases and at the end of the volatilization, the working frequency stabilizes gradually. As mentioned at the beginning, it is the coefficient α between the SWF shift Δfr and the MRF shift Δf0 what we exploit. It is easy to notice that in this result, the initial frequency, which corresponds to a vapor concentration of 0, is different. Such an offset of the initial frequency is caused the by the offset of the initial phase. During the adjustment of the phase compensator before a close-loop test, an error of operation is inevitable as the signal superposition is depending on subjective judgment. Therefore, an offset of initial phase was introduced, noted as φoffset. The frequency offset foffset is thereby created. The frequency offset is 38 Hz in this experiment.
As φoffset is the phase offset at f = f0 in an open-loop adjustment, from (2) we have:
φ a ( f 0 ) + φ b ( f 0 ) = 2 π + φ offset
When the system works in a close-loop mode, we also have:
φ a ( f 0 + f offset ) + φ b ( f 0 + f offset ) = 2 π
Developing (6) at f = f0, using a Taylor development:
φ a ( f 0 ) + φ a ( f 0 ) f f offset + 1 2 ! φ a 2 ( f 0 ) f 2 f offset 2 + ... + φ b ( f 0 ) + φ b ( f 0 ) f f offset + 1 2 ! φ b 2 ( f 0 ) f 2 f offset 2 + ... = 2 π
As the linearity shown in Figure 4, in the interval of this development, the high order derivatives are far less than the first order derivative. Thus:
φ a ( f 0 ) + φ a ( f 0 ) f f offset + φ b ( f 0 ) + φ b ( f 0 ) f f offset = 2 π
From (5) and (7), we have:
f offset = φ offset φ a ( f 0 ) f + φ b ( f 0 ) f
It is seen from (8) how φoffset causes foffset, the offset of the initial frequency in experimental results. The literature expression of α later on will show that this coefficient will not be affected by the initial frequency offset.

3.2. Differences between MRF Shifts and SWF Shifts

Figure 9 shows the MRF and the SWF as a function of gas concentration. After six sequential ethanol injections of equal quantity, the MRF decreases from 65,670 Hz to 65,455 Hz, and the SWF decreases from 65,708 Hz to 65,498 Hz. A result of simulation of SWF calculated from the MRF measured is also presented for comparison.
It is the ratio of the slope of two experimental curves that represents the coefficient α. The linear expression is obtained by interpolating the curves using the least-square principle. As the close-loop has a higher resolution power, its data shows a better correlation:
f 0 : y = 36.393 x + 65676 , R 2 = 0.9966 f r : y = 35.179 x + 65710 , R 2 = 0.9999 f r s i m u l a t e : y = 35.089 x + 65710 , R 2 = 0.9969
Both experimental result and simulation result demonstrate the existence of a coefficient α which makes Δfr = α·Δf0. The value of α is also calculated: αreal = 0.9666, αsimulate = 0.9642.
We are going to explore this a little bit further. As a frequency offset foffset is observed, we are going to demonstrate it has no influence to α. Suppose that the SWF of the SAD circuit loop decreases from fr to f r , when the MRF decreases form f0 to f 0 . Transform the definition f r f r = α ( f 0 f 0 ) to:
( 1 α ) ( f 0 f 0 ) = f offset f offset
After mass loading, at the new resonance frequency f 0 , a new phase offset φ offset exists:
φ a ( f 0 ) + φ b ( f 0 ) = 2 π + φ offset
From the discussion of Figure 4, there is:
φ a ( f ) = φ a ( f Δ f 0 )
A deduction follows the routine of (8), we have:
f offset = φ offset φ a ( f 0 ) f + φ b ( f 0 ) f
For the first term of the denominator, use (11), we have:
φ a ( f 0 ) f = φ a ( f 0 Δ f 0 ) f = φ a ( f 0 ) f
For the second term of the denominator, develop φ b ( f 0 ) f using Taylor development at f 0 :
φ b ( f 0 ) f = φ b ( f 0 ) f + 1 2 ! 2 φ b ( f 0 ) f 2 Δ f 0 + 1 3 ! 3 φ b ( f 0 ) f 3 Δ f 0 2 + ...
Notice the linearity of φb(f), we have:
φ b ( f 0 ) f = φ b ( f 0 ) f
Substitute the denominator of (12) by (13) and (15), using (8), we have:
f offset f offset = φ offset φ offset φ a ( f 0 ) f + φ b ( f 0 ) f
Using (5) and (10), we have:
φ offset φ offset = φ b ( f 0 ) f Δ f 0
Thus:
f offset f offset = φ b ( f 0 ) f φ a ( f 0 ) f + φ b ( f 0 ) f Δ f 0
From (18) α can be written as:
α = φ a ( f 0 ) f φ a ( f 0 ) f + φ b ( f 0 ) f
The literature expression shows that α is independent with φoffset. It is the phase-frequency curve of frequency selection network and of feedback module that determine the coefficient. As shown in Figure 5, at the point of f0, the slope of the phase-frequency curve φa (f) is much larger than that of φb (f), thereby the coefficient α has a value close to 1.
To calculate the coefficient α from Equation (19), we draw the phase-frequency curves in Figure 10, where the experimental data is presented together with the simulation one. For the frequency selection network, both the results from experiment and simulation show a good linearity at the neighborhood of f0, whereas at the point far away from f0, the slopes of the two curves deviate. The deviation is likely to be a result of non-ideal properties of the micro-cantilever. As for the feedback module, both two curves show a good linearity, and an approximate slope.
Based on literature result in (19), the coefficient α can be determined from the slopes of the phase-frequency curves: αreal = 0.9659, αsimulate = 0.9640. These results confirm our findings in gas experiments. Furthermore, as the coefficient α is related to the phase-frequency curves of the frequency selection network and of the feedback module, it is a typical value for each realization of a SAD circuit loop.

4. Conclusions

Gas experiments have been performed on a MEMS gas sensing system. The system is based on a silicon beam with a piezoelectric zinc oxide layer, and a SAD circuit loop is built to actuate the cantilever and detect its frequency shifts. The experimental results confirm, as predicted by another paper, the existence of a coefficient α, relating the mechanical resonance frequency shift and system working frequency shift. It was also discussed the relation between this coefficient and the phase-frequency curves of the SAD circuit loop.
As analyzed in this paper, the mass loading of the micro-cantilever modifies its MRF as well as the phase-frequency curves of the SAD circuit loop, and therefore the SWF shifts to meet the phase condition of Formula (2) so that the SAD circuit loop could continue its close-loop working. This paper shows that the SWF shift observed in operations is not an appropriate indicator of the gas concentration. To characterize the concentration of analyte, the SWF shift should be revised with a coefficient, to trace the MRF shift. And as such a coefficient is derived from the Barkhausen conditions, it should be taken into account in all kinds of sensing system working with an actuating-detecting circuit loop.

Acknowledgments

The authors would like to acknowledge the support by National 863 Project of China and Science Fund of China Post Doctoral Scientists.

References

  1. Zougagh, M; Rios, A. Micro-electromechanical sensors in the analytical field. Analyst 2009, 134, 1274–1290. [Google Scholar]
  2. Lange, D; Hagleitner, C. Complementary metal oxide Semiconductor cantilever arrays on a single chip: mass-sensitive detection of volatile organic compounds. Anal. Chem 2002, 74, 3084–3095. [Google Scholar]
  3. Zhou, Q. Study on MEMS chemical gas sensor based on resonating microcantilever. M.S. thesis; Tsinghua University: Beijing, China, 2006. [Google Scholar]
  4. Li, P; Zhao, J. Resonating frequency of a SAD circuit loop and inner micro-cantilever in a gas sensor. IEEE Sens. J. (In press)
  5. Datar, R; Kim, S; Jeon, S; Hesketh, P; Manalis, S; Boisen, A; Thundat, T. Cantilever sensors: Nanomechanical tools for diagnostics. MRS Bull 2009, 34, 449–454. [Google Scholar]
  6. Finot, E; Passian, A; Thundat, T. Measurement of mechanical properties of cantilever shaped materials. Sensors 2008, 8, 3497–3541. [Google Scholar]
  7. Fadel, L; Lochon, F. Chemical sensing: millimeter size resonant microcantilever performance. J. Micromech. Microeng 2004, 14, S23–S30. [Google Scholar]
  8. Bedair, S; Fedder, G. CMOS MEMS oscillator for gas chemical detection. Proc. IEEE 2004, 2, 955–958. [Google Scholar]
  9. Vancura, C; Ruegg, M. Magnetically actuated complementary metal oxide semiconductor resonant cantilever gas sensor systems. Anal. Chem 2005, 77, 2690–2699. [Google Scholar]
  10. Li, Y; Vancura, C. Monolithic CMOS multi-tranducer gas sensor microsystem for organic and inorganic analytes. Sens. Actuat 2007, B126, 431–440. [Google Scholar]
  11. Lu, J; Ikehara, T. Characterization and improvement on quality factor of microcantilevers with self-actuation and self-sensing capability. Microelectron. Eng 2009, 86, 1208–1211. [Google Scholar]
  12. Mutyala, MSK; Bandhanadham, D. Mechanical and electronic approaches to improve the sensitivity of microcantilever sensors. Acta Mech. Sin 2009, 25, 1–12. [Google Scholar]
  13. Yu, H; Li, X. Resonant-cantilever bio/chemical sensors with an integrated heater for both resonance exciting optimization and sensing repeatability enhancement. J. Micromech. Microeng 2009, 19, 045023. [Google Scholar]
  14. Rogers, B; Manning, L; Jones, M; Sulchek, T; Murray, K; Beneschott, B; Adames, JD; Hu, Z; Thundat, T; Cavazos, H; Minne, SC. Mercury vapor detection with a self-sensing, resonating piezoelectric cantilever. Rev. Sci. Instr 2003, 74, 4899–4901. [Google Scholar]
  15. Feng, G. Resonant sensing theory and devices; Tsinghua University Publishing House: Beijing, China, 2008; pp. 25–26. [Google Scholar]
  16. Gao, W. Study on resonant MEMS gas sensor based on thermally excited micro cantilever; M.S. Thesis; Tsinghua University: Beijing, China, 2008. [Google Scholar]
Figure 1. Electrical model of piezoelectric micro-cantilever.
Figure 1. Electrical model of piezoelectric micro-cantilever.
Sensors 10 07044f1
Figure 2. Impedance of the piezoelectric micro-cantilever.
Figure 2. Impedance of the piezoelectric micro-cantilever.
Sensors 10 07044f2
Figure 3. Frequency selection network module.
Figure 3. Frequency selection network module.
Sensors 10 07044f3
Figure 4. Schematic illustration of a SAD circuit.
Figure 4. Schematic illustration of a SAD circuit.
Sensors 10 07044f4
Figure 5. Simulation of phase-frequency curves of the SAD circuit loop, with AB and A’B’ the linear interval. Where f0 = 65,670 Hz and Δf0 = 200 Hz.
Figure 5. Simulation of phase-frequency curves of the SAD circuit loop, with AB and A’B’ the linear interval. Where f0 = 65,670 Hz and Δf0 = 200 Hz.
Sensors 10 07044f5
Figure 6. Scanning electron micrographs of the micro-cantilever: (a) Before the deposition of PEO sensitive material. (b) After the deposition of PEO sensitive material.
Figure 6. Scanning electron micrographs of the micro-cantilever: (a) Before the deposition of PEO sensitive material. (b) After the deposition of PEO sensitive material.
Sensors 10 07044f6
Figure 7. Schematic drawing of the experimental setup of the open-loop tests [3].
Figure 7. Schematic drawing of the experimental setup of the open-loop tests [3].
Sensors 10 07044f7
Figure 8. Typical raw data of the frequency shift of a SAD circuit loop and the internal micro-cantilever. To improve the visibility of the data points, the data are sampled with a time step of 40 seconds. The abscissa is time and each unit represents 40 seconds.
Figure 8. Typical raw data of the frequency shift of a SAD circuit loop and the internal micro-cantilever. To improve the visibility of the data points, the data are sampled with a time step of 40 seconds. The abscissa is time and each unit represents 40 seconds.
Sensors 10 07044f8
Figure 9. Frequency shift in response to the injection of ethanol. (a) MRF of the micro-cantilever in open-loop test [ Sensors 10 07044i1]. (b) SWF in close-loop test [ Sensors 10 07044i2]. (c) Simulation results of SWF [ Sensors 10 07044i3].
Figure 9. Frequency shift in response to the injection of ethanol. (a) MRF of the micro-cantilever in open-loop test [ Sensors 10 07044i1]. (b) SWF in close-loop test [ Sensors 10 07044i2]. (c) Simulation results of SWF [ Sensors 10 07044i3].
Sensors 10 07044f9
Figure 10. Experimental phase-frequency curves in compare with simulation results at MRF f0 = 65,670 Hz and SWF fr = 65,708 Hz. (a) Simulation results for SAD[ Sensors 10 07044i4], feedback module[ Sensors 10 07044i5] and frequency selection network[ Sensors 10 07044i6]. (b) Experimental measurements of SAD[ Sensors 10 07044i7], feedback module[ Sensors 10 07044i8] and frequency selection network[ Sensors 10 07044i9].
Figure 10. Experimental phase-frequency curves in compare with simulation results at MRF f0 = 65,670 Hz and SWF fr = 65,708 Hz. (a) Simulation results for SAD[ Sensors 10 07044i4], feedback module[ Sensors 10 07044i5] and frequency selection network[ Sensors 10 07044i6]. (b) Experimental measurements of SAD[ Sensors 10 07044i7], feedback module[ Sensors 10 07044i8] and frequency selection network[ Sensors 10 07044i9].
Sensors 10 07044f10

Share and Cite

MDPI and ACS Style

Guan, L.; Zhao, J.; Yu, S.; Li, P.; You, Z. Investigation of the Frequency Shift of a SAD Circuit Loop and the Internal Micro-Cantilever in a Gas Sensor. Sensors 2010, 10, 7044-7056. https://doi.org/10.3390/s100707044

AMA Style

Guan L, Zhao J, Yu S, Li P, You Z. Investigation of the Frequency Shift of a SAD Circuit Loop and the Internal Micro-Cantilever in a Gas Sensor. Sensors. 2010; 10(7):7044-7056. https://doi.org/10.3390/s100707044

Chicago/Turabian Style

Guan, Liu, Jiahao Zhao, Shijie Yu, Peng Li, and Zheng You. 2010. "Investigation of the Frequency Shift of a SAD Circuit Loop and the Internal Micro-Cantilever in a Gas Sensor" Sensors 10, no. 7: 7044-7056. https://doi.org/10.3390/s100707044

Article Metrics

Back to TopTop