Mechanical on-chip microwave circulator.

Nonreciprocal circuit elements form an integral part of modern measurement and communication systems. Mathematically they require breaking of time-reversal symmetry, typically achieved using magnetic materials and more recently using the quantum Hall effect, parametric permittivity modulation or Josephson nonlinearities. Here we demonstrate an on-chip magnetic-free circulator based on reservoir-engineered electromechanic interactions. Directional circulation is achieved with controlled phase-sensitive interference of six distinct electro-mechanical signal conversion paths. The presented circulator is compact, its silicon-on-insulator platform is compatible with both superconducting qubits and silicon photonics, and its noise performance is close to the quantum limit. With a high dynamic range, a tunable bandwidth of up to 30 MHz and an in situ reconfigurability as beam splitter or wavelength converter, it could pave the way for superconducting qubit processors with multiplexed on-chip signal processing and readout.Nonreciprocal optical elements often require magnetic materials in order to break time-reversal symmetry. Here, Barzanjeh et al. demonstrate a magnetic-free on-chip microwave circulator that utilizes the interference from six electro-mechanical signal paths.

Introduction

Nonreciprocal devices are quintessential tools to suppress spurious modes, interferences and unwanted signal paths[1]. More generally, circulators can be used to realize chiral networks[2] in systems where directional matter-light coupling is not easily accessible. In circuit quantum electrodynamics[3] circulators are used for single-port coupling or as isolators to protect the vulnerable cavity and qubit states from electromagnetic noise and strong parametric amplifier drive tones. State-of-the-art passive microwave circulators are based on magneto-optic effects that require sizable magnetic fields[4,5], incompatible with ultra-low loss superconducting circuits. Due to the design principle their size is at least on the order of the wavelength and during manufacturing they need to be tuned and optimized one by one. Commercial circulators can therefore not be integrated on-chip causing additional losses and forming a major roadblock towards a fully integrated quantum processor based on superconducting qubits.

Many recent theoretical and experimental efforts have been devoted to overcome these limitations both in the optical[6-8] and microwave regimes[9-18]. In parallel, the rapidly growing field of optomechanical and electromechanical systems has shown promising potential for applications in quantum information processing and communication, in particular for microwave to optical conversion[19,20] and amplification[21]. Very recently, several theoretical proposals[22-24] have pointed out that reservoir-engineered optomechanical systems[25] can lead to nonreciprocity and first isolators have just been demonstrated in the optical domain[26-28].

Here we present an on-chip microwave circulator using a frequency tuneable silicon-on-insulator electromechanical system[29] that is compatible with superconducting qubits[30]. The device can be reconfigured in situ as a filter, splitter, mixer, isolator or circulator. For the first mechanically mediated microwave circulator, we achieve an isolation of up to 24–38 dB, compared to total losses of 4.7–8.5 dB and at most 4–7 added noise quanta over an instantaneous bandwidth of 630 Hz.

Results

Device characterization and tunability

The main elements of the microchip circulator device are shown in Fig. 1a, b. The circuit comprises three high-impedance spiral inductors (L i) capacitively coupled to the in-plane vibrational modes of a dielectric nanostring mechanical resonator. The nanostring oscillator consists of two thin silicon beams that are connected by two symmetric tethers and fabricated from a high resistivity silicon-on-insulator device layer[29]. Four aluminium electrodes are aligned and evaporated on top of the two nanostrings, forming one half of the vacuum gap capacitors that are coupled to three microwave resonators and one DC voltage bias line as shown schematically in Fig. 1c (see Supplementary Tables 1 and 2 for details). The device is mounted on the mixing chamber plate of a cryogen-free dilution refrigerator at a temperature of T f = 10 mK and all incoming lines are strongly filtered and attenuated to suppress Johnson and phase noises.

Fig. 1

Microchip circulator and tunability. a Scanning electron micrograph of the electromechanical device including three microwave resonators, two physical ports labelled 1 and 2, one voltage bias input (V dc) on the top left. The inset shows the spiral inductor cross-overs in the green dashed boxed area. The total device area is 0.3 mm by 0.45 mm. Scale bar (yellow) corresponds to 100 μm. b Enlarged view of the silicon nanostring mechanical oscillator with four vacuum-gap capacitors coupled to the three coil inductors and one voltage bias input. Insets show details of the nanobeam as indicated by the dashed and dotted rectangles. Scale bar (yellow) corresponds to 10 μm. c Electrode design and electrical circuit diagram of the device. The input modes a couple inductively to the microwave resonators with inductances L , coil capacitances C , additional stray capacitances C and the motional capacitances C . The reflected tones a pass through a separate chain of amplifiers each, and are measured at room temperature using a phase-locked spectrum analyzer (not shown). The simulated displacement of the lowest frequency in-plane flexural modes of the nanostring are shown in the two insets. Colour indicates relative displacement. d Resonator reflection measurement of the three microwave resonators of an identical device, as a function of the applied bias voltage and a fit (dashed lines) to Δω = α 1 V 2 + α 2 V 4 with the tunabilties α 1/2π = 0.53 MHz/V2 and α 2/2π = 0.05 MHz/V4 with a total tunable bandwidth of 30 MHz for resonator 2 at 9.8 GHz

