Next Article in Journal
Kinematic Measurements of Novel Chaotic Micromixers to Enhance Mixing Performances at Low Reynolds Numbers: Comparative Study
Next Article in Special Issue
Directional and Fast Photoluminescence from CsPbI3 Nanocrystals Coupled to Dielectric Circular Bragg Gratings
Previous Article in Journal
Development of a Wearable Glove System with Multiple Sensors for Hand Kinematics Assessment
Previous Article in Special Issue
In-Fiber BaTiO3 Microsphere Resonator for High-Sensitivity Temperature Measurement
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Battery-Less Wireless Respiratory Sensor Using Micro-Machined Thin-Film Piezoelectric Resonators

Department of Electrical and Computer Engineering, University of Central Florida, Orlando, FL 32828, USA
*
Author to whom correspondence should be addressed.
Micromachines 2021, 12(4), 363; https://doi.org/10.3390/mi12040363
Submission received: 3 March 2021 / Revised: 24 March 2021 / Accepted: 25 March 2021 / Published: 27 March 2021
(This article belongs to the Special Issue Micro-Nano Science and Engineering)

Abstract

:
In this work, we present a battery-less wireless Micro-Electro-Mechanical (MEMS)-based respiration sensor capable of measuring the respiration profile of a human subject from up to 2 m distance from the transceiver unit for a mean excitation power of 80 µW and a measured SNR of 124.8 dB at 0.5 m measurement distance. The sensor with a footprint of ~10 cm2 is designed to be inexpensive, maximize user mobility, and cater to applications where disposability is desirable to minimize the sanitation burden. The sensing system is composed of a custom UHF RFID antenna, a low-loss piezoelectric MEMS resonator with two modes within the frequency range of interest, and a base transceiver unit. The difference in temperature and moisture content of inhaled and exhaled air modulates the resonance frequency of the MEMS resonator which in turn is used to monitor respiration. To detect changes in the resonance frequency of the MEMS devices, the sensor is excited by a pulsed sinusoidal signal received through an external antenna directly coupled to the device. The signal reflected from the device through the antenna is then analyzed via Fast Fourier Transform (FFT) to extract and monitor the resonance frequency of the resonator. By tracking the resonance frequency over time, the respiration profile of a patient is tracked. A compensation method for the removal of motion-induced artifacts and drift is proposed and implemented using the difference in the resonance frequency of two resonance modes of the same resonator.

1. Introduction

It is estimated that from the hundreds of millions burdened with chronic respiratory conditions globally, annually four million die prematurely [1]. This immense humanitarian and economic toll can be alleviated with proper care and monitoring of respiratory function in at-risk patients [1,2].
Respiration rate is an excellent predictor for the need to admit patients at risk of cardiac arrest into intensive care units [3]. Studies suggest that the relative changes in respiration rate are a more significant indicator of at-risk patients than systolic blood pressure in the case of physiologically unstable patients [4]. Regardless of the accuracy of diagnostic tools, limited information gathering at in-clinic visits provides a narrow snapshot of the patient’s physiology that is not measured under natural conditions [5]. This creates an inherent limitation in the quality of the clinical evaluation based on the results of diagnostic tools.
In addition to diagnostic applications, studies have shown the benefit of long-term continuous respiration monitoring for patient’s recovery [6]. As such, research results emphasize the need for respiration monitoring systems that are capable of continuously measuring both respiration rate and respiration depth [7].
Despite the stated benefits of monitoring the respiratory rate, research shows that, even in cases where the patient’s main complication is a respiratory condition, proper measurement is routinely neglected [8]. This could be the result of the relative complexity of conventional respiratory monitoring diagnostic tools that necessitate the patient be connected to a stationary measurement unit with cables and tubes. This cumbersome setup renders continuous long-term monitoring impractical for most patients.
Monitoring systems designed to fill this gap in healthcare can be categorized according to the physiological signatures that are monitored in proxy of respiration [9]. The main categories are airflow temperature [10], moisture content [11], flow rate [12], acoustic signature [13], and chest movement [2]. While these sensors lend themselves to non-intrusive continuous long-term monitoring, typically their on-sensor circuitry and batteries increase sensor unit cost, footprint, and patient discomfort and hinder continuous long-term monitoring.
To address the aforementioned needs in our previous work, a resonant wireless respiration monitoring sensor capable of measuring respiration rate without the need for onboard batteries or circuits was demonstrated for the first time [14]. With a weight of only 9 g, the sensor is composed of a Micro-Electro-Mechanical (MEMS) resonator connected to a commercial RFID antenna operating in the ISM band (@902 MHz) and was wirelessly energized from a transmitter.
Although this sensor offered a compact form-factor (~9 cm2), negligible thickness, and low unit cost, the relatively limited measurement range (~25 cm) was less than desired. Here, we build on this work by significantly improving the measurement range from a maximum of 25 to 2 m in this paper, enhancing resilience to motion-induced artifacts, and improving SNR and thus sensor accuracy without sacrificing the sensor’s small footprint and unit cost. To achieve this, we utilize a dual-mode low-loss, high Q, MEMS resonator with a custom-designed RFID antenna. We illustrated that our sensor is capable of wirelessly measuring human respiratory profiles by proxy of moisture condensation and heat exchange.
The structure of this work is as follows. After a brief overview of the operation of the sensing system, we introduce individual components of the sensor, the impact of exposure to respiratory flow on the MEMS resonator as it relates to sensing, and the details of wireless power delivery to the sensor. The results of respiratory monitoring using the proposed sensing system are presented in Section 3. Finally, the main findings of this work are briefly described in the conclusion.

2. Systems, Components and Methods

