Figure of Merit Enhancement of Laterally Vibrating RF-MEMS Resonators via Energy-Preserving Addendum Frame

This paper examines a new technique to improve the figure of merit of laterally vibrating RF-MEMS resonators through an energy-preserving suspended addendum frame structure using finite element analysis. The proposed suspended addendum frame on the sides of the resonant plate helps as a mechanical vibration isolator from the supporting substrate. This enables the resonator to have a low acoustic energy loss, resulting in a higher quality factor. The simulated attenuation characteristics of the suspended addendum frame are up to an order of magnitude larger than those achieved with the conventional structure. Even though the deployed technique does not have a significant impact on increasing the effective electromechanical coupling coefficient, due to a gigantic improvement in the unloaded quality factor, from 4106 to 51,136, the resonator with the suspended frame achieved an 11-folds improvement in the figure of merit compared to that of the conventional resonator. Moreover, the insertion loss was improved from 5 dB down to a value as low as 0.7 dB. Furthermore, a method of suppressing spurious mode is demonstrated to remove the one incurred by the reflected waves due to the proposed energy-preserving structure.


Introduction
RF-MEMS resonators have been used for many applications, ranging from acoustic filters used for communication systems to timing application, owing to their excellent power handling capability and low motional resistance [1][2][3][4][5]. These merits are associated with piezoelectric transduction, which offers strong electromechanical coupling, and the use of substrate materials that have a low intrinsic damping and high energy density, such as single-crystal silicon (SCS) [1].
Amidst the most attractive properties of RF-MEMS resonators are their compatibility with existing CMOS technology, making them highly attractive for monolithic integration with CMOS circuitry [6][7][8][9][10]. In most of their applications, resonators with a high Q are of great importance in fulfilling the stringent requirements of the device. Accordingly, in applications for oscillators, Q sets the phase noise performance [8,11]. In filters, a high Q reduces the insertion loss while narrowing the bandwidth [12][13][14], whereas, in sensing applications it benefits the sensing resolution [15][16][17]. However, in order to obtain the aforementioned requirements, the Q of piezoelectric resonators still needs to be improved. Another important parameter which expresses the electromechanical transduction efficiency of a resonator is the effective electromechanical coupling coefficient (k 2 e f f ) [18][19][20]. The product of these two parameters, as given in Equation (1) and used as a figure of merit (FoM), is an important parameter for judging how good a resonator is. For this, designing resonators with a higher figure of merit is a hot topic for research. Hence, numerous studies have made efforts to enhance the FoM of a resonator by improving one or both of Q and k 2 e f f .
Even though there are many intrinsic and extrinsic loss mechanisms, such as dielectric loss, piezoelectric loss, surface loss, and Ohmic loss, it has been widely reported that anchor loss is the most dominant loss mechanism of MEMS resonators severely limiting the quality factor effect in laterally vibrating RF MEMS resonators [21][22][23][24][25][26]. This loss occurs due to the elastic wave propagating out from the resonant cavity to the supporting substrate through the supporting tether. The primary design approach to reduce this loss is to locate the tether at the nodal points of the resonator where there is a minimum vibration, reducing the width of the tether and using a tether length of quarter an acoustic wavelength [27][28][29][30]. However, a limit to reduce the tether dimension for practical feasibility is reported as the pitfall of these approaches. In Ref. [31], our group demonstrated the effect of swastika hole-shaped phononic crystals (SW-PnC) with significantly high acoustic band gaps; these, on their placement on the anchoring boundaries, improve the quality factor by over 1.4-fold. Furthermore, in [32] the effect of placing 2 D-PnC on the tether and/or the anchoring boundaries of the resonator is demonstrated, and an improvement of about 3-fold in Q is reported. Another approach reported recently is deploying acoustic reflectors on the anchoring boundaries of the resonator, aiming to reflect some of the radiated elastic waves back to the resonant plate [33]. The major pitfall of these techniques is the limited enhancement in Q, as the reflection of acoustic waves is finished once it reaches the anchoring boundaries or at least after it passes the supporting tether. Furthermore, a MEMS support design to reduce the anchor loss of PZT-on-Si disc resonators by using a mechanical impedance transformer design on the tether is reported in [34]; this achieved a significant improvement in reducing the anchor loss.
This paper aims to introduce a new approach for improving FoM by using a suspended energy-preserving addendum frame. We firstly describe the layout of suspended energy-preserving addendum frames along with the mechanical equivalent model and transmissibility analysis in comparison with that of conventional resonator structure. Subsequently, the design and simulation of resonators with the proposed energy-preserving features and the potential improvements in the Q are explored through a series of simulations. Afterwards, frequency domain analysis is performed for all resonators, and the electrical responses are thoroughly discussed, followed by a description of the spurious mode mitigation technique introduced in this study. The finite element package COMSOL Multiphysics was used to calculate the eigen modes and frequency responses of all the resonators designed. A loss factor (η εS ) of 0.0003 was applied in all the simulations for the piezoelectric material (AlN) used for this work.

