Resonators based on two-dimensional (2D) materials have exceptional properties for application as nanomechanical sensors, which allows them to operate at high frequencies with high sensitivity. However, their performance as nanomechanical sensors is currently limited by their low quality ( Q)-factor. Here, we make use of micro-electromechanical systems (MEMS) to apply pure in-plane mechanical strain, enhancing both their resonance frequency and Q-factor. In contrast to earlier work, the 2D material resonators are fabricated on the MEMS actuators without any wet processing steps using a dry-transfer method. A platinum clamp, which is deposited by electron beam-induced deposition, is shown to be effective in fixing the 2D membrane to the MEMS and preventing slippage. By in-plane straining the membranes in a purely mechanical fashion, we increase the tensile energy, thereby diluting dissipation. This way, we show how dissipation dilution can increase the Q-factor of 2D material resonators by 91%. The presented MEMS actuated dissipation dilution method does not only pave the way toward higher Q-factors in resonators based on 2D materials, but also provides a route toward studies of the intrinsic loss mechanisms of 2D materials in the monolayer limit.

Nanomechanical resonators made of two-dimensional (2D) materials are the subject of intensive research due to their remarkable properties. Their low mass, combined with their high Young’s modulus, leads to resonance frequencies that are typically a few tens of MHz.1 Yet, their extreme flexibility in the out-of-plane direction enhances the sensitivity to external stimuli and makes them promising for various applications, including mass,1–4 force,5–7 pressure,6,8–11 and temperature sensing.6,12,13

The performance of nanomechanical resonant sensors and clocks is generally limited by their dissipation per cycle ( Δ W). A low Δ W results in a high quality ( Q-)factor, which is the ratio of stored energy W to Δ W over a single oscillation cycle ( Q = 2 π W / Δ W). A low Δ W and, thus, a high Q-factor physically insulates the resonator from external noise sources, allowing long-term coherent oscillations while minimizing energy dissipation to the environment,14,15 thus enabling low phase-noise oscillators16 and high-performance noise-rejection filters.15,17,18

Increasing the tensile stress can be an effective strategy to realize high- Q resonators.19,20 The tension leads to an increase in the stored energy W without significantly affecting losses, thereby increasing the ratio W / Δ W and the Q-factor,21 an effect that is commonly known as dissipation dilution. Since this strategy was very successful in realizing high- Q resonators in SiN that can be grown with high intrinsic tensile stress,20 it was also considered as a method for increasing the quality factor of 2D materials that could not be grown with high intrinsic stress such that they might become serious contenders for high- Q sensors and high f Q resonant quantum devices. Out-of-plane electrostatic and thermal forces were used to achieve this quality factor tuning.22 However, attempts to increase the Q-factor through out-of-plane electrostatic gating of 2D material membranes typically result in a reduction of the Q-factor instead of an increase by dissipation dilution.4,5,23 This reduction is due to the voltage-dependent electronic Joule dissipation of the displacement current within the resonator.23,24 Furthermore, the out-of-plane electrostatic pulling force increases sidewall adhesion, which facilitates dissipation through coupling with the substrate.22,25,26 Thermal expansion-based tuning strategies to increase Q22 have the drawback of making it difficult to distinguish tension effects from other thermal effects on Q. For example, the change in the membrane temperature changes material parameters that can increase damping via, e.g., the thermoelastic dissipation mechanism.15,27

Considering these drawbacks, it would be ideal to apply pure in-plane mechanical stress to increase the Q of 2D resonators. To reach this goal, special micro-electromechanical system (MEMS) actuators were developed in earlier works.28–30 However, all these approaches have a step in which liquids, such as glues or etchants, are in contact with the 2D material, thereby affecting the device properties and Q. Moreover, the involved liquids and the complex fabrication process result in a low device yield, making it difficult to show the clamping efficiency of 2D materials on MEMS actuators. The field is, thus, in need of a simplified, dry fabrication method for MEMS actuators with integrated and clamped 2D materials.

In this paper, we achieve this goal and provide evidence of dissipation dilution in suspended 2D material resonators that are controllably tensioned in-plane using a micro-electromechanical systems (MEMS) actuator, as illustrated in Fig. 1. To enable these experiments, we introduce a dry-transfer method to precisely suspend 2D materials over MEMS gaps and rigidly clamp them with a layer of platinum using electron-beam-induced deposition (EBID). We actuate the membrane resonances optothermally and record the resulting motion using an interferometry setup, from which we extract the Q-factor ( Q) and the resonance frequency ( f 0). By applying strain with the MEMS actuator to the 2D material resonator in a mechanical and controllable fashion, we find an increase in the resonance frequency as well as the Q-factor, which are consistent with dissipation dilution. Our findings, thus, provide a new way for enhancing the Q-factor via dissipation dilution in 2D material resonators realized with a dry fabrication method.

FIG. 1.

Tunable 2D material resonator via MEMS-induced tension. (a) Side view of the final device illustrating the 2D material, EBID clamps, drive- and readout lasers, and the force F cd induced by the MEMS actuator. Dimensions are in micrometers. (b) SEM image under an angle of the device. (c) Top view of an empty MEMS actuator. Scale bar: 50  μm. Voltage over comb-drive fingers ( V cd). (d) Detailed view of the comb drive indicating the direction of the comb-drive force ( F cd).

FIG. 1.

Tunable 2D material resonator via MEMS-induced tension. (a) Side view of the final device illustrating the 2D material, EBID clamps, drive- and readout lasers, and the force F cd induced by the MEMS actuator. Dimensions are in micrometers. (b) SEM image under an angle of the device. (c) Top view of an empty MEMS actuator. Scale bar: 50  μm. Voltage over comb-drive fingers ( V cd). (d) Detailed view of the comb drive indicating the direction of the comb-drive force ( F cd).

Close modal