The voltage bias line can be used to generate an attractive force that pulls the nanobeam and tunes the operating point frequencies of the device[31]. Figure 1d shows the measured resonance frequency change as a function of the applied bias voltage V dc. As expected, resonators 1 and 3 are tuned to higher frequency due to an increased vacuum gap, whereas resonator 2 is tuned to lower frequency. At a large tunable bandwidth of up to 30 MHz, as obtained for resonator 2, the ability to excite the motion directly and to modulate the electromechanical coupling in situ represents an important step towards new optomechanical experiments and more practicable on-chip reciprocal and nonreciprocal devices.

As a first step, we carefully calibrate and characterize the individual electromechanical couplings and noise properties, similar to ref. [32]. The measured thermalization temperature of the two mechanical modes are 18 and 25 mK and the final occupancy for the reported device reached as low as 0.6 and 2 quanta for standard motional sideband cooling. We then measure the bidirectional frequency conversion between two microwave resonator modes as mediated by one mechanical mode[33]. The incoming signal photons can also be distributed to two ports with varying probability as a function of the parametric drive strength and in direct analogy to a tunable beam splitter. We present the relevant sample parameters in Supplementary Tables 1 and 2, the theoretical analysis in Supplementary Note 2 and experimental results of this bidirectional frequency conversion process in Supplementary Fig. 1.

Theoretical model

Directionality is achieved by engaging the second mechanical mode, a method that was developed in parallel to this work[34,35] for demonstrating nonreciprocity in single-port electromechanical systems. In fact, creating a parametric coupling between the two electromagnetic and the two mechanical modes by four microwave pumps with frequencies slightly detuned from the lower motional sidebands of the resonances, creates two paths for exchanging photons, as shown in Fig. 2a. Nonreciprocity appears when these paths interfere destructively, leading to breaking the symmetry between the two directions. For a detailed description, we begin with the theoretical model describing two microwave cavities with resonance frequencies ω and total linewidths κ with i = 1, 2 parametrically coupled to two distinct modes of a mechanical resonator with resonance frequencies ω and damping rates γ with j = 1, 2. To establish the parametric coupling, we apply four microwave tones, with frequencies detuned by δj from the lower motional sidebands of the resonances, as shown in Fig. 2a. In a reference frame rotating at the frequencies ω and ω  + δj, the linearized Hamiltonian in the resolved sideband regime (ω     κ 1, κ 2) is given by (ħ = 1)where a (b ) is the annihilation operator for the cavity i (mechanics j), and g 0ij are the effective and vacuum electromechanical coupling rates between the mechanical mode j and cavity i, respectively, while n is the total number of photons inside the cavity i due to the drive with detuning Δ, and ϕ is the relative phase set by drives. Here, Δ11 = Δ21 = ω  + δ 1 and Δ22 = Δ12 = ω  + δ 2 are the detunings of the drive tones with respect to the cavities and H off describes the time-dependent coupling of the mechanical modes to the cavity fields due to the off-resonant drive tones. These additional coupling terms create cross-damping[36] and renormalize the mechanical modes, and can only be neglected in the weak coupling regime for .

Fig. 2

Electromechanical isolator. a Mode coupling diagram for electromechanically induced nonreciprocity. Two microwave cavities (C 1 and C 2) are coupled to two mechanical modes (M 1 and M 2) with the electromechanical coupling rates G (where i, j = 1, 2), inducing two distinct signal conversion paths. Power spectral density (PSD) of the two microwave cavities and arrows indicating the frequency of the four microwave pump tones slightly detuned by δ from the lower motional sidebands of the resonances. All four pumps are phase-locked while the signal tone is applied. Only one of the microwave source phases is varied to find the optimal interference condition for directional transmission between ports 1 and 2. b Measured power transmission (dots) in forward |S 21|2 (cavity 1 → cavity 2) and backward directions |S 12|2 (cavity 2 → cavity 1) as a function of probe detuning for two different phases ϕ = ±102.6º. The solid lines show the results of the coupled-mode theory model discussed in the text. Grey shaded areas denote cavity loss and green shaded areas insertion loss. c Experimental data and theoretical model of measured transmission coefficients |S 12|2 and |S 21|2 as a function of signal detuning and pump phase ϕ. Dashed lines indicate the line plots for the two phases ϕ = ±102.6º as shown in b