The proposed sensor, shown in Figure 1, essentially consists of a piezoelectric-based micro-machined resonator that is excited wirelessly by a pulsed sinusoidal signal transmitted from a transmitter antenna. During the off period of the excitation signal, the resonant sensor returns a decaying sinusoidal signal at the natural resonance frequency of the resonator [15] assuming that the excitation signal frequency was close to the natural resonance frequency of the MEMS resonator. The sensor’s response waveform is recovered by a receiver antenna and transferred to a high-frequency oscilloscope. Similar to the technique used in [15], time gating is utilized to separate the response of the sensor from the excitation signal at the receiver and the resonance frequencies of the sensor are extracted from the response using Fast Fourier Transform (FFT) analysis.
Maximum sensor stimulation occurs when the frequency of the excitation signal matches that of the natural frequencies of the MEMS resonators. As such, using a single signal generator we successively target the resonance frequency of the MEMS resonators. Frequency modulation of the resonance modes stemming from exposure to respiratory flow necessitates that the excitation signal frequency closely follows the resonance frequencies [16]. For this reason, a dynamic frequency locking mechanism is incorporated into the signal acquisition/generation software.
In the following section, we detail how exposure to respiratory flow impacts the MEMS device and we describe the criteria used to design and optimize these devices for passive wireless sensing.

2.1. Low-Loss, High-Quality-Factor Piezoelectric MEMS Resonator

The first rendition of the sensor in this work is made up of two RF MEMS resonators wire bonded to a custom-designed flexible dipole antenna. As will be detailed in the following sections, the goal of this multi-resonator operation is the elimination of the environmental effect (such as temperature) on the sensor.
Here, we utilize thin-film piezoelectric-on-substrate (TPoS) MEMS resonators operating near the ISM band (902–928 MHz). TPoS MEMS resonators were chosen based on their low insertion loss and high quality-factor [17]. The MEMS used here, schematically shown in Figure 1, is composed of a piezoelectric layer sandwiched between two metal layers stacked on top of a nanocrystalline diamond substrate. AlN with a thickness of 500nm forms the piezoelectric layer and the nanocrystalline substrate has a thickness of 3 µm. Details of the operation of the MEMS device, including the design and fabrication process, can be found in [18,19].
To model the sensor, an equivalent circuit for the MEMS resonator is developed. Here, the well-established modified Butterworth-Van Dyke (mBVD) model is used [20], which models the effective mass, stiffness, and loss of a single resonance mode as a series inductor ( L m ) , capacitor ( C m ) , and resistor ( R m ) , respectively. In addition, the energy stored in the electrical field formed between the ground and signal electrodes is represented as a capacitor ( C 0 ) parallel with the series R L C mechanical branch. Furthermore, a resistor modeling ohmic losses of the resonator is placed in series with the other elements. The antenna was modeled as a complex load in series with a voltage source in transmission and receiving mode, respectively.
The complete model of the sensor including the antenna and the calculated mBVD model values for the two MEMS resonators used in this paper is shown in Figure 2. Notice that MEMS#1 has two series R L C branches with each corresponding to a different resonance mode.
The impedance responses of the two resonators intended to be used for sensing are shown in Figure 3. The modes of interest were measured to be at 874.7 MHz, 879.9 MHz, and 888 MHz with an unloaded quality factor of ~1212, ~2750, and ~1750, respectively. The proximity of the resonance frequencies of the two MEMS devices assists in the design of a dipole antenna with acceptable efficiency at both frequencies. However, if the resonance frequencies are too close, we risk misidentification of the modes during measurement.
The impedances of the two resonators, as modeled by the mBVD model (shown in Figure 2), were plotted alongside the measured impedance and are shown to be in good agreement.

2.2. Sensing Mechanism

The response of the sensor to respiratory flow can be thought of as a reaction to the modulation of the MEMS resonators’ temperature and surface vapor condensation due to the periodic expulsion of air from the lungs [9]. Through disparate coupling pathways, both surface condensation and increase in temperature result in a reversible drop in MEMS device’s resonance frequency. Critical for our application, however, the resulting frequency shift due to both mechanisms is of the same polarity (negative).
Heat transferred from hot exhalatory flow to the resonator has the principal effect of lowering the effective stiffness of the resonator which in turn decreases the resonance frequency [19]. The temperature-frequency characteristic of an MEMS resonator is dependent on the properties of the material that constitutes the resonators. For a given respiration-induced temperature change of Δ T , the frequency shift can be calculated as:
Δ f T = T C F × f 0 × Δ T
Here, T C F is the temperature coefficient of frequency and has been measured to be effectively constant (~−9.6 ppm) in the temperature range of interest for the targeted lateral extensional resonance mode of the TPoS MEMS resonators [18] used here, and f 0 is the initial resonance frequency of the same mode.
The frequency shift as a function of temperature was simulated using COMSOL for all three modes measured in Figure 2. The resonance frequency for a given temperature is simulated from eigenmode analysis. To achieve this, first, the Temperature Coefficients of Elasticity (TCE) of the AlN and Mo layers were extracted from [21,22], respectively. Next, using these values we extract the TCE of the resonators’ nanocrystalline diamond layer by fitting the simulated temperature-frequency profile to the measured profile. Knowing this value, we can simulate the temperature-frequency profile of all 3 modes, as shown in Figure 4.
As exhaled air with high moisture concentration passes over the MEMS resonator, a thin layer of water vapor condenses on the surface of the device. This thin layer increases the effective mass of the resonator proportional to its coverage and thickness. Modeling an MEMS resonator as a mass-spring-damper system, the resonance frequency can be calculated as:
f 0 = 1 2 π k m 0
Here, m 0 is the initial effective mass and k is effective stiffness. Assuming a small mass loading term ( Δ m ), the frequency shift is calculated to be proportional to the effective mass of the adsorbed water vapor:
Δ f m = f 0 2 m 0 Δ m
In addition to mass loading, adsorption of moisture results in higher acoustic energy loss stemming from both acoustic mismatch between condensation and the micromechanical structure and also high acoustic loss in the condensation layer [23,24].
Both changes in temperature and vapor mass loading are periodic time-varying parameters that are coupled through heat transfer occurring during vapor condensation and evaporation. Therefore, the total sensor frequency shift corresponding to a single resonance mode resulting from exposure to respiratory flow can be written as:
Δ f ( t ) = f 0 2 m 0 Δ m ( t ) T C F × f 0 × Δ T ( t )
The same sensing pathways that facilitate respiration sensing can also act as a source of measurement error considering environmental variability. As described in Equations (3) and (4) (and illustrated in Figure 4), any change in temperature, including environmental temperature, results in a change in the resonance frequency of the resonators. If this change is sudden (such as in the case of exposure to sunlight and opening of refrigerator door), it can be incorrectly categorized as a respiration event. As shown in Figure 4, the temperature-frequency profile of the 879.9 and 888 MHz modes are almost identical in the 20–50 °C range. In effect, by monitoring the difference in resonance frequencies we can distinguish frequency change stemming from environmental temperature variation from respiration.
Δ f 888 ( t ) Δ f 879.9 ( t ) = [ f 0 2 m 0 Δ m ( t ) ] 888 [ f 0 2 m 0 Δ m ( t ) ] 879.9
The RFID antenna exhibits slightly different impedance characteristics as a function of environmental electromagnetic loading. As a result, the sensor is susceptible to motion-induced artifacts stemming from spatial and temporal changes in surrounding objects. The change in sensor frequency resulting from motional artifacts is expected to be similar for different modes; therefore, by monitoring the difference in the resonance frequencies of two different modes, we can discern respiration from such artifacts. The viability of the two compensation methods is studied experimentally in Section 3 of this work.

