Non-Degenerate Multimode Optomechanics

Non-Degenerate Multimode Optomechanics

L. F. Buchmann Department of Physics, University of California, Berkeley, California 94720, USA    D. M. Stamper-Kurn Department of Physics, University of California, Berkeley, California 94720, USA Materials Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA

We theoretically investigate interactions between non-degenerate mechanical oscillators mediated by a time-dependent cavity field. We obtain a reduced master equation valid for all optomechanical systems operating in the weak coupling regime. This master equation includes all forms of decoherence and back-action due to the dissipation of the field mediating the interaction. We apply the master equation to study two resonant coupling schemes within a rotating-wave approximation: the beam splitter Hamiltonian and the two-mode parametric amplifier. In both cases, the effective unitary interaction can be made arbitrarily strong compared to the decoherence due to dissipation of the mediating field by choosing appropriate detunings.

42.50.Pq 03.65.-w 06.20.-f 43.50.+y

I Introduction

The ability to control a physical system at the quantum level has, whenever it was achieved, profoundly influenced physics and many of its applications. Two prominent examples, the development of lasers and their application to cool atomic gases, have allowed to use light as a tool to perform a wide range of tasks and lead to improvements of precision measurements by orders of magnitude. In the last few years, successes in cavity optomechanics have lead to similar control over motional degrees of freedom of mesoscopic and microscopic oscillators review ().

Experiments have successfully demonstrated quantum effects involving a single mechanical oscillator OMsqueezing (); Cleland (); ultracold (); sidebandasym (); forcesensitivity (); palomaki () and theoretical studies point out the suitability of optomechanical systems to observe truly macroscopic non-classical behavior nonclassicalOM (); MAQRO (); tunneling (); nonclassicalOM2 (); Yanbei (). While the vast majority of the experimental results involve a single mechanical mode, typical solid-state resonators exhibit a multitude of resonances over a wide range of frequencies regalmultimode (); mukundmultimode (). In systems where the mechanical modes consist of the collective motion of a trapped atomic cloud ultracold (); Esslinger () one can engineer non-degenerate optomechanical modes by trapping several clouds in an optical superlattice coldatomsmultimode (). Numerous non-degenerate mechanical modes are thus a common feature of optomechanical systems, which feature, however, remains largely unexploited to date as most experimental results multimodeOM5 (); multimodeOM6 (); multimodeOM3 (); regalmultimode () involve nearly degenerate mechanical modes.

While multimode optomechanics has been studied theoretically for specific applications mattash (); multimodeOM1 (); multimodeOM2 (); multimodeOM3 (); multimodeOM4 (); phasecon (); mediatedinteractions (); borkje (), there has been no general formalism to analyze the possibilities of tuning non-degenerate modes into resonance. The aim of the present work is to provide that formalism, including the effects of quantum noise and measurement back-action of the mediating field. The resulting equations are applicable over the whole range of optomechanical systems, from large scale interferometers for gravitational wave detection gravwave () to oscillations of ultracold atomic clouds ultracold () or optomechanical microcavities review (); microcavity ().

The effective interaction between non-degenerate mechanical modes is suppressed by their frequency difference. Together with the typically small optomechanical couplings this prohibits effective multimode physics involving modes with vastly different frequencies. However, other fields of physics successfully use auxiliary modes to resonantly couple non-degenerate states. In atomic physics, for instance, two laser beams appropriately detuned from a third state can drive a two-photon resonance between two non-degenerate atomic states textbook (). In non-linear optics, two optical fields can be tuned into resonance by an appropriate pump creating a standing wave in a non-linear crystal textbook (). The optomechanical analogue of these techniques is an oscillating intensity in a cavity field coupled to two non-degenerate mechanical modes and is the subject of this Article.

A time dependent intracavity intensity is created by any anharmonic pump, the simplest example being the field created by two harmonic drives with distinct frequencies. The resulting cavity field can mediate a resonant coupling irrespective of the frequency difference of the mechanical modes. Recent theoretical studies have applied specific examples of this very effect for specialized applications, such as the cooling of a low-lying mechanical mode borkje (), mechanical phase-conjugation phasecon () or preparation of non-classical states mattash (). The price for resonant mediated interaction between non-degenerate modes is additional decoherence due to the dissipative nature of the mediating field. The cavity field that provides mechanical coupling also carries information about the mechanical state out of the cavity and this loss of unitarity is necessarily accompanied by noise entering the mechanical system mediatedinteractions ().

