Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Strong opto-electro-mechanical coupling in a silicon photonic crystal cavity

Open Access Open Access

Abstract

We fabricate and characterize a microscale silicon opto-electromechanical system whose mechanical motion is coupled capacitively to an electrical circuit and optically via radiation pressure to a photonic crystal cavity. To achieve large electromechanical interaction strength, we implement an inverse shadow mask fabrication scheme which obtains capacitor gaps as small as 30 nm while maintaining a silicon surface quality necessary for minimizing optical loss. Using the sensitive optical read-out of the photonic crystal cavity, we characterize the linear and nonlinear capacitive coupling to the fundamental ωm/2π = 63 MHz in-plane flexural motion of the structure, showing that the large electromechanical coupling in such devices may be suitable for realizing efficient microwave-to-optical signal conversion.

© 2015 Optical Society of America

1. Introduction

Microelectromechanical systems (MEMS) are a widespread technology platform with a vast number of applications. MEMS devices are found, for instance, in a variety of hand-held electronic devices, often as accelerometers [1], microphones [2] or pressure sensors [3]. Recently, MEMS have been proposed for energy harvesting applications [4], ultra-high resolution mass sensors [5] and as a suitable candidate for the development of biological sensors in lab-on-a-chip technologies [6]. New generation of MEMS have critical dimensions down below the microscale, and into the nanoscale, opening the possibility to integrate these devices with other nanotechnologies. In the case of nanophotonics, the emerging field of cavity optomechanics uses the radiation pressure force to probe and control the state of a mechanical actuator embedded in an optical cavity. Optomechanical nanophotonic devices, such as photonic crystal “zipper” cavities [7] and optomechanical crystals (OMCs) [8], have been proven effective for near quantum-limited position read-out [9, 10] and strong back-action effects, as shown by the cooling of a mechanical mode to its quantum ground state of motion [11].

The integration of MEMS with optomechanical devices may be useful both for microphotonic circuits, where MEMS may be employed to tune the optical properties of devices, as well as for MEMS sensors, where optomechanical devices may be used for shot-noise-limited readout and back-action modification of the mechanical response. Moreover, an integrated MEMS-optomechanics technology could allow for up-conversion of low-frequency electrical signals to an optical carrier, mediated by an intermediate mechanical transducer. The ultimate goal of such a conversion scheme would see the realization of coherent, quantum frequency translation between an optical and microwave cavity which shares the same mechanical resonator [1216]. Different approaches to realize such frequency translation include the use of a silicon nitride membrane vertically stacked within an electronic and optical cavity [17, 18], and the creation of piezoelectric nanobeam optomechanical crystals [19].

Our approach to integrating optomechanical and electromechanical devices utilizes a silicon (Si) optomechanical photonic crystal whose mechanical degree of freedom is shared with an electrical capacitor [20]. Using a nanoslotted planar photonic crystal slab, it is possible to localize the optical field to the center slot of the slab and the capacitor to the other outer edges of the slab, with both optical and elecrostatic fields connected to the same slab motion. This avoids losses in the optical path, and in the case of superconducting circuits, avoids optically-induced electrical losses from the breaking of Cooper pairs. In this kind of planar device, the capacitive element is effectively one dimensional, being formed by two parallel metal wires. A strong coupling between the mechanical motion and the optical or electrical mode (in what follows we refer to the lower frequency mode as the electrical mode) can be realized by making the electromagnetic mode volume commensurate with the acoustic wavelength of the mechanical resonance. In the case of near-infrared optics and GHz mechanical resonances, one finds there is a common wavelength scale, whereas for radio-frequency or microwave electrical modes, the mechanical and electrical length scales are vastly different. In the microwave frequency range this requires decreasing the capacitor gap size to tens of nanometers, where metallic or superconducting boundaries are used to effectively ”squeeze” the electric field into a small volume [21, 22].

In this Article, we push the fabrication limits of the silicon-on-insulator (SOI) cavity-electro-optomechanics platform first presented in Ref. [20], to achieve large electromechanical coupling by engineering a narrow electrode gap while retaining a high optical Q-factor. To this end, we use an inverse shadow-mask evaporation which allows us to fabricate electrode gaps as small as d ≈ 30 nm. Compared to other fabrication techniques, such as focused ion beam milling [22], this method maintains a pristine semiconductor surface, avoiding damage to the optical resonator [23]. An outline of the paper is as follows. We begin with a review of the optical, electrical, and mechanical design of the structure, followed by a description of the methods used for device fabrication. Optical and mechanical characterization of the device are then presented. This is followed by a characterization of both the linear and nonlinear electromechanical coupling. We conclude by discussing the potential application of these devices for efficient and noise-free microwave-to-optical signal conversion.

2. Design and Fabrication

As shown in Fig. 1, the electro-optomechanical device studied here is based around a silicon thin-film photonic crystal in which a linear waveguide is formed around a central nanoscale air slot (a so-called W1 slotted waveguide) [20, 24]. An optical resonant cavity is formed from the waveguide by creating a defect along the axial length of the waveguide in which the parameters of the waveguide are slowly modified. This results in an optical mode confined in the s ∼ 80 nm air slot and localized to the defect region, as shown in the finite-element-method (FEM) simulation of Fig. 1(a). Mechanical motion of the structure is allowed by undercutting, and suspending the Si device layer. Two additional gaps are fabricated on the outer edge of the two photonic crystal slabs to accommodate capacitor electrodes which connect the mechanical motion of the slabs to an electrical circuit (see Fig. 1c). The whole slab structure is clamped on the ends to the underlying SiO2 (BOX) layer and Si substrate, resulting in a fundamental in-plane mode with simulated frequency of ωm/2π = 67 MHz. The deformation profile, |Q(r)|, of the differential motion of the two slabs is shown in Fig. 1(b). A spatial separation between light and metals, as in our device, has a twofold benefit. From the photonic side it makes metallic-related losses negligible, and from the electrical side it avoids any stray light on the electrodes, which can be crucial if superconducting materials are used.

 figure: Fig. 1

