High quality factor graphene-based 2D heterostructure mechanical resonator

High quality factor graphene-based 2D heterostructure mechanical resonator

M. Will JARA-FIT and 2nd Institute of Physics, RWTH Aachen University, 52074 Aachen, Germany    M. Hamer School of Physics and Astronomy and Manchester Centre for Mesoscience and Nanotechnology, University of Manchester, Oxford Road, Manchester M13 9PL, United Kingdom    M. Müller JARA-FIT and 2nd Institute of Physics, RWTH Aachen University, 52074 Aachen, Germany    A. Noury    P. Weber    A. Bachtold ICFO-Institut de Ciencies Fotoniques, The Barcelona Institute of Science and Technology, 08860 Castelldefels, Barcelona, Spain    R. V. Gorbachev School of Physics and Astronomy and Manchester Centre for Mesoscience and Nanotechnology, University of Manchester, Oxford Road, Manchester M13 9PL, United Kingdom    C. Stampfer stampfer@physik.rwth-aachen.de JARA-FIT and 2nd Institute of Physics, RWTH Aachen University, 52074 Aachen, Germany    J. Güttinger JARA-FIT and 2nd Institute of Physics, RWTH Aachen University, 52074 Aachen, Germany

Ultralight mechanical resonators based on low-dimensional materials are well suited as exceptional transducers of minuscule forces or mass changes. However, the low dimensionality also provides a challenge to minimize resistive losses and heating. Here, we report on a novel approach that aims to combine different 2D materials to tackle this challenge. We fabricated a heterostructure mechanical resonator consisting of few layers of niobium diselenide (NbSe) encapsulated by two graphene sheets. The hybrid membrane shows high quality factors up to 245’000 at low temperatures, comparable to the best few-layer graphene mechanical resonators. In contrast to few-layer graphene resonators, the device shows reduced electrical losses attributed to the lower resistivity of the NbSe layer. The peculiar low temperature dependence of the intrinsic quality factor points to dissipation over two-level systems which in turn relax over the electronic system. Our high sensitivity readout is enabled by coupling the membrane to a superconducting cavity which allows for the integration of the hybrid mechanical resonator as a sensitive and low loss transducer in future quantum circuits.

mechanical resonator, graphene, NbSe2, 2D heterostructures, NEMS, cavity readout

Peter Grünberg Institute (PGI-9), Forschungszentrum Jülich, 52425 Jülich, Germany \abbreviations2D,Nb,NbSe,PMMA, TLS

Keywords: mechanical resonator, graphene, NbSe2, 2D heterostructures, NEMS, cavity readout.

Van-der-Waals heterostructures based on complementary two-dimensional (2D) materials are a topic of intense research 1, 2 as they can be used to fabricate tailored electrical and optical devices with superior properties 3, 4, 5, 6. While mechanical devices based on individual 2D materials have shown tunable mechanical frequency and high quality factors 7, 8, 9, 10, 11, 12, 13, 14, the suitability of 2D heterostructures for mechanical resonators has not been explored so far. Another interesting question is how mechanical vibrations will interact with exotic states in encapsulated 2D materials such as 2D superconductors 6, 15 or 2D magnets 16, 17.

In order to use heterostrucutres for mechanical applications, it is crucial that the mechanical quality factor is not significantly degraded by interlayer friction forces between the materials. So far, a detailed understanding of the energy dissipation in 2D mechanical resonators down to low temperatures is missing, despite various calculations 18, 19, 20 and experiments 21, 22, 23, 24. On the one hand calculations on nanometer sized few-layer graphene resonators suggest the importance of interlayer friction forces 19, on the other hand the highest quality () factors were reported in multilayer graphene resonators 10, 14. Experimental limitations in studying the intrinsic mechanical dissipation are imposed by parasitic edge modes in doubly clamped devices 22, 24, the dependence of spectral measurements on amplitude fluctuations and frequency fluctuations 25, 14 and the challenge of a non-invasive readout technique down to milikelvin temperatures 14.

In this work we demonstrate a high quality van-der-Waals heterostructure mechanical resonator based on few-layer niobium diselenide (NbSe) encapsulated by graphene. The membrane motion is sensitively probed by capacitive coupling to a superconducting microwave cavity. This readout technique is ideal to study the energy decay as it allows for time resolved ring-down measurements which are independent of frequency fluctuations 14, 26. We demonstrate a high mechanical quality factor over 240’000 which compares favorably to the best multilayer graphene resonators 7, 9, 11, 14. Furthermore, by pulling the membrane electrostatically, the mechanical quality factor changes less compared with pure graphene resonators which can be explained by a reduced electrical resistance. By analyzing the low temperature dependence of the intrinsic energy decay we find an anomalous temperature dependence in the highest Q resonators of both types. This finding can be well understood by modeling the dissipation with two-level fluctuators that mediate the dissipation of mechanical energy to the electronic bath 27, 28, 29.