In this Article, we derive a general and complete master equation for the behavior of two non-degenerate mechanical modes interacting via weak coupling to a shared cavity field. Our treatment is complete in the sense that all sources of decoherence, back-action and noise that are a consequence of the mediating field’s dissipation are included. Taking these effects into account is crucial for the operation of any optomechanical device at the quantum level and a prerequisite to the observation of non-classical effects. The resulting equation is general in the sense that we do not assume any specific detuning or rely on the system’s being in any particular parameter regime, apart from the weak coupling assumption which is accessible to virtually all optomechanical systems to date review (). The derivation of our master equation does not rely on the mediating field’s being the mode of an electromagnetic resonator. With a few simple modifications it will also describe interactions mediated by bulk-phononic modes mukundmed (). Similarly it is straightforward to add individual, Markovian heat baths for the mechanical oscillators and while we have neglected them for the discussion of the individual couplings, their effect is discussed in the final section.

After deriving the general master equation, we will study the two resonant effective interactions described by the beam splitter Hamiltonian


and the two-mode parametric amplifier


where denotes the annihilation operator for phonons in oscillator . The former allows the mapping of a quantum state from one oscillator to the other while the latter leads to non-classical correlations between the two oscillators, akin to parametric down-conversion in quantum optics. Most results obtained in single-mode optomechanics consist of these two Hamiltonians resonantly coupling the cavity-mode with the mechanical mode via the matching of mechanical side bands with the cavity resonance. Optomechanical cooling, for instance, is a realization of the beam splitter Hamiltonian together with the large dissipation of the electromagnetic field modes, while parametric driving is a realization of the parametric amplifier involving mechanical and electromagnetic modes palomaki (). Realizing these couplings between two mechanical modes opens the door to more customized and precise sensing of forces and fields, the application of quantum computation algorithms to mechanical systems computing () and the study of macroscopic decoherence diosi ().

Of crucial importance is the ratio between the strength of the mediated unitary interactions and the decoherence rates due to the mediating field being emitted from the cavity. This quantity is a measure for the ability of the effective interaction to create or maintain non-classical correlations between the two modes and thus a gauge for the classicality of the interactions mediatedinteractions (); kafri1 (); kafri2 (). We find that the relationship between unitary coupling strengths and total decoherence rates depends on the system parameters and the detunings of the pumps with respect to the cavity resonance, such that different choices of detunings turn out to be ideal for a given application. Most importantly, we find that even for systems where the cavity emission rate of the mediating field exceeds all of the involved mechanical frequencies there are regimes where the engineered interaction allows non-classical motional states or the transfer of nonclassical states between oscillators.

Ii Master Equation

ii.1 System Hamiltonian

Consider an optomechanical system consisting of two mechanical modes with natural frequencies and coupled to the field of an electromagnetic resonator with single photon/single phonon coupling strengths and respectively.A sketch of a possible setup involving spatially separate membranes coupled to the field of an optical Fabry-Pérot resonator is given in Fig. 1. The mechanical modes need not be located on distinct oscillators and the mediating field does not have to be optical. The system could also consist of two non-degenerate drumhead modes of a capacitor in a superconducting LC circuit palomaki () with photons in the microwave regime or two samples of ultracold atoms trapped in an optical cavity coldatomsmultimode ().

Figure 1: (Color online) (a) Sketch of a possible experimental realization involving separate optomechanical membranes. The mechanical modes could also be on the same resonator. (b) A scattering process contributing to the coherent exchange of phonons in the resonant beam splitter interaction. The three other contributing processes result from time reversal and pump photons with frequency scattering with oscillator 2. (c) A scattering process contributing to the parametric amplification of the mechanical modes. The three other processes are obtained in the same way as in panel (b).

Neglecting all other modes of the mechanical system and the cavity, the Hamiltonian in the lab frame is given by




and accounts for the intensity of the intracavity field being emitted to a reservoir at zero temperature at a rate 111The frequencies of electromechanical systems working in the microwave regime cannot be assumed to feature a vanishing thermal photon occupation. However, typical devices are operated cryogenically, where a zero temperature reservoir is reasonable.. We have set , dropped constant terms and defined the input coherent field driving the cavity with the rate . The coupling of each oscillator to its thermal reservoir is neglected in favor of a focus on decoherence due to the dissipative nature of the field that mediates the interaction. Thermal effects are explicitly discussed later.

To tune a desired coupling into resonance, we need the cavity intensity to oscillate in time. This is achieved by driving the cavity with a bi-chromatic pump,


