## Abstract

When a trichromatic laser field is applied to a cavity optomechanical system within the single-photon strong-coupling regime, we find that the motion of mirror can evolve into a dark state such that the cavity field mode cannot absorb energy from the external field. Via tuning three components of the pumping field to be resonant to the carrier, red-sideband and blue-sideband transitions in the displaced representation respectively, the state of mirror motion can exhibit non-classical properties, such as that in the Lamb-Dicke limit, the state evolves into a squeezed coherent state, and beyond the limit, the state can become a squeezed non-Gaussian state.

© 2014 Optical Society of America

## 1. Introduction

A great surge of interest in macroscopic mechanical systems covers manipulation of the center-of-mass motion [1–4], preparation of nonclassical states of mechanical oscillators [5] and light fields [6, 7], etc. In special, preparations of the nonclassical states of mechanical oscillators can demonstrate the quantum control at large scales [8], which possess potential applications in such as precision measurement [9], long-distance quantum communication and networking [10–12]. Nonclassical states such as squeezed state [13], Schrödinger ‘cats’ [14] and Fock states [15] have been already investigated in various optomechanical systems. In these studies linearized Langevin equation or master equation to quantum fluctuations around classical steady states approach is employed to find system observables, in which the position of the mechanical oscillator linearly modulates the cavity frequency while the cavity field is driven by a strong optical field [16, 17].

With the rapid technical advancement, the single-photon strong-coupling regime has been experimentally realized in ultracold atomic gas system [6, 18, 19], where the presence of a single photon displaces the mechanical oscillator by more than its zero-point uncertainty. There are several micro-optomechanical systems [20, 21] currently approaching to the regime in optical domain, and the microwave optomechanical circuit system also gives the possibility to approach the elusive single-photon strong coupling limit [22, 23]. In this regime the single-photon coupling strength becomes comparable with the frequency of mechanical oscillator and larger than the damping rate of the cavity mode, and then the inherently nonlinear interaction becomes much more significant. The nonlinearity is much helpful to create more general, possibly more interesting and useful states at the level of single photons and phonons, especially non-Gaussian quantum states with negative Wigner funcitons, which indicate strong nonclassicality and genuine quantum behavior and have advantages over Gaussian ones in quantum information protocols [24]. For example, Ren *et al.* [25] proposed a heralded probabilistic scheme to generate mechanical NOON state via manipulating photon transport in the single-photon strong-coupling optomechanical systems, Rabl [26] investigated the photon blockade effect in the single-photon optomechanical system, which may provide an essential ingredient for potential applications as a quantum nonlinear device, and Akram *et al.* [27] proposed a scheme to prepare entangled mechanical ‘cat’ states via conditional single-photon optomechanics.

In this paper we study the preparation of non-classical states of mirror motion when a trichromatic field is applied to a cavity optomechanical system within the single-photon strong-coupling regime. Via utilizing the quantum interference effects [28], which are much useful in the optomechanical system such as manipulation of electromagnetically induced transparency (EIT) to produce desired transmission [29–32], the motion of the mirror may evolve into a dark state such that the cavity-field mode cannot absorb energy from the external field, and accordingly, the optomechanical cavity is decoupled from the external field. However, compared with [28] where a bichromatic laser field is applied to couple the carrier and red-sideband transitions in the displaced representation, in our system when a third frequency component couples to blue-sideband transition, we will demonstrate that more interesting dark states can be generated. We resort to variance of quadrature operator and Wigner function to analyze the non-classical properties of the states and we find that in the Lamb-Dicke limit, i.e., the single-photon coupling strength much smaller than the mechanical frequency, the obtained state can be a squeezed coherent state, while beyond the limit, i.e., in the strong-coupling regime, the obtained state can become a squeezed non-Gaussian state with the negativity of Wigner function. As a consequence, for Gaussian and non-Gaussian optomechanical system the use of an extra frequency component to drive the blue-sideband transition can lead to the squeezed state of mechanical motion [33–35]. Moreover, in contrast to earlier generation of non-classical mechanical states relying on conditional measurement [36] and non-Gaussian input light [5], we can achieve the non-Gaussian mechanical state only with use of a trichromatic laser field.

The organization of the paper is as follows. In Sec. 2, the single-photon optomechanical system is introduced. In Sec. 3 the preparation of non-classical states of mirror motion and the validity of the optomechanical state are discussed, and in the last the conclusion is drawn.

## 2. Description of the single-photon optomechanical system

We consider a general optomechanical system, in which the cavity field couples to a movable end mirror via radiation pressure. The cavity field is driven by a trichromatic laser field with respective frequencies *ω _{i}* and amplitudes Ω

*(*

_{i}*i*= 1, 2, 3). The setup is schematically shown in Fig. 1, and the Hamiltonian of the system is given by (

*ħ*= 1)