Fig. 1 (a) Simulated Ey(r) electric field component of the optical mode. (b) Total deformation |Q(r)| of the fundamental planar flexural mode of the Si slab. (c) Scanning electron microscope (SEM) image of the top-view of a typical realized device. The scale bar is 5 μm. The zoomed-in view to the right depicts the electrode gap, whose size is d < 30 nm. (d) Illustration of the inverse shadow mask technique. A mask is created with ZEP resist to define the electrical contacts, while the capacitor electrodes are defined using the etched Si edges. After lift-off and full release, the metal within the gap is removed together with the SiO2 buffer layer. The scale bar is 5 μm. The purpose of the small trench on the top of the picture is to redistribute the residual silicon stress in order to avoid out-of-plane bending of the device.

Download Full Size | PDF

The main difference between the device studied here and the one reported in Ref. [20] resides in the electrode gap size, which has been reduced from 250 nm to less than 30 nm using an inverted shadow-mask evaporation technique in a two-layer lithography process (see Fig. 1d). In the first step of the device fabrication, a ZEP resist mask is defined with electron beam lithography. The 220 nm Si device layer of the SOI wafer used in this work is dry etched with an inductively-coupled SF6/C4F8 plasma reactive-ion etching process. After a H2SO4/H2O2 Piranha clean, the device layer is partially undercut in a 15 second hydrofluoric acid (HF) wet etching step. The electrical circuit is then defined in a ZEP lift-off process using an aligned electron beam lithography step, followed by deposition of chromium (5 nm) and gold (50 nm) layers and then lift-off. In this last step, the mask for the capacitor electrodes is defined using the partially underetched silicon slot, resulting in a metal gap size defined by the etched slot itself. After lift-off in ZDMAC, the entire optomechanical device is released in HF, whereby the metal deposited between the capacitor electrodes is removed along with the underlying BOX layer. A SEM of the final device is shown in Fig. 1(c). As can be seen in the enlarged view, our fabrication technique enables the realization of clean and narrow vacuum electrode gaps of width d < 30 nm.

3. Optomechanical characterization

The experimental set-up used to probe the optical, mechanical, and electrical properties of the fabricated device is shown in Fig. 2(a). The device fundamental optical resonance is probed through optical transmission measurements using a tunable external cavity semiconductor diode laser (Newfocus Velocity series) whose frequency tuning is calibrated with an unbalanced fiber Mach-Zender interferometer. Optical coupling to a given device is achieved via a tapered and dimpled optical fiber probe, which when placed in near-field of the photonic crystal cavity allows for evanescent coupling of light between the fiber and cavity [25]. An optical transmission scan, shown in Fig. 2(b) for a typical device (device A in what follows), shows a fundamental optical resonance with center wavelength λc = 1522.94 nm and intrinsic quality factor Qc = 8.9 × 104.

 figure: Fig. 2

Fig. 2 (a) Experimental setup for sample characterization. The tunable diode laser is calibrated with a high-finesse Mach Zender interferometer (MZ) [27] and can be locked using a wavemeter (WM) for feedback. After being attenuated (VOA) and polarization controlled (FPC), the evanescent field of a tapered optical fiber is used to couple light into the photonic crystal cavity (PCC). The optical transmission is then measured with a low bandwidth photodetector (D1) in the case of optical spectroscopy, or via a high-speed photodetector (D2) connected to a Spectrum Analyzer (SA) for analysis of the mechanical spectrum. In addition, voltage sources (VDC), a signal generator (VAC) or a Network analyzer (NA) were used to measure the electric response of the device. (b) Measured normalized optical transmission through the photonic crystal of device A. The red line is the best Lorentzian fit of the resonance, showing an intrinsic optical Q-factor of Qc ∼ 9 × 104. (c) Power spectral density of the optically transduced thermal Brownian motion of the structure of device A. The strong resonance peak at ωm/2π = 63 MHz is the fundamental in-plane mechanical mode of the slab structure, whereas the smaller peaks represent weakly coupled (predominantly) out-of-plane modes of the slab. The red line is a multi-Lorentzian fit to the transduced mechanical spectrum. (d) Simulated capacitance as a function of gap size, d, for the capacitor on one side of the slab structure. The data fits with a power law (red line) Cm = [6.933/d0.533] fF, with d in units of nanometers. Inset: simulated Ey electrostatic field amplitude within the capacitor. The field isolines are shown in black.

Download Full Size | PDF

The mechanical mode studied in this work is not the fundamental in-plane differential slab mode, but rather the fundamental in-plane mode of only a single slab. This is an artifact of the measurement technique we employ, in which to mechanically stabilize the optical fiber taper we place it in direct contact with one of the photonic crystal slabs. This effectively decouples the one slab from the other, and thus we only measure and actuate the motion of the free slab without the taper on it. The fundamental in-plane mode of the free slab still modifies the air slot gap size in the center of the photonic crystal, and thus induces a frequency shift of the optical cavity resonance which we quantify by an optomechanical coupling parameter gOM (defined below) in units of GHz/nm. The transmitted optical power for a laser, frequency locked on the side of the optical cavity resonance, carries a signal corresponding to the thermal Brownian motion of the structure, as shown in Fig. 2(c) for device A. The measured fundamental mechanical mode of this device is centered at ωm/2π = 63 MHz with a corresponding quality factor of Qm = 150, limited by atmospheric pressure squeeze-film damping [26]. The smaller peaks in the optically-transduced mechanical spectrum are (predominantly) out-of-plane slab modes which are weakly coupled to the optical cavity resonance.

The optomechanical coupling parameter, gOM, is defined as the fractional change in the energy stored in the optical resonance per unit displacement of the mechanical resonance. The fractional energy shift (δŪE) can be numerically calculated with a perturbative approach, evaluating an integral of the electric and displacement fields of the optical resonance over the dielectric boundaries of the structure (here we assume no change in the magnetic energy due to the fact that the magnetic permeability of most materials is approximately unity at optical frequencies) [8, 28]:

