Next Article in Journal
Design and Analysis of a Compact UWB Band Notch Antenna for Wireless Communication
Previous Article in Journal
Electronic Systems and Offsite Touristic Activities Based on Geological Concepts: A Speculative Discussion
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Proceeding Paper

A System-Level Modelling of Noise in Coupled Resonating MEMS Sensors †

by
Vinayak Pachkawade
Independent scientific consultant and a researcher, 4000 Liège, Belgium
Presented at the 1st International Electronic Conference—Futuristic Applications on Electronics, 1–30 November 2020; Available online: https://iec2020.sciforum.net/.
Published: 30 October 2020

Abstract

:
This paper presents realistic system-level modeling of effective noise sources in a coupled resonating mode-localized MEMS sensors. A governing set of differential equations are used to build a numerical model of a mechanical noise source in a coupled-resonator sensor and an effective thermo-mechanical noise is quantified through the simulation performed via SIMULINK. On a similar note, an effective noise that stems from the electronic readout used for the coupled resonating MEMS sensors is also quantified. Various noise sources in electronic readout are identified and the contribution of each is quantified. A comparison between an effective mechanical and electronic noise in a sensor system aids in identifying the dominant noise source in a sensor system. A method to optimize the system noise floor for an amplitude-based readout is presented. The proposed models present a variety of operating conditions, such as finite quality factor, varying coupled electrostatic spring strength, and operation with in-phase and out-of-phase mode. The proposed models aim to study the impact of fundamental noise processes that govern the ultimate resolution into a coupled resonating system used for various sensing applications.

1. Introduction

In a typical MEMS resonator system, electronic circuits (typically a current to voltage converter also called trans-impedance amplifier) is used for readout and processing of the electrical signal (output motional current from a resonator) provided by the sensing element. A total noise level of the system is therefore due to the combined effect of the mechanical–thermal noise in the mechanical domain, the electrical noise of the (resistive) mechanical sensing element [1] and the input-referred noise of the readout circuits [1,2,3,4]. In a MEMS sensor system, noise from a readout circuit dominates the noise performance [5]. However, the details of the mechanical–thermal noise should also be considered relevant, and a noise source (mechanical or electronic) that determines the detection limit of the sensor should be analyzed.
In this paper, an attempt is made to model and quantify an effective noise into the two coupled-resonator (WCR) system used for sensing applications. Through mathematical analysis and simulations obtained via the Simulink model, a source of the dominant noise source is identified.An effective noise floor in the design is quantified. Following this, the finest possible resolution for the design is proposed.

2. Noise Floor in a Coupled Resonator Sensor

In order to understand the impact of system noise on the performance of a coupled-resonator (CR) sensor, a noise source was added to the model. Output response was analyzed to determine the ultimate detection limit, i.e., resolution of a CR sensor.

Thermo-Mechanical Noise in CR Sensor Design