Mechanical Equivalent Model
As is described prior, the main source of dissipation in RF-MEMS resonators is the propagation of acoustic energy to the anchoring boundaries through the supporting tether, which leads to a decrease in Q. Aiming to isolate the acoustic vibration between the resonant plate and the supporting substrate and taking the idea of mechanical vibration isolators as a motivation, a suspended addendum frame is added to the conventional structure along with a novel meandering tether support. Figure 1a,b show the quarter plate of the conventional MEMS resonator and the modified resonator with the suspended addendum frame, hereafter named the energy-preserving frame (EPF), designed to vibrate at a frequency of 140.32 MHZ in the first fundamental mode. Figure 1c,d show the schematic view of a mass-spring-damper equivalent model for the conventional and EPF MEMS resonators, respectively. The conventional resonator can be represented by a simple massspring-damper system of a single mass: where m rp , x 1 , c 1 , k 1 , and F 0 stand for the mass of the resonant plate, the total equivalent displacement, the damping coefficient, the stiffness, and the applied force, respectively.
Assuming that the quality factor of the resonator is mainly dependent on the anchor loss, it is reasonable to neglect the effect of damping. In that case, the displacement x 1 can be expressed as: where r = ω ω n and ω n = k 1 /m rp . The force transmitted to the fixed support can be expressed as:

Transmissibility Analysis
To further demonstrate the acoustic wave attenuation capability of the proposed method, a transmission line simulation was performed for the EPF and conventional resonators. As can be seen in Figure 1a,b, an input frequency-dependent harmonic displacement excitation ( ) with an amplitude of 1 nm was applied at the boundary probe 1. Then, the vibrational transmissibility could be sensed at boundary probe 2 as a harmonic displacement of ( ). The frequency-dependent transmission (dB) could be obtained by the relation [30]: As can be seen from Figure 2, at the resonance frequency of 140.32 MHz the simulated attenuation for the conventional and proposed EPF resonators were 20 dB and 124 dB, As a result, the ratio of transmitted force to the excitation force is called transmissibility, T, which is expressed as: On the other hand, the resonator EPF can be represented by the coupling of two masses: one for the resonant plate (m rp ) and the other for the suspended addendum frame (m sa ).
where x i , k i , and c i stand for the displacement, stiffness, and damping of conventional and EPF resonators. With the same context to the single mass system, neglecting the damping effect, the motions of the masses m rp and m sa to an input force F eq are: The transmissibility of the energy isolator can then be expressed as: Taking the law of conservation of energy into consideration and assuming that excitation forces are equal (F 0 = F eq ) and of an equal stiffness (k 1 = k 2 ), as the vibration of the suspended mass is caused by the resonant body, the displacement x 2 is smaller than x 1 . As a result, the transmissibility of the resonator with a suspended frame is expected to be much lower than that of a conventional resonator (i.e., higher isolation capacity).