δU¯E=V(Q˜(r)n)(Δε|E|2Δε1|D|2)d2rVε|E|2d3r,
where ε is the dielectric constant, Δε ≡ (ε1 − ε2) is the difference in the dielectric constant across the boundary between material 1 and material 2, Δε−1 ≡ (1/ε1 − 1/ε2), and E (D) is the parallel (perpendicular) component with respect to the boundary, ∂V, of the electric (displacement) field. Here, a generalized coordinate for the mechanical resonance of u = max(|Q(r)|) is chosen, corresponding to a normalized displacement field of = Q/max(|Q|) in eq. (1). The optomechanical coupling, representing the optical resonance frequency shift per unit displacement amplitude u of the mechanical resonance, is given by gOM = (1/2)ωcδŪE, where ωc is the optical resonance frequency and the factor of 1/2 accounts for the energy in the magnetic field which is decoupled from the mechanical motion. For the fundamental optical mode coupled to the fundamental in-plane mechanical mode of a single slab, the optomechanical coupling is evaluated to be gOM/2π = 36 GHz/nm. This value is in good agreement with the experimentally measured value of 37 ± 6 GHz/nm for device A, determined from the radiation pressure induced mechanical frequency shift (see for example Ref. [20]). Note that we have implicitly chosen a positive amplitude u to correspond to outward motion of the photonic crystal slabs as shown in Fig. 1(b), resulting in a reduced capacitor gap, d, and an increased central slot width, s. In what follows we continue to use the same generalized mechanical coordinate, u = max(|Q|), in order to be self-consistent with the calculated value of gOM.

4. Linear electromechanical coupling

The motion of the fundamental in-plane mechanical mode can also be detected and actuated via the capacitor on the outer edge of the photonic crystal slab. The simulated capacitance (for a single side of the slab) is found to be in the range 0.5–2 fF, scaling with the gap d as shown in Fig. 2(d). In analogy with the optomechanical coupling, it is possible to quantify the capacitive, linear electromechanical interaction through the coupling parameter Ge,1 ≡ (1/Cm)(∂Cm/∂u) [17], where Cm is the capacitance modulated by the mechanical mode and u is again the chosen generalized coordinate of the mechanical displacement. Ge,1 can be evaluated starting from the electrostatic energy in a capacitor biased with a fixed charge Q, Ues = Q2/2Cm. For fixed charge, the fractional electrostatic energy shift due to motion of the mechanical resonance can be numerically calculated for an arbitrary mechanical displacement profile from the unperturbed fields in the capacitor using the same integral as for the calculation of the optomechanical coupling in eq. (1), with Ge,1 = −δŪE. Here we assume perfectly conducting boundary conditions and zero fields within the metal wires, a good approximation at microwave frequencies and below. Typically, one would like to couple the mechanical motion to an electrical resonant circuit. In such a case, the coupling capacitance is closed by an inductor (L) to form an LC resonant circuit of frequency ωLC/2π = (LC)−1/2. Assuming capacitive coupling only to the mechanical resonator, the corresponding electromechanical coupling is given by, gEM = (−ωLC/2)Ge,1, in direct analogy to the optomechanical coupling. Addition of the inductor usually accompanies an unwanted parasitic capacitance, Cs, which is not coupled to the mechanical resonator and which reduces Ge,1 by a participation factor ηe = Cm/(Cm + Cs). Here we will concern ourselves primarily with the coupling parameter Ge,1, however, in the conclusion we will further discuss coupling to a microwave LC resonator.

To measure the linear electromechanical coupling parameter we apply a voltage to the electrodes and measure the resulting mechanical displacement using optical read-out. This is done for both a static voltage (VDC) and for a small modulated voltage (VAC) at half the resonance frequency of the fundamental in-plane differential mode. For example, Fig. 3(a) shows a waterfall plot of the transmission through the optical cavity for an applied bias voltage of VDC = 0 to 5.6 V for device A. In order to determine the electromechanical coupling from this tuning data, we first consider the force exerted on the fundamental in-plane mechanical mode for a potential difference V across the capacitor,

Fcap=12CmV2Ge,1,
where again Ge,1 is specific to the amplitude coordinate u of a given mechanical resonance. For a static voltage, the resulting mechanical displacement amplitude u due to Fcap is inversely proportional to the mode effective spring constant, keff=meffωm2. Here, meff (= 10.4 pg) is the effective motional mass of the fundamental in-plane mechanical mode of a single slab defined as:
meff=Q*(r)ρ(r)Q(r)d3rmax(|Q|2),
where ρ is the mass density of Si. This definition of meff is consistent with our choice of definition of u corresponding to the maximum amplitude of the mechanical displacement profile. The effective spring constant, obtained by combining the simulated meff and the measured mechanical frequency, is equal to keff = 1.8 kN/m for device A. This agrees within 18% of the full numerical simulation of the deformation under a constant load applied to the center point of the mechanical mode. The mechanical deformation shifts the optical resonance frequency according to,
Δωc,DC(V)=gOMCmGe,12keffVDC2=αVDC2,
where α ≡ (gOMCmGe,1)/(2keff) is the optical tunability of the structure.

 figure: Fig. 3

Fig. 3 (a) Waterfall plot of the optical transmission spectra for the fundamental optical cavity resonance as a function of DC voltage (VDC = 0 to 5.6 V, 0.2 V steps). The x-axis has been normalized according to the optical center wavelength at zero applied bias, λc,0. (b) Measured optical resonance wavelength (open circles) versus VDC2. Solid black curve corresponds to a linear fit to the measured data. (c) Linear coupling parameter extracted from different experiments (DC (circles), AC (triangles)) and for different samples (A (red) and B (blue)). The solid black line is the numerically simulated linear electromechanical coupling using eq. (1). The dashed green line is the linear electromechanical coupling calculated using the functional dependence of the capacitance versus gap d shown in Fig. 2(d).

