## Abstract

We propose a scheme for generation of a nonlinear coherent state (NCS) of a mechanical resonator (MR) in an optomechanical micro-cavity, in which a two-level quantum dot (QD) and the microcavity are respectively driven by a strong laser and a weak laser. This microcavity can be engineered within a photonic band-gap (PBG) material. By properly tuning the frequency of the weak driving field, two-photon blockade phenomenon occurs. The QD-cavity subsystem can evolve into a dark state due to the damping of the microcavity and the elimination of the decay rate of the QD at selected frequencies in the PBG material. In this situation, the phonon mode of the MR can be prepared into a NCS, which is a non-classical state and possesses the sub-Poisson statistics. We also demonstrate the Wigner function of the NCS, which negativity implies its non-classicality.

© 2016 Optical Society of America

## 1. Introduction

In the emerging field of optomechanical systems, numerous rapid progresses with the radiation pressure have recently been obtained during measuring and cooling the mechanical motion [1,2]. It is reported that the quantum control of mechanical motion and the control of light or microwave field in cavity optomechanical systems can be achieved to some extent [3, 4]. The cooling of a mechanical resonator (MR) to its ground state of motional degree of freedom is an important procedure [5, 6]. Experimentally, it has been realized that micro- and nanomechanical resonator can be cooled to or near its quantum ground state in various mechanical systems either through direct cryogenic cooling [7], or laser cooling using microwave [8, 9] and optical cavity fields [10]. When the MR is cooled to the ground state, it is possible to generate and detect the non-classical states of the MR, such as the squeezed state [3, 11, 12] and the entangled state [13]. The squeezed state of a movable mirror can be obtained via feeding broadband squeezed vacuum light into the cavity in dissipative optomechanics [11]. By driving an optomechanical cavity with two controllable lasers with different amplitudes, the dissipative mechanism with the driven cavity acts as an engineered reservoir, the squeezed state of a MR can also be obtained [14]. Recently, a scheme for generation of the non-classical state of the MR is proposed, in which the MR is nonlinearly coupled to a two-level system [15]. The non-classical states of the MR can also be produced by applying the generation of motional dark states [16].

The nonlinear coherent state (NCS) as another concerned non-classical states has also attracted many interests in view of its promising application in quantum algebras [17–19], which is based on noncommutative space rather than the usual space [20]. In the description of the motion of a trapped ion, Vogel *et al.* put forward the definition of the NCS as the eigenstate of the non-Hermitian operator *f*(*n̂*)*a* [21,22], where the deformation function *f*(*n̂*) is an operator-valued function of number operator *n̂* = *a*^{†}*a*. In this case, the commutator [*f*(*n*)*a*, *a*^{†} *f*(*n*)] is not a constant or a linear function of the generators of the Lie algebra but nonlinear in the generators. It is found that the NCS exhibits non-classical features such as amplitude squeezing and self-splitting, which are accompanied by pronounced quantum interference effects [21,23]. The generation schemes of the NCS are explored in some other literatures, e.g. a single-atom laser [17], a micromaser under the intensity-dependent Jaynes-Cummings model [24], an exciton in a wide quantum dot [25], and a trapped ion [21,26]. The stationary state of a single-atom laser is shown to be a phase-averaged NCS and the interaction-induced transition probabilities of the single-atom laser effectively depend on the field intensity [17]. The quantized motion of an exciton in a wide quantum dot interacting with two laser beams can be prepared in a NCS and the exciton states can be utilized for quantum information processes [25]. The NCS is employed to discuss the optimization of quantum information which is important in binary (or multibinary) communication [27]. Usually, the discussion of NCS is related with the non-Gaussian quantum states. The occurrence of the negativity of the Wigner function (WF) for some non-Gaussian states implies that they possess the non-classical properties. The non-Gaussian state plays a significant role in entanglement distillation [28, 29], quantum computation with cluster states [30], loophole-free Bell tests with the continuous variables [31, 32], and quantum error correction [33]. Besides, the non-Gaussian state can be exploited to improve the quantum communication protocols in quantum information processing [34] and improve the optimal estimation of losses at the ultimate quantum limit [35]. Recent years, it is found that the optomechanical systems possess the intrinsic nonlinearity thus hold promising applications in quantum information processing [36] and ultra-precision measurement [37]. So how to prepare the NCS of a MR in the optomechanical system is a topic of great significance.

The generation of a non-classical light source is one of the major requirements of quantum information processor associated with modern semiconductor fabrication technology [38]. For this purpose, employing photon blockade phenomenon is one of methods. This phenomenon is first proposed by Imamoḡlu [39], which occurs when the absorption of a first input photon by an optical device blocks the transmission of a second one, thereby leading to non-classical output photon statistics. The inherent non-classical property at single photon level provides a way to control single photon, which is commonly used in quantum key distribution system [40] and plays a fundamental role in proposed devices for quantum information processing [38]. The photon blockade is first observed in cavity quantum electrodynamics(QED) systems with the strong atom-cavity coupling [41], and then has been demonstrated in different solid-state schemes, e.g. in cavity QED systems with a quantum dot (QD) in photonic crystal cavities [42], with a QD strongly coupled to a photonic crystal resonator [43] and other circuit-QED systems [44, 45]. Recently, much attention has been focused on multi-photon blockade [46–48]. The photon pairs generated by the blockade are often required for generation of the entangled light. Two-photon blockade in the strong-coupling qubit-cavity system is investigated using a modified Lindblad master equation [47]. Two- and three- photon blockade phenomena observed respectively in the standard cavity QED and circuit QED systems are attributed to photon-induced tunneling [46]. It is demonstrated that multi-photon blockade may inhibit additional photon being absorbed, which is caused by the nonlinearity of the Jaynes-Cummings (JC) model [48]. In a word, it makes sense to propose a scheme in the optomechanical system for preparing a non-classical state of the MR via the multi-photon blockade.

In this paper, we propose a scheme for generating the NCS of the phonons via the two-photon blockade and the dark state of the system, then investigate its non-classical properties. A two-level QD is driven by a strong laser and the microcavity is engineered within a photonic band-gap (PBG) material. In the dressed state basis induced by the strong field, resonant excitation of two-photon can be implemented by the weak driving field. In this case, the two-photon blockade can be realized in the optomechanical system by tuning the frequency of the weak driving field. The optomechanical system can evolve into a dark state due to the damping of the microcavity and the decay of the dressed QD at selected frequencies in the PBG material. As a result, the phonon mode of the MR can be prepared in a NCS, which is a non-classical state and possesses the properties of the sub-Poisson statistics. We also demonstrate the Wigner functions (WF) of NCS, which negativity implies its non-classicality.

The paper is organized as follows. In section 2, we show the physical model for the driven QD inside an optomechanical microcavity coupled to a weak laser. The whole system is transformed into the polaron frame, and the master equation of the optomechanical system is obtained using the rotating wave, secular and Born-Markov approximations. The preparation of the NCS and its properties are investigated in section 3. Finally, we draw our conclusions in section 4.

## 2. System and master equations

#### 2.1. Description of the system

The system under consideration consists of a high-Q optomechanical microcavity, which is technically engineered inside a PBG material and coupled with an excitonic two-level QD. The optomechanical microcavity with the frequency *ω _{c}* can be coupled with the MR of frequency

*ω*via a combination of radiation pressure and photostriction [49]. The cavity is also driven by a weak laser field with frequency

_{m}*ω*and amplitude

_{p}*α*, as shown in Fig. 1. The two-level QD is represented by the usual Pauli spin- $\frac{1}{2}$ operators

_{p}*σ*

_{+}and

*σ*

_{−}, which satisfy the commutation relations [

*σ*

_{+},

*σ*

_{−}] =

*σ*and [

_{z}*σ*,

_{z}*σ*

_{±}] = ±2

*σ*

_{±}. It is placed inside the microcavity and coupled to another strong laser field with frequency

*ω*and Rabi frequency Ω. This model can be realized by a optomechanical microcaviy, which includes a defective photonic crystal nanobeam and a MR with the “breathing” mode [50], and a QD ironing the surface of the beam (see Fig. 1). The motion of the MR manifests as a successive of contractions and expansions along the

_{L}*x*axis and is driven by a strong laser. The microcavity is driven another weak laser light. In the rotating frame at the frequency of

*ω*, the total Hamiltonian of the system is given by (setting

_{L}*h̄*= 1 throughout the paper)

*a*(

*a*

^{†}),