The beat-note between the two drives will provide the necessary energy to tune either Hamiltonian (1) or (2) into resonance. The density operator of the system satisfies the master equation


with the Lindblad superoperator


Without loss of generality we will assume and . The next step is to unitarily displace the cavity field operator by a time-dependent complex function: . If we choose




where , all source terms for cancel and it only carries quantum fluctuations of the cavity field. Without loss of generality, we may choose the phases of such that both are real. The displacement operation introduces the term driving the mechanical degrees of freedom. The time independent part of this term is a constant force that will shift the equilibrium position of each oscillator. Additionally, the beat-note between the two pumps will drive each oscillator off-resonantly. The constant force does not affect the dynamics of the system. For the small mechanical linewidths typical of optomechanics review () the off-resonant terms drive oscillator into a coherent state with phonon number of the order of . For weak couplings his small oscillation can be neglected for all effects that are considered in the remainder of this paper.

Similarly, we can also neglect the interactions quadratic in and thus find in a co-rotating frame for the optical operators the Hamiltonian


ii.2 Tracing out the cavity field

Tracing out the electromagnetic system will yield an effective master equation coupling the mechanical modes. To perform the partial trace over electromagnetic degrees of freedom we restrict the optical Hilbert space to the optical vacuum and first-order coherences. This is formally a projection of the optical Hilbert space unto a lower-dimensional subspace as described below. Introducing Fock states of the optical subsystem and operators


the reduction is formalized as


By virtue of the displacement operation this is a reasonable approximation. The reduced density operator of the mechanical system becomes


Note, however, that this identification can only be made after taking the trace over optical degrees of freedom in Eq. (8), as there are non-vanishing contributions from the one-photon sector. Taking this into consideration, we find the equation of motion


with denoting the Hermitian conjugate. The evolution of the operator is governed by


Assuming that changes slowly on a timescale we can solve this equation by elementary transformations. The result is


This expression can be plugged into Eq. (II.2) to give the reduced master equation for the mechanical system. Introducing the vectors , where denotes the transpose, we obtain


Here are matrices




and .

This is the general form of the master equation in the lab frame, including all off-resonant terms. It is easily reduced to the system of a single harmonic oscillator (with and or ) or a single harmonic pump (with and or ). As sketched in Fig. 1 (b) and (c), mediated interactions involve at least two photon-phonon scattering processes. Thus all effective amplitudes appearing in Eqs. (21) are of second order in the single phonon/single photon coupling strength.

The unitary evolution described by Eq. (19) is conceptually simple. Each pump shifts the resonance of the oscillators by the frequency


where we have separated the time-independent spring shifts from time-dependent contributions easily found from Eqs. (21). The explicit form of the spring shifts is


The mechanical resonances oscillate with small amplitudes around their means . Due to the weakness of optomechanical couplings we have and the oscillation of their resonance frequency will not affect the dynamics of the system for sufficiently non-degenerate modes, quantified by . Additionally, the light-field also establishes an optical spring between the oscillators, with spring constant given by


where is the zero-point displacement of oscillator . The re-scaling by the two oscillator lengths serves to give the units of force per length but does not appear in the equations of motion, where this length has been absorbed into the definition of the operators . This spring too consists of constant and time-dependent contributions. In contrast to the single-mode optical springs given in Eq. (22), it is now the time-dependent contributions that can create resonant couplings between the non-degenerate oscillators.

Iii Resonant Couplings

In the following section we will focus on the cases where the time-dependence of the optical spring between the oscillator resonates with eigenmodes of the two-oscillator system. We will focus on a regime where the frequency difference between the two oscillators is large compared to the effective coupling strengths,


For such systems, a monochromatic pump does not lead to multimode behavior and the system can be treated as two separate single-mode optomechanical oscillators. In our case of two harmonic pumps, the beat-note between the two drives can tune a desired interaction into resonance. For such cases condition (25) assures off-resonant contributions to be suppressed such that they can be neglected within a rotating-wave approximation.

The beat-note between the two optical pumps has to hit the resonances of the optically shifted mechanical frequencies . This turns the resonance conditions into non-linear equations with respect to the optical detunings. The correction of the modified resonance conditions compared to their bare values are small, but we have computed them for completeness and they are given in the appendix. From here on, we assume that these corrections are taken into account and take them into the definitions of the mechanical frequencies, i.e.


This substitution has to be carried through in Eqs. (21) to be self consistent and does not lead to qualitatively different dynamics in the weak-coupling regime.