Download Full Size | PDF

Alternatively, if a modulated voltage VAC = V0 cos(ωACt) is applied to the capacitor, the resulting mechanical displacement is filtered by the mechanical response function, with the maximum displacement being enhanced by the mechanical Q-factor for an on-resonance capacitive force, ωAC = ωm/2. In this case, the time-average of the optical transmission spectrum assumes a double-dip lineshape (see Ref. [20]) with a separation between the minima given by,

Δωc,AC=αQmV02.

Using eqs. (4) and (5) we can extract Ge,1 from the experimental static and modulated voltage tuning curves, respectively. As expected from eq. (4), the applied static voltage blue shifts the optical resonance as shown in Fig. 3(a) for device A with capacitor gap (from SEM images) of d = 28 nm. The optical resonance wavelength versus VDC2 is shown in Fig. 3(b) using the same color scale, from which a linear fit yields an optical tunability of α/2π = 0.31 GHz/V2. AC measurement of α for device A at a drive frequency of ωAC = ωm/2 yields a small (20%) difference compared to the static voltage measurement, attributable to the impedance mismatch in our AC drive circuit. A second device (device B) with capacitor gap d = 46 nm was also measured using the AC tuning method. Using the measured optomechanical coupling (gOM) along with the simulated capacitance (Cm) and motional mass (meff), we show a plot of the inferred linear electromechanical coupling parameter Ge,1 from DC and AC tuning measurements in Fig. 3(c) for both device A and device B. The measured linear electromechanical coupling agrees well with the numerically simulated curve using eq. (1) (solid black curve).

5. Nonlinear electromechanical coupling

The small capacitor gaps in our devices makes nonlinear terms in the electromechanical interaction relevant. Electromechanical nonlinearities can be used to generate squeezing [2931] and mechanical parametric amplification [32, 33], as well as for logic operation in the classical [34] and quantum regime [35]. Assuming the voltage across the capacitor electrodes follows the applied voltage during the mechanical motion of the electrodes (This assumption requires that the drive circuit be able to provide the current necessary to the electrodes such that as the electrodes undergo mechanical motion the capacitor voltage follows the applied voltage. This feedback from the drive circuit modifies the electrostatic force curve similar to the optical feedback in cavity-optomechanics, resulting in a dynamic spring effect. A wholly different force curve and nonlinear coupling parameter results if we assume a closed capacitor system with fixed charge on the electrodes.), and expanding the capacitive force to linear order in amplitude u, we can define the second order nonlinear coupling parameter as Ge,2 = (1/Cm)(∂2Cm/∂u2). The effective spring constant of the mechanical system is ktot = keff − (1/2)CmV2Ge,2 [29, 36], resulting in an electrical shift of the mechanical resonance frequency given by,

Δωm=(ωmCmGe,24keff)V2=βV2,
where β ≡ −(ωmCmGe,2)/(4keff) is defined as the mechanical tunability. From the mechanical frequency shift versus applied voltage, we can obtain the nonlinear coupling parameter Ge,2 in a similar fashion to that used to determine the linear parameter Ge,1.

Before proceeding to measurements of the nonlinear electromechanical coupling, it is instructive to consider again the capacitor force expression. If a voltage signal with mixed AC and DC components is fed to the capacitor an additional resonant term appears:

FcapV2=VDC2+VAC2+2VDCVAC.
The third, mixed term, explicitly given by 2VDCV0 cos(ωACt), has a maximum mechanical response at ωAC = ωm and can be suppressed or enhanced by controlling the DC bias. This is a useful and well known property, which allows control of electromechanical nonlinearities [37, 38]. The resonant nature of this mixed term is also useful to perform homodyne detection of the coherent mechanical oscillations induced by an AC drive. To this end, we use a network analyzer (NA) of which port 1 is connected to the capacitor electrodes while port 2 is connected to the optical photodetector used to read-out the mechanical motion (see Fig. 2(a)). The S-parameter S21 in such a scheme will therefore carry the amplitude and phase response of the mechanical resonator to the electrical driving force.

Measurements of the nonlinear electromechanical response were performed on device B, whose fundamental in-plane mechanical mode has a resonance frequency of ωm/2π ≈ 49 MHz, intrinsic mechanical Q-factor Qm = 344, and spring constant keff = 1.09 kN/m. In order to avoid spurious transduction of large amplitude mechanical motion in these measurements, a fixed probe laser wavelength detuned ∼ 50 nm from the optical resonance of device B is used. Measurement of the electrostatic modification to the mechanical frequency is first measured, using a weak AC drive voltage of V0 <1 mVrms and DC voltages varying from 0.2 to 6 V. The frequency of VAC is swept using port 1 of the NA and the photodetected signal of the mechanical response is measured on port 2. For these drive levels (VDCV0), the mechanical frequency shift in eq. (6) is dominated by the VDC2 term. As can be seen in the waterfall plot of Fig. 4(a), the mechanical resonance frequency of device B red-shifts with applied DC voltage, corresponding to electrostatic softening. The extracted resonance frequency scales quadratically with the applied DC bias, yielding a mechanical tunability parameter of β/2π = −3.989 kHz/V2.

 figure: Fig. 4

Fig. 4 Nonlinear electromechanical response of device B. (a) Measured amplitude response (|S21|) for a weak AC voltage probe (V0 <1 mVrms) and DC bias ranging from 0.2 V (blue curve) to 6 V (red curve) in 0.2 V steps. Successive spectra are offset from the VDC = 0.2 V spectrum in 1 dB increments. (b) Measured amplitude response (|S21|) of the optically transduced mechanical motion to a strong near-resonant electric drive (ωAC ∼ ωm) with fixed VDC=6 V and V0 ranging from 277–870 mVrms. Increasing the drive amplitude, the mechanical Lorentzian peak (red) starts to assume the form of a Duffing oscillator (green). After an instability threshold, the mechanics enters a frequency entrainment regime in which it is completely in phase with the AC drive (blue). (c) Measured relative phase response (ϕ(S21)) of the optically transduced mechanical motion for the same drive conditions as in (b). Note that the phase slope is about the same inside and outside the frequency entrainment region indicating that the phase of the mechanical response is locked to the phase of the drive. Successive spectra in (b) and (c) are offset from the VAC = 277 mVrms spectrum in 1 dB increments.