*â*and

*b̂*are the annihilation operators of cavity and mechanical modes and their resonant frequencies are

*ω*and

_{c}*ω*respectively. Single-photon coupling strength is

_{m}*g*=

*ω*

_{c}x_{0}/

*L*, where ${x}_{0}=\sqrt{\hslash /(2M{\omega}_{m})}$ is the zero-point fluctuation of the mirror with its mass

*M*, and cavity’s rest length is

*L*. In a frame rotating at the cavity’s frequency

*ω*, the transformed Hamiltonian reads

_{c}*=*

_{i}*ω*−

_{i}*ω*are the detunings of the laser components from the cavity frequency.

_{c}The first two terms of *Ĥ′* correspond to a Hamiltonian *Ĥ*_{0} without driving, which conserves the photon number, i.e., [*â*^{†}*â*, *Ĥ*_{0}] = 0, and it can be diagonalized by introducing a displaced oscillator basis as

*Ĥ*

_{0}is

*n*〉

*and |*

_{c}*p*〉

*are the number states of the cavity and mirror, and $D(ng/{\omega}_{m})=\text{exp}\left[\frac{ng}{{\omega}_{m}}({\widehat{b}}^{\u2020}-\widehat{b})\right]$ is the displacement operator. The state |*

_{m}*p̃*(

*n*)〉

*is an*

_{m}*n*-photon displaced number state [37]. The corresponding eigenvalue is

*ε*=

_{n,p}*pω*−

_{m}*n*

^{2}

*g*

^{2}/

*ω*, where the energy of phonon state shifts

_{m}*g*

^{2}/

*ω*induced by a single-photon radiation pressure.

_{m}Within the eigenbasis of *Ĥ*_{0}, the cavity annihilation operator *â* is expressed in the form

*ξ*=

*g/ω*indicates the oscillator displacement in units of 2

_{m}*x*

_{0}that is related to effects of the radiation pressure of a single photon [38], and ${L}_{r}^{s}(x)$ are associated Laguerre polynomials. Then, from Eq. (5) the Hamiltonian

*Ĥ′*can be reexpressed as

*κ*is much smaller than the oscillating frequency

_{c}*ω*of mirror, the energy levels are clearly separated. Due to the strong single-photon radiation pressure interaction, the eigenstates of the system are anharmonic and then through appropriately choosing the detunings Δ

_{m}*the different non-classical states of the movable mirror may be generated.*

_{i}## 3. Preparation of non-classical states for the mirror motion

Now we turn to the problem of how to effectively reduce the Hamiltonian (8) to a desired form, which can be shown to operate in a finite-dimensional Hilbert space. In the case that the frequency components *ω _{i}* of trichromatic field are equally spaced, the detunings Δ

*satisfy the conditions:*

_{i}*Ĥ*

_{eff}describes the resonant interactions and is given in the form

*V̂*describes the off-resonance transitions

*V̂*excludes those terms contained in

*Ĥ*

_{eff}, and the off-resonance detunings are

*δ*(

_{i}*n*,

*p*,

*p′*) are sufficiently large compared with the effective coupling strength ${A}_{p,{p}^{\prime}}^{(n)}{\mathrm{\Omega}}_{i}$, the transitions described by

*V̂*can be neglected in the rotating-wave approximation. Consequently, if the cavity is initially in the vacuum state, the system will be confined within the zero- and single-photon states. This result can be also interpreted as the photon blockade effect [26], which is based on the fact that the excitation to two-photon space from one-photon states is far off-resonance, i.e., ${\delta}_{i}(2,p,{p}^{\prime})\gg {A}_{p,{p}^{\prime}}^{(2)}{\mathrm{\Omega}}_{i}$. Accordingly, the excitation to higher photon number space is much weaker and negligible.

The Hamiltonian *Ĥ*_{eff} describes the resonant transitions between zero- and one-photon states, which is sketched in Fig. 2. In there the frequency element *ω*_{1} resonantly drives the carrier transitions |*ψ*_{0,}* _{p}*〉 ↔ |

*ψ*

_{1,}

*〉 , element*

_{p}*ω*

_{2}resonantly couples to the red-sideband transitions |

*ψ*

_{1,}

*〉 ↔ |*

_{p}*ψ*

_{0,p+1}〉 and element

*ω*

_{3}resonantly couples to the blue-sideband transitions |

*ψ*

_{0,}

*〉 ↔ |*

_{p}*ψ*

_{1,p+1}〉. In the system there exists a class of dark states |

*D*〉 for the effective Hamiltonian in Eq. (11), which is the eigenvector of zero eigenvalue, i.e.,

*Ĥ*

_{eff}|

*D*〉 = 0 [28]. The explicit expression for the dark state is given by

*β*

_{0}= 1 and