Resonant coupling arises when the modulation frequency of the spring connecting the two oscillators matches the difference of the two individual oscillator resonances,


where we have introduced for the difference of the mechanical frequencies. The resulting interaction – governed by Hamiltonian (1) – can be used to exchange the quantum states of the two mechanical modes irrespective of their frequency difference. A diagram for one particular set of scattering processes contributing to this interaction is given in Fig. 1 (b). The two classical pumps driving the cavity establish the conditions for resonant exchange of phonons between the mechanical modes via scattering of intracavity photons. Alternatively, the modulation of the coupling spring can match the sum of the oscillator frequencies,


with being the average frequency of the bare mechanical modes. In this case the resulting interaction is a parametric drive of a collective mode of the two-oscillator system described by Hamiltonian (2). An example for a process contributing coherently to this effective interaction is given in Fig. 1 (c).

In the following, we will study the two resonant cases in detail. To this end, it is useful to transform into a co-rotating frame for the mechanical operators, achieved by applying the unitary


to the master equation (19). The transformed Hamiltonian does not contain the free mechanical Hamiltonians , but the frequencies now modulate the off-diagonal elements of the matrices appearing in the master equation,


iii.1 Beam Splitter

In this subsection we will study the setting where the effective interaction is given by the beam splitter Hamiltonian (1). This interaction leads to a coherent exchange of quantum states between the two mechanical modes at a rate . Applications include the cooling of a low-frequency mode by means of an auxiliary high-frequency mode, the study of mediated interactions mediatedinteractions () or quantum state mapping. In the time-domain the resonant coupling corresponds to a modulation of the strength of a spring connecting the two oscillators. It is instructive to discuss this condition in the frequency domain.

Figure 2: (Color online) Sketch of the intracavity photon power spectral density neglecting the cavity profile for clarity. (a) and (b) show the contributions of the two pumps in red (solid) and blue (dashed) respectively. The difference between the pump frequencies is arranged in the beam splitter configuration specified by Eq. (27) with the pumps indicated by arrows representing delta-peaks at and . Each pump creates a pair of sidebands for each mechanical oscillator, as indicated in the sketch. The overlapping sidebands are emphasized by the dashed lines and formally captured by Eqs. (33). The exact behavior of this quantity including the cavity profile can be obtained from Eq. (18)

The optomechanical interaction leads to side bands on the cavity field, identifiable in Eq. (18) and associated with photon-phonon scattering processes. Each pump creates two pairs of side bands, one for each oscillator. Their optical frequencies are


The resonance condition (27) implies


where we have introduced the central detuning


These identities represent overlapping side bands of different oscillators, as sketched in Fig. 2. When a photon coming from one pump mode leaves the cavity in the other pump mode, it will have interacted with both oscillators and mediated the coherent transfer of a phonon between oscillators. The distinct overlapping side bands in Eqs. (33) provide two distinguishable pathways for phonon exchange and the coherent effective coupling rate is the sum of the amplitudes for each pathway. Additionally, there are four unmatched side bands with frequencies


These unmatched side bands contain photons which have scattered phonons off a single mechanical mode. They carry information about the oscillators out of the cavity and are sources of dissipation and back-action.

The resonant coupling terms are found in the matrices ; using identities (33) we find


where the star denotes complex conjugation. The coherent processes are contained in the imaginary parts of these coupling terms,


where the two terms correspond to amplitudes for the scattering being mediated by the matched side bands on the red and blue side of the pumps respectively. The real parts of Eq. (36) together with resonant terms in describe the dissipative dynamics. We find


Explicitly separating unitary and dissipative terms yields after some algebra the master equation


with the effective Hamiltonian


and the collective mode


where we have assumed and introduced for notational brevity. The chosen normalization of the collective mode yields non-canonical commutation relations but simplifies the coefficients appearing in the master equation. They are


The master equation (39) describes two resonantly coupled bosonic modes connected to three independent heat baths. The two heat baths described by coupling rates and effective thermal occupations affect each oscillator separately, while the bath described by with effective occupation acts on a collective mode of the system.

iii.1.1 Unitary Evolution

For given side band resolution and optomechanical couplings, the maximum exchange rate between the two states is reached for . For typical experimental parameters the coupling is around  review () and can be tuned via the pump strengths. The coupling strength as a function of central detuning changes its behavior depending on the side band resolution as plotted in Fig. 3.