Transmissibility Analysis
To further demonstrate the acoustic wave attenuation capability of the proposed method, a transmission line simulation was performed for the EPF and conventional resonators. As can be seen in Figure 1a,b, an input frequency-dependent harmonic displacement excitation U in ( f ) with an amplitude of 1 nm was applied at the boundary probe 1. Then, the vibrational transmissibility could be sensed at boundary probe 2 as a harmonic displacement of U out ( f ). The frequency-dependent transmission (dB) could be obtained by the relation [30]: As can be seen from Figure 2, at the resonance frequency of 140.32 MHz the simulated attenuation for the conventional and proposed EPF resonators were 20 dB and 124 dB, respectively. This implies that the amount of energy transmitted to probe 2 (which is on the anchoring substrate of the resonators) for the conventional resonator was orders of magnitude greater than that of the EPF resonator. The large improvement in the attenuation of the EPF resonator was consistent with the displacement profile of the resonators in the anchoring boundaries seen in Figure 1a,b, among which the conventional resonator had a larger displacement, implying a larger amount of energy leaking from the resonant plate compared to the one with the EPF resonator.

Device Design and Simulation
Cognizant to the resonators designed above for transmissibility analysis, all the resonators designed henceforth were set to vibrate in lateral mode with a frequency of about 140 MHz. The resonant plate was constructed from a piezoelectric film of aluminum nitride (AlN) sandwiched with a 10 µm anisotropic single-crystal silicon (SCS) substrate, which is also used as a common ground, and 1 µm-thick aluminum (Al) used as an interdigitated electrode for the input and output terminals. Furthermore, each of the resonators were designed to be transduced in the seventh order width-extensional (WE) mode. The seventh-order resonant frequency ( ) of the width extension mode resonators was determined by the center-to-center distance of the electrodes (Pitch), which corresponds to 1/7 of the total resonant plate width ( ), as can be seen from Figure 3. In other words, the nth order resonant frequency ( ) is n times the fundamental mode exhibited by the overall resonant plate.
where ν is the velocity of the acoustic wave and n is the harmonic number that usually represents the number of IDT fingers (i.e., unless some portion of the resonant plate is left inactive-for example, such as reducing the loading effect of electrodes). Considering the much thicker single-crystal silicon used in this work compared to that of the piezoelectric material (AlN) and the electrode metal (Al), the vibration property of the device is mainly dominated by the single-crystal silicon. In this example, an acoustic velocity of 8500 m/s is used for analysis (i.e., acoustic velocity in anisotropic single-crystal silicon, <100> orientation). Despite the addition of some suspended structures (i.e., a frame) for the sake of improving the performance of the resonators, as is briefly discussed later, the dimensions of the resonant plates were consistent across all the resonators designed. In this regard, the length of the resonators ( ), and the width of resonators ( ), were designed with the proportion of 3.5 (i.e., =3.5× ).
It is noteworthy that, in order to reduce the computational time for the 3D FEA simulation of the resonators, only a quarter of the resonators were simulated and the symmetric condition was applied at the symmetric sides, as seen in Figure 3, to ensure that the eigen modes were solved for the whole resonant plate. Apart from this, a perfectly matched layer (PML) with a width of 3λ was applied to absorb the energy that escaped

