Next Article in Journal
Experimental Study on Damage Detection in ECC-Concrete Composite Beams Using Piezoelectric Transducers
Next Article in Special Issue
Multipurpose Polymer Bragg Grating-Based Optomechanical Sensor Pad
Previous Article in Journal
A Framework for Constructing a Secure Domain of Sensor Nodes
Previous Article in Special Issue
Highly Sensitive Charge Sensor Based on Atom-Assisted High-Order Sideband Generation in a Hybrid Optomechanical System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New Design of an MOEMS Gyroscope Based on a WGM Microdisk Resonator

Key Laboratory of Micro-Inertial Instrument and Advanced Navigation Technology, Ministry of Education, School of Instrument Science and Engineering, Southeast University, Nanjing 210096, China
*
Author to whom correspondence should be addressed.
Sensors 2019, 19(12), 2798; https://doi.org/10.3390/s19122798
Submission received: 6 May 2019 / Revised: 30 May 2019 / Accepted: 3 June 2019 / Published: 21 June 2019
(This article belongs to the Special Issue Optomechanical Sensors)

Abstract

:
In this paper, we present a new design for a micro-opto-electro-mechanical (MOEMS) gyroscope based on a whispering-gallery mode (WGM) microdisk resonator and MEMS resonator. The mechanical characteristics, frequency split, and quality factor (Q) of the MEMS resonator; the optical characteristics, Q value, and coupling regimes of the WGM resonator; and the coupling between the two resonators were analyzed. Its operation principle—the transformation process from angular velocity to the resonance wavelength of the WGM resonator—is presented at same time. Next, the analysis conclusions were validated with the help of simulations in ANSYS and FDTD (Finite-Difference Time-Domain) Solutions. Afterwards, some key specifications were estimated based on the results of simulations. Lastly, the fabrication process is detailed.

1. Introduction

Gyroscopes are essential for a variety of applications ranging from inertial navigation to motion sensing in video games. The typical gyroscope can be realized using capacitive, piezoelectric, electric resistance, or optical methods. In the last few decades, micro-electro-mechanical (MEMS) gyroscopes have been widely used due to their advantages in both volume and the cost of their manufacturing. However, the drawbacks to using an MEMS gyroscope using the capacitive method, which is the main type currently in the market, include suboptimal precision and anti-electromagnetic interference. Thus, they are not suitable for use in military weapons and commercial production, which demand the use of a low-cost gyroscope with high precision. Meanwhile, gyroscopes using optical detection methods do not suffer from the same disadvantages as the MEMS gyroscope and provide superior displacement resolution [1], resilience to electromagnetic interference, and a long-range readout [2]. Furthermore, trends toward the integration of micro-photonics and microelectronics in a silicon-on-insulator (SOI) platform have made a lot of progress during recent years [3]; this makes the micro-opto-electro-mechanical (MOEMS) gyroscope an attractive option.
There are many different kinds of MOEMS inertial sensors according to their micro-photonic components, including grating, an optical cavity, and a wave guide. Among the different micro-photonic components, microdisk and microring resonators have been actively pursued for chip-scale silicon photonics. The high quality factors (Q) and small mode volumes that have been achieved in micrometer-scale devices result in a very long photon life time and strong electric field enhancement inside these resonators. The unique properties of these microdisk resonators make them ideal candidates for many applications including the filtering, sensing, delay, and enhancement of light–matter interactions [4]. Moreover, since the microdisk is a traveling-wave resonator, it is possible to realize a very high transmission efficiency of energy by appropriately designing the waveguide–resonator coupling [5]. Hence, microring and microdisk resonators on different substrates have been used previously as displacement sensors, biosensors, ultrasound sensors, and chemical sensors and in optomechanical systems [6,7,8,9,10,11,12,13]. Therefore, we propose a new type of gyroscope based on the use of whispering-gallery mode (WGM) resonators fabricated on an SOI wafer; the WGM resonator was chosen for the kernel sensing components of the gyroscope.
The gyroscope presented in this paper can be divided into two parts: the MEMS disk resonator and the microdisk optical resonator. The sense mode of the MEMS disk resonator, when operated in the drive mode under the input signal from the electrodes, will be activated when it is rotating. The WGM resonator was used to measure the deformation of the MEMS disk resonator in the sense mode, which is related to the angular velocity. By measuring the shift of the transmission spectrum of the WGM resonator, the angular velocity could be determined.

2. Structure of the MOEMS Gyroscope and Operational Principle of the MEMS Disk Resonator

2.1. Structure of the New MOEMS Gyroscope

In this paper, we designed a new type of MOEMS gyroscope based on both mechanical and optical resonators.
The entire structure of the gyroscope is shown in Figure 1a. The gyroscope was fabricated in an SOI wafer; its components are shown in Figure 1b. Figure 1 (b-1) to (b-8) shows the glass gap, the metal sealed wall, the WGM resonator, the metal pad used for electric signal transmission, the SiO2 layer, the MEMS disk resonator, the buried oxide layer, and the handle layer, respectively. Figure 1c,d shows the partial enlarged drawing, in which (c) shows the optical parts of the WGM resonator including the grating couplers, tapers, waveguide, and disk resonant cavities, and (d) shows a more detailed depiction of these optical components.
Figure 1e is a cutaway drawing of the gyroscopes; information about the parts can be found in Table 1. The operational principle and fabrication process are described in the later chapters.

2.2. Operational Principle of the MEMS Disk Resonator