We use a MEMS actuator that is designed and fabricated in the commercially available XMB10 process from X-FAB.31,32 The resulting device (see Fig. 1) consists of a moving shuttle with 38 comb fingers that is held suspended by four serpentine flexures. The flexures are connected to fixed anchors that have aluminum bond pads for making electrical contact by wire bonding. The crystalline silicon shuttle has a thickness of 15  μm, a length of 520  μm, a comb finger length of 103  μm, and an asymmetric finger spacing of 2.0 and 4.0  μm. The membranes are suspended over a 6  μm trench between the fixed anchor and the moving shuttle. Inside the trench, at a depth of 5  μm below the shuttle surface, a suspended silicon beam acts as a mirror for interferometric readout of the membrane motion. In order to transfer the membranes onto the MEMS actuator, we first mechanically exfoliate 2D materials33 onto a 5 × 5 mm 2 PDMS sheet on a microscope slide. Next, we use a microscope to select membranes on the PDMS sheet with a minimum length of 20  μm such that they can cover the suspended trench as well as parts of the fixed anchor and the moving shuttle. Membranes are selected based on their flatness and uniformity. Once we find a suitable membrane, we use a dome-shaped PDMS stamp covered with a sacrificial polypropylene carbonate (PPC) film on a microscope slide to pick it up from the PDMS sheet.34 The utilization of a PDMS dome results in a smaller contact area with the MEMS actuator, approximately 350  μm in diameter, which allows the precise positioning of a membrane using a manual transfer system of HQ2D while minimizing contaminations. The microscope slides ensure that we can see the flake during transfer and position it accurately. We then bring the membrane on the PDMS dome in contact with the MEMS actuator and heat the stage to 110  °C (above the melting point of PPC). This causes the PPC film to melt and ensures the transfer of the membrane onto the designated area [see Fig. 2(b)]. After the transfer, we wire bond the MEMS actuator and connect all terminals to a common ground; this crucial step prevents any electrostatic force-induced movement during the rest of the fabrication process. Next, we remove the PPC film from the membrane and the MEMS actuator through annealing in a high vacuum oven at a pressure below 10 5 mbar for 3 h at a temperature of 300  °C (exceeding the decomposition temperature of PPC35). After annealing, we inspect the sample optically to confirm the PPC removal [see Fig. 2(c)]. Finally, we clamp the membrane with a layer of platinum using EBID.36 Other materials or composites for clamping might be considered as well for the ease of fabrication or enhancing the mechanical properties of the final device, e.g., by maximizing the stiffness of the clamps.37 We used platinum as it was readily available and is chemically inert, yet compatible with EBID. Detailed information on the fabrication procedure is available in the supplementary material S1.

FIG. 2.

Fabrication procedure: drawings, optical and SEM micrographs. (a) Device straining platform. Scale bar: 50  μm. (b) Device after the membrane transfer using a sacrificial PPC layer. (c) Device after annealing. (d) Illustration of the EBID process, an optical image of the device after clamping, and a false-colored SEM image highlighting the platinum clamps.

FIG. 2.

Fabrication procedure: drawings, optical and SEM micrographs. (a) Device straining platform. Scale bar: 50  μm. (b) Device after the membrane transfer using a sacrificial PPC layer. (c) Device after annealing. (d) Illustration of the EBID process, an optical image of the device after clamping, and a false-colored SEM image highlighting the platinum clamps.

Close modal

In total, we fabricated four devices (D1–D4) with different 2D materials using the method outlined above (see Table I). Figure 1(a) shows a schematic cross section of a device, including the dimensions. All membranes have a suspended length of 6  μm, width w, and thickness t, as determined with a white light interferometer (see Table I).

TABLE I.

Characteristics of fabricated devices D1–D4, including 2D material, clamped using EBID of platinum, width w, thickness t, and fitted parameters Q int / f plate 2 and Qext, where Qint is the intrinsic quality factor due to bending and elongation losses, fplate is the resonance frequency due to the bending rigidity, and Qext is damping from extrinsic dissipation sources [see Eq. (3)]. The membrane thickness is determined using AFM. The width is obtained by measuring an optical or SEM image.

DeviceMaterialClampedw (μm)t (nm) Q int / f plate 2 (Hz−2)QintQext
D1 Graphene No 7.2 9.8 ± 0.33 7.084 × 10−12 <71 182 
D2 Graphene No 16.7 22.47 ± 0.15 2.106 × 10−12 <127 >1000 
D3 MoS2 Yes 19.9 77.79± 0.33 1.362 × 10−12 <64 180 
D4 WS2 Yes 16.0 94.39 ± 0.23 2.515 × 10−12 <203 >1000 
DeviceMaterialClampedw (μm)t (nm) Q int / f plate 2 (Hz−2)QintQext
D1 Graphene No 7.2 9.8 ± 0.33 7.084 × 10−12 <71 182 
D2 Graphene No 16.7 22.47 ± 0.15 2.106 × 10−12 <127 >1000 
D3 MoS2 Yes 19.9 77.79± 0.33 1.362 × 10−12 <64 180 
D4 WS2 Yes 16.0 94.39 ± 0.23 2.515 × 10−12 <203 >1000 

To mechanically tension 2D materials with the MEMS actuators [see Figs. 1(b)1(d)], we apply a potential difference V cd between the comb fingers. Due to the asymmetric placement of the comb fingers [see Fig. 1(d)], a force F cd = 1 2 C x V cd 2 will act on the moving shuttle and will tension the membrane and the four serpentine flexures. Here, C x is the change in capacitance C between the comb fingers with respect to a change in position x of the moving shuttle.

