Cavity optomechanics mediated by a quantum two-level system

This is an electronic reprint of the original article. This reprint may differ from the original in pagination and typographic detail. Author(s): Title: Year: Version: Please cite the original version: All material supplied via JYX is protected by copyright and other intellectual property rights, and duplication or sale of all or part of any of the repository collections is not permitted, except that material may be duplicated by you for your research use or educational purposes in electronic or print form. You must obtain permission for any other use. Electronic or print copies may not be offered, whether for sale or otherwise to anyone who is not an authorised user. Cavity optomechanics mediated by a quantum two-level system Pirkkalainen, J.-M.; Cho, S.U.; Massel, Francesco; Tuorila, J.; Heikkilä, Tero; Hakonen, P.J.; Sillanpää, M.A. Pirkkalainen, J.-M., Cho, S.U., Massel, F., Tuorila, J., Heikkilä, T., Hakonen, P.J., & Sillanpää, M.A. (2015). Cavity optomechanics mediated by a quantum two-level system. Nature Communications, 6, Article 6981. https://doi.org/10.1038/ncomms7981 2015 ARTICLE Received 12 Feb 2015 |Accepted 19 Mar 2015 |Published 27 Apr 2015 Cavity optomechanics mediated by a quantum two-level system J.-M. Pirkkalainen1, S.U. Cho2,w, F. Massel3, J. Tuorila4, T.T. Heikkila ¨3, P.J. Hakonen2& M.A. Sillanpa ¨a ¨1 Coupling electromagnetic waves in a cavity and mechanical vibrations via the radiation pressure of photons is a promising platform for investigations of quantum–mechanical properties of motion. A drawback is that the effect of one photon tends to be tiny, and hence one of the pressing challenges is to substantially increase the interaction strength. A novel scenario is to introduce into the setup a quantum two-level system (qubit), which, besides strengthening the coupling, allows for rich physics via strongly enhanced nonlinearities. Here we present a design of cavity optomechanics in the microwave frequency regime involving a Josephson junction qubit. We demonstrate boosting of the radiation–pressure interaction by six orders of magnitude, allowing to approach the strong coupling regime. We observe nonlinear phenomena at single-photon energies, such as an enhanced damping attributed to the qubit. This work opens up nonlinear cavity optomechanics as a plausible tool for the study of quantum properties of motion. DOI: 10.1038/ncomms7981 OPEN 1Department of Applied Physics, Aalto University, PO Box 11100, FI-00076 Aalto, Finland. 2Department of Applied Physics, Low Temperature Laboratory, Aalto University, PO Box 15100, FI-00076 Aalto, Finland. 3University of Jyva ¨skyla ¨, Department of Physics, Nanoscience Center, University of Jyva ¨skyla ¨,PO Box 35 (YFL) FI-40014 University of Jyva ¨skyla ¨, Jyva ¨skyla ¨, Finland. 4Department of Physics, University of Oulu, FI-90014 Oulu, Finland. wPresent address: Department of Physics and Astronomy, Seoul National University, Seoul 151-747, Korea. Correspondence and requests for materials should be addressed to M.A.S. (email: Mika.Sillanpaa@aalto.fi). NATURE COMMUNICATIONS | 6:6981 | DOI: 10.1038/ncomms7981 | www.nature.com/naturecommunications 1 &2015 Macmillan Publishers Limited. All rights reserved. The basic setup of cavity optomechanics1,2 involves an otherwise ordinary optical Fabry–Perot cavity resonator, but having one of the mirrors compliant for mechanical vibrations. The oscillatory displacement xof the mirror couples to the length that corresponds to the frequency o c of the cavity. This interaction can also be interpreted as that resulting from the radiation pressure caused by the cavity photons and characterized by the coupling energy g 0 (qo c /qx)x zp . Here, x zp is the zeropoint motion amplitude. Cavity optomechanics has recently received a lot of attention, for it has turned out as an excellent test system for investigations on quantum–mechanical properties of the motion of nearly macroscopic bodies3,4. Several different physical systems have been demonstrated, which realize the radiation-pressure interaction between electromagnetic and mechanical modes. In particular, the frequency of the electromagnetic cavity need not be optical. One of the most successful realizations in the recent past has been the use of microwave-regime superconducting cavities made with lithographic techniques3,5,6. Here, the displacement changes the effective capacitance C(x) of the cavity and hence its resonance frequency o c , and the coupling energy is g 0 ¼(o c /2C)(qC/qx)x zp . In all the physical realizations, except in the case of microscopic motional degrees of freedom7,8, the radiationpressure coupling energy is much smaller than the cavity decay rate k. Therefore, the cavity has to be strongly irradiated up to a high photon number n c c1 in order to enhance the intrinsically small interaction. This, however, linearizes the interaction and creates two linearly coupled harmonic oscillators, which is a somewhat limited system. Hence, one nowadays strives to enhance the magnitude of the interaction to exploit the intrinsic nonlinearities9–14. The next goal is to reach a strong coupling threshold g 0 ]kwhere some quantum phenomena become observable9,11,13,15,16. The largest radiation-pressure coupling values g 0 /2pC1 MHz have been demonstrated in the optical frequency domain (ref. 4). These systems also possess the largest ratio of the coupling versus cavity decay k, which is, g 0 /kB0.2%. This is still two and half orders of magnitude away from strong coupling. Here we build on the scheme proposed in ref. 17. The scheme is motivated by several successful experiments coupling superconducting Josephson junction quantum circuits to micromechanical resonators18–22, although canonical radiation– pressure phenomena have not been observed in those setups. We demonstrate a microwave-regime optomechanical device that allows for high values of g 0 /k, and at the same time, displays phenomena not previously observed in cavity optomechanics either in the microwave or the optical frequency range. Results The basic idea. The device consists of a tripartite system made with patterned thin films of superconducting aluminium on a silicon chip, containing first of all a microwave resonator (called cavity in the following). Another part is a micromechanical resonator and the third is a charge quantum two-level system (qubit), which consists of small Josephson tunnel junctions and mediates the interaction between the former two parts (see Fig. 1 and Supplementary Fig. 1). The bare cavity and the mechanical modes have the frequencies o0 cand o0 m, and their dynamical variables are the creation and annihilation operators aw,aand bw, b, respectively. The qubit is represented by the Hamiltonian H QB ¼B x s x /2 B y s y /2 B z s z /2. Here, s j are Pauli matrices and B j are the controllable pseudomagnetic fields, which also depend on the cavity and mechanical coordinates. The level spacing of the qubit is B¼ffiffiffiffiffiffiffiffiffiffiffiffi PjB2 j q. More discussion of the Hamiltonian is found in Supplementary Methods. A fundamental reason to expect rich optomechanical physics in the scheme owes to the nonlinearity of the quantum two-level system17–25. An important ingredient here is the large d.c. voltage bias V g applied to the mechanical resonator represented as a movable capacitance C g (x). If coupled to a two-level system (qubit), this creates a longitudinal coupling g m (bwþb)s z in the natural basis of the charge qubit, with a maximal effect on the energy of the qubit by the mechanical resonator. When coupling this qubit to a cavity with the transverse interaction g c (awþa)s x , the cavity frequency will experience a Stark shift g2 c=o0 cB  . This Stark shift contains a contribution from the vibrations due to their effect on qubit level spacing. The system then closely approximates the radiation-pressure interaction between the effective cavity mode (frequency o c ) and the effective mechanical mode (frequency o m ), which both are Stark shifted from the bare frequencies. The radiation-pressure coupling strength between these modes can be amplified in this setting by a factor of the order of CV g /(2e), which amounts to B4to6 orders of magnitude for typical experimental parameters (ref. 26). Moreover, there are intriguing nonlinear corrections to the basic setup, first of all, a cross-Kerr type interaction in the next higher order. Most importantly, the damping of the mechanics becomes affected by the nonlinearities. The physical device (see Fig. 1c and Supplementary Fig. 2) is fabricated by standard electron-beam lithography. The cavity is a 5-mm long superconducting microstrip resonating at the (bare) cavity frequency o0 c=2p¼4:93 GHz. The charge qubit has the island capacitance C S ¼C g þC 1 þC 2 , which gives the charging energy of a single electron E C ¼e2/2C S . It is of the same order as x Vg CL Φ Cc Cg(x) Vin Vout EJ1,C1 EJ2,C2 gc gm To cavity ω P Figure 1 | Optomechanics with a qubit. (a) Illustration showing the idea for introducing a two-level system (qubit or artificial atom) into a cavity optomechanical setup. Inside the cavity (blue) there is a quantum two-level system (green), which is mechanically compliant (red). (b) Equivalent electrical circuit of the microwave design with the same colour codes as in a. The parameters in the equivalent circuit are CC0.33 pF, LC3.2 nH, C g C1.8 fF, C S C4.9 fF, E J ¼E J1 þE J2 C10.5 GHz, E C C5.3 GHz. The reflected pump microwave V out is detected outside the cryostat. The Josephson junctions forming the qubit are marked by crossed boxes, and the mechanical resonator is coupled to the qubit island via a movable capacitance C g (x). (c) Micrographs of the device showing the qubit and the mechanical resonator. The Josephson junctions are marked by arrows. Scale bar, 2 mm. ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms7981 2NATURE COMMUNICATIONS | 6:6981 | DOI: 10.1038/ncomms7981 | www.nature.com/naturecommunications &2015 Macmillan Publishers Limited. All rights reserved. the Josephson energy E J , which is E J /E C C2.0. A strong qubitcavity coupling with the coupling energy up to B500 MHz (depending on the qubit bias) sets the system near the ultrastrong coupling27 of circuit QED, and hence gives rise to a pronounced Stark shift of the cavity. The micromechanical resonator is made as a 2-mm wide bridge suspended across a 50 nm vacuum gap atop the qubit island (ref. 21) and resonating at the frequency o m /2pC65 MHz of the lowest flexural mode. The qubitmechanics coupling varies here between g m /2pC80 MHzy 160 MHz depending on the value of V g . We note that in spite of the large coupling exceeding the bare mechanical frequency, the Schrieffer–Wolff approach of deriving the radiation-pressure interaction is justified (see Supplementary Methods). Measurements. In all the measurements of the mechanical motion we use the basic idea of cavity optomechanics shown in Fig. 2a, where the cavity is irradiated by the pump microwave at a frequency o p , which usually resides near the red sideband o c o m . The pump detuning D(o p o c )/o m is hence around DB1. The motional sidebands at ±o m about the pump give information about the mechanics. In the following, the sideband closer to the cavity is studied. The main result of the radiationpressure interaction is the ‘optical spring effect’ which modifies both the mechanical frequency and damping. The changes exerted are denoted as doopt m, and gopt m4g2 0nc=k, respectively. The experiments are performed in a dilution refrigerator at 20 mK temperature. We first identify the mechanical resonance by pumping the cavity very strongly up to n c B109such that the Josephson effect, which has its influence only in the low-photon regime n c B1, becomes insignificant, and what is left behind is a linear cavity with the bare frequency o0 c. The estimated bare optomechanical coupling between the bare cavity and the mechanics is only g 0 /2pB1 Hz. Such a small coupling does not yet allow for observing the small thermal motion of the mechanical resonator at a low temperature. We hence actuate the motion via the gate line, up to about 10-pm displacement amplitude, obtaining the resonance curve shown in Supplementary Fig. 3. This way we determine the intrinsic mechanical frequency o m /2pC65 MHz and linewidth that is g m /2p¼15 kHz. Notice that o m decreases when V g deviates from zero due to the well-known electrostatic softening of the mechanical spring. Next, we decrease the pump power by about 10 orders of magnitude down to P p B1 fW. Although power levels are drastically reduced, the detection sensitivity is enhanced and we can observe the sideband peak corresponding to the thermal motion, shown in Fig. 2b, at 20 mK temperature, amounting to about six phonons at equilibrium. The peak frequency is somewhat shifted from the calibration experiment, and the linewidth of this peak is about twice the intrinsic linewidth. As discussed below, this results from the qubit-induced nonlinearities. We now turn into careful characterization of the novel type of optomechanical interaction. We first show that the device displays the basic physics due to the radiation-pressure interaction, namely, the optical spring effect. As discussed below, the latter causes the mechanical frequency and damping to differ from the intrinsic values. We now characterize the basic behaviour of the microwave cavity. As seen in Fig. 2c, the frequency of the effective cavity depends periodically on the gate charge n g ¼C g V g /2e. Each period corresponds to the tunnelling of one extra Cooper pair onto the qubit island. Here we have a large dc voltage V g B5y10 V, which is needed to establish the electromechanical interaction, hence n g 4104. However, because of periodicity, we henceforth refer to n g via its non-integer part for simplicity. On top of the large d.c. voltage, we use a fine tuning voltage to adjust the radiation-pressure physics by tuning the gate charge within one period. The period of the cavity response is one electron (0.5) instead of two as expected. This is due to a quasiparticle slowly jumping on and off the island21,28. We thus see a double image, shifted by one electron, of the cavity transition. These are marked by the red and blue dashed lines in Fig. 2c. At a given n g , there are hence two different cavity frequencies. The system switches between these states while the quasiparticle randomly tunnels on and off the island at the typical rates well below a MHz (refs 29,30). These rates are much slower than the mechanical Vout/Vin (dB) 4.8 4.9 5.0 f – fP (MHz) Spectral density (arb) 65.1 65.2 65.3 65.4 65.5 65.6 m c m pp + m 0p – m ng fP (GHz) 4.8 4.85 4.9 0.5 0.0 –0.5 ng 0.5 0.0 –0.5 0 –10 –5 –15 –20 65.2 65.4 f – fP (MHz) Power (arb) 0.5 0.6  Figure 2 | Basic characterization. (a) The frequencies involved: the vertical axis gives the spectrum of either the cavity (o c ) or the mechanical resonator (o m ). The pump (frequency o p ) is applied around the detuning DC1 below the cavity frequency. (b) Thermal motion measured at 20 mK. The parameters are DC1, n c C0.05, n g ¼0.25, V g ¼6.5 V. (c) Gate charge modulation of the cavity resonance absorption. The resonance o c (n g ) of the effective cavity appears as the periodic black-blue line. The dashed red and blue lines are theoretical fits to two of the branches, obtained by using the qubit-cavity Hamiltonian17 in numerical diagonalization. The flux bias was F/F 0 C0.39. The vertical green and grey lines represent the frequencies o p and the sideband o p þo m , respectively, used to obtain the data in Fig. 3a,b. (d) Spectral density of the emission from the cavity around the mechanical sideband. NATURE COMMUNICATIONS | DOI: 10.1038/ncomms7981 ARTICLE NATURE COMMUNICATIONS | 6:6981 | DOI: 10.1038/ncomms7981 | www.nature.com/naturecommunications 3 &2015 Macmillan Publishers Limited. All rights reserved. frequency, and hence in the long-time average we see responses from both cavity states superimposed on top of each other. From the geometry we find qC g /qxC24 nF/m and x zp C6 fm. We determine the coupling from the experimental data using g 0 ¼(qo c /qx)x zp , where (qo c /qx)¼V g /2e(qo c /qn g )(qC g /qx) can be related to the slope of the gate charge modulation of the cavity frequency seen in Fig. 2c. For n g ¼0.25 where the cavity is reasonably well visible, we have (qo c /qn g )C(2p)250 MHz and hence g 0 /2pC0.5 MHz when V g ¼4.6 V. Since |g 0 (n g ¼0.25)| ¼ |g 0 (n g ¼0.75)|, this value is the same for both cavity branches. Gate charge dependence. One can characterize the radiationpressure physics as a function of, for example, the value of the gate charge n g . Changes in n g change the value of g 0 , which is proportional to the slope (qo c /qn g ). However, a change in n g has also other consequences. With a given pump frequency, the pump detuning Dchanges according to the effective cavity frequency. Similarly, also the effective mechanical frequency o m will change (due to Stark shift) with n g . Finally, these effective modes exhibit the radiation-pressure interaction visible as further changes in frequency and damping due to the optical spring effect. We denote the resulting total frequency of the mechanics as otot m¼omþdoopt m. The total damping is given by gtot m¼gmþgopt mþgf m, where the last term gf mis due to the qubit and is discussed below. In Fig. 2d we show a result of this kind of a measurement, allowing us to connect the thermal emission peak with the cavity response. A maximum emission occurs at n g C0.25, which here coincides with the sideband resonance and hence the expected maximum signal. The two crossing branches are associated with the two cavity branches. We next look in detail at the frequency and linewidth of the thermal motion peak in the measurement as in Fig. 2d. The extracted parameters are displayed in Fig. 3a,b revealing three distinct sets of peaks. To understand this data, we first calculate the Stark-shifted o m and o c , as well as g 0 from the diagonalization of the qubit-cavity Hamiltonian (see Supplementary Methods) and then apply the results in refs 17,31 to obtain a prediction for the optical spring. Each of the two cavity branches causes its own sideband response. These arise via the detuning and coupling shown in Fig. 3c,d. The blue curves in reality originate from n g values between 0.25...0.5, but reflect to the left about n g ¼0.25 due to the quasiparticle. The data and theory are marked with the compatible red and blue colour codes in Figs 2 and 3 and Supplementary Fig. 4. We obtain a good agreement between the experiment and theory for the mentioned two branches. However, the peaks visible around n g ¼0...0.15, marked in black, are not directly explained by this model. We attribute them to that the qubit spends some part of the time in the excited state as observed in this kind of measurements21. The excited state gives rise to another cavity frequency above o c , having the highest qo c /qn g around n g ¼0. The qubit being excited, however, is beyond the current model and we cannot make a quantitative comparison to theory. Damping due to the qubit. In the rest of the paper, we select a fixed value n g ¼0.25, which simplifies the analysis since there is only one cavity branch and a constant g 0 . Another way to examine the radiation-pressure interaction is to change the pump detuning (Fig. 3e). The total damping is maximized close to the sideband resonance in agreement with theory. The total damping should increase linearly with the pump photon number n c owing to the contribution gopt m, moreover, it should depend quadratically on g 0 at a given n c . We test this in Fig. 4 by varying the photon number at a few fixed values of g 0 , which are set by changing the gate voltage V g . The data points are extracted from the data sets in Supplementary Fig. 5, which we here show as an example of thermal motion peaks. As observed in Fig. 4, the damping towards n c -0 is clearly higher than the intrinsic g m . The extra damping is attributed to that the qubit opens another dissipation channel due to the hybridization of the qubit and the mechanics17, analogously to the Stark shift, and is given by gf m’gQBg2 mB2 x B3om :ð1Þ fm tot (MHz) m tot (kHz) m tot (kHz) 0 0.05 0.1 0.15 0.2 0.25 0 2 4 6 x 10 5 g0 (Hz) n g n g 0 0.05 0.1 0.15 0.2 0.25 65.1 65.2 65.3 65.4 65.5 0 0.05 0.1 0.15 0.2 0.25 0 10 20 30 40 50 60 70 80 n g 8 4.76 4.78 4.8 4.82 0 20 40 60 80 –1 0 1 2 0 0.05 0.1 0.15 0.2 0.25 n g Δ/ωm f P (GHz) Figure 3 | Optomechanics enhanced by a qubit. (a,b) Optical spring effect as a function of gate charge. The solid lines are a fit to theory and are plotted in the region where detuning and coupling are such that the peaks are visible. The red and blue colours refer to the two cavity branches as in Fig. 2. The black circles present data that presumably originate from higher excited levels of the qubit (see main text). The pump frequency o p /2p¼4.79 GHz, and the sideband (o p þo m )/2p¼4.855 GHz. (c) Illustration of the pump detuning for the two branches for the measurements in aand b.(d) Radiation-pressure coupling extracted from o c (n g ). (e) Optical spring as a function of pump detuning, at fixed gate charge n g ¼0.25. The solid line is a fit to theory. In all data, the cavity photon number is in the range 0.1...1. This is estimated based on the input power and cable attenuation, and in the end fine-tuned as an adjustable parameter. In all the panels, V g ¼4.6 V. ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms7981 4NATURE COMMUNICATIONS | 6:6981 | DOI: 10.1038/ncomms7981 | www.nature.com/naturecommunications &2015 Macmillan Publishers Limited. All rights reserved. Here, g QB is the qubit relaxation rate. We could not independently measure it, but a good fit of the data is obtained with a reasonable value g QB C(2p)60 MHz. Although the damping due to the basic cavity cooling gopt m  gives rise to back-action sideband cooling of the mechanical resonator, further theoretical work is needed to understand the effect of the qubit contribution gf mas well as the nonlinearities such as the cross-Kerr effect on the final temperature. In Fig. 4 we see that the damping grows first linearly up to about n c B1 and then saturates. The saturation of the radiationpressure physics around n c B1 is due to reaching the limits of the linear regime of the effective cavity. This is set by the requirement that the total phase excursions ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ncþ1=2 pfzp should be small as compared to 2pin order not to hit the nonlinearity of the Josephson cosine function. This is analogous to the transition of a trapped ion out from the Lamb–Dicke regime (see Supplementary Methods). Here, the phase zero-point fluctuation is f zp B0.2, which agrees with the scale we observed. Discussion The largest g 0 /2pC1.6 MHz obtained in this work is achieved by half of a flux quantum flux bias, and setting a high gate voltage V g ¼9.5 V. This value is the largest ever reported in a cavity optomechanical setting. This was achieved by integrating, for the first time, a quantum two-level system into such a setup, allowing boosting of the radiation-pressure coupling by six orders of magnitude. The ratio of coupling to the cavity linewidth, B4%, can likely be further increased by designing the coherence properties of the qubit. We estimate that boosting factors up to eight orders of magnitude, giving a radiation-pressure coupling up to 100 MHz, are in principle possible with an optimized device (see Supplementary Fig. 6). The setup introduces long-sought nonlinearities into the framework of cavity optomechanics and enables their utilization, for example, in Fock-state measurements of the phonon states, or for tests of the foundations of quantum mechanics32–34. Methods Chip layout.A photograph of the chip is shown in Supplementary Fig. 2. The ends of the meandering inductance are close to each other such that the double-junction charge qubit can be placed to connect the ends of the meander, which is required for the scheme. The device is fabricated on a regular silicon substrate using a similar process as in ref. 21. The device was designed to not have a too small cavity capacitance for the system to reside in the Lamb–Dicke limit for qubit-cavity coupling (see discussion in Supplementary Methods). Most of the capacitance is provided by a bonding pad-like structure visible between the coupling capacitance and the meandering inductance, marked Cin Supplementary Fig. 2. Experimental setup.The measurements are done at the base temperature of a dilution refrigerator using the cabling setup as shown in Supplementary Fig. 7. The pump microwave reflects from the cavity and is amplified at a 4K stage using a low-noise amplifier. The system noise temperature is in the range of 3...4 K. Filtering the room temperature noise away from reaching the sample via the cable connected to the mechanical resonator (that is, qubit gate) is challenging because the typical solution (resistive attenuators) cannot be used owing to the high voltages involved. We hence use reflective filtering. From the qubit point of view as well, the filter needs to be reactive, because dissipative 50 Oenvironment could cause too high spontaneous decay. The filtering task is carried out by a combined bias-T and low-pass filter, marked yellow in Supplementary Fig. 7. Without the filter, we observe substantial excitation of the qubit to the higher levels, such that the spectrum of the effective cavity changes qualitatively and cannot be modelled by ground-state transitions. With the filter installed, the excitation is modest, and cavity response well understood. The filter provides enough transmission (33 dB) at frequencies around o m such that we can actuate the mechanics if needed. Cavity characterization.The cavity frequency has a strong dependence on the magnetic flux Fapplied through the loop, as shown in Supplementary Fig. 8, which displays the microwave reflection response quite similar to Fig. 2c in the main text. The cavity resonance appears as the yellow-blue coloured dip varying between 4.99 and 4.84 GHz, periodically in flux with the period of one flux quantum F 0 .Inthe measurement, we applied noise to the gate such that the response is a weighted average over all the possible gate n g values. In the range B0.4yn g y0.6 the resonance is not well visible because of the strong gate dispersion. The theoretical curves plotted on top of the data in Supplementary Fig. 8 are obtained by numerical diagonalization of Supplementary Equation 9, but disregarding the mechanics (that is, n g (x)-n g ), which here is insignificant. The overlaid curves indicate a good match between the experiment and the modelling. The cavity linewidth khas a strong dependence on the qubit gate charge, as displayed in Supplementary Fig. 9. While the changes are due to the qubit, we cannot make a definite conclusion of whether it is caused by qubit relaxation or dephasing. We also checked that simulating the inhomogenous broadening by establishing the theory with bare k, and selecting a ensemble of f c , and averaging over the ensemble does not markedly change the results. References 1. Kippenberg, T. J. & Vahala, K. J. Cavity optomechanics: Back-action at the mesoscale. Science 321, 1172–1176 (2008). 2. Aspelmeyer, M., Kippenberg, T. J. & Marquardt, F. Cavity optomechanics. Rev. Mod. Phys. 86, 1391–1452 (2014). 3. Teufel, J. D. et al. Sideband cooling of micromechanical motion to the quantum ground state. Nature 475, 359–363 (2011). 4. Chan, J. et al. Laser cooling of a nanomechanical oscillator into its quantum ground state. Nature 478, 89–92 (2011). 5. Regal, C. A., Teufel, J. D. & Lehnert, K. W. Measuring nanomechanical motion with a microwave cavity interferometer. Nat. Phys. 4, 555–560 (2008). 6. Massel, F. et al. Microwave amplification with nanomechanical resonators. Nature 480, 351–354 (2011). nc 0 1 2 3 4 5 6 0 50 100 150 200 250 300 350 Vg= 6.5 V Vg= 4.6 V Vg= 9.5 V 0 0.5 1 1.5 2 2.5 0 20 40 60 80 100 nc m tot (kHz) m tot (kHz) m Figure 4 | Nonlinear cavity damping. (a) The total (tot) mechanical linewidth as a function of cavity photon number. Red circles: V g ¼6.5 V, g 0 /2pC0.7 MHz. Black circles: V g ¼9.5 V yielding higher g 0 /2pC1.0 MHz. Squares: smaller g 0 /2pC0.5 MHz obtained with V g ¼4.6 V. The solid lines are expectations from the basic linear model, that is gtot m¼gmþgopt m.The dashed lines include the qubit-induced extra linewidth, gtot m¼gmþgopt mþgf m.(b) Zoom-in of the linear regime of a. In all the data, DC1, F/F 0 C0.39. NATURE COMMUNICATIONS | DOI: 10.1038/ncomms7981 ARTICLE NATURE COMMUNICATIONS | 6:6981 | DOI: 10.1038/ncomms7981 | www.nature.com/naturecommunications 5 &2015 Macmillan Publishers Limited. All rights reserved. 7. Brennecke, F., Ritter, S., Donner, T. & Esslinger, T. Cavity optomechanics with a Bose-Einstein condensate. Science 322, 235–238 (2008). 8. Purdy, T. P. et al. Tunable cavity optomechanics with ultracold atoms. Phys. Rev. Lett. 105, 133602 (2010). 9. Rabl, P. Photon blockade effect in optomechanical systems. Phys. Rev. Lett. 107, 063601 (2011). 10. Nunnenkamp, A., Børkje, K. & Girvin, S. M. Single-photon optomechanics. Phys. Rev. Lett. 107, 063602 (2011). 11. Qian, J., Clerk, A. A., Hammerer, K. & Marquardt, F. Quantum signatures of the optomechanical instability. Phys. Rev. Lett. 109, 253601 (2012). 12. Xuereb, A., Genes, C. & Dantan, A. Strong coupling and long-range collective interactions in optomechanical arrays. Phys. Rev. Lett. 109, 223601 (2012). 13. Ludwig, M., Safavi-Naeini, A. H., Painter, O. & Marquardt, F. Enhanced quantum nonlinearities in a two-mode optomechanical system. Phys. Rev. Lett. 109, 063601 (2012). 14. Lemonde, M.-A., Didier, N. & Clerk, A. A. Nonlinear interaction effects in a strongly driven optomechanical cavity. Phys. Rev. Lett. 111, 053602 (2013). 15. Børkje, K., Nunnenkamp, A., Teufel, J. D. & Girvin, S. M. Signatures of nonlinear cavity optomechanics in the weak coupling regime. Phys. Rev. Lett. 111, 053603 (2013). 16. Liao, J.-Q. & Nori, F. Photon blockade in quadratically coupled optomechanical systems. Phys. Rev. A 88, 023853 (2013). 17. Heikkila ¨, T. T., Massel, F., Tuorila, J., Khan, R. & Sillanpa ¨a ¨, M. A. Enhancing optomechanical coupling via the Josephson effect. Phys. Rev. Lett. 112, 203603 (2014). 18. Etaki, S. et al. Motion detection of a micromechanical resonator embedded in a d.c. SQUID. Nat. Phys. 4, 785–788 (2008). 19. LaHaye, M. D., Suh, J., Echternach, P. M., Schwab, K. C. & Roukes, M. L. Nanomechanical measurements of a superconducting qubit. Nature 459, 960–964 (2009). 20. O’Connell, A. D. et al. Quantum ground state and single-phonon control of a mechanical resonator. Nature 464, 697–703 (2010). 21. Pirkkalainen, J.-M. et al. Hybrid circuit cavity quantum electrodynamics with a micromechanical resonator. Nature 494, 211–215 (2013). 22. Gustafsson, M. V. et al. Propagating phonons coupled to an artificial atom. Science 346, 207–211 (2014). 23. Ramos, T., Sudhir, V., Stannigel, K., Zoller, P. & Kippenberg, T. J. Nonlinear quantum optomechanics via individual intrinsic two-level defects. Phys. Rev. Lett. 110, 193602 (2013). 24. Pflanzer, A. C., Romero-Isart, O. & Cirac, J. I. Optomechanics assisted by a qubit: from dissipative state preparation to many-partite systems. Phys. Rev. A 88, 033804 (2013). 25. Jo ¨ckel, A. et al. Sympathetic cooling of a membrane oscillator in a hybrid mechanical-atomic system. Nature Nanotechnology 10, 55–59 (2014). 26. Sillanpa ¨a ¨, M., Roschier, L. & Hakonen, P. The inductive single-electron transistor (L-SET). Phys. Rev. Lett. 93, 066805 (2004). 27. Niemczyk, T. et al. Circuit quantum electrodynamics in the ultrastrongcoupling regime. Nat. Phys. 6, 772–776 (2010). 28. Schreier, J. A. et al. Suppressing charge noise decoherence in superconducting charge qubits. Phys. Rev. B 77, 180502 (2008). 29. Naaman, O. & Aumentado, J. Time-domain measurements of quasiparticle tunneling rates in a single-cooper-pair transistor. Phys. Rev. B 73, 172504 (2006). 30. Shaw, M. D., Lutchyn, R. M., Delsing, P. & Echternach, P. M. Kinetics of nonequilibrium quasiparticle tunneling in superconducting charge qubits. Phys. Rev. B 78, 024503 (2008). 31. Marquardt, F., Chen, J. P., Clerk, A. A. & Girvin, S. M. Quantum theory of cavity-assisted sideband cooling of mechanical motion. Phys. Rev. Lett. 99, 093902 (2007). 32. van Wezel, J. & Oosterkamp, T. H. A nanoscale experiment measuring gravity’s role in breaking the unitarity of quantum dynamics. Proc. R. Soc. A 468, 35–56 (2012). 33. Pikovski, I., Vanner, M. R., Aspelmeyer, M., Kim, M. S. & Brukner, C. Probing Planck-scale physics with quantum optics. Nat. Phys. 8, 393–397 (2012). 34. Carlisle, A., Kwon, H., Jeong, H., Ferraro, A. & Paternostro, M. On the limitations of an optomechanical route to quantum macroscopicity of superposition states. preprint at http://arxiv.org/abs/arXiv:1411.7980 (2014). Acknowledgements This work was supported by the Academy of Finland (contract 250280, CoE LTQ) and by the European Research Council (240387-NEMSQED, 240362-Heattronics, 615755CAVITYQPD), and by FP7 grant 323924 iQUOEMS. The work benefited from the facilities at the Micronova Nanofabrication Center and at the Low Temperature Laboratory infrastructure. Author contributions J.-M.P. and S.U.C. designed and fabricated the devices. J.-M.P. carried out the measurements and a part of the data analysis. F.M. and J.T. and T.T.H. provided the theoretical results. M.A.S. conceived the work and designed the experimental setup with P.J.H. Additional information Supplementary Information accompanies this paper at http://www.nature.com/ naturecommunications Competing financial interests: The authors declare no competing financial interests. Reprints and permission information is available at http://npg.nature.com/ reprintsandpermissions/ How to cite this article: Pirkkalainen, J.-M. et al. Cavity optomechanics mediated by a quantum two-level system. Nat. Commun. 6:6981 doi: 10.1038/ncomms7981 (2015). This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/ ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms7981 6NATURE COMMUNICATIONS | 6:6981 | DOI: 10.1038/ncomms7981 | www.nature.com/naturecommunications &2015 Macmillan Publishers Limited. All rights reserved.