Self-Calibration and Performance Control of MEMS with Applications for IoT

A systemic problem for microelectromechanical systems (MEMS) has been the large gap between their predicted and actual performances. Due to process variations, no two MEMS have been able to perform identically. In-factory calibration is often required, which can represent as much as three-fourths of the manufacturing costs. Such issues are challenges for microsensors that require higher accuracy and lower cost. Towards addressing these issues, this paper describes how microscale attributes may be used to enable MEMS to accurately calibrate themselves without external references, or enable actual devices to match their predicted performances. Previously, we validated how MEMS with comb drives can be used to autonomously self-measure their change in geometry in going from layout to manufactured, and we verified how MEMS can be made to increase or decrease their effective mass, damping, and or stiffness in real-time to match desired specifications. Here, we present how self-calibration and performance control may be used to accurately sense and extend the capabilities of a variety of sensing applications for the Internet of things (IoT). Discussions of IoT applications include: (1) measuring absolute temperature due to thermally-induced vibrations; (2) measuring the stiffness of atomic force microscope or biosensor cantilevers; (3) MEMS weighing scales; (4) MEMS gravimeters and altimeters; (5) inertial measurement units that can measure all four non-inertial forces; (6) self-calibrating implantable pressure sensors; (7) diagnostic chips for quality control; (8) closing the gap from experiment to simulation; (9) control of the value of resonance frequency to counter drift or to match modes; (10) control of the value of the quality factor; and (11) low-amplitude Duffing nonlinearity for wideband high-Q resonance.


Introduction
Attributes of microelectromechanical systems (MEMS) include a wide variety of transduction capabilities in packages that are small in size, weight, power, and cost. Such attributes facilitate the Internet of Things (IoT), where physically-sensed information about the environment is passed between objects (things) and the Internet. In addition to monitoring the performance of critical processes, MEMS that are widely and ubiquitously distributed will enable the Internet to sense the physical world. Examples of IoT application areas include: (1) society-health care, media and entertainment, smart environments, intelligent transportation, smart retail, security and surveillance; (2) environment-smart agriculture, disaster management, pollution control, and smart power plants; and (3) industry-supply chain management, aerospace and aviation, transportation and logistics, smart metering, warehouse and storage [1][2][3][4]. IoT will generate an enormous amount of data. The quality of sensed data will be based on the performance of MEMS. Although MEMS are very precise, they have accuracy issues.
There are many IoT applications that do to not require MEMS sensors with high accuracy, precision, or performance control, but for present or future applications that cannot afford wide margins of uncertainty, then highly accurate, precise, and extended behavioral control may make