In this section, a model to quantify a magnitude of a thermal-mechanical noise is presented. Small moving parts in MEMS are especially susceptible to mechanical noise resulting from molecular agitation [1]. In ultra-low level signal detection, mechanical-thermal noise plays an important role in setting up the effective noise floor of a sensor system, and, thus a minimum detection limit. In the context of a CR sensor, it has been concluded that it is the thermal-mechanical noise that governs the ultimate detection limit of a CR sensor [2]. In this paper, a mechanical-thermal noise source in 2 DoF CR sensor is modeled and its impact on the resolution (i.e., lowest possible detectable physical quantity) is quantified for the amplitude-based output of the design.
Consider a schematic representation of a 2-DoF CR sensor system as shown in Figure 1a. It consists of effective mass, Mi, spring constant, Ki and damping constant, ci, (i = 1, 2). A coupling spring, Kc and damping constant, cc that may exist between the two resonators are also shown. Design conditions are: M1 = M2 = M, Km1 = Km2 = Km, c1 = c2 = cc = c and Kc << Km. In the design, it is stated that Xji is the amplitude variations of a j-th resonator (j = 1, 2) at the i-th mode of frequency response (i = 1, 2) due to the noise-induced into the system. Therefore, an expression of the output/s (considering a system involving an intrinsic (mechanical) and extrinsic (electronic) noise sources) may be obtained as follows: (i) a minimum resolvable amplitude ratio (AR) shift is ( X 1 X 2 ) min j i and (ii) a minimum resolvable amplitude shift is ( X ) min j i . A power spectral density (PSD) of an amplitude displacement noise of j-th resonator (j = 1, 2) at the i-th mode of the frequency response (i = 1, 2) can be investigated (through the Simulink model). An effect of a mechanical noise was modeled by adding a force term into the governing set of equations of motion (for 2 DoF) as follows
M X ¨ 1 + ( c + c c ) X ˙ 1 + ( K m + K c ) X 1 c X ˙ 2 K c X 2 = F n o i s e _ r m s
M X ¨ 2 + ( c + c c ) X ˙ 2 + ( K m + K c ) X 2 c X ˙ 1 K c X 1 = F n o i s e _ r m s
A spectral density of a noise forcing term is given as F n o i s e _ d e n s i t y = 4   k B T   c N/Hz0.5 [1,5]. Here, kB is the Boltzmann constant (≈1.380 × 10−23 Joule/Kelvin), T is the temperature (300 Kelvin) and c is the damping coefficient (c = 0.0031 Ns/m in our case). As observed, the spectral density of a mechanical noise force depends on temperature and the magnitude of mechanical damping. A noise forcing term with the calculated average value, 5.136 × 10−22 N2 was added as an excitation force into the Simulink model shown in Figure 1b.

3. Results

A noise PSD (in dB/Hz) for displacements X1 and X2 was plotted as shown in Figure 2. Simulations were run for the varying strength of a coupling spring, Kc between the two resonators in our design. Condition for the simulation are as follows: c = 0.0031 Ns/m, (Q ≈ 2547), ∆k = 0, κ = −0.0032, excitation force applied to mass M1. ∆k is the perturbation into the stiffness of one of the resonators in 2-DoF CR sensor system.
A resultant displacement noise PSD (due to noise forcing term) is −582.2 dB/Hz (at mode 1) and −572.2 dB/Hz (at mode 2) (refer to Figure 2a) which correspond to an equivalent magnitude | X j | of 7.76 × 10−30 m2/Hz (at mode 1) and 2.45 × 10−29 m2/Hz (at mode 2) respectively. Assuming a measurement bandwidth of 10 Hz around the resonant mode frequencies, ωi (i = 1, 2), an average value of a mechanical thermal noise of jth resonator at the ith mode (j, i = 1, 2) can be computed as X j i _ a v g (refer Table A1). All tables are placed in the appendix section.
Theoretical estimation of effective noise current for our design (due to intrinsic mechanical-thermal noise) is given as i m o t X j i = η ω i X j i , where, η   ω i , and Xji are transduction factor, angular frequency and the maximum displacement amplitude of the jth resonator (j = 1, 2) at the ith mode of the frequency response (i = 1, 2), respectively [6]. By determining the X j i the effective mechanical-thermal noise current (i.e., corresponding variations in the motional currents, i m o t X j i due to noise forcing term) was quantified (refer Table A1). Since the thermal noise amplitude of mode 1 is relatively lower (refer to Figure 2), the computation for the best case is given.

3.1. Modelling the Noise from the Electronic Readout