*C*is a normalization constant. It should be noted here the cavity mode is in vacuum and in the following explicit discussions we can numerically prove that the cavity mode is empty. Moreover, when we consider the damping rate of the cavity field, which is represented by the Liouvillian operator

*𝒦̂ρ̂*=

*κ*/2(2

*âρ̂â*

^{†}−

*â*

^{†}

*âρ̂*−

*ρ̂â*

^{†}

*â*), the dark state is also the eigenvector of operator

*𝒦̂*with zero eigenvalue. Thus the cavity damping does not influence the dark state and may even help to generate the dark state. By substituting the expression into the equation, we obtain a chained set of recursive relation

_{3}= 0, i.e., blue-sideband transition is not driven, and the recursive relation in Eq. (15) is consistent with the result therein. However, with use of the frequency element

*ω*

_{3}to couple to the blue-sideband transition, we will display that the obtained dark state can exhibit squeezing in both the weak- and strong-coupling regime. Especially in the strong-coupling regime, the mechanical dark state can possess both the squeezing and non-Gaussian quantum properties.

From the Hamiltonian (11), we can see that for a given positive integer *N*, the Hilbert space can be confined within *N* phonon number states if one sets the transition coefficients to the *N* +1 phonon number states
${\mathrm{\Omega}}_{2}{A}_{N+1,N}^{(1)}$ and
${\mathrm{\Omega}}_{3}{A}_{N,N+1}^{(1)}$ zero. Thus, there will be none excitations from *N* phonon states to *N* + 1 phonon states, as shown by the dashed lines in Fig. 2. The explicit expressions for
${A}_{N+1,N}^{(1)}$ and
${A}_{N,N+1}^{(1)}$ read

*g*=

*g*. Here we take

_{N}*g*the smallest positive value, and the phonon state can be limited within

_{N}*N*phonon number space.

With the special value of *g _{N}*, the recursive relations in Eq. (15) become a group of

*N*+ 1 linear equations. With the condition

*β*

_{0}= 1, all the coefficients

*β*can be calculated from the first

_{p}*N*recursive equations, and to obtain the self-consistent results the last equation

_{1}from the values of Ω

_{2}and Ω

_{3}. In order to facilitate the mathematical treatment, we write the chained set of finite linear recursive equations in Eq. (15) into a matrix form where

*N*+ 1) × 1 column vector and

*T*denotes the transpose of a vector, and the coefficient matrix

*M*of linear equations is a (

*N*+ 1) × (

*N*+ 1) tridiagonal matrix and written as

Because of *β*_{0} = 1, there are *N* terms to calculate. Therefore, in order to obtain the nontrivial solutions of the equations, the rank of coefficient matrix *M* should be *N*. The condition of rank(*M*) = *N* can be achieved via utilizing the recurrence relation of continuant [39], and choosing an appropriate value of Ω_{1} for a given group of Ω_{2} and Ω_{3} together with the coupling strength *g _{N}* given from Eqs. (16).

#### 3.1. Squeezed coherent mechanical state within the Lamb-Dikce Limit

In the following we first consider the dark state of the oscillator motion when the system is well within the Lamb-Dicke limit *g* ≪ *ω _{m}*, i.e.,

*ξ*≪ 1. From the expressions of the associated Laguerre polynomials, ${A}_{p-1,p}^{(1)}$, ${A}_{p,p}^{(1)}$ and ${A}_{p+1,p}^{(1)}$ in the Lamb-Dicke limit become

*H*(

_{p}*z*) are the Hermite polynomials [40]. To fulfill the relation

*g*≪

*ω*, we can choose the

_{m}*g*= 0.1916

_{N}*ω*with the mirror’s motion state confined within phonon number

_{m}*N*= 99, where the value

*ξ*

^{2}= 0.037 well falls in the approximation in Eq. (21). In Fig. 3, we plot the numerical result of the phonon number distribution

*P*calculated from Eq. (18) compared with that of squeezed coherent state obtained from Eq. (24) with the parameters: Ω

_{m}_{1}= 0.0038

*ω*, Ω

_{m}_{2}= 0.0225

*ω*, Ω

_{m}_{3}= −0.015

*ω*, where the populations are hardly distributed above the phonon number

_{m}*n*= 20, and we plot the distributions within phonon number

*n*= 20 to clearly show the similarity. The quadrature operator of phonon mode is $\widehat{X}=({\widehat{b}}^{\u2020}+\widehat{b})/\sqrt{2}$ and the variance of

*X̂*is determined by the relation The variance calculated from the analytical coefficient in Eq. (24) is (Δ

*X̂*)

^{2}= 0.5(

*ν*−

*μ*)

^{2}= 0.1(7dB), and from the numerical results we can obtain the variance (Δ

*X̂*)