*b*(

*b*

^{†}) and ${a}_{\lambda}\left({a}_{\lambda}^{\u2020}\right)$ are the annihilation (creation) operators of the microcavity field, the mechanical mode and the photonic crystal vacuum reservoir, respectively. Δ

*=*

_{j}*ω*−

_{j}*ω*(

_{L}*j*=

*a*,

*c*,

*λ*) are the detunings of the QD’s resonance frequency, of the microcavity-mode frequency and of the crystal reservoir frequency from the driving field frequency, respectively. Under the electric dipole and the rotating-wave approximations, the interactions between the driven QD, the microcavity, the mechanical mode and independent bathes can be expressed as [51]

*=*

_{p}*ω*−

_{p}*ω*is the frequency difference between the weak (

_{L}*ω*) and strong (

_{p}*ω*) driving field. The coefficients

_{L}*g*

_{1}and

*g*

_{2}describe the coupling strengths of the microcavity with the QD and phonon mode of the MR. We have assumed that the coupling coefficient

*g*

_{1}is independent on frequency.

*g*

_{2}=

*ω*

_{c}x_{0}/

*L*in Eq. (3) is the single-photon coupling strength, where ${x}_{0}=\sqrt{\frac{1}{2M{\omega}_{m}}}$ is the zero-point fluctuation of the MR with its mass

*M*, and microcavitys rest length is

*L*.

In Eq. (3), the information about the frequency dependent of the photonic crystal is stored in the coupling strength *g _{λ}* (

*ω*) which can be written as

_{λ}*g*(

_{λ}*ω*) =

_{λ}*g*(

_{λ}D*ω*) with

_{λ}*g*being a constant proportional to the dipole moment of the QD and

_{λ}*D*(

*ω*) corresponding to the transfer function of the photonic crystal reservoir. Experimentally, when the length of the waveguide in a 3D heterostructure is 12

_{λ}*μm*, the local photon density of states (LDOS) is demonstrated to have a jump of an order of 100 within a range of relative frequencies about 10

^{−4}

*ω*, where

_{b}*ω*is the cutoff frequency of the waveguide mode [52]. The large discontinuity with step-shaped LDOS can be provided by a waveguide cutoff mode within a 3D-2D PGB heterostructure [52,53]. Besides, a forklike shape LDOS can also be engineered using a trimodal waveguide architecture in a PBG microchip using suitable cutoff frequencies for two of the three waveguide modes [52]. In this paper, the large discontinuity in the LDOS and the relative positions of the relevant frequencies considered in our study are presented in Fig 2, such that |

_{b}*D*(

*ω*)|

_{λ}^{2}=

*u*(

*ω*−

_{λ}*ω*),

_{b}*u*(

*ω*−

_{λ}*ω*) being a unit step function. One valuable feature of waveguide structure is that near the waveguide cutoff frequency, the field is almost the same for rods along the waveguide direction. Therefore, the QD which is close to rods and along the waveguide direction will experience the same field. The cavity resonant frequency can be engineered to be within a PBG or high-LDOS region and near the LDOS discontinuity [52, 54]. Since the QD is driven by a strong laser field, three Mollow bands appear in the fluorescent scattering spectrum [55]. When the Mollow components at different frequencies (especially at lower and higher frequencies) fall within the high LDOS region or the low LDOS region, then they will experience strongly different mode densities [53]. It will strongly influence the dynamics of QD and results in new emission effect, i.e. |

_{b}*D*(

*ω*)|

_{λ}^{2}= 1 for

*ω*>

_{λ}*ω*and |

_{b}*D*(

*ω*)|

_{λ}^{2}= 0 for

*ω*<

_{λ}*ω*. This technique supplies a way for us to eliminate spontaneous transitions at selected frequencies of the QD-cavity system. In our scheme, this property of the PBG material is very important to generate the NCS of the MR, because the generation of the dark state of the system requires the elimination of the spontaneous transitions at both lower and central frequencies of the QD-cavity system, as will be discussed later.

_{b}The strong driving field can be viewed as a dressing field for the QD. Therefore, we diagonalize the first two terms of the Hamiltonian *H*_{0}, i.e., the Hamiltonian of the QD and the coherent interaction with the strong driving field, to find the eigenstates of the coherently driven QD [56]

*H*

_{0}and the interaction Hamiltonian

*H*

_{I}in Eq. (1) are rewritten as

*R*= |

_{ij}*i*〉 〈

*j*|(

*i*,

*j*= +, −) are the dressed-state transition operators,

*R*

_{3}=

*R*

_{++}−

*R*

_{−−}is the population difference operator.

In order to conveniently study the dynamics of the system including QD, microcavity and the MR, we transform the total system into the polaron frame [57] and introduce a renormalized MR-microcavity coupling strength. The polaron transformation of the total system can be written as

where*S*=

*ηa*

^{†}

*a*(

*b*

^{†}−

*b*), and

*η*=

*g*

_{2}/

*ω*. The Lamb-Dicke parameter

_{m}*η*measures the localization of the motional ground state of the MR relative to the (effective) wavelength of the transition. After the transformation, the Hamiltonian of the whole system can be written as with

*H′*

_{I1}in Eq. (10) represents the interaction among the QD, the microcavity and the MR,

*H′*

_{I2}characterizes the interaciton of the microcavity with the MR and the weak laser field. It is indicated from Eq. (11) that the absorption processes of a cavity photon are inseparable from the transition processes of the QD and the generation or annihilation of phonons. In the following section, this will be employed to generate the NCS of the MR. The term

*H′*

_{I3}is the interaction of the QD with the photonic crystal reservoir.

#### 2.2. Derivation of the effective master equation

By use of the Baker-Campbell-Hausdorff Relation [58, 59], the exponential *e*^{−η(b−b†)} can be expanded as

*= 2Ω̄ and*

_{c}*ω*≈ 4Ω̄. With this condition, in the interaction picture defined by

_{m}*H′*

_{0}, we employ the rotating wave approximation to neglect the terms which oscillate rapidly and keep the nearly resonant oscillating terms. The term

*H′*

_{I1}in Eq. (10) is given by

*δ*= Δ

_{m}*+ 2Ω̄ −*

_{c}*ω*is the detuning of the MR with the microcavity and the Rabi frequency of the strong driving field, which in our concerned case satisfies the condition

_{m}*δ*≪ Ω̄

_{m}*,*

*ω*. The operator-valued function

_{m}*f*

_{1}(

*b*

^{†}

*b*) in Eq (14) is of the form

*H′*

_{I2}in Eq. (11) is of form

*δ*= Δ

_{cp}*− Δ*

_{c}*is the detuning of the microcavity from the weak laser field. In Eq. (16), we have used the condition*

_{p}*δ*≪

_{cp}*ω*to employ the rotating wave approximation. This condition will be discuss later. Another operator-valued function

_{m}*f*

_{2}(

*b*

^{†}

*b*) in Eq. (16) is

*H′*

_{I3}is of the form

As mentioned above, in order to explore the dynamics of the system with a three-particle direct interaction, we find that there is one time-independent Hamiltonian in *H̃*_{I1}(*t*), i.e., *H _{J}* =

*g′*

_{1}

*c*

^{2}(

*a*

^{†}

*R*

_{−+}+

*aR*

_{+−}), which corresponds to the Jaynes-Cummings(JC) interaction Hamiltonian describing the interaction between the dressed QD and the microcavity. The JC model has been known to describe an extensive range of coherent phenomena including the spontaneous collapse, the periodic Rabi-flopping and revivals of a two-level atom in a resonant cavity [60]. The non-classical character is investigated by the transmitted light in a driven cavity-atom system both in theoretical [61] and experimental [41] literatures. It is found that remarkable nonlinearities occurs when the coupling between the cavity mode and the atom is strong.

*H*can be diagonalized, and the eigenstates consist of a ground state |0〉 = |−, 0〉 with energy

_{J}*E*

_{0}= 0 and excite states $|{E}_{\pm}^{(n)}\u3009=\frac{1}{\sqrt{2}}\left(|-,n\u3009\pm |+,n-1\u3009\right)$ (

*n*= 1, 2,···) with corresponding energies ${E}_{\pm}^{(n)}=\pm \sqrt{n}{g}_{1}^{\prime}{c}^{2}$. |+(−),

*n*〉 denotes the higher- (lower-) energy eigenstate with

*n*excitations. The splitting between the energy eigenstates in each manifold has a nonlinear dependence on