Based on their geometry, the MEMS actuators used here can controllably strain membranes up to 11%, which is set by the maximum in-plane displacement of 0.67 μm (1/3 of the 2 μm actuation gap) over a 6  μm suspended length. As the maximum in-plane displacement is the limiting factor, pre-deformations and wrinkles in 2D materials22,26,38 can reduce the maximum achievable strain. Moreover, the bare MEMS actuator has a pull-in voltage of 13.5 ± 0.5  V (see supplementary material S2). As a result, if the stiffness of the 2D membrane stiffness is low, e.g., because it is wrinkled, the MEMS actuator will collapse at this voltage. However, for certain devices, the stiffness of the 2D membrane was high such that it increased the total device stiffness significantly with respect to that of the MEMS springs and, as a consequence, led to an increase in the pull-in voltage. We were able to apply voltages up to 60 V for device D4 without a pull-in, which demonstrates that the force and stiffness provided by the 2D material were substantial.

We measure the dynamics of the membranes using an optical interferometry setup [see Fig. 3(a)]. The devices are placed in a vacuum chamber with a pressure below 10 5 mbar and are actuated using a blue diode laser ( λ = 405 nm) that is power-modulated through a vector network analyzer (VNA). We use a red He–Ne laser ( λ = 632 nm) to measure the motion of the membrane as its reflected intensity highly depends on the distance between the membrane and the mirror [see Fig. 3(b)]. The intensity of the reflected red laser light is detected using a photodiode and further processed by the VNA. Figure 3(c) shows a typical response recorded by the VNA. We observe multiple peaks in the spectrum that we identify as resonance frequencies. To analyze these further, we fit the frequency response M of the fundamental resonance to the well-known harmonic oscillator model,15,39,40
(1)
where f 0 is the resonance frequency in Hz, Q is the quality factor, and A is the peak amplitude. Figure 3(c) shows this fit on the frequency response of device D4 at V cd = 2 V, from which we extract a resonance frequency of 8.22 MHz. We then repeat the measurement for different V cd [see Fig. 3(d)] to analyze the dependence of f 0 and Q on the applied force F cd. As Fig. 3(d) shows, we see an increase in f 0 with applied V cd due to the increase in tension in the membrane, which is in agreement with results reported in the literature.30,37,41,42 We observe a similar increase in the resonance frequency with applied V cd for the higher mode between 11 and 12 MHz, as depicted in supplementary material S3.
FIG. 3.

Device measurements. (a) Fabry–Perot interferometry setup for measuring the resonance frequencies. Vector Network Analyzer (VNA), V cd and stage control (control), Beam expander (BE), Beam splitter (BS), Quarter wave plate (QWP), Dichroic mirror (DM), White light source (WL), Photodetector (PD), Red ( λ = 632 nm) He-Ne Laser (RL), and Blue ( λ = 405 nm) Diode Laser (BL). (b) Side view of the device inside the high vacuum chamber. (c) Fitting the harmonic oscillator equation on the experimental data of the VNA for device D4 at V cd = 2 V to extract the resonance frequency ( f 0) and the quality factor (Q). (d) Experimental data of the VNA for device D4 for V cd = 0, 2, 4, 6, and 8 V (DC).

FIG. 3.

Device measurements. (a) Fabry–Perot interferometry setup for measuring the resonance frequencies. Vector Network Analyzer (VNA), V cd and stage control (control), Beam expander (BE), Beam splitter (BS), Quarter wave plate (QWP), Dichroic mirror (DM), White light source (WL), Photodetector (PD), Red ( λ = 632 nm) He-Ne Laser (RL), and Blue ( λ = 405 nm) Diode Laser (BL). (b) Side view of the device inside the high vacuum chamber. (c) Fitting the harmonic oscillator equation on the experimental data of the VNA for device D4 at V cd = 2 V to extract the resonance frequency ( f 0) and the quality factor (Q). (d) Experimental data of the VNA for device D4 for V cd = 0, 2, 4, 6, and 8 V (DC).

Close modal

To evaluate the effectiveness of a deposited layer of platinum using EBID in clamping 2D material resonators, we compare the response of an unclamped device (D1) to the response of a clamped one (D3). We vary the comb-drive force F cd for both devices by varying V cd. In a single sweep, we start at V cd = 0 V, increase the voltage to a maximum V max, return to 0 V, decrease the voltage to V max, and finally return to 0 V. For each subsequent sweep, we increase V max by 1 V up to a maximum of 5 V for device D1, to gradually increase the maximum force exerted on the membrane, and monitor the shift in the fundamental resonance. Figure 4 shows f 0 as a function of V cd extracted from one such experiment for devices D1 and D3. For the unclamped device [see Fig. 4(a)], we observe a significant decrease in f 0 at V cd = 0 V after each voltage sweep. At the last measurement point at 0 V, f 0 shows an irreversible reduction from 4.30 to 3.73 MHz, which corresponds to a decrease of approximately 13% with respect to the very first measurement on this membrane. The permanent reduction of f 0 might be attributed to irreversible slippage, unsticking, or ironing out the wrinkles that increase the effective length of the membrane, reduce its tension, and, thus, resonance frequency. This is in contrast to recent observations of a reversible sliding scenario where a closed f 0 vs V cd loop is expected.43 

FIG. 4.

Comparison of the resonance frequency ( f 0) vs the comb-drive voltage ( V cd). F cd indicates the pulling direction of the suspended shuttle by the comb-drive actuator. Scale bars: 10  μm. (a) Unclamped device D1, V max = 2 V (blue), 3 V (orange), 4 V (yellow), and 5V (purple), false-colored optical image (blue: 2D membrane) (b) Device D3 clamped with a platinum layer deposited by EBID, V max = 30 V (blue), false-colored SEM image (blue: 2D membrane, purple: platinum clamps).

FIG. 4.

Comparison of the resonance frequency ( f 0) vs the comb-drive voltage ( V cd). F cd indicates the pulling direction of the suspended shuttle by the comb-drive actuator. Scale bars: 10  μm. (a) Unclamped device D1, V max = 2 V (blue), 3 V (orange), 4 V (yellow), and 5V (purple), false-colored optical image (blue: 2D membrane) (b) Device D3 clamped with a platinum layer deposited by EBID, V max = 30 V (blue), false-colored SEM image (blue: 2D membrane, purple: platinum clamps).