Figure 1: (a) Microscope image of the graphene/NbSe2/graphene assembly on the polymethyl methacrylate (PMMA) membrane before transfer to the microwave circuit. The single-layer nature of the graphene sheets is confirmed by Raman spectroscopy (see SI Section 2). The NbSe flake consists of - layers according to optical contrast measurements. (b) Optical microscope image of the hybrid membrane after transfer onto the niobium (Nb) electrodes (light gray). The membrane is fixed to the electrodes by crosslinked PMMA stripes and a ring shaped structure around the suspended device. (c) Schematic cross-section of the mechanical resonator along the black arrows in (b). The gate electrode allows to apply a potential between the membrane and the cavity electrode, displacing the membrane. The initial distance between the cavity electrode and the membrane is . (d) Optical contrast image of the entire chip, showing the superconducting cavity including a weak and a strong coupling port and the two electrodes contacting the membrane. Inset: Schematic of the measurement setup. A microwave source with frequency is connected to the weak port and the cavity field is detected over the strong port at the cavity resonance frequency . A static membrane-cavity voltage V and a AC mechanical drive voltage at frequency is applied at the membrane. (e) Exemplary spectral measurement of the mechanical amplitude as a function of with the extracted mechanical resonance frequency . (f) Resonant frequency of the hybrid membrane as a function of V. For the measurement we drive with a constant force (see SI Section 2). The inset shows that we pump the cavity with a red detuned pump tone such that .

The optical image of the device in Figure 1 shows the graphene-NbSe hybrid membrane, which is capacitively coupled to a superconducting niobium (Nb) cavity. The heterostructure is assembled in an argon atmosphere to avoid contamination and degradation of the NbSe layer 6. The membrane consists of few layers of NbSe that are encapsulated by single-layer graphene flakes on each side (see Figure 1a). Alternatively, it is also possible to use hBN crystals or other stable 2D materials for the encapsulation. The main reason to use graphene is the high mechanical strength of the material and the minimal added effective mass by monolayer encapsulation. Another beneficial aspect of graphene is the ease to find single-layer graphene. We expect the NbSe layer to induce superconductivity in the graphene due to the proximity effect. However, as the right contact showed an infinite resistance, we were not able to observe superconductivity in transport measurements. As previously demonstrated with few-layer graphene membranes 9, the encapsulated membrane is transferred on top of a pre-patterned Nb cavity and clamped by cross-linking the polymethyl methacrylate (PMMA) membrane that was used to transfer the 2D membrane (Figure 1b,c). The membrane is contacted over the left electrode which allows to apply a static gate voltage and a resonant driving voltage with respect to the grounded cavity electrode (Figure 1c,d). The as-fabricated device has an initial separation of , which is weakly reduced by for the maximum applied voltage of in this work. This static displacement is inferred from the shift of the cavity frequency as a function of due to the change of the cavity capacitance 26. Tunable mechanical resonance frequencies, low mass and high mechanical quality factors are key assets of mechanical resonators based on graphene. We show that this quantities can be maintained also in graphene-based heterostructure membranes. Figure 1e shows the spectral response of the hybrid resonator to a driving tone with a resonance frequency  MHz. This frequency is comparable to drums with graphene membranes of similar radius to this . The mechanical vibration is probed by the capacitive coupling between the electromagnetic pump field in the cavity and the mechanical motion 9. The pump field is injected into the cavity over a weakly coupled port and the scattered field is detected over a strongly coupled port and amplified by a low noise amplifier at the stage of the cryostat (Figure 1d). For all measurements we use a red-detuned pump tone such that the anti-Stokes scattered field becomes resonant with the cavity frequency (see inset Figure 1f). The coupling is characterized by the so called single photon-phonon coupling constant , which we estimate in SI Section 2. We use sufficiently low pump fields such that we can neglect optomechanical backaction. In order to evaluate the mass of the membrane we measure the dependence of (see Figure 1f). The decrease of with increasing is a hallmark of an electro-mechanical resonator under tension and can be attributed to a softening of the mechanical potential by the electric field. For the fundamental mode of a circular resonator the gate voltage dependent mechanical resonance frequency can be modeled by 9


Here is the maximum resonance frequency at the mechanical charge neutrality point (), is the electric constant, the effective mass of the fundamental mode, the two dimensional mass density of the membrane and the radius of the gate electrode. From a fit to the data we obtain and . Given the radius of the membrane, this effective mass corresponds to the mass of roughly 20 layers of graphene with . Considering that one layer of NbSe () has the mass of 5.2 graphene layers, the mass is in agreement with a heterostructure composed of three layers of NbSe encapsulated by two layers of graphene and some extra mass that might be attributed to polymer residues from the fabrication. From we can estimate the initial tension of the circular membrane 9