*n*. We assume that the coupling coefficient

*g′*

_{1}

*c*

^{2}is much larger than the damping rate

*κ*of the microcavity. By cooperatively modulating the frequencies

*ω*,

_{p}*ω*of the driven fields and the strength Ω, two-photon transition processes can be tuned to occur resonantly and the other excitations are far off resonance. In this case, one can truncate the Hilbert space and describe the QD-cavity subsystem with the following four states

_{L}For the convenience of expression, in Eq. (19) we have denoted |0〉,
$|{E}_{-}^{(1)}\u3009$,
$|{E}_{+}^{(1)}\u3009$ and
$|{E}_{-}^{(2)}\u3009$ as |0〉, |1〉, |2〉 and |3〉 respectively, with corresponding energies *E*_{0} = 0, *E*_{1} = −*g′*_{1}*c*^{2}, *E*_{2} = *g′*_{1}*c*^{2} and
${E}_{3}=-\sqrt{2}{g}_{1}^{\prime}{c}^{2}$. In our scheme, the decay rate *κ* and *γ* of the coupled QD-cavity system are much larger than that of the MR. We declare here that the four states concerned in Eq. (19) are established in the subspace of the reduced QD-cavity subsystem, which involves the strong driving field *ω _{L}*. Moreover, the resonant frequency of the cavity

*ω*is at least a factor of 10

_{c}^{3}with the frequency of

*ω*[8–10,50,62–66]. In this situation, the energy interval between the subspaces of the reduced QD-cavity system is much larger than the energy interval of the MR Fock space. Therefore, we suppose that these above four states are enough to describe the considered system even when taking 15 Fock states of the MR into account as shall be discussed in section 3. The energy diagram of the above states of QD-cavity subsystem is shown in Fig 3. As a result, the time-dependent Hamiltonian

_{m}*H*are diagonalized as

_{J}In this new basis of Eq. (19), we can rewrite the photon creation operator of the cavity and the transition operator of the QD as

*H̃*

_{I1},

*H̃*

_{I2}and

*H̃*

_{I3}. After that, we work in the interaction picture defined by

*H*of Eq. (20) and obtain the interaction Hamiltonian of the system in the truncated Hilbert space as where

_{J}*V*

_{2}in Eq. (22) is of the form

_{1}=

*δ*−

_{cp}*g′*

_{1}

*c*

^{2}, Δ

_{2}=

*δ*+

_{cp}*g′*

_{1}

*c*

^{2}, ${\mathrm{\Delta}}_{3}={\delta}_{\mathit{cp}}-\left(\sqrt{2}-1\right){g}_{1}^{\prime}{c}^{2}$, and ${\mathrm{\Delta}}_{4}={\delta}_{\mathit{cp}}-\left(\sqrt{2}+1\right){g}_{1}^{\prime}{c}^{2}$.

*V*

_{2}represents the interaction with larger detuning between the microcavity and the MR when we assume the detuning as ${\delta}_{\mathit{cp}}=\frac{\sqrt{2}}{2}{g}_{1}^{\prime}{c}^{2}\gg {\alpha}_{p}$. It is obvious that

*δ*and

_{m}*δ*are the same order of magnitude. As stated before, since we are only interested in the case Δ

_{cp}*= 2Ω̄ and*

_{c}*ω*≈ 4Ω̄, the parameter

_{m}*δ*= Δ

_{m}*+ 2Ω̄ −*

_{c}*ω*satisfies the condition

_{m}*δ*≪ Ω̄

_{m}*,*

*ω*. Consequently, the condition

_{m}*δ*≪

_{cp}*ω*is satisfied automatically, which is the condition we have used during the calculation of Eq. (16). The two states |1〉 and |2〉 in the Eq. (24) are intermediate states for the two-photon transition and they are largely detuned from the weak probe field. Therefore, we choose ${\delta}_{\mathit{cp}}=\frac{\sqrt{2}}{2}{g}_{1}^{\prime}{c}^{2}$ to adiabatically eliminate |1〉 and |2〉 and calculate the effective Hamiltonian of

_{m}*V*

_{2}as In the derivation of the effective Hamiltonian, we also neglect some constant terms which correspond to the energy shift of the system. As a result, the effective Hamiltonian

*V′*

_{2}in the interaction picture takes the form

*ω*from the weak field are present, as shown in Fig 3. It also should be mentioned that once this subsystem absorbs two photons with frequency

_{p}*ω*, the energy level is shifted far-off resonance with two-photon resonance frequency. It prevents another two photons from entering the microcavity until the previous two photons leave, which is known as two-photon blockade.

_{p}Taking the similar processes, we now calculate the dissipative terms. For the dissipative contribution originating from the interaction of the QD with the photonic crystal reservoir as depicted by *H̃*_{I3} in Eq. (18), we assume that the coupling between the QD and photonic crystal vacuum reservoir is weak and its change on the reservoir is negligible. Also, a fast time scale for the decay of the reservoir correlation is assumed such that the secular [67] and Born-Markov [58,59] approximations can be used. In quantum optics, these standard approximations have been applied with considerable success to a great deal of atom-field interaction problems in photonic crystals [55, 68]. A similar method with details has been introduced in other literature [69]. We apply the second-order perturbation theory to trace over the reservoir degrees of freedom from the following equation

_{R}denotes tracing over the photonic crystal vacuum reservoir variables, and

*H̃*(

_{int}*t*) is the interaction of the QD with photonic crystal in the interaction picture, i.e.,

*H̃*(

_{int}*t*) =

*H̃*

_{I3}(

*t*). Under the Born-Markov approximation, the operator

*ρ*(

_{T}*t′*) in Eq. (27) can be replaced by

*ρ*(

*t*) ⊗

*ρ*(0), where

_{λ}*ρ*(

*t*) is the reduced density operator for the three-particle (QD, cavity and MR) system and

*ρ*(0) is the initial density operator of the photonic crystal vacuum reservoir. After straightforward calculations, we obtain the dissipative term

_{λ}*ℒ*induced by the interaction of the QD with photonic crystal vacuum reservoir, which has Liouvillian form as

_{λ}ρ*γ*

_{0,±}in Eq. (28) are responsible for the damping among the dressed states of the dressed-QD system. The parameter ${\gamma}_{0}=2\pi {c}^{2}{s}^{2}{\sum}_{\lambda}{g}_{\lambda}^{2}{\left|D({\omega}_{\lambda})\right|}^{2}\delta ({\omega}_{\lambda}-{\omega}_{L})$ is the spontaneous emission rate with which the electron of the QD transits at central frequency of the dressed states. The parameter ${\gamma}_{+}=2\pi {c}^{4}{\sum}_{\lambda}{g}_{\lambda}^{2}{\left|D({\omega}_{\lambda})\right|}^{2}\delta ({\omega}_{\lambda}-{\omega}_{L}-2\overline{\mathrm{\Omega}})$ represents spontaneous emission rate of electron which occurs at high frequency

*ω*

_{+}=

*ω*+ 2Ω̄ of the dressed states from the upper dressed state |+〉 of one manifold to the lower dressed state |−〉 of the manifold below, whereas the parameter ${\gamma}_{-}=2\pi {s}^{4}{\sum}_{\lambda}{g}_{\lambda}^{2}{\left|D({\omega}_{\lambda})\right|}^{2}\delta ({\omega}_{{\lambda}_{-}}-{\omega}_{L}+2\overline{\mathrm{\Omega}})$ describes spontaneous emission which occurs at lower frequency

_{L}*ω*

_{−}=

*ω*− 2Ω̄ of the dressed states from the lower dressed state of one manifold to the upper dressed state of the manifold below. However, in our paper, the microcavity is engineered in a PBG material. As mentioned before, by appropriately engineering the system topology, the QD coupled to the cavity field can experiences the discontinuity of the LDOS [54]. It strongly influences the dynamics of QD and results in new emission effects, i.e. |

_{L}*D*(

*ω*)|

_{λ}^{2}=

*u*(

*ω*−

_{λ}*ω*) = 0 for

_{b}*ω*<

_{λ}*ω*and |

_{b}*D*(

*ω*)|

_{λ}^{2}= 1 for

*ω*>

_{λ}*ω*. For simplicity, we also assume that the photonic density of modes is constant over the spectral regions surrounding the dressed-state resonant frequencies

_{b}*ω*and

_{L}*ω*± 2Ω. In this way, the spontaneous transitions of the QD at frequencies