To see how the nonreciprocity arises, we use the quantum Langevin equations of motion along with the input–output theorem to express the scattering matrix S of the system described by the Hamiltonian (1), and relating the input photons a in,(ω ) at port i to the output photons a out,(ω ) at port j via with i = 1, 2. The dynamics of the four-mode system described by Hamiltonian (1) is fully captured by a set of linear equations of motion as verified in Supplementary Notes 3–6. Solving these equations in the frequency domain, using the input–output relations, and setting ϕ 22 = ϕ, ϕ 11 = ϕ 21 = ϕ 12 = 0, the ratio of backward to forward transmission readsHere, Σ = 1 + 2i[(−1) δ − ω]/γ is the inverse of the mechanical susceptibility divided by the mechanical linewidth γ and is the electromechanical cooperativity. Note that, in Eq. (2) we assume the device satisfies the impedance-matching condition on resonance, i.e., S (ω = 0) = 0, which can be achieved in the high-cooperativity limit ().

Inspection of Eq. (2) reveals the crucial role of the relative phase between the drive tones ϕ and the detuning δ to obtain nonreciprocal transmission. When the cooperativities for all four electromechanical couplings are equal () then perfect isolation, i.e. λ = 0, occurs forEquation (3) shows that on resonance (ω = 0) tan[ϕ]   δ, highlighting the importance of the detuning δ to obtain nonreciprocity. Tuning all four drives to the exact red sideband frequencies (δ = 0) results in bidirectional behaviour (λ = 1). At the optimum phase ϕ given by Eq. (3), ω = 0, and for two mechanical modes with identical decay rates (γ  = γ  = γ) the transmission in forward direction is given bywhere η 1(2) = κ ext,1(2)/κ 1(2) is the resonator coupling ratio and κ  = κ int, + κ ext, is the total damping rate. Here κ int, denotes the internal loss rate and κ ext, the loss rate due to the cavity to waveguide coupling. Equation (4) shows that the maximum of the transmission in forward direction, , occurs when and for large cooperativities . These conditions, as implemented in our experiment, enable the observation of asymmetric frequency conversion with strong isolation in the backward direction and small insertion loss in forward direction.

Bidirectional wavelength conversion

For bidirectional wavelength conversion, higher cooperativity enhances the bandwidth, as shown in Supplementary Note 2. In contrast, the bandwidth of the nonreciprocal conversion is independent of cooperativity and set only by the intrinsic mechanical linewidths γ , which can be seen in Eq. (2). This highlights the fact that the isolation appears when the entire signal energy is dissipated in the mechanical environment, a lossy bath that can be engineered effectively[25]. In the present case it is the off-resonant coupling between the resonators and the mechanical oscillator that modifies this bath. The applied drives create an effective interaction between the mechanical modes, where one mode acts as a reservoir for the other and vice versa. This changes both the damping rates and the eigenfrequencies of the mechanical modes.

It, therefore, increases the instantaneous bandwidth of the conversion and automatically introduces the needed detuning, which is fully taken into account in the theory.

Two-port microwave isolator

Using the on-chip electromechanical microwave circuit shown in Fig. 1a, we experimentally realize directional wavelength conversion between two superconducting coil resonators at (ω 1,ω 2)/2π = (9.55, 9.82) GHz coupled to two different physical waveguide ports and measurement lines with (η 1, η 2) = (0.74, 0.86). Here we use the two lowest-frequency vibrational in-plane modes of the mechanical resonator at (ω , ω )/2π = (4.34, 5.64) MHz with intrinsic damping rates (γ , γ )/2π = (4, 8) Hz. The vacuum electromechanical coupling strengths for these mode combinations are (g 0,11, g 0,12, g 0,21, g 0,22)/2π = (33, 34, 13, 31) Hz. The microwave resonators are driven with four coherent microwave sources with powers (P 11, P 12, P 21, P 22) = (−73.3, −68.7, −66.9, −67.4) dBm at the device inputs that correspond to the single cavity–single mechanical cooperativities (C 11, C 12, C 21, C 22) = (47, 43.8, 41.9, 56.9). Transmission parameters are measured by using a weak probe signal with a signal power of only −117 dBm at the device inputs.