2.3. Antenna Design

A wearable wireless sensor that is designed to enable patient mobility requires a compact, lightweight, and omnidirectional antenna. While many commercial RFID antennas meet the above criteria, they are specifically designed for certain RFID tag chips and therefore cannot be used for our MEMS sensors. It should be noted that the impedance mismatch between the MEMS sensors and antenna is a major source of efficiency loss and thus must be kept at a minimum.
While it is feasible to incorporate an on-sensor matching network, this additional component would both increase sensor loss and complicate the development of the sensor. To this end, we developed a custom dipole antenna with lateral dimensions of ~2 × 6.5 cm2, negligible thickness, and a weight of only ~10 g. Based on [25], this simple uniplanar antenna was designed specifically for RFID tag applications and has the advantage of relatively straightforward impedance tunability whilst retaining a compact design.
To design the antenna, first, the S-parameters of the MEMS resonators were characterized using Rhode & Schwarz ZNB-8 Vector Network Analyzer (VNA). Given that the two MEMS devices are electrically connected in parallel, Keysight Advanced Design System (ADS) was used to calculate the equivalent impedance.
The major challenge in designing this antenna is achieving acceptable impedance matching between the antenna and MEMS sensor despite the difference in impedance at their respective resonance frequencies. Rather than matching the antenna perfectly at only one of the MEMS resonance modes, our solution was to match the regression line of the equivalent MEMS sensor’s impedance at the resonance modes. Whilst the power lost to impedance mismatch is not eliminated at either of the MEMS sensor’s resonance modes, this design achieves the overall minimum impedance mismatch. It should be noted that this antenna was originally designed with only 879.9 and 880 MHz modes in mind.
To realize this design, first, the antenna dimensions are estimated such that the effective length of the antenna is equivalent to a half wavelength at the antenna’s center frequency. In the next step, ANSYS High-Frequency Structure Simulator (HFSS) is utilized to fine-tune the dimensions to achieve matching to the equivalent impedance of the two MEMS. This fine-tuning process is illustrated in the HFSS simulation results shown in Figure 5. We observed that the real and imaginary parts of the antenna impedance are mainly a function of L 1 and L 3 , respectively.
The impedance profile of the finalized design with dimensions of W = 3 mm, L 1 = 35 mm, L 2   = 20 mm, and L 3   = 37.2 mm is shown in Figure 6. Using this design, the power transfer efficiency between the antenna and 3 MEMS resonance modes is shown in Table 1. The efficiency was found to be noticeably higher for the 879.9 and 888 MHz modes at 77.6% and 81.7%, respectively, as compared to the 874.7 MHz mode at 20%. This difference can be mainly attributed to the significant difference in impedance at 874.7 MHz compared to the other two modes.
The antenna was fabricated using the print-screen method [25] on a flexible Pyralux copper laminate sheet with a copper and polyimide thickness of 35 and 25 µm, respectively. In the first step of fabrication, the paper mask is printed and laminated onto the Pyralux sheet by applying heat. Next, the sheet is placed in a copper etchant bath for ~4 min.

2.4. Power Delivery to Sensor