^{2}= 0.0952(7.2dB). The variance is smaller than the fluctuation of coherent state 0.5, and thus the mechanical motion is squeezed. Moreover, to quantitatively describe how close of the analytical state and the numerical result, we resort to fidelity

*F*= |〈

*ϕ*

_{num}|

*ϕ*

_{sqcoh}〉|

^{2}, where |

*ϕ*

_{num}〉 denotes the numerical result and |

*ϕ*

_{sqcoh}〉 denotes squeezed coherent state. Numerical calculation of fidelity between the two states is

*F*= 99.9% with the parameters in Fig. 3, which can display the consistency between two results.

Thus in the Lamb-Dicke limit, the dark state of the mirror motion is well described by a squeezed coherent state. As compared with the case of two-tone laser applied on the cavity [28], with use of a third frequency component to drive the blue-sideband transition, the dark state becomes squeezed. In physics, it is analogous to the technique of reservoir engineering in the linearized optomechanical system in [41], where an additional field is applied to couple the blue-sideband transition based on the red-sideband cooling model, and the motion of the oscillator evolves into a dark state which is a squeezed vacuum state. Moreover, if we remove the frequency element *ω*_{1}, i.e. Ω_{1} = 0, the model in Fig. 2 is similar to the model in [41], and the recursive relation (23) becomes the recursive equation of squeezed vacuum state [40]. The motion of the mirror is in a squeezed vacuum state. And when a trichromatic field is applied, the motion of the mirror is in a squeezed coherent state.

In order to gain insight of quantum distribution of the dark state in phase space, we use the Wigner function to present the differences between the classical and quantum behaviors in phase space. The Wigner function is defined as

*ρ̂*is the density operator and

*D̂*(

*α*) = exp(

*αb̂*

^{†}−

*α*

^{*}

*b̂*) is the displacement operator [42]. Wigner function of the dark state is plotted in phase space in Fig. 4, where

*X*and

*Y*are position and momentum. The distribution is suppressed along the

*X̂*quadrature component and presents crescent shape, which can provide a clear insight: the state is close to a squeezed coherent state [43]. Thus in the Lamb-Dicke limit, the dark state of mechanical motion is Gaussian. According to Hudson’s theorem [44] for pure states, the border between Gaussian and non-Gaussian states coincides exactly with the one between states with positive and negative Wigner functions. Obviously there is no negative region in Fig. 4. However, we know that negative Wigner functions of non-Gaussian quantum states indicate strong nonclassicality and genuine quantum behavior, and non-Gaussian states have advantages over Gaussian ones. The use of mechanical nonlinearity will help generate the non-Gaussian state, and thus we will investigate the mechanical behavior in the single-photon strong-coupling regime.

#### 3.2. Non-Gaussian mechanical state beyond Lamb-Dicke limit

Further, when beyond the Lamb-Dicke limit, the nonlinear optomechanical effects become significant and may play a vital role in preparation of non-Gaussian steady state of the mechanical oscillator [38], such as a strongly negative Wigner density and photon-photon correlation function if the optomechanical coupling is comparable to or larger than the optical decay rate and the mechanical frequency [45–47]. It is also possible to extend the squeezing beyond squeezed coherent state and consider the squeezed non-Gaussian state, which can provide a basis for several proposals of quantum information protocols [48,49]. The squeezed non-Gaussian state has been investigated using the cross-Kerr nonlinearity and the corresponding Wigner function has a specific crescent shape as well as negative values [50].

In the paper we also resort to Wigner function to show the effects of nonlinear coupling on preparation of nonclassical state of the oscillator motion, since for pure states, according to Hudson’s theorem [44] the presence of negative regions in the Wigner distribution is an authentic sign of quantum character. Beyond the Lamb-Dicke limit, such as that when the motion state is confined within the number state *N* = 13, the single-photon coupling strength is *g _{N}* = 0.5124

*ω*, and Wigner function of the dark state in Fig. 5 presents the crescent shape of the squeezed state as well as negative values. Here we choose Ω

_{m}_{2}= 0.0225

*ω*, Ω

_{m}_{3}= −0.015

*ω*that govern the degree of squeezing in the Lamb-Dicke regime, and then numerically determine Ω

_{m}_{1}= 0.0029

*ω*from Eq. (20) with the requirement of rank(

_{m}*M*) =

*N*. The variance of the quadrature operator

*X̂*is numerically obtained, which is equal to 0.386(1.12dB) and below the fluctuation of the coherent state. Thus the mechanical state is squeezed even in the single-photon strong-coupling limit, and the multi-tone driving approaches to generate squeezing in [51] can be also useful beyond the linearized optomechanical system. Meanwhile, the minimum of the Wigner function is −0.1. Therefore, in the strong-coupling optomechanical system driven by a trichromatic field, the mirror motion can be in a squeezed non-Gaussian state.

