substantial progress in fabricating high-@xmath0 mechanical resonators with high frequencies , as well as recent success in cooling them close to the motional ground state @xcite , inaugurate a new research field of manifold fundamental interest @xcite .
the regime of very low temperature , where a quantum mechanical description predicts only few quanta of mechanical motion , promises potential insight into some fundamental questions of decoherence , as well as various technical applications that make use the of the expected quantum behavior @xcite .
fundamental questions , concerning the border between the classical ( macroscopic ) and the quantum ( microscopic ) worlds @xcite , trigger a natural interest in preparing quantum states of `` as large as possible '' objects and demonstrating their distinct quantum behavior by appropriate measurements . see @xcite for a recent review .    regarding this major goal , it is important to stress that the dynamics of a purely harmonic quantum system is analogous to its classical dynamics , in the sense that expectation values of canonical observables follow the classical equations of motion @xcite .
therefore , it is a common approach to introduce nonlinearities in a quantum system in order to detect quantum behavior . while there may be the possibility to achieve the strong optomechanical coupling regime @xcite and make use of the nonlinear nature of the standard optomechanical coupling , or to couple to a nonlinear ancilla system @xcite , we propose here a different approach : the use of an optoelectromechanical system featuring a tunable _ mechanical _ nonlinearity per phonon .
the latter originates from the intrinsic geometric nonlinearity of elastic systems @xcite and its amount per motional quanta is enhanced with the help of electrostatic fields .
this has the advantage that the linear optomechanical coupling is preserved as a control channel providing techniques as , for example , the sideband driving technique used in @xcite .
the regime of large mechanical nonlinearity then enables new means to control the mechanical motion at the quantum level , if combined with the coupling to a high - q optical cavity mode , as well as the application of suitable gradient fields .
the intrinsic anharmonicity in the mechanical motion of micro- and nanomechanical resonators is usually small and therefore only relevant in the regime of large oscillation amplitudes . in order to render the anharmonicity relevant for displacements at the scale of the quantum mechanical zero point motion
, we propose to use electrostatic gradient forces to enhance the latter @xcite .
these forces result from the dielectric properties of the resonator material when an inhomogeneous external electric field is applied .
they can be used to effectively reduce the resonator s stiffness and therefore its resonance frequencies . in turn , this has the effect that the zero point deflection is enhanced up to an extent that the nonlinear contribution becomes important .    using this technique ,
the nonlinearity per phonon can be made large enough , that distinct transitions in the mechanical spectrum can be resonantly addressed while interacting with other quantum systems .
examples are the selective sideband driving of transitions in the mechanical spectrum @xcite , or the resonant exchange of excitations within an array of nanoresonators via a common cavity mode @xcite .    in this paper , we explicitly derive the fundamental mode properties of a nonlinear mechanical resonator , subject to the aforementioned gradient forces , to obtain a suitable model for the mechanical degree of freedom .
we then derive the specifics of the optomechanical coupling to a high finesse cavity and analyze possible source of damping and decoherence in detail
. we also summarize different control schemes , associated with suitable laser drives for the cavity and gradient fields from the tip electrodes .
two applications of these control mechanisms have been proposed previously @xcite , considering state of the art experimental components .
we finally discuss a readout scheme for the nonlinear mechanical resonators considered .
the remainder of the paper is organized as follows . in section [ sec : elasticity ] we describe the nonlinear dynamics of thin rods starting from elasticity theory and derive the resulting fundamental mode hamiltonian . in section [ sec : softening ] we describe our approach to enhance the intrinsic mechanical nonlinearity by mode softening with gradient forces .
then we introduce the optomechanical model with several laser driven cavity modes in section [ sec : optomech ] . after introducing a possible setup for an implementation in section [ sec : setup ] ,
we quantitatively discuss its central physical properties , namely the optomechanical coupling mechanism in section [ sec : coupling ] , as well as potential setup specific losses in section [ sec : losses ] . in section [ sec : controlmech ] we review some mechanisms that can be used to control the mechanical motion at the level of single phonons and which have been applied in previous works @xcite .
finally , in section [ sec : measurement ] we introduce methods to obtain information on the mechanical state from the output spectrum of a probe laser .
the harmonic description of the transverse motion of thin rods is based on only considering the bending energy for small deflections @xcite .
we consider here thin rods , which means that the cross - sectional dimensions , such as the radius for circular cross sections or width and depth for rectangular cross sections , are much smaller than the length @xmath1 .
we also consider the rod to be homogeneous along the longitudinal axis , here parametrized by @xmath2 , with constant mass line density @xmath3 and use thin rod elasticity theory .
the planar deflection in the transverse direction is described by a field @xmath4 and we consider a bridge geometry where the end points at @xmath5 and @xmath6 are fixed , i.e.  @xmath7 and @xmath8 .
the lagrangian within this approximation reads @xmath9\ , , \label{lagrangian}\ ] ] with a kinetic part as well as the bending energy @xmath10=\frac{1}{2}\int\limits \mathcal{f}\tilde{\kappa}^2(y'')^2{\rm d}x\ , .
\label{bending_energy}\ ] ] here , @xmath11 is the linear modulus of the rod given by the young modulus @xmath12 of the material times the cross - section area @xmath13 , and @xmath14 is the ratio between the bending and compressional rigidities and depends on the cross - sectional geometry , where @xmath15 is the in - plane coordinate within the cross - section that is directed along the deflection with origin at the neutral line , see figure [ fig : bending_fig1 ] . for a rectangular cross - section of thickness @xmath16 , @xmath17 , whereas for a circular cross - section with radius @xmath18 , @xmath19 . for a cylindrical shell like a nanotube one finds @xmath20 .
the energy ( [ bending_energy ] ) results from the fact that for small curvature @xmath21 , the local strain inside the rod is linear with respect to the distance @xmath15 to the neutral surface ( cf .
1 ) while the free energy density is quadratic with respect to the strain .     depends on the transverse coordinate @xmath15 and on the local curvature @xmath21 and determines the energy density @xmath22 . ]
the lagrangian ( [ lagrangian ] ) leads to the equation of motion @xmath23 as this equation is linear in @xmath24 and its derivatives , it leads to harmonic dynamics characterized by the following eigenmodes @xmath25}\,,\nonumber\end{aligned}\ ] ] with frequencies @xmath26 , where @xmath27 is the phase velocity of compressional phonons along the rod .
the @xmath28 are the roots of the transcendental equation @xmath29 , with @xmath30 the smallest one .
the @xmath31 are normalization constants chosen such that @xmath32 .
we choose this normalization so that the coefficients in a mode expansion represent the maximum amplitudes of the deflection associated to each mode .
introducing now the canonical momentum @xmath33 , as well as the expansion of the field into the modes @xmath34 yields the hamilton function of a harmonic oscillator for each mode @xmath35 with the deflection @xmath36 and mode momentum @xmath37 for the @xmath38-th mode , as well as the effective mode masses @xmath39 .
corrections to this harmonic description that lead to nonlinearities originate from a stretching effect that occurs due to the deflection if the end points of the rod are fixed @xcite .
the resulting strain leads to an additional energy @xmath40 where the stretched length is @xmath41 with @xmath1 being the zero deflection length . including this streching energy leads to a nonlinear extension of the hamiltonian , which after inserting the modes given in eq .
leads to @xmath42 where @xmath43 .
we quantize this model by introducing bosonic mode operators @xmath44 and @xmath45 , given by @xmath46 where we introduced the zero point motion amplitudes @xmath47 for each mode .
this leads to the description @xmath48 with nonlinearity @xmath49 as in a rigorous elasticity treatment this description arises from an adiabatic elimination of the stretching modes , the indices in eq .
should run up to an @xmath50 corresponding to an `` ultraviolet '' cutoff @xmath51 .
the terms involving higher order modes induce small frequency shifts and nonlinear mode coupling .
however , the later is found to be negligible for the parameters considered ( see appendix [ appendix b ] ) and the shift of the fundamental mode can be ignored given the additional tunable electrostatic contribution ( see section [ sec : softening ] ) .
therefore we can restrict our description to the fundamental mode with @xmath52 and drop this label to get the usual hamiltonian for the duffing oscillator @xmath53 where we have introduced the fundamental frequency @xmath54 and the effective mass of the fundamental mode @xmath55 .
the anharmonicity is given by @xmath56 in terms of phonon creation and annihilation operators @xmath57 and @xmath58 this hamiltonian reads , @xmath59 with the nonlinearity parameter @xmath60 . here
, the frequencies @xmath61 and @xmath62 refer to fundamental mode properties that result only from the intrinsic elastic forces in the absence of any externally applied forces on the rod .
as we will describe in the next section , external forces can be used to tune the resonance frequency , @xmath63 , which will in turn change the zero point motion amplitude , @xmath64 , and hence the nonlinearity , @xmath65 .
in this section we describe how electric fields generated by tip electrodes that are placed near the center of the doubly clamped nanobeam , see figure [ fig : bending_fig2 ] , can be employed to control its dynamical properties . in particular ,
inhomogeneous gradient fields can be used to enhance the nonlinearity per phonon .
the later scales as @xmath66 and can thus be enhanced by lowering the harmonic oscillation frequency @xmath61 .
one way to change the mode frequency that has been discussed previously , is to add an additional external force along the rod axis , that causes compressive or tensile strain @xcite .
an alternative approach , that promises better control but yields the same potential for the fundamental mode , is to use a static electric field that is strongly inhomogeneous in the direction of deflection .
if the rod material shows suitable dielectric properties , this produces an additional , inverted square potential with respect to the deflection , see figure [ fig : bending_fig2 ] .
here we consider tip electrodes with suitable applied voltages at each side of the nanoresonator that generate an electrostatic field . in our calculation
we model the electrodes by point charges @xmath67 and @xmath68 , which is valid given that the relevant tip radii are much smaller than the gap between the electrodes ( cf .  section [ sec : setup ] ) . the electrostatic energy associated to the dielectric per unit length along the rod can be described by @xmath69\,,\ ] ] where @xmath70 are the co - ordinates along the resonator axis and the direction of its deflection .
@xmath71 are external field components parallel and perpendicular to the beam axis and @xmath72 the respective screened polarizabilities
. we can expand @xmath73 to second order in the displacement @xmath24 and get an additional contribution to the hamiltonian of the nanobeam that reads @xmath74{\rm d}x\,,\nonumber\end{aligned}\ ] ] where we dropped the displacement independent constant @xmath75 which is irrelevant for the dynamics . inserting the modes defined in equation we get @xmath76 with @xmath77 the contributions @xmath78 will cause weak interactions between the modes , as @xmath78 is not a diagonal matrix . by diagonalizing the hamiltonian of the rod in the presence of electrostatic fields
, one can find new normal modes .
for the parameters considered in this work , the induced corrections to the mode shapes are however found to be negligibly small .
we thus focus on the fundamental mode contributions @xmath79 and @xmath80 , where @xmath81 .    .
the field profile ( red ) leads to an inverted parabola for the dielectric potential ( blue ) around the equilibrium position . ]
one may also consider time dependent electric fields and it is convenient to separate between static and time dependent contributions , @xmath82 whereas the time dependent contributions @xmath83 and @xmath84 generate drives applied to the nanobeam and will be considered in section [ sec : controlmech ] , we now focus on the constant contributions @xmath85 and @xmath86 , which can be employed to tune the nanobeam .
the electrostatic force @xmath85 causes a static deflection of the nanobeam and shifts its equilibrium position .
however if the rod interacts with the photon fields of a nearby cavity , those fields will also cause a deflecting force . for convenience
we choose @xmath85 such that these two forces compensate and the equilibrium position remains unshifted ( see section [ sec : optomech ] ) .
the electrostatic potential associated to @xmath86 in turn is an inverted harmonic potential that lowers the harmonic oscillation frequency of the nanobeam .
therefore , we consider @xmath87 as the `` tuned '' mechanical hamiltonian with a reduced frequency @xmath88 . in a phononic description
this hamiltonian reads @xmath89 where the nonlinearity per phonon @xmath90 is now increased by a factor @xmath91 , cf .
as an example , for maximum applied fields at the tube @xmath92vm@xmath93 and @xmath94vm@xmath93 and a gap between the electrodes of size @xmath95 nm ( cf .
3 ) , typical parameters discussed in section [ sec : setup ] , yield @xmath96 allowing to boost @xmath97 by at least an order of magnitude .    for further calculations , it is convenient to express all observables in the energy eigenbasis of the hamiltonian @xmath98 , so that @xmath99 where the energy eigenstates @xmath100 and energy levels @xmath101 , as well as the displacement matrix elements @xmath102 have to be determined numerically . for small nonlinearites
@xmath103 analytical expression can be obtained as the hamiltonian is approximately diagonal in fock basis since one may apply a rotating wave approximation in the nonlinear part @xmath104 where @xmath105 , @xmath106 and the eigen - energies are given by @xmath107 .
the time dependent contributions @xmath108 and @xmath109 in turn can be used to drive or temporarily detune the resonator .
this has been used in @xcite to perform local gate operations on nanobeams acting as qubits .
in the following section we will now discuss an optomechanical interaction between the nano resonator described by the hamiltonian in equation ( [ hm1 ] ) and the resonance modes of a high finesse cavity .
we consider a typical optomechanical setup with a micro toroid cavity coupled to several nanomechanical resonators , where the displacement of the latter modifies the frequencies of the cavity modes .
the cavity modes with frequencies @xmath110 that contribute to the dynamics are described by photon creation and annihilation operators @xmath111 and @xmath112 .
they are each driven by a classical laser of input power @xmath113 and frequency @xmath114 . the coupling strength between cavity mode @xmath115 and nanobeam @xmath116
is given by @xmath117 , where @xmath118 is the optical frequency shift per deflection @xmath119 .
thus , in a frame rotating with the laser field modes , the hamiltonian describing the coupled system reads ( @xmath120 ) @xmath121+\sum\limits_jh_{\mathrm{m},j}\nonumber\\ & + \sum\limits_{ij}g_{0,ij}{a_i}^{\dag}a_i{\ensuremath{\mathcal{x}}}_j\ , .
\label{hamiltonian}\end{aligned}\ ] ] here we have introduced the laser detunings denoted by @xmath122 and drive amplitudes @xmath123 , with @xmath124 being the cavity decay rates into the associated outgoing electromagnetic modes .
both the field inside the cavity and the mechanical motion are subject to damping , which in the regime of weak optomechanical coupling , small nonlinearity and low mechanical occupation is well described by a master equation with lindblad form damping terms . with the decay rates for the cavity modes @xmath125 and the mechanical damping rates @xmath126 , this master equation reads @xmath127&+\sum\limits_i\frac{\kappa_i}{2}{\ensuremath{\mathcal{d}}}(a_i)\rho\,+\\   \nonumber&+\sum\limits_j\frac{\gamma_j}{2}\left[(\overline n_j+1){\ensuremath{\mathcal{d}}}(b_j)\rho+\overline n_j { \ensuremath{\mathcal{d}}}({\ensuremath{b^{\dag}}}_j)\rho\right]\ , .
\end{aligned}\ ] ] here we introduced the lindblad form dissipator @xmath128 , as well as the bose occupation number @xmath129^{-1}$ ] of the phonon bath mode with frequency @xmath130 at temperature @xmath131 . a more precise treatment of the mechanical damping accounts for the fact that , due to the mechanical nonlinearity , there is more than one bath mode resonantly coupling to the mechanical mode
however , for small nonlinearity , equation ( [ mg ] ) proofs to be sufficiently accurate .    as usual , we expand the cavity field operators around their steady state values and adopt a shifted representation @xmath132 , with @xmath133 , in which the master equation has the same form as in equation ( [ mg ] ) but with the shifted system hamiltonian @xmath134 with @xmath135 and where we have dropped the nonlinear terms @xmath136 in the coupling , which is valid for @xmath137 .
we have also assumed that the static electric fields for each beam are chosen such that @xmath138 , so that their equilibrium positions are undeflected .
we now turn to discuss a possible experimental setup that would allow to explore the physics described by the model presented in equation .
d ) schematic cross - sectional view of the nems - chip and the cavity rim ( yellow ) . the nems - chip is mounted in the evanescent cavity field at a distance @xmath16 to the cavity rim surface .
the nanoresonator is displaced from the closest point to allow for an optomechanical coupling that is linear in the resonator s deflection.,title="fig : " ] +    for the experimental realization of the model in equation , we envisage a setup as shown in fig .  [
fig : setup ] , comprising a nems chip containing the nanobeam resonators and a high finesse toroidal microcavity @xcite .
each nanobeam resonator consists of a suspended single - walled carbon nanotube ( cnt ) with radius @xmath18 ( @xmath140ms@xmath93 , @xmath141kgm@xmath142 ) @xcite , and has electrodes in its vicinity that generate the electric fields for controlling and driving it .
furthermore , the nanoresonator interacts with the evanescent field of the microcavity via optical gradient forces ( see section [ sec : coupling ] ) . given the state of the art
, a cnt is a favorable system to implement the proposed optoelectromechanical scheme @xcite . in particular , the intrinsic nonlinearity per phonon scales as @xmath143 , which favors small transverse dimensions and masses .
additionally , carbon nanotubes show ultra - low dissipation @xcite , that is expected to decrease further in the regime of small amplitudes @xcite .
it is convenient to use ( @xmath144,@xmath145 ) cnts with radius @xmath146 , as these tubes show relatively large static polarizabilities with @xmath147 , @xmath148 @xcite , where @xmath149 ( @xmath150 ) is the polarizability parallel ( perpendicular ) to the tube axis . to obtain nonlinearities that are large enough
, it is convenient to use cnt lengths below @xmath151 for @xmath152 we have @xmath153 and @xmath154 .
the cavity is a silica ( @xmath155 ) microtoroid with resonant wavelengths @xmath156 , circumference @xmath157 , finesse @xmath158 @xcite and mode volume @xmath159 .
these parameters correspond to @xmath160 , @xmath161 and a decay length of the evanescent field @xmath162 [ cf .
eq .   in section [ sec : coupling ] ] .
the nems chip is placed at a distance @xmath163 from the cavity rim .
the cnt - resonator is displaced from the closest point to allow for a linear coupling , as the resonator moves in the plane of the chip s surface ( cf .
section [ sec : coupling ] ) , and the electrodes are aligned parallel to the rim of the cavity to minimize additional cavity losses they might induce , see section [ sec : losses ] .    in the following two sections we analyze important practical aspects of the envisioned implementation with a carbon nanotube coupled to a toroidal microcavity in more detail .
thus , readers who are only interested in the results of the discussed mechanism may directly turn to section [ sec : controlmech ] .
in this section we derive an estimate for the coupling between the mechanical displacement and the cavity field .
the coupling arises as the energy of the dielectric oscillator in the evanescent electric field depends on the field strength at the location of the oscillator . as the evanescent field decays with distance to the cavity rim , altering the oscillator - cavity distance by displacing the oscillator results in a change of energy .
thus , the coupling part of the hamiltonian is given by @xmath164 where the polarization is @xmath165 with @xmath166 the screened polarizability tensor and @xmath167 the cavity field .
the integration is taken over the oscillator volume @xmath168 and without loss of generality we consider a single cavity mode with resonant frequency @xmath169 .
given the dimensions of the toroidal microcavity , its torus can be locally modeled as a cylindrical waveguide of radius @xmath170 note that this differs from the definition of @xmath170 used in ref .   by a factor of @xmath171 [ cf .  below eq .
we introduce cylindrical coordinates @xmath172 with the @xmath173-direction along the waveguide axis ( cf .
[ fig : setup]c ) and consider te@xmath174 modes , as a transverse electric field is advantageous to suppress loss mechanisms that are discussed in section [ sec : losses ] .
the corresponding transverse fields are given by @xcite @xmath175 for the field inside the waveguide , @xmath176 , and @xmath177 outside the waveguide , @xmath178 .
the axial field reads @xmath179 with the modified bessel function @xmath180 and the bessel function of the first kind @xmath181 .
here @xmath182 , @xmath183 is the wavevector component parallel to the waveguide axis and @xmath184 and @xmath185 are the transverse wave vectors outside and inside the waveguide , respectively . henceforth , we assume a refractive index such that @xmath186 and a frequency @xmath169 well above cutoff , i.e.  @xmath187 where @xmath188 and @xmath189 is the first zero of @xmath190 .
these assumptions imply that @xmath191 , @xmath192 , @xmath193 , @xmath194 , and @xmath195 . within these approximations ,
the ratio of the axial magnetic field at @xmath196 to its value at the origin is given by @xmath197 , and the evanescent field can be written as @xmath198 this field will later be used to estimate losses induced by the electrodes . in order to determine the optomechanical coupling strength , we write the electric field in its quantized form @xmath199 with photon creation ( annihilation ) operators @xmath200 and where @xmath201 is the corresponding normalized eigenmode , satisfying @xmath202 we consider the following definition for the mode volume @xmath203 from eqs .   and and properties of the bessel functions
, we obtain that the maximum electric field inside the cavity is given by @xmath204 where @xmath205 is the first positive root of @xmath206 .
thus , neglecting the small contributions of the evanescent part to the integrations in eqs .   and which are higher order in @xmath207 , we obtain from eqs .  - @xmath208 for @xmath178 . here , @xmath209 and @xmath210 denotes the ratio of the electric field at the waveguide s surface to the maximum field @xmath211 .
the zero point motion @xmath64 of the oscillator is small compared to the decay length @xmath212 of the evanescent field and the same holds for the transverse dimensions of the oscillator .
thus , we can linearize @xmath213 around the equilibrium position of the nanoresonator and assume that the electric field is constant everywhere inside the resonator volume @xmath168 . subsequently , by comparing the hamiltonians and , and using eqs .   and
, we find for the opto - mechanical coupling rate , @xmath214 here , we have neglected the contribution of the perpendicular polarizability since , for carbon nanotubes , the perpendicular polarizability is typically one order of magnitude smaller than the parallel one , and again used @xmath215 so that only the derivative of the exponential factor is relevant . the geometry of the setup we consider here is illustrated in figure [ fig : setup ] and the dependence on the alignment and positioning of the nanotube is accounted for in the correction factor @xmath216 . for a te@xmath174 mode of the cavity field , the electric field
is directed along @xmath217 , i.e. tangential to the cavity rim , while the nanotube is aligned along @xmath218 .
in addition , the deflection of the nanotube in the direction @xmath219 is not aligned with the interaction - energy gradient , which is approximately along @xmath220 .
finally , if @xmath16 is the distance of the chip to the cavity rim , the actual distance of the oscillator to the rim is @xmath221 .
taking into account these various issues , we find for the correction factor @xmath222 this is maximized for @xmath223 and @xmath224 , where we consider the leading order in the small parameter @xmath225 .
these optimal angles yield @xmath226 , resulting in @xmath227 for the parameters introduced in section [ sec : setup ] ( i.e.  @xmath155 and @xmath228 ) .
finally , for those typical values we obtain @xmath229 .
in addition to the well known loss mechanisms of photon losses from the cavity and intrinsic phonon losses of the mechanical resonator @xcite , there can be further sources of loss in our setup due to the presence of the tip electrodes . in this section
we show that these additional loss mechanisms are negligible for the parameters we envision .      exploiting the nonlinearity of the nanoresonators to control their dynamics
demands a high-@xmath0 optical cavity as the linewidth @xmath230 needs to be at most comparable to the mechanical nonlinearity @xmath97 .
using conventional metallic electrodes to generate the inhomogeneous control fields can potentially increase the cavity losses . to achieve the necessary low losses it is crucial to have deep subwavelength transverse dimensions for the electrodes .
we now give an estimate of the photon losses that may be induced by such electrodes and show that these are negligible .
to do so , we model the electrodes as metallic cylinders and assume that their radius @xmath231 is much smaller than the decay length of the evanescent cavity field , @xmath232 .
we assume the electrodes to be parallel to the waveguide representing the cavity rim , with a small misalignment angle @xmath233 . for the relevant te@xmath234 modes , the losses result solely from the misalignment , as in the small radius regime considered they arise only from the field along the electrode axis which vanishes for @xmath235 .
the resulting finesse @xmath236 for the cavity can be obtained from the finesses @xmath237 associated to different decay channels , which can be assumed to be independent such that @xmath238 for each loss channel , the finesse @xmath237 can be determined from the ratio of the power circulating in the cavity @xmath239 to the fraction of power that it lost through the respective channel @xmath113 , using that the time - averaged stored energy and free spectral range are given , respectively , by @xmath240 and @xmath241 .
thus , we arrive at @xmath242 the loss channels that will be considered subsequently are ( 1 ) scattering by the `` bulk '' of the electrodes modeled as a single metallic cylinder , ( 2 ) scattering by the `` gap '' between the electrodes and ( 3 ) absorption .
the independence between the contributions ( 1 ) and ( 2 ) assumed in eq .   amounts to neglecting the interference between them which is permissible when estimating an upper bound .
[ [ incident - field - and - circulating - power ] ] incident field and circulating power + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    for our calculations , we introduce new cylindrical co - ordinates @xmath243 for the electrode with the @xmath244-direction along its axis .
we first express the cavity field that is incident on the electrode in these primed coordinates .
the result will later be used to determine the scattered and absorbed fractions of the incident power .    to this end
, we consider the projection of the electric field determined in section [ sec : coupling ] onto the electrode s axis , see fig .
[ fig : coordinates_electrode_wg ] , @xmath245 which completely determines the losses in the small radius regime considered .
the origin of the primed axis lies at @xmath246 and the relevant unit vectors are given by @xmath247 for points on the @xmath244-axis we have @xmath248 for our further calculations it is convenient to express the incident field via its fourier transform @xmath249 . using this and equations ( [ coordinate_rel1])-([coordinate_rel2 ] ) in equations and we arrive at @xmath250 where we have substituted @xmath251 and @xmath252 .
one can find an approximation for the integral @xmath253 by applying the method of steepest descents , using @xmath254 and @xmath255 , which for @xmath256 yields @xmath257     and waveguide coordinates @xmath258 .
the misalignment is determined by a small angle @xmath233.,scaledwidth=40.0% ]    finally , we determine the circulating power to leading order in @xmath207 , which is given by @xmath259 where we have used eqs .   and , that @xmath260 , and the time - averaged poynting vector @xmath261 with the vacuum permeability @xmath262 .
[ [ scattering - losses ] ] scattering losses + + + + + + + + + + + + + + + + +    here , we model the electrodes as a single metallic cylinder which for simplicity is assumed to be perfectly conducting since this maximizes the scattering and , thus , provides an estimate of an upper bound to the corresponding losses that is independent of material properties naturally for the transparent electrode scenario considered below , in section [ electrode losses].d , these losses would be substantially smaller than this upper bound .
we expand the scattered field into solutions of the wave equation in cylindrical coordinates for the electrode .
as the radius of the electrode @xmath231 is much smaller than the wavelength of the incident field @xmath263 , all contributions to the scattered power are suppressed at least like @xmath264 , except for s - wave scattering of tm modes , for which the suppression is only logarithmic .
this can be understood in terms of the taylor expansions of the corresponding cylindrical harmonics and the incident field . in turn , to determine the tm s - wave scattering to leading order in @xmath265 , the incident field can be assumed to be constant for a given cross section and determined by the field at the electrode s
center @xmath266 .
we neglect multiple scattering between the waveguide and the electrode and ignore the dielectric substrate of the latter .
thus , the scattered field @xmath267 is determined from the homogeneous boundary condition at the surface of the electrode @xmath268 , and an outgoing - wave boundary condition at infinity for @xmath267 here @xmath269 is the evanescent contribution .
the transverse fields of an outgoing tm solution with @xmath244-dependence @xmath270 are given by @xmath271 with @xmath272 .
for @xmath273 the solution is evanescent and does not contribute to the scattered power . for s - wave scattering @xmath274 , where @xmath275 , and to leading order in @xmath265 , the scattered field can be written as @xmath276 where @xmath277 is the fourier transform of the incident field @xmath266 and we have used the approximation @xmath278 for @xmath279 .
we calculate the scattered power by integrating the energy flux across a cylinder coaxial with the electrode with radius @xmath280 .
thus , from eqs .
- we obtain for the scattered power @xmath281}}_{g } \label{scatpower}\end{aligned}\ ] ] where we have used @xmath282 and @xmath283 $ ] for @xmath284 .
the integral @xmath285 can be estimated by performing the substitution @xmath286/\ln(2kr)$ ] and considering @xmath287 , which yields @xmath288 .
hence , from eqs .  , - , and , and using @xmath289 and @xmath290 for @xmath256 , and @xmath291 we arrive at @xmath292}}{n_{\mathrm{c}}\left(k{\ensuremath{a_{\mathrm{c}}}}\sqrt{n_{\mathrm{c}}^2 - 1}\right)^3|\ln2kr'|}\,.\ ] ] if we now use that for relevant parameters @xmath293 , @xmath294 and @xmath295 , we find a lower bound for the finesse associated to scattering losses , @xmath296 which is independent of the ratio @xmath297 .
thus , for @xmath298 , @xmath299 , @xmath155 and @xmath300 , we find @xmath301 .
hence losses due to scattering are clearly sufficiently suppressed for electrodes of subwavelength radius that are approximately aligned with the cavity rim .
[ [ gap - contribution - to - scattering - losses ] ] gap contribution to scattering losses + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    so far we have considered scattering from one single continuous electrode .
actually , our setup comprises two electrodes , see figure [ fig : setup ] , separated by a gap in which the nanomechanical oscillator is positioned .
we now estimate an upper bound to the additional losses that this gap may induce . here
we assume as before @xmath302 so that the gap is much smaller than the distance over which substantial currents are induced in the electrodes and one may consider @xmath235 for this estimate . denoting the direction joining the electrodes by @xmath303
, we focus on the relevant regime @xmath304 so that one can assume that only the incident field component along @xmath303 is screened by them .
we model the gap of size @xmath305 as a perfectly conducting sphere with radius @xmath306 subject to an external field along @xmath303 determined by @xmath307 .
this should provide an upper bound for the magnitude of the total induced dipole @xmath308 , which in turn determines the leading contribution in @xmath309 to the scattered power .
thus we find @xmath310 @xcite for optimal placement , which yields for the scattered power @xmath311 ^ 2    { \ensuremath{\mathrm{e}^}}{-2\kappa_{\bot}d}\cos^2\theta'_*\,.\ ] ] here we have used eq .  , @xmath312 and @xmath313 . along the same lines as before , using eqs .   and and @xmath223 , the associated finesse reads @xmath314 which yields @xmath315 for the parameters discussed in section [ sec : setup ] and @xmath95 nm . thus , additional scattering losses due to the gap between the electrodes are also negligible .    [ [ absorption - losses ] ] absorption losses + + + + + + + + + + + + + + + + +    we assume here transparent electrodes afforded by cylindrical shells with 2d conductivity @xmath316 that absorb the power @xmath317 since @xmath318 induces a current @xmath319 in each electrode here , as in [ electrode losses].b .
, we neglect the small gap . from eqs .   and
- we obtain for the absorbed power @xmath320 ^ 2 { \ensuremath{\mathrm{e}^}}{-2\kappa_{\bot}d}|b_{z}(0)|^2\,\times\\      & \times\underbrace{\int \limits_{-\infty}^{\infty}{\rm d}x\frac{{\ensuremath{\mathrm{e}^}}{-2\kappa_{\bot}(d+{\ensuremath{a_{\mathrm{c } } } } ) ( \sqrt{1+x^2}-1)}}{(1+x^2)^{3/2}}}_{j}{\nonumber}\end{aligned}\ ] ] where we have again substituted @xmath251 . using the method of steepest descents we estimate @xmath321 for @xmath322 . thus using @xmath323 and eqs .  , and , yields for the finesse associated to absorption @xmath324 ^ 2\,.\ ] ] we consider now the specific case where the transparent electrodes are provided by a pair of nanotubes . the latter exhibit
a maximum in the conductivity @xmath325 @xcite . by assuming an off - resonant @xmath326 with @xmath327
, we get @xmath328 with the fine structure constant @xmath329 .
if we consider the same values as before except that now @xmath330 and @xmath331 and assuming @xmath332 we find @xmath333 . hence even though absorption losses clearly dominate over scattering losses , their effect can still be neglected for electrode radii @xmath334 and alignment angles @xmath330 .
here we give estimates of decoherence rates for the nanoresonator induced by noise in the inhomogeneous electric fields .
such noise might originate from voltage fluctuations @xmath335 due to the electrodes resistance ( johnson - nyquist noise ) or from moving charges on the chip surface ( @xmath336-noise ) .
we calculate the respective single - phonon decoherence rates @xmath337 and @xmath338 from the corresponding noise spectra @xmath339 and @xmath340 using the relation @xmath341 with @xmath342 where @xmath343 is the force fluctuation acting on the resonator . the electrostatic gradient force acting on a resonator
can be expressed by @xmath344 where we have estimated the field gradient at a distance @xmath345 from the electrode by @xmath346 and used the fact that the field mainly acts on the nanotube in a region of length @xmath347 . for a field with fluctuations associated to different independent sources @xmath348 , the force fluctuations are then given by @xmath349 thus , the resulting decoherence rates read @xmath350 where the @xmath351 are the noise spectra for the different electric field fluctuations .
[ [ johnson - nyquist - noise ] ] johnson - nyquist noise + + + + + + + + + + + + + + + + + + + + +    for johnson - nyquist noise @xcite , we have fluctuating voltages @xmath335 with @xmath352 for an ambient temperature @xmath131 and an internal resistance @xmath353 . for our setup
we find @xmath354 at @xmath355 , which is well below the relevant mechanical decoherence rate @xmath356 for relevant resistances @xmath357 .
[ [ f - noise ] ] @xmath336-noise + + + + + + + + + + + + + + + + + + + + + + +    the origin of @xmath336-noise is usually associated with surface charge fluctuations in the device . an electric field noise density
@xmath358 has been measured at @xmath359 and at a distance of @xmath360 between a charged resonator and a gold surface @xcite . for a scaling @xmath361 @xcite
this corresponds to @xmath362 for our conditions with @xmath355 and @xmath363 .
thus we expect for the associated decoherence rate @xmath364 , which is again well below the mechanical decoherence rate @xmath365 .
these results are also corroborated by recent estimates that were obtained for a related setup @xcite .
the hamiltonian of the full optomechanical system with tuned nanobeams potentially leads to complex dynamics for photons and phonons . here
, we focus on scenarios where driven cavity modes and suitable electric gradient fields are used to control the dynamics of one or several nanobeams .
we summarize the basic principles for three conceptually different schemes , namely ( i ) the selective addressing of transitions in the mechanical spectrum by cavity sideband driving , ( ii ) the coherent interaction between several nonlinear nanobeams mediated by a common driven cavity mode and ( iii ) the manipulation of a single resonator s state by time dependent gradient fields .
the first scheme represents a suitable extension of the standard sideband cooling technique @xcite to nonlinear resonators . here
, the detuning of a red ( blue ) detuned laser drive is only resonant with one specific transition @xmath366 ( @xmath367 ) in the nonlinear mechanical spectrum , see figure [ fig : schemes]a .
therefore , if the mechanical nonlinearity is resolved by the cavity linewidth , @xmath368 , appropriate laser drives can lead to highly nonclassical steady states for the mechanical motion . for a single nanobeam for example , this allows for the preparation of stationary fock states with high fidelity @xcite .
the sideband driving technique could potentially also be applied to more complicated level structures , for example the collective modes of several interacting nanobeams , and thus constitutes a versatile control mechanism .
the second scheme uses the cavity to mediate a coherent coupling between several nanobeams that all couple to the same photon mode . here
, the photon mode is driven with a large detuning to be off - resonant to any mechanical transition frequency . the coherent photon background field that builds up inside the cavity leads to an effective interaction @xmath369 between any pairs @xmath370 of nanobeams .
in order to exchange excitations via this coupling , proper resonance conditions have to be met . by tuning each of the nanobeams using their respective electrodes , interactions between desired pairs of beams
can be realized @xcite .
furthermore , due to the nonlinear spectra , it is possible to restrict the dynamics of the nanobeams to the `` qubit '' subspace built up by the states @xmath371 , c.f .
equation , for each resonator .
finally , beside the static tuning capability , see equation , the gradient fields provided by the tip electrodes can be used to perform coherent operations on any nanobeam .
this becomes most obvious if one considers the qubit subspace @xmath372 for one nanobeam . here
, a drive @xmath373 , where @xmath374 is the qubit transition frequency , implements a @xmath375 rotation , see figure [ fig : schemes]b .
a temporary shift of the qubit transition frequency @xmath376 can be achieved by a temporary @xmath84 contribution , which corresponds to a @xmath377 rotation , see figure [ fig : schemes]c .
note that the drive @xmath378 associated with the tip electrodes is a coherent drive , while the cavity sideband driving technique constitutes a stochastic drive .
together with the coherent coupling of several nanobeams , the time dependent gradient fields can for example be employed to build up a universal set of quantum gates for quantum information processing @xcite .     +
the steady state of a nanoresonator can be probed with an additional laser , weakly driving one cavity mode on resonance , i.e. with @xmath379 .
then , information about the state of the mechanical resonator can be extracted from the sideband structures in the power spectrum , which correspond to photons that have been up- or down converted during the interaction with the mechanical motion .
the intensity of the sideband peaks depends on the population of the mechanical energy eigenstates .
thus , the power spectrum only provides information about the diagonal entries @xmath380 of the density matrix describing the mechanical resonator represented in the basis formed by eigenstates of the hamiltonian @xmath98 , see eq .
( [ eq : diaghm ] ) . for a probe laser resonantly driving a cavity mode at frequency @xmath381 , the spectrum shows a sideband structure @xmath382 with lorentzian sideband peaks determined by @xmath383 ^ 2+({\ensuremath{\gamma_{\mathrm{eff}}}}^{nm}/2)^2}\,,\ ] ] with @xmath384 as given in equation ( [ gammaeff ] ) , see appendix [ appendix b ] for details . here , @xmath385 is the optomechanical coupling strength associated with the probe laser and the @xmath386 denote the mechanical transition frequencies .
the peaks appear in groups with @xmath387 ; see figure [ fig : spectrum ] .
the occupation probabilities @xmath388 for the eigenstates @xmath389 can be extracted from the peak intensities within the main sidebands with @xmath390 @xcite .        for the readout of mechanical qubits
as discussed in @xcite , a shelving technique can be used to determine whether a qubit is in state @xmath391 or in state @xmath392 . here
, a balanced cycling transition between @xmath393 using a cooling laser on @xmath394 and a coherent rf - drive with local gradient fields on those two levels causes a continuous stream of up converted photons only if the resonator is found in the state @xmath392 .
this can be detected by measuring the corresponding sideband spectrum . here
, a large enough number of photons has to be collected before external damping destroys the intermediate state , which requires @xmath395 .
we have introduced a scheme to access a novel regime of optomechanics where the motion of the nanomechanical resonator becomes anharmonic and thus allows to explore genuine quantum dynamics . in our approach inhomogeneous electrostatic fields
are applied to the nanomechanical resonator to enhance its anharmonicity per phonon until it becomes comparable to the linewidth of a high finesse optical cavity . for realistic experimental conditions ,
sufficiently large optomechanical couplings can be realized and losses induced by the tip electrodes can be suppressed to a negligible level .
furthermore populations of the energy eigenstates of such nonlinear mechanical oscillators can be extracted from the output spectrum of a probe laser .
the approach thus paves the way towards exploring nonclassical dynamics of nanomechanical oscillators at the single - phonon level .
this work is part of the emmy noether project ha 5593/1 - 1 and was supported the crc 631 , both funded by the german science foundation ( dfg ) .
iwr acknowledges support by the german science foundation ( dfg ) via the grant wi 3859/1 - 1 and nanosystems initiative munich ( nim ) .
the authors thank j.  kotthaus , t.j .
kippenberg , and a.  schliesser for enlightening discussions .
in order to estimate the strength of the mode coupling , we rewrite the nonlinearity as @xmath396\ , , \label{bracket}\ ] ] where @xmath397 is the physical mass of the rod . the term in brackets solely depends on the mode shape for the doubly clamped boundary conditions and is independent of the parameters @xmath398 , @xmath3 , @xmath1 , of the oscillator .
this can be seen from substituting @xmath399 , which yields @xmath400 table [ couple_vals ] shows some numerically obtained values for the bracket in equation that are relevant for the fundamental mode . in the case of an electrostatically tuned resonator ,
the modified nonlinear couplings read @xmath401 where @xmath402 is the factor by which the frequency of mode @xmath115 is reduced due to the presence of the gradient fields .
while this factor is usually intended to be larger than unity for the fundamental mode @xmath403 , it remains close to unity for the higher modes .
phonon transfer between modes is strongly suppressed because of resonance mismatches , as @xmath404 for processes where the phonon number in each mode is not preserved .
one should note that the relevant frequency ratios scale as @xmath405 and the dominant processes of this type affecting the fundamental mode involve its coupling to the next higher mode with the same symmetry , which is the third mode .
in addition , the fundamental mode experiences a modification of its rigidity due to the thermal and quantum fluctuations of higher order modes with @xmath406 .
this effect however can be taken into account by a proper redefinition of the fundamental mode s rigidity .
only pairs of modes @xmath370 with the same parity yield a finite coupling of this type . [ cols="^,^,^,^,^,^",options="header " , ]
a steady state of the mechanical resonator can be probed via a resonant laser drive on an additional cavity mode .
the quantum motion of the nanoresonator is described by the reduced master equation @xcite , @xmath408+\frac{1}{2}\sum\limits_{nm , j}a^{nm}_j{\ensuremath{\mathcal{d}}}({\ensuremath{|n\rangle}}{\ensuremath{\langlem|}})\mu+{\ensuremath{\mathcal{d}}}_{\gamma}\mu\ , .
\label{eq : reduced_mg}\ ] ] here , @xmath409 includes the external mechanical damping via standard lindblad terms and the influence of the lasers is given by the rates @xmath410 here , for example , @xmath411 labels the probe laser with @xmath412 and the other lasers with @xmath413 are used for the steady state preparation , see section [ sec : controlmech ] .
the rates @xmath414 have to be small , i.e. @xmath415 , to assure a weak measurement .
the output power spectrum is given by @xmath416 where the output fields @xmath417 are related to the intracavity fields @xmath418 via the standard input - output relation @xcite , @xmath419 here , we only focus on the output for the probe field and thus drop the index @xmath116 .
the dynamics of the intra cavity field can be described by a quantum langevin equation @xmath420 where @xmath421 and @xmath422 are the fluctuations of the input field in the laser mode and the other bath modes . defining @xmath423 to be a solution for @xmath424
, we can integrate formally and apply a dyson series type expansion to first order in the optomechanical coupling strength @xmath385 , to find for the motion of the cavity modes , cf .
@xcite , @xmath425\\   & \approx   a^{(0)}(t)-{\ensuremath{\mathrm{i}}}\frac{g_{\mathrm{m}}}{2}\sum\limits_{n < m}x_{nm}\left[\frac{{\ensuremath{|n\rangle } } { \ensuremath{\langlem|}}(t)}{{\ensuremath{\mathrm{i}}}\delta_{nm}-\kappa/2}+\text{h.c.}\right]\,.\label{dyson}\end{aligned}\ ] ] the contribution of the fluctuations is included in the free field solution @xmath426 and the input field operators have already been written in a shifted representation @xmath427 , which splits off the coherent part of the input .
we substitute and into and concentrate on the contributions to the sidebands , which to lowest order in @xmath428 are given by the first order two - time correlations of the mechanical motion .
the latter can be calculated from the reduced master equation using the quantum regression theorem .
we find that @xmath429 where @xmath430 are the probabilities to find the resonator in the eigenstate @xmath389 , are the only nonvanishing contributions .
thus the sideband spectrum around the probe laser frequency reads @xmath431 the resulting peak linewidths @xmath432 satisfy @xmath433 , where the nonlinearity @xmath97 is the typical peak distance within the fine structure of one sideband , see figure [ fig : spectrum ] .
t.  rocheleau , t.  ndukum , c.  macklin , j.b .
hertzberg , a.a .
clerk , and k.c .
schwab , nature * 463 * , 72 ( 2010 ) .
oconnell , m.  hofheinz , m.  ansmann , r.c .
bialczak , m.  lenander , e.  lucero , m.  neeley , d.  sank , h.  wang , m.  weides , j.  wenner , j.m .
martinis and a.n .
cleland , nature * 464 * , 697 ( 2010 ) .
teufel , t.  donner , dale li , j.h .
harlow , m.s .
allman , k.  cicak , a.j .
sirois , j.d .
whittaker , k.w .
lehnert , and r.w .
simmonds , nature * 475 * , 359 ( 2011 ) .
j.  chan , t.p .
mayer alegre , a.h .
safavi - naeini , j.t .
hill , a.  krause , s.  grblacher , m.  aspelmeyer , and o.  painter , nature  * 478 * , 89 ( 2011 ) .
e.  verhagen , s.  deleglise , s.  weis , a.  schliesser , and t.  j.  kippenberg , nature * 482 , * 63 ( 2012 ) . t.  j. kippenberg and k.  j.  vahala , science * 321 * , 1172 ( 2008 ) a.d .
armour , m.p .
blencowe , and k.c .
schwab , phys .
lett .  * 88 * , 148301 ( 2002 ) .
f.  marquardt and s.m .
girvin , physics * 2 * , 40 ( 2009 ) .
teufel , t.  donner , m.a .
castellanos - beltran , j.w .
harlow , and k.w .
lehnert , nature nanotech .  * 4 * , 820 ( 2009 ) .
a.   schliesser , r.  rivire , g.  anetsberger , o.  arcizet , t.j .
kippenberg , nat .
phys .  * 5 * , 509 ( 2009 ) .
g.  anetsberger , e.  gavartin , o.  arcizet , q.p .
unterreithmeier , e.m .
weig , m.l .
gorodetsky , j.p .
kotthaus , t.j .
kippenberg , phys .
a * 82 * , 061804 ( 2010 ) .
adler and a.  bassi , science * 325 * , 275 ( 2009 ) .
m. aspelmeyer , t. j. kippenberg , f. marquardt , arxiv:1303.0733 ( 2013 ) .
ehrenfest , z.  phys .
a * 45 * , 455 ( 1927 ) .
p. rabl , phys .
lett . * 107 * , 063601 ( 2011 ) .
a. nunnenkamp , k. borkje , and s. m. girvin , phys .
lett . * 107 * , 063602 ( 2011 ) .
j.  c.  sankey , c.  yang , b.  m.  zwickl , a.  m.  jayich and j.  g.  e.  harris , nat .
* 6 * 707 ( 2010 ) k.  jacobs , phys .
lett .  * 99 * , 117203 ( 2007 ) . s.m .
carr , w.e .
lawrence , and m.n .
wybourne , phys .
b * 64 * , 220101 ( 2001 ) .
v.  peano and m.  thorwart , new j.  phys .
* 8 * , 21 ( 2006 ) .
e.  babourina - brooks , a.  doherty and g.j .
milburn , new j.  phys .  * 10 * , 105020 ( 2008 ) .
i.  katz , a.  retzker , r.  straub , and r.  lifshitz , phys .
* 99 * , 040404 ( 2007 ) . m.j .
hartmann and m.b .
plenio , phys .
* 101 * , 200503 ( 2008 ) .
dykman , m.  marthaler , and v.  peano , phys . rev . a * 83
* , 052115 ( 2011 ) . q.p .
unterreithmeier , e.m .
weig , and j.p .
kotthaus , nature * 458 * , 1001 ( 2009 ) .
landau and e.m .
lifshitz , _ theory of elasticity _ , pergamon , oxford ( 1986 ) .
p.  werner and w.  zwerger , europhys .
lett.*65 * ( 2004 ) .
g.  anetsberger , o.  arcizet , q.p .
unterreithmeier , r.  rivire , a.  schliesser , e.m .
weig , j.p .
kotthaus , and t.j .
kippenberg , nat .
phys .  * 5 * , 909 ( 2009 ) .
m.  pllinger , d.  oshea , f.  warken and a.  rauschenbeutel , phys .
lett .  * 103 * , 053901 ( 2009 ) .
k. h. lee , t. g. mcrae , g.i .
harris , j. knittel , and w. p. bowen , phys .
lett . * 104 * , 123604 ( 2010 ) .
steele , a.k .
httel , b.  witkamp , m.  poot , h.b .
meerwaldt , l.p .
kouwenhoven , and h.s.j .
van der zant , science * 325 * , 1103 ( 2009 ) .
b.  lassagne , y.  tarakanov , j.  kinaret , d.  garcia - sanchez , and a.  bachtold , science * 325 * , 1107 ( 2009 ) . m.  s.  hofmann , j.  t.  glckert , j.  noe , c.  bourjau , r.  dehmel , and a.  hogele , nature nanotech .
* 8 , * 502 ( 2013 ) .
j. waissman , m. honig , s. pecker , a. benyamini , a. hamo , and s. ilani , nature nanotech . * 8 , * 569 ( 2013 ) .
s. stapfner , l. ost , d. hunger , j. reichel , i. favero , and e. m. weig , appl .
phys . lett . *
102 * , 151910 ( 2013 ) .
httel , g.a .
steele , b.  witkamp , m.  poot , l.p .
kouwenhoven and h.s.j .
van der zant , nano lett .
* 9 * , 2547 ( 2009 ) .
a.  eichler , j.  moser , j.  chaste , m.  zdrojek , i.  wilson - rae , and a.  bachtold , nature nanotech .
* 6 * , 339 ( 2011 ) .
b.  kozinsky and n.  marzari , phys .
* 96 * , 166801 ( 2006 ) .
t.  j.  kippenberg , s.  m.  spillane , and k.  j.  vahala , app .
. lett . * 85 * , 6113 ( 2004 ) j.  d.  jackson , _ classical electrodynamics _
, wiley , new york ( 1998 ) i. wilson - rae , phys . rev .
b * 77 * , 245418 ( 2008 ) . l.  g.  remus , m.  p.  blencowe , and y.  tanaka , phys .  rev .
b * 80 * , 174103 ( 2009 ) .
joh , j.  kinder , l.h .
herman , s .- y .
ju , m.a .  segal , j.n .
johnson , g.k .- l .  chan and j.  park , nature nanotech .
* 6 * , 51 ( 2011 ) .
h.  nyquist , phys .
lett . * 32 * , 110 ( 1928 ) b.  c.  stipe , h.  j.  mamin , t.  d.  stowe , t.  w.  kenny , and d. rugar , phys .
* 87 * , 096801 ( 2001 ) .
p.  rabl , s.  j.  kolkowitz , f.  h.  l. koppens , j.  g.  e. harris , p.  zoller , and m.  d.  lukin , nat . phys .
* 6 * , 602 ( 2010 ) .
i.  wilson - rae , c.  galland , w.  zwerger , and a.  imamolu , new j. phys .
* 14 * , 115003 ( 2012 ) i. wilson - rae , n. nooshi , w. zwerger , and t.j .
kippenberg , phys .
* 99 * , 093901 ( 2007 ) .
f.  marquardt , j.p .
chen , a.a .
clerk , and s.m .
girvin , phys .
lett .  * 99 * , 093902 ( 2007 ) . c.  w.  gardiner and m.  j.  collett , phys .
a * 31 * , 3761 ( 1985 )