Close modal

In contrast to the unclamped device D1, device D3 contains a clamped membrane, as depicted in the inset of Fig. 4(b). As Fig. 4(b) shows, we observe a notably different f 0 vs V cd response when compared to the unclamped device. Initially, f 0 measures 7.49 MHz, ascending to 13.67 MHz at 30 V, and decreasing back to 7.42 MHz at 0 V, a change of less than 1%. In subsequent sweeps, the resonance frequency at V cd = 0 V remains stable within a range of 0.02 MHz. These measurements show that the resonance frequency f 0 of the device D3 is much more stable than that of the unclamped device, even at substantially higher actuation voltages, resulting in a 36 times larger force ( F cd V cd 2) than applied in measurements on device D1. This difference between the clamped and unclamped devices indicates that the deposited layer of platinum using EBID is effective in preventing permanent tension reduction during MEMS actuation. Since the EBID clamps are separated by a few micrometers from the edge of the trench and, thus, cannot significantly affect unwrinkling and adhesion mechanisms inside the trench, we also conclude that the irreversible changes in f 0 in Fig. 4(a) are most likely due to slippage of a large part of the flake over the silicon surface.

Furthermore, it is worth noting that the maximum voltage of 30 V significantly surpasses the pull-in voltage V PI of the bare MEMS actuator, which was determined to be 13.5 ± 0.5 V (see supplementary material S2), providing evidence that the clamped membrane generates a substantial force on the shuttle, thus preventing its pull-in and collapse. Simultaneously, according to Newton’s third law, we conclude that the MEMS actuator effectively pulls and strains the clamped membranes. An equivalent observation for unclamped device D2 and clamped device D4 can be found in supplementary material S4.

By fitting resonance peaks, such as in Fig. 3(d), we extract both the quality factor Q and resonance frequency f 0. To study the effectiveness of the MEMS actuator in tuning the quality factor by dissipation dilution, we plot Q vs f 0 for both unclamped and clamped devices over the full actuation voltage range in Fig. 5. It is seen that the Q-factor of all devices increases with f 0, as expected in a dissipation dilution scenario. For instance, device D4 experiences a 30% increase in f 0 and a 91% increase in Q. Interestingly, despite the presence of slipping effects, this trend even seems to hold for the unclamped devices D1 and D2 shown in Figs. 5(a) and 5(b), respectively.

FIG. 5.

Dissipation dilution in unclamped devices D1–D2 and clamped devices D3–D4. Quality factor ( Q) vs the resonance frequency ( f 0). The unclamped devices exhibit a decline in f 0 and Q due to slippage, in contrast with the clamped devices, which show consistent values for f 0 and Q at V cd = 0 V. (a) Experimental data from unclamped device D1 from the V max = 4 V-cycle. The device starts at f 0 = 4.17 MHz and Q = 73.5 and stops at f 0 = 3.96 MHz and Q = 69.0 due to slippage. (b) Experimental data from unclamped device D2 from the V max = 12 V-cycle. The device starts at f 0 = 8.45 MHz and Q = 150 and stops at f 0 = 7.95 MHz and Q = 113 due to slippage. (c) Experimental data of clamped device D3 from the V max = 30 V-cycle. The device maintains a consistent f 0 and Q of 7.50 MHz and 63.8, respectively. (d) Experimental data of clamped device D4 from the V max = 60 V-cycle. The device maintains a consistent f 0 and Q of 9.70 MHz and 202, respectively.

FIG. 5.

Dissipation dilution in unclamped devices D1–D2 and clamped devices D3–D4. Quality factor ( Q) vs the resonance frequency ( f 0). The unclamped devices exhibit a decline in f 0 and Q due to slippage, in contrast with the clamped devices, which show consistent values for f 0 and Q at V cd = 0 V. (a) Experimental data from unclamped device D1 from the V max = 4 V-cycle. The device starts at f 0 = 4.17 MHz and Q = 73.5 and stops at f 0 = 3.96 MHz and Q = 69.0 due to slippage. (b) Experimental data from unclamped device D2 from the V max = 12 V-cycle. The device starts at f 0 = 8.45 MHz and Q = 150 and stops at f 0 = 7.95 MHz and Q = 113 due to slippage. (c) Experimental data of clamped device D3 from the V max = 30 V-cycle. The device maintains a consistent f 0 and Q of 7.50 MHz and 63.8, respectively. (d) Experimental data of clamped device D4 from the V max = 60 V-cycle. The device maintains a consistent f 0 and Q of 9.70 MHz and 202, respectively.

Close modal
We now compare the experimental relation between f 0 and Q with theory. From the literature,14 the quality factor Q D of a resonator in the presence of dissipation dilution is given by (see supplementary material S5)
(2)
where f 0 is the measured resonance frequency, f plate is the frequency in the presence of bending rigidity, and Q int is the intrinsic dissipation. In the presence of other external dissipation mechanisms with quality factor Q ext, the total Q-factor can be further reduced15,19 to a value Q,
(3)

We fit Eq. (3) to the experimental data in Fig. 5 with Q int / f plate 2 and Q ext as fit parameters (orange lines). We provide the fitting values in Table I. Note that the fitted lines are almost linear because they span only a small frequency range. For the fit parameter Q int / f plate 2, we find a similar order of magnitude for all devices. As we assume that f 0 equals f plate when the membranes are unstrained ( V cd = 0 V, see supplementary material S5), we estimate Q int by multiplying Q int / f plate 2 with f 0 ( 0 V ) 2. The large values of Q int (see Table I) support the underpinning hypothesis of Eq. 2 that the storage modulus of these materials is much larger than their loss modulus. For devices D1 and D3, we fitted an extrinsic damping contribution of Q ext 180. For devices D2 and D4, the data could be fitted well without assuming an external dissipation source Q ext, and therefore, we took 1 / Q ext = 0. The good agreement between the experimental data and the fitted curves using Eq. (3) provides evidence that dissipation dilution can account for the observations.