Download Full Size | PDF

This corresponds to a nonlinear coupling parameter of Ge,2 = 3.86 × 10−4 nm−2, in reasonable agreement with the numerically simulated value of 4.3 × 10−4 nm−2.

Measurements of device B under a strong near resonant AC drive (V0 = 277 to 870 mVrms) and fixed DC bias of VDC = 6 V are shown in Fig. 4(b). With increasing AC drive amplitude the mechanical resonance evolves from a symmetric Lorentzian (red curve), to an asymmetric lineshape characteristic of a Duffing oscillator (green curve), and finally into an instability regime (blue curve) in which the mechanical resonance becomes entrained by the driving tone. The amplitude in this instability regime becomes roughly constant versus drive frequency, and the phase response assumes a constant slope versus frequency in the entrained region as shown in Fig. 4(c). Such instabilities are known to occur for parametric driving at frequencies near 2ωm/n, n being an integer corresponding to the nth instability “tongue” [39]. As we are driving near resonance, the parametric instability in this case would correspond to the second tongue (n = 2), with a theoretical threshold drive amplitude given by [39]:

V0,th=2keffQm1/2CmGe,2VDC.
The estimated parametric instability threshold using the Ge,2 value from the mechanical softening measurement, is V0,th = 39 Vrms, far larger than the measured instability threshold of only V0,th ≈ 600 mVrms. This large discrepancy likely indicates that the influence of the large resonant driving force (which is present in addition to the resonant parametric drive) cannot be ignored, and that other mechanical nonlinearities such as the cubic Duffing term also play a role in the onset of the instability.

6. Conclusion

To conclude, we consider application of the demonstrated electro-opto-mechanical device to quantum conversion between electrical and optical signals [13, 14]. In particular, we consider a system in which the capacitor electrodes (Cm ∼ 1.2 fF) of the current device are connected together through a wire inductor, forming a lumped element LC resonator in the microwave frequency range. For such a resonant circuit, the linear electromechanical coupling can be quantified by g0,egEMxzpf, where gEM = (−ωLC/2)Ge,1 and xzpf = (ħ/2meffωm)1/2 is the zero-point amplitude of the mechanical resonator. A similar relation exists on the optomechanical side, with g0,ogOMxzpf. Physically, g0,e and g0,o represent the shift in the microwave and optical resonance for an amplitude of motion equal to the zero-point amplitude of the mechanical resonator, respectively. For the optomechanical device studied here, xzpf = 3.48 fm and g0,o/2π = 125 kHz. For an LC resonator of frequency ωLC = 10 GHz, compatible with current superconducting quantum circuits [40], the naively estimated electromechanical coupling rate would be g0,e/2π = 250 Hz. However, one must also consider the stray capacitance associated with adding a large inductance (for the lumped element LC resonators of interest here, the dominant parasitic capacitance arises from the proximity of different regions of the meandered or coiled wire of the inductor [41,42]). Using a planar rectangular spiral inductor with large inductance per wind [42], numerical simulations indicate that a 10 GHz microwave resonance can be realized with a 54 nH spiral inductor formed on a Si slab (thickness 220 nm) with wire widths of 500 nm and inter-wire spacing of 500 nm. Such an inductor fits within a 81 μm square, and has a stray capacitance estimated to be Cs = 3.26 fF. Here we have included a 3 μm vacuum layer, corresponding to the undercutting of the BOX layer, between the Si device layer and a 725 μm Si handle wafer. This reduces the estimated electromechanical coupling by a factor of ηe = Cm/(Cm + Cs) = 0.27, to a value of g0,e/2π = 70 Hz for the device studied here.

The above analysis should be compared against the aluminum superconducting LC resonators of Ref. [42], which employ a vertically-layered vacuum gap capacitor with drumhead mechanical modes around 10 MHz, and have been used to realize strong electromechanical back-action sufficient to cool the mechanical mode to its quantum ground-state and realize efficient microwave-to-optical signal conversion [18]. The electromechanical coupling in such devices is in the g0,e/2π ∼ 200 Hz range, comparable to the estimated value for the nanoslot devices of this work. Remaining challenges to realizing efficient microwave to optical signal conversion in planar Si devices similar to those studied here, include the attainable mechanical and microwave Q-factor at milliKelvin temperatures. Cryogenic temperatures are required both to reduce the thermal noise in the mechanics, as well as to limit the loss in the microwave circuit. Efficient optical coupling to similar Si optomechanical devices at milliKelvin temperatures has recently been demonstrated [43], as have mechanical Q-factors greater than 105 at temperatures of T ∼ 10 K [44]. A far more significant challenge will likely be the realization of low-loss superconducting resonators on the SOI wafer platform, in which a high resistivity (> 1 kOhm cm) Si device layer must be used [45] and the Si surface must be appropriately passivated [46]. Finally, in order to realize microwave-to-optical conversion free of quantum noise (spontaneous Stokes scattering), one must operate in both the sideband resolved limit for the microwave resonator (a relatively straightforward requirement with MHz mechanical resonators) and the optical resonator. Attaining sideband resolution in the optical domain with the photonic crystal nanocavities of this work, and their accompanying large optical linewidth (κ/2π ≳ 200 MHz [47]), requires mechanical resonances with frequency ≳ 100 MHz. The devices presented here represent a step toward this limit (ωm/2π ∼ 60 MHz), but further consideration of the mechanical design, including optomechanical crystal cavities [8], is currently underway.

Acknowledgments