Device Design and Simulation
Cognizant to the resonators designed above for transmissibility analysis, all the resonators designed henceforth were set to vibrate in lateral mode with a frequency of about 140 MHz. The resonant plate was constructed from a piezoelectric film of aluminum nitride (AlN) sandwiched with a 10 µm anisotropic single-crystal silicon (SCS) substrate, which is also used as a common ground, and 1 µm-thick aluminum (Al) used as an interdigitated electrode for the input and output terminals. Furthermore, each of the resonators were designed to be transduced in the seventh order width-extensional (WE) mode. The seventhorder resonant frequency ( f r ) of the width extension mode resonators was determined by the center-to-center distance of the electrodes (Pitch), which corresponds to 1/7 of the total resonant plate width (W r ), as can be seen from Figure 3. In other words, the n-th order resonant frequency ( f r ) is n times the fundamental mode exhibited by the overall resonant plate.
where ν is the velocity of the acoustic wave and n is the harmonic number that usually represents the number of IDT fingers (i.e., unless some portion of the resonant plate is left inactive-for example, such as reducing the loading effect of electrodes). Considering the much thicker single-crystal silicon used in this work compared to that of the piezoelectric material (AlN) and the electrode metal (Al), the vibration property of the device is mainly dominated by the single-crystal silicon. In this example, an acoustic velocity of 8500 m/s is used for analysis (i.e., acoustic velocity in anisotropic single-crystal silicon, <100> orientation). Despite the addition of some suspended structures (i.e., a frame) for the sake of improving the performance of the resonators, as is briefly discussed later, the dimensions of the resonant plates were consistent across all the resonators designed. In this regard, the length of the resonators (L r ), and the width of resonators (W r ), were designed with the proportion of 3.5 (i.e., L r = 3.5 × W r ).
Micromachines 2022, 13, x FOR PEER REVIEW 6 of 14 from the resonator and avoid reflections at the boundaries. Table 1 lists the values of basic design parameters for the resonators.

AlN
Si AL  To fully demonstrate the ability of the proposed method to meet the stated objectives, two different devices were designed. Device one was a conventional resonator, as shown in Figure 4a, along with the corresponding displacement mode shape in the y-field. The second device was a resonator with a suspended energy-preserving addendum frame (EPF), as shown in Figure 4b. It can be seen in Figure 4a that the rectangular resonant plate of a conventional resonator was attached to the undercut by supporting tethers, which in turn were directly attached to the PML. Whereas, in the case of the EPF resonator, the resonant plate was not directly connected to the undercut. Rather, it was connected through the suspended addendum via a meandering tether with a width of ( ), which It is noteworthy that, in order to reduce the computational time for the 3D FEA simulation of the resonators, only a quarter of the resonators were simulated and the symmetric condition was applied at the symmetric sides, as seen in Figure 3, to ensure that the eigen modes were solved for the whole resonant plate. Apart from this, a perfectly matched layer (PML) with a width of 3λ was applied to absorb the energy that escaped from the resonator and avoid reflections at the boundaries. Table 1 lists the values of basic design parameters for the resonators. To fully demonstrate the ability of the proposed method to meet the stated objectives, two different devices were designed. Device one was a conventional resonator, as shown in Figure 4a, along with the corresponding displacement mode shape in the y-field. The second device was a resonator with a suspended energy-preserving addendum frame (EPF), as shown in Figure 4b. It can be seen in Figure 4a that the rectangular resonant plate of a conventional resonator was attached to the undercut by supporting tethers, which in turn were directly attached to the PML. Whereas, in the case of the EPF resonator, the resonant plate was not directly connected to the undercut. Rather, it was connected through the suspended addendum via a meandering tether with a width of (W t2 ), which is half of the width of the tether (W t1 ) (i.e., W t1 = W t /2, Figure 1). The width of the addendum (W sa ) was tuned with FEA simulation, with an optimum value equal to two times the pitch width. Moreover, the suspended addendum frame was placed in an equal distance (i.e., gap of addendum (G a )) from the resonant plate and the undercut. is half of the width of the tether ( ) (i.e., = 2 ⁄ , Figure 1). The width of the addendum ( ) was tuned with FEA simulation, with an optimum value equal to two times the pitch width. Moreover, the suspended addendum frame was placed in an equal distance (i.e., gap of addendum ( )) from the resonant plate and the undercut.     Figure 4c presents the vibration modes of a resonator in its fundamental mode and three overtones (third, fifth, and seventh overtones) for the conventional and EPF resonators. For the sake of understanding the effect of vibration mode on the quality factor of the resonator, the same rectangular plate size was used. The fundamental mode along with three consecutive overtones with frequencies of 20 MHz, 60 MHz, 99.8 MHz, and 140.3 MHz, respectively, were analyzed with FEA modal analysis. Amongst the conventional resonators, the largest overtone in the list (seventh overtone) experienced the largest loss, resulting in the minimum Q anchor (6740). Whereas, it was from the seventh overtone of the EPF resonator that the maximum Q anchor (97,575) was achieved. Consequently, further analysis was conducted for this mode in order to find the impact of the proposed structure on the performance parameters.