In Secs. IIV, we studied the effect of in-plane stress on the Q-factor and the resonance frequency of membranes made of multi-layered 2D materials. The results show that the quality factor can be enhanced by generating tension on the 2D material membranes, and the obtained data match reasonably well with a dissipation dilution model. We note that this observation does not rule out other mechanisms that may affect the quality factors. These include (i) mode coupling to other modes of the membrane;44–46 (ii) a change in dissipation rate Δ W with applied tension;26,47,48 (iii) changes in the membrane’s geometry with increasing tension, which could be attributed to suppressing wrinkles;38,49 (iv) the release of edge adhesion during membrane straining,50 which facilitates the transfer of energy from the resonating membrane to the substrate, leading to acoustic radiation losses;51 or (v) material fatigue behavior of both the 2D material and the platinum clamps under repeated mechanical cycling. Future work could focus on systematically studying these effects as a function of applied tension.

The main advantage of using a MEMS actuator to apply the tension is that it provides a pure in-plane force, in contrast to thermal, electrostatic gating, or gas pressure-based approaches, where the force is accompanied by other effects that might modify Q. The presented method also has some limitations. First of all, actuator device fabrication, design, and wire bonding are time-consuming. Moreover, the actuation range is limited to 1/3 of the actuator gap, although other actuator designs can provide a more extensive range. Finally, transferring membranes to a movable shuttle is a complicated task and can affect membrane suspending and adhesion to the MEMS actuator.

A key prospect of the presented work is the ability to apply MEMS tensioning to significantly raise the quality factor of 2D materials, with the aim to eventually scale it down to the monolayer limit and outperform record silicon nitride/carbide devices.52–54 Although currently, Qs of 2D materials seem to be limited by their low stress and by fabrication artifacts, such as wrinkles, in the monolayer limit, the intrinsic Q int of 2D materials might be much higher than that of multilayer materials since interlayer dissipation mechanisms and thermoelastic damping are largely eliminated.55 Moreover, for sensing applications, such as resonant mass sensors, the large aspect ratio and low mass of 2D materials can increase sensitivity. In this respect, the high ultimate tensile stress of 2D materials, such as graphene, which was demonstrated to exceed 10 GPa,56 can potentially outperform the tension limitations of even silicon carbides.52 

Moreover, it is of interest to evaluate the minimum measured Q in more detail. For device D4, we found a Q of 202. Disregarding the small effect of pre-stress, we expect this Q to be roughly equal to the intrinsic Q since the thickness of the membrane is relatively high and the intrinsic Q is primarily governed by the bending contribution. This intrinsic Q is over 32 times lower than that of silicon nitride, for which a value of 69 t nm 1 was found,57 which for a thickness equal to that of device D4 of 94.39 nm evaluates to Q int = 6513. A question for further study is if the low intrinsic Q in device D4 is due to the intrinsic crystal properties of WS 2 compared to SiN or if the value is limited by imperfections, such as wrinkles. It is clear that the MEMS implementation of dissipation dilution demonstrated in this work, which results in Q enhancement up to 91%, is only a first step toward Q-factor engineering of 2D materials and by no means is yet as efficient as that in SiN resonators, where Q / Q int factors of over 10 4 are reached.53 

To conclude, we show signatures of dissipation dilution in 2D material resonators. We use a MEMS actuator to strain the membrane uniaxially and, thus, tune the resonance frequency and the Q-factor. To induce dissipation dilution, we developed a device fabrication method using dry transfer of membranes on the MEMS actuator and a clamping technique using EBID of platinum to prevent edge-slippage. The MEMS platform can also be used to study slipping and sliding effects, where slipping was observed to reduce both the Q-factor and f 0, in line with the dissipation dilution mechanism that increases Q with the MEMS actuation force. By pulling on the membranes, we control the in-plane tension of the membranes, which resulted in an increase in resonance frequency f 0 of 30% and an accompanying increase of 91% in the Q factor. Our results enable a leap in developing higher- Q resonators based on 2D materials with potential applications in sensing, time-keeping, and information processing.

See the supplementary material for detailed fabrication instructions (S1), MEMS stiffness characterization (S2), frequency tuning of a higher mode (S3), f0 vs V measurements for devices D2 and D4 (S4), and derivation of the dissipation dilution model (S5).

We would like to thank the Plantenna Research Program funded by the 4TU Federation, the European Union’s Horizon 2020 Research and Innovation Program under ERC Consolidator Grant NCANTO (No. 101125458) for financial support, and Europractice for providing the design tools and multi-project wafer service.

The authors have no conflicts to disclose.

Michiel P. F. Wopereis: Data curation (equal); Formal analysis (equal); Investigation (equal); Methodology (equal); Writing – original draft (equal); Writing – review & editing (equal). Niels Bouman: Conceptualization (equal); Data curation (equal); Formal analysis (equal); Investigation (equal); Methodology (equal). Satadal Dutta: Conceptualization (equal); Funding acquisition (equal). Peter G. Steeneken: Conceptualization (equal); Data curation (equal); Funding acquisition (equal); Supervision (equal); Writing – original draft (equal); Writing – review & editing (equal). Farbod Alijani: Conceptualization (equal); Data curation (equal); Supervision (equal); Writing – original draft (equal); Writing – review & editing (equal). Gerard J. Verbiest: Conceptualization (equal); Data curation (equal); Funding acquisition (equal); Supervision (equal); Visualization (equal); Writing – original draft (equal); Writing – review & editing (equal).

