Abstract
We propose a scheme for preparing nanomechanical oscillators in nonclassical steady states, characterized by a pronounced negative Wigner function. In our optomechanical approach, the mechanical oscillator couples to multiple laser-driven resonances of an optical cavity. By lowering the resonance frequency of the oscillator via an inhomogeneous electrostatic field, we significantly enhance its intrinsic geometric nonlinearity per phonon. This causes the motional sidebands to split into separate spectral lines for each phonon number and transitions between individual phonon Fock states can be selectively addressed. We show that this enables the preparation of the nanomechanical oscillator in a single-phonon Fock state. Our scheme can, for example, be implemented with a carbon nanotube dispersively coupled to the evanescent field of a state of the art whispering gallery mode microcavity.
GENERAL SCIENTIFIC SUMMARY Introduction and background. The growing field of optomechanics aims to explore the quantum physics of mechanically oscillating structures. Via their clamping, such mechanical oscillators are typically coupled to a thermal environment which thermalises their motion. To access their quantum properties, researchers have therefore cooled small mechanical oscillators to their ground states. However, the ground state of purely harmonic oscillators is a Gaussian state and does not show any non-classical features.
Main results. In our paper we introduce a scheme to prepare nanomechanical oscillators in manifestly non-classical states. Our approach electrostatically enhances the intrinsic nonlinearity of doubly clamped oscillators, which typically only becomes relevant at large deflection amplitudes, to such an extent that the nonlinearity per phonon becomes comparable to the optical linewidth of a high-finesse cavity. When combined with optomechanical coupling to the cavity, this approach opens up a variety of possibilities for manipulating the mechanical oscillator by driving the cavity modes with lasers. As one application, we show that suitably adjusted laser fields prepare the mechanical oscillator in a stationary phonon number state for which its quantum behaviour becomes most obvious since its Wigner function reaches the lowest possible value.
Wider implications. Our approach enables a large variety of possibilities for the manipulation of optomechanical systems that we expect to be of considerable use in potential applications for nanomechanical structures in quantum information processing. An observation of the quantum effects for mesoscopic mechanical objects that we predict would furthermore allow to experimentally test potential limitations of quantum mechanics.
1. Introduction
The manipulation of mechanical degrees of freedom by light scattering, at present known as optomechanics, has attracted increasing interest in recent years [1, 2]. Current research in this direction is largely driven by substantial progress in cooling mechanical oscillators to their ground states [3–8] and carrying out displacement measurements close to the standard quantum limit [9–11]. These developments set the stage for investigations of the quantum regime of macroscopic mechanical degrees of freedom, an endeavour that could lead to experimental tests of potential limitations of quantum mechanics [12]. Here, we introduce an optoelectromechanical technique for preparing a mechanical oscillator in a nonclassical steady state with a negative Wigner function by enhancing the anharmonicity of its motion.
As has been known since the initial formulation of quantum mechanics, the dynamics of a purely harmonic quantum system is hard to distinguish from its classical counterpart [13]. As an example, electrical circuits in the superconducting regime do not display quantum behaviour unless they feature a nonlinear element such as a Josephson junction. Along these lines, nanomechanical oscillators that are coupled to nonlinear ancilla systems have been considered [6, 14, 15]. Here, in contrast, we introduce a framework that only makes use of intrinsic properties of the mechanical oscillator. This approach allows one to prepare the latter in nonclassical states that are stationary states of its dissipative motion. The lifetime of these states is thus not limited by the lifetime of the oscillator's excitations.
In fact, doubly clamped mechanical resonators such as nanobeams naturally feature an intrinsic, geometric nonlinear contribution to their elastic energy that gives rise to a Duffing nonlinearity [16–19] and typically only becomes relevant for large deflection amplitudes. To investigate quantum effects of mechanical degrees of freedom [20–22] it is thus of key importance to devise means for enhancing the nonlinearity of nanomechanical structures.
Here we consider inhomogeneous electrostatic fields to decrease the harmonic oscillation frequency of a nanomechanical oscillator [23]. This procedure enhances the amplitude of its zero point motion and therefore leads to an amplification of its nonlinearity per phonon, which eventually becomes comparable to the optical linewidth of a high-finesse cavity. When combined with an optomechanical coupling to the cavity, this approach opens up a variety of possibilities for manipulating the mechanical oscillator by driving the cavity modes with lasers. In particular, the motional sidebands split into separated lines for each phonon number n, cf figure 1. This feature allows us to selectively transfer population from the phonon Fock state |n〉 to |n + 1〉 (|n'〉 to |n' − 1〉) by driving the spectral line corresponding to n (n') phonons in the blue (red) sideband.
Figure 1. Level diagram of the coupled cavity-resonator system in a shifted representation. |nm1j〉 denotes a single-photon state of cavity mode j with n phonons in the mechanical mode, while |0〉c denotes the vacuum (i.e. the steady state of the decoupled cavity). Due to the nonlinearity, the mechanical energy level spacings depend on n and can be addressed individually. Only resonant processes are sketched with a blue detuned laser acting on the transition |0〉m → |1〉m and red detuned lasers acting on transitions |2〉m → |1〉m and |3〉m → |2〉m.
Download figure:
Standard imageTo demonstrate the capabilities of our approach, we show that suitably adjusted laser fields prepare the mechanical oscillator in a stationary state close to a phonon Fock state, i.e. in a state for which the quantum behaviour of an oscillator is most obvious [24]. The highly nonclassical nature of this state manifests itself in a Wigner function that reaches the lowest possible value of −2/π. This can be verified via the power spectrum of an optical probe field [25, 26]. Here, the stationary nature of the motional states considerably simplifies their measurement as compared to phonon Fock states that might be prepared by projective measurements [27–29].
Our scheme can, for example, be implemented with state of the art nanoelectromechanical systems (NEMS) comprising carbon-based mechanical resonators [30, 31] dispersively coupled to high-finesse toroidal [32] or fibre-based [33] microcavities (see figure 2). Carbon-based oscillators are particularly well suited as they combine very small transverse dimensions with ultra-low dissipation [34].
Figure 2. (a) Possible setup with a CNT dispersively coupled to the evanescent field of a microtoroid and tip electrodes used for the electrostatic softening of the mechanical frequency (blue cylinders). (b) NEMS chip carrying the nanotube and the electrodes. The oscillation direction of the nanotube is indicated by the red arrows. In (a), the nanotube and electrodes are not drawn to scale with the toroid and the NEMS chip is not shown for illustration purposes. The placement of the nanotube and electrodes minimizes the degradation of the cavity finesse and maximizes optomechanical coupling.
Download figure:
Standard imageThe remainder of the paper is organized as follows. In section 2, we discuss the Hamiltonian of the oscillating nanobeam, including its geometric nonlinearity. In section 3, we then introduce the optomechanical coupling of the beam to a high-finesse cavity and present the master equation that describes the dynamics of our system including the relevant damping mechanisms. In section 4, we provide an analysis of how driven cavity modes can drive the nanomechanical resonator into a nonclassical steady state by adiabatically eliminating the photon degrees of freedom. In what follows, we analyse important aspects of an implementation in more detail. In section 5, we describe the electrostatic technique that lowers the harmonic oscillation frequency of the beam and hence increases its nonlinearity per phonon. Section 6 describes a possible feasible experimental setup including a carbon nanotube (CNT) coupled to a microtoroidal cavity. We then present examples of states that could be prepared in this setup in section 7. Finally, we discuss a technique to measure the prepared state (section 8), and end the paper with the conclusions and an outlook.
2. The oscillating nanobeam
The dynamics of the fundamental flexural mode of a thin beam can be described by the Hamiltonian [35]
where X is the deflection, P its conjugate momentum, ωm the mode frequency and m* the effective mass. We assume clamped–clamped boundary conditions and apply thin beam theory [36]. This yields for the mode frequency without applying elastic stress or other external forces, and for the nonlinearity which results from the stretching induced by the deflection. Here, L is the suspended length, cs is the phase velocity of compressional phonons and is given by the relevant transverse dimension of the beam times a geometric prefactor.
The single-mode model (1) is valid for sufficiently low phonon occupancies and weak nonlinearities such that phonon–phonon interactions involving higher harmonics can be neglected. Thus, we focus on the latter regime and write Hm in terms of phonon creation (annihilation) operators b† (b) to obtain
with nonlinearity , where is the amplitude of the oscillator's zero-point motion.
Note that λ∝ωm−2 can be significantly enhanced by 'softening' the mode, i.e. lowering ωm. While this can be achieved by applying compressive stress along the beam [35, 37], we follow here an alternative approach to allow for better tunability. We consider an in-plane flexural resonance (cf figure 2(b)) and use tip electrodes to apply a strong inhomogeneous static electric field that polarizes the beam. As a dielectric body is attracted towards stronger fields, this adds an inverted parabola to the potential that counteracts the harmonic contribution due to elasticity. Thereby, as will be borne out in detail in section 5, the oscillation frequency is reduced [23, 38, 39].
3. Optomechanics
We consider a typical optomechanical setup where the displacement X = xZPM(b + b†) of the nanomechanical oscillator dispersively couples to several laser-driven cavity resonances. In a frame rotating with the frequencies of the laser fields ωL,j the system Hamiltonian reads
where Δj = ωL,j − ωj is the detuning of a given laser j from the closest cavity mode with frequency ωj and photon annihilation operator aj, is the driving strength associated with that laser with input power Pin,j and κex is the corresponding external decay rate. The coupling strength corresponds to the shift of the respective cavity resonance as per the oscillator's deflection [40].
We incorporate the damping of the cavity field at rate κ and mechanical dissipation at rate γm in a master equation and shift the cavity normal coordinates to their steady-state value aj → aj + αj, with
Neglecting higher-order terms in the photon–phonon coupling for the regime where 〈a†jaj〉 ≪ |αj|2 [3, 4, 40] leads to the linearized Hamiltonian
where gm,j = 2αjxZPMG0,j is the enhanced optomechanical coupling. The term ∼|α|2(b† + b), causing a static shift of the mechanical mode due to radiation pressure, is omitted here since it can be compensated for by a suitable choice of electrostatic softening fields (see section 5). The master equation for the shifted system now reads
where
with the Bose number at the environment temperature T. In equation (8), we have applied a rotating wave approximation (RWA) that is justified for the high mechanical quality factors we are interested in [40].
Equation (6) describes the dynamics of our system. However, the behaviour of the mechanical mode becomes more apparent after adiabatically eliminating the photon degrees of freedom, as we do in the next section.
4. Reduced master equation
For |gm,j| ≪ κ, we now derive a reduced master equation (RME) for the mechanical motion by adiabatically eliminating the cavity degrees of freedom [40]. We start by applying an RWA to the mechanical Hamiltonian which yields
with the shifted frequency ωm' = ωm + λ' and the nonlinearity per phonon λ' = 6λ. Here, the RWA is warranted provided the relevant phonon numbers n fulfil n2λ' ≪ 6ωm. Within the RWA, the mechanical energy eigenstates are Fock states and the large nonlinearity, , we are interested in naturally leads to Fock number resolved dynamics. In an interaction picture with respect to Hm' transition operators between neighbouring phonon Fock states rotate at frequencies corresponding to the respective mechanical energy level spacings, as
with
and
We further focus on the resolved sideband regime for the softened frequency κ ≪ ωm, which is suitable for efficient back action cooling, as well as , to ensure that the influence of the cavity modes on the mechanical motion overwhelms environmental heating. The resulting RME for the reduced density operator of the nanomechanical resonator, , reads
where the driven cavity modes act as an additional bath for each number transition with damping terms
Note that the frequency shifts
are negligible compared to ωm',λ'. The cavity-induced transition rates between neighbouring Fock states
however, can exceed the mechanical damping , which leads to a significant impact on the mechanical steady state.
The steady state of the RME (12) is diagonal in the Fock basis,
Hence the relations
determine all occupation probabilities Pn and therefore the mechanical state μ. The relations (18) show the central feature of the driven nonlinear nanoresonator. Here implies . For a detuning Δj = ± δn we find that An±,j > An∓,j and the corresponding laser drive therefore enhances (decreases) Pn/Pn−1.
This feature can be used to approach phonon Fock states |n〉, as illustrated in figure 1 for n = 1. Fock states are highly nonclassical and for n odd they exhibit a maximally negative Wigner function at the origin. Ensuring that each laser drive preferentially addresses only one transition |n〉 → |n ± 1〉 naturally requires λ'⪆κ. This can be achieved by electrostatic mode softening as we describe in the next section.
5. Dielectric softening
The harmonic frequency ωm of the nanoresonator can be reduced by adding an electrostatic potential Ves∝ − X2 to the elastic potential, where X is the deflection of the resonator. This is achieved by placing charged electrodes at both sides of the nanoresonator (see figure 2), which generate a strongly inhomogeneous electrostatic field. To estimate the potential that softens the beam, we consider its electrostatic energy line density
Here, y and x are the co-ordinates along the beam and the direction of its deflection, E||,⊥ are the electrostatic field components parallel and perpendicular to the beam and α||,⊥ the respective screened polarizabilities. The electrostatic energy as a functional of the deflection x(y) thus reads
We expand to second order in the transverse deflection x(y) from its equilibrium configuration x(y) = 0. The leading contributions for the fundamental mode with profile ϕ0(y) can be found by considering x(y) = ϕ0(y)X. They read
where Ves,0 is an irrelevant constant and
The linear term Ves,1X shifts the beam's equilibrium position and can be employed to compensate for shifts induced by the dispersive coupling to the cavity fields by satisfying . The quadratic term Ves,2X2 in turn is negative and reduces the fundamental frequency as
At , the buckling instability [35, 37] occurs and by tuning the electrostatic fields below it, one can reduce ωm by any desired factor ζ = ωm,0/ωm. Note that the softening increases the zero point motion amplitude and, hence, also the optomechanical coupling per photon G0xZPM. A particularly useful feature of this dielectric approach to soften the nanobeam is that the electrostatic potential and hence ωm, G0xZPM and λ can be tuned in situ. To illustrate the feasibility of our approach, we now turn to the description of a feasible experimental setup for its implementation.
6. The carbon nanotube and microtoroid
As a suitable device for realizing our scheme, we envisage a setup with a single-walled CNT as the mechanical oscillator (cs = 21 000 ms−1) that interacts with the evanescent field of a whispering gallery mode of a microtoroid cavity [32], cf figure 2. A CNT is a favourable candidate since without softening (ωm = ωm,0) one obtains , indicating that a low effective mass m* and a small transverse dimension are desirable. For a CNT of radius R, one has . We consider a wide-band-gap semiconducting CNT that for the infrared laser wavelengths involved will behave as a dielectric [41]. Taking into account that the transverse optical polarizability is negligible (i.e. depolarization) and approximating the longitudinal one by its static value, the optomechanical coupling can be estimated as (cf [32])
where Vc = πa2cLc is the mode volume of the whispering gallery mode with cavity length Lc and effective mode waist 2ac. ξ is the ratio of the field at the cavity surface to the maximum field inside, d the distance between the chip carrying the nanotube and the microtoroid's surface and κ−1⊥ the decay length of the evanescent field. The correction factor C depends on the type of whispering gallery mode and the placement of the CNT. We assume for simplicity that the NEMS and photonic chips are perpendicular to each other and find for a transverse electric (TE) mode and optimal placement . This optimum involves a small displacement of the CNT midpoint from the equatorial plane of the toroid to ensure coupling to the in-plane motion and an angle of ∼50° between the CNT axis and the equatorial plane of the toroid to circumvent depolarization.
In this setup, the tip electrodes used to lower the mechanical oscillation frequency need to be positioned close to the cavity surface to ensure a sufficiently inhomogeneous electrostatic field while maintaining the optomechanical coupling high enough. Hence, a concomitant degradation of the cavity's finesse due to photon scattering or absorption at the tip electrodes is a natural concern. We find that cavity losses are minimized if one uses TE whispering gallery modes and tip electrodes of sub-wavelength diameter del that are aligned approximately parallel to the equatorial plane of the toroid (cf figures 2(a) and (b)). For this arrangement the dominant loss mechanism is absorption. The ratio of absorbed power Pa to power circulating in the cavity Pc can be estimated to be (cf [42], chapter 8)
where σ is the 2D optical conductivity of the electrodes and θ their misalignment with respect to the equatorial plane of the toroid. The resulting finesse is given by
where is the bare cavity finesse. For relevant cavity parameters we find if we assume θ ⩽ 1°, del ⩽ 10 nm and . Metallic CNTs, for example, fulfil the latter two requirements provided their resonances do not match the relevant cavity frequencies ωj [43].
In the following, we now give some examples of steady states for the mechanical oscillator that can be generated in the setup described here.
7. Examples
To illustrate the potential of our scheme, we show in figure 3 the results for a (10,0) CNT characterized by R = 0.39 nm, α∥ = 142(4π0 Å2) and α⊥ = 10.9(4π
0 Å2) [44]. The nanotube has a length of L = 1.0 μm and is coupled to whispering gallery modes with wavelengths 2πc/ωj ≈ 1.1 μm, ac = 1.4 μm and ξ = 0.4 of a silica microtoroid with finesse [45] and circumference Lc = 1.35 mm. The applied voltages are chosen so that tip electrodes at a distance of 20 nm from the tube induce maximal fields at its axis of E∥ ≈ 1.20 × 107 V m−1 and E⊥ ≈ 1.778 × 106 V m−1 and lower its mechanical frequency by a factor ζ = 4.0. For the optomechanical couplings, we consider gm,j/2π = 21.0 kHz for all j which can be achieved with d = 50 nm, κex/κ = 0.1 and a launched power Pin,j = 1.2 W. The above parameters lead to ωm/2π = 5.23 MHz, λ/2π = 209 kHz and κ/2π = 52.3 kHz. We assume CNT Q-values of Qm = 5 × 106 and Qm = 5 × 105 2 and an environmental temperature T = 20 mK.
Figure 3. Steady state of an oscillating (10,0) CNT with L = 1.0 μm and Qm = 5 × 106 for T = 20 mK, and Lc = 1.35 mm. (a) Fock state occupation probabilities Pn. (b) The Wigner function W. These results are obtained from a numerical solution of equation (6).
Download figure:
Standard imageFor these parameters and for laser detunings Δ1 ≈ δ1 and Δ2,3 ≈ −δ2,3 a single phonon Fock state can be prepared with 91% fidelity (for Qm = 5 × 106), as shown in figure 3. This steady state of the mechanical oscillator has a Wigner function that is clearly negative at the origin, W(0,0) = −0.53 for Qm = 5 × 106 and W(0,0) = −0.30 for Qm = 5 × 105.
The results shown in figure 3 have been obtained from a numerical solution of equation (6) with numerically optimized detunings Δj. These numerical results do not rely on the rotating wave approximation used in equation (9). Here the fast rotating terms lead to corrections to the energy level spacings δn that can be comparable to κ. Hence to match the laser detunings with sufficient accuracy to the energy level spacings δn, the latter have been determined numerically.
In turn, for the more moderate parameters: L = 1.7 μm, , Lc = 1.80 mm (κ/2π = 64.2 kHz), ζ = 3.3 (which imply ωm/2π = 2.13 MHz and λ/2π = 85.5 kHz), launched laser powers Pin,1/2 = 22 mW and Pin,3 = 44 mW, Qm = 1.5 × 106 and T = 30 mK, we still find a significant negative peak of depth Wmin = −0.15 (see figure 4).
Figure 4. Steady state of an oscillating (10,0) CNT with L = 1.7 μm and Qm = 1.5 × 106 for T = 30 mK, and Lc = 1.8 mm. (a) Fock state occupation probabilities Pn. (b) The Wigner function W. These results are obtained from a numerical solution of equation (6).
Download figure:
Standard imageWe note that the deep resolved sideband and undercoupled conditions used here, where the frequencies of the driving lasers are detuned by several cavity linewidths from the optical resonances of the cavity and losses into the coupling fibre are only a small fraction of the total losses, imply that the power absorbed inside the cavity is much smaller than the launched power (below 30 μW in the example of figure 3 and below 7 μW in the example of figure 4). Whereas an intracavity absorption ∼1 μW has already been shown to be compatible with a cryogenic environment [46], even higher finesses and lower absorption could be attained by using crystalline resonators [47]. Additionally, the CNT polarizabilities could be enhanced well above the static values we have assumed by profiting from the corresponding excitonic optical resonances [48], allowing for a substantial reduction in the laser powers required.
Furthermore, in both examples, the voltages applied to soften the frequency of the CNT nanomechanical resonator generate electric fields at the CNT that are well below its threshold for dielectric breakdown and are comparable to fields applied in recent experiments [23].
Finally, we note that our approach can be generalized to the preparation of Fock states with n > 1.
8. Measurements
To verify whether a nonclassical steady state with a negative Wigner function W has been successfully prepared, one can measure sidebands of the output power spectrum S(ω) of an additional cavity mode weakly driven by a probe laser with frequency ωL. Within our quantum noise approach [40] and under the same conditions for which equation (12) holds, we find for the sidebands (i.e. where ω ≠ ωL) (figure 5)
with sideband linewidths
Figure 5. Power spectrum S(ω) of a resonant probe laser with the sidebands split into well-resolved lines corresponding to the phonon numbers n. This example shows the result of equation (28) for the parameters used in figure 3.
Download figure:
Standard imageSince the negativity of the Wigner function is in our scheme expected to be maximal at the origin where , the nonclassicality of the state can be read off from S(ω). For a probe laser that is resonant with a cavity mode, one has An+ = An− and the ratios of peak intensities S(ωL + δn)/S(ωL − δn) ≈ Pn/Pn−1 together with allow us to determine the Pn assuming that occupations above a certain Fock number are negligible.
9. Conclusions and outlook
We have shown that a dielectrically enhanced nonlinearity in nanomechanical oscillators enables controlling their motion at a Fock state resolved level. The scheme set forth opens new avenues for optomechanical manipulations of nanobeams, in particular by exploiting the in situ tunability of the electrostatic softening approach introduced.
Acknowledgments
This work is part of the DFG Emmy Noether project HA 5593/1-1 and was supported by Nanosystems Initiative Munich. The authors thank J Kotthaus, T J Kippenberg and A Schliesser for useful discussions.
Footnotes
- 2