Figure 3: (Color online) Normalized coupling strength as a function of central detuning and side band resolution. The line gives parameters resulting in zero coupling, . Other parameters are .

For there are two non-zero values of resulting in no effective coupling, indicated by a line in the figures. This cancellation is a temporal analogue of a double-slit interference. The processes of phonon exchange via either matched side band (see Fig. 2 and Eqs. (33)) can coherently cancel. This is, however, only possible if the photon coherence time exceeds the inverse separation of the two processes .

iii.1.2 Dissipative Dynamics

The single-mode heat baths, characterized by Eqs. (42b),(42c),(42e) and (42f), correspond to the unmatched side bands of the two pumps, see Fig. 2. These photons only carry information about the positions of single oscillators and affect the mechanical modes separately. Photons in the overlapping side bands partially contribute to the coherent exchange of excitations through , but also leak out of the cavity, affecting the collective mode . Effective temperatures are defined by the occupations and . Parameter regimes where or become negative are to be understood in the sense of statistical mechanics and indicate anti-damping of the mode in question.

For pumps that are arranged symmetrically around the cavity resonance () there is no effective spring coupling the two oscillators, . The light field only measures the mechanical system and its dissipation is pure measurement back-action of the light field on the mechanical oscillators 222The effective temperature diverges, but the quantity appearing in the master equation remains finite. Cavities with a dissipation rate smaller than the frequency difference of the oscillators can resolve each heat bath independently, as plotted in Fig. 4. In such systems it is possible to cool or amplify each mode or the collective mode selectively and it is sufficient to consider only the dominant heat bath.

iii.1.3 Limits to Coherence

The detuning resulting in the strongest coupling strength is not necessarily the ideal choice to couple two mechanical modes. For coherent manipulations on a quantum level it is necessary to minimize the influence of decoherence. The total decoherence rate of the system due to the dissipation of the mediating field is given by


and is plotted in Fig. 4 as a function of and .

Figure 4: (Color online) Scaled total decoherence rate as a function of detuning and optical dissipation rate. Other parameters are the same as in Figs. 3.

Each contribution is a Lorentzian with width centered on the blue detuned side of the cavity resonance at and . For large , the ratio is typically smaller than unity for in the vicinity of and thus one would have to conclude that in this regime high-fidelity transfer of quantum states between oscillators is not possible. However, for this ratio becomes a monotonically increasing function of and can in principle be made arbitrary large 333The limiting factor is given by the available power in the pumps and the weak coupling regime.. This remains true irrespective of the cavity’s quality factor and is analogous to atomic physics, where the decoherence due to the spontaneous emission from an excited state can be circumvented by coupling two lower lying states via a far detuned two-photon transition.

iii.2 Two-mode Parametric Amplification

Let us now turn to the case where the resonant coupling results in two-mode parametric amplification as described by Hamiltonian (2). This interaction leads to fast growth of the amplitudes of both modes. The phase of each oscillator separately becomes uniformly distributed, while the phase of their collective motion becomes increasingly well defined. In the absence of dissipation, the resulting two-mode state is not separable and exhibits squeezing of the uncertainty of one collective quadratures of the system below its zero-point value textbook (). Additionally this systems represents a pair of positive and negative mass oscillators containing a quantum-mechanics free subsystem negmassosc (). Even for degenerate modes, the interaction (2) needs auxiliary fields to become resonant. For the condition (28) we easily find


A contributing resonant scattering process is given in Fig. 1 (c) and the side band configuration in frequency space is sketched in panel (b) of Fig. 5.

Figure 5: (Color online) Sketch of the photon power spectral density neglecting the cavity profile. (a) and (b) show the power spectral density due to the two pumps in red (solid) and blue (dashed) respectively. The two pumps, indicated as arrows at and , are arranged in the parametric amplifier configuration specified by Eq. (28). We emphasized the overlapping sidebands captured by Eqs. (44) with dashed lines.

The pumps are arranged such that a blue side band of oscillator 1 is overlapped with the red side band of oscillator 2 and vice versa. Similar “crossing of side bands” was found in other set-ups to lead to entanglement huatang (), as is typically associated with parametric amplification.
The resonant coupling terms are found in the off-diagonal elements of the matrices and . Using Eqs. (44), we find


The resonant terms in the matrices contribute to dissipation


Writing out the master equation using these expressions and separating unitary from dissipative evolution, we arrive at


where the effective Hamiltonian is




and the collective mode is now