_{L}*ω*

_{−}and

*ω*can be selectively eliminated, i.e.,

_{L}*γ*

_{0}=

*γ*

_{−}= 0.

The decay of the cavity *ℒ _{c}ρ* =

*κ𝒟*[

*a*] is phenomenologically added with the definition $\mathcal{D}[\mathcal{O}]=\frac{1}{2}\left(2\mathcal{O}\rho {\mathcal{O}}^{\u2020}-\rho {\mathcal{O}}^{\u2020}\mathcal{O}-{\mathcal{O}}^{\u2020}\mathcal{O}\rho \right)$.

*κ*is the decay rate of the microcavity which may lead to a trend towards vacuum mode of the cavity. In fact,

*ℒ*=

_{c}ρ*κ𝒟*[

*a*] can also be obtained following above procedures of

*ℒ*. After the polaron transformation of Eq. (7), the decay of the cavity is of the form

_{λ}ρ*ℒ*and

_{c}ρ*ℒ*into the truncated subspace of the dressed QD-cavity system as illustrated in Eq. (19), we finally obtain the master equation for the density operator

_{λ}ρ*ρ*of the three-particle system as the form

*ℒ*and

_{λ}ρ*ℒ*in Eq. (30) being the form

_{c}ρ*η*

^{m+l}(−

*η*)

^{n+k}/(

*m*!

*n*!

*k*!

*l*!), ${B}_{j}=\frac{3}{2}-\sqrt{2}{\left(-1\right)}^{j}$. In the derivations of

*ℒ*and

_{c}ρ*ℒ*, the intermediate states |1〉 and |2〉 (see figure 3) are not abandoned because they are sometimes occupied during the dissipation processes.

_{λ}ρThe purpose of this paper is to prepare the NCS of the MR’s phonon mode. In the next section, we will demonstrate the generation of the NCS applying the dark state of the total system. The properties of this non-classical state is also explored.

## 3. Preparation of NCS using dark state

Up to now, we have obtained the effective master equation which is able to describe the dynamics of three-particle system including the MR, the microcavity and the dressed QD. We also mentioned the occurrence of two-photon blockade. However, another significant object of this paper is to prepare the NCS of the phonon. In the following, on the basis of the effective Hamiltonian, we investigate the dark state of the system and propose a scheme to generate the NCS of the MR. Some properties of the NCS are also discussed.

We first consider the case that the interaction Hamiltonian *V*_{1} and *V′*_{2} can be neglected, i.e., *V*_{1} = *V′*_{2} = 0. In this case, the dynamic of the whole system is determined by the dissipation term *ℒρ* = *ℒ _{λ}ρ* +

*ℒ*. The steady-state behavior of the system can be obtained via the steady equation

_{c}ρ*ℒρ*(∞) = 0 with $\rho (\infty )={\sum}_{j=0}^{3}{\rho}_{jj}|j\u3009\hspace{0.17em}\u3008j|\otimes {\rho}_{bj}(\infty )$, where

*ρ*(∞) is the reduced density operator of the MR in the steady state when the QD-cavity subsystem is in state |

_{bj}*j*〉.

*ρ*= 〈

_{jj}*j*|

*ρ*|

*j*〉 are elements of the reduced density matrix of the QD-cavity subsystem. In the steady-state solution, we find that

*ρ*

_{00}= [

*γ*

_{+}/(

*γ*

_{+}+

*γ*

_{−})]

^{2}. It should be stressed here that if the spontaneous emission rate can be eliminated

*γ*

_{−}= 0, the reduced QD-cavity subsystem evolves into the state |0〉 = |−, 0〉 as depicted in Eq. (19), i.e.,

*ρ*

_{00}= 1,

*ρ*= 0, (

_{jj}*j*= 1, 2, 3). This means that the microcavity is in the vacuum state and the QD is in the lower dressed state |−〉, corresponding to

*ρ*(∞) = |0〉 〈0| ⊗

*ρ*(∞). In fact, this condition

_{b}*γ*

_{−}= 0 can be realized by technically engineering the microcavity into the PBG material. As mention before, by properly engineering the system topology, the QD coupled to the microcavity can experiences a discontinuity of the LDOS. It strongly influences the dynamics of QD and results in new emission effects, i.e. |

*D*(

*ω*)|

_{λ}^{2}= 0 for

*ω*<

_{λ}*ω*and |

_{b}*D*(

*ω*)|

_{λ}^{2}= 1 for

*ω*>

_{λ}*ω*. In this way, the spontaneous emissions occurring at frequencies

_{b}*ω*

_{−}and

*ω*are selectively eliminated, i.e., both the spontaneous emission rates

_{L}*γ*

_{−}and

*γ*

_{0}are zero. Only the spontaneous emission occurring at the frequency

*ω*

_{+}=

*ω*+ 2Ω̄ is allowed to happen with the spontaneous emission rate

_{L}*γ*

_{+}. In our scheme, a nonlinear coherent state of the MR mode is generated based on the presence of a dark state. Additionally, if the PBG material were to be moved such that none of the spontaneous transitions are eliminated, the steady state of the system can not evolve into such a dark state, then the scheme for generation of the target state may be the other project which is out of our interest.

When taking the couplings of the microcavity with the QD and the weak laser field into account, i.e., *V*_{1} ≠ 0, *V′*_{2} ≠ 0, we devote to finding the solution which simultaneously satisfy the equations −*i*[*V*_{1} + *V′*_{2}, *ρ*] = 0 and *ℒρ* = 0. Because the solution of the latter one is that in which the QD-cavity subsystem is in the state |0〉 as discussed above. Therefore the steady state of the whole system can be described by |Ψ〉 = |0, *φ _{b}*〉, (

*ρ*= |0,

*φ*〉 〈0,

_{b}*φ*|), where |

_{b}*φ*〉 is the motional state of the mechanical mode. When the total system is in such a state, the QD-microcavity system stops to fluoresce and the motional state |

_{b}*φ*〉 can be obtained. Therefore, it can be said that the system is in a motional dark state [21], which has various types. This target state can be obtained from equation

_{b}*φ*〉 satisfying We call the eigenstates of Eq. (33) with the boson annihilation operator

_{b}*b*and nonlinear functions of number operator

*N̂*=

*b*

^{†}

*b*as NCSs [21]. At this point, we can see that the state |

*φ*〉 is a NCS of the MR which is obtained by using the dark state properties of the system.

_{b}In order to explore the properties of |*φ _{b}*〉, we expand the state |

*φ*〉 as

_{b}*c*is the probability amplitude of the Fock state |

_{n}*n*〉. In the Fock representation, we use the definitions of

*f*(

_{i}*b*

^{†}

*b*), (

*i*= 1, 2) in Eq. (15) and Eq. (17) to expand Eq. (33). After some straightforward calculations, we find that the coefficients of the Fock states satisfy a recursion relation as

*η*

_{0}such that ${L}_{n}^{(0)}\left({\eta}^{2}\right)=0$ for

*n*=

*N*. 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

*c*

_{n+1}zero. This recursion relation ensures that all coefficients

*c*vanish for

_{n}*n*>

*N*. Thus, there will be none excitations from

*N*phonon states to

*N*+ 1 phonon states.

So far, by designing a suitable coherent Hamiltonian and the interaction between the phonon bath and the coupled QD-cavity system, the mechanical mode can be prepared to the NCS |*φ _{b}*〉 which is a non-classical motional dark state. It should be highlighted here that the prepared state |

*φ*〉 of the mechanical mode is completely decoupled from the QD and the cavity. This can be understood in the polaron frame described by Eq. (7), because the whole system can be described by |Ψ〉 = |0

_{b}*, φ*〉, where the state |

_{b}*φ*〉 of mechanical mode is decoupled from the QD and the cavity. However, when going back from this polaron frame, we obtain

_{b}*S*is defined below Eq. (7), i.e.

*S*=

*ηa*

^{†}

*a*(

*b*

^{†}−

*b*). |0〉 = |−, 0〉 is the dark state of the reduced QD-cavity subsystem, in which the cavity is in the vacuum state and the QD is in the lower dressed state |−〉. As a result, |Ψ〉 = |0,

*φ*〉 keeps the same form under the inverse polaron transformation. Therefore, the prepared NCS of the MR is completely decoupled from the QD and the cavity.