Power delivery is a fundamental limitation of passive sensing. In our respiratory sensor, the distance between the patient and transceiver unit, and hence the quality of the mobility offered to the patient, is directly limited by the strength of the signal delivered to and received from the sensor. As such, understanding the relationship between sensing distance and wireless signal strength is crucial. For this reason, we have developed a framework that connects the transmitted power to antenna characteristics, sensor-transceiver distance, sensor characteristics, signal acquisition settings, and ultimately the strength of the sensor’s signal, as measured by the receiving unit.
For a given excitation signal, the power delivered to the passive wireless sensor is determined by the antenna gain, sensor distance, and sensor loss and can be calculated using the Friis Equation [26].
P S P T x = G S G T x λ 2 L S ( 4 π ) 2 D T S 2
Here, P S and P T x are the power of the signal delivered to the sensor and transmitted by the transmitter antenna, respectively, and G T x and G S are the gain of the transmitter and sensor antenna, respectively. D T S represents the transmitter-sensor distance. λ is the EM wavelength at the excitation signal frequency (assuming the frequency of excitation signal matches either MEMS resonance mode) and L s represents the total power lost in the sensor. Sources of loss include impedance mismatch between the resonator and antenna, Ohmic loss, dielectric loss, and losses occurring in the MEMS resonator [27]. Modifying the Friis equation for ringdown sensor measured using FFT windowing, the power of the sensor’s response as received by the oscilloscope can be calculated as:
P R x P T x = G R x G T x G S 2 λ 4 L S 2 ( 4 π ) 4 D R S 2 D T S 2 × ( 1 e T w ω / Q ) T w ω / Q
Here, T w is the FFT time window, ω is the angular frequency of the MEMS resonance mode, and Q is the effective quality factor at the receiver, as measured by the ringdown method. It should be noted that while theoretically increasing the transmitted power can significantly improve measurement range, in practice transmitted power is often limited by the RF front end (amplifier linearity and noise figure) and also concerns patient safety from prolonged exposure to high-power RF radiation. This highlights the importance of minimizing sources of energy loss in the sensor to maximize measurement distance. This equation is applied to the respiratory sensor in the following section and is utilized to calculate the actual total power loss in the sensor after all the different components are integrated.

3. Results and Discussion

The measurement setup and the acquisition flow chart used here are shown in Figure 7 and Figure 8, respectively. The transmitter and receiver antennas were connected to a Rohde & Schwarz, Munich, Germany, SMC100A signal generator and RTO 1024 oscilloscope, respectively. The signal generator is externally triggered by the oscilloscope and both are controlled and recorded through their LAN port from a remote PC unit using National Instruments, Austin, TX, USA, LabVIEW. The excitation signal alternatively targets the frequency of the resonance modes and their response in the form of a decaying sinusoidal signal is recorded by the oscilloscope and subsequently passed on to the PC for FFT analysis.
The high quality-factor of the MEMS resonance modes combined with the frequency-modulated sensing mechanism necessitates a dynamic frequency tracking system whereby the frequency of the excitation signal locks onto that of the resonance modes throughout the measurement. For this purpose, the frequency of the excitation signal targeting a given mode was set as the moving average of the preceding measured resonance frequencies of that mode.
In sensing applications, typically the size of the moving average window is determined as a compromise between responsiveness and SNR; however, in this case, given that a slow response to shifts in resonance frequency results in lower SNR, a small window is preferred. A moving average with a larger window is applied in the post-processing step. The window size for the tracking system and post-processing step was 2 and 8, respectively.
The time resolution (sampling rate) of the respiratory profile determines how faithfully the patient’s respiration is reconstructed and is a critical parameter as it sets the maximum detectable feature frequency through the Nyquist–Shannon sampling theorem. Given that the typical respiratory cycle of an adult human lasts ~4 s, our measurement system’s time resolution of ~130 ms fulfills the Nyquist frequency condition and can theoretically resolve features up to ~3.85 Hz. It should be noted that, in our measurement system, the transfer of signals acquired by the oscilloscope to the PC and their processing by LabVIEW is the determining factor for time resolution.

3.1. Resilience to Environmental Interference

The goal of incorporating more than one MEMS resonator into the sensor was to distinguish between the effects of respiratory flow and environmental variability on the sensor. The two major sources of environmental variability focused in this work are (1) discerning changes stemming from respiratory flow from environmental temperature variation and (2) negating motional artifacts stemming from small changes in the performance of the antenna as a function of spatial and temporal variability in electromagnetic environmental loading.
Given that a MEMS resonator’s resonance frequency is a function of temperature, any sudden change in environmental temperature can be miscategorized as a respiratory event. As such, by incorporating two MEMS devices with similar TCF (as shown in Figure 4), we sought to remove frequency modulation originating from environmental temperature change from the respiratory profile. In other words, we endeavored to decouple temperature-induced frequency shift from mass loading frequency shift.
The viability of this decoupling mechanism was experimentally studied by flowing hot air over the sensor and monitoring the frequency shift. This experiment is designed to replicate a sudden change in environmental temperature in the absence of any significant mass loading. As shown in Figure 9, the frequency shift exhibited by the 879.9 and 888 MHz modes is not consistent with the assumption of equal MEMS resonator temperature that underpins the decoupling mechanism. We hypothesize that this temperature gradient is in part the result of unequal heat transfer due to uncontrolled airflow at the surface of the sensor.
The RFID antenna exhibits a slightly different impedance profile depending on electromagnetic loading from its environment. The impact of the orientation of the sensor concerning the excitation and receiver antennas was tested by four separate 90° rotations of the sensor. This experiment was designed to replicate the rotation of the patient’s body. The measurement was performed for a peak and mean power of 80 mW and 80 µW at a distance of 50 cm in an anechoic chamber. The results shown in Figure 10 illustrate that sensor orientation significantly impacts the frequency corresponding to both modes, as seen by the receiver antenna.
Without a compensation mechanism, the sensing system would categorize the rotation of a patient’s head as a respiration event. However, as shown in Figure 10, by tracking the difference between the frequencies of the two modes, we can effectively identify changes in frequency associated with motion-induced artifacts. It should be noted that the same compensation mechanism was not observed between either of the two resonance modes of MEMS#1 (874.6 MHz and 879.9 MHz) and that of MEMS#2 (888 MHz). The spike-like drop in frequencies is due to the physical handling of the sensor by the experimenter’s hand.
Given the results of the two experiments, we concluded that while the sensor is capable of motion-induced artifact removal, in its current format, distinguishing rapid changes in environmental temperature from the respiratory flow is infeasible without utilizing a software-based solution. As a result, the 888 MHz resonance mode was not monitored in the following measurement results.
Using the measurement setup described above, the respiratory profile of a healthy adult was measured with a mean transmitted power of 80 µW and peak transmitted power of 80 mW. The response of the sensor, as seen by the receiver antenna, is shown in Figure 11.
As shown in Figure 12, the sensor has measured three distinct respiratory patterns, irregular, uniform-normal, and shallow, each of which is interrupted by a brief period in which the patient ceased breathing. This pattern of brief cessation of breathing events is commonly found in patients suffering from sleep apnea.