The structure of the MEMS disk resonator is shown in Figure 1. From this figure, it can be seen that the resonator consists of two glass substrates with electrical connection lines and bonding pads, an electrode, a disk resonator, and 16 discrete pillar electrodes that encircle the disk and are evenly distributed. Meanwhile, there are 16 cone-shaped feed-through holes located on the top glass substrate patterned with 16 electrodes and a through hole located on the bottom glass for the DC bias connection. The operational principle of the disk resonator is shown in Figure 2; it is driven into oscillation along the drive axis, and then the Coriolis force caused by the rotation in the detecting axis, which is vertical to the plane, gives rise to the motion in its sense axis.
The angle between the drive axis and the sense axis is 45°. For the sake of a preliminary analysis of the structure of the MEMS disk resonator, the simplified equations of the motion are expressed as follows [14]:
x ¨ + 2 ξ x ω x x ˙ + ω x 2 x = F 0 sin ( ω x t ) m e f f y ¨ + 2 ξ y ω y y ˙ + ω y 2 y = 4 A g Ω x ˙
where x and y are the displacement of the drive axis and sense axis, ξ is the damping ratio, ω is the resonant frequency, meff is the effective mass, Ag is the angular gain, Ω is the angular rate, and F0 is the driving force which is defined by the input signal.
Thus, the equation of motion in the drive axis can be expressed as
x ( t ) = Q x F 0 k x sin ( ω x t π 2 ) = x 0 sin ( ω x t π 2 )
where Qx is the Q value of the drive mode, Qx = 1/2ξx, x0 is the amplitude of the driving axis, and kx is the stiffness of the drive axis. Assuming that the gyroscope is in mode matching, the mechanical sensitivity can be written as
{ | y 0 Ω | = 4 A g Q x F 0 ω x m e f f 2 ( ω y 2 ω x 2 ) 2 + ( ω y ω x Q y ) 2 = 4 A g Q x x 0 ω 0 ( k x = k y , ω x = ω y ) m e f f = V ρ ( φ x 1 2 + φ y 1 2 + φ z 1 2 ) d V = V ρ m ( φ x 2 2 + φ y 2 2 + φ z 2 2 ) d V A g = V ρ ( φ x 1 φ y 2 φ x 2 φ y 1 ) d V 2 m e f f
where (ϕx1, ϕy1, ϕz1, ϕx2, ϕy2, ϕz2) are the shape functions of the disk resonator and ρm is the density of the material.
From Equation (3), it can be concluded that the higher mechanical sensitivity can result in a larger deformation in the sense axis of the disk resonator at the same angular rate. Meanwhile, the deformation of the MEMS disk resonator in the sense axis will be transmitted to the WGM resonator and change its optic characteristics. In the next few chapters, the output value of optical detection is found to be positively correlated with the deformation of the WGM resonator. Hence, it is very important to improve the mechanical sensitivity.
Based on Equation (3), there are three possible ways to increase the mechanical sensitivity: (1) by increasing the driving force; (2) by decreasing the frequency split between the drive and sense modes through optimal design of the gyroscope; and (3) by improving the Q value in the drive and sense modes.
Based on results we have previously published [14,15], the optimized parameters of the disk resonator are shown in Table 2. The frequency split and Q value of the resonator are 13 Hz and 90,000.

3. Analysis and Simulation of the WGM Resonator

3.1. Operational Principle of the WGM Resonator

Another component of the gyroscope is the WGM resonator. The resonator comprised a straight waveguide and a disk optical cavity, which were placed close to each other so that light coupled from one to the other. When the circumference of a disk cavity is an integer number of wavelengths, the disk cavity is resonant to the wavelength and the light power stored in the disk builds up; this is referred to as the whispering-gallery mode (WGM). The light in WGM resonators can orbit thousands of times before escaping; therefore, the detection sensitivity of this method is expected to be greatly enhanced. The applied force on the resonator can then be detected by measuring the frequency shift of the WGM, which is caused by the physical elongation of the disk.
Therefore, once an applied angle velocity exerts the Coriolis force on the sense axis of the MEMS disk resonator, the WGM resonator will experience the strain, which causes the resonant wavelength to shift; the angle velocity can then be determined by measuring the frequency shift of the WGM. This indicates that the WGM resonator is the key device in the sensor, and it is essential to analyze its characteristics and optimization method.

3.2. Approximation of the Modes in a WGM Disk Resonator

Considering a WGM disk resonator with a radius R, light is confined inside the resonator by total internal reflection (TIR) because the angles between the beams and the normal directions of the boundary satisfy Snell’s law. The propagation path of light inside the resonator can be approximated as a polygon, as shown in Figure 3a, and the smaller the resonance wavelength, the higher the order of the polygon. The resonators usually satisfy the condition of R >> λ0 so that the distance of the travelling path of photons inside the optical disk resonator is equal to its equatorial circumference [16]. Assuming that the circumference of the cavity is L (L = 2πR, R is the radius of the cavity), the resonant wavelength of the disk can be expressed as
λ = n e f f L l = 2 π R n e f f l
where λ is the resonant wavelength, neff is the effective index of refraction of the optical mode in the disk cavity, and l is an integer representing the resonance order.
While the geometrical optics view gives an intuitive picture, a more quantitative description of the WGM is given by Maxwell’s Equations.
× E = i ω μ H × H = i ω ε E
In Maxwell’s equation, an isotropic medium with constant scalar permittivity and permeability is free of charge and current in order to investigate the electromagnetic fields in the disk resonators.
In the disk resonator, due to the confinement in the vertical direction, the photon is restricted in a plane. Hence, there are two polarizations: TE (E field parallel to the disk plane) and TM (E field perpendicular to the disk plane). Thus, Equation (9) becomes scalar in the z direction and Fz corresponds to Hz (Ez) for the TE (TM) modes. Assuming a separable solution for Fz(ρ,φ,z), we then have the following:
Φ Z 2 Ψ ρ 2 + Φ Z 1 ρ Ψ ρ + Ψ Z 1 ρ 2 2 Φ ϕ 2 + Ψ Φ 2 Z z 2 = β 2 Ψ Φ Z 1 Ψ 2 Ψ ρ 2 + 1 Ψ 1 ρ Ψ ρ + 1 Φ 1 ρ 2 2 Φ ϕ 2 + 1 Z 2 Z z 2 = β 2
If we let
2 Z z 2 = β z 2 Z 2 Φ ϕ 2 = l 2 Φ β ρ 2 + β z 2 = β 2
then we get a classic Bessel Differential Equation:
ρ 2 2 Ψ ρ 2 + ρ Ψ ρ + [ ( β ρ ρ ) 2 l 2 ] Ψ = 0 .
In the disk resonator, the general solution of Bessel functions of the first order is as follows:
Ψ = A J l ( β ρ ρ ) + B Y l ( β ρ ρ ) .
Jl(βρρ) and Yl(βρρ) represent the Bessel functions of the first and second kind, in which l is the angular mode number. Because Yl (βρρ) is divergent when ρ tends to 0, and the field should be finite at R = 0, the solution can be expressed as
Ψ = A J l ( β ρ ρ ) ( ρ R )
where A is a constant dependent on the boundary conditions.
Meanwhile, the transient field exponentially attenuates along the radial direction:
Ψ = J l ( β ρ R ) exp ( α ( ρ R ) ) ( ρ > R )
where α = 2 π n e f f 2 n 0 2 / λ .
The solution of the first line in Equation (7) follows the standard slab mode calculation:
Z ( z ) = cos ( β z z )
where β = mπ/t, t is the thickness of the disk resonator, and m is the longitudinal mode number.
With the help of the FDTD Solution, the first- and second-order whispering-gallery modes of the resonator were determined. Through changing the disk resonator radii from 200 nm to 600 nm in the simulation, we found that the resonator has the same resonance wavelength of 1.55 µm with radii of 418 nm and 428 nm, while the excited WGM modes are of the first and second order, respectively, as shown in Figure 4.

3.3. Quality Factor Q of the WGM Resonator

The loss of a resonator is an important parameter and is used to determine the sensitivity of the gyroscope. Resonator loss is commonly expressed in terms of the quality factor Q. The quality factor Q is related to the cavity photon lifetime τ, as represented in the following equation [16]:
Q = λ 0 Δ λ = ω π
where ω is the optical frequency (ω = 2πc/λ).
The quality factor value Q primarily depends on several loss mechanisms in the optical WGM resonators. The overall Q value can be expressed as
Q 1 = Q c o u p l e 1 + Q b e n d 1 + Q s s 1 + Q b u l k 1 + Q w 1
where Qcouple is the loss induced by the coupling of the straight waveguide and the optical resonator, which will be discussed in next section; Qbend is due to the bending loss of the waveguide and the WGM resonator; Qss is the scattering losses from the residual surface in homogeneities; Qw is the absorption losses due to contamination on the surface of the resonator; and Qbulk is the bulk absorption in the material.

3.4. The Model of the Theoretical Approach of Single Waveguides Coupled to WGM Resonators

It is impossible for a WGM resonator to perform measurement without an input–output interface. Hence, a typical resonator having a microdisk as its resonant cavity that is closely coupled with one straight bus waveguide, which serves as the input and output, was designed, as shown in Figure 5. Based on the optical field analysis performed in the previous sections, the light from the waveguide coupled to the disk was found to be an evanescent wave through the narrow gap, which is illustrated in Figure 5b.
Based on the coupled mode theory [17], the coupled mode equation is given as
[ E 2 E 3 ] = [ t j κ j κ t ] [ E 1 E 4 ]
where t is the self-coupling coefficient of the waveguide and disk, while κ is the mutual coupling coefficient for both of them. t2 + κ2 = (1 − τ)2, where τ is the loss factor.
Solving the above equations yields the expression for the transmission factor in relation to the input waveguide:
T = | E 2 E 1 | 2 = α 2 + | t | 2 2 α | t | cos ( θ + ϕ t ) 1 + α 2 | t | 2 2 α | t | cos ( θ + ϕ t )
where α and θ represent the loss and the phase shift, respectively, after one circulation. α = eτL and θ = βL where β is the propagation constant, which depends on the effective refractive index of the resonator as β = 2πneff/λ; ϕ t is the additional phase shift in the resonator and the length L = 2πR.
When in the resonance mode, the coupling can be considered as one of three regimes:
(1)
Under-coupled regime: When the waveguide is far from the resonator, the coupling is very weak. When the waveguide moves close to the resonator, the overlap of the waveguide mode and WGMs will increase. However, when |t| > α, the transmission T decreases continuously from unity and approaches zero.
(2)
Critical coupled regime: When the waveguide gets closer to the resonator such that |t| = α, the normalized transmission T is zero, which means that all of the input power is coupled to the resonator.
(3)
Over-coupled regime: When the gap between the waveguide and the resonator is further decreased, the overlap of the modes becomes larger and the normalized transmission T becomes greater than 0 again; this is called the over-coupled regime.
The classical normalized transmission spectra of the three different coupling regimes are shown in Figure 6.
For the above scenarios, we use power normalization so that |E1|2 and |E2|2 are the respective traveling wave powers. We will, without experiencing a loss of generality, take the incident power |E1|2 to be unity. At resonance, θ + ϕ t = 2 π n , n is an integer, and
| E 2 | 2 = ( α | t | ) 2 ( α α | t | ) 2 .
The relationship between the light intensity and the self-coupling coefficient t at resonant frequency when α = 0.9 and 0.99 is plotted in Figure 7.
For high-Q resonators, when α is close to unity, the portion of the curve to the right of the critical coupling point is extremely steep. Small changes in α for a given t, or vice versa, can control the transmitted power E2 between unity and zero. An understanding of the relationship between the deformation of the resonator caused by the rotation and α and/or t will allow for the identification of a method to optimize the performance of the gyroscope.

4. Analysis and Simulation of the MOEMS Gyroscope

In the last two chapters, the two main elements of the MOEMS gyroscope were introduced. The mechanical properties of the MEMS resonator, which converts angular velocity to displacement of the sense axis in the disk, were analyzed. At the same time, the optical properties of the WGM resonator were detailed.
Based on that, the optional principle of the MOEMS gyroscope can be concluded. When an applied angular velocity exerts the Coriolis force on the MOEMS structure, the force on the disk resonator causes its rings to become deformed. Meanwhile, the WGM resonator will experience the same deformation for it is manufactured on the ring. The deformation of the WGM resonator also means a change in its circumference or the optical distance L. Based on Equation (4), the optical resonator wavelength is related to L. Hence, it can be found that the applied angular velocity can be measured according to the shift in the optical resonator frequency.
In summary, the function of the MOEMS gyroscope is based on the transformation process from angular velocity to deformation of the MEMS resonator, to deformation of the WGM resonator, to a shift in the resonator frequency. In this section, the whole process will be discussed.

4.1. The Relationship between Frequency Shift and Deformation of the WGM Resonator

Due to physical elongation, the shape of the disk is modified. Assuming that the path length L of the disk is changed by ΔL, the shift in its resonant wavelength Δλ can be expressed as follows:
λ m = L m n e f f l = = ( L + Δ L ) n e f f l Δ λ = λ m λ = L m n e f f l L n e f f l = Δ L n e f f l
When the disk is strained, the diameter becomes larger in the direction of the strain and smaller in the perpendicular direction, and a circular disk becomes an ellipse. The deformation of the WGM resonator is related to the strain of the MEMS resonator for it is on the outer ring of the disk. Hence, the equation can be written as
Δ λ = λ m λ = L m n e f f l L n e f f l = { [ 2 π b + 4 ( a b ) ] 2 π r 0 } n e f f l
where r0 is the initial radius of the optical disk resonator and a, b are the semimajor axis and semiminor axis of the elliptical WGM resonator.
Because the WGM resonator is on the MEMS resonator and its thickness is much smaller (less than 1%), the deformation of the WGM resonator is exactly the same as the strain of the silicon element below the corresponding position. Based on the above formula and the definition of strain, the following results are obtained:
ε = Δ L L = [ 2 π b + 4 ( a b ) ] 2 π R 2 π R = Δ λ λ Δ λ = ε λ

4.2. The Relationship between Deformation of the WGM Resonator and Stress on the MEMS Resonator

From the last section, it is known that the resonator wavelength shift of the WGM resonator is proportional to its strain. If the relational expression between deformation of the optical resonator and stress on the disk resonator is known, the wavelength shift caused by the angular velocity can be calculated combined with Equations (3) and (20).
In order to facilitate the analysis, a simplified MOEMS gyro deformation model is presented, as shown in Figure 8. The radii of the optical resonators discussed in this paper are smaller than the width of the ring in the MEMS resonator, so the deformation of optical cavity is closely related to the deformation of the outermost ring caused by strain in the ring. With the help of the Timoshenko Formula, the distribution of strain in the ring can be obtained:
M ( θ ) = 1 2 F × R a ( cos θ 2 π ) σ θ θ ( θ ) = ± 3 2 M ( θ ) H δ 2 = ± 3 4 1 H δ 2 F × R a ( cos θ 2 π )
where F is the force in the radial direction, Ra is the center diameter of the ring, θ is the azimuthal angle, M(θ) is the bending moment of the cross section at angle θ, σθθ(θ) is the hoop force at angle θ, H is the thickness of the ring, and δ is the displacement of the ring.
Based on other researchers’ work, the complete stress distribution is [18,19]
σ θ θ ( α , θ ) = 2 F π R 0 [ 1 2 ( 1 α cos θ ) sin 2 θ ( 1 + α 2 2 α cos θ ) 2 ( 1 + α cos θ ) sin 2 θ ( 1 + α 2 + 2 α cos θ ) 2 ] σ r r ( α , θ ) = 2 F π R 0 [ 1 2 ( 1 α cos θ ) ( cos θ α ) 2 ( 1 + α 2 2 α cos θ ) 2 ( 1 + α cos θ ) ( cos θ + α ) 2 ( 1 + α 2 + 2 α cos θ ) 2 ] σ r θ ( α , θ ) = 2 F π R 0 [ ( 1 α cos θ ) ( cos θ α ) sin θ ( 1 + α 2 2 α cos θ ) 2 + ( 1 + α cos θ ) ( cos θ + α ) sin θ ( 1 + α 2 + 2 α cos θ ) 2 ]
where R0 is the outer radius of the ring; α is the ratio of the distance from a point to the center to the outer radius (α = Ra/R0); and σθθ(α,θ), σrr(α,θ), and σ(α,θ) are the hoop force, radial force, and shear force at angle θ with ratio α, respectively.
Thus, the main strain force τmax is
{ p ( α , θ ) = 1 2 [ σ θ θ ( α , θ ) + σ r r ( α , θ ) ] + 1 2 [ σ θ θ ( α , θ ) σ r r ( α , θ ) ] 2 + 4 σ r θ 2 ( α , θ ) q ( α , θ ) = 1 2 [ σ θ θ ( α , θ ) + σ r r ( α , θ ) ] 1 2 [ σ θ θ ( α , θ ) σ r r ( α , θ ) ] 2 + 4 σ r θ 2 ( α , θ ) τ max ( α , θ ) = p q 2 = 1 2 [ σ θ θ ( α , θ ) σ r r ( α , θ ) ] 2 + 4 σ r θ 2 ( α , θ )
The strain ε can be calculated by
ε = σ E = S f [ α , θ , τ max ( α , θ ) ] d δ E
where E is the Young’s modulus of the ring, σ is the combined stress on the element corresponding to the location of the WGM resonator, and f(α,θ,τ) is the surface equation of the strain force. The ring coordinates (α,θ) are transformed from the orthogonal coordinates (x,y):
{ α = x 2 + y 2 R 0 θ = arctan x y
Considering the complexity of the disk resonator, it is impossible to get the analytic solutions directly, so the relationships between deformation of the disk resonator and the optical cavity with different radii were analyzed with the help of the finite element simulation software ANSYS Workbench. The results of the static simulation are shown in Figure 9.
Using the simulation results, more credible results were obtained in support of the theoretical analysis above:
{ Δ x = K Y x Y Δ y = K Y y Y
where △x and △y are the variations of the WGM resonator radius in the tangential and radial directions of the disk resonator, respectively; Y is the maximum displacement of the disk resonator in the force direction; and KYx and KYy are coefficients concluded from the simulation results. When r0 = 5 µm or 3 µm, KYx = 2.301 × 10−4 or 1.4616 × 10−4 and KYx = 3.807 × 10−6 or 2.438 × 10−6.

5. Simulation of the MOEMS Gyroscope

In the last section, the transformation processes of deformation to a resonant wavelength shift in the WGM resonator and stress on the MEMS resonator to deformation of the WGM resonator were discussed. In this section, a simulation test was used to prove the working principle of the MOEMS gyroscope and calculate its sensitivity and cross-axis sensitivity.

5.1. Simulation Analysis of the Resonant Wavelength Shift of a Deformed WGM Resonator

Based on the conclusion of Section 4, the semimajor axis and semiminor axis of the elliptical WGM resonator are rx = R ± △x and ry = R ± △y. Because △x >> △y, the variation △y can be ignored, and the maximum of wavelength shift is
Δ λ max λ = Δ L L [ 2 π b + 4 ( a b ) ] 2 π R 2 π R ± 4 Δ x 2 π r 0 = ± 2 Δ x π r 0 = ± k x λ Δ x r 0 .
Assume that r0 = 3 µm, 5 µm, or 10 µm and that △x = (10−4 ~ 0.01)*r0; with the help of FDTD Solutions, the transmission spectrum of the WGM resonators can be calculated. With △x as the transverse axis and △λ as the longitudinal axis, the results are arranged as shown in Figure 10.
According to the curves on the right-hand side in Figure 10a–c, the resonant wavelength shift △λ is proportional to the deformation △x in WGM resonators with different radii.
Dividing the deformation △x by the initial radius r0 and comparing this with the resonance wavelength shift △λ divided by the initial resonance wavelength λ0 gives the results shown in Figure 11.
Comparing this with Equation (27), it can be seen that the relative variation of the resonance wavelength λ/λ0 is proportional to the strain x/r0, and the scale factor k is independent of the initial radius ( Δ λ λ = k x λ Δ x r 0 , k ≈ 0.4486). Hence, the results obtained by simulation and by theoretical calculation are similar.

5.2. Key Specifications of the MOEMS Gyroscope

In Section 4.2, the relationship between deformations of the MEMS resonator and the optical resonator was preliminarily concluded from static structure simulations. With the help of harmonic response simulation, the WGM resonator deformation at an applied angular velocity can be calculated. By combining all the conclusions above, some key specifications of the gyroscope were estimated.
(1) Sensitivity
Assuming that the drive force of the gyroscope is 4.6 × 10−6 N, the resulting structural parameters of the MEMS resonator are given in Table 2; the resonant frequencies of the drive and sense modes are 8324.8 Hz and 8326.2 Hz, and the results of the harmonic response simulation are shown in Figure 12. It can be seen that the maximum drive amplitude of the MEMS resonator was 2.9329 µm (X0 = 2.9329 µm) when it was at the resonance point.
Combining Equations (3), (26), and (27), we get
{ | Y Ω | = 4 A g Q x X 0 ω 0 Δ x = K Y x Y Δ λ max λ = Δ L L ± k x λ Δ x r 0
Hence, the sensitivity of the MOEMS gyroscope is
S M O E M S = | Δ λ Ω | = | Δ L Ω L λ | | [ 2 π b + 4 ( a b ) ] 2 π R 2 π R Ω λ | = k x λ K Y x r 0 4 A g Q x 0 ω 0 λ
Considering the Q value and the complexity of the transmission spectrum, a WGM resonator with a 5 µm radius was selected in the MOEMS gyroscope; consequently, λ0 = 1.59732 µm, Kyx = 2.301 × 10−4, and K = 0.4886. Substituting all these parameters into Equation (29), we get SMOEMS = 3.002 pm/(°·s−1).
(2) Measurement Range
Once an angular velocity is applied to the MOEMS gyroscope, its transmission spectra will drift due to the resonance wavelength shift. In order to avoid inaccuracies, the drift (△λrange) should be less than the separation between successive resonances; this separation is also called the free spectral range (FSR). Meanwhile, for the sake of minimizing the interference from transmission spectra overlapping, we set △λrange = 0.5*FSR. Based on the results of the FDTD simulation, as shown in Figure 13, it can be seen that FSR = 0.01004 µm and △λrange = 5.02 nm. If the watch window is set in a proper range, the drifting is clear and the shift in the resonance wavelength can be picked up easily.
Knowing △λrange, the range of measurable angle velocity can be calculated by Ω r a n g e = ± 1 2 | Δ λ r a n g e S M O E M S | = ± 836.11 ° / s .
(3) Nonlinear error
Within the measurement range of the MOEMS gyroscope, the maximum shift in the resonance wavelength is ±2.51 nm; based on Equation (27), the radius variation of the WGM resonator, △x, is about ±16.4679 nm. Hence, by comparing the linear and the nonideal results, as shown in Figure 14, the nonlinear error is calculated by following equation.
δ = d ( Δ λ ) max Δ λ r a n g e / 2 × 100 %
Thus, the nonlinear error of the MOEMS gyroscope in measurement is about 1.67% based on the simulation results.
(4) Cross-Axis Coupling and Equivalent Signals
In the gyroscope, due to fabrication imperfections disrupting the isotropy of the MEMS disk resonator, causing frequency mismatch and mode coupling [20], we consider the cross-axis coupling. Equation (1) can be re-written as
m e f f x ¨ + 2 m e f f ξ x ω x x ˙ + ω x 2 m e f f x + k y x y = F 0 sin ( ω x t ) y ¨ + 2 ξ y ω y y ˙ + ω y 2 y = k x y x + 4 A g Ω x ˙
where kxy and kyx are the equivalent cross-axis-coupling stiffness.
According to the method in Section 4.2, the simulation was used to analyze the cross-axis coupling. The position of the WGM resonator and applied force is shown in Figure 15a, while the results of the optical resonator deformation under variable force are given in Figure 15b.
Thus, we can get the relationship between the deformations of the two elements:
{ Δ x 0 Δ y = K X y , q X
where KXy,q is the coefficient concluded from the simulation results. When r0 = 5 µm, KYx = 2.28 × 10−8.
Referring to Equation (33), the cross-axis coupling can be calculated as follows.
{ Δ y = K X y , q X 0 Δ λ q = | Δ L L λ | k x λ Δ y r 0 λ
Thus, the wavelength shift caused by cross-axis coupling is about 0.01 pm. Based on the sensitivity of the gyroscope, the orthogonal signal which is equivalent to the angular velocity is about 0.01°/s when only the cross-axis coupling is considered.

6. Process and Result of Fabrication

The gyroscope can be fabricated in a special silicon-on-insulator (SOI) wafer with five layers. From the top to the bottom of the wafer is an SiO2 layer 2 µm thick, a device layer (Si) 60 µm thick, a buried oxide layer 1µm thick, a handle layer 400 µm thick, and an SiO2 layer 2 µm thick. The device layer and the handle layer are made on low-resistivity, boron-doped (P-type) <100> single crystal silicon wafer. BF33 borosilicate glass was chosen as the material for the substrates and for the cap of the gyroscope due to its close coefficient of thermal expansion to that of silicon. The main process flows are shown in Figure 16:
(1)
A 2 µm plasma-enhanced chemical vapor deposition (PECVD) silicon nitride (SiNx) layer was deposited on the top and bottom sides of the SOI wafer.
(2)
The photoresist was put on the bottom side of the wafer through spin coating and patterned by mask1 to define the holes in the handle layer.
(3)
The through-holes in the handle layer were formed through deep reactive ion etching (DRIE), and then part of the buried oxide layer was etched by the input KOH vapor from the holes.
(4)
The bottom SiO2 layer and SiNx layer were polished after Step (3), then a layer of PMMA was applied on the top SiNx layer. After that, the pattern of the optical disk resonator and the wave guide was defined using e-beam lithography.
(5)
The remaining PMMA was washed off after the WGM resonator and the wave guide were created through reactive ion etching (RIE), and then a new PMMA layer was applied on the same side. Next, the pattern of the grating coupler was defined on the PMMA using e-beam lithography, and the coupler was made with the help of inductively coupled plasma (ICP) etching.
(6)
The PMMA was cleaned off after Step (5), and then another SiO2 layer was deposited on the front of the wafer to protect the optical devices made in the last step.
(7)
The photoresist was spun on the deposited SiO2 layer, and then the pattern of metal pads was defined on the photoresist using mask2.
(8)
With the pattern from Step (7), the position of the pad in the device layer was exposed after wet etching, and then 30/300 nm thick Cr/Au was deposited in the pattern to form the pads through a lift-off process. After that, the remaining photoresist was cleaned off the wafer.
(9)
We spun the photoresist again and defined the pattern of the MEMS disk resonator in the photoresist using mask3.
(10)
The resonator structure and the electrodes were simultaneously released using the Bosch ICP.
(11)
The pattern of the bonding point and the metal envelope wall in the glass, which would be used for gold silicon bonding, was created. As in Steps (7) and (8), the pattern was defined using mask4 and the object was formed through a lift-off process.
(12)
We spun the photoresist on the glass and then defined the pattern of the cavity using mask5, then wet etched the cavity in the glass.
(13)
We spun the photoresist on another side of the glass and defined the pattern of the through holes using mask6, then wet etched the holes in the glass.
(14)
The glass cap was bonded with the SOI through gold–silicon bonding with the bonding point and the metal envelope wall. Then, another glass wafer was bonded on the bottom of the SOI wafer using the anodic bonding method.

7. Conclusions

The results show that microdisk WGM resonators will experience optical resonance shifts when exposed to the deformation force caused by an MEMS resonator under an external angular velocity. The analysis and simulation results explained the operational principle of the optical gyroscope and gave us a few different ways to improve the performance of the gyroscope through mechanical and optical elements. With these results, some key specifications including sensitivity, measurement range, nonlinear error, and bias instability were calculated and analyzed. The fabrication process was also discussed. These results are only the first step towards the design and development of angular velocity sensors that are based on this principle; our future work will focus on manufacturing, experiments, and improvement of the sensor.

Author Contributions

D.X. and L.Z. conceived and designed the study; L.H. finalized the simulation; L.H. analyzed the results; all authors wrote the paper.

Acknowledgments

This work was supported in part by the National Natural Science Foundation (No. 61571127, No.61871125), the Joint Fund of Ministry of Education for Equipment and Pre-research (6141A02022333), and the Pre-research Fund (61405170103).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Krishnamoorthy, U.; Olsson, R.H.; Bogart, G.R.; Baker, M.S.; Carr, D.W.; Swiler, T.P.; Clews, P.J. In-plane MEMS-based nano-g accelerometer with sub-wavelength optical resonant sensor. Sens. Actuators A Phys. 2008, 145–146, 283–290. [Google Scholar] [CrossRef]
  2. Krause, A.G.; Winger, M.; Blasius, T.D.; Lin, Q.; Painter, O. A high-resolution microchip optomechanical accelerometer. Nat. Photonics 2012, 6, 768–772. [Google Scholar] [CrossRef] [Green Version]
  3. Soltani, M.; Yegnanarayanan, S.; Adibi, A. Ultra-high Q planar silicon microdisk resonators for chip-scale silicon photonics. Opt. Express 2007, 15, 4694–4704. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Barclay, P.E.; Srinivasan, K.; Painter, O.; Lev, B.; Mabuchi, H. Integration of fiber-coupled high-Q SiNx microdisks with atom chips. Appl. Phys. Lett. 2006, 89, 131108. [Google Scholar] [CrossRef]
  5. Niehusmann, J.; Vörckel, A.; Bolivar, P.H.; Wahlbrink, T.; Henschel, W.; Kurz, H. Ultrahigh-quality-factor silicon-on-insulator microring resonator. Opt. Lett. 2004, 29, 2861–2863. [Google Scholar] [CrossRef] [PubMed]
  6. Li, J.J.; Zhu, K.D. All-optical mass sensing with coupled mechanical resonator systems. Phys. Rep. 2013, 525, 223–254. [Google Scholar] [CrossRef]
  7. Bhola, B.; Steier, W.H. A Novel Optical Microring Resonator Accelerometer. IEEE Sens. J. 2007, 7, 1759–1766. [Google Scholar] [CrossRef]
  8. Chao, C.Y.; Guo, L.J. Design and optimization of microring resonators in biochemical sensing applications. J. Lightwave Technol. 2006, 24, 1395–1402. [Google Scholar] [CrossRef] [Green Version]
  9. Cho, S.Y.; Jokerst, N.M. A Polymer Microdisk Photonic Sensor Integrated onto Silicon. IEEE Photonic Technol. Lett. 2006, 18, 2096–2098. [Google Scholar] [CrossRef]
  10. Morand, A.; Zhang, Y.; Martin, B.; Phan Huy, K.; Amans, D.; Benech, P.; Verbert, J.; Hadji, E.; Fédéli, J. Ultra-compact microdisk resonator filters on SOI substrate. Opt. Express 2006, 14, 12814–12821. [Google Scholar] [CrossRef] [PubMed]
  11. Kiyat, I.; Kocabas, C.; Aydinli, A. Integrated micro ring resonator displacement sensor for scanning probe microscopies. J. Micromech. Microeng. 2004, 14, 374. [Google Scholar] [CrossRef]
  12. Fang, W.; Buchholz, D.B.; Bailey, R.C.; Hupp, J.T.; Chang, R.P.H.; Cao, H. Detection of chemical species using ultraviolet microdisk lasers. Appl. Phys. Lett. 2004, 85, 3666–3668. [Google Scholar] [CrossRef]
  13. Krioukov, E.; Klunder, D.J.; Driessen, A.; Greve, J.; Otto, C. Sensor based on an integrated optical microcavity. Opt. Lett. 2002, 27, 512–514. [Google Scholar] [CrossRef] [PubMed]
  14. Xia, D.; Huang, L.; Xu, L.; Gao, H. Structural Analysis of Disk Resonance Gyroscope. Micromachines 2017, 8, 296. [Google Scholar] [CrossRef] [PubMed]
  15. Xia, D.; Hu, Y.; Kong, L.; Chang, C. Design of a Digitalized Microgyroscope System Using ΣΔ Modulation Technology. IEEE Sens. J. 2015, 15, 3793–3806. [Google Scholar] [CrossRef]
  16. Le, T.M. Integration of Practical High Sensitivity Whispering Gallery Mode Resonator Sensors. Bachelor’s Thesis, USC, Los Angeles, CA, USA, 2012. [Google Scholar]
  17. Yariv, A. Critical coupling and its control in optical waveguide-ring resonator systems. IEEE Photonic Technol. Lett. 2002, 14, 483–485. [Google Scholar] [CrossRef] [Green Version]
  18. Chianese, R.B.; Erdlac, R.J. The general solution to the distribution of stresses in a circular ring compressed by two forces acting along a diameter. Q. J. Mech. Appl. Math. 1988, 41, 239–247. [Google Scholar] [CrossRef]
  19. Muskhelishvili, N.I. Some Basic Problems of the Mathematical Theory of Elasticity; Springer: Dordrecht, The Netherlands, 1977. [Google Scholar]
  20. Painter, C.C.; Shkel, A.M. Identification of anisoelasticity for electrostatic trimming of rate-integrating gyroscopes. In Proceedings of the Smart Structures and Materials 2002: Smart Electronics, MEMS, and Nanotechnology, San Diego, CA, USA, 11 July 2002; SPIE: New York, NY, USA, 2002; pp. 157–168. [Google Scholar]
Figure 1. Diagrams of the structure of the micro-opto-electro-mechanical (MOEMS) gyroscope: (a) Diagram of the entire gyroscope; (b) Exploded view; (c,d) Partial enlarged drawing; (e) Cutaway drawing.
Figure 1. Diagrams of the structure of the micro-opto-electro-mechanical (MOEMS) gyroscope: (a) Diagram of the entire gyroscope; (b) Exploded view; (c,d) Partial enlarged drawing; (e) Cutaway drawing.
Sensors 19 02798 g001
Figure 2. Diagram of the micro-electro-mechanical (MEMS) disk resonator modes: (a) Drive mode; (b) Sense mode.
Figure 2. Diagram of the micro-electro-mechanical (MEMS) disk resonator modes: (a) Drive mode; (b) Sense mode.
Sensors 19 02798 g002
Figure 3. (a) Geometric schematic of the whispering-gallery mode (WGM) resonance; (b) Electromagnetic field distribution of the WGM resonance.
Figure 3. (a) Geometric schematic of the whispering-gallery mode (WGM) resonance; (b) Electromagnetic field distribution of the WGM resonance.
Sensors 19 02798 g003
Figure 4. (a) Magnetic field of the WGM resonator in the first order (R = 418 nm); (b) Magnetic field of the WGM resonator in the second order (R = 428 nm).
Figure 4. (a) Magnetic field of the WGM resonator in the first order (R = 418 nm); (b) Magnetic field of the WGM resonator in the second order (R = 428 nm).
Sensors 19 02798 g004
Figure 5. Model of the WGM resonator: (a) geometric optical mode; (b) optical field distribution.
Figure 5. Model of the WGM resonator: (a) geometric optical mode; (b) optical field distribution.
Sensors 19 02798 g005
Figure 6. Normalized transmission output in different coupling regimes.
Figure 6. Normalized transmission output in different coupling regimes.
Sensors 19 02798 g006
Figure 7. The universal transmission plot for the configuration in Figure 5.
Figure 7. The universal transmission plot for the configuration in Figure 5.
Sensors 19 02798 g007
Figure 8. A simplified model of disk resonator deformation.
Figure 8. A simplified model of disk resonator deformation.
Sensors 19 02798 g008
Figure 9. (a) Deformation of the MEMS disk resonator under static force; (b) Radius variation of the optical resonator; (c) Ratios of radii variation in optical and MEMS resonators.
Figure 9. (a) Deformation of the MEMS disk resonator under static force; (b) Radius variation of the optical resonator; (c) Ratios of radii variation in optical and MEMS resonators.
Sensors 19 02798 g009
Figure 10. FDTD simulation results of the relationship between △x and △λ in WGM resonators with different radii: (a) r0 = 3 µm; (b) r0 = 5 µm; (c) r0 = 10 µm.
Figure 10. FDTD simulation results of the relationship between △x and △λ in WGM resonators with different radii: (a) r0 = 3 µm; (b) r0 = 5 µm; (c) r0 = 10 µm.
Sensors 19 02798 g010
Figure 11. Relationship between △x/r0 and △λ/λ0 in WGM resonators with different initial radii.
Figure 11. Relationship between △x/r0 and △λ/λ0 in WGM resonators with different initial radii.
Sensors 19 02798 g011
Figure 12. Harmonic response simulation results of the MEMS resonator.
Figure 12. Harmonic response simulation results of the MEMS resonator.
Sensors 19 02798 g012
Figure 13. Transmission spectra drift of the MOEMS gyroscope.
Figure 13. Transmission spectra drift of the MOEMS gyroscope.
Sensors 19 02798 g013
Figure 14. The output of the MOEMS gyroscope in linear and nonideal conditions.
Figure 14. The output of the MOEMS gyroscope in linear and nonideal conditions.
Sensors 19 02798 g014
Figure 15. (a) A simplified model about cross-axis coupling in disk resonator. (b) The comparison of cross-axis coupling signal and drive signal.
Figure 15. (a) A simplified model about cross-axis coupling in disk resonator. (b) The comparison of cross-axis coupling signal and drive signal.
Sensors 19 02798 g015
Figure 16. The fabrication process of the MOEMS gyroscope.
Figure 16. The fabrication process of the MOEMS gyroscope.
Sensors 19 02798 g016
Table 1. Information about each layer in Figure 1e.
Table 1. Information about each layer in Figure 1e.
No.NameMaterialThickness(µm)
e-1Gap layerGlass100
e-2Sealing layerAu/Cr/Sn1
e-3Optical components layerSiNx0.35
e-4Metal pads layerAu/Cr0.3
e-5Optical substrate layerSiO22
e-6Device layerDoped silicon60
e-7Buried oxide layerSiO22
e-8Handle layerMonocrystalline silicon 500
Table 2. The parameters of the selected MEMS disk resonator.
Table 2. The parameters of the selected MEMS disk resonator.
ParameterValue
Spoke number16
Ring number60
Spoke width20 µm
Spoke length10 µm
Ring width20 µm
Electrode gap1.5 µm

Share and Cite

MDPI and ACS Style

Xia, D.; Huang, L.; Zhao, L. A New Design of an MOEMS Gyroscope Based on a WGM Microdisk Resonator. Sensors 2019, 19, 2798. https://doi.org/10.3390/s19122798

AMA Style

Xia D, Huang L, Zhao L. A New Design of an MOEMS Gyroscope Based on a WGM Microdisk Resonator. Sensors. 2019; 19(12):2798. https://doi.org/10.3390/s19122798

Chicago/Turabian Style

Xia, Dunzhu, Lingchao Huang, and Liye Zhao. 2019. "A New Design of an MOEMS Gyroscope Based on a WGM Microdisk Resonator" Sensors 19, no. 12: 2798. https://doi.org/10.3390/s19122798

APA Style

Xia, D., Huang, L., & Zhao, L. (2019). A New Design of an MOEMS Gyroscope Based on a WGM Microdisk Resonator. Sensors, 19(12), 2798. https://doi.org/10.3390/s19122798

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