Max
From the mode shapes (displacement distribution and y-field) of the resonators presented in Figure 4a,b, it is apparent that the displacement in the supporting tether of the conventional resonator was much larger than that of the other three resonators, suggesting that the energy dissipation was higher (i.e., the color scale represents maximum displacement field with dark red and minimum displacement with dark blue).
Furthermore, the total displacement plots of the resonant plate at the A-A' and C-C' lines (see Figure 3), as demonstrated in Figure 5a,c, reveal that the conventional resonator had the minimum displacement. This implies that the amount of energy on the resonant body is lower than the EPF resonator. On the other hand, the plots shown in Figure 5b,d at the lines B-B' and D-D' (see Figure 3), which are at the tether and the anchoring substrate, respectively, revealed that the conventional resonator has significantly higher displacements. This indicates that the energy leaking into the supporting substrate is higher than that in the EPF resonators. On top of that, the FEA simulated Q anchor results of 97,575 for the EPF resonator were attained, which account for an improvement as high as 14.5 fold that of the conventional resonator's result of 6740. loss, resulting in the minimum (6740). Whereas, it was from the seventh overtone of the EPF resonator that the maximum (97,575) was achieved. Consequently, further analysis was conducted for this mode in order to find the impact of the proposed structure on the performance parameters.
From the mode shapes (displacement distribution and y-field) of the resonators presented in Figure 4a,b, it is apparent that the displacement in the supporting tether of the conventional resonator was much larger than that of the other three resonators, suggesting that the energy dissipation was higher (i.e., the color scale represents maximum displacement field with dark red and minimum displacement with dark blue).
Furthermore, the total displacement plots of the resonant plate at the A-A' and C-C' lines (see Figure 3), as demonstrated in Figure 5a,c, reveal that the conventional resonator had the minimum displacement. This implies that the amount of energy on the resonant body is lower than the EPF resonator. On the other hand, the plots shown in Figure 5b,d at the lines B-B' and D-D' (see Figure 3), which are at the tether and the anchoring substrate, respectively, revealed that the conventional resonator has significantly higher displacements. This indicates that the energy leaking into the supporting substrate is higher than that in the EPF resonators. On top of that, the FEA simulated results of 97,575 for the EPF resonator were attained, which account for an improvement as high as 14.5 fold that of the conventional resonator's result of 6740.

Discussion
The FEA-simulated admittance response of the conventional resonator along with the EPF resonator is plotted in Figure 6a,b. The summaries of the results obtained are listed in Table 2, from which it can be observed that a k 2 e f f of 0.13% is obtained from the EPF resonator and one of 0.14% is obtained from the conventional resonator. This indicates that the employed techniques do not have a significant effect on the coupling coefficient. On the other hand, the simulated transmission S21 (dB) responses plotted in Figure 7a,b reveal a significant improvement in the insertion loss (IL) from 5 dB to 0.7 dB. The unloaded Q (Q u ) also shows a significant improvement of about 12.5-fold. On top of that, as shown in Table 2, the EPF resonator demonstrated a Q u * k 2 e f f product (i.e FoM) of 66) which shows an 11-fold improvement compared to the conventional resonator with a value of about 6.