This work was supported by the DARPA MESO program, the AFOSR Hybrid Nanophotonics MURI, the Institute for Quantum Information and Matter, an NSF Physics Frontiers Center with support of the Gordon and Betty Moore Foundation, and the Kavli Nanoscience Institute at Caltech. AP gratefully acknowledge funding from EU through Marie Curie Actions, project NEMO (GA 298861). AT acknowledges partial financial support from the ERC through the advanced grant SoulMan.

References and links

1. N. Yazdi, F. Avazi, and K. Najafi, “Micromachined inertial sensors,” Proc. IEEE 86, 1640–1659 (1998). [CrossRef]  

2. T. Tajima, T. Nishiguchi, S. Chiba, A. Morita, M. Abe, K. Tanioka, N. Saito, and M. Esashi, “High-performance ultra-small single crystalline silicon microphone of an integrated structure,” Microelectron. Eng. 67–68, 508–519 (2003). [CrossRef]  

3. W. P. Eatony and J. H. Smith, “Micromachined pressure sensors: review and recent developments,” Smart Mater. Struct. 6, 530–539 (1997). [CrossRef]  

4. K. A. Cook-Chennault, N. Thambi, and A. M. Sastry, “Powering MEMS portable devices: a review of non-regenerative and regenerative power supply systems with special emphasis on piezoelectric energy harvesting systems,” Smart Mater. Struct. 17, 043001 (2008). [CrossRef]  

5. M. S. Hanay, S. Kelber, A. K. Naik, D. Chi, S. Hentz, E. C. Bullard, E. Colinet, L. Duraffourg, and M. L. Roukes, “Single-protein nanomechanical mass spectrometry in real time,” Nature Nanotech. 7, 602–608 (2012). [CrossRef]  

6. A. Manz, N. Graber, and H. M. Widmer, “Miniaturized total chemical analysis systems: a novel concept for chemical sensing,” Sens. Actuat. B: Chem. 1, 244–248 (1990). [CrossRef]  

7. M. Eichenfield, R. M. Camacho, J. Chan, K. J. Vahala, and O. Painter, “A picogram- and nanometer-scale photonic crystal opto-mechanical cavity,” Nature 459, 550–555 (2009). [CrossRef]   [PubMed]  

8. M. Eichenfield, J. Chan, R. M. Camacho, K. J. Vahala, and O. Painter, “Optomechanical crystals,” Nature 462, 78–82 (2009). [CrossRef]   [PubMed]  

9. J. D. Cohen and S. M. Meenehan, “Optical coupling to nanoscale optomechanical cavities for near quantum-limited motion transduction,” Opt. Express 21, 11227 (2013). [CrossRef]   [PubMed]  

10. G. Anetsberger, E. Gavartin, O. Arcizet, Q. P. Unterreithmeier, E. M. Weig, M. L. Gorodetsky, J. P. Kotthaus, and T. J. Kippenberg, “Measuring nanomechanical motion with an imprecision below the standard quantum limit,” Phys. Rev. A 82, 061804 (2010). [CrossRef]  

11. J. Chan, T. P. M. Alegre, A. H. Safavi-Naeini, J. T. Hill, A. Krause, S. Groblacher, M. Aspelmeyer, and O. Painter, “Laser cooling of a nanomechanical oscillator into its quantum ground state,” Nature 478, 89–92 (2011). [CrossRef]   [PubMed]  

12. K. Stannigel, P. Rabl, A. S. Sørensen, P. Zoller, and M. D. Lukin, “Optomechanical Transducers for Long-Distance Quantum Communication,” Phys. Rev. Lett. 105, 220501 (2010). [CrossRef]  

13. A. H. Safavi-Naeini and O. Painter, “Proposal for an optomechanical traveling wave phononphoton translator,” New J. Phys. 13, 013017 (2011). [CrossRef]  

14. C. A. Regal and K. W. Lehnert, “From cavity electromechanics to cavity optomechanics,” J. Phys.: Conf. Ser. 264, 012025 (2011).

15. S. Barzanjeh, M. Abdi, G. J. Milburn, P. Tombesi, and D. Vitali, “Reversible Optical-to-Microwave Quantum Interface,” Phys. Rev. Lett. 109, 130503 (2012). [CrossRef]   [PubMed]  

16. Y.-D. Wang and A. A. Clerk, “Using dark modes for high-fidelity optomechanical quantum state transfer,” New J. Phys. 14, 105010 (2012). [CrossRef]  

17. T. Bagci, A. Simonsen, S. Schmid, L. G. Villanueva, J. Zeuthen, E. Appel, J. M. Taylor, A. Srensen, K. Usami, A. Sorrenson, and E. S. Polzik, “Optical detection of radio waves through a nanomechanical transducer,” Nature 507, 81–85 (2014). [CrossRef]   [PubMed]  

18. R. W. Andrews, R. W. Peterson, T. P. Purdy, K. Cicak, R. W. Simmonds, C. A. Regal, and K. W. Lehnert, “Bidirectional and efficient conversion between microwave and optical light,” Nat. Phys. 10, 321–326 (2014). [CrossRef]  

19. J. Bochmann, A. Vainsencher, D. D. Awschalom, and A. N. Cleland, “Nanomechanical coupling between microwave and optical photons,” Nat. Phys. 9, 712–716 (2013). [CrossRef]  

20. M. Winger, T. D. Blasius, T. P. M. Alegre, A. H. Safavi-Naeini, S. Meenehan, J. Cohen, S. Stobbe, and O. Painter, “A chip-scale integrated cavity-electro-optomechanics platform,” Opt. Express 19, 24905–24921 (2011). [CrossRef]  

21. K. Cicak, M. Allman, J. Strong, K. Osborn, and R. Simmonds, “Vacuum-Gap Capacitors for Low-Loss Super-conducting Resonant Circuits,” IEEE Trans. Appl. Supercond. 19, 948–952 (2009). [CrossRef]  

22. J. Sulkko, M. A. Sillanpaa, P. Hakkinen, L. Lechner, M. Helle, A. Fefferman, J. Parpia, and P. J. Hakonen, “Strong Gate Coupling of High-Q Nanomechanical Resonators,” Nano Lett. 10, 4884–4889 (2010). [CrossRef]   [PubMed]  