Figure 3a represents the schematic representation of transimpedance amplifier that may be realized using OPA 381 [7] integrated circuit (IC) as an interface circuit to the CR design. Two such ICs may be deployed for motional current pick-up into each output channel of the CR output. Figure 3b is the equivalent noise model for the interface circuit used (transimpedance amplifier). A noise contribution from each of these sources is estimated by superposition and nodal analysis.
An amplifier (OPA 381) [7] used in the analysis has an input current noise density, in ≈ 20 fA/Hz0.5, and an input voltage noise density, vn ≈ 70 nV/Hz0.5. In a CR sensor system, in the presence of an applied ac signal vac, an impact of noise on the motional current output i m o t j i of j-th resonator (j = 1, 2) at the i-th mode of response (i = 1, 2) can be evaluated (refer Figure 1b). This ac voltage signal together with the DC voltage, Vdc exerts a net electrostatic actuation force at the resonant frequency. In order to evaluate the impact of the electronic preamplifier (OPA 381) noise on the minimum shifts into the amplitude-based output/s (and thus quantify measurand, ∆k or δ k min (normalized)), all noise components were re-expressed as equivalent current as given in Table A2. In Table A2, i j n o i s e (rms) is the noise current from the j-th resonator at the i-th mode of response. Rx is the motional resistance, R x = d 4 k e f f M V d c 2 ε 2 A 2 Q [8] (Rx = 4 MΩ for a Q ≈ 2547 with other parameters being constant, in our case), Rf is the effective value of a feedback resistor used in the preamplifier (1 MΩ). The term kB is the Boltzmann constant (≈1.380 × 10−23 Joule/Kelvin) and T is the temperature (300 Kelvin). Term B is the integration bandwidth of 10 Hz around the resonator’s mode frequencies, f_ip and f_op) taken for computations for optimum noise estimation. From Table A2, an effective theoretical noise floor (i.e., input-referred electronic noise current, i j n o i s e due to electronic readout) of our design is ≈1.56 × 10−13 Arms. A total noise (mechanical resonator + electronics) in our design is then estimated by the vector sum of the uncorrelated noise sources i.e., It2 = I12 + I22, where, I1 is effective mechanical noise current and I2 is effective electronic readout noise current. It is calculated to be ≈1.56 × 10−13 Arms.
An electrostatic forcing term was applied to the Simulink model shown in Figure 1b. An effective system noise was added to the Simulink model. Other model conditions are the same as given in Figure 2. As a result, the following computations are obtained (refer to Table A3). Calculated values for motional current amplitudes are in agreement with the simulated values.

3.2. Resolving the Lowest Possible Shifts (i.e., Resolution)

The minimum resolvable shift in the voltage amplitudes of our 2-DoF CR sensor was derived as V j i n o i s e V o u t j i [6]. Here, V j i n o i s e is the output refereed noise voltage of the j-th coupled-resonator at the i-th mode. V o u t j i is the noiseless deterministic output voltages (determined as V o u t j i = i m o t j i × R f , where, i m o t j i = η ω i x j i ). Therefore, for the j-th resonator lowest possible resolvable mode amplitude shift, V j 1 n o i s e V o u t j 1 is 156 × 10 9 0.135 1.15 × 10 6 , (dimensionless) for mode 1, and V j 2 n o i s e V o u t j 2 is 156 × 10 9 0.271 5.756 × 10 7 (dimensionless) for mode 2. Similarly, the minimum resolvable shift in the amplitude ratio readout, ARi at the i-th mode (i = 1, 2) is ( V 11 n o i s e V o u t 11 ) 2 + ( V 21 n o i s e V o u t 21 ) 2 1.481 × 10 6 (dimensionless) for mode 1 and ( V 12 n o i s e V o u t 12 ) 2 + ( V 22 n o i s e V o u t 22 ) 2 3.89 × 10 7 (dimensionless) for mode 2. Since rms amplitudes of output voltage (essentially a motional current) at the mode 1 for j-th resonator is relatively closer (worst-case signal-to-noise (S/N) ratio) to the rms noise amplitude, it gives the possibility to determine the worst-case lowest possible shift (thus resolution) in the design. The effective resolution (theoretical) of our design for AR-based readout is ≈3.89 × 10−7 which may be compared to the theoretically calculated resolution (≈6 × 10−3) in MEMS coupled-resonator design in [6].