with the combined two-dimensional Young’s modulus. Here = 2 (3) is the number of graphene (NbSe) layers and is the elastic stiffness of graphene (NbSe30, 31.

Figure 2: (a) Energy decay measurements at various . The resonator drive is stopped at which leads to a subsequent decay of the mechanical amplitude. The plotted traces are the results of 10000 averaged energy decay measurements. The highest quality factors are obtained at low gate voltage with values up to at . (b) dependence of the inverse quality factor of the membrane (blue). Electronic dissipation is induced over capacitive displacement currents to compensate for the mechanical motion (schematic inset). The fitting parameter for the electronic dissipation is an effective resistance with (black line, see text). Every point is the average of three measurements.

In order to probe the mechanical energy dissipation independent of frequency noise and mechanical nonlinearities, we perform ring-down measurements (see Figure 2). For this measurement, the mechanical resonator is driven to a constant amplitude until the drive is stopped at time and the decay of the amplitude is recorded as a function of time (Figure 2a). From the exponential energy decay the decay rate is extracted at V. The corresponding quality factor is among the highest measured so far in mechanical resonators based on 2D materials 11, 14.

Thanks to reduced electrical loss in our device, a high mechanical Q factor of is still maintained at increased gate voltage in contrast to few-layer graphene mechanical resonators 11, 26. In Figure 2b the inverse mechanical quality factor obtained from energy decay measurements is plotted as a function of gate voltage . The decrease of the quality factor can be attributed to electrical loss induced by capacitive displacement currents due to the motion of the membrane 23, 26. By increasing more charges have to flow to compensate the change of the capacitance induced by the mechanical motion. The resulting loss of mechanical energy can be described by


with an effective resistance of the membrane and the derivative of the capacitance with respect to the deflection , which is obtained over the measured shift of the cavity resonance frequency as a function of . The data agrees well with with with an intrinsic dissipation and (black line in Figure 2b). This effective resistance is comparable to the sheet resistance of a three layer NbSe flake in the normal conducting state 6.
Despite the apparent absence of superconductivity in the membrane the Joule dissipation is reduced by at least a factor of four in the heterostructure device compared with graphene-based mechanical resonators (see Table 1).

Membrane [layers]
30 graphene
6 graphene
2 graphene, 3 NbSe
Table 1: Summary of gate dependent dissipation in different resonators (from Refs. 26, 14) in comparison with our heterostructure membrane regarding effective resistance and internal quality factor.
Figure 3: (a) Inverse quality factor as a function of the cryostat temperature . The color shading highlights three different temperature regimes (yellow, orange, red). The temperature of the mixing chamber (mc) has been measured with a calibrated Magnicon noise thermometer. (b) Relative change of the resonance frequency as a function of cryostat temperature with respect to the resonance frequency at . The inset depicts a sketch of a TLS with asymmetry and tunneling coupling . In panels (a) and (b) the different black lines correspond to the processes listed in Table 2. The blue line is the sum of the contributions from processes I and II. Note that in panel (a) the grey dot-dashed line with is used to fit the data in the yellow regime.

By analyzing the temperature dependence of and in Figure 3 we find evidence that the quality factor is limited by the coupling to two-level systems (TLS). Figure 3a shows the inverse quality factor as a function of cryostat temperature. The gate voltage was set to in order to minimize Joule heating losses and probe the intrinsic dissipation. Three different temperature dependencies are identified in the measurements indicated by the yellow, orange and red shading in Figure 3a. We will focus first on the anomalous crossover at , which cannot be explained by adding two independent dissipation channels with and to . This is because the stronger temperature dependence would dominate at higher temperatures. It is important to note that this behavior is not related to the heterostructure type of our membrane as a very similar dependence has been observed in a few-layer graphene resonator (see SI Section 3). Looking at previous low temperature measurements in graphene or nanotube-based resonators 32, 33, 34, 24 a proportionality of to or has been observed individually. However, an increase of the temperature dependence at lower temperature has not been reported before. We attribute this difference to the higher precision and reduced heating involved in the detection with a superconducting cavity combined with ring-down measurements. Next we will show that these features can be understood in the framework of two-level fluctuators that are relaxed by the interaction with conduction electrons.

We now recall the model of a two-level fluctuator or tunneling system, which has been proposed to describe low temperature effects in a wide range of isolating and metallic glasses and crystalline resonators 35, 36, 37, 27, 28, 38, 39, 40, 41, 42, 43, 34. In its simplest form, the tunneling system is defined by an energy asymmetry and a matrix element  28 for tunneling between the two states (see inset in Figure 3b). For the description of the physics it is often helpful to characterize two-level systems over its energy splitting and relaxation rate . Since the resonator interacts with an ensemble of two-level systems, a specific distribution for both and is assumed depending on the exact modeling 41. The mechanical motion interacts with the TLS over a periodic modulation of the potential. Thereby the TLS can absorb mechanical energy which is subsequently emitted in the phonon or electron bath providing a mechanical energy decay channel. The absorption process of the TLS can be either resonant (I) or over relaxation absorption (II).

(I) Resonant absorption is dominant at low temperatures where . The standard TLS model predicts and a resonance frequency shift for this process 28. The effects of resonant absorption can be saturated once the mechanical pumping rate gets stronger than the relaxation rate.

(II) Absorption relaxation is dominant for . In this process the mechanical strain fluctuations are affecting the TLS energy mainly over changes in the asymmetry energy  28. The resulting modulation of influences the interaction of the TLS with thermal phonons and electrons. This process gives rise to the peculiar kink in the temperature dependence at , where the fastest relaxation rate equals the mechanical frequency at . In the following we will discuss this process for both temperatures smaller (IIa) and larger (IIb) than . (IIa) For where there is a strong temperature dependence of the dissipation with for phonon relaxation and a linear dependence for relaxation over conduction electrons 27. Electronic relaxation leads to a greatly reduced relaxation time which lowers from a few Kelvin to the sub Kelvin regime for MHz resonators 28. The associated contribution to the frequency shift is negligible, so that the dependence of is given by resonant absorption processes involving TLS with , leading to . (IIb) At where the relaxation is faster than the mechanical modulation and the dissipation rate is proportional to . The resulting damping is . The resonance frequency shift due to IIb is negative with . The combined effect of I and II predicts a reduced but still positive frequency dependence. The positive dependence is not visible in the measurement and might be counteracted by a negative temperature dependence e.g. due to the negative thermal expansion of graphene (see SI Figure S5). Note that above a few Kelvin the temperature dependence is expected to become negative due to the onset of absorption relaxation by phonons 28.

Quantity Process Standard TLS model
I (saturable)
Table 2: Predictions for and based on the standard TLS theory in metallic systems 28. For both quantities the contribution from resonant (I) and relaxation absorption (II) add up. The crossover temperature separates the different contribution of relaxation absorption for (IIa) and for (IIb). depends on the frequency offset where and is related to the crossing temperature . The parameter is related to the TLS density and the coupling strength of deformations to the TLS energy splitting. The parameter is proportional to the conduction electron density and their coupling to TLS.

A quantitative comparison between theory and experiment is generally difficult and will serve here only to get a qualitative understanding of the parameters. We compare our measurements with the standard TLS theory 28, 41. We assume that electronic relaxation dominates over phonon relaxation in the measured temperature regime due to the low (fast relaxation time) and below . From the frequency dependence we extract K from the kink position and from the slope below (see formulas in Table 2). Note that the frequency offset has been adjusted such that the temperature corresponding to equals . From the saturation of for we extract . The difference between and might arise from changes in the density of the contributing TLS at the different temperatures. 40. We extract a TLS electron coupling constant of (dotted line) from the fit of the linear temperature dependence of contribution IIa. The temperature dependence in a high Q graphene device yields comparable values with and (see SI Section 3). These values are also in the same range that has been measured previously in metallic glasses and polycrystalline metals with and  41, 44. At the lowest temperatures the dissipation has a negative temperature dependence as predicted by theory (see also SI Figure S3). However, the measurement shows a reduced temperature dependence compared to the model (I), which might be explained by insufficient thermalization of the membrane to the mixing chamber plate at the lowest temperatures (mK, see Figure 3b). Please note that the fit merely emphasizes the change of sign in temperature dependence and is not meant to be universally valid. In the low temperature dependence of the graphene device we don’t observe any signature from process I. This might be explained with a complete saturation of the resonant TLS at the lowest drive amplitudes studied. Furthermore, in contrast to Ref. 45 we do not observe a saturation of the damping with increasing drive power but rather an increase in damping attributed to additional decay channels 14. A deviation from the theoretical prediction of the standard model is not unexpected as it has been observed before in metallic glasses 27. A nonsaturable attenuation might arise from more asymmetric TLS or a modified scattering cross-section 44. Additionally, the effect of individual TLS with and deviating from the assumed contribution is enhanced due to the tiny volume of the membrane 46.

The microscopic mechanism of TLS damping in 2D membranes is unclear so far. A previous calculation on dissipation in graphene and carbon nanotube resonators investigated electrically active TLS in the surrounding of the vibrating structure and concluded that the resulting dissipation will be negligible 18. However, residues from fabrication or imperfections in the crystal might still contain mechanically activated TLS 34.

TLS damping is not a universal limitation in 2D material based mechanical resonators. In resonators with lower factors we observe usually a saturation or reduced temperature dependence at low temperatures. This might be related to clamping induced loss 18 or loss over parasitic modes 22, 24 which masks the effect of the TLS.

We have not found any indication for superconductivity in the measured heterostructure device. First, the electric resistance extracted from the dissipation as a function of gate voltage is compatible with normal conducting NbSe. Second, we would expect the superconductivity to suppress the coupling of the TLS to electrons below  28, 42, which would lead to a higher kink temperature and increased temperature dependence compared to the measured data.

In conclusion, we presented a mechanical resonator based on van-der-Waals heterostructures of two-dimensional materials. The high mechanical quality, frequency tunability and low electrical resistivity measured in a graphene/NbSe/graphene resonator exemplifies the opportunities of 2D-heterostructures for mechanical resonators. The enhanced and partially linear temperature dependence of the dissipation is understood with mechanically mediated relaxation of two-level systems by conduction electrons. It might be possible to further improve the mechanical quality by using electrically isolating boron nitride for the encapsulation or reducing polymer residues in the fabrication. In future experiments, the combination of superconductivity and low mass in a 2D membrane opens the potential for ultra low loss and high sensitivity mechanical resonators. This is because a superconducting mechanical resonator suppresses electrical losses induced by mechanical motion and allows for higher pumping fields in a cavity based readout. Both are important for the exploration of quantum motion and improvements of force sensitivity in atomically thin mechanical resonators 9, 10, 11, 26.


We thank Uwe Wichmann for help on the electronics and for fruitful discussions. Support by the Helmholtz Nano Facility (HNF) 47 at the Forschungszentrum Jülich, the Excellence Initiative (RWTH Start-Up grant), the ERC-carbonNEMS and the ERC-GQEMS (GA-Nr. 280140) are gratefully acknowledged. The work is further supported by the ERC advanced grant 692876, the Foundation Cellex, Severo Ochoa (SEV-2015-0522), the grant FIS2015-69831-P of MINECO, and the Fondo Europeo de Desarrollo Regional (FEDER).

  • Supplement: Membrane characterization by Raman, calculation of resonator amplitude, data on temperature dependent dissipation in additional devices, effect of thermal expansion on resonance frequency and dependence of dissipation on mechanical drive voltage.


  • Novoselov \latinet al. 2016 Novoselov, K. S.; Mishchenko, A.; Carvalho, A.; Castro Neto, A. H. Science 2016, 353, 461–472.
  • Liu \latinet al. 2016 Liu, Y.; Weiss, N. O.; Duan, X.; Cheng, H.-C.; Huang, Y.; Duan, X. Nat. Rev. Mats. 2016, 1, 16042.
  • Dean \latinet al. 2010 Dean, C.; Young, A.; Meric, I.; Lee, C.; Wang, L.; Sorgenfrei, S.; Watanabe, K.; Taniguchi, T.; Kim, P.; Shepard, K.; J., H. Nat. Nanotech. 2010, 5, 722–726.
  • Britnell \latinet al. 2013 Britnell, L.; Ribeiro, R. M.; Eckmann, A.; Jalil, R.; Belle, B. D.; Mishchenko, A.; Kim, Y.-J.; Gorbachev, R. V.; Georgiou, T.; Morozov, S. V.; Grigorenko, A. N.; Geim, A. K.; Casiraghi, C.; Neto, A. H. C.; Novoselov, K. S. Science 2013, 340, 1311–1314.
  • Koppens \latinet al. 2014 Koppens, F.; Mueller, T.; Avouris, P.; Ferrari, A.; Vitiello, M.; Polini, M. Nat. Nanotech. 2014, 9, 780–793.
  • Cao \latinet al. 2015 Cao, Y. \latinet al. Nano Lett. 2015, 15, 4914–4921.
  • Eichler \latinet al. 2011 Eichler, A.; Moser, J.; Chaste, J.; Zdrojek, M.; Wilson-Rae, I.; Bachtold, A. Nat. Nanotech. 2011, 6, 339–342.
  • Lee \latinet al. 2013 Lee, J.; Wang, Z.; He, K.; Shan, J.; Feng, P. X.-L. ACS Nano 2013, 7, 6086–6091.
  • Weber \latinet al. 2014 Weber, P.; Güttinger, J.; Tsioutsios, I.; Chang, D. E.; Bachtold, A. Nano Lett. 2014, 14, 2854–2860.
  • Song \latinet al. 2014 Song, X.; Oksanen, M.; Li, J.; Hakonen, P. J.; Sillanpää, M. A. Phys. Rev. Lett. 2014, 113, 027404.
  • Singh \latinet al. 2014 Singh, V.; Bosman, S. J.; Schneider, B. H.; Blanter, Y. M.; Castellanos-Gomez, A.; Steele, G. A. Nat. Nanotech. 2014, 9, 1–5.
  • Castellanos-Gomez \latinet al. 2015 Castellanos-Gomez, A.; Singh, V.; van der Zant, H. S. J.; Steele, G. A. Ann. Phys. 2015, 527, 27–44.
  • Morell \latinet al. 2016 Morell, N.; Reserbat-Plantey, A.; Tsioutsios, I.; Schädler, K. G.; Dubin, F.; Koppens, F. H. L.; Bachtold, A. Nano Lett. 2016, 16, 5102–5108.
  • Güttinger \latinet al. 2017 Güttinger, J.; Noury, A.; Weber, P.; Eriksson, A. M.; Lagoin, C.; Moser, J.; Eichler, C.; Wallraff, A.; Isacsson, A.; Bachtold, A. Nat. Nanotech. 2017, 12, 631–636.
  • Xi \latinet al. 2015 Xi, X.; Zhao, L.; Wang, Z.; Berger, H.; Forro, L.; Shan, J.; Mak, K. F. Nat. Nanotech. 2015, 10, 765–769.
  • Lee \latinet al. 2016 Lee, J.-U.; Lee, S.; Ryoo, J. H.; Kang, S.; Kim, T. Y.; Kim, P.; Park, C.-H.; Park, J.-G.; Cheong, H. Nano Lett. 2016, 16, 7433–7438.
  • Chittari \latinet al. 2016 Chittari, B. L.; Park, Y.; Lee, D.; Han, M.; MacDonald, A. H.; Hwang, E.; Jung, J. Phys. Rev. B: Condens. Matter Mater. Phys. 2016, 94, 184428.
  • Seoánez \latinet al. 2007 Seoánez, C.; Guinea, F.; Castro Neto, A. H. Phys. Rev. B: Condens. Matter Mater. Phys. 2007, 76, 125427.
  • Kim and Park 2009 Kim, S. Y.; Park, H. S. App. Phys. Lett. 2009, 94, 101918.
  • von Oppen \latinet al. 2009 von Oppen, F.; Guinea, F.; Mariani, E. Phys. Rev. B: Condens. Matter Mater. Phys. 2009, 80, 075420.
  • Chen \latinet al. 2009 Chen, C.; Rosenblatt, S.; Bolotin, K.; Kalb, W.; Kim, P.; Kymissis, I.; Stormer, H.; Heinz, T.; Hone, J. Nat. Nanotech. 2009, 4, 861–867.
  • Zande \latinet al. 2010 Zande, A. M. v. d.; Barton, R. A.; Alden, J. S.; Ruiz-Vargas, C. S.; Whitney, W. S.; Pham, P. H. Q.; Park, J.; Parpia, J. M.; Craighead, H. G.; McEuen, P. L. Nano Lett. 2010, 10, 4869–4873.
  • Song \latinet al. 2012 Song, X.; Oksanen, M.; Sillanpää, M. A.; Craighead, H. G.; Parpia, J. M.; Hakonen, P. J. Nano Lett. 2012, 12, 198–202.
  • Takamura \latinet al. 2016 Takamura, M.; Okamoto, H.; Furukawa, K.; Yamaguchi, H.; Hibino, H. Micromachines 2016, 7, No. 158.
  • Miao \latinet al. 2014 Miao, T.; Yeom, S.; Wang, P.; Standley, B.; Bockrath, M. Nano Lett. 2014, 14, 2982–2987.
  • Weber \latinet al. 2016 Weber, P.; Güttinger, J.; Noury, A.; Vergara-Cruz, J.; Bachtold, A. Nat. Commun. 2016, 7, 12496.
  • Golding \latinet al. 1978 Golding, B.; Graebner, J. E.; Kane, A. B.; Black, J. L. Phys. Rev. Lett. 1978, 41, 1487–1491.
  • Phillips 1987 Phillips, W. A. Rep. Prog. Phys. 1987, 50, 1657 – 1708.
  • Fefferman \latinet al. 2008 Fefferman, A. D.; Pohl, R. O.; Zehnder, A. T.; Parpia, J. M. Phys. Rev. Lett. 2008, 100, 195501.
  • Lee \latinet al. 2008 Lee, C.; Wei, X.; Kysar, J. W.; Hone, J. Science 2008, 321, 385–388.
  • Barmatz \latinet al. 1975 Barmatz, M.; Testardi, L. R.; Di Salvo, F. J. Phys. Rev. B: Condens. Matter Mater. Phys. 1975, 12, 4367–4376.
  • Jiang \latinet al. 2004 Jiang, H.; Yu, M.-F.; Liu, B.; Huang, Y. Phys. Rev. Lett. 2004, 93, 185501.
  • Hüttel \latinet al. 2009 Hüttel, A. K.; Steele, G. A.; Witkamp, B.; Poot, M.; Kouwenhoven, L. P.; van der Zant, H. S. J. Nano Lett. 2009, 9, 2547–2552.
  • Imboden and Mohanty 2014 Imboden, M.; Mohanty, P. Physics Reports 2014, 534, 89 – 146.
  • Anderson \latinet al. 1972 Anderson, P. W.; Halperin, B. I.; Varma, C. M. Philos. Mag. 1972, 25, 1–9.
  • Phillips 1972 Phillips, W. A. J. Low Temp. Phys. 1972, 7, 351–360.
  • Jäckle 1972 Jäckle, J. Z. Phys. 1972, 257, 212–223.
  • Kleiman \latinet al. 1987 Kleiman, R. N.; Agnolet, G.; Bishop, D. J. Phys. Rev. Lett. 1987, 59, 2079–2082.
  • Venkatesan \latinet al. 2010 Venkatesan, A.; Lulla, K. J.; Patton, M. J.; Armour, A. D.; Mellor, C. J.; Owers-Bradley, J. R. Phys. Rev. B: Condens. Matter Mater. Phys. 2010, 81, 073410.
  • Rivière \latinet al. 2011 Rivière, R.; Deléglise, S.; Weis, S.; Gavartin, E.; Arcizet, O.; Schliesser, A.; Kippenberg, T. J. Phys. Rev. A: At., Mol., Opt. Phys. 2011, 83, 063835.
  • Esquinazi 2013 Esquinazi, P. Tunneling systems in amorphous and crystalline solids; Springer Science & Business Media: Berlin, 2013; p 187.
  • Lulla \latinet al. 2013 Lulla, K. J.; Defoort, M.; Blanc, C.; Bourgeois, O.; Collin, E. Phys. Rev. Lett. 2013, 110, 177206.
  • Faust \latinet al. 2014 Faust, T.; Rieger, J.; Seitner, M. J.; Kotthaus, J. P.; Weig, E. M. Phys. Rev. B: Condens. Matter Mater. Phys. 2014, 89, 100102.
  • Coppersmith and Golding 1993 Coppersmith, S. N.; Golding, B. Phys. Rev. B: Condens. Matter Mater. Phys. 1993, 47, 4922–4936.
  • Singh \latinet al. 2016 Singh, V.; Shevchuk, O.; Blanter, Y. M.; Steele, G. A. Phys. Rev. B: Condens. Matter Mater. Phys. 2016, 93, 245407.
  • Remus \latinet al. 2009 Remus, L. G.; Blencowe, M. P.; Tanaka, Y. Phys. Rev. B: Condens. Matter Mater. Phys. 2009, 80, 174103.
  • Forschungszentrum Jülich GmbH. HNF - Helmholtz Nano Facility 2017 Forschungszentrum Jülich GmbH. HNF - Helmholtz Nano Facility, Journal of large-scale research facilities 2017, 3, A112.

Appendix A Supporting information

a.1 Characterization of the membrane by Raman

Figure 4: (a) Image of the device with the red dot marking the position of the Raman measurement, where both the top and bottom graphene flake are present. (b) Raman measurement of the membrane with single Lorentzian fits of the 2D and G line showing that top and bottom graphene flake are single layer.

Raman measurements are conducted to determine the number of layers for the encapsulating graphene after the measurement in the dilution refrigerator. Figure 4a shows the position of the suspended region that is investigated with Raman. In Figure 4b the single Lorentzian fit of the 2D peak to the Raman spectrum with a at indicates single layer graphene. A Raman map of the entire membrane shows that the drum collapsed after the measurement in the dilution refrigerator with only the region around the position marked by the red dot still suspended. The Raman spectrum of the NbSe layer in the drum area is also investigated. However, the resolution of the Raman signal is insufficient to determine the width of the shear mode, which perishes in the Rayleigh peak. The soft mode, A and E can be determined and show a Raman shift in the expected region around (not shown) 48.

a.2 Calculation of the membrane amplitude

The amplitude of the membrane can be estimated from the optomechanical coupling as shown in 49, 50:


Here and are the coupling constants of the cavity, are calculated from the amplifiers and attenuators, is the power of the pump wave and is the change of the cavity resonance frequency with displacement of the membrane. The output power is measured, whereas the other quantities are extracted from other calibration measurements as explained in the following.

The of the measurement setup inside the cryostat system is . Adding this to the whole measurement setup outside of the cryostat with a total equals to a total . However, for the measurement with the weak port no bandpass filters with a combined attenuation of are needed resulting in a total . This attenuation also leads to a reduction of the drive amplitude. The drive RMS voltage at the sample is adjusted to the applied back gate voltage to avoid high forces resulting in a movement of the membrane. We modify the applied drive voltage at the sample so that the drive force is independent of the applied back gate voltage .

The read out technique utilizes red sideband pumping of a microwave cavity. The cavity linewidth is determined by the capacitive coupling of the weak and strong port for feed in and read out of the pump frequencies 50. The total cavity linewidth of is composed of the capacitive coupling to the two ports , and additional internal losses . The coupling constants are extracted from fits to transmission measurements. All measurements are conducted in a dilution refrigerator with a base temperature of well below the critical temperatures of the Nb cavity ().

The single phonon-photon coupling rate is estimated by 50


Here is the cavity resonance frequency, is the total capacitance of the cavity, z is the zero point motion of the resonator and is the derivative of the capacitance between the mechanical resonator and the gate with respect to the center position of the membrane . For we obtain from a plate capacitor estimation , yielding in close agreement with () extracted from the change of the cavity resonance frequency by pulling on the membrane with .

Due to the limited coupling we are not able to confirm the estimation of with thermal calibration measurements. However, the coupling is compatible with - from previous samples with similar geometry 49, 50.

a.3 Temperature dependence of dissipation in a high Q few-layer graphene and an hBN encapsulated graphene device

Figure 5: Temperature dependence of and in a high Q few-layer graphene device (a, b) and an encapsulated graphene device (c, d). (a, c) The dotted dashed line corresponds to with extracted from a linear fit between and (a) and from a fit between and (c). (b, d) Temperature dependence of measured with the ring-down method. The dashed line in (b) corresponds to the saturation expected in the standard tunneling model and the dotted line is a fit between and to . (d) The dotted line represents a fit between and .

Figure 5 shows the temperature dependence of and in a high Q few-layer graphene resonator previously studied in Refs. 49, 51 and an hBN/graphene/hBN resonator. The resonance frequency change in Figure 5a is extracted from thermal motion measurements because the frequency shift is not resolved in the ring-down measurements. Above the data follows a logarithmic temperature dependence expected from resonant interaction with TLS. The deviation below is at least partially attributed to heating of the membrane due to the increased cavity pump power needed to resolve the thermal motion 49. The shift in resonance frequency shown in Figure 5c also shows a deviation from the otherwise logarithmic behaviour below which is connected to an insufficient thermalization of the device to the mixing chamber plate as mentioned in the main text. The extracted value for is in the expected range.
The temperature dependence in Figure 5b shows again a crossover from a stronger to a weaker temperature dependence upon increasing the temperature. The dashed line corresponds to the saturation level expected from the standard TLS model with from Figure 5a. However, as with the heterostructure device this plateau does not describe the data well. Besides having a too small value the measurement does not show an immediate saturation but rather a dependence of above . Such a dependence has been observed in several graphene based resonators above  52. To still find an estimate for it was calculated using the previously extracted value for yielding . Another difference compared to the heterostructure device in the main text is the absence of a clear sign for a resonant dissipation contribution to at the lowest temperatures. In the hBN/graphene/hBN device the effect of saturation due to resonant absorption and insufficient thermalization in combination with a lower resonance frequency of around , which should shift the position of the change in sign to lower temperatures, is not distinguishable.
It will be interesting to investigate this regime in more devices and for higher mechanical modes to gain more insight.
The inverse Quality factor dependence of the encapsulated device (Figure 5d) can be described by a fit for which still indicates, that the dissipation is connected to electrons rather than phonons.

a.4 Dependence of the dissipation on mechanical drive voltage

In the measured range the heterostructure device does not show the expected saturation of the dissipation when increasing the drive voltage at but rather shows a decrease of as shown in Figure 7. As pointed out in the main text non saturable attenuation might be caused by heavily asymmetric TLS and modified scattering cross-sections. The influence of individual TLS with deviating energies and tunneling splittings from the assumed contribution might also be increased due to the tiny membrane volume. The unexpected behaviour around in the temperature dependent data (Figure 6) might be an additional indication for that.

Figure 6: The measured as function of different drive voltages and measurement times. Both the measurements are consistent within the error bars and were taken coming from low temperature and heating up the sample. Both measurements show data at and . However, the first measurement (blue data points) uses a higher drive of and faster time constant of . The second measurements (red data points) have been taken at and an acquisition time constant of .
Figure 7: measured at base temperature as function of drive voltage. No saturation of the damping with increased drive voltage is observed.

a.5 Effect of thermal expansion on the mechanical resonance frequency

One possible explanation for the observed discrepancy of the frequency shift between the TLS theory and the measured data above in favor of lower frequencies is the influence of the negative thermal expansion of the graphene layers. Assuming that the mechanical properties of the multilayer are dominated by the graphene based on the higher Youngs modulus one can roughly estimate the frequency shift due to the change in strain of the membrane. Using a value of /K for the thermal expansion coefficient (slightly below room temperature value 53) and results in the red curve in Fig. 8. We point out that neither the expansion of the exact resonator drum was regarded nor a deeper investigation in the actual size of this contribution was conducted and the analysis serves only to give a rough estimation of the effect.

Figure 8: Temperature dependence of the relative shift in mechanical resonance frequency with temperature. The red curve illustrates the influence of an additional contribution to the frequency shift caused by variation of the strain in the membrane due to the negative thermal expansion of the graphene layers.


  • Xi \latinet al. 2015 Xi, X.; Zhao, L.; Wang, Z.; Berger, H.; Forro, L.; Shan, J.; Mak, K. F. Nat. Nanotech. 2015, 10, 765–769.
  • Weber \latinet al. 2016 Weber, P.; Güttinger, J.; Noury, A.; Vergara-Cruz, J.; Bachtold, A. Nat. Commun. 2016, 7, 12496.
  • Weber \latinet al. 2014 Weber, P.; Güttinger, J.; Tsioutsios, I.; Chang, D. E.; Bachtold, A. Nano Lett. 2014, 14, 2854–2860.
  • Güttinger \latinet al. 2017 Güttinger, J.; Noury, A.; Weber, P.; Eriksson, A. M.; Lagoin, C.; Moser, J.; Eichler, C.; Wallraff, A.; Isacsson, A.; Bachtold, A. Nat. Nanotech. 2017, 12, 631–636.
  • Imboden and Mohanty 2014 Imboden, M.; Mohanty, P. Physics Reports 2014, 534, 89 – 146.
  • Yoon \latinet al. 2011 Yoon, D.; Son, Y.-W.; Cheong, H. Nano Lett. 2011, 11, 3227–3231.
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description