The data that support the findings of this study are available from the corresponding author upon reasonable request.

1.
A.
Sakhaee-Pour
,
M. T.
Ahmadian
, and
A.
Vafai
, “
Applications of single-layered graphene sheets as mass sensors and atomistic dust detectors
,”
Solid State Commun.
145
,
168
(
2008
).
2.
B.
Lassagne
,
D.
Garcia-Sanchez
,
A.
Aguasca
, and
A.
Bachtold
, “
Ultrasensitive mass sensing with a nanotube electromechanical resonator
,”
Nano Lett.
8
,
3735
(
2008
).
3.
J.
Atalaya
,
J. M.
Kinaret
, and
A.
Isacsson
, “
Nanomechanical mass measurement using nonlinear response of a graphene membrane
,”
Europhys. Lett.
91
,
48001
(
2010
).
4.
V.
Singh
,
S.
Sengupta
,
H. S.
Solanki
,
R.
Dhall
,
A.
Allain
,
S.
Dhara
,
P.
Pant
, and
M. M.
Deshmukh
, “
Probing thermal expansion of graphene and modal dispersion at low-temperature using graphene nanoelectromechanical systems resonators
,”
Nanotechnology
21
,
165204
(
2010
).
5.
P.
Weber
,
J.
Guttinger
,
A.
Noury
,
J.
Vergara-Cruz
, and
A.
Bachtold
, “
Force sensitivity of multilayer graphene optomechanical devices
,”
Nat. Commun.
7
,
12496
(
2016
).
6.
M. C.
Lemme
,
S.
Wagner
,
K.
Lee
,
X.
Fan
,
G. J.
Verbiest
,
S.
Wittmann
,
S.
Lukas
,
R. J.
Dolleman
,
F.
Niklaus
,
H. S. J.
van der Zant
,
G. S.
Duesberg
, and
P. G.
Steeneken
, “
Nanoelectromechanical sensors based on suspended 2D materials
,”
Research
2020
,
8748602
.
7.
I. E.
Rosłoń
,
A.
Japaridze
,
P. G.
Steeneken
,
C.
Dekker
, and
F.
Alijani
, “
Probing nanomotion of single bacteria with graphene drums
,”
Nat. Nanotechnol.
17
,
637
(
2022
).
8.
J.
Romijn
,
S.
Vollebregt
,
R. J.
Dolleman
,
M.
Singh
,
H. S.
Van Der Zant
,
P. G.
Steeneken
, and
P. M.
Sarro
, “A miniaturized low power Pirani pressure sensor based on suspended graphene,” in 2018 IEEE 13th Annual International Conference on Nano/Micro Engineered and Molecular Systems (NEMS) (IEEE, 2018).
9.
R. J.
Dolleman
,
D.
Davidovikj
,
S. J.
Cartamil-Bueno
,
H. S. J.
Van Der Zant
, and
P. G.
Steeneken
, “
Graphene squeeze-film pressure sensors
,”
Nano Lett.
16
,
568
(
2016
).
10.
D.
Shahdeo
,
A.
Roberts
,
N.
Abbineni
, and
S.
Gandhi
, “Graphene based sensors,” in Analytical Applications of Graphene for Comprehensive Analytical Chemistry, Comprehensive Analytical Chemistry (
Elsevier
, 2020) pp. 175–199.
11.
J. H.
Choi
,
J.
Lee
,
M.
Byeon
,
T. E.
Hong
,
H.
Park
, and
C. Y.
Lee
, “
Graphene-based gas sensors with high sensitivity and minimal sensor-to-sensor variation
,”
ACS Appl. Nano Mater.
3
,
2257
(
2020
).
12.
R. J.
Dolleman
,
G. J.
Verbiest
,
Y. M.
Blanter
,
H. S. J.
Van Der Zant
, and
P. G.
Steeneken
, “
Nonequilibrium thermodynamics of acoustic phonons in suspended graphene
,”
Phys. Rev. Res.
2
,
012058
(
2020
).
13.
H.
Liu
, et al., “Enhanced photothermal response near the buckling bifurcation in 2d nanomechanical resonators,”
2D Mater.
11,
025028
(2023).
14.
L.
Sementilli
,
E.
Romero
, and
W. P.
Bowen
, “
Nanomechanical dissipation and strain engineering
,”
Adv. Funct. Mater.
32
,
2105247
(
2022
).
15.
J. M. L.
Miller
,
A.
Ansari
,
D. B.
Heinz
,
Y.
Chen
,
I. B.
Flader
,
D. D.
Shin
,
L. G.
Villanueva
, and
T. W.
Kenny
, “
Effective quality factor tuning mechanisms in micromechanical resonators
,”
Appl. Phys. Rev.
5
,
041307
(
2018
).
16.
D.
Leeson
, “
A simple model of feedback oscillator noise spectrum
,”
Proc. IEEE
54
,
329
(
1966
).
17.
T. B.
Gabrielson
, “
Fundamental noise limits for miniature acoustic and vibration sensors
,”
J. Vib. Acoust.
117
,
405
(
1995
).
18.
T.
Gabrielson
, “
Mechanical-thermal noise in micromachined acoustic and vibration sensors
,”
IEEE Trans. Electron Devices
40
,
903
(
1993
).
19.
S.
Schmid
,
Fundamentals of Nanomechanical Resonators
(
Springer
,
2016
).
20.
S. S.
Verbridge
,
J. M.
Parpia
,
R. B.
Reichenbach
,
L. M.
Bellan
, and
H. G.
Craighead
, “
High quality factor resonance at room temperature with nanostrings under high tensile stress
,”
J. Appl. Phys.
99
,
124304
(
2006
).
21.
S. A.
Fedorov
,
N. J.
Engelsen
,
A. H.
Ghadimi
,
M. J.
Bereyhi
,
R.
Schilling
,
D. J.
Wilson
, and
T. J.
Kippenberg
, “
Generalized dissipation dilution in strained mechanical resonators
,”
Phys. Rev. B
99
,
054107
(
2019
).
22.
D.
Davidovikj
,
M.
Poot
,
S. J.
Cartamil-Bueno
,
H. S. J.
Van Der Zant
, and
P. G.
Steeneken
, “
On-chip heaters for tension tuning of graphene nanodrums
,”
Nano Lett.
18
,
2852
(
2018
).
23.
N.
Morell
,
A.
Reserbat-Plantey
,
I.
Tsioutsios
,
K. G.
Schädler
,
F.
Dubin
,
F. H. L.
Koppens
, and
A.
Bachtold
, “
High quality factor mechanical resonators based on WSe 2 monolayers
,”
Nano Lett.
16
,
5102
(
2016
).
24.
X.
Song
,
M.
Oksanen
,
M. A.
Sillanpää
,
H. G.
Craighead
,
J. M.
Parpia
, and
P. J.
Hakonen
, “
Stamp transferred suspended graphene mechanical resonators for radio frequency electrical readout
,”
Nano Lett.
12
,
198
(
2012
).
25.
J.
Rieger
,
A.
Isacsson
,
M. J.
Seitner
,
J. P.
Kotthaus
, and
E. M.
Weig
, “
Energy losses of nanomechanical resonators induced by atomic force microscopy-controlled mechanical impedance mismatching
,”
Nat. Commun.
5
,
3345
(
2014
).
26.
P. G.
Steeneken
,
R. J.
Dolleman
,
D.
Davidovikj
,
F.
Alijani
, and
H. S. J.
van der Zant
, “
Dynamics of 2D material membranes
,”
2D Mater.
8
,
042001
(
2021
).
27.
R.
Lifshitz
and
M. L.
Roukes
, “
Thermoelastic damping in micro- and nanomechanical systems
,”
Phys. Rev. B
61
,
5600
(
2000
).
28.
G. J.
Verbiest
,
D.
Xu
,
M.
Goldsche
,
T.
Khodkov
,
S.
Barzanjeh
,
N.
Von Den Driesch
,
D.
Buca
, and
C.
Stampfer
, “
Tunable mechanical coupling between driven microelectromechanical resonators
,”
Appl. Phys. Lett.
109
,
143507
(
2016
).
29.
H. H.
Perez-Garza
,
E. W.
Kievit
,
G. F.
Schneider
, and
U.
Staufer
, “
Highly strained graphene samples of varying thickness and comparison of their behaviour
,”
Nanotechnology
25
,
465708
(
2014
).
30.
Y.
Xie
,
J.
Lee
,
Y.
Wang
, and
P. X.
Feng
, “
Straining and tuning atomic layer nanoelectromechanical resonators via comb–drive MEMS actuators
,”
Adv. Mater. Technol.
6
,
2000794
(
2021
).
31.
X.
Zou
,
S.
Ahmed
,
N.
Jaber
, and
H.
Fariborzi
, “A compact high-sensitivity temperature sensor using an encapsulated clamped-clamped MEMS beam resonator,” in 2021 21st International Conference on Solid-State Sensors, Actuators and Microsystems (Transducers) (IEEE, 2021).
32.
P.
Streit
,
R.
Forke
,
S.
Voigt
,
U.
Schwarz
,
R.
Ziegenhardt
,
S.
Weidlich
,
D.
Billep
,
M.
Gaitzsch
, and
H.
Kuhn
, “Vibration sensors with a high bandwidth and low SNR, enhanced with post processing gap reduction,” in 2022 23rd International Conference on Thermal, Mechanical and Multi-Physics Simulation and Experiments in Microelectronics and Microsystems (EuroSimE) (IEEE, 2022).
33.
K. S.
Novoselov
,
D.
Jiang
,
F.
Schedin
,
T. J.
Booth
,
V. V.
Khotkevich
,
S. V.
Morozov
, and
A. K.
Geim
, “
Two-dimensional atomic crystals
,”
Proc. Natl. Acad. Sci. U.S.A.
102
,
10451
(
2005
).
34.
K.
Kinoshita
,
R.
Moriya
,
M.
Onodera
,
Y.
Wakafuji
,
S.
Masubuchi
,
K.
Watanabe
,
T.
Taniguchi
, and
T.
Machida
, “
Dry release transfer of graphene and few-layer h-bn by utilizing thermoplasticity of polypropylene carbonate
,”
npj 2D Mater. Appl.
3
,
22
(
2019
).
35.
S.
Wang
,
Y.
Huang
,
B.
Liao
,
G.
Lin
,
G.
Cong
, and
L.
Chen
, “
Structure and properties of poly(propylene carbonate)
,”
Int. J. Polym. Anal. Charact.
3
,
131
(
1997
).
36.
M.
Lee
,
D.
Davidovikj
,
B.
Sajadi
,
M.
Siskins
,
F.
Alijani
,
H. S. J.
van der Zant
, and
P. G.
Steeneken
, “
Sealing graphene nanodrums
,”
Nano Lett.
19
,
5313
(
2019
).
37.
M.
Goldsche
,
G. J.
Verbiest
,
T.
Khodkov
,
J.
Sonntag
,
N. V. D.
Driesch
,
D.
Buca
, and
C.
Stampfer
, “
Fabrication of comb-drive actuators for straining nanostructured suspended graphene
,”
Nanotechnology
29
,
375301
(
2018
).
38.
R. J.
Nicholl
,
H. J.
Conley
,
N. V.
Lavrik
,
I.
Vlassiouk
,
Y. S.
Puzyrev
,
V. P.
Sreenivas
,
S. T.
Pantelides
, and
K. I.
Bolotin
, “
The effect of intrinsic crumpling on the mechanics of free-standing graphene
,”
Nat. Commun.
6
,
8789
(
2015
).
39.
R. M.
Schmidt
,
G.
Schitter
, and
A.
Rankers
,
The Design of High Performance Mechatronics: High-Tech Functionality by Multidisciplinary System Integration
(
IOS Press
,
2020
).
40.
R. S.
Figliola
and
D. E.
Beasley
,
Theory and Design for Mechanical Measurements
(
John Wiley & Sons
,
2014
).
41.
G. J.
Verbiest
,
M.
Goldsche
,
J.
Sonntag
,
T.
Khodkov
,
N.
von den Driesch
,
D.
Buca
, and
C.
Stampfer
, “
Tunable coupling of two mechanical resonators by a graphene membrane
,”
2D Mater.
8
,
035039
(
2021
).
42.
J.
Sonntag
,
M.
Goldsche
,
T.
Khodkov
,
G.
Verbiest
,
S.
Reichardt
,
N. V.
Den Driesch
,
D.
Buca
, and
C.
Stampfer
, “Engineering tunable strain fields in suspended graphene by microelectromechanical systems,” in 2019 20th International Conference on Solid-State Sensors, Actuators and Microsystems & Eurosensors XXXIII (TRANSDUCERS & EUROSENSORS XXXIII) (IEEE, 2019).
43.
Y.
Ying
,
Z.-Z.
Zhang
,
J.
Moser
,
Z.-J.
Su
,
X.-X.
Song
, and
G.-P.
Guo
, “
Sliding nanomechanical resonators
,”
Nat. Commun.
13
,
6392
(
2022
).
44.
X.
Wang
,
D.
Zhu
,
X.
Yang
,
L.
Yuan
,
H.
Li
,
J.
Wang
,
M.
Chen
,
G.
Deng
,
W.
Liang
,
Q.
Li
,
S.
Fan
,
G.
Guo
, and
K.
Jiang
, “
Stressed carbon nanotube devices for high tunability, high quality factor, single mode GHz resonators
,”
Nano Res.
11
,
5812
(
2018
).
45.
A.
Keşkekler
,
O.
Shoshani
,
M.
Lee
,
H. S. J.
Van Der Zant
,
P. G.
Steeneken
, and
F.
Alijani
, “
Tuning nonlinear damping in graphene nanoresonators by parametric–direct internal resonance
,”
Nat. Commun.
12
,
1099
(
2021
).
46.
A.
Keşkekler
,
V.
Bos
,
A. M.
Aragón
,
P. G.
Steeneken
, and
F.
Alijani
, “
Multimode nonlinear dynamics of graphene resonators
,”
Phys. Rev. Appl.
20
,
064020
(
2023
).
47.
C.
Zener
, “
Internal friction in solids. I. Theory of internal friction in reeds
,”
Phys. Rev.
52
,
230
(
1937
).
48.
S.
Schmid
,
K. D.
Jensen
,
K. H.
Nielsen
, and
A.
Boisen
, “
Damping mechanisms in high-q micro and nanomechanical string resonators
,”
Phys. Rev. B
84
,
165307
(
2011
).
49.
R. J.
Nicholl
,
N. V.
Lavrik
,
I.
Vlassiouk
,
B. R.
Srijanto
, and
K. I.
Bolotin
, “
Hidden area and mechanical nonlinearities in freestanding graphene
,”
Phys. Rev. Lett.
118
,
266101
(
2017
).
50.
R. J.
Dolleman
,
Y. M.
Blanter
,
H. S. J.
Van Der Zant
,
P. G.
Steeneken
, and
G. J.
Verbiest
, “
Phonon scattering at kinks in suspended graphene
,”
Phys. Rev. B
101
,
115411
(
2020
).
51.
P.
Steeneken
,
J.
Ruigrok
,
S.
Kang
,
J.
Van Beek
,
J.
Bontemps
, and
J.-J.
Koning
, “Parameter extraction and support-loss in MEMS resonators,” arXiv:1304.7953 (2013).
52.
M.
Xu
,
D.
Shin
,
P. M.
Sberna
,
R.
Van Der Kolk
,
A.
Cupertino
,
M. A.
Bessa
, and
R. A.
Norte
, “
High–strength amorphous silicon carbide for nanomechanics
,”
Adv. Mater.
36
,
2306513
(
2023
).
53.
D.
Shin
,
A.
Cupertino
,
M. H. J.
de Jong
,
P. G.
Steeneken
,
M. A.
Bessa
, and
R. A.
Norte
, “
Spiderweb nanomechanical resonators via Bayesian optimization: Inspired by nature and guided by machine learning
,”
Adv. Mater.
34
,
2106248
(
2022
).
54.
A.
Cupertino
,
D.
Shin
,
L.
Guo
,
P. G.
Steeneken
,
M.
Bessa
, and
R.
Norte
,
Centimeter-scale Nanomechanical Resonators with low Dissipation
(
2023
).
55.
N.
Lindahl
,
D.
Midtvedt
,
J.
Svensson
,
O. A.
Nerushev
,
N.
Lindvall
,
A.
Isacsson
, and
E. E. B.
Campbell
, “
Determination of the bending rigidity of graphene via electrostatic actuation of buckled membranes
,”
Nano Lett.
12
,
3526
(
2012
).
56.
M.
Goldsche
,
J.
Sonntag
,
T.
Khodkov
,
G. J.
Verbiest
,
S.
Reichardt
,
C.
Neumann
,
T.
Ouaj
,
N.
Von Den Driesch
,
D.
Buca
, and
C.
Stampfer
, “
Tailoring mechanically tunable strain fields in graphene
,”
Nano Lett.
18
,
1707
(
2018
).
57.
L.
Villanueva
and
S.
Schmid
, “
Evidence of surface loss as ubiquitous limiting damping mechanism in SiN micro- and nanomechanical resonators
,”
Phys. Rev. Lett.
113
,
227201
(
2014
).