3.2. Power Delivered to the Sensor

The response of the sensor at the 874.7 and 879.9 MHz modes were measured as a function of distance from the base unit for a peak and mean transmitted power of 80 mW and 80 µW, respectively. By knowing the performance of the sensor for a given input power, we can use the modified Friis equation, Equation (7), for sensor, transmitter, and receiver antenna gains of 1.4 dBi, 8 dBi, and 8 dBi, respectively, to predict signal strength as a function of sensor distance. The simulated and measured signal power is shown in Figure 13.
By fitting Equation (5) to the measured power at the receiver, the total effective loss of the sensor for the 874.7 and 879.9 MHz modes was found to be ~15 and ~13.4 dB, respectively. The difference between these values and the simulation results in Table 1 can be attributed to discrepancies between HFSS simulation of the antenna and the fabricated antenna, wire-bond inductance, and degradation of MEMS resonator due to aging effects.
The precision of the sensing system was evaluated by calculating the measurement variance for a relatively long duration (~45 min) at a distance of 0.5 m from the receiver and excitation antenna in an anechoic chamber. Similar to the previous experiment, the peak and mean transmitted power were 80 mW and 80 µW, respectively. Following measurement, frequency drift stemming from environmental temperature variation was removed in post-processing. As shown in Figure 14, by plotting the histogram of the measurement and fitting it to a normal distribution, the standard deviation was calculated to be ~503 Hz. Using this value, we can calculate the SNR = μ/σ to be 124.8 dB.

4. Conclusions

A passive wireless respiratory sensing system based on the frequency modulation of a TPoS MEMS resonator was developed. With a size of only ~10 cm2, the sensor was shown to be capable of measuring a respiratory profile from up to a distance of 2 m from the transceiver unit for mean transmitted power of 80 µW and peak transmitted power of 80 mW. The SNR was measured to be 124.8 dB at a 50 cm sensor distance with the same transmitted power.
The respiratory sensor is composed of an antenna wire bonded to two low-loss piezoelectric MEMS resonators and is excited via a transmitter antenna and its response is recorded by a second receiver antenna. We showed that an increase in temperature and mass loading from surface vapor condensation, both stemming from exposure to exhaled airflow, results in a substantial and reversible drop in the MEMS resonance frequencies. FFT analysis is then used to extract the shift in MEMS devices’ resonance frequencies and is plotted over time to construct the patient’s respiratory profile.
To minimize energy lost to impedance mismatch between the MEMS resonators and sensor antenna, a custom RFID antenna was designed to achieve optimal impedance matching with two MEMS resonators. This was achieved despite their significant difference in impedance at their respective resonance modes.
We showed that the Friis equation, widely used to calculate the maximum operating distance for wireless passive devices, is inadequate for ringdown-type sensors. A modified version of the equation, taking into account the FFT window and sensor quality factor, was presented and used to measure the effective loss of the sensor at different operational modes.
The sensor was shown to be capable of clearly distinguishing shallow and irregular breathing from normal respiration and was also utilized to detect brief cessation of breathing events, important for studies of sleep apnea.
The sensing system utilizes a novel method for canceling the effect of motion-induced artifacts and drift on the response of the sensor. It was shown that by incorporating an MEMS resonator with two modes near the frequency of interest and monitoring the difference between their resonance frequencies, frequency shifts stemming from variation in the antenna’s position can be discerned from respiration.

Author Contributions

S.M. and R.A. developed the idea, and S.M. generated the simulation results. S.M. and P.A. collected the experimental data and wrote the paper; J.H. and X.G. developed RFID antenna, and R.A. revised the manuscript; R.A. and S.M. analyzed the results. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Science Foundation, grant number 1711632.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing not applicable.

Acknowledgments