_{b}Recently, when a two-level system and a mechanical mode are both strongly coupled to a cavity, the preparation of the mechanical modes in a non-Gaussian dark state has also been studied [70]. Employing additional auxiliary cavity mode to construct two Liouvillians, an approximated entangled state are prepared in the weak coupling limit, i.e., *ζ* = (*g _{q}/g_{m}*)

^{2}≤ 1, where

*g*and

_{q}*g*are coupling strengthes of the cavity with a two-level system and a mechanical mode respectively. When ignoring the spontaneous emission of the two-level system and measuring the two-level system, the mechanical oscillator can be prepared in a superposition of Fock states with a fidelity of

_{m}*F*≈ 0.83. However, it should be emphasized that the literature [70] is distinct from this paper. In this paper, when properly tuning the frequency of the weak driving field, two-photon blockade phenomenon occurs. The QD-cavity subsystem can evolve into a dark state due to the damping of the microcavity and the elimination of the decay rate of the QD at selected frequencies in the photonic band-gap (PBG) material. By exploiting the interaction between the phonon bath and the coupled QD-cavity system, a suitable coherent Hamiltonian is designed so that the mechanical mode can naturally be prepared into a NCS. The prepared NCS of the MR in this paper is a definite pure state. Besides, the prepared NCS is completely decoupled from the QD and the cavity. The intrinsic non-linearity and non-classical behavior of the state will be analysed via phonon number statistics, second-order correlation function, and Wigner function.

#### 3.1. Phonon number statistics

The phonon statistics of a MR are the “finger-print” of its quantum state, from which a number of useful measures of non-classicality may be inferred. The phonon number distribution (PND) of the dark states is complicated by analytical ways with the correlated polynomials. In this paper, we obtain the results by numerically solving the equation involving Laguerre polynomial. For simplicity, we factitiously choose the parameter *η* = *g*_{2}/*ω _{m}* = 0.305463 to satisfy the equation

*N*> 15, the transition coefficients to the

*N*+ 1 phonon number states vanish.

We should note that *η* = *g*_{2}/*ω _{m}* = 0.305463 means the system approaches the single-photon strong coupling regime. If the single-photon optomechanical coupling strength

*g*

_{0}is an appreciable fraction of the mechanical resonance frequency

*ω*and the combined nonlinearity parameter satisfies the relation ${g}_{0}^{2}/\kappa {\omega}_{m}\ge 1$, the system is in the single-photon ultra-strong coupling regime characterized by cavity photon blockade [71, 72]. This challenging idea has inspired emerging attentions from experimental field. The nonlinearity of the phonon mode of the MR originating from the strong coupling between the microcavity and the MR has been discussed based on capacitive [73] or magnetic [74] coupling in the microwave-regime optomechanical devices. By enhancing the effect of a single-photon in the optical cavity, the strong-coupling regime is reached with the auxiliary of the squeezing of a cavity mode [4]. A single-photon ultra-strong coupling regime is proposed to be reached in an optomechanical system [75]. For a system with a quantum dot strongly coupled to a photonic-crystal nanocavity, J. Vučković et. al have studied non-classical higher-order photon correlations [76]. The strong coupling of a single two-level solid-state system with a photon is realized with a single quantum dot in a semiconductor microcavity [77]. Here, we supply some main parameters to implement our scheme with

_{m}*ω*/2

_{c}*π*= 195THz and

*ω*/2

_{m}*π*= 3.68GHz,

*κ*/2

*π*= 150MHz,

*g*

_{1}= 5

*κ*, which are similar to that in the literature [10]. The condition

*η*=

*g*

_{2}/

*ω*= 0.3 indicates that a value of

_{m}*g*

_{2}/2

*π*= 1.1GHz is needed to be achieved. At present, the optomechanical coupling strength

*g*

_{2}has realized

*g*

_{2}/2

*π*= 11.5MHz, and is predicted to approach

*g*

_{2}/2

*π*= 26MHz [65]. It is still inadequate for realizing our proposal. However, very recently, it has been proposed that applying a large-amplitude, strongly detuned mechanical parametric drive to amplify mechanical zero-point fluctuations, the radiation pressure interaction is enhanced and thus the overall coupling strength is dramatically increased, allowing us to completely enter the quantum-enabled, strong-coupling regime [78]. The strong radiation pressure forces from squeezed light on a mechanical oscillator is observed [79]. These inspiring and promising progresses enable us to believe that our scheme can be achieved in the near future.

In Eq. (37), the fact *c*_{n+1} = 0 for *n* > *N* = 15 indicates that the MR is trapped into the Fock superposition state via the cavity damping and selectively elimination the spontaneous emission *γ*_{−} in the PBG material. The specific expression of *c _{n}* is

*χ*. The PND is defined as where |

*n*〉 is the eigenstate of the phonon number operator

*b*

^{†}

*b*. The probability of finding

*n*phonons in the field is given by

*P*(

*n*) = |

*c*|

_{n}^{2}. Note that for the NCS, the probability of the phonon number depends on the nonlinearity functions and varies with parameter

*α*

_{p}/g_{1}. In figure 4, we show the PND of the MR as a function of

*n*for the different values of

*α*

_{p}/g_{1}. From this figure it is seen that the probability is nearly close to zero for

*n*≥ 7. As apparently shown in the figure 4, the peak of

*P*depends on the parameter

_{n}*α*

_{p}/g_{1}, and which position moves forward to the larger number distribution with increasing of the parameter

*α*

_{p}/g_{1}, meanwhile value of the peak is reduced. With the increasing of the parameter

*α*

_{p}/g_{1}, the mean phonon number in the NCS may increase, which means that multi-phonon processes are enhanced.

#### 3.2. Second-order correlation function(SOCF)

In order to extract more information about the NCS of the MR, we calculate the corresponding second-order correlation function(SOCF) which is correlated with the coherence and statistical properties a field. The normalized SOCF of the phonon mode of the MR in the steady state can be defined as [80]

*s*describes the steady state of the field and

*b*

^{†}(

*b*) is the creation (annihilation) operator of the MR. This function can give the influences of detecting one phonon of the source on the result of detecting another one. The SOCF

*g*

^{(2)}(0) can be applied to distinguish whether the statistical properties of the field is super-Poissonian (

*g*

^{(2)}(0) > 1), Poissonian (

*g*

^{(2)}(0) = 1), or sub-Poissonian (

*g*

^{(2)}(0) < 1). The sub-Poissonian statistics is usually correlated to non-classical state [81, 82].

In order to get the insight into the non-classical properties of the phonon, we numerically calculate the SOCF in the Fock state representation with the following expression

*α*

_{p}/g_{1}in figure 5. It is clearly shown from figure 5 that the value of SOCF

*g*

^{(2)}(0) decreases with the increasing of

*α*

_{p}/g_{1}and approaches the vicinity of

*g*

^{(2)}(0) < 1 when the parameter

*α*

_{p}/g_{1}approximately takes the value 0.1225. The fact

*g*

^{(2)}(0) < 1 indicates the sub-Poissonian statistics and the non-classicality of the generated phonon state. It also implies that the NCS of the MR |

*φ*〉 exhibits the sub-Poissonian behaviour and non-classical effect.

_{b}#### 3.3. Wigner function

The Wigner function (WF) [83] was first introduced by Wigner in 1932 to calculate quantum correction to a classic distribution function of a quantum-mechanical system. It is a quasi-probability distribution which fully describes the states of a quantum system in phase space, (either the position-momentum space for an harmonic oscillator or, equivalently, the space spanned by two orthogonal quadratures of the electromagnetic fields for a single-mode state of light), in the same fashion as a probability distribution (non-negative by definition) characterizes a classical system. It has now become a useful tool to study the non-classical properties of quantum states, because the partial negativity of the WF [83,84] is indeed a good indication of the highly non-classical character of the state. In order to further unveil the non-classicality of the NCS of phonon, in this subsection, we supply the expression of WF for the phonon in Fock state representation and numerically illustrate the corresponding property.

For a single-mode system, the WF of the density operator *ρ* in Fock state representation |*n*〉 is expressed as

*ρ*is the reduced density operator of the MR and

*D*(

*β*) = exp(

*βb*

^{†}−

*β*

^{*}

*b*) is the displacement operator. The WF of the NCS |

*φ*〉 is depicted in figure 6 in phase space, where

_{b}*X̂*and

*Ŷ*correspond to the position and the momentum respectively. The distribution is suppressed along