3.3. Influence of Coupling Strength on the Effective Noise Floor

This section shows if electrostatic spring strength, Kc, coupling the two resonators influence the effective noise floor. Figure 4 shows a simulated frequency spectrum of the j-th resonator (j = 1, 2) resonators for two variants of coupling spring, Kc. As observed, with the given set of operating conditions, the output signal from resonator 2 offers a relatively lower noise floor as compared to resonator 1 and is independent of the coupling strength (refer to Figure 4a,b). In addition, it is useful to utilize the output of the j-th resonator at mode 2 as it offers a higher amplitude and thus maximum shifts (sensitivity). From Figure 4, it is inferred that smaller coupling strength, Kc leads to reduced noise floor (about −12 dB gain in the noise floor thus enhancing the ultimate detection limit of the measurand) for the resonating output of the resonator 2.

4. Discussion and Conclusions

This work attempted to model and quantify the noise in the design of a 2-DoF MEMS CR sensor. An expected variation in the sensor output arising from the intrinsic (mechanical) and extrinsic (electronic) noise sources were individually evaluated via theory and system-level simulations. From Table A2, it is seen that though mechanical-thermal noise sets the ultimate noise floor in this design, it is an electronic noise that dominates the overall performance (electronics noise exceeds mechanical noise by two orders of magnitude). Through quantifying the system noise in the system, a system-level model aids in simulating the sensor output response. The ultimate detection limit in the 2-DoF CR sensor design is determined. A method to optimize the system noise/improve detection limit is also proposed (via tuning the coupling spring).

Appendix A