The authors would like to thank the National Science Foundation (NSF) for the generous support of this project (grant no. 1711632) and also Rhode and Schwarz for the generous supply of equipment used in this work.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Ferkol, T.; Schraufnagel, D. The Global Burden of Respiratory Disease. Ann. Am. Thorac. Soc. 2014, 11, 404–406. [Google Scholar] [CrossRef] [PubMed]
  2. Chu, M.; Nguyen, T.; Pandey, V.; Zhou, Y.; Pham, H.N.; Bar-Yoseph, R.; Radom-Aizik, S.; Jain, R.; Cooper, D.M.; Khine, M. Respiration rate and volume measurements using wearable strain sensors. NPJ Digit. Med. 2019, 2, 1–9. [Google Scholar] [CrossRef] [PubMed]
  3. Fieselmann, J.F.; Hendryx, M.S.; Helms, C.M.; Wakefield, D.S. Respiratory rate predicts cardiopulmonary arrest for internal medicine inpatients. J. Gen. Intern. Med. 1993, 8, 354–360. [Google Scholar] [CrossRef] [PubMed]
  4. Subbe, C.P.; Davies, R.G.; Williams, E.; Rutherford, P.; Gemmell, L. Effect of introducing the Modified Early Warning score on clinical outcomes, cardio-pulmonary arrests and intensive care utilisation in acute medical admissions*. Anaesthesia 2003, 58, 797–802. [Google Scholar] [CrossRef] [PubMed]
  5. Aliverti, A. Wearable technology: Role in respiratory health and disease. Breathe 2017, 13, e27–e36. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Hmeidi, H.; Motamedi-Fakhr, S.; Chadwick, E.K.; Gilchrist, F.J.; Lenney, W.; Iles, R.; Wilson, R.C.; Alexander, J. Tidal breathing parameters measured by structured light plethysmography in children aged 2–12 years recovering from acute asthma/wheeze compared with healthy children. Physiol. Rep. 2018, 6, e13752. [Google Scholar] [CrossRef]
  7. Folke, M.; Cernerud, L.; Ekstrome, M.; Hök, B. Critical review of non-invasive respiratory monitoring in medical care. Med. Biol. Eng. Comput. 2003, 41, 377–383. [Google Scholar] [CrossRef] [PubMed]
  8. Cretikos, A.M.; Bellomo, R.; Hillman, K.; Chen, J.; Finfer, S.; Flabouris, A. Respiratory rate: The neglected vital sign. Med. J. Aust. 2008, 188, 657–659. [Google Scholar] [CrossRef]
  9. Massaroni, C.; Nicolò, A.; Presti, D.L.; Sacchetti, M.; Silvestri, S.; Schena, E. Contact-based methods for measuring respir-atory rate. Sensors 2019, 19, 908. [Google Scholar] [CrossRef] [Green Version]
  10. Marks, M.K.; South, M.; Carter, B.G. Measurement of respiratory rate and timing using a nasal thermocouple. J. Clin. Monit. 1995, 11, 159–164. [Google Scholar] [CrossRef]
  11. Kano, S.; Kim, K.; Fujii, M. Fast-Response and Flexible Nanocrystal-Based Humidity Sensor for Monitoring Human Respiration and Water Evaporation on Skin. ACS Sens. 2017, 2, 828–833. [Google Scholar] [CrossRef]
  12. Jiang, P.; Zhao, S.; Zhu, R. Smart Sensing Strip Using Monolithically Integrated Flexible Flow Sensor for Noninvasively Monitoring Respiratory Flow. Sensors 2015, 15, 31738–31750. [Google Scholar] [CrossRef]
  13. Li, S.-H.; Lin, B.-S.; Tsai, C.-H.; Yang, C.-T. Design of Wearable Breathing Sound Monitoring System for Real-Time Wheeze Detection. Sensors 2017, 17, 171. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Moradian, S.; Abdolvand, R. MEMS-based passive wireless respiration profile sensor. In Proceedings of the 2016 IEEE SENSORS, Orlando, FL, USA, 30 October–3 November 2016. [Google Scholar]
  15. Polunin, P.M.; Yang, Y.; Dykman, M.I.; Kenny, T.W.; Shaw, S.W. Characterization of MEMS Resonator Nonlinearities Using the Ringdown Response. J. Microelectromech. Syst. 2016, 25, 297–303. [Google Scholar] [CrossRef]
  16. Carr, A.R.; Patel, Y.H.; Neff, C.R.; Charkhabi, S.; Kallmyer, N.E.; Angus, H.F.; Reuel, N.F. Sweat monitoring beneath garments using passive, wireless resonant sensors interfaced with laser-ablated microfluidics. NPJ Digit. Med. 2020, 3, 1–9. [Google Scholar] [CrossRef]
  17. Abdolvand, R.; Lavasani, H.M.; Ho, G.K.; Ayazi, F. Thin-film piezoelectric-on-silicon resonators for high-frequency ref-erence oscillator applications. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2008, 55, 2596–2606. [Google Scholar] [CrossRef] [PubMed]
  18. Fatemi, H.; Zeng, H.; Carlisle, J.A.; Abdolvand, R. High-frequency thin-film AlN-on-diamond lateral–extensional reso-nators. J. microelectromech. Syst. 2013, 22, 678–686. [Google Scholar] [CrossRef]
  19. Fatemi, H.; Abdolvand, R. Low-loss lateral-extensional piezoelectric filters on ultrananocrystalline diamond. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2013, 60, 1978–1988. [Google Scholar] [CrossRef]
  20. Larson, J.; Bradley, P.; Wartenberg, S.; Ruby, R. Modified Butterworth-Van Dyke circuit for FBAR resonators and automated measurement system. In Proceedings of the 2000 IEEE Ultrasonics Symposium—An International Symposium (Cat. No.00CH37121), San Juan, PR, USA, 22–25 October 2000; Volume 1, pp. 863–868. [Google Scholar] [CrossRef]
  21. Lin, C.-M.; Yen, T.-T.; Felmetsger, V.V.; Hopcroft, M.A.; Kuypers, J.H.; Pisano, A.P. Thermally compensated aluminum nitride Lamb wave resonators for high temperature applications. Appl. Phys. Lett. 2010, 97, 083501. [Google Scholar] [CrossRef]
  22. Farraro, R.; McLellan, R.B. Temperature dependence of the Young’s modulus and shear modulus of pure nickel, platinum, and molybdenum. Met. Mater. Trans. A 1977, 8, 1563–1565. [Google Scholar] [CrossRef]
  23. Mahdavi, M.; Mehdizadeh, E.; Pourkamali, S. Piezoelectric MEMS resonant dew point meters. Sensors Actuators A Phys. 2018, 276, 52–61. [Google Scholar] [CrossRef]
  24. Nguyen, C.C.; Ngo, V.K.T.; Le, H.Q.; Li, W.L. Influences of relative humidity on the quality factors of MEMS cantilever resonators in gas rarefaction. Microsyst. Technol. 2018, 25, 2767–2782. [Google Scholar] [CrossRef]
  25. Chen, J.-P.; Hsu, P. A Compact Strip Dipole Coupled Split-Ring Resonator Antenna for RFID Tags. IEEE Trans. Antennas Propag. 2013, 61, 5372–5376. [Google Scholar] [CrossRef]
  26. Balanis, C.A. Antenna Theory: Analysis and Design; John Wiley & Sons: Hoboken, NJ, USA, 2016. [Google Scholar]
  27. Abdolvand, R.; Fatemi, H.; Moradian, S. Quality Factor and Coupling in Piezoelectric MEMS Resonators. In Stretchable Bioelectronics for Medical Devices and Systems; Springer: Berlin, Germany, 2017. [Google Scholar]