Again, the commutation relation of this mode depends on the relative coupling strengths and vanishes for , indicating a quantum-mechanics free subsystem negmassosc (). Since involves annihilation and creation operators, one cannot unambiguously assign an effective temperature to its reservoir. Instead, we will discuss the dynamics resulting from these dissipative terms from a numerical simulation later on. The parameters characterizing the dissipative evolution are


In contrast to the beam splitter case, the collective dynamics contained in Eqs. (51) and (53c) are independent of and depend on the frequency difference only.

iii.2.1 Unitary Evolution

The parameter is proportional to the gain of the parametric amplifier and formally equivalent to with the substitution . It reaches its maximum at and its qualitative behavior is plotted in Fig. 6.

Figure 6: (Color online) Normalized parametric driving strength as a function of central detuning and side band resolution. The line gives parameters with . Other parameters are . Note the formal equivalence up to a different normalization compared to Fig. 3.

The single-quanta processes leading to parametric amplification can be found from Figs. 1 (c) and 5. A photon is transferred from one pump to the other by the creation (or annihilation) of a phonon in each mechanical mode. This also clarifies the creation of entanglement by this interaction, as phonons in the two oscillators are created in pairs and the noise affecting each oscillator becomes correlated. Just as in the beam splitter configuration there are necessarily two distinct pathways for the coherent process separated in frequency by . Consequently systems with photon coherence times longer than the inverse of the frequency difference, , exhibit non-zero detunings where the pathways interfere destructively and lead to . These detunings are drawn as lines in Fig. 6.

For non-vanishing the unitary interaction correlates the motion of both oscillators while increasing their energy exponentially. Whether the mediated interaction can entangle the two mechanical modes depends on the ratio of to the total decoherence rate in the system.

iii.2.2 Dissipative Dynamics

For the parametric amplifier the dissipative behavior is somewhat less intuitive than in the beam splitter situation. While the effect of the two pseudo-thermal, single-mode reservoirs leads to decoherence of each oscillator’s quantum state with rate towards a thermal state with occupations , the collective mode is coupled to a non-trivial reservoir. To understand its phenomenology we solved the master equation


numerically for the initial vacuum state. The probability distributions for the two single mode oscillators and one of their collective modes after are depicted in Figs. 7 together with half-maximum contours compared to their initial state. The dissipation term (54) increases the variance of both quadratures for each oscillator separately, albeit at a different rate. On the collective mode it heats only a single quadrature, leaving the orthogonal one unchanged. These dynamics are the result of a back-action evading measurement on the collective mode . For , this evolution is accompanied by the diffusion of the orthogonal quadrature due to the non-vanishing value of in addition to the unitary evolution of parametric amplification.

Figure 7: (Color online) Phase-space probability distributions (Wigner functions) for: (a): oscillator 1, (b): oscillator 2, (c): collective mode under evolution of given by Eq.(54). The dashed black line is the contour of half maximum, while the thick white one marks the same contour for a vacuum state.

Just as in the beam splitter case, a cavity with small dissipation can resolve the different dissipation channels. If the cavity decay rate is smaller than the sum of the mechanical frequencies, , one can choose whether the system couples to the single-mode reservoirs through (53a) and (53b) or to the collective, non-trivial reservoir (53c). If the decay rate is furthermore smaller than the frequency difference, , either one of the two different collective dynamics governed by Eq. (53c) can be chosen separately.

iii.2.3 Entanglement

The orientation of the squeezed collective quadrature in phase-space is not aligned with the squeezing from two-mode parametric amplification of the effective Hamiltonian (2), which occurs at an angle of with respect to the axes chosen in Fig. 7 (c)textbook (). In general this dissipative behavior will therefore counteract the entanglement generated by the unitary evolution. The relevant quantity which determines if the two oscillators can be driven into an inseparable state is the ratio


Two oscillators can be driven into an entangled state for . We plotted this quantity in dependence of and in Figs. 8 for different parameter regimes.

Figure 8: (Color online) Effective coherences (as defined in Eq. (55)) as a function of central detuning and cavity decay rate .(a) and (c) are for widely separated oscillators, , while (b) and (d) are for . The red line denotes the value . Other parameters are .

The behavior of is qualitatively similar to the quantum/classical transition of linearly coupled oscillators found in mediatedinteractions (). Perhaps most importantly, scales linearly with the central detuning for high optical decay rates and can in principle be made arbitrarily large. For very large detunings this requires high pump powers, since the effective coupling decreases due to smaller intracavity photon numbers, , and weaker effective coupling , see Eq. (51).

