The present invention relates to MEMS resonators, more particularly to MEMS piezoresistive resonators, and the driving thereof.
Micro-Electro Mechanical System (MEMS) resonators are known. MEMS resonators are formed from, for example, doped silicon, and are considered as a possible replacement for quartz crystal oscillators in applications of tuneable oscillators such as mobile telephones.
One particular type of MEMS resonator is a MEMS piezoresistive resonator. A MEMS piezoresistive resonator is a longitudinal mode resonator whose excitation is induced by electrostatic actuation and sensing of which is performed by sensing the piezoresistive effect of the doped material, typically silicon, forming the resonator. Such a MEMS piezoresistive resonator is described in WO 2004/053431, the contents of which are incorporated herein by reference.
The way in which a MEMS resonator functions is usually considered as corresponding to that of a lumped mass and a massless spring.
Conventionally, MEMS piezoresistive resonators are operated at their fundamental eigenmode. The fundamental eigenmode provides the fundamental resonant frequency f0 when considering only longitudinal vibration, i.e. in the direction of actuation, rather than bending. This ties in with the conventional understanding of the above mentioned consideration of the way a MEMS resonator functions, since a resonator comprising a lumped mass and a massless spring only has one resonant frequency, the fundamental resonant frequency f0.
The present inventor has realised that the physical effects involved in the operation of specifically a MEMS piezoresistive resonator allow the possibility to use higher order eigenmodes, i.e. additional eigenmodes that occur at higher frequencies than the frequency at which the fundamental eigenmode occurs. Such higher order eigenmodes may be considered as, or called, higher harmonics. The present inventor has further realised that analysis of the characteristics of the higher order eigenmodes allows a possibility for surprising performance gains to be achieved by determining or specifying various design features of an MEMS piezoresistive resonator in relation to the characteristics of the higher order eigenmodes.
In a first aspect, the present invention provides a method of operating a MEMS piezoresistive resonator; the method comprising driving the MEMS piezoresistive resonator at a higher order eigenmode than the fundamental eigenmode.
The method may further comprise arranging a route of flow of a sense current in relation to a characteristic of the higher order eigenmode.
The route of flow of the sense current may be arranged in relation to the characteristic of the higher order eigenmode by being at a point of maximum displacement of the higher order eigenmode.
The route of flow of the sense current may be arranged in relation to the characteristic of the higher order eigenmode by being at a point of maximum rate of change with respect to distance of displacement of the higher order eigenmode.
The route of flow of the sense current may be arranged by locating the end point of a trench of the MEMS piezoresistive resonator at a position determined in relation to the characteristic of the higher order eigenmode, the trench being formed between two beams of the MEMS piezoresistive resonator.
In a further aspect, the present invention provides a method of fabricating a MEMS piezoresistive resonator; the method comprising fabricating the MEMS piezoresistive resonator such that the route of flow of a sense current in operation is arranged in relation to a characteristic of a higher order eigenmode of the MEMS piezoresistive resonator compared to the fundamental eigenmode of the MEMS piezoresistive resonator.
The route of flow of the sense current may be arranged in relation to the characteristic of the higher order eigenmode by being at a point of maximum displacement of the higher order eigenmode.
The route of flow of the sense current may be arranged in relation to the characteristic of the higher order eigenmode by being at a point of maximum rate of change with respect to distance of displacement of the higher order eigenmode.
The route of flow of the sense current may be arranged by fabricating the MEMS piezoresistive resonator with a trench formed between two beams of the MEMS piezoresistive resonator, the end of the trench being located at a position determined in relation to the characteristic of the higher order eigenmode.
In a further aspect, the present invention provides a MEMS piezoresistive resonator comprising a trench formed between two beams, the end of the trench located at a position of maximum displacement or maximum rate of change with respect to distance of displacement of a higher order eigenmode of the MEMS piezoresistive resonator.
The MEMS piezoresistive resonator may further comprise a head that is wider than the width of the overall footprint of the two beams, and the end of the trench may be located at a point along the beams that is before where the beams reach the head.
In a further aspect, the present invention provides a MEMS piezoresistive resonator system; the system comprising: a MEMS piezoresistive resonator; and a driver for the MEMS piezoresistive resonator arranged to drive the MEMS piezoresistive resonator at a frequency of a higher order eigenmode than the fundamental eigenmode of the MEMS piezoresistive resonator.
Embodiments of the present invention will now be described, by way of example, with reference to the accompanying drawings, in which:
The beams and heads together form a space therebetween which we shall hereinafter call a trench 15.
The MEMS piezoresistive resonator 8 further comprises two actuation electrodes, namely a first electrode 18 and a second electrode 19. The first electrode 18 is positioned at the end of the first head 16 that is furthest from the beams, with a first gap 20 between the first electrode 18 and the end of the first head 16. The second electrode 19 is positioned at the end of the second head 17 that is furthest from the beams, with a second gap 21 between the second electrode 19 and the end of the second head 17.
In operation, excitation of the resonator is induced by electrostatic actuation and sensing is done by piezoresistive effect of the doped silicon beams. In more detail, both an AC voltage and a DC voltage are applied to the electrodes 18, 19, to drive the resonator structure into resonant vibration. A sense current 22 is sent through the beams 11, 12, 13, 14 which form the resonator arms, via the anchors 9, 10, as shown in
It can be seen in
The following dimensions are labelled. W indicates the width of a beam 11, 12, 13, 14. L indicates the length or distance between the anchors 9, 10 and the heads 16, 17, which is therefore the length the trench 15 would be (starting from the anchors) if it extended all the way to the head. R indicates the distance from the anchors 9, 10 to the end of the trench 15. The length of the heads 16, 17 is indicated by “a”. The width of the heads 16, 17 is indicated by “b”. The variable “x” indicates the distance from the anchors in the direction of the head.
In this embodiment the above mentioned dimensions of the MEMS piezoresistive resonator 8 have the following values:
The width of the trench 15, i.e. the spacing between the beams, is 1 μm. The thickness of the material of the MEMS piezoresistive resonator 8 is 1.5 um, the material being silicon (on insulator). Relevant material properties of the silicon forming the MEMS piezoresistive resonator 8 are a Young's modulus of 181 GPa and a density of 2329 kg/m3.
The frequency of the fundamental eigenmode of the above described MEMS piezoresistive resonator 8 is 56.5 MHz. However, in this embodiment, rather than being operated at its fundamental eigenmode, the MEMS piezoresistive resonator 8 is operated at a higher order eigenmode (higher harmonic), in particular at its second eigenmode, which in the case of the above described MEMS piezoresistive resonator 8 occurs at a frequency of 211.14 MHz. In an oscillator, the MEMS piezoresistive resonator 8 is driven at the required frequency, being the determined frequency of the selected higher eigenmode, by any suitable means, for example by suitable choice of the feedback amplifier's characteristics, as can be routinely implemented by the skilled person in the field of oscillator circuit design. An overview of such higher order eigenmode (harmonic) driving, and in particular how the choice of the above dimensions is inter-related with this, will now be given, including by reference to
The resonator provided by MEMS piezoresistive resonator 8 can be considered to be a clamped rod with varying cross-sectional area, A(x). For such a rod, compressional (as opposed to bending) eigenfrequencies can be calculated. Standing waves are found, that match the boundary conditions of this system. Other names for these vibrations are bulk acoustic waves or standing pressure waves. In the case of solid material, these are the same as the longitudinal eigenvibrations of the resonator. Using a one-dimensional, continuous model where only design lengths and material properties are used, the eigenfrequencies and corresponding modeshapes can be found analytically.
In
The first eigenmode 31, i.e. at the fundamental frequency f0, occurs at a frequency of 56.5 MHz. The second eigenmode 32, i.e. the second harmonic, occurs at a frequency of 211.14 MHz; and the third eigenmode 33, i.e. the third harmonic, occurs at a frequency of 365.55 MHz. A vertical line 38 at approximately x=17.75 μm indicates the interface between the beams and the head, corresponding to the fact that L=17.75 μm as described earlier above with reference to
The fundamental mode 31 describes, analogous to the mass-spring model, a linear slope of u(x) in the spring-part 40 and a much flatter u(x) in the head part 42. The strain in the material of the MEMS piezoresistive resonator 8 follows from the displacement, i.e. strain=du/dx. A massless spring would stretch with a constant strain all over the spring. A typical standing wave would be a non-constant strain.
The higher-order eigenmodes, i.e. harmonics, show that alternating compression and extension occur in both the spring part 40 and in the mass part 42. The spatial period of u(x) for any given eigenmode is the same in the spring part 40 as in the mass part 42, since the material and hence the material properties are the same in the spring part 40 as in the mass part 42. However, the amplitude of u(x) of any given eigenmode is different in the spring part 40 compared to in the mass part 42. It is believed by the inventor that this may be explained by the concept that in effect the heavier mass part 42 “swings-up” the lighter spring part 40. It is further believed by the inventor that the different eigenmodes each in effect balance harmonic deformation energy in both the spring part 40 and the head part 42.
The inventor has considered the different eigenmodes in terms of their potential for improved operation of the MEMS piezoresistive resonator 8. In particular, the inventor has considered various specific displacements shown in
Such increased displacement is used for improved performance, particularly improved gain. In particular, this is achieved in this embodiment by specifying the length of the trench 15 to fit to the x position where the maximum displacement 50 occurs, so that the sense current 22 passes across the MEMS piezoresistive resonator 8 at this location, and by performing resistive read-out at the frequency=211.14 MHz of the second eigenmode 32 i.e. the second harmonic.
A further refinement is that instead of considering simply displacement, including maximum displacement, one may consider slope of the displacement (with respect to x). Further details of this possibility are given below immediately following description of a further embodiment that includes a “straight MEMS piezoresistive resonator”.
The above described first embodiment is merely one example of MEMS piezoresistive resonator 8 shape and dimensions, and choice of eigenmodes, and location of trench etc. that can be used. It will be appreciated by the skilled person that a wide range of other specific embodiments may be arrived at by considering and analysing different shapes of MEMS piezoresistive resonator 8 and carrying out trial or modelling predictions to assess the use of different higher orders of eigenmodes and/or different trench lengths etc. Further background considerations and analysis will now be described with reference to
To calculate the frequencies at which stable standing pressure waves will occur in the dogbone resonator, the resonator is modelled as a rod with varying cross-section. The cross-sectional area A(x) leads to both the mass-per-meter m(x) and the stiffness EA(x) of the rod being discontinuous functions. The rod can be modelled as two connected rods, with conditions at both outer ends and conditions at the interface of the two. First, we derive the general differential equation that governs the standing waves in any rod.
At resonance, every single part of the material is displaced periodically, where a displacement U(x,t) can be expressed as
U(x,t)=u(x)g(t), with g(t)=cos(ωt) (1)
A general displacement profile will cause a strain
so that at any cross section a force
F(x,t)=EA(x)∈(x,t), (3)
where E is Young's modulus of elasticity, is acting on the surface. Cutting the rod at locations x2 and x1 therefore results in a resultant force acting on the isolated mass from x1 to x2. By writing
the ‘forces at two surfaces’ can be seen to be the result of integrating the derivative of their cause. This will prove to be useful when the force balance is completed.
A resultant force will accelerate a piece of mass, according to Fres=ma, where m is a mass and a is the acceleration. For an isolated piece of material in our vibrating rod, in between x1 and x2, mass and acceleration are defined as the integral
where the material density ρ is used. The expressions (4) and (5) have to balance, resulting in
This balance should hold for any part of the rod-in-resonance, so we can drop the integration limits to arrive at
on which separation of the variables t and x can be performed. By using U=u(x) g(t) and by noting that A(x) is non-constant, so that the product rule applies, we can write (by denoting the derivative in x with “'” and by denoting the derivative in t with “.”) that
since when
g=cos(ωt), then {umlaut over (g)}=−ω2 cos(ωt) (9)
Dividing both sides of Equation (8) by g results in
E[AU′]′=−ω2ρAU, (10)
where ω is still unknown.
The general formulation of standing waves or axial vibration of a rod with varying cross section is now written as
from which eigenfrequencies ω and corresponding modal shapes u(x) can be found. This second-order differential equation holds for the rod as a whole, as well as for parts of it. Because of the discontinuous A(x), see
Two second order equations require four boundary conditions, but since ω is also unknown, we need five. The conditions are summarized as
The last boundary condition follows from the fact that for the rod as a whole, the equation on both sides has to be continuous and differentiable. The mentioned scaling property stems from the fact that if u is a solution, then also 2u or 3u are solutions. To find the correct relative amplitudes of all found modes—there are infinite values ω and infinite functions u(x) that satisfy the equation—the modes will be normalized to the mass-function, according to
resulting in a normalization factor ni for the mode with index i.
For the two parts of the rod with cross sections A1 and A2, the general equation can be written as
Our solutions before applying the boundary conditions are
u
1(x)=c1 cos(βx)+c2 sin(βx), and (16)
u
2(x)=c3 cos(βx)+c4 sin(βx). (17)
The first boundary condition implies that
u
1(O)=O→c1=O, (18)
The scaling property described before is used to state
c2=1 (19)
since one of the coefficients or amplitudes can be chosen. The strain-free end implies
−c3 sin [β(L+a)]β+c4 cos [β(L+a)]β=0, so
c
4
=c
3 tan [β(L+a)]. (20)
This is the first of three relations between c3, c4 and β. Using this expression for c4, the continuity boundary condition can be expressed as
c
3=cos [β(L+a)] sec [aβ] sin [Lβ]. (21)
Finally, the condition of differentiability, together with the expressions for c3, and c4, allows us to find values for β by
Equation (22) can be re-written in other periodic functions, but the general statement can be made, that solutions for the eigenfrequencies ω correspond to the values of β at which two harmonic functions of different period intersect.
Using a found value of β, the values for c3, and c4, can be calculated. With c1=0 and c2=1, we can now present the modes corresponding to the eigenfrequencies. To scale these eigenfunctions properly with respect to each other, a normalization is undertaken, according to
This results in normalization factors n(i) being
[n(i)]2=2β/dp[−ω(−2Lβ+sin [2Lβ])+b sin [Lβ]2(aβsec[aβ]2+tan [aβ])], (24)
where d is the resonator thickness. The modeshapes are
where c3, and c4 follow from Equations (21) and (20) respectively.
Analytical expressions for the eigenfrequencies as well as the eigenfunctions are derived. In the extreme cases where either a or L is reduced to zero, and the resonator reduces to a conventional rod of constant cross section, or when simply A1=A2, all expressions remain valid. For such a rod, literature provides expressions that have the same results for ω and u(x) as derived here.
The mode shapes can also be normalized to unity-at-the-edge. This normalization is defined as
q
(i)
c
3 cos [β(i)(L+a)]+c4 sin [β(i)(L+a)]=1, (26)
which is not necessary for the uniform rod, since these modes are all already found to end at unity. The reason for this normalization, however, lies in the practical field of determining the actuation gap size through the effective spring constant. Based on a found value of β, the normalization factors are
g
(f)=cos [aβ] csc [Lβ]. (27)
With the normalization q(i) known, an expression for the effective spring constant keff can be derived.
Consider the fundamental resonance at ω0 with corresponding mode shape u(x). At resonance, Rayleigh's quotient (see L. Meirovitch, “Elements of vibration analysis”, McGraw-Hill, 1986) is defined as
where Vmax is the maximum potential energy and T the reference kinetic energy. These values scale with vibration amplitude. A mass-spring system with a force F acting on the outer edge of the rigid mass, would induce a potential energy of k
Rayleigh's quotient R(u) is written as
where u(L+a)=1. The kinetic energy is by definition
where again u(L+a)=1. In order to have the found functions u1 and u2 meet this unit amplitude at the end, they are multiplied by the normalization factor q(i), in this case q0. The energy balance of Equation (29) can now be written as
so that an effective spring constant can be calculated. The dimensions are checked to be [Nm−1] as for a spring. For the general dogbone, the integral is split in two by writing
resulting in
The effective mass meff follows from the natural frequency and the found spring constant, using q0 based on β0 as in Equation (27).
The resonant frequencies, the shapes of the standing waves, the normalization factors and the effective spring constant all rely on found values of β. As stated, these values can be found from
For this equation, there is no mathematical expression to provide all intersections analytically. Next, we will present two different approximations to find only the first intersection.
The first approximation is a very accurate expression. However, this accurate expression is valid only when A1/A2≈1, but was tested to hold up to a ratio of 4, or ¼ alternatively. For most dogbone resonators, this is sufficient, but the generality for making either L or a zero, is lost, unfortunately. The first intersection, hence the value of β0 is found at
Basically, a tan-function is fitted, of which the squared function approximates the product of the two different tan-functions in the original equation.
The equation for the intersections can also be written as
which can be seen as a long-period cosine starting just below 1, intersecting a short-period cosine starting just above 1, as depicted in
The value β0 can be approximated again from this expression, but cannot be found analytically. When both cosine-functions are approximated by a Taylor series expansion around β=0, the second order (or parabola) expressions yield an intersection at
which is very similar to the expressions found in other studies. The first resonant frequency is then
with ω0, in radians per second.
Although exact expressions are not available, the two given approximations are useful. The first gives a very good estimate for a not-too-extreme dogbone and the second gives an estimate that is worse, but in formulation very convenient to use in e.g. designing a resonator for a certain operating frequency. Since for both approximations the un-approximated functional description is available, very good approximations can be found after a few iteration steps of a Newton-Rhapson root finding procedure.
Given an estimate xn, an improved estimate xn+1 can be found, according to
Here, function f(x) is a differentiable function of which the zero-crossings or roots are sought. Close to a root, this iteration will, generally, double the amount of correct digits in the estimate each time.
The resonator is part of an oscillator. To oscillate, the output of the resonator is fed back with amplification to drive it again. The actuation is electrostatic, but the measured signal is obtained by using the piezoresistive properties of the material. Due to the piezoresistive effect, the material experiences a relative change in resistance, according to
Where K is the so-called gauge factor, which is a material property. A DC current through the spring-arms of the resonator will provide a voltage at the output that is modulated due to a changing ∈(t) at resonance.
Modeshapes u(x) for the dogbone resonator have been derived. The strain is ∈(x)=du/dx, so that the piezoresistive effect, related to the total encountered strain, for the measurement current in the dogbone resonator is
where x′ is the point of return, i.e. the point as explained in the above embodiment where the sense current is made to cross the MEMS piezoresistive resonator, e.g. by location of the end of the trench at that point. If the current follows a path that is a result of some implantation, the multiplication by a factor of 2 can be possibly be increased to 4 or 6. But for now, as the point x′ is defined by the length of the trench, we can observe that for the present dog bone, the values
u
(1)(L)=n(1) sin [β(1)L] and u(2)(L)=n(2) sin [β(2)L], (42)
which are the normalized displacements at x′=L for modes 1 and 2, differ. This means that the to be measured resistive effect is larger, when the present dogbone resonator is excited at the second resonance.
When an optimized dogbone resonator is manufactured, the resistive effect can be amplified even further by reducing the length of the trench. The value for x′, which is the length of the trench, so that 0<x′<L, is found as the point where u(x) has the first maximum. It can be found as the first zero of the derivative, resulting in
Neglecting the amplitude, the value for x′ is found to be
For each resonant mode, the optimum x′ can be calculated.
The function ∈(x)=du/dx shows that the highest stresses, since σ=E∈, will occur in the first part of the resonator, since, generally, the amplitudes of u(x) are higher in the first part. The period of the standing wave is related to β and is equal in both parts of the resonator. Because u(0)=0, the highest stress can be stated to occur at x=0. Driving the resonator in higher modes than the fundamental mode, causes higher stresses, as the slopes in u(x) get larger. The maximum stress will be larger, compared to the fundamental mode, by a factor of
where the cos-factor indicates the use of the derivative of u1.
For the dogbone under consideration, while assuming that the trench would have an optimized length, this could provide a nearly 4 times higher strain signal, when exciting it in the third mode rather than in the first one. For this mode the trench length should be x′=0.53 L. The signal gain of a factor of 4 comes with a gain of nearly 10 in the maximum stress that occurs.
This has now completed the above mentioned theoretical account provided with reference to
It will be appreciated from the above description of the first embodiment and also the theoretical account that the skilled person may analyse and select different eigenmodes as required, according to the details of a specific MEMS piezoresistive resonator under consideration. For example, although in the above described first embodiment the second eigenmode gave the largest gain, for other MEMS piezoresistive resonators of different shapes and dimensions other eigenmodes may provide the largest gain.
Furthermore, the skilled person may decide to implement other embodiments where a given eigenmode is used even if it is not the eigenmode that provides the highest gain. Yet further, the skilled person may decide to implement embodiments where the sense current is made to cross the MEMS piezoresistive resonator at a point that does not represent the highest gain for the chosen eigenmode.
The skilled person may use any appropriate analysis approach (or combination of analysis approaches), for example employing trial and error, to arrive at a chosen eigenmode and/or current crossing location. For example, mathematical modelling may be employed, or actual MEMS piezoresistive resonators may be built and driven at different frequencies etc. on a trial and error basis.
In the above embodiments, by virtue of locating the end of the trench before the start of the head, the remaining part of the beams comprises in fact a single wide portion comprising the width of the two beams plus the width of the trench. However, this need not be the case. For example,
In the above embodiments the MEMS piezoresistive resonator is a so-called dogbone resonator, i.e. is shaped with a head (mass) part at each end that is wider than the beam (spring) part. However, in other embodiments other shapes of resonator may be used, and their higher order eigenmodes exploited. For example, a straight MEMS piezoresistive resonator may be used, i.e. one in which the overall footprint taken up by the beams is the same width as that of the heads. In fact, with such a shape, there is effectively no distinction between a mass and a spring as such; longitudinal vibration modes are continuous standing pressure waves of uniform amplitude throughout the body. Further details of a suitable structure of a straight MEMS piezoresistive resonator are given in WO2004/053431, the contents of which are incorporated herein by reference.
In the above description, the aspect of improved displacement has been concentrated on. However, another possibility is to consider slope of the displacement (with respect to x). The derivative of displacement is strain, and increasing this increases the output signal, i.e. gain. This applies to both dogbone MEMS piezoresistive resonators and straight MEMS piezoresistive resonators. That is, straight MEMS piezoresistive resonators are capable of higher gain using eigenmodes other than the fundamental eigenmode even though it would not particularly be expected to exhibit higher maximum displacement along the body compared to its fundamental mode, since it does nevertheless exhibit steeper derivatives when using eigenmodes other than the fundamental eigenmode.
Such consideration of derivative of displacement leads to further possibilities for optimising or otherwise choosing the position where the sense current should cross the MEMS piezoresistive resonator, in addition or as an alternative to the above described choice made to optimise displacement. Up until this position the derivative is positive, hence a contribution to the resistive readout is made along the walk from origin to the above selected sense current crossing point. However, one may further consider that besides wanting to run the current through regions of high derivative (strain), it is also desirable to keep the current path as short as possible, to reduce the overall resistance along the sense path. If this is of interest, then one possibility for defining the best sense current crossing point may be considered as “in between the origin and the point of maximum displacement”, since the strain gradually decreases to be zero at the point of maximum displacement. Note how this differs to the prior art use of fundamental eigenmode where the strain is constant all along the beams, and hence the choice of sense current crossing point has no effect.
In the above embodiments the trench 15 (and where included the further trench space 15′) comprises a low-pressure gas or vacuum, the MEMS piezoresistive resonator being encapsulated. Another possibility is that the trench 15 (and where included the further trench space 15′) is air. In other embodiments the trench 15 (and where included the further trench space 15′) may be any suitable material, such as a dielectric material, whose electrical conductivity is lower than that of the beams, such that the sense current follows a route that excludes the trench.
In further embodiments, the materials and form of the MEMS piezoresistive resonator may be varied compared to the above described embodiments. For example, the MEMS piezoresistive resonator need not be “dogbone” shaped nor straight in that any other suitable shapes may be used. Furthermore, the beams and anchor may be replaced by alternative structures. Yet further, a MEMS piezoresistive resonator may be used that has less symmetry than the one described above, for example there may be a head and beams on only one side of the anchors.
Number | Date | Country | Kind |
---|---|---|---|
07109409.8 | Jun 2007 | EP | regional |
Filing Document | Filing Date | Country | Kind | 371c Date |
---|---|---|---|---|
PCT/IB08/52088 | 5/28/2008 | WO | 00 | 11/20/2009 |