Figure 2b shows the measured transmission of the wavelength conversion in the forward |S 21|2 and backward directions |S 12|2 as a function of probe detuning for two different phases as set by one out of the four phase-locked microwave drives. At ϕ = −102.6º and over a bandwidth of 518 Hz, we measure high transmission from cavity 1 to 2 with an insertion loss of 2.4 dB due to finite input matching and a resonator loss of 1.9 dB due to finite intrinsic resonator linewidths. In the backward direction, the transmission is suppressed by up to 40.4 dB. Likewise, at the positive phase of ϕ = 102.6º the transmission from cavity 1 to 2 is suppressed while the transmission from cavity 2 to 1 is high. In both cases, we observe excellent agreement with theory (solid lines). Figure 2c shows the S parameters for the whole range of phases ϕ, which are symmetric and bidirectional around ϕ = 0. We find excellent agreement with theory over the full range of measured phases with <10% deviation to independently calibrated drive photon numbers and without any other free parameters.

Extension to a microwave circulator

The described two-port isolator can be extended to an effective three-port device by parametrically coupling the third microwave resonator capacitively to the dielectric nanostring, as shown in Fig. 1a. The third resonator at a resonance frequency of ω 3/2π = 11.30 GHz is coupled to the waveguide with η 3 = 0.52 and to the two in-plane mechanical modes with (g 0,31, g 0,32)/2π = (22, 45) Hz. Similar to the isolator, we establish a parametric coupling between cavity and mechanical modes using six microwave pumps with frequencies slightly detuned from the lower motional sidebands of the resonances, which for certain pump phase combinations can operate as a three-port circulator for microwave photons, see Fig. 3a. The pump powers at the device inputs are (P 11, P 12, P 21, P 22, P 31, P 32) = (−72.5, −69, −67.5, −68, −69, −70) corresponding to the single resonator cooperativities (C 11, C 12, C 21, C 22, C 31, C 32) = (56.5, 40.9, 35, 49.6, 99.2, 49, 6). Using an additional microwave source as a weak probe signal with a signal power of only −117 dBm at the device inputs, we measure the power transmission between all ports and directions as shown in Fig. 3b for a single fixed phase of ϕ = −54º, optimized experimentally for forward circulation.

Fig. 3

Electromechanical circulator. a Mode coupling diagram describing the coupling between three microwave cavities (C 1, C 2 and C 3) and two mechanical modes (M 1 and M 2) with electromechanical coupling rates G (where i = 1, 2, 3 and j = 1, 2), creating a circulatory frequency conversion between the three cavity modes, as indicated by the grey arrow. b Measured power transmission in forward |S 21|2, |S 32|2 and |S 13|2 (red dots) and backward directions |S 12|2, |S 23|2 and |S 31|2 (blue dots) as a function of probe detuning for a pump phase ϕ = −54º. The solid lines show the prediction of the coupled-mode theory model discussed in the text. The inset shows the signal propagation between the three resonator modes and the black arrow indicates the circulator direction. Grey shaded areas denote cavity loss and green shaded areas insertion loss. c Measured S parameters and theoretical model as a function of detuning and pump phase. Dashed lines indicate the line plot positions shown in b

At this phase, we see high transmission in the forward direction S 21,32,13 with an insertion loss of (3.8, 3.8, 4.4) dB, due to imperfect input matching, and an isolation in the backward direction S 12,23,31 of up to (18.5, 23, 23) dB over a bandwidth of 628.5 Hz. The full dependence of the circulator scattering parameters on the drive phase is shown in Fig. 3c, where we see excellent agreement with theory. The added noise photon number of the device is found to be (n add,21, n add,32, n add,13) = (4, 6.5, 3.6) in the forward direction and (n add,12, n add,23, n add,31) = (4, 4, 5.5) in the backward direction, limited by the thermal occupation of the mechanical modes and discussed in more detail in Supplementary Fig. 3.

Discussion

In conclusion, we demonstrated a frequency tunable and in situ reconfigurable signal processing device that can act as a filter, wavelength converter, beam splitter, isolator or circulator for microwave photons. The circulator is highly directional and operates with relatively low loss and added noise. Improvements of the circuit design and fabrication will help to increase the instantaneous bandwidth and decrease the insertion losses. Compared to the so far significantly higher bandwidth Josephson devices, a mechanical approach is insensitive to magnetic field noise and offsets, mechanical devices have higher dynamical range due to the smaller nonlinearity, and well-confined mechanical modes are typically less prone to parasitic coupling when integrated in larger systems. In addition, mechanical systems have the potential for hybrid microwave and optical signal processing, in particular for non-reciprocity between microwave and optical propagating fields. The presented external voltage bias offers new ways to achieve directional amplification and squeezing of microwave fields in the near future.

Data availability

The data that support the findings of this study are available from the corresponding author upon reasonable request.