Just as in the beam splitter case the ratio of unitary coupling strength to decoherence rate is analogous to two-photon transitions in atomic physics textbook (). For large cavity-detunings the lossy photonic state becomes only virtually excited and thus decoherence from its dissipation can be made arbitrarily small. This allows to create non-classical correlations in a cavity with a low quality factor for sufficiently large central detunings and sufficiently high pump-powers.

The dissipation rate in the denominator of Eq. (55) assumes that information contained in the photons emitted by the cavity is discarded. Using the outcomes of measurements performed on these photons to post-select data – or alternatively using these photons directly to control a feed-back mechanism – would allow to create entangled mechanical modes even in regimes where . A detailed study of such mechanisms is beyond the scope of this work.

iii.3 Thermal Heat bath

Throughout this article we have neglected the independent thermal heat baths for each oscillator. The derivation of the master-equation relies on the coupling strength of each oscillator to its bath being much smaller than the cavity decay rate ; a condition typically satisfied for optomechanical systems review (). The discussed coherent effects are negligible if thermal decoherence rates dominate the system,

where denotes the occupation of oscillator ths thermal reservoir. Even for such thermal systems, however, the resonant coupling can significantly alter the dynamics if the effective coupling is at least of the order of the reservoir couplings,

These conditions are not necessarily restrictive because each can be tuned via the intracavity photon numbers .

Iv Summary

We investigated a scheme to couple non-degenerate mechanical modes resonantly using optomechanical coupling to a single cavity mode. Applying a bi-chromatic pump leads to time-dependent optical springs that can tune couplings between non-degenerate modes into resonance. We derived the reduced master equation for the two oscillators containing frequency shifts, unitary coupling and non-unitary effects such as dissipation and back-action. It is generally applicable for optomechanical systems in the weak coupling regime, ranging from macroscopic oscillators to ultracold atomic systems. We applied the master equation to discuss two resonant cases in depth: a beam splitter and a parametric amplifier. While the former allows more general cooling schemes and quantum state transfer between different frequency regimes, the latter could realize entanglement of spatially separated mechanical modes. We showed that these two fundamental quantum operations can be executed coherently irrespective of the cavity quality factor. Additionally, both cases exhibit mechanical quantum interference for sufficiently low cavity decay rates.

The coupling scheme allows customized interactions among distinct elements in arrays consisting of more than two mechanical modes array () and the derived master equation is easily generalized to include more modes. Beyond mechanical excitations, the presented theory can be adapted to describe a variety of hybrid quantum systems, for instance ensembles of spins interacting with a single cavity mode spinoptodyn () or spatially separated nitrogen vacancy centers in diamond coupled by a high-frequency phonon mode nvdiamonds ().