In physics, because of nonlinear coupling in the single-photon strong-coupling limit, the effective coupling strength between different phonon states varies, which means that the coupling strength is dependent on phonon number state. For example, we can write the coupling coefficient in the form

*f*(

*p*) is nonlinear function of phonon number

*p*. Then the recursive relation in Eq. (15) becomes nonlinear. The nonlinear coupling between the mechanical states can modify the coherence between different mechanical Fock states and has pronounced quantum interference effects, which is of benefit to the generation of non-Gaussian state [50, 52].

In order to gain a comprehensive insight into the changes of dark states of mechanical motion with the coupling strengths, we choose another group of parameters: *g _{N}* = 0.3678

*ω*(

_{m}*N*= 26), Ω

_{2}= 0.0225

*ω*, Ω

_{m}_{3}= −0.015

*ω*, and Ω

_{m}_{1}= 0.00963

*ω*fulfilling the requirement of rank(

_{m}*M*) =

*N*, where the coupling strength is between the values of well within and beyond the Lamb-Dicke regimes. The Wigner function and phonon number distribution are shown in Fig. 6, where the minimum of Wigner function is −0.05 and the variance of the quadrature operator

*X̂*is 0.142(5.5dB). The dark state is also a squeezed non-Gaussian state. However the negativity of Wigner function is smaller than that of strong single-photon regime while the degree of squeezing is smaller than the squeezed coherent state in Lamb-Dicke regime. With the increasing single-photon coupling strength, the degree of squeezing decreases and negativity of Wigner function increases. It seems that although the squeezing and negativity of Wigner function are the signatures for the non-classical mechanical motional state, they do not possess the correspondence. This is because Wigner function is connected to all orders of correlations, while squeezing is just connected to the second-order correlation, and they do not have a one-to-one match [53].

Moreover, if we change the field amplitudes Ω_{2} and Ω_{3}, in the Lamb-Dicke limit, the dark state of mechanical motion is still a squeezed coherent state except that the degree of squeezing varies followed the ratio of Ω_{2} and Ω_{3}. We now discuss the changes of dark states in the single-photon strong-coupling regime. We choose a group of parameters: *g _{N}* = 0.5124

*ω*(

_{m}*N*= 13), Ω

_{2}= 0.016

*ω*, Ω

_{m}_{3}= −0.01

*ω*, and Ω

_{m}_{1}= 0.002

*ω*fulfilling the requirement of rank(

_{m}*M*) =

*N*. The Wigner function and phonon number distribution are shown in Fig. 7, where the minimum of Wigner function is −0.069 and the variance of the quadrature operator

*X̂*is 0.239(3.21dB). The dark state is a squeezed non-Gaussian state with the smaller negativity of Wigner function and larger degree of squeezing compared with those in Fig. 5.