23. J. Tian, W. Yan, Y. Liu, J. Luo, D. Zhang, Z. Li, and M. Qiu, “Optical quality improvement of Si photonic devices fabricated by focused-ion-beam milling,” J. Lightwave Technol. 27, 4306–4310 (2009). [CrossRef]  

24. A. H. Safavi-Naeini, T. P. M. Alegre, M. Winger, and O. Painter, “Optomechanics in an ultrahigh-Q two-dimensional photonic crystal cavity,” Appl. Phys. Lett. 97, 181106 (2010). [CrossRef]  

25. C. P. Michael, M. Borselli, T. J. Johnson, C. Chrystal, and O. Painter, “An optical fiber-taper probe for wafer-scale microphotonic device characterization,” Opt. Express 15, 4745–4752 (2007). [CrossRef]   [PubMed]  

26. S. S. Verbridge, H. G. Craighead, and J. M. Parpia, “A megahertz nanomechanical resonator with room temperature quality factor over a million,” Appl. Phys. Lett. 92, 013112 (2008). [CrossRef]  

27. M. Borselli, High-Q microresonators as lasing elements for silicon photonics, Ph.D. thesis, California Institute of Technology (2006).

28. S. G. Johnson, M. Ibanescu, M. A. Skorobogatiy, O. Weisberg, J. D. Joannopoulos, and Y. Fink, “Perturbation theory and Maxwell’s equations with shifting material boundaries,” Phys. Rev. E 65, 066611 (2002). [CrossRef]  

29. D. Rugar and P. Grütter, “Mechanical parametric amplification and thermomechanical noise squeezing,” Phys. Rev. Lett. 67, 699–702 (1991). [CrossRef]   [PubMed]  

30. R. Almog, S. Zaitsev, O. Shtempluck, and E. Buks, “Noise Squeezing in a Nanomechanical Duffing Resonator,” Phys. Rev. Lett. 98, 078103 (2007). [CrossRef]   [PubMed]  

31. M. Poot and H. X. Tang, “Phase shifting, strong coupling, and parametric feedback squeezing of an opto-electromechanical device,” in Proceedings of CLEO: 2013 p. CW3F.7 (2013).

32. C. C. Wu and Z. Zhong, “Capacitive Spring Softening in Single-Walled Carbon Nanotube Nanoelectromechanical Resonators,” Nano Lett. 11, 1448–1451 (2011). [CrossRef]   [PubMed]  

33. A. Szorkovszky, A. A. Clerk, A. C. Doherty, and W. P. Bowen, “Detuned mechanical parametric amplification as a quantum non-demolition measurement,” New J. Phys. 16, 043023 (2014). [CrossRef]  

34. I. Mahboob and H. Yamaguchi, “Bit storage and bit flip operations in an electromechanical oscillator,” Nat. Nano 3, 275–279 (2008). [CrossRef]  

35. S. Rips and M. J. Hartmann, “Quantum Information Processing with Nanomechanical Qubits,” Phys. Rev. Lett. 110, 120503 (2013). [CrossRef]   [PubMed]  

36. Q. P. Unterreithmeier, E. M. Weig, and J. P. Kotthaus, “Universal transduction scheme for nanomechanical systems based on dielectric forces,” Nature 458, 1001–1004 (2009). [CrossRef]   [PubMed]  

37. M. Agarwal, K.-K. Park, R. Candler, B. Kim, M. Hopcroft, S. A. Chandorkar, C. Jha, R. Melamud, T. Kenny, and B. Murmann, “Nonlinear Characterization of Electrostatic MEMS Resonators,” in International Frequency Control Symposium and Exposition, 2006 IEEE pp. 209–212 (2006). [CrossRef]  

38. V. Kaajakari, T. Mattila, A. Lipsanen, and A. Oja, “Nonlinear mechanical effects in silicon longitudinal mode beam resonators,” Sens. Actua. A 120, 64–70 (2005). [CrossRef]  

39. R. Lifshitz and M. C. Cross, “Response of parametrically driven nonlinear coupled oscillators with application to micromechanical and nanomechanical resonator arrays,” Phys. Rev. B 67, 134302 (2003). [CrossRef]  

40. S. M. Girvin, M. H. Devoret, and R. J. Schoelkopf, “Circuit QED and engineering charge-based superconducting qubits,” Phys. Scr. T137, 014012 (2009). [CrossRef]  

41. K. Cicak, D. Li, J. A. Strong, M. S. Allman, F. Altomare, A. J. Sirois, J. D. Whittaker, J. D. Teufel, and R. W. Simmonds, “Low-loss superconducting resonant circuits using vacuum-gap-based microwave components,” App. Phys. Lett. 96, 093502 (2010). [CrossRef]  

42. J. D. Teufel, T. Donner, D. Li, J. W. Harlow, M. S. Allman, K. Cicak, A. J. Sirois, J. D. Whittaker, K. W. Lehnert, and R. W. Simmonds, “Sideband cooling of micromechanical motion to the quantum ground state,” Nature 475, 359–363 (2011). [CrossRef]   [PubMed]  

43. S. M. Meenehan, J. D. Cohen, S. Groeblacher, J. T. Hill, A. H. Safavi-Naeini, M. Aspelmeyer, and O. Painter, “Thermalization properties at mK temperatures of a nanoscale optomechanical resonator with acoustic-bandgap shield,” http://arxiv.org/abs/1403.3703 (2014).

44. A. H. Safavi-Naeini, S. Groeblacher, J. T. Hill, J. Chan, M. Aspelmeyer, and O. Painter, “Squeezed light from a silicon micromechanical resonator,” Nature 500, 185–189 (2013). [CrossRef]   [PubMed]  

