Dissipative synthesis of mechanical Fock-like states

TL;DR

This paper proposes a scheme to stabilize and generate mechanical Fock-like quantum states in an optomechanical system by combining dissipative squeezing with mechanical nonlinearity, enabling near-perfect Fock state approximations.

Contribution

It introduces a novel method to produce and stabilize arbitrary mechanical Fock-like states using dissipative processes and nonlinearities in an optomechanical setup.

Findings

Achieves near-unit fidelity in approximating Fock states of any number.

Demonstrates stabilization of nonclassical mechanical states via dissipation.

Provides a tunable scheme using three control lasers for state preparation.

Abstract

The observation of genuine quantum features of nano-mechanical motion is a key goal for both fundamental and applied quantum science. To this end, a promising approach is the stabilization of nonclassical features in the presence of dissipation, by means of the tunable coupling with a photonic environment. Here we present a scheme that combines dissipative squeezing with a mechanical nonlinearity to stabilize arbitrary approximations of (displaced) mechanical Fock state of any number. We consider an optomechanical system driven by three control lasers--at the cavity resonance and at the two mechanical sidebands--that couple the amplitude of the cavity field to the resonator's position and position squared. When the amplitude of the resonant drive is tuned to some specific values, the mechanical steady state is found in a (displaced) superposition of a finite number of Fock states, which for large enough squeezing achieves near-unit fidelity with a (displaced) Fock state of any desired number.

Equations

$$\hat{H}={\omega }_c{\hat{a}}^{†}\hat{a}+{\omega }_m{\hat{b}}^{†}\hat{b}-g_0^{(1)}{\hat{a}}^{†}\hat{a}(\hat{b}+{\hat{b}}^{†})-g_0^{(2)}{\hat{a}}^{†}\hat{a}(\hat{b}+{\hat{b}}^{†}{)}^2+{\hat{H}}_{\mathrm{drive}},$$

$$\hat{H}_{\mathrm{drive}}=\hat{a}^{\dagger}\bigl{(}\varepsilon_{-}e^{-i\omega_{-}t}+\varepsilon_{0}e^{-i\omega_{0}t}+\varepsilon_{+}e^{-i\omega_{+}t}\bigr{)}+\mathrm{H.c.}\,,$$

$${\hat{H}}_{\mathrm{RWA}}=-{\hat{d}}^{†}(G_{-}\hat{b}+G_{+}{\hat{b}}^{†}+G_0\{\hat{b},{\hat{b}}^{†}\})+H.c.,$$

$$(G_{-}\hat{b}+G_{+}{\hat{b}}^{†}+G_0\{\hat{b},{\hat{b}}^{†}\})|\phi ⟩=0.$$

$$G_0=\frac{G}{\sqrt{2(2n+1)}},$$

$${\hat{b}}^{†}\hat{b}\hat{D}\left(\sqrt{n+\frac{1}{2}}\right)|\phi ⟩=n\hat{D}\left(\sqrt{n+\frac{1}{2}}\right)|\phi ⟩,$$

$$|{\phi }_n⟩=\hat{D}\left(-\sqrt{n+\frac{1}{2}}\right)|n⟩$$

$$G_0=\sqrt{\frac{G_{+}{G}_{-}}{2(2n+1)}},$$

$${\phi }_n(q)\propto e^{-q\sqrt{\zeta (1+2n)}}{e}^{-\frac{q^2}{2}}{H}_n\left(q+\frac{(1+\zeta )\sqrt{\zeta (1+2n)}}{2\zeta }\right),$$

$${\phi }_n(q)\propto e^{-\frac{1}{2}(q-{\xi }_n{)}^2}{H}_n\left(q-{\xi }_n+\frac{(1-\zeta )\sqrt{\zeta (1+2n)}}{2\zeta }\right),$$

$${\phi }_n(q)\propto \sum _{k=0}^n\left(\genfrac{}{}{0ex}{}{n}{k}\right)c_n^{-k}{e}^{-\frac{1}{2}(q-{\xi }_n{)}^2}{H}_k\left(q-{\xi }_n\right),$$

$$|\varphi_{n}\rangle=\mathcal{N}_{n}\hat{D}\bigl{(}\xi_{n}/\sqrt{2}\bigr{)}\sum_{k=0}^{n}\binom{n}{k}c_{n}^{\,-k}|k\rangle\,,$$