Discussion
The FEA-simulated admittance response of the conventional resonator along with the EPF resonator is plotted in Figure 6a,b. The summaries of the results obtained are listed in Table 2, from which it can be observed that a of 0.13% is obtained from the EPF resonator and one of 0.14% is obtained from the conventional resonator. This indicates that the employed techniques do not have a significant effect on the coupling coefficient. On the other hand, the simulated transmission S21 (dB) responses plotted in Figure  7a,b reveal a significant improvement in the insertion loss (IL) from 5 dB to 0.7 dB. The unloaded Q ( ) also shows a significant improvement of about 12.5-fold. On top of that, as shown in Table 2, the EPF resonator demonstrated a * product (i.e FoM) of 66) which shows an 11-fold improvement compared to the conventional resonator with a value of about 6.

Discussion
The FEA-simulated admittance response of the conventional resonator along w the EPF resonator is plotted in Figure 6a,b. The summaries of the results obtained listed in Table 2, from which it can be observed that a of 0.13% is obtained from t EPF resonator and one of 0.14% is obtained from the conventional resonator. This in cates that the employed techniques do not have a significant effect on the coupling coe cient. On the other hand, the simulated transmission S21 (dB) responses plotted in Figu 7a,b reveal a significant improvement in the insertion loss (IL) from 5 dB to 0.7 dB. T unloaded Q ( ) also shows a significant improvement of about 12.5-fold. On top of th as shown in Table 2, the EPF resonator demonstrated a * product (i.e FoM) of which shows an 11-fold improvement compared to the conventional resonator with value of about 6.    The main pitfall encountered by the proposed FoM enhancement technique is that it introduces undesirable in-band and out-of-band spurious modes which mainly result from the reflected transverse waves. As these spurious modes highly affect the main mode, the study of the mitigation technique is necessary. In this regard, a new method of spurious reduction is introduced, as presented in the next section.

Spurious Mode Suppression
Spurious modes are considered as a major bottleneck in achieving high performance with RF-MEMS resonators. These unwanted spurious modes are mainly caused by reflections of acoustic waves by the acoustic boundaries and energy-preserving structures, in this particular case. As can be seen from the frequency response in Figure 6b, a significantly large amount of spurious mode appears before and after the main mode, which is not the case for the conventional resonator's frequency response, as shown in Figure 6a. That being so, an approach for suppressing spurious mode is devised by adding metal patches on the edges of the electrodes, as is illustrated in Figure 8.
The main pitfall encountered by the proposed FoM enhancement technique is that it introduces undesirable in-band and out-of-band spurious modes which mainly result from the reflected transverse waves. As these spurious modes highly affect the main mode, the study of the mitigation technique is necessary. In this regard, a new method of spurious reduction is introduced, as presented in the next section.