Self-Calibration
This section describes how packaged and deployed MEMS devices with comb drives can be made to accurately measure their own mass, damping, stiffness, displacement, and force without the need for in-factory calibration or external references. The extension of this metrology to other types of MEMS sensors is discussed below. The method is accurate, repeatable, and reliable. Here, accuracy is a measure of the difference between the average value and the true value; uncertainty is the order of the most significant uncertain digit due to the totality of measurement noise; repeatability is the ability of a measurement method to obtain the same value after reassembling the experiment; and reliability is the ability of the same measurement method to achieve the same results by using different testing equipment (i.e., test equipment by a different manufacturer). We define verification as a test of how well different methods of analyses agree with each other (such as analytical theory versus simulation); and we define validation as a test of how well a measurement method agrees with the true value.
Validation. The lack of accurate and traceable measurement standards for MEMS makes direct validation difficult because the true value is unknown. However, a condition that must hold true in the validation of any metrological method is that two or more sensors that use that method must agree on their measured values whenever identical disturbances are applied. For instance, although the amount of an applied force may be unknown, two different sensors that use the same metrological method to measure that force should agree on their measurements of that force.
Assumptions. We assume that: (1) Geometric features that are within close-proximity to each other undergo identical process variations. This assumption is, of course, measurable through validation tests. (2) The following analysis requires that comb drives operate within their linear range of motion.
Sequence of measurements. We previously presented the electrically-probed measurement of the gap between structures and flexure width. The measurements were validated against scanning electron microscopy [21,22]. In either case, a relationship was necessarily identified that expressed the mechanical quantity as a function of electrical measurands only. Doing so enables measurements to be performed completely within a packaged chip. Since most MEMS relations are coupled to each other by one or more other mechanical quantities, then a sequence of quantity extraction exists. That is, relationships between unknown mechanical quantities m i (such as stiffness, width, mass, gap, damping, force, Young's modulus, density, etc.) and measurable electrical quantities e i (such as capacitance, voltage, frequency, etc.) can be expressed by the following electro-micro-metrology (EMM) functions G i : m 1 = G * 1 (e 1 ) m 2 = G 2 (m 1 , e 2 ) = G 1 (e 1 , e 2 ) m 3 = G 3 (m 2 , e 3 ) = G 2 (e 1 , e 2 , e 3 ) . . . m N = G N (m N−1 , e N ) = G N (e 1 , . . . , e N ) (1) where the first seed relation G * 1 is special in that it is not a function of any mechanical quantity. All subsequent measurements of mechanical quantities m i , i > 1, depend on one or more prior measurements m j , j < i. However, each m i quantity is ultimately a function of electrical measurands only. In all, relations in Equation (1) form a ordered sequence of quantity extractions. Due to interdependency between mechanical quantities, a large number of mechanical quantities can be determined from a much smaller number of electrical measurands.
Quantifiable uncertainty. The uncertainty of an electrical measurand, δe j , is easily determined by, for example, the order of the most significant uncertain digit on its readout meter. Subsequently, the uncertainty of the mechanical quantity, δm i , of the MEMS device can be determined by a 1st-order multivariate Taylor expansion about e j as shown in Equation (2).
In Equation (2), the small-valued electrical uncertainties δe j are multiplied by the large-valued electromechanical sensitivities ∂G i /δe j . It is not uncommon to have uncertainties on the order of O 10 −18 with corresponding sensitivities on the order of O 10 12 . As length scale decreases, so too does the sensitivity, which is a reason why this metrology method becomes tractable at the microscale and below: Example. Let's apply the above EMM concepts to the type of MEMS shown in Figure 1. The device consists of two pairs of comb drives for actuation and sensing, a pair of asymmetric gaps, and a pair of folded flexure spring supports. Upon fabrication, packaging, and deployment, the device would be subjected to process variations, packaging stress, and environmental changes, yielding the following unknown properties: Young's modulus, Poisson's ratio, viscosity, geometric overetch, layer thickness, material density, permittivity, fringing field factor, curvature of radius for the fillets located at all vertices, gaps between comb fingers, etc. As a consequence, the quantities that depend on such properties are unknown as well. This includes mass, damping, stiffness, force, displacement, etc. Simulated illustration of an electrically-probed calibration of a MEMS device. One probe applies a voltage V1 to close one of two asymmetric gaps, while the differential capacitance to traverse gap1 is measured. Many other aspects of the design, such as the type flexures, are irrelevant. By closing the gaps and resonating the structure, the device's mass, damping, stiffness, and state can be accurately and precisely measured after packaging and long-term dormancy. This is, it is capable of self-calibration.
The first step is to find an EMM expression * 1 G that relates an unknown mechanical quantity to electrical quantities. For our testcase, we choose our first seed relation to be overetch. A variety of methods can be used to close gaps, depending on gap size, such as DC voltage, resonance, mechanical force, etc. Overetch can then be used to determine displacement (G2), which will be used to determine force (G3), which will be used to determine stiffness (G4), and so on. As derived in Reference [2], geometric overcut can be determined by measuring the change in capacitances ( 1 C  and 2 C  ) required to close the manufactured gaps (gap1 and gap2 in Figure 1). That is:    cancels out the parasitic capacitance 1 P C , which enables this method to be reliable because parasitic capacitance will be different for each comb drive and between each test setup (e.g., dielectric charging, probe contact area, cable orientation, etc.). The welldefined gap-stops enable the method to be repeatable. And utilizing precise electrical measurands enables high accuracy and low uncertainty. For instance, Analog Devices (ADI, Norwood, MA, USA) reports to have sensed an average comb drive displacement of 100 femtometers (a thousandth of the diameter of a hydrogen atom) due to a change in comb drive capacitance of a zeptofarad [23]. Such a sensitivity ratio (comb drive displacement to change in capacitance) agrees with our results below. The result from Equation (3) can be used to measure the next mechanical quantity, displacement. Since the comb drives will be operating within their linear range, then the ratio of change in capacitance to the distance traversed is the comb drive constant 1 , the deflection can be measured as: Figure 1. Simulated illustration of an electrically-probed calibration of a MEMS device. One probe applies a voltage V 1 to close one of two asymmetric gaps, while the differential capacitance to traverse gap 1 is measured. Many other aspects of the design, such as the type flexures, are irrelevant. By closing the gaps and resonating the structure, the device's mass, damping, stiffness, and state can be accurately and precisely measured after packaging and long-term dormancy. This is, it is capable of self-calibration.
The first step is to find an EMM expression G * 1 that relates an unknown mechanical quantity to electrical quantities. For our testcase, we choose our first seed relation to be overetch. A variety of methods can be used to close gaps, depending on gap size, such as DC voltage, resonance, mechanical force, etc. Overetch can then be used to determine displacement (G 2 ), which will be used to determine force (G 3 ), which will be used to determine stiffness (G 4 ), and so on. As derived in Reference [2], geometric overcut can be determined by measuring the change in capacitances (∆C 1 and ∆C 2 ) required to close the manufactured gaps (gap 1 and gap 2 in Figure 1). That is: where ∆gap [m] is overcut, i.e., the difference in going from layout gap gap 1,Layout to manufactured gap gap 1 ; and n = 1 is the ratio between the gaps, i.e., gap 2,Layout = n gap 1,Layout . This first measurement of a mechanical quantity is a function of electrical measurements and exactly-known design parameters.
It is important to note that the change in capacitance ∆C 1 = C comb cancels out the parasitic capacitance C P 1 , which enables this method to be reliable because parasitic capacitance will be different for each comb drive and between each test setup (e.g., dielectric charging, probe contact area, cable orientation, etc.). The well-defined gap-stops enable the method to be repeatable. And utilizing precise electrical measurands enables high accuracy and low uncertainty. For instance, Analog Devices (ADI, Norwood, MA, USA) reports to have sensed an average comb drive displacement of 100 femtometers (a thousandth of the diameter of a hydrogen atom) due to a change in comb drive capacitance of a zeptofarad [23]. Such a sensitivity ratio (comb drive displacement to change in capacitance) agrees with our results below.
The result from Equation (3) can be used to measure the next mechanical quantity, displacement. Since the comb drives will be operating within their linear range, then the ratio of change in capacitance to the distance traversed is the comb drive constant Ψ 1 = ∆C 1 /gap 1 [F/m]. Therefore, for any intermediate deflection, x < gap 1 , that produces a change in capacitance ∆C(x), the deflection can be measured as: The linear operating range also permits us to equate ∂C 1 /∂x= ∆C 1 /∆x = Ψ 1 . Therefore, the drive force, generated by an applied voltage V 1 , which includes variable comb finger gaps and fringing fields, is: The ratio of the change in force to change in displacement is a measure of stiffness, which includes nonlinearities, defects, fillets, non-rectangular cross sections, course sidewalls, variable widths, variable Young's modulus, anchor compliance, etc. The stiffness of the system is: Once stiffness is known, a measure of system mass can be obtained by measuring undamped displacement resonance frequency. However, since a high vacuum is necessary for undamped displacement-resonance frequency, a much more practical measure can be done by measuring the velocity-resonance frequency ω .
x max = ω 0 = √ K/M, which is independent of damping, unlike displacement-resonance. Since mass, stiffness, and velocity-resonance is related by M = K/ω 2 0 , from Equation (6) we have: Damping D is related to mass by D = 2Mγ 2 , where the exponential decay rate is γ = ω 2 0 − ω 2 r /2 and ω r is the damped displacement-resonance frequency. We therefore have: At this point in our sequence of quantity extractions, we now have quantities that define the microsystem's equation of motion along the x-direction: Continuing on, the sequence of quantity extractions leads to relations for displacement amplitude, quality factor, Young's modulus, material density, viscosity, layer thickness, etc. For example, amplitude and qualify factor are: and: where ω d = ω 2 0 + ω 2 r /2 is the exponentially decaying oscillation frequency without drive excitation (F dr = 0).
So instead of being limited to frequency response, this metrological method enables accurate measurements of mass, damping, stiffness, force, and state of the system. In particular, this facilitates accurate sensing of force and displacement due to external disturbances, which can be quite useful since nearly all physical quantities can be traced to relations that depend on force and displacement. As Equation (2) describes, measurements of uncertainties for each of the extracted quantities of (3) to (11) can be determined by a truncated multivariate Taylor expansion, where each electrical measurand is expressed in terms of its nominal and uncertain quantities; i.e., e measured → e ± δe . For instance, for capacitance and voltage, we have: and similarly: where it is assumed that δC and δV represent the totality of all noise, which can be conservatively determined by the decimal place of the most significant uncertain digit on, say, the capacitance or voltage meter's readout display. For example, the uncertainty expressions for displacement (4) and force (4) are: where the curly-bracketed expressions are the sensitivities. For the device in Figure 1, the sensitivity of δx is on the order of O 10 9 meters/Farad for the device in Figure 1, so an uncertainty for δx that would be on the order of an angstrom requires that δC = O 10 −10 [F].
Off-the-shelf capacitance meters can have precisions of δC = O 10 −12 F, or a pF [24][25][26]. The most precise capacitance meters to date have δC = O 10 −21 F, or zF [27,28]. Similarly, the sensitivities of δF are O 10 6 and O 10 −9 for the first and second terms in Equation (15). The order of precision for voltage source/meters typically range from millivolts to nanovolts. Such a range of precisions for capacitance and voltage are expected to result in a range of force uncertainties from micronewtons to femtonewtons. However, there are also other ways to improve precision. In Reference [21] it is shown that the order of EMM uncertainties also depends on design parameters, which affect the order of sensitivities for δx and δF. That is, there will be a decrease in sensitivity if there is a decrease in comb finger gap g f , or if there is an increase in the number of comb fingers N f , gap ratio n, layer thickness h, or gap-stop size gap 1 . So by choice of device design parameters and meter precisions, a desired precision magnitude may be achieved. Tangible examples of force magnitudes are provided in Table 1.

Performance Control
The previous section described how system mass, damping, stiffness, force, and state can be accurately and precisely measured using electro micro metrology (EMM), this section describes how the effective mass, damping, and stiffness of a MEMS device can be changed, such that the microsystem can modify its performance on demand (POD). EMM can be used to recalibrate to the modified effective mass, damping, and stiffness of a PODMEMS device.
With control over effective mass, damping, and stiffness, PODMEMS devices are expected to be able to compensate for performance variations due to processing, packaging, temperature, noise, damping, or extend dynamic range, or be used to accommodate multiple application modes. Such control is achieved by continuously monitoring the state of the proof mass, and feeding back forces onto the proof mass that are proportional to displacement, velocity, and or acceleration. A system-level view of the feedback for a PODMEMS device is shown in Figure 2. With the addition of feedback forces fb F , the equation of motion (9) becomes: where fb

Performance Control
The previous section described how system mass, damping, stiffness, force, and state can be accurately and precisely measured using electro micro metrology (EMM), this section describes how the effective mass, damping, and stiffness of a MEMS device can be changed, such that the microsystem can modify its performance on demand (POD). EMM can be used to recalibrate to the modified effective mass, damping, and stiffness of a PODMEMS device.
With control over effective mass, damping, and stiffness, PODMEMS devices are expected to be able to compensate for performance variations due to processing, packaging, temperature, noise, damping, or extend dynamic range, or be used to accommodate multiple application modes. Such control is achieved by continuously monitoring the state of the proof mass, and feeding back forces onto the proof mass that are proportional to displacement, velocity, and or acceleration. A system-level view of the feedback for a PODMEMS device is shown in Figure 2. With the addition of feedback forces fb F , the equation of motion (9) becomes: where fb

Performance Control
The previous section described how system mass, damping, stiffness, force, and state can be accurately and precisely measured using electro micro metrology (EMM), this section describes how the effective mass, damping, and stiffness of a MEMS device can be changed, such that the microsystem can modify its performance on demand (POD). EMM can be used to recalibrate to the modified effective mass, damping, and stiffness of a PODMEMS device.
With control over effective mass, damping, and stiffness, PODMEMS devices are expected to be able to compensate for performance variations due to processing, packaging, temperature, noise, damping, or extend dynamic range, or be used to accommodate multiple application modes. Such control is achieved by continuously monitoring the state of the proof mass, and feeding back forces onto the proof mass that are proportional to displacement, velocity, and or acceleration. A system-level view of the feedback for a PODMEMS device is shown in Figure 2. With the addition of feedback forces fb F , the equation of motion (9) becomes: where fb

Performance Control
The previous section described how system mass, damping, stiffness, force, and state can be accurately and precisely measured using electro micro metrology (EMM), this section describes how the effective mass, damping, and stiffness of a MEMS device can be changed, such that the microsystem can modify its performance on demand (POD). EMM can be used to recalibrate to the modified effective mass, damping, and stiffness of a PODMEMS device.
With control over effective mass, damping, and stiffness, PODMEMS devices are expected to be able to compensate for performance variations due to processing, packaging, temperature, noise, damping, or extend dynamic range, or be used to accommodate multiple application modes. Such control is achieved by continuously monitoring the state of the proof mass, and feeding back forces onto the proof mass that are proportional to displacement, velocity, and or acceleration. A system-level view of the feedback for a PODMEMS device is shown in Figure 2. With the addition of feedback forces fb F , the equation of motion (9) becomes: where fb

Performance Control
The previous section described how system mass, damping, stiffness, force, and state can be accurately and precisely measured using electro micro metrology (EMM), this section describes how the effective mass, damping, and stiffness of a MEMS device can be changed, such that the microsystem can modify its performance on demand (POD). EMM can be used to recalibrate to the modified effective mass, damping, and stiffness of a PODMEMS device.
With control over effective mass, damping, and stiffness, PODMEMS devices are expected to be able to compensate for performance variations due to processing, packaging, temperature, noise, damping, or extend dynamic range, or be used to accommodate multiple application modes. Such control is achieved by continuously monitoring the state of the proof mass, and feeding back forces onto the proof mass that are proportional to displacement, velocity, and or acceleration. A system-level view of the feedback for a PODMEMS device is shown in Figure 2. With the addition of feedback forces fb F , the equation of motion (9) becomes: where fb

Performance Control
The previous section described how system mass, damping, stiffness, force, and state can be accurately and precisely measured using electro micro metrology (EMM), this section describes how the effective mass, damping, and stiffness of a MEMS device can be changed, such that the microsystem can modify its performance on demand (POD). EMM can be used to recalibrate to the modified effective mass, damping, and stiffness of a PODMEMS device.
With control over effective mass, damping, and stiffness, PODMEMS devices are expected to be able to compensate for performance variations due to processing, packaging, temperature, noise, damping, or extend dynamic range, or be used to accommodate multiple application modes. Such control is achieved by continuously monitoring the state of the proof mass, and feeding back forces ..
x τ onto the proof mass that are proportional to displacement, velocity, and or acceleration.
A system-level view of the feedback for a PODMEMS device is shown in Figure 2. With the addition of feedback forces F f b , the equation of motion (9) becomes: .. where x τ , and F M = M e ..
x(t − τ) of the proof mass are subject to feedback delay τ due to the RC time constant of the analog feedback circuit. The quantities K e , D e , and M e are the electrically-generated proportionality constants that can increase or decrease the effective mass, damping, or stiffness system. That is, Equation (16) may be rewritten as: ..
x + D e .
x + D e f f .
x where the effective mass, damping, and stiffness may be expressed as: assuming low latency. For large latencies, exact expressions for M e f f , D e f f , and K e f f are derived in [29]. As seen in Equation (18), there is a small amount of crosstalk within the effective mass and damping, i.e., D e τ, K e τ, M e ω 2 τ, where the amount of crosstalk is reduced by feedback delay.
circuit. The quantities Ke, De, and Me are the electrically-generated proportionality constants that can increase or decrease the effective mass, damping, or stiffness system. That is, Equation (16) where the effective mass, damping, and stiffness may be expressed as:  [29]. As seen in Equation (18), there is a small amount of crosstalk within the effective mass and damping, i.e., e D  , e K  , 2 e M   , where the amount of crosstalk is reduced by feedback delay.
. Figure 2. MEMS + performance controller [29]. Symmetric feedback components on the right and left (not shown) sides of the structure operate 180° out of phase for continuous feedback response throughout each cycle. Feedback forces are applied to left and right combs according to position, velocity, and acceleration.
The dynamical range of PODMEMS is much greater than what can be achieved by MEMS that do not exploit feedback. Due to feedback latency, instability is also possible. The dynamic range and stability of a PODMEMS device can be visually determined by Figure 3a-d as follows (exact expressions for stability are derived in Reference [29]). The vertical and horizontal axes in Figure 3 are Ke/K and De/D. The family of curves (or lobes) indicate the bounds of instability, where PODMEMS is stable within the lobe, which is bounded between the curve and dashed line at Ke/K = −1. PODMEMS is unstable outside a given lobe. The family of lobes are parameterized by Me in Figure 3a  The dynamical range of PODMEMS is much greater than what can be achieved by MEMS that do not exploit feedback. Due to feedback latency, instability is also possible. The dynamic range and stability of a PODMEMS device can be visually determined by Figure 3a-d as follows (exact expressions for stability are derived in Reference [29]). The vertical and horizontal axes in Figure 3 are K e /K and D e /D. The family of curves (or lobes) indicate the bounds of instability, where PODMEMS is stable within the lobe, which is bounded between the curve and dashed line at K e /K = −1. PODMEMS is unstable outside a given lobe. The family of lobes are parameterized by M e in Figure 3a (2) as shown in Figure 3d, there are regions where the overall damping can be negative while stability is maintained; (3) the electrical mass, damping, and stiffness are allowed to be several orders greater in magnitude than their purely mechanical counterparts; and (4) due to the small length scale and time scale of MEMS, the domain of stability for any PODMEMS is very sensitive to feedback latency.  (17), the characteristic effective exponential decay rate eff  , quality factor eff Q , amplitude at displacement resonance max x , and velocity resonance 0  are:

With a system described by Equation
, where, instead of such characteristic behaviors being constant as they are for passive MEMS, such characteristics for PODMEMS are greatly modifiable, controlled by electrical feedback forces that depend on real-time state monitoring.
Let's apply the above PODMEMS concepts to the MEMS device shown in Figure 2. The design parameters for the device are: pairs of folded flexures having width w = 2 µm, thickness h = 20 µm, and length L = 294.7 µm; pairs of 100-finger comb drive arrays with finger length Lf = 20 µm, finger width  (17), the characteristic effective exponential decay rate γ e f f , quality factor Q e f f , amplitude at displacement resonance x max , and velocity resonance ω 0 are:

With a system described by Equation
Q e f f = ω d,e f f 2γ e f f (20)
When designing feedback control systems, it can be useful to express maximal mechanical forces in terms of equivalently-applied voltages; i.e., the voltage necessary to generate the equivalent feedback force. For simple harmonic motion, we have: where the applied control voltages V K , V D , and V M , correspond to feedback forces F K , F D , and F M . For example, using the PODMEMS device from Figure 2 where ω r = 50 krad/ sec and Ψ ≈

Possible IoT Applications
This section envisions a few possible applications of self-calibration and performance control that may be useful for IoT and other application areas.

Metrology
The benefits of self-calibration are expected to: (1) reduce the cost of MEMS devices since the costly expense of in-factory calibration can be reduced or eliminated, which should also increase manufacturing throughput; (2) greatly extend the usefulness of sensors since devices will be able to re-calibrate after long-term dormancy or after harsh environmental change; (3) improve the quality of data being analyzed in terms of increased accuracy and reduced uncertainty; and (4) close the large gap in going from experiment back to simulation, to build experimentally-accurate predictive computer models of IoT sensors and of their environments. Examples follow: Temperature sensors. Accurate measurements of temperature can be useful for improved: weather predictions, ecological health monitoring, infrared sensing, sensing of factory equipment or chemical processes, implantable health monitoring, continuous measurement of engine efficiency, heat energy, pressure, volume, thermal efficiency, coefficient of performance, entropy production, etc. Due to the length scale, most microdevices are subject to measurable thermally-induced vibrations. From the equipartition theorem, mechanical potential energy is related to thermal energy by 1 2 K x 2 = 1 2 k B T [30] where k B is the Boltzmann constant, T is absolute temperature, and x 2 is the mean-square of displacement due to thermally-induced vibrations, and other sources. Since EMM enables accurate measurements of flexure stiffness K and displacement x, then by the equipartition theorem, the absolute temperature T of the MEMS device can be accurately measured.
Atomic force microscopy and biosensors. Once a MEMS device is accurately calibrated, it becomes an accurate sensor for calibrating other devices, such as cantilevers for atomic force microscope (AFM), cantilevers for biosensors, the self-calibration of AFM-on-chips, mass sensors, etc. See Table 1. Due to the inaccuracies and large uncertainties (>10%) of conventional cantilever stiffness measurement methods [30][31][32][33][34], it is difficult for biotechnologists to discern targeted from nonspecific bindings, and most AFM users do not use the AFM to measure force [35][36][37]. However, by pressing the tip of a cantilever against the sidewall of a calibrated MEMS device (Figure 4), the MEMS device can apply a known force (5) to deflect the cantilever by measured amount (4). The force balance relation would be F = (K + K AFM )x, where K AFM is the unknown AFM cantilever stiffness to be determined. Since x AFM = x, then the photodiode can be simultaneously calibrated. Atomic force microscopy and biosensors. Once a MEMS device is accurately calibrated, it becomes an accurate sensor for calibrating other devices, such as cantilevers for atomic force microscope (AFM), cantilevers for biosensors, the self-calibration of AFM-on-chips, mass sensors, etc. See Table  1. Due to the inaccuracies and large uncertainties (>10%) of conventional cantilever stiffness measurement methods [30][31][32][33][34], it is difficult for biotechnologists to discern targeted from nonspecific bindings, and most AFM users do not use the AFM to measure force [35][36][37]. However, by pressing the tip of a cantilever against the sidewall of a calibrated MEMS device (Figure 4), the MEMS device can apply a known force (5) to deflect the cantilever by measured amount (4). The force balance relation would be is the unknown AFM cantilever stiffness to be determined. Since AFM x x  , then the photodiode can be simultaneously calibrated. Mass sensors. Depending on application, a mass m can be adsorbed onto a calibrated vibrating diaphragm, or m can be measured within a microfluidic channel embedded within flexures and proof mass. Although nano-to micro-mechanical resonant mass sensors have been able to detect a change in frequency due to additional mass m, the accurate measurement of that added mass has been difficult due to the inability to accurately measure the system stiffness K and mass M [38][39][40]. The ability of EMM to accurately measure MEMS stiffness K and mass M can be used to measure the added mass m of a resonating device; i.e., where m is the unknown to be determined. Gravimetry. Gravimeters are often used to measure changes in the earth's gravitational field due to large underground deposits of resources, human activities, or seismology. The accuracy of relative or absolute MEMS-based gravimeters [41,42] depend on accurate calibration of stiffness and mass. By tilting an EMM-calibrated MEMS device from being perpendicular to parallel to the direction of gravitational acceleration, then absolute gravity can be measured as  is the unknown gravitational acceleration being measured.
Altimetry. Conventional MEMS altimeters work by measuring atmospheric pressure [43,44]. The measurement of height can be an additional dimension added to GPS. As a MEMS device changes in height hg, there is a small change in gravity. The change in gravitational force on a MEMS proof mass Mass sensors. Depending on application, a mass m can be adsorbed onto a calibrated vibrating diaphragm, or m can be measured within a microfluidic channel embedded within flexures and proof mass. Although nano-to micro-mechanical resonant mass sensors have been able to detect a change in frequency due to additional mass m, the accurate measurement of that added mass has been difficult due to the inability to accurately measure the system stiffness K and mass M [38][39][40]. The ability of EMM to accurately measure MEMS stiffness K and mass M can be used to measure the added mass m of a resonating device; i.e., ω 0 = K/(M + m), where m is the unknown to be determined.
Gravimetry. Gravimeters are often used to measure changes in the earth's gravitational field due to large underground deposits of resources, human activities, or seismology. The accuracy of relative or absolute MEMS-based gravimeters [41,42] depend on accurate calibration of stiffness and mass. By tilting an EMM-calibrated MEMS device from being perpendicular to parallel to the direction of gravitational acceleration, then absolute gravity can be measured as K x = M ..
x g is the unknown gravitational acceleration being measured.
Altimetry. Conventional MEMS altimeters work by measuring atmospheric pressure [43,44]. The measurement of height can be an additional dimension added to GPS. As a MEMS device changes in height h g , there is a small change in gravity. The change in gravitational force on a MEMS proof mass of, say, M= 10 × 10 −9 kg, by the height of a person, the Empire State Building, or the cruising altitude of a passenger jet is 4 fN, 1 pN, or 26 pN, as determined by Newton's law of gravitation: where G = 6.67 × 10 −11 m 3 · kg −1 · s −2 , M Earth = 5.9 × 10 24 kg, R Earth = 6.37 × 10 6 m. Or, if gravity is assumed to be constant, the change in height ∆h g can be determined by equating gravitational potential energy to the mechanical potential energy stored in the flexure: M ..
x g ∆h g = 1 2 K x 2 − x 2 re f erence . Inertial measurement units. More accurate translational and or rotational reckoning may be required by mobile IoT devices such as drones, self-driving vehicles, robots, precision handheld tools, automated surgical tools, manufacturing machines, refined GPS (<<1 m), air/marine/space craft, etc. Present IMUs have difficulty with drift due to thermal expansion from changes in temperature [45]. EMM can be used to recalibrate IMUs as temperatures change, and because the equation of motion is known, the forces due to movement within a non-inertial reference frame can be accurately measured. These forces include: (1) the Coriolis force due to the proof mass M moving with velocity . r = .
x in a frame that is rotating with frequency vector ω: F Coriolis = −2M ω ×

R.
Since a self-calibrated IMU would be able to directly measure these non-inertial forces, then the respective unknowns (ω, . ω, r, ..

R) become measurable.
Pressure sensors/Sound sensors. Conventional MEMS pressure sensors or microphones are usually based on piezo or capacitive transduction [46][47][48]. Pressure sensors typically consist of a diaphragm that is exposed to external pressures and substances. The diaphragm of implantable pressure sensors for blood pressure monitoring can become coated with biomatter, which changes the diaphragm's effective stiffness. Using a vertical comb drive to sense diaphragm deflection, the diagram may be re-calibrated to its new stiffness prior to measurement. Additionally, if an implantable sensor consisted of a pair of calibrated pressure sensors at a known distance apart, then in addition to absolute pressure, the viscosity and velocity of blood flow may also be measured.
Diagnostic structure/Quality control. Examples of prior efforts in MEMS diagnostics include using large arrays of test structures [49] and coupled measurement plus computer modeling [50]. Once a device's mass, damping, and stiffness are measured, continuing on with the sequence of property extractions, geometric and material properties can be determined as well. In this way, an EMM diagnostic chip can accompany each process run for quality control purposes. Since EMM can be electrically probed, the testing and corresponding modification of the fabrication recipe can be automated. The diagnostic chip may also be placed at different locations about the wafer to examine how, say, Young's modulus, varies across the wafer, from wafer to wafer, from run to run. For MEMS devices that are not amenable to direct EMM calibration, an EMM test structure device can be manufactured alongside the main device, within close-proximity, such that the main IoT sensor and the EMM test structure would be subjected to closely-matched geometric and material property variations.
Experimentally-accurate modeling. The performance of a fabricated MEMS device diverges significantly from its originating CAD model due process variations and other nonidealities. By substituting a device's measured parameters into the CAD model, higher-order effects can be studied for a better understanding of the device behavior. Improved understandings can enable the causes of subtle disturbances to be more easily identified or can pinpoint areas of improvements for the next design generation to shorten design cycles. This can ultimately extend to the refined modeling of distributed IoT phenomena such as environmental models.

Performance Control
The benefits of performance control are expected to: (1) close the large gap in going from simulation to experiment, whereby the device corrects for variations in processing, packaging, and environment; (2) optimize sensing behavior to environmental conditions; (3) extend utility and dynamic range beyond the limits of structural micromachining constraints; and (4) maintain or improve manufacturing yield while improving performance.
Resonance control. MEMS sensors often utilize resonance for stronger signal to noise ratios, filtering, signal processing, timing, inertial navigation, communications, etc. As previously mentioned, prescribing a particular resonance is problematic due to process variations and requires tuning [51]. As illustrated in Figure 5a, a wafer of identically laid out resonators that were predicted to have a particular resonance frequency within CAD will differ in resonance upon fabrication. However, by adjusting the effective mass or stiffness by electrical force feedback, the desired resonance can be obtained. Moreover, the resonances of all devices can change to a different value on demand (see Figure 5b). For example, the resonance can be reduced by a factor of two by reducing the effective stiffness by a factor of four. That is, by setting the feedback electrical stiffness to be K e = −

Performance Control
The benefits of performance control are expected to: (1) close the large gap in going from simulation to experiment, whereby the device corrects for variations in processing, packaging, and environment; (2) optimize sensing behavior to environmental conditions; (3) extend utility and dynamic range beyond the limits of structural micromachining constraints; and (4) maintain or improve manufacturing yield while improving performance.
Resonance control. MEMS sensors often utilize resonance for stronger signal to noise ratios, filtering, signal processing, timing, inertial navigation, communications, etc. As previously mentioned, prescribing a particular resonance is problematic due to process variations and requires tuning [51]. As illustrated in Figure 5a, a wafer of identically laid out resonators that were predicted to have a particular resonance frequency within CAD will differ in resonance upon fabrication. However, by adjusting the effective mass or stiffness by electrical force feedback, the desired resonance can be obtained. Moreover, the resonances of all devices can change to a different value on demand (see Figure 5b). For example, the resonance can be reduced by a factor of two by reducing the effective stiffness by a factor of four. That is, by setting the feedback electrical stiffness to be , then the new resonance becomes   (a) (b) Frequency locking. Drift is a significant challenge for devices that rely on a particular resonance frequency, such as IMUs, clocks, and notch filters. Small shifts in resonance frequency from drift can result in significant amplitude attenuation due to extremely narrow bandwidths. Besides tracking frequency [52,53], it may be possible to counter drift by correcting changes in the effective stiffness in real time to maintain a constant resonant frequency. For instance, since MEMS resonance (i.e., stiffness) drifts with temperature,   K T , the MEMS device can be regarded as the slave-oscillator of a phase locked loop, where the constant driving frequency is the master-oscillator. A phase detector measures the phase between the drifting MEMS frequency and the constant driving frequency. The phase signal is then amplified, filtered, and can be used to increase or decrease the electrical feedback stiffness , PLL e K which counters the drifted frequency of the MEMS device, maintaining a desired Frequency locking. Drift is a significant challenge for devices that rely on a particular resonance frequency, such as IMUs, clocks, and notch filters. Small shifts in resonance frequency from drift can result in significant amplitude attenuation due to extremely narrow bandwidths. Besides tracking frequency [52,53], it may be possible to counter drift by correcting changes in the effective stiffness in real time to maintain a constant resonant frequency. For instance, since MEMS resonance (i.e., stiffness) drifts with temperature, K(T), the MEMS device can be regarded as the slave-oscillator of a phase locked loop, where the constant driving frequency is the master-oscillator. A phase detector measures the phase between the drifting MEMS frequency and the constant driving frequency. The phase signal is then amplified, filtered, and can be used to increase or decrease the electrical feedback stiffness K PLL e , which counters the drifted frequency of the MEMS device, maintaining a desired frequency ω 0 = [K(T) ± K PLL e ]/M = const., which effectively locks the device to a constant frequency that is independent of temperature.
Mode matching. The ability to match or follow frequencies can be beneficial. For instance, the small Coriolis force that is proportional to the velocity of the primary mode in vibratory gyros is what drives the secondary mode to amplitudes that are often on the order of nanometers. The size of the amplitude depends on how well the primary and secondary resonance modes are matched. Matching can also eliminate blind spots. Blinds spots in vibratory gyros happen when the velocity in the primary mode nears the turning points of oscillation, when velocity becomes too small to produce a significant Coriolis force for the secondary mode, F Coriolis = −2M ω × . r 1 . That is, F Coriolis → 0 as . r 1 → 0 at the turning points of oscillation, regardless of the size of the input disturbance ω during that moment. Vibratory gyros experience such blindness twice per period, which contributes to inaccurate results. However, a pair of matched gyros that maintain a 90 • phase between them should able to eliminate blind spots, whereby one gyro will have a peak velocity whenever the other has a zero velocity, such that either velocity is significant at all times.
Nonlinear dynamics. Geometric nonlinearity in stiffness affects the linearity of sensors. Such nonlinearity can be reduced or increased by a feedback force of the Duffing-type F κ = ±κ e (x τ ) 3 . Such control can be used to increase the linearity of sensors or to introduce nonlinearity to linear (small-deflection) devices [20]. The early onset of nonlinearity can be used to maintain a large resonance amplitude during frequency shifts by spring hardening or softening. For example, a feedback voltage of V = 4.5 V at an amplitude of x = 2.2 × 10 −7 m is equivalent to κ e = 4 × 10 12 N/m 3 , which is two orders greater than the purely mechanical geometric stiffness (see Figure 6). If the driving frequency is detuned to ω = 51.5 krad/s, instead of driving at ω = ω 0 = 50.0 krad/s, then as the response curve translates to the left due to an increase in temperature of, say, 100 o C, then the displacement amplitude would only change by about 0.02%. However, if the response curve was that of a linear high-Q device (vertical pole with an extremely narrow bandwidth), then a slight shift in resonance would experience significant attenuation in displacement amplitude. , which effectively locks the device to a constant frequency that is independent of temperature. Mode matching. The ability to match or follow frequencies can be beneficial. For instance, the small Coriolis force that is proportional to the velocity of the primary mode in vibratory gyros is what drives the secondary mode to amplitudes that are often on the order of nanometers. The size of the amplitude depends on how well the primary and secondary resonance modes are matched. Matching can also eliminate blind spots. Blinds spots in vibratory gyros happen when the velocity in the primary mode nears the turning points of oscillation, when velocity becomes too small to produce a significant Coriolis force for the secondary mode,  is equivalent to 12 3 4 10 N m , e    which is two orders greater than the purely mechanical geometric stiffness (see Figure 6). If the driving frequency is detuned to 51.5   krad/s, instead of driving at  x + D .
x + Kx + (κ + κ e )x 3 = F 0 cos(ωt), where M = 8 × 10 −10 kg, D = 1.55 × 10 −7 Ns/m, K = 2 N/m, κ = 4 × 10 10 N/m 3 , κ e = 4 × 10 12 N/m 3 , and F 0 = 1.7 × 10 −9 N due to a driving voltage of V 0 = 0.9 V. Note that κ e is two orders larger than the purely mechanical κ, which enables the device to experience high nonlinearity at a much smaller amplitude, which would otherwise result in a linear response. The benefit of an early onset of nonlinearity enables large amplitudes to be maintained during resonance drift; i.e., left (or right) translation of the response curve due to an increase (or decrease) in temperature.
Displacement noise reduction. An increase in effective damping can be useful for static measurements, reduction of transients, or nanoscale manipulation or positioning beyond the thermal noise limit [17,18]. For static measurements, systems that are effectively critically-damped offer the fastest route to static equilibrium by setting D e f f to be 4M e f f K e f f . This is done by feeding back a force that is proportional and opposite to sensed velocity. A reduction in displacement noise by damping could also benefit applications that would otherwise require a significant amount of averaging.
High Q. Quality factors Q are critically important for gyro sensitivity (bias • /h). Applications for gyro sensitivity range from low-end tactical (1 • /h < bias < 15 • /h, 1000 < Q < 10,000), tactical (0.1 • /h < bias < 1 • /h, 10,000 < Q < 40,000), to navigational (0.001 • /h < bias < 0.1 • /h, 40,000 < Q < 10M). As Q increases from 1000 to 10 M, cost increases from $10 to $100k. A decrease in effective damping D eff should be useful for achieving a higher effective quality factor Q eff , whereby the energy lost per cycle is fed back into the system. Achieving high-Q eff through feedback should be less costly. Usually, the higher the Q, the more difficult it is to match driving frequency with device resonance frequency due to drifting temperatures due to a narrowing vertical bandwidth. However, as modeled in Figure 6, nonlinear feedback may be viewed as effectively increasing the bandwidth, not by increasing the width of the curve, but by controllably bending the curve to the right by spring hardening or bending the curve to the left by spring softening. Such nonlinearity effectively increases the frequency range at which the driving frequency results in a large displacement amplitude. This enables high-Q resonators to more easily maintain a large displacement amplitude while the resonance drifts by using a detuned driving frequency. Previously, spring softening has been achieved by using small gap-closing electrodes with limited displacement amplitude, and spring hardening has been achieved by exploiting mechanical geometric stiffness due to large displacement amplitudes [54].
Modularity. Devices that can modify their effective mass, damping, and or stiffness should be amenable to achieving modularity, where a new device should be able to mimic the performance of the old device that it is replacing. Or conversely, a distribution of old sensors may be required to update their performances to that of new sensing specifications.
Beyond micromachining limits. It should be possible to significantly reduce the effective mass or stiffness along a particular degree of freedom to that of, say, a nanoscale device. That is, while there are micromachining limits to the flexures and proof masses for MEMS, the effective stiffness or mass can be reduced in the x-direction to that of a much smaller device. Since this method of performance enhancement does not require approaching the fragile limits of micromachining, more robust devices can be designed that maintain high manufacturing yields.

Conclusions
The small size and low power requirements of MEMS sensors are favorable attributes for IoT. IoT applications that require accurate sensing analysis with low uncertainty require the sensed input data to be also be accurate with low uncertainty. Presently, most MEMS devices are inaccurate with large uncertainties. However, we have shown that is it possible to develop MEMS devices that can autonomously calibrate themselves to achieve much higher accuracy and much lower uncertainty. And since uncomprehensive in-factory calibration is time-consuming and bottlenecks throughput, it can increase the cost MEMS devices by as much as 300%. However, self-calibration may be able to reduce this cost by reducing or eliminating the need for in-factory calibration; and since self-calibration is repeatable and reliable, the method may lead to international standardization. Another beneficial attribute of MEMS devices for IoT is their exceptional performance, where they can, for instance, achieve much higher quality factors than purely electrical systems. But performance depends on structural design, which has been constrained by manufacturing limits. However, our results suggest that such constraints may be bypassed by feedback performance control, where performances that were previously intractable could be easily achieved by controlling the effective mass, damping, and stiffness of the system. While these initial studies of self-calibration and performance control utilize the comb drive sensor/actuator, further investigations are needed to extend these methods to the other types of microtransducers. The ability of future microsensors to change their performance behaviors on demand, or to accurately self-characterize themselves, are expected to greatly extend the applications and utility of MEMS for IoT, as well as other areas.