Figure 1. (a) A schematic representation of the sensing system. The resonant sensor is excited by a pulsed sinusoidal signal (blue) from a transmitter and its response (orange) as seen by a receiver is analyzed via Fast Fourier Transform (FFT). (b) Sample of a sensor used in this study. (c) The downward frequency shift of resonators from exposure to respiratory flow. (d) Respiratory profile of the healthy subject. Each fall and rise in frequency corresponds to a single respiratory cycle.
Figure 1. (a) A schematic representation of the sensing system. The resonant sensor is excited by a pulsed sinusoidal signal (blue) from a transmitter and its response (orange) as seen by a receiver is analyzed via Fast Fourier Transform (FFT). (b) Sample of a sensor used in this study. (c) The downward frequency shift of resonators from exposure to respiratory flow. (d) Respiratory profile of the healthy subject. Each fall and rise in frequency corresponds to a single respiratory cycle.
Micromachines 12 00363 g001
Figure 2. The equivalent circuit model of the sensor in receive mode. The model of the sensor (MEMS#1) is outlined in orange with the mechanical branches corresponding to each of the two modes outlined in green and gray. The model of the reference device (MEMS#2) is outlined in light blue. The antenna is modeled as a voltage source in series with a complex load is outlined in dark blue.
Figure 2. The equivalent circuit model of the sensor in receive mode. The model of the sensor (MEMS#1) is outlined in orange with the mechanical branches corresponding to each of the two modes outlined in green and gray. The model of the reference device (MEMS#2) is outlined in light blue. The antenna is modeled as a voltage source in series with a complex load is outlined in dark blue.
Micromachines 12 00363 g002
Figure 3. The measured and modeled impedance magnitude and phase of MEMS#1 (a) and MEMS#2 (b) resonators used in the sensor with the mode shapes (as simulated using COMSOL) corresponding to (a) 874.7 MHz and 879.9 MHz and (b) 888 MHz modes placed as insets.
Figure 3. The measured and modeled impedance magnitude and phase of MEMS#1 (a) and MEMS#2 (b) resonators used in the sensor with the mode shapes (as simulated using COMSOL) corresponding to (a) 874.7 MHz and 879.9 MHz and (b) 888 MHz modes placed as insets.
Micromachines 12 00363 g003
Figure 4. Simulated temperature-frequency profile of 874.7 and 879.9 MHz modes of thin-film piezoelectric-on-substrate (TPoS) MEMS resonator used in the sensor. This resonator is fabricated on a 3 µm nanocrystalline diamond-on-silicon substrate covered by 500 nm of sputtered AlN. The frequency change profile of the 879.9 and 888 MHz modes is very similar because of their common operating mode, and the small variation can be attributed to the difference in center frequency.
Figure 4. Simulated temperature-frequency profile of 874.7 and 879.9 MHz modes of thin-film piezoelectric-on-substrate (TPoS) MEMS resonator used in the sensor. This resonator is fabricated on a 3 µm nanocrystalline diamond-on-silicon substrate covered by 500 nm of sputtered AlN. The frequency change profile of the 879.9 and 888 MHz modes is very similar because of their common operating mode, and the small variation can be attributed to the difference in center frequency.
Micromachines 12 00363 g004
Figure 5. (a) Schematic of the sensor with the critical design dimensions of the RFID antenna highlighted. (b) Simulated antenna impedance for L 3   sweep and W = 3 mm, L 1   = 33 mm, and L 2   = 33 mm. (c) Simulated antenna impedance for L 1   sweep and W = 3 mm, L 2   = 33 mm, and L 3   = 30 mm.
Figure 5. (a) Schematic of the sensor with the critical design dimensions of the RFID antenna highlighted. (b) Simulated antenna impedance for L 3   sweep and W = 3 mm, L 1   = 33 mm, and L 2   = 33 mm. (c) Simulated antenna impedance for L 1   sweep and W = 3 mm, L 2   = 33 mm, and L 3   = 30 mm.
Micromachines 12 00363 g005
Figure 6. Real and imaginary impedances of the final antenna design with the fabricated antenna inserted as inset.
Figure 6. Real and imaginary impedances of the final antenna design with the fabricated antenna inserted as inset.
Micromachines 12 00363 g006
Figure 7. Measurement setup in an anechoic chamber for studying power at the receiver end as a function of sensor distance. Both the Signal generator and oscilloscope are controlled via LabVIEW from a PC placed outside the chamber.
Figure 7. Measurement setup in an anechoic chamber for studying power at the receiver end as a function of sensor distance. Both the Signal generator and oscilloscope are controlled via LabVIEW from a PC placed outside the chamber.
Micromachines 12 00363 g007
Figure 8. Signal acquisition flow chart highlighting the major components/stages of the respiratory sensing system. The Power Amplifier (PA) component is strictly optional. Here, the profile derived from the reference MEMS is not displayed in the output.
Figure 8. Signal acquisition flow chart highlighting the major components/stages of the respiratory sensing system. The Power Amplifier (PA) component is strictly optional. Here, the profile derived from the reference MEMS is not displayed in the output.
Micromachines 12 00363 g008
Figure 9. Measured temperature shift in 879.9 and 888 MHz modes in response to two separate heating periods. Before the first heating event, the temperature of the resonators (each correlating to an individual mode) was identical, and, following the second cool-down period, the temperature of the modes converged.
Figure 9. Measured temperature shift in 879.9 and 888 MHz modes in response to two separate heating periods. Before the first heating event, the temperature of the resonators (each correlating to an individual mode) was identical, and, following the second cool-down period, the temperature of the modes converged.
Micromachines 12 00363 g009
Figure 10. (ac) Measured frequency corresponding to each of the sensor’s 3 resonance modes at 50 cm distance and peak and mean powers of 80 mW and 80 µW, respectively. The sensor was rotated by 90° (highlighted in red) to simulate the sudden rotation of the patient’s body. (d) For the 874.7 and 879.9 MHz modes, this frequency shift can be canceled by subtracting the 879.9 and 874.7 MHz modes, whilst the same compensation method was found to be ineffective between 888 and 879.9 MHz (e), 888 MHz, and 874.7 MHz (f).
Figure 10. (ac) Measured frequency corresponding to each of the sensor’s 3 resonance modes at 50 cm distance and peak and mean powers of 80 mW and 80 µW, respectively. The sensor was rotated by 90° (highlighted in red) to simulate the sudden rotation of the patient’s body. (d) For the 874.7 and 879.9 MHz modes, this frequency shift can be canceled by subtracting the 879.9 and 874.7 MHz modes, whilst the same compensation method was found to be ineffective between 888 and 879.9 MHz (e), 888 MHz, and 874.7 MHz (f).
Micromachines 12 00363 g010
Figure 11. Excitation signal and sensor response, as measured by the receiving antenna, for the 874.7 (top) and 879.7 MHz (middle) modes. The FFT of the two signals are superimposed and shown in the bottom panel.
Figure 11. Excitation signal and sensor response, as measured by the receiving antenna, for the 874.7 (top) and 879.7 MHz (middle) modes. The FFT of the two signals are superimposed and shown in the bottom panel.
Micromachines 12 00363 g011
Figure 12. (a) Human respiratory measurement setup, illustrating the sensor placed on the patient’s upper lip region. (b) Respiratory profile of a human subject demonstrating non-uniform, uniform moderate, and shallow respiration in that order. Each drop in the frequency and subsequent recovery constitutes a single exhalation and inhalation cycle. Periods in which the patient ceased to breathe are highlighted in red.
Figure 12. (a) Human respiratory measurement setup, illustrating the sensor placed on the patient’s upper lip region. (b) Respiratory profile of a human subject demonstrating non-uniform, uniform moderate, and shallow respiration in that order. Each drop in the frequency and subsequent recovery constitutes a single exhalation and inhalation cycle. Periods in which the patient ceased to breathe are highlighted in red.
Micromachines 12 00363 g012
Figure 13. Signal power at receiver for a peak transmitter power of 19 dBm. The difference in measured and simulated values can be largely attributed to the impact of the directionality of the transmitter and receiver antennae at 50 cm.
Figure 13. Signal power at receiver for a peak transmitter power of 19 dBm. The difference in measured and simulated values can be largely attributed to the impact of the directionality of the transmitter and receiver antennae at 50 cm.
Micromachines 12 00363 g013
Figure 14. Top: Sensor measurement for 45 min for a peak transmitted power of 80 mW and mean transmitted power of 80 µW. Bottom: Histogram of sensor response fitted with Gaussian distribution with a standard deviation of ~503 Hz. Each count corresponds to a single measured response.
Figure 14. Top: Sensor measurement for 45 min for a peak transmitted power of 80 mW and mean transmitted power of 80 µW. Bottom: Histogram of sensor response fitted with Gaussian distribution with a standard deviation of ~503 Hz. Each count corresponds to a single measured response.
Micromachines 12 00363 g014
Table 1. Simulated power transfer efficiency between the final antenna design and 3 MEMS resonance modes. MEMS1 and MEMS2 are the sensor and reference devices, respectively. MEMS and port numbers are detailed in Figure 4.
Table 1. Simulated power transfer efficiency between the final antenna design and 3 MEMS resonance modes. MEMS1 and MEMS2 are the sensor and reference devices, respectively. MEMS and port numbers are detailed in Figure 4.
Resonance Mode (MHz)Zantenna
(Port 1)
ZMEMS1
(Port 2)
ZMEMS2
(Port 3)
S11 (dB)S21 (dB)S31 (dB)
874.715.6 + 5.8j24.109−44.57j 21.91−24.3j−3.5−3.5−9.6
879.915.9 + 9.9j24.3−19.9j9.9−60.9j−11.3−0.8−10.5
88816.4 + 16.9j9.1−62.4j21.9−24.3j−14.1−10.2−0.6
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Moradian, S.; Akhkandi, P.; Huang, J.; Gong, X.; Abdolvand, R. A Battery-Less Wireless Respiratory Sensor Using Micro-Machined Thin-Film Piezoelectric Resonators. Micromachines 2021, 12, 363. https://doi.org/10.3390/mi12040363

AMA Style

Moradian S, Akhkandi P, Huang J, Gong X, Abdolvand R. A Battery-Less Wireless Respiratory Sensor Using Micro-Machined Thin-Film Piezoelectric Resonators. Micromachines. 2021; 12(4):363. https://doi.org/10.3390/mi12040363

Chicago/Turabian Style

Moradian, Sina, Parvin Akhkandi, Junyi Huang, Xun Gong, and Reza Abdolvand. 2021. "A Battery-Less Wireless Respiratory Sensor Using Micro-Machined Thin-Film Piezoelectric Resonators" Micromachines 12, no. 4: 363. https://doi.org/10.3390/mi12040363

APA Style

Moradian, S., Akhkandi, P., Huang, J., Gong, X., & Abdolvand, R. (2021). A Battery-Less Wireless Respiratory Sensor Using Micro-Machined Thin-Film Piezoelectric Resonators. Micromachines, 12(4), 363. https://doi.org/10.3390/mi12040363

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

Article Metrics

Back to TopTop