45. A. D. O’Connell, M. Ansmann, R. C. Bialczak, M. Hofheinz, N. Katz, E. Lucero, C. McKenney, M. Neeley, H. Wang, E. M. Weig, A. N. Cleland, and J. M. Martinis, “Microwave dielectric loss at single photon energies and millikelvin temperatures,” Appl. Phys. Lett. 92, 112903 (2008). [CrossRef]  

46. Y. Shim, J.-P. Raskin, C. R. Neve, and M. Rais-Zadeh, “RF MEMS Passives on High-Resistivity Silicon Substrates,” IEEE Microw. Wirel. Compon. Lett. 23, 632–634 (2013). [CrossRef]  

47. J. Chan, A. H. Safavi-Naeini, J. T. Hill, S. Meenehan, and O. Painter, “Optimized optomechanical crystal cavity with acoustic radiation shield,” Appl. Phys. Lett. 101, 081115 (2012). [CrossRef]  

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (4)

Fig. 1
Fig. 1 (a) Simulated Ey(r) electric field component of the optical mode. (b) Total deformation |Q(r)| of the fundamental planar flexural mode of the Si slab. (c) Scanning electron microscope (SEM) image of the top-view of a typical realized device. The scale bar is 5 μm. The zoomed-in view to the right depicts the electrode gap, whose size is d < 30 nm. (d) Illustration of the inverse shadow mask technique. A mask is created with ZEP resist to define the electrical contacts, while the capacitor electrodes are defined using the etched Si edges. After lift-off and full release, the metal within the gap is removed together with the SiO2 buffer layer. The scale bar is 5 μm. The purpose of the small trench on the top of the picture is to redistribute the residual silicon stress in order to avoid out-of-plane bending of the device.
Fig. 2
Fig. 2 (a) Experimental setup for sample characterization. The tunable diode laser is calibrated with a high-finesse Mach Zender interferometer (MZ) [27] and can be locked using a wavemeter (WM) for feedback. After being attenuated (VOA) and polarization controlled (FPC), the evanescent field of a tapered optical fiber is used to couple light into the photonic crystal cavity (PCC). The optical transmission is then measured with a low bandwidth photodetector (D1) in the case of optical spectroscopy, or via a high-speed photodetector (D2) connected to a Spectrum Analyzer (SA) for analysis of the mechanical spectrum. In addition, voltage sources (VDC), a signal generator (VAC) or a Network analyzer (NA) were used to measure the electric response of the device. (b) Measured normalized optical transmission through the photonic crystal of device A. The red line is the best Lorentzian fit of the resonance, showing an intrinsic optical Q-factor of Qc ∼ 9 × 104. (c) Power spectral density of the optically transduced thermal Brownian motion of the structure of device A. The strong resonance peak at ωm/2π = 63 MHz is the fundamental in-plane mechanical mode of the slab structure, whereas the smaller peaks represent weakly coupled (predominantly) out-of-plane modes of the slab. The red line is a multi-Lorentzian fit to the transduced mechanical spectrum. (d) Simulated capacitance as a function of gap size, d, for the capacitor on one side of the slab structure. The data fits with a power law (red line) Cm = [6.933/d0.533] fF, with d in units of nanometers. Inset: simulated Ey electrostatic field amplitude within the capacitor. The field isolines are shown in black.
Fig. 3
Fig. 3 (a) Waterfall plot of the optical transmission spectra for the fundamental optical cavity resonance as a function of DC voltage (VDC = 0 to 5.6 V, 0.2 V steps). The x-axis has been normalized according to the optical center wavelength at zero applied bias, λc,0. (b) Measured optical resonance wavelength (open circles) versus V DC 2. Solid black curve corresponds to a linear fit to the measured data. (c) Linear coupling parameter extracted from different experiments (DC (circles), AC (triangles)) and for different samples (A (red) and B (blue)). The solid black line is the numerically simulated linear electromechanical coupling using eq. (1). The dashed green line is the linear electromechanical coupling calculated using the functional dependence of the capacitance versus gap d shown in Fig. 2(d).
Fig. 4
Fig. 4 Nonlinear electromechanical response of device B. (a) Measured amplitude response (|S21|) for a weak AC voltage probe (V0 <1 mVrms) and DC bias ranging from 0.2 V (blue curve) to 6 V (red curve) in 0.2 V steps. Successive spectra are offset from the VDC = 0.2 V spectrum in 1 dB increments. (b) Measured amplitude response (|S21|) of the optically transduced mechanical motion to a strong near-resonant electric drive (ωAC ∼ ωm) with fixed VDC=6 V and V0 ranging from 277–870 mVrms. Increasing the drive amplitude, the mechanical Lorentzian peak (red) starts to assume the form of a Duffing oscillator (green). After an instability threshold, the mechanics enters a frequency entrainment regime in which it is completely in phase with the AC drive (blue). (c) Measured relative phase response (ϕ(S21)) of the optically transduced mechanical motion for the same drive conditions as in (b). Note that the phase slope is about the same inside and outside the frequency entrainment region indicating that the phase of the mechanical response is locked to the phase of the drive. Successive spectra in (b) and (c) are offset from the VAC = 277 mVrms spectrum in 1 dB increments.

Equations (8)

Equations on this page are rendered with MathJax. Learn more.

δ U ¯ E = V ( Q ˜ ( r ) n ) ( Δ ε | E | 2 Δ ε 1 | D | 2 ) d 2 r V ε | E | 2 d 3 r ,
F cap = 1 2 C m V 2 G e , 1 ,
m eff = Q * ( r ) ρ ( r ) Q ( r ) d 3 r max ( | Q | 2 ) ,
Δ ω c , DC ( V ) = g OM C m G e , 1 2 k eff V DC 2 = α V DC 2 ,
Δ ω c , AC = α Q m V 0 2 .
Δ ω m = ( ω m C m G e , 2 4 k eff ) V 2 = β V 2 ,
F cap V 2 = V DC 2 + V AC 2 + 2 V DC V AC .
V 0 , th = 2 k eff Q m 1 / 2 C m G e , 2 V DC .
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.