This work was supported by the Air Force Office of Scientific Research, NSF and a grant from the Swiss National Science Foundation.


  • (1) P. Meystre, Annalen der Physik 525, 215 (2013),
    T. J. Kippenberg and K. J. Vahala, Science 29, 1172 (2008).
    M. Aspelmeyer, T. J. Kippenberg, F. Marquardt, arXiv:1303.0733 (2013),
    D. M. Stamper-Kurn, arXiv:1204.4351 (2012).
  • (2) D. Brooks et al. Nature 488, 476 (2012).
  • (3) A. D. O‘Connell et al. Nature 464, 697 (2010).
  • (4) N. Brahms, T. Botter, S. Schreppler, D. W. C. Brooks and D. M. Stamper-Kurn, Phys. Rev. Lett. 108 , 133601 (2012).
  • (5) A. H. Safavi-Naeini et. al. Phys. Rev. Lett. 108, 033602 (2012).
  • (6) S. Schreppler et. al. Science 344, 1486 (2014).
  • (7) T. A. Palomaki, J. D. Teufel, R. W. Simmonds and K. W. Lehnert, Science 8, 710 (2013).
  • (8) O. Romero-Isart et. al. Phys. Rev. Lett. 107, 020405 (2011).
  • (9) R. Kaltenbaek et. al. Exp. Astron. 34, 123 (2012).
  • (10) L. F. Buchmann, L. Zhang, A. Chiruvelli and P. Meystre, Phys. Rev. Lett. 108, 210403 (2012).
  • (11) H. Tan, F. Bariani, G. Li, and P. Meystre, Phys. Rev. A 88, 023817 (2013).
  • (12) Y. Chen, J. Phys. B: At. Mol. Opt. Phys. 46, 104001 (2013).
  • (13) P.-L. Yu, T. P. Purdy, C. A. Regal, Phys. Rev. Lett. 108, 083603 (2012).
  • (14) S. Chakram, Y. S. Patil, L. Chang, M. Vengalattore Phys. Rev. Lett. 112, 127201 (2014).
  • (15) Ferdinand Brennecke, Stephan Ritter, Tobias Donner, Tilman Esslinger, Science 322, 235 (2008).
  • (16) T. Botter, D. W. C. Brooks, S. Schreppler, N. Brahms, and D. M. Stamper-Kurn, Phys. Rev. Lett. 110, 153001 (2013).
  • (17) M. J. Woolley, A. A. Clerk, Phys. Rev. A 89, 063805 (2014).
  • (18) S. Chakram, Y. S. Patil, M. Vengalattore, arXiv:1412.8536 (2014).
  • (19) H. Seok, L. F. Buchmann, E. M. Wright, and P. Meystre, Phys. Rev. A 88, 063850 (2013)
  • (20) G. Heinrich, F. Marquardt Europhys. Lett. 93, 18003 (2011).
  • (21) F. Massel et al. Nature Communications 3, 987 (2012).
  • (22) M. J. Woolley and A. A. Clerk Phys. Rev. A 87, 063846 (2013).
  • (23) L. F. Buchmann and D. M. Stamper-Kurn, Annalen der Physik 527, 156 (2015).
  • (24) L. F. Buchmann, E. M. Wright, and P. Meystre, Phys. Rev. A 88, 041801(R) (2013).
  • (25) T. Ojanen and K. Borkje, Phys. Rev. A. 90, 013824 (2014).
  • (26) A. B. Shkarin et al. Phys. Rev. Lett. 112, 013602 (2014)
  • (27) Q. Lin et al. Nature Photonics 4, 236 (2010).
  • (28) Bassan, Massimo (Ed.), Advanced Interferometers and the Search for Gravitational Waves, (Springer, 2014).
  • (29) Hu, Yi-Wen and Xiao, Yun-Feng and Liu, Yong-Chun and Gong, Qihuang, Frontiers of Physics 8, 475 (2013).
  • (30) P. Meystre, M. Sargent III, Elements of Quantum Optics, 4th Ed. (Springer, 2010).
  • (31) M. Schmidt, Max Ludwig and Florian Marquardt, New J. Phys. 14, 125005 (2012).
  • (32) L. Diosi, Phys. Rev. A 40, 1165 (1989).
  • (33) Dvir Kafri and J. M. Taylor, arXiv:1311.4558 (2013).
  • (34) Kafri, D., J. M. Taylor, and G. J. Milburn, arXiv:1401.0946 (2014).
  • (35) H. Tan, L. F. Buchmann, H. Seok and G. Li, Phys. Rev. A 87, 022318 (2013).
  • (36) M. Tsang and C. M. Caves, Phys. Rev. X 2, 031016 (2012).
  • (37) D E Chang, A H Safavi-Naeini, M Hafezi and O Painter, New J. Phys. 13 023003 (2011).
  • (38) N. Brahms and D. M. Stamper-Kurn, Phys. Rev. A 82, 041804(R) (2010).
  • (39) P. Ovartchaiyapong, K. W. Lee, B. A. Myers and A. C. Bleszynski Jayich, Nature Communications 5, 4429 (2014).
Figure 9: (Color online) Corrections of the resonant pump detunings from their bare value. Top row for the beam splitter setting , bottom row for parametric amplification, . All plots have and the optical dissipation rate is for (a) and (d), for (b) and (e) and for (c) and (f).

Appendix A Nonlinear corrections to the resonance conditions

Adopting the notation of Sec. II, with and denoting the bare mechanical frequencies and their optical spring shift respectively, the two modified resonance conditions read


These resonance conditions are nonlinear equations with respect to the two pump frequencies because the individual springs depend on the position of each pump with respect to the cavity resonance. While the individual shifts are subtracted from each other in the case of a beam splitter coupling, they add up for the parametric amplifier. For mechanical modes with identical couplings they remain small corrections to the bare frequency difference/sum, but since they provide limits to the fidelity of coherent state transfer or the gain of parametric amplification and the resulting entanglement, these effects have to be taken into account for applications in the quantum regime. The extent of the correction is measured by the dimensionless quantities


for the beam splitter and parametric amplifier respectively. We have calculated these expressions numerically and plotted them as a function of the bare mechanical frequency difference and the average of the pump detunings in Fig. 9.

Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description