*X̂*quadrature component and presents crescent shape. According to Hudson’s theorem [85], the presence of negative regions in the Wigner distribution is an authentic sign of quantum character. Obviously, there is a negative region in figure 6, which reveals the strong non-classicality and genuine quantum behavior.

As mentioned at the beginning of this subsection, the partial negativity of the WF [83, 84] is a good indication of the highly non-classical character of a quantum state. In order to quantitatively explore the non-classicality of the state, we study the absolute value *P _{W}* of the total negative probability of the WF, which can be defined by [84, 86]

*ϖ*is the negative region of WF distribution. Since the negative region

*W*(

*x*,

*y*) of WF has been already obtained. Then, the value of

*P*can be numerically calculated according to the above definition. One should note in Eq. (43) that only the negative value of

_{W}*W*(

*x*,

*y*) is integrated. During the calculation, the value of

*P*is set to be zero when

_{W}*W*(

*x*,

*y*) > 0.

We illustrate *P _{W}* of the total negative probability of the WF as a function of the parameter

*α*

_{p}/g_{1}in figure 7. It can be seen from figure 7 that with the increasing of

*α*

_{p}/g_{1},

*P*takes a definite value when the parameter

_{W}*α*

_{p}/g_{1}approaches the region

*α*

_{p}/g_{1}> 0.09. The peak of

*P*occurs at

_{W}*α*

_{p}/g_{1}≈ 0.145. Whereas, as found from figure 5 that

*g*

^{(2)}(0) takes a smaller value than 1 in the region

*α*

_{p}/g_{1}≥ 0.1225. In this region 0.09 <

*α*

_{p}/g_{1}< 0.1225,

*g*

^{(2)}(0) takes a larger value than 1, being unable to tell whether or not the non-classical effect appears. But in the same region

*P*takes a definite value, which is the indication of non-classicality. The region for

_{W}*g*

^{(2)}(0) < 1 is narrower than that for

*P*> 0. In other words, the WF of the phonon state has a negative region and the state possesses the non-classicality, but it does not show the sub-poissonian phonon statistics [87]. Therefore, we emphasize here that for the prepared NCS of the phonon mode in this paper, the absolute value

_{W}*P*may be a better choice for quantifying the non-classicality than the SOCF.

_{W}Up to now, we have prepared the NCS of the phonon mode of MR by using the dark state of the system. Both the facts that *g*^{(2)}(0) < 1 and the WF approaches a negative region are convincing indications of the sub-Poissonian behaviour and non-classicality of the NCS. However, we should point out here that in all the calculations, the damping of the MR is neglected. Because for a MR with high quality factor, the mechanical damping *γ _{m}* could be much smaller than the cavity dissipation rate

*κ*. The mechanical damping can be neglected during the time needed to achieve the steady states. As mentioned before, in our scheme, the four states concerned in Eq. (19) are established in the subspace of the reduced QD-cavity subsystem, which involves the strong driving field

*ω*. The resonant frequency of the cavity

_{L}*ω*is at least 10

_{c}^{3}times the frequency of

*ω*[8–10, 50, 62–66]. In this situation, the energy interval between the subspaces of the reduced QD-cavity system is much larger than the energy interval of the MR Fock space. Therefore, the four states of the QD-cavity subsystem and 15 Fock states of the MR are valid to describe the considered system. Here, we should mention that the increasing of

_{m}*ω*corresponds to a larger quality factor, which is estimated to 1.1 × 10

_{m}^{6}when

*ω*= 2

_{m}*π*× 5.3GHz at cryogenic temperature [50]. A larger

*ω*than the decay of the cavity

_{m}*κ*means the appearance of the sideband excitations. Then, the sideband cooling method can be employed. For a given cavity matched with the QD, the increasing of

*ω*may dissatisfy the condition ${g}_{0}^{2}/\kappa {\omega}_{m}\ge 1$ thus lead to difficult achievement of strong coupling regime, where

_{m}*g*

_{0}is the optomechanical coupling strength. In this case, as discussed before, it has been proposed recently that applying a large-amplitude, strongly detuned mechanical parametric drive to amplify mechanical zero-point fluctuations, the radiation pressure interaction is enhanced and thus the overall coupling strength is dramatically increased, allowing one to completely enter the quantum-enabled, strong-coupling regime [78]. The strong radiation pressure forces from squeezed light on a mechanical oscillator is observed [79]. These inspiring and promising progresses enable us to believe that our scheme can be achieved in the near future.

As is known that non-Gaussian states might be more appropriate than Gaussian state in the sense of the robustness in the applications of teleportation and the cloning [88–90] of quantum states. The nonlinear coherent state generated in the solid-state optomechanical system may benefit for the optomechanical quantum information processing with photons and phonons [36].

## 4. Conclusion

In this paper, we have proposed a scheme to prepare a NCS for the MR by using the dark state of the optomechanical system, in which a two-level QD located inside a microscopic cavity is engineered in a PBG material. The effective master equation of the optomechanical system is obtained in the polaron frame. Two-photon transition can be realized in the dressed states of the QD-microcavity subsystem via adjusting the discrepancy of the frequencies between the micro-cavity and the weak laser field. By exploiting the interaction between the phonon bath and the coupled QD-cavity system, a suitable coherent Hamiltonian is designed so that the mechanical mode can naturally be prepared into a NCS, which is a definite pure state and is completely decoupled from the QD and the cavity. The NCS of the phonon exhibits strong non-classicality since its SOCF shows the sub-Poissonian distribution and its Wigner function approaches the negative region. Our proposal provides a way to prepare a NCS with nonclassicality in an optomechanical crystal system. We hope this could stimulate discussions about the applications of non-Gaussian non-classical states in the quantum information processing [36] and quantum teleportation [91] in the optomechanical systems. The NCS could contribute to the application in optimizing quantum information in the near future, which is important in binary (or multibinary) communication [27].

## Acknowledgments

This work was supported by the National Natural Science Foundation of China (NSFC) (Grants No. 61275123, No.11474119, No.11404103 and No. 11304024) and the National Basic Research Program of China (Grant No. 2012CB921602).

## References and links

**1. **A. H. Safavi-Naeini, T. P. M. 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]

**2. **E. Verhagen, S. Deléglise, S. Weis, A. Schliesser, and T. J. Kippenberg, “Quantum-coherent coupling of a mechanical oscillator to an optical cavity mode,” Nature **482**, 63–67 (2012). [CrossRef] [PubMed]

**3. **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]

**4. **X. Y. Lü, Y. Wu, J. R. Johansson, H. Jing, J. Zhang, and F. Nori, “Squeezed optomechanics with phase-matched amplification and dissipation,” Phys. Rev. Lett. **114**, 093602 (2015). [CrossRef] [PubMed]

**5. **Y. C. Liu, Y. F. Xiao, X. S. Luan, and C. W. Wong, “Dynamic dissipative cooling of a mechanical resonator in strong coupling optomechanics,” Phys. Rev. Lett. **110**, 153606 (2013). [CrossRef] [PubMed]

**6. **W. J. Gu and G. X. Li, “Quantum interference effects on ground-state optomechanical cooling,” Phys. Rev. A **87**, 025804 (2013). [CrossRef]

**7. **A. D. O’Connell, M. Hofheinz, M. Ansmann, R. C. Bialczak, M. Lenander, E. Lucero, M. Neeley, D. Sank, H. Wang, M. Weides, J. Wenner, J. M. Martinis, and A. N. Cleland, “Quantum ground state and single-phonon control of a mechanical resonator,” Nature **464**, 697–703 (2010). [CrossRef]

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

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

**10. **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]

**11. **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]

**12. **J. Q. Liao and C. K. Law, “Parametric generation of quadrature squeezing of mirrors in cavity optomechanics,” Phys. Rev. A **83**, 033820 (2011). [CrossRef]

**13. **X. W. Xu, Y. J. Zhao, and Y. X. Liu, “Entangled-state engineering of vibrational modes in a multimembrane optomechanical system,” Phys. Rev. A **88**, 022325 (2013). [CrossRef]

**14. **A. Kronwald, F. Marquardt, and A. A. Clerk, “Arbitrarily large steady-state bosonic squeezing via dissipation,” Phys. Rev. A **88**, 063833 (2013). [CrossRef]

**15. **T. Ramos, V. Sudhir, K. Stannigel, P. Zoller, and T. J. Kippenberg, “Nonlinear quantum optomechanics via individual intrinsic two-level defects,” Phys. Rev. Lett. **110**, 193602 (2013). [CrossRef] [PubMed]