In addition, one should test the validity of the dark state as well as the effective Hamiltonian (11). We will solve the evolution of the system state |Ψ(*t*)〉 which obeys the Schrödinger equation
$i\frac{d}{dt}|\mathrm{\Psi}(t)\u3009={\widehat{H}}^{\prime}(t)|\mathrm{\Psi}(t)\u3009$ with the Hamiltonian (2) without the rotating-wave approximation. Here we consider the time interval *t* ≫ 1/*κ*, where 1/*κ* represents the readout timescale of the cavity field. For example, a cavity-Cooper pair transistor (cCPT) mechanical resonator, *ω _{m}/κ* ∼ 10 is achievable. The time-dependent dark state governed by effective Hamiltonian (11) is given by

*ε*

_{0,p}=

*pω*. In particular, we choose a dark state as the initial system state |Ψ(0)〉 = |

_{m}*D*(0)〉 and calculate the fidelity

*F*= |〈

*D*(

*t*)|Ψ(

*t*)〉|

^{2}between |Ψ(

*t*)〉 and |

*D*(

*t*)〉 to verify the rotating-wave approximation, where system state $|\mathrm{\Psi}(t)\u3009=|\mathrm{\Psi}(0)\u3009-i{\int}_{0}^{t}d{t}^{\prime}{\widehat{H}}^{\prime}({t}^{\prime})|\mathrm{\Psi}({t}^{\prime})\u3009$. In Fig. 8 we plot the fidelity as a function of time. Fidelity is

*F*≈ 1 for a long time, i.e.

*t*≫

*κ*

^{−1}, and thus the system state |Ψ(

*t*)〉 can be well described by the dark state |

*D*(

*t*)〉 and the effective Hamiltonian (11) is valid. The high frequency oscillation of the fidelity is caused by the off-resonant transitions in Eq. (12), and the phenomena of collapse and revival can be physically understood from the summation of different Rabi oscillations which are related to different values of detunings and the coupling strengths that vary for different phonon number states. However, due to the too many oscillations in the system the revival is hard to reach for the unity.

#### 3.3. Effects of cavity and mechanical dissipations on the mechanical dark state

In the following we will discuss the generation of the dark state for the mirror motion under the help of the cavity dissipation together with the influence of mechanical damping. In the frame of Hamiltonian given in Eq. (2), the evolution of the system with the cavity and phonon damping is governed by the master equation

*ρ̂*is the density matrix of the photon-mirror system,

*γ*is the mechanical damping rate and

*n̄*is the phonon number in the thermal environments. It should be noted here that the dark state in Eq. (14) is also the eigenvector of Liouvillian operator of cavity damping with zero eigenvalue, which can help to generate the dark state. The mechanical quality factor can reach

*Q*∼ 10

^{5}, and the damping rate

*γ*∼ 10

^{−5}

*ω*. Around ground-state cooling of mechanical oscillator and in the timescale

_{m}*κ*

^{−1}≪

*t*≪

*γ*

^{−1}, we can neglect the affect of the mechanical damping on the generation of dark state.

We numerically solve the master equation (29) with the initial ground state of system, and to measure to probability of the system in the dark state we also resort to the fidelity

In the Fig. 9 we plot the fidelity that increases with time and at time*T*= 2500

*ω*

^{−1}the fidelity can reach the value 0.91. It should be noted here the line in Fig. 9 also has small oscillations as there in Fig. 8, but much smaller compared to the degree scale in Y-axis. Thus the dark state can be generated from the ground state of the cavity-mirror system with the help of the cavity damping, and the Wigner function also presents negative values and the minimal negativity is about −0.055.

Finally, there are a number of techniques capable of detecting tiny mechanical vibrations [54], and also various proposals to perform quantum state reconstruction (QSR) of mechanical motion. For example, mechanical quadrature tomography using back-action-evading interactions can provide the position or momentum distributions of mechanical state and Wigner function for the typical optomechanical system [55]. Thus the squeezed non-Gaussian state can be visualized with this approaching techniques.

## 4. Conclusion

In summary, we have investigated the dark state of mirror motion in a single-photon optomechanical system when a trichromatic field is applied to drive a cavity mode. In the displaced representation, three frequency components couple to the carrier, red-sideband and blue-sideband transitions respectively, the mirror motion can evolve into a dark state which arises from the quantum interference effects. Moreover, via carefully tuning single-photon coupling strength, the phonon state can be confined within a low number. In the Lamb-Dicke limit, the mirror motion can evolve into a squeezed coherent state, and beyond the Lamb-Dicke limit, the state evolves into a squeezed non-Gaussian state, where Wigner function of mirror motion state presents the crescent shape of the squeezed state as well as negative values. Finally, we apply the fidelity to exhibit the validity of the non-classical dark state of mirror motion in the system and discuss the effects of cavity damping on the generation of dark state.

## Acknowledgments

This work is supported by the National Natural Science Foundation of China (Grant No. 61275123), the National Basic Research Program of China (Grant No. 2012CB921602), and the Key Laboratory of Advanced Micro-Structure of Tongji University.

## References and links

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

**2. **Z. Q. Yin, “Phase noise and laser-cooling limits of optomechanical oscillators,” Phys. Rev. A **80**, 033821 (2009). [CrossRef]

**3. **J. Zhang, Y. X. Liu, and F. Nori, “Cooling and squeezing the fluctuations of a nanomechanical beam by indirect quantum feedback control,” Phys. Rev. A **79**, 052102 (2009). [CrossRef]

**4. **Y. Li, L. A. Wu, and Z. D. Wang, “Fast ground-state cooling of mechanical resonators with time-dependent optical cavities,” Phys. Rev. A **83**, 043804 (2011). [CrossRef]

**5. **F. Khalili, S. Danilishin, H. X. Miao, H. Müller-Ebhardt, H. Yang, and Y. B. Chen, “Preparing a mechanical oscillator in non-gaussian quantum states,” Phys. Rev. Lett. **105**, 070403 (2010). [CrossRef] [PubMed]

**6. **D. W. C. Brooks, T. Botter, S. Schreppler, T. P. Purdy, N. Brahms, and D. M. Stamper-Kurn, “Non-classical light generated by quantum-noise-driven cavity optomechanics,” Nature **488**, 476–480 (2012). [CrossRef] [PubMed]

**7. **Y. D. Wang and A. A. Clerk, “Reservoir-engineered entanglement in optomechanical systems,” Phys. Rev. Lett. **110**, 253601 (2013). [CrossRef] [PubMed]

**8. **T. J. Kippenberg and K. J. Vahala, “Cavity opto-mechanics,” Opt. Express **15**, 17172–17205 (2007). [CrossRef] [PubMed]

**9. **C. Caves, “Quantum-mechanical noise in an interferometer,” Phys. Rev. D **23**, 1693–1708 (1981). [CrossRef]

**10. **H. J. Kimble, “The quantum internet,” Nature **453**, 1023–1030 (2008). [CrossRef] [PubMed]

**11. **B. Kraus and J. I. Cirac, “Discrete entanglement distribution with squeezed light,” Phys. Rev. Lett. **92**, 013602 (2004). [CrossRef] [PubMed]

**12. **M. Paternostro, W. Son, and M. S. Kim, “Complete conditions for entanglement transfer,” Phys. Rev. Lett. **92**, 197901 (2004). [CrossRef] [PubMed]

**13. **W. J. Gu, G. X. Li, and Y. P. Yang, “Generation of squeezed states in a movable mirror via dissipative optomechanical coupling,” Phys. Rev. A **88**, 013835 (2013). [CrossRef]

**14. **S. Bose, K. Jacobs, and P. L. Knight, “Preparation of nonclassical states in cavities with a moving mirror,” Phys. Rev. A **56**, 4175–4186 (1997). [CrossRef]

**15. **S. Rips, M. Kiffner, I. Wilson-Rae, and M. J. Hartmann, “Steady-state negative Wigner functions of nonlinear nanomechanical oscillators,” New J. Phys. **14**, 023042 (2012). [CrossRef]

**16. **I. Wilson-Rae, N. Nooshi, W. Zwerger, and T. J. Kippenberg, “Theory of ground state cooling of a mechanical oscillator using dynamical backaction,” Phys. Rev. Lett. **99**, 093901 (2007). [CrossRef] [PubMed]

**17. **F. Marquardt, J. P. Chen, A. A. Clerk, and S. M. Girvin, “Quantum theory of cavity-assisted sideband cooling of mechanical motion,” Phys. Rev. Lett. **99**, 093902 (2007). [CrossRef] [PubMed]

**18. **K. W. Murch, K. L. Moore, S. Gupta, and D. M. Stamper-Kurn, “Observation of quantum-measurement backaction with an ultracold atomic gas,” Nat. Phys. **4**, 561–564 (2008). [CrossRef]

**19. **F. Brennecke, S. Ritter, T. Donner, and T. Esslinger, “Cavity optomechanics with a Bose-Einstein Condensate,” Science **322**, 235–238 (2008). [CrossRef] [PubMed]

**20. **A. H. Safavi-Naeini, J. Chan, J. T. Hill, Thiago P. Mayer Alegre, A. Krause, and O. Painter, “Observation of quantum motion of a nanomechanical resonator,” Phys. Rev. Lett. **108**, 033602 (2012). [CrossRef] [PubMed]

**21. **S. Weis, R. Rivière, S. Deléglise, E. Gavartin, O. Arcizet, A. Schliesser, and T. J. Kippenberg, “Optomechanically induced transparency,” Science **330**, 1520–1523 (2010). [CrossRef] [PubMed]

**22. **T. T. Heikkilä, F. Massel, J. Tuorila, R. Khan, and M. A. Sillanpää, “Enhancing optomechanical coupling via the Josephson effect,” arXiv:1311.3802.

**23. **A. J. Rimberg, M. P. Blencowe, A. D. Armour, and P. D. Nation, “A cavity-Cooper pair transistor scheme for investigating quantum optomechanics in the ultra-strong coupling regime,” New J. Phys. **16**, 055008 (2014). [CrossRef]

**24. **J. Eisert, S. Scheel, and M. B. Plenio, “Distilling gaussian states with Gaussian operations is impossible,” Phys. Rev. Lett. **89**, 137903 (2002). [CrossRef] [PubMed]

**25. **X. X. Ren, H. K. Li, M. Y. Yan, Y. C. Liu, Y. F. Xiao, and Q. H. Gong, “Single-photon transport and mechanical NOON-state generation in microcavity optomechanics,” Phys. Rev. A **87**, 033807 (2013). [CrossRef]

**26. **P. Rabl, “Photon blockade effect in optomechanical systems,” Phys. Rev. Lett. **107**, 063601 (2011). [CrossRef] [PubMed]

**27. **U. Akram, W. P. Bowen, and G. J. Milburn, “Entangled mechanical cat states via conditional single photon optomechanics,” New J. Phys. **15**, 093007 (2013). [CrossRef]

**28. **G. F. Xu and C. K. Law, “Dark states of a moving mirror in the single-photon strong-coupling regime,” Phys. Rev. A **87**, 053849 (2013). [CrossRef]

**29. **G. S. Agarwal and Sumei Huang, “Electromagnetically induced transparency in mechanical effects of light,” Phys. Rev. A **81**, 041803(R) (2010). [CrossRef]

**30. **H. Xiong, L. G. Si, A. S. Zheng, X. X. Yang, and Y. Wu, “Higher-order sidebands in optomechanically induced transparency,” Phys. Rev. A **86**, 013815 (2012). [CrossRef]

**31. **S. Huang and G. S. Agarwal, “Electromagnetically induced transparency with quantized fields in optocavity mechanics,” Phys. Rev. A **83**, 043826 (2011). [CrossRef]

**32. **A. H. Safavi-Naeini, T. P. Mayer Alegre, J. Chan, M. Eichenfield, M. Winger, Q. Lin, J. T. Hill, D. E. Chang, and O. Painter, “Electromagnetically induced transparency and slow light with optomechanics,” Nature **472**, 69–73 (2011). [CrossRef] [PubMed]

**33. **A. Mari and J. Eisert, “Gently modulating optomechanical systems,” Phys. Rev. Lett. **103**, 213603 (2009). [CrossRef]

**34. **A. Farace and V. Giovannetti, “Enhancing quantum effects via periodic modulations in optomechanical systems,” Phys. Rev. A **86**, 013820 (2012). [CrossRef]

**35. **W. J. Gu and G. X. Li, “Squeezing of the mirror motion via periodic modulations in a dissipative optomechanical system,” Opt. Express **21**, 020423 (2013). [CrossRef]

**36. **M. R. Vanner, “Selective linear or quadratic optomechanical coupling via measurement,” Phys. Rev. X **1**, 021011 (2011).

**37. **F. A. M. de Oliveira, M. S. Kim, P. L. Knight, and V. Buzek, “Properties of displaced number states,” Phys. Rev. A **41**, 2645–2652 (1990). [CrossRef] [PubMed]

**38. **A. Nunnenkamp, K. Børkje, and S. M. Girvin, “Single-photon optomechanics,” Phys. Rev. Lett. **107**, 063602 (2011). [CrossRef] [PubMed]

**39. **M. E. A. El-Mikkawy, “On the inverse of a general tridiagonal matrix,” Appl. Math. Comput. **150**, 669–679 (2004). [CrossRef]

**40. **J. S. Peng and G. X. Li, *Introduction to Modern Quantum Optics* (World Scientific, 1998).

**41. **P. Rabl, A. Shnirman, and P. Zoller, “Generation of squeezed states of nanomechanical resonators by reservoir engineering,” Phys. Rev. B **70**, 205304 (2004). [CrossRef]

**42. **S. M. Barnett and P. M. Radmore, *Methods in Theoretical Quantum Optics* (Clarendon, 1997).

**43. **U. Leonhardt, *Measuring the Quantum State of Light* (Cambridge University, 1997).

**44. **R. L. Hudson, “When is the wigner quasi-probability density non-negative?” Rep. Math. Phys. **6**, 249–252 (1974). [CrossRef]

**45. **J. Qian, A. A. Clerk, K. Hammerer, and F. Marquardt, “Quantum signatures of the optomechanical instability,” Phys. Rev. Lett. **109**, 253601 (2012). [CrossRef]

**46. **P. D. Nation, “Nonclassical mechanical states in an optomechanical micromaser analog,” Phys. Rev. A **88**, 053828 (2013). [CrossRef]

**47. **N. Lörch, J. Qian, A. Clerk, F. Marquardt, and K. Hammerer, “Laser theory for optomechanics: limit cycles in the quantum regime,” Phys. Rev. X **4**, 011015 (2014).

**48. **W. J. Munro, K. Nemoto, R. G. Beausoleil, and T. P. Spiller, “High-efficiency quantum-nondemolition single-photon-number-resolving detector,” Phys. Rev. A **71**, 033819 (2005). [CrossRef]

**49. **J. Fiurášek, L. Mišta, and R. Filip, “Entanglement concentration of continuous-variable quantum states,” Phys. Rev. A **67**, 022304 (2003). [CrossRef]

**50. **T. Tyc and N. Korolkova, “Highly non-Gaussian states created via cross-Kerr nonlinearity,” New. J. Phys. **10**, 023041 (2008). [CrossRef]

**51. **M. J. Woolley and A. A. Clerk, “Two-mode squeezed states in cavity optomechanics via engineering of a single reservoir,” Phys. Rev. A **89**, 063805 (2014). [CrossRef]

**52. **R. L. de, M. Filho, and W. Vogel, “Nonliear coherent states,” Phys. Rev. A **54**, 4560–4563 (1996). [CrossRef]

**53. **E. V. Shchukin and W. Vogel, “Nonclassical moments and their measurement,” Phys. Rev. A **72**, 043808 (2005). [CrossRef]

**54. **Y. Greenberg, Y. A. Pashkin, and E. II’ichev, “Nanomechanical resonators,” Phys.-Usp. **55**, 382 (2012). [CrossRef]

**55. **M. R. Vanner, I. Pikovski, and M. S. Kim, “Towards optomechanical quantum state reconstruction of mechanical motion,” arXiv:1406.1013.