Table A1. Table showing the computation of the effective noise forcing term, displacement noise amplitudes and the resultant noise current due to the mechanical noise forcing term.
Table A1. Table showing the computation of the effective noise forcing term, displacement noise amplitudes and the resultant noise current due to the mechanical noise forcing term.
TermValueExpression
X j 1 _ a v g 7.762 × 10 29 m 2 X j i _ a v g = | X j i | × d f  1
X j i _ r m s 8.81 × 10 15 m X j 1 _ r m s = X j 1 _ a v g  2
i m o t X j 1 4.29 × 10 15 A r m s i m o t X j i = η ω i X j i  3
F n o i s e _ P S D 5.136 × 10 23 N 2 / H z ( F n o i s e _ d e n s i t y ) 2 = 4   k B T   c *
F n o i s e _ a v g 5.136 × 10 22 N 2 F n o i s e _ a v g = 0 10 ( 4   k B T   c ) 2 d f **
F n o i s e _ r m s 2.26 × 10 11 N F n o i s e _ r m s = F n o i s e _ a v g ***
1 A best-case estimate for average noise power for mode 1. 2 An effective rms value of a mechanical-thermal noise for mode 1. 3 An effective noise current due to the mechanical noise forcing term. * Power spectral density (PSD) of mechanical noise force generator. ** Average (mean square) value of a mechanical thermal noise force generator, assumed measurement bandwidth is 10 Hz. *** A effective rms amplitude of mechanical noise forcing term (mean square) value of mechanical thermal noise of j-th resonator at the i-th mode (j, i = 1, 2) will be X j i _ a v g = | X j i | × d f . Refer Table A1 for the computations. As observed from Figure 2, noise magnitude is lower for mode 1.
Table A2. Table showing the contribution of the noise current due to the several noise terms used in the electronic readout for the mechanical CR sensor.
Table A2. Table showing the contribution of the noise current due to the several noise terms used in the electronic readout for the mechanical CR sensor.
Noise SourcesEffective Noise Current (Arms)Expression
Mechanical thermal-noise 14.9 × 10−15 i X j i = η ω i X j i
Feedback resistance4.06 × 10−13 4 k B T × B R f
Input voltage noise of a pre-amplifier4.34 × 10−13 v n 2 ( 1 + R x R f ) 2 / R x 2 × B × 1.57
Input current noise of a pre-amplifier 29.92 × 10−14 i n 2 × B × 1.57
input-referred electronic noise current, i j n o i s e 1.56 × 10−13 Arms i j n o i s e = ( i j n o i s e ) 2 = [ i n 2 + v n 2 ( 1 + R x R f ) 2 / R x 2 + 4 k B T R f ]
Total system noise (mechanical + electronic)≈1.56 × 10−13 ArmsIt2 = I12 + I22
1 Best-case calculation for the out-of-phase mode as it shows the lower noise amplitude. 2 Factor 1.57 is the roll-off rate of a filter (1-pole) [9].
Table A3. Computations and simulated results in response to the electrostatic forcing term used in the Simulink model. Noise due to mechanical and electronic terms was taken into account to simulate the output response.
Table A3. Computations and simulated results in response to the electrostatic forcing term used in the Simulink model. Noise due to mechanical and electronic terms was taken into account to simulate the output response.
TermValue (Mode 1)Value (Mode 2)Expression
maximum displacement 1, | x j i | xj1 ≈ 0.419 µm xj2 ≈ 0.836 µm N/A
motional current amplitudes 2, i m o t j i ≈ 192 nA≈ 384 nA i m o t j 1 η × 2 π f _ o p × x j 1 , i m o t j 2 η × 2 π f _ i p × x j 2
Effective output voltage 3, V o u t j i = 0.135 Vrms= 0.271 Vrms V o u t j 1 = i m o t j 1 × R f 0.192 V max , V o u t j 2 = i m o t j 2 × R f 0.384 V max
output referred noise voltage 2, V j i n o i s e ≈156 × 10−9 VrmsN/A i j 1 n o i s e × Rf
lowest possible resolvable mode amplitude shift, V j i n o i s e V o u t j i 156 × 10 9 0.135 1.15 × 10 6 (dimensionless) 156 × 10 9 0.271 5.756 × 10 7 (dimensionless) V j 1 n o i s e V o u t j 1 , V j 2 n o i s e V o u t j 2
minimum resolvable shift in the amplitude ratio readout, ARi 1.481 × 10 6 (dimensionless) 3.89 × 10 7 (dimensionless) ( V 11 n o i s e V o u t 11 ) 2 + ( V 21 n o i s e V o u t 21 ) 2 , ( V 12 n o i s e V o u t 12 ) 2 + ( V 22 n o i s e V o u t 22 ) 2
1 Derived from the Simulink. 2 Theory. 3 Amplification factor, Rf = 1 MΩ, as a feedback resistor in preamplifier OPA 381.

References

  1. Gabrielson, T.B. Mechanical-Thermal Noise in Micromachined Acoustic and Vibration Sensors. IEEE Trans. Electron. Devices 1993, 40, 903–909. [Google Scholar] [CrossRef]
  2. Juillard, J.; Prache, P.; Ferreira, P.M.; Barniol, N. Ultimate limits of differential resonant MEMS sensors based on two coupled linear resonators. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2018, 65, 2440–2448. [Google Scholar] [CrossRef] [PubMed]
  3. Juillard, J.; Prache, P.; Ferreira, P.M.; Barniol, N. Impact of output metric on the resolution of mode-localized MEMS resonant sensors. In Proceedings of the 2017 Joint Conference of the European Frequency and Time Forum and IEEE International Frequency Control Symposium (EFTF/IFCS), Besancon, France, 9–13 July 2017. [Google Scholar] [CrossRef]
  4. Pandit, M.; Zhao, C.; Sobreviela, G.; Mustafazade, A.; Du, S.; Zou, X.; Seshia, A.A. Closed-Loop Characterization of Noise and Stability in a Mode-Localized Resonant MEMS Sensor. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2019. [Google Scholar] [CrossRef] [PubMed]
  5. Djurić, Z. Mechanisms of noise sources in microelectromechanical systems. Microelectron. Reliab. 2002, 40, 919–932. [Google Scholar] [CrossRef]
  6. Thiruvenkatanathan, P.; Woodhouse, J.; Yan, J.; Seshia, A.A. Limits to mode-localized sensing using micro- and nanomechanical resonator arrays. J. Appl. Phys. 2011, 109, 104903. [Google Scholar] [CrossRef]
  7. Precision, Low Power, 18 MHz Transimpedance Amplifier. Available online: www.ti.com (accessed on 23 June 2019).
  8. Hung, L.; Nguyen, C. High-Q Low-Impedance MEMS Resonators; Technical Report No. UCB/EECS-2012-218; University of California: Berkeley, CA, USA, December 2012. [Google Scholar]
  9. TI Precision Labs—Op Amps: Noise 1|TI.com Video. Available online: https://training.ti.com/ti-precision-labs-op-amps-noise-1 (accessed on 23 June 2019).