$$\hat{f}=\mathcal{G}\hat{\beta }+\sqrt{\frac{\cosh r\sinh r}{2(2n+1)}}\{{\hat{b}}^{†},\hat{b}\},$$

Full text

\pubvolume

xx \issuenum1 \articlenumber5

\historyReceived: date; Accepted: date; Published: date

\TitleDissipative synthesis of mechanical Fock-like states \AuthorMatteo Brunelli 1 and Oussama Houhou 2,3 \AuthorNamesMatteo Brunelli and Oussama Houhou

\corresCorrespondence: mb2236@cam.ac.uk

1 Introduction

The motional state of atomic or mechanical degrees of freedom can be manipulated via the interaction with the electromagnetic field confined in a cavity. Such a possibility is best illustrated by cavity cooling, which has been successfully applied to single atoms AtomCooling , ions IonCooling , and micro- and nano-mechanical resonators Chan2011 ; Teufel2011 ; Kippen2012 . Recent breakthroughs in the dissipative preparation of mechanical squeezed states MechSqueezing1 ; MechSqueezing2 ; MechSqueezing3 ; ResEngIons3 , where a cavity-assisted scheme is designed to cool the target system directly into a squeezed state of motion, can be thought of as a powerful development of this paradigm ResEng1 ; Kronwald ; Wang ; Woolley ; JieLi . However, for many applications, ranging from fundamental tests of quantum mechanics to quantum information precessing, the stabilization of highly pure states with non-Gaussian features is needed instead. In cavity optomechanics, the quadratic optomechanical coupling has been exploited for the dissipative preparation of Schrödinger cat states OptoCat1 ; OptoCat2 , but the existence of multiple steady states requires the unpractical initialization of the system in a state of definite parity. Recently we have shown that a tunable optomechanical coupling which has both a linear and quadratic component enables the stabilization of pure non-Gaussian states without requiring any initialization Me1 ; Me2 . For specific values of the amplitude of the laser drives new families of nonclassical states can be stabilized, which correspond to (squeezed and displaced) superpositions of a finite number of Fock states. Here we focus on a specific instance, namely on one such (displaced) finite superposition that approximates—in principle with arbitrary fidelity—any number state in the harmonic ladder (modulo a displacement).

2 Results

We consider an optomechanical system where the frequency of a cavity mode parametrically couples to the displacement and squared displacement of a mechanical resonator. The Hamiltonian is given by (we set $\hbar=1$ throughout)

[TABLE]

where $\hat{a}$ ($\hat{b}$) is the annihilation operator of the cavity (mechanical) mode of frequency $\omega_{c}$ ($\omega_{m}$) and $g_{0}^{(1)}$, $g_{0}^{(2)}$ respectively quantifies the linear and quadratic single-photon coupling. Such linear-and-quadratic coupling can be realized in membrane-in-the-middle setups NonLin ; NonLin1 ; NonLin2 , cold atoms NonLin3 , microdisk resonators NonLin4 and photonic crystal cavities NonLin5 ; NonLin6 . The cavity has a decay rate $\kappa$ and is driven with three tones

[TABLE]

applied on the cavity resonance ($\omega_{0}=\omega_{c}$), and on the lower and upper mechanical sideband $(\omega_{\pm}=\omega_{c}\pm\omega_{m})$. After standard linearization (we dub $\hat{d}$ the fluctuation operator of the cavity field), moving into a frame rotating with the free cavity and mechanical Hamiltonian, and focusing on the good cavity limit ($\kappa\ll\omega_{m}$) we get

[TABLE]

where we set $G_{\pm}=g_{0}^{(1)}\alpha_{\pm}$, $G_{0}=g_{0}^{(2)}\alpha_{0}$, and $\alpha_{\pm,0}$ are the steady values of the cavity amplitude at each frequency component; we will assume these couplings to be real and positive without loss of generality. After a transient time the cavity field is found in the vacuum while the mechanical resonator in a pure state $|\varphi\rangle$ that satisfies the condition

[TABLE]

Note that when the nonlinear term in absent, namely $\nobreak{G_{0}\equiv 0}$, we recover dissipative squeezing with a squeezing degree $r=\tanh^{-1}(G_{+}/G_{-})$ Kronwald .