**16. **W. J. Gu, G. X. Li, S. P. Wu, and Y. P. Yang, “Generation of non-classical states of mirror motion in the single-photon strong-coupling regime,” Opt. Express **22**, 18254–18267 (2014). [CrossRef] [PubMed]

**17. **S. Ya. Kilin and A. B. Mikhalychev, “Single-atom laser generates nonlinear coherent states,” Phys. Rev. A **85**, 063817 (2012). [CrossRef]

**18. **H. Fakhri and A. Hashemi, “Nonclassical properties of the q-coherent and q-cat states of the Biedenharn-Macfarlane q oscillator with q>1,” Phys. Rev. A **93**, 013802 (2016). [CrossRef]

**19. **S. Dey, “Q-deformed noncommutative cat states and their nonclassical properties,” Phys. Rev. D **91**, 044024 (2015). [CrossRef]

**20. **S. Dey and A. Fring, “Squeezed coherent states for noncommutative spaces with minimal length uncertainty relations,” Phys. Rev. D **86**, 064038 (2012). [CrossRef]

**21. **R. L. de Matos Filho and W. Vogel, “Nonlinear coherent states,” Phys. Rev. A **54**, 4560 (1996). [CrossRef] [PubMed]

**22. **V. I. Manko, G. Marmo, and E. C. G. Sudarshan, “F-oscillators and nonlinear coherent states,” Phys. Scr. **55**, 528 (1997). [CrossRef]

**23. **W. Vogel and R. L. de Matos Filho, “Nonlinear Jaynes-Cummings dynamics of a trapped ion,” Phys. Rev. A **52**, 4214 (1995). [CrossRef] [PubMed]

**24. **M. H. Naderi, M. Soltanolkotabi, and R. Roknizadeh, “A theoretical scheme for generation of nonlinear coherent states in a micromaser under intensity-dependent Jaynes-Cummings model,” Eur. Phys. J. D **32**, 397–408 (2005). [CrossRef]

**25. **M. B. Harouni, R. Roknizadeh, and M. H. Naderi, “Nonlinear coherent state of an exciton in a wide quantum dot,” J. Phys. B: At. Mol. Opt. Phys. **41**, 225501 (2008). [CrossRef]

**26. **Z. Kis, W. Vogel, and L. Davidovich, “Nonlinear coherent states of trapped-atom motion,” Phys. Rev. A **64**, 033401 (2001). [CrossRef]

**27. **E. M. F. Curado, J-P Gazeau, and Ligia M. C. S. Rodrigues, “Nonlinear coherent states for optimizing quantum information,” Phys. Scr. **82**, 038108 (2010). [CrossRef]

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

**29. **G. Giedke and J. I. Cirac, “Characterization of Gaussian operations and distillation of Gaussian states,” Phys. Rev. A **66**, 032316 (2002). [CrossRef]

**30. **A. P. Lund, T. C. Ralph, and H. L. Haselgrove, “Fault-tolerant linear optical quantum computing with small-amplitude coherent states,” Phys. Rev. Lett. **100**, 030503 (2008). [CrossRef] [PubMed]

**31. **H. Nha and H. J. Carmichael, “Proposed test of quantum nonlocality for continuous variables,” Phys. Rev. Lett. **93**, 020401 (2004). [CrossRef] [PubMed]

**32. **F. L. Kien, S. D. Gupta, V. I. Balykin, and K. Hakuta, “Spontaneous emission of a cesium atom near a nanofiber: efficient coupling of light to guided modes,” Phys. Rev. A **72**, 032509 (2005). [CrossRef]

**33. **J. Niset, J. Fiurasek, and N. J. Cerf, “No-go theorem for Gaussian quantum error correction,” Phys. Rev. Lett. **102**, 120501 (2009). [CrossRef] [PubMed]

**34. **M. G. Genoni and M. G. A. Paris, “Quantifying non-Gaussianity for quantum information,” Phys. Rev. A , **82**, 052341 (2010). [CrossRef]

**35. **G. Adesso, F. DellAnno, S. De Siena, F. Illuminati, and L. A. M. Souza, “Optimal estimation of losses at the ultimate quantum limit with non-Gaussian states,” Phys. Rev. A **79**, 040305 (2009). [CrossRef]

**36. **K. Stannigel, P. Komar, S. J. M. Habraken, S. D. Bennett, M. D. Lukin, P. Zoller, and P. Rabl, “Optomechanical quantum information processing with photons and phonons,” Phys. Rev. Lett. **109**, 013603 (2012). [CrossRef] [PubMed]

**37. **M. Li, H. X. Tang, and M. L. Roukes, “Ultra-sensitive NEMS-based cantilevers for sensing, scanned probe and very high-frequency applications,” Nat. Nanotecnol. **2**, 114–120 (2007). [CrossRef]

**38. **E. Knill, R. Laflamme, and G. J. Milburn, “A scheme for efficient quantum computation with linear optics,” Nature **409**, 46–52 (2001). [CrossRef] [PubMed]

**39. **A. Imamoglu, H. Schmidt, G. Woods, and M. Deutsch, “Strongly interacting photons in a nonlinear cavity,” Phys. Rev. Lett. **79**, 1467 (1997). [CrossRef]

**40. **R Alléaume, F Treussart, G Messin, Y Dumeige, J-F Roch, A Beveratos, R Brouri-Tualle, J-P Poizat, and P Grangier, “Experimental open-air quantum key distribution with a single-photon source,” New J. Phys. **6**, 92 (2004). [CrossRef]

**41. **K. M. Birnbaum, A. Boca, R. Miller, A. D. Boozer, T. E. Northup, and H. J. Kimble, “Photon blockade in an optical cavity with one trapped atom,” Nature **436**, 87–90 (2005). [CrossRef] [PubMed]

**42. **A. Faraon, I. Fushman, D. Englund, N. Stoltz, P. Petroff, and J. Vuckovic, “Coherent generation of non-classical light on a chip via photon-induced tunnelling and blockade,” Nature Phys. **4**, 859–893 (2008). [CrossRef]

**43. **D. Englund, A. Majumdar, A. Faraon, M. Toishi, N. Stoltz, P. Petroff, and J. Vuckovic, “Resonant excitation of a quantum dot strongly coupled to a photonic crystal nanocavity,” Phys. Rev. Lett. **104**, 073904 (2010). [CrossRef] [PubMed]

**44. **C. Lang, D. Bozyigit, C. Eichler, L. Steffen, J. M. Fink, A. A. Abdumalikov, M. Baur, S. Filipp, M. P. da Silva, A. Blais, and A. Wallraff, “Observation of resonant photon blockade at microwave frequencies using correlation function measurements,” Phys. Rev. Lett. **106**, 243601 (2011). [CrossRef] [PubMed]

**45. **A. J. Hoffman, S. J. Srinivasan, S. Schmidt, L. Spietz, J. Aumentado, H. E. Türeci, and A. A. Houck, “Dispersive photon blockade in a superconducting circuit,” Phys. Rev. Lett. **107**, 053602 (2011). [CrossRef] [PubMed]

**46. **A. Miranowicz, M. Paprzycka, Y. X. Liu, J. Bajer, and F. Nori, “Two-photon and three-photon blockades in driven nonlinear systems,” Phys. Rev. A **87**, 023809 (2013). [CrossRef]

**47. **W. W. Deng, G. X. Li, and H. Qin, “Enhancement of the two-photon blockade in a strong-coupling qubit-cavity system,” Phys. Rev. A **91**, 043831 (2015). [CrossRef]

**48. **S. S. Shamailov, A. S. Parkins, M. J. Collett, and H.J. Carmichael, “Multi-photon blockade and dressing of the dressed states,” Opt. Commun. **283**, 766–772 (2010). [CrossRef]

**49. **M. Aspelmeyer, T. J. Kippenberg, and F. Marquardt, “Cavity optomechanics,” Rev. Mod. Phys. **86**, 1391–1452 (2014). [CrossRef]

**50. **R. Riedinger, S. Hong, R. A. Norte, J. A. Slater, J. Shang, A. G. Krause, V. Anant, M. Aspelmeyer, and S. Gröblacher, “Non-classical correlations between single photons and phonons from a mechanical oscillator,” Nature **530**, 313–316 (2016). [CrossRef] [PubMed]

**51. **S. Carlig and M. A. Macovei, “Long-time correlated quantum dynamics of phonon cooling,” Phys. Rev. A **90**, 013817 (2014). [CrossRef]