Figure 1. (a) Schematic representation of a 2-DoF coupled-resonator (CR) sensor system. Proof mass, Mi, springs, Ki and dampers, ci (i = 1, 2) are shown (b) Simulink block diagram representing the system-level model.
Figure 1. (a) Schematic representation of a 2-DoF coupled-resonator (CR) sensor system. Proof mass, Mi, springs, Ki and dampers, ci (i = 1, 2) are shown (b) Simulink block diagram representing the system-level model.
Engproc 03 00005 g001
Figure 2. Simulated power spectrum density (PSD) of a displacement noise of resonator 1 and 2 subject to mechanical-thermal noise forcing term on resonator 1. (a) Kc = −393.5 N/m and (b) Kc = −1000 N/m. Simulation conditions are as follows: c = 0.0031 Ns/m, (Q ≈ 2547), ∆k = 0, excitation force applied to mass M1. A normalized coupling factor, κ is given by Kc/Keff. ∆k is the induced perturbation into the stiffness of one of the resonators in a 2-DoF CR sensor system.
Figure 2. Simulated power spectrum density (PSD) of a displacement noise of resonator 1 and 2 subject to mechanical-thermal noise forcing term on resonator 1. (a) Kc = −393.5 N/m and (b) Kc = −1000 N/m. Simulation conditions are as follows: c = 0.0031 Ns/m, (Q ≈ 2547), ∆k = 0, excitation force applied to mass M1. A normalized coupling factor, κ is given by Kc/Keff. ∆k is the induced perturbation into the stiffness of one of the resonators in a 2-DoF CR sensor system.
Engproc 03 00005 g002
Figure 3. (a) OPA 381 interface for a 2-DoF CR sensor design (b) equivalent noise circuit to evaluate circuit dominant noise.
Figure 3. (a) OPA 381 interface for a 2-DoF CR sensor design (b) equivalent noise circuit to evaluate circuit dominant noise.
Engproc 03 00005 g003
Figure 4. Simulated PSD of the motional current output signal in two weakly coupled-resonators for the following operating condition, Q = 2547, δk = 0, F = 1 N.
Figure 4. Simulated PSD of the motional current output signal in two weakly coupled-resonators for the following operating condition, Q = 2547, δk = 0, F = 1 N.
Engproc 03 00005 g004
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Pachkawade, V. A System-Level Modelling of Noise in Coupled Resonating MEMS Sensors. Eng. Proc. 2020, 3, 5. https://0-doi-org.brum.beds.ac.uk/10.3390/IEC2020-06971

AMA Style

Pachkawade V. A System-Level Modelling of Noise in Coupled Resonating MEMS Sensors. Engineering Proceedings. 2020; 3(1):5. https://0-doi-org.brum.beds.ac.uk/10.3390/IEC2020-06971

Chicago/Turabian Style

Pachkawade, Vinayak. 2020. "A System-Level Modelling of Noise in Coupled Resonating MEMS Sensors" Engineering Proceedings 3, no. 1: 5. https://0-doi-org.brum.beds.ac.uk/10.3390/IEC2020-06971

Article Metrics

Back to TopTop