In order to characterize the steady state $|\varphi\rangle$, let us first assume that the amplitudes at the two mechanical sidebands are equal, i.e., $G_{\pm}=G$. In this case it is enough to notice that for the following values of the resonant coupling

[TABLE]

the condition expressed in Eq. (4) becomes

[TABLE]

where $\hat{D}$ is the displacement operator and $n\in\mathbb{N}$ a non-negative integer (to stress this dependence we set $|\varphi_{n}\rangle\equiv|\varphi\rangle$ from now on). This is in turn equivalent to

[TABLE]

and proves that the steady state is indeed a displaced Fock state. In particular, by tuning the amplitude of the resonant drive in Eq. (5) any state in the Fock state ladder can be stabilized.

The class of steady states obtained in Eq. (7) turns out to be unstable Me2 . However, it can be seen as the limit $G_{+}\rightarrow G_{-}$ of the more general case $G_{+}\neq G_{-}$ with

[TABLE]

which is guaranteed to be stable as long as $G_{+}<G_{-}$. In order to find the new steady state, we can project Eq. (4) onto the position eigenstate $|q\rangle$ and obtain a differential equation for the associated wave function $\varphi_{n}(q)$. The solution of such equation reads

[TABLE]

where we set $\zeta=\tanh r\in[0,1)$. Note that the integer order of the Hermite polynomial is determined by the resonant coupling in Eq. (8). By completing the square in the exponent we get

[TABLE]

where $\xi_{n}=-\sqrt{\zeta(1+2n)}$. Note that for $\zeta\rightarrow 1$ we correctly recover the wave function of a displaced quantum harmonic oscillator. We now exploit the following property of the Hermite polynomials, $H_{n}(x+y)=\sum_{k=0}^{n}\binom{n}{k}H_{k}(x)(2y)^{n-k}\,,$ which leads us to

[TABLE]

with $c_{n}=-\frac{(1-\zeta)}{4\zeta}\xi_{n}$. From the last line we can finally read the explicit expression of the state

[TABLE]

where the normalization factor is given by $\mathcal{N}_{n}=\left[{}_{2}F_{1}\left(-n,-n;1;c_{n}^{-2}\right)\right]^{-1/2}$. The steady state is now given by the action of a $n$-dependent displacement on a superposition of a finite number ($n+1$) of elements. It is easily checked that in the limit $\zeta\rightarrow 1$ the superposition collapses to the single element of Eq. (7). On the other hand, for any non-zero value of the squeezing parameter the state $|\varphi_{n}\rangle$ displays negativity in the Wigner distribution and the larger the amount squeezing the closer the resemblance with a Fock state. This feature is clear from Fig. 1, where we show the Wigner distribution for a given $n$ ($n=5$) and different values of the squeezing parameter $\zeta$. We clearly see that the distribution, which for lower values of $\zeta$ is skewed toward one side, progressively straightens to approach that of a Fock state. We can thus think of $|\varphi_{n}\rangle$ as a state that approximates any given displace Fock state, to an extent that improves with the amount of squeezing available. Mechanical dissipation—not considered here—sets a limit on the precision of such approximation. Yet, one can show that it is still possible to approximate with near-unit fidelity any Fock state Me2 .

Coming back to Eq. (4), we notice that $|\varphi_{n}\rangle$ is the state uniquely annihilated by the nonlinear operator

[TABLE]

where $\hat{\beta}=\cosh r\hat{b}+\sinh r\hat{b}^{{\dagger}}$ is a Bogoliubov mode and $\mathcal{G}=\sqrt{G_{-}^{2}-G_{+}^{2}}$. The nonlinear contribution added to the Bogoliubov transformation makes the nature of $\hat{f}$ non bosonic.

3 Discussion

We presented an exactly solvable model to augment dissipative squeezing by means of a quadratic nonlinearity. The model can be implemented in optomechanical cavity and the states stabilized by our protocol approximate displaced multi-phonon Fock state of any desired number.

Acknowledgements.

M. B. is supported by the European Union’s Horizon 2020 research and innovation programme under grant agreement No 732894 (FET Proactive HOT). O. H. acknowledges support from the SFI-DfE Investigator programme (grant 15/IA/2864), the EU Horizon2020 Collaborative Project TEQ (grant agreement No 766900) and from the EPSRC project EP/P00282X/1. \conflictsofinterestThe authors declare no conflict of interest. The funding sponsor had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results. \reftitleReferences