**52. **R. Wang and S. John, “Engineering the electromagnetic vacuum for controlling light with light in a photonic-band-gap microchip,” Phys. Rev. A **70**, 043805 (2004). [CrossRef]

**53. **X. Ma and S. John, “Switching dynamics and ultrafast inversion control of quantum dots for on-chip optical information processing,” Phys. Rev. A **80**, 063810 (2009). [CrossRef]

**54. **L. Florescu, S. John, T. Quang, and R. Wang, “Theory of a one-atom laser in a photonic band-gap microchip,” Phys. Rev. A **69**, 013816 (2004). [CrossRef]

**55. **S. John and T. Quang, “Spontaneous emission near the edge of a photonic band gap,” Phys. Rev. A **50**, 1764 (1994). [CrossRef] [PubMed]

**56. **C. Cohen-Tannoudji, J. Dupont-Roc, and G. Grynberg, *Atom-photon interactions* (Wiley, 1992).

**57. **C. Roy and S. Hughes, “Phonon-dressed Mollow triplet in the regime of cavity quantum electrodynamics: excitation-induced dephasing and nonperturbative cavity feeding effects,” Phys. Rev. Lett. **106**, 247403 (2011). [CrossRef] [PubMed]

**58. **M. O. Scully and M. S. Zubairy, *Quantum optics* (Cambridge, 1997).

**59. **J. S. Peng and G. X. Li, *Introduction to modern quantum optics* (World Scientific, 1998).

**60. **N. B. Narozhny, J. J. Sanchez-Mondragon, and J. H. Eberly, “Coherence versus incoherence: collapse and revival in a simple quantum model,” Phys. Rev. A **23**, 236 (1981). [CrossRef]

**61. **A. Auffèves-Garnier, C. Simon, J. M. Gérard, and J. P. Poizat, “Giant optical nonlinearity induced by a single two-level system interacting with a cavity in the Purcell regime,” Phys. Rev. A **75**, 053823 (2007). [CrossRef]

**62. **J. Restrepo, C. Ciuti, and I. Favero, “Single-polariton optomechanics,” Phys. Rev. Lett. **112**, 013601 (2014). [CrossRef] [PubMed]

**63. **X. W. Xu, H. Wang, J. Zhang, and Y. X. Liu, “Engineering of nonclassical motional states in optomechanical systems,” Phys. Rev. A **88**, 063819 (2013). [CrossRef]

**64. **R. Jacob, S. Winnerl, M. Fehrenbacher, J. Bhattacharyya, H. Schneider, M. T. Wenzel, H. G. von Ribbeck, L. M. Eng, P. Atkinson, O. G. Schmidt, and M. Helm, “Intersublevel spectroscopy on single InAs-quantum dots by terahertz near-field microscopy,” Nano Lett. **12**, 4336–4340 (2012). [CrossRef] [PubMed]

**65. **R. Leijssen and E. Verhagen, “Strong optomechanical interactions in a sliced photonic crystal nanobeam,” Sci. Rep. **5**, 15974 (2015). [CrossRef] [PubMed]

**66. **F. Lecocq, J. D. Teufel, J. Aumentado, and R. W. Simmonds, “Resolving the vacuum fluctuations of an optomechanical system using an artificial atom,” Nature Phys. **11**, 635–639 (2015). [CrossRef]

**67. **G. S. Agarwal, L. M. Narducci, D. H. Feng, and R. Gilmore, “Intensity correlations of a cooperative system,” Phys. Rev. Lett. **42**, 1260 (1979). [CrossRef]

**68. **M. Florescu and S. John, “Single-atom switching in photonic crystals,” Phys. Rev. A **64**, 033801 (2001). [CrossRef]

**69. **C. W. Gardiner and P. Zoller, *Quantum noise* (Springer, 2004).

**70. **A. C. Pflanzer, O. Romero-Isart, and J. I. Cirac, “Optomechanics assisted by a qubit: from dissipative state preparation to many-partite systems,” Phys. Rev. A **88**, 033804 (2014). [CrossRef]

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

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

**73. **J. M. Pirkkalainen, S. U. Cho, F. Massel, J. Tuorila, T. T. Heikkilä, P. J. Hakonen, and M. A. Sillanpää, “Cavity optomechanics mediated by a quantum two-level system,” Nat. Commun. **6**, 6981 (2015). [CrossRef] [PubMed]

**74. **G. Via, G. Kirchmair, and O. Romero-Isart, “Strong single-photon coupling in superconducting quantum magnetomechanics,” Phys. Rev. Lett. **114**, 143602 (2015). [CrossRef] [PubMed]

**75. **P. D. Nation, J. Suh, and M. P. Blencowe, “Ultrastrong optomechanics incorporating the dynamical Casimir effect,” Phys. Rev. A **93**, 022510 (2016). [CrossRef]

**76. **A. Rundquist, M. Bajcsy, A. Majumdar, T. Sarmiento, K. Fischer, K. G. Lagoudakis, S. Buckley, A. Y. Piggott, and J. Vučković, “Nonclassical higher-order photon correlations with a quantum dot strongly coupled to a photonic-crystal nanocavity,” Phys. Rev. A **90**, 023846 (2014). [CrossRef]

**77. **J. P. Reithmaier, G. Sȩk, A. Löffler, C. Hofmann, S. Kuhn, S. Reitzenstein, L. V. Keldysh, V. D. Kulakovskii, T. L. Reinecke, and A. Forchel, “Strong coupling in a single quantum dot-semiconductor microcavity system,” Nature **432**, 197–200 (2004). [CrossRef] [PubMed]

**78. **M. A. Lemonde, N. Didier, and A. A. Clerk, “Enhanced nonlinear interactions in quantum optomechanics via mechanical amplification,” Nat. Commun. **7**, 11338 (2016). [CrossRef] [PubMed]

**79. **J. B. Clark, F. Lecocq, R. W. Simmonds, J Aumentado, and J. D. Teufel, “Observation of Strong Radiation Pressure Forces from Squeezed Light on a Mechanical Oscillator,” http://arxiv.org/abs/1601.02689 (2016).

**80. **R. J. Glauber, “The quantum theory of optical coherence,” Phys. Rev. **130**, 2529 (1963). [CrossRef]

**81. **B. M. Ann, Y. Song, J. Kim, D. Yang, and K. An, “Correction for the detector-dead-time effect on the second-order correlation of stationary sub-Poissonian light in a two-detector configuration,” Phys. Rev. A **92**, 023830 (2015). [CrossRef]

**82. **Y. H. Zhou, H. Z. Shen, and X. X. Yi, “Unconventional photon blockade with second-order nonlinearity,” Phys. Rev. A **92**, 023838 (2015). [CrossRef]

**83. **E. Wigner, “On the quantum correction for thermodynamic equilibrium,” Phys. Rev. **40**, 749 (1932). [CrossRef]

**84. **A. Kenfack and K. Zyczkowski, “Negativity of the Wigner function as an indicator of non-classicality,” J. Opt. B: Quantum Semiclass. Opt. **6**, 396 (2004). [CrossRef]

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

**86. **S. B. Li, X. B. Zou, and G. C. Guo, “Nonclassicality of quantum excitation of classical coherent field in photon-loss channel,” Phys. Rev. A **75**, 045801 (2007). [CrossRef]

**87. **A. Biswas and G. S. Agarwal, “Nonclassicality and decoherence of photon-subtracted squeezed states,” Phys. Rev. A **75**, 032104 (2007). [CrossRef]

**88. **P. T. Cochrane, T. C. Ralph, and G. J. Milburn, “Teleportation improvement by conditional measurements on the two-mode squeezed vacuum,” Phys. Rev. A **65**, 062306 (2002). [CrossRef]

**89. **S. Olivares and M. G. A. Paris, “De-Gaussification by inconclusive photon subtraction,” Laser Phys. **16**, 1533–1550 (2006). [CrossRef]

**90. **N. J. Cerf, O. Kruger, P. Navez, R. F. Werner, and M. M. Wolf, “Non-Gaussian cloning of quantum coherent states is optimal,” Phys. Rev. Lett. **95**, 070501 (2005). [CrossRef] [PubMed]

**91. **F. DellAnno, S. De Siena, L. Albano, and F. Illuminati, “Continuous-variable quantum teleportation with non-Gaussian resources,” Phys. Rev. A **76**, 022301 (2007). [CrossRef]