Spurious Mode Suppression
Spurious modes are considered as a major bottleneck in achieving high performance with RF-MEMS resonators. These unwanted spurious modes are mainly caused by reflections of acoustic waves by the acoustic boundaries and energy-preserving structures, in this particular case. As can be seen from the frequency response in Figure 6b, a significantly large amount of spurious mode appears before and after the main mode, which is not the case for the conventional resonator's frequency response, as shown in Figure 6a. That being so, an approach for suppressing spurious mode is devised by adding metal patches on the edges of the electrodes, as is illustrated in Figure 8. The metal patch on the edge of an electrode creates different acoustic impedance zones wherein the transverse spurious waves could be significantly influenced into having a vanishing coupling, resulting in spurious free responses. The metal patch installed has patch depth ( ) and patch length ( ) as indicated in Figure 8. Moreover, the patch width is fixed to be identical with that of the IDT finger (electrode) width ( ) for the sake of simplicity and its insignificant effect on the results.
The first investigation is performed on the selection mechanism of materials to be used as a patch. In this regard, five metals (i.e., Ag, Cu, Au, Pt and W) were investigated. According to the resulting admittance, Y11 (dB), and conductance, G (dB), responses presented in Figure 9a,b, the performance is in direct correlation with the acoustic impedance (Z) of the patch metals used, which can also be seen from Table 3 and the frequency responses. Consequently, tungsten (W) shows superior performance, and is used in subsequent analysis to determine the other parameters of the patch. The optimum dimensions of the patch for the best suppression of spurious modes are tuned based on FEA simulations as it is demonstrated by the admittance, Y11 (dB) and conductance, G (dB) responses in Figure 10a  The metal patch on the edge of an electrode creates different acoustic impedance zones wherein the transverse spurious waves could be significantly influenced into having a vanishing coupling, resulting in spurious free responses. The metal patch installed has patch depth (H patch ) and patch length (L patch ) as indicated in Figure 8. Moreover, the patch width is fixed to be identical with that of the IDT finger (electrode) width (W e ) for the sake of simplicity and its insignificant effect on the results.
The first investigation is performed on the selection mechanism of materials to be used as a patch. In this regard, five metals (i.e., Ag, Cu, Au, Pt and W) were investigated. According to the resulting admittance, Y11 (dB), and conductance, G (dB), responses presented in Figure 9a,b, the performance is in direct correlation with the acoustic impedance (Z) of the patch metals used, which can also be seen from Table 3 and the frequency responses. Consequently, tungsten (W) shows superior performance, and is used in subsequent analysis to determine the other parameters of the patch. The optimum dimensions of the patch for the best suppression of spurious modes are tuned based on FEA simulations as it is demonstrated by the admittance, Y11 (dB) and conductance, G (dB) responses in Figures 10a,b and 11a        (a) (b) Figure 10. Illustration of simulated electrical responses using different patch lengths (i.e., 4 ⁄ , 3 ⁄ , 2 ⁄ and ). The patch depth ( ) of 2 ⁄ and patch metal of tungsten (W) is used for the simulation. (a) Admittance, Y11 (dB), and (b) conductance, G (dB).  (a) (b) Figure 11. Illustration of simulated electrical responses using different patch depths (i.e., 4 ⁄ , 2 ⁄ , 3 4 ⁄ , and ). The patch length ( ) of 3 ⁄ and the patch metal of tungsten (W) are used for the simulation. (a) Admittance, Y11 (dB), and (b) conductance, G (dB).

Conclusions
In conclusion, the proposed suspended energy-preserving addendum frame lends itself to acoustic energy confinement applications, as it offered a significant improvement in the FoM of about 11-fold compared to the conventional resonator. Using FEA-based modeling techniques, we found a design that led to an improvement in the unloaded quality factor ( ) as high as 51,136, approximately 12.5-fold that of the conventional resonator, which offered a maximum value of about 4106. The best insertion loss value recorded was 0.7 dB for the resonator with the proposed technique, which is significantly better than the 5 dB insertion loss value obtained for the conventional resonator. Aside from this, the proposed method does not have a significant effect on the effective electromechanical coupling coefficient. Moreover, a new method for suppressing spurious modes is introduced, from which the FEA simulation results demonstrate a superior performance in mitigating the spurious modes. Conflicts of Interest: The authors declare no conflict of interest.

Conclusions
In conclusion, the proposed suspended energy-preserving addendum frame lends itself to acoustic energy confinement applications, as it offered a significant improvement in the FoM of about 11-fold compared to the conventional resonator. Using FEA-based modeling techniques, we found a design that led to an improvement in the unloaded quality factor (Q u ) as high as 51,136, approximately 12.5-fold that of the conventional resonator, which offered a maximum value of about 4106. The best insertion loss value recorded was 0.7 dB for the resonator with the proposed technique, which is significantly better than the 5 dB insertion loss value obtained for the conventional resonator. Aside from this, the proposed method does not have a significant effect on the effective electromechanical coupling coefficient. Moreover, a new method for suppressing spurious modes is introduced, from which the FEA simulation results demonstrate a superior performance in mitigating the spurious modes.