# Low-power phonon lasing through position-modulated Kerr-type nonlinearity

P. Djorwe, Yan Pennec, Bahram Djafari-Rouhani

► To cite this version:

P. Djorwe, Yan Pennec, Bahram Djafari-Rouhani. Low-power phonon lasing through position-modulated Kerr-type nonlinearity. Scientific Reports, 2019, 9 (1), pp.1684. 10.1038/s41598-019-38578-8. hal-03146998

HAL Id: hal-03146998

https://hal.science/hal-03146998

Submitted on 19 Feb 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
SCIENTIFIC REPORTS

OPEN

# Low-power phonon lasing through position-modulated Kerr-type nonlinearity

P. Dhorwe, Y. Pennec & B. Djafari-Rouhani

Received: 24 July 2018

Accepted: 14 December 2018

Published online: 08 February 2019

We demonstrate low-power amplification process in cavity optomechanics (COM). This operation is based on the nonlinear position-modulated self-Kerr interaction. Owing to this nonlinear term, the effective coupling highly scales with the photon number, resulting in a giant enhancement of the cooperativity. Even for small nonlinearity, the system reaches the amplification threshold for weak driving strength, leading to low-power phonon lasing. This amplifier can be phase-preserving and provides a practical advantage related to the power consumption issues. This work opens up new avenues to perform low-power and efficient amplifiers in optomechanics and related fields.

Cavity optomechanics (COM), which is devoted to explore interaction between electromagnetic radiation and mechanical object, provides a platform to perform phonon lasing action¹,² based on backaction amplification³. At the threshold of parametric instability, backaction-induced mechanical gain overcomes mechanical loss, resulting to an amplification process that leads to coherent phonon oscillations¹. Similar to stimulated emission of photon lasing in cavity with a gain medium, backaction amplification induces stimulated phonon lasing from the parametric instability's threshold⁴. Both phonon lasing and amplification are actively studied in optomechanics, and other fields, for technological purposes.

Single-phonon Fock state detection has been recently performed in⁵,⁶, and that paves a way to variety of quantum state engineering tasks including, quantum information processing⁷ and quantum entanglement of remote mechanical elements⁸. The development of single phonon source reveals also a process towards technologies for precision sensing⁹. Moreover, intensive researches are carried out in order to improve these achievements, for instance to lower the power need for phonon lasing. Such performance has been recently realized using *PT*-symmetry optomechanics¹⁰, and the polarization of light in coupled optomechanical devices¹¹.

To use quantum signals, those involving only a few quanta, they need to be amplified and must display some degree of purity to carry certain amount of information. This issue can be handled by characterizing the amplification process as the system gets closer to the phonon lasing threshold. These characteristics include phase-preserving amplification¹²⁻¹⁴, the power gain, the added noise, and the gain-bandwidth product¹⁵⁻¹⁸. The purpose of a phase-preserving amplifier is to make a weak signal large, regardless of its phase, so that it can be perceived by devices unable to resolve the original signal, while sacrificing as little as possible in signal-to-noise ratio¹⁴. Engineering amplifier having large power gain with quantum-limited added noise¹⁷,¹⁸, and without limitation on the gain-bandwidth product¹⁷ are useful for applications such as sensing³,⁹ and quantum information processing⁷. Another interesting feature of amplifiers is their ability to reach high gain amplification for low input power. Such an amplifier has been recently realized, and has shown high gain for low-power operation and low quantum limit noise performance¹⁹.

Recently, nonlinear position-modulated self-Kerr interaction has been engineered in COM, and it was found that it leads to an effective coupling that scales with the square of the photon number²⁰. Such nonlinearity can be derived in the situation where Kerr nonlinear coefficient is modulated by the mechanical position. Cavity polaritons have been recently proposed as good candidates for such interactions in optomechanical systems²¹. Position modulation of Kerr nonlinearity can be possible in optical cavity, where the entire space between the two mirrors (where one of them can move) of the cavity is filled by a $\chi^{(3)}$ medium. However, the resulting position-dependent Kerr term is weak, since it is inversely proportional to the cavity length²⁰. Therefore, the promising systems to implement this kind of interaction are superconducting microwave resonators, where a large optical Kerr nonlinearity at microwave frequencies can be generated²². Beside of microwave systems, atom-optical systems can

Institut d'Electronique, de Microélectronique et Nanotechnologie, UMR CNRS 8520 Université de Lille, Sciences et technologies, Villeneuve d'Ascq, 59652, France. Correspondence and requests for materials should be addressed to P.D. (email: philippe.djorwe@univ-lille1.fr)
be also used for such nonlinear couplings. Position-dependent Kerr nonlinearity has recently led to a strong coupling, with a magnitude exceeding the strength of the Kerr coefficient, even for low driving strength$^{20,21}$. In the red sideband regime, these authors have demonstrated motional cooling, mode splitting and multistability for low-power$^{20}$. Owing to these interesting nonlinear phenomena, our aim here is to perform low-power phonon laser at the blue sideband, resulting from low-power amplification induced by this nonlinearity. To understand this low-power operation, we have derived the cooperativity that highly scales with the photon number. This pushes the system fastly near to the amplification's threshold. We have characterized this amplifier, which shows high gain and phase-preserving close to the phonon lasing threshold. This work opens up promising ways to develop low-power amplifiers based on the position-modulated self-Kerr interaction in COMs and superconducting microwave setups$^{22,23}$.

## Results

**Hamiltonian and dynamical equations.** Position-modulated self-Kerr nonlinearity can be engineered in COMs$^{24–26}$ or in the superconducting microwave setups exhibiting a giant Kerr nonlinearity$^{22}$. The idea of this engineering is based on the fact that, the Kerr nonlinear coefficient is modulated by the position of the mechanical resonator connected to the system. Such nonlinear interaction, has been recently investigated in$^{20,21}$. This has led to an effective coupling that highly scales with the photon number. In the red sideband, this effective coupling leads to low powers required for motional cooling, the emergence of multistability, and other interesting nonlinear features. Owing to these exciting nonlinear effects, here we move to the blue sideband and investigate on phonon lasing and amplification phenomenon. In the frame rotating at the driving field frequency $\omega_p$, the Hamiltonian describing the generic system (with $\hbar = 1$) is given by,

$$H = -\Delta_0 a^\dagger a + \omega_m b^\dagger b - g_l a^\dagger a(b^\dagger + b) - g_{nl} a^\dagger a a(b^\dagger + b) + E(a + a^\dagger). \quad (1)$$

In this Hamiltonian, $a(b)$ is the annihilation bosonic operator for the intracavity field (mechanical resonator), $H_{l,int} = -g_l a^\dagger a(b^\dagger + b)$ and $H_{nl,int} = -g_{nl} a^\dagger a^\dagger aa(b^\dagger + b)$ describe the linear and nonlinear interactions. The first two terms represent the cavity and mechanical free energy respectively, while the last term stands for the driving energy. The parameters $\omega_m$ and $\Delta_0 = \omega_p - \omega_{cav}$ are the mechanical frequency of the resonator and the detuning between the driving ($\omega_p$) and the cavity eigenfrequency ($\omega_{cav}$). The linear and nonlinear optomechanical couplings are denoted by $g_l$ and $g_{nl}$, respectively. The mechanical displacement is defined as $x = x_{ZPF}(b + b^\dagger)$, where $x_{ZPF} = \sqrt{\hbar/(2m\omega_m)}$ is the zero-point fluctuation amplitude of the mechanical resonator, with $m$ its effective mass. From the Hamiltonian given in Eq. (1), the Nonlinear Langevin Equations (NLEs), including cavity ($\kappa$) and mechanical ($\gamma_m$) dissipations, can be derived as,

$$\begin{cases} \dot{a} = \left[ i(\Delta_0 + g_l(b^\dagger + b) + 2g_{nl}(b^\dagger + b)a^\dagger a) - \frac{\kappa}{2} \right] a - iE, \\ \dot{b} = -\left( i\omega_m + \frac{\gamma_m}{2} \right) b + ig_l a^\dagger a + ig_{nl} a^\dagger a a. \end{cases} \quad (2)$$

Throughout the work, we assume the hierarchy of parameters $\gamma_m$, $g_l \ll \kappa \ll \omega_m$, similar to the experiments carried out in the resolved sideband regime$^{5,6}$. Our numerical and analytical investigations will be done at the sideband $\Delta_0 = \omega_m$. To get insight of the phonon lasing phenomenon, we linearize Eq. (2) in the limit of large driving field. In this case, both intracavity field (a) and mechanical degrees of freedom (b) can be split into their average fields ($\alpha(t), \beta(t)$) with some amount of fluctuations ($\delta\alpha(t), \delta\beta(t)$) as follows,

$$\begin{cases} \delta\alpha(t) = a(t) - \alpha(t), \\ \delta\beta(t) = b(t) - \beta(t). \end{cases} \quad (3)$$

Using Eq. (3) in Eq. (2), leads to the steady-state dynamics,

$$\begin{cases} \dot{\alpha} = \left( i\Delta - \frac{\kappa}{2} \right) \alpha + \sqrt{\kappa} \delta\alpha^{in}, \\ \dot{\beta} = - \left( i\omega_m + \frac{\gamma_m}{2} \right) \beta + ig_l |\alpha|^2 + ig_{nl} |\alpha|^4, \end{cases} \quad (4)$$

with the corresponding lowest order fluctuations dynamics, including noises,

$$\begin{cases} \delta\dot{\alpha} = \left(i\Delta - \frac{\kappa}{2}\right)\delta\alpha + i\chi(\delta\beta^* + \delta\beta) + i\eta(\delta\alpha^* + \delta\alpha) + \sqrt{\kappa}\delta\alpha^{in}, \\ \delta\dot{\beta} = -\left(i\omega_m + \frac{\gamma_m}{2}\right)\delta\beta + i\chi(\delta\alpha^* + \delta\alpha) + \sqrt{\gamma_m}\delta\beta^{in}. \end{cases} \quad (5)$$

In Eq. (4) and Eq. (5), we have set for convenience $-iE = \sqrt{\kappa} \alpha^{in}$ where the driving strength $\alpha^{in}$ is related to the input power $P_m$ as $\alpha^{in} = \sqrt{\frac{P_m}{\hbar\omega_p}}$. We have defined the effective coupling as $\chi = \chi_0 \alpha$ with $\chi_0 = (g_l + 2g_{nl}|\alpha|^2)$, $\eta = 4g_{nl}|\alpha|^2 \text{Re}(\beta)$ and $\Delta = \Delta_0 + 2\chi_0 \text{Re}(\beta)$. The noise operators are characterized by the noncommuting properties, i.e., $\langle \delta\alpha^{in}(t) \rangle = 0$, $\langle \delta\alpha^{in\dagger}(t')\delta\alpha^{in}(t) \rangle = n_\alpha \delta(t' - t)$, and $\langle \delta\alpha^{in}(t')\delta\alpha^{in\dagger}(t) \rangle = (n_\alpha + 1)\delta(t' - t)$ for the input
field and $\langle\delta\beta^{in}(t)\rangle = 0$, $\langle\delta\beta^{in\dagger}(t')\delta\beta^{in}(t)\rangle = n_{th}\delta(t' - t)$, and $\langle\delta\beta^{in}(t')\delta\beta^{in\dagger}(t)\rangle = (n_{th} + 1)\delta(t' - t)$ for the thermal bath. The quantities $n_\alpha$ and $n_{th}$ are the equilibrium occupation numbers for the input field and mechanical oscillator, respectively.

**Stability.** The steady-state solutions $\alpha_s$ and $\beta_s$ are derived from Eq. (4), providing $\dot{\alpha} = 0$ and $\dot{\beta} = 0$. These solutions are given by,

$$ \left\{ \begin{aligned} \alpha_s &= \frac{\sqrt{\kappa} \alpha^{in}}{\left(\frac{\kappa}{2} - i[\Delta_0 + 2\chi_0 \text{Re}(\beta_s)]\right)}, \\ \beta_s &= (g_l + g_{nl} |\alpha_s|^2) \frac{|\alpha_s|^2}{\omega_m - i \frac{\gamma_m}{2}}. \end{aligned} \right. \qquad (6) $$

By setting the intracavity intensity as $I = |\alpha_s|^2$, one shows from Eq. (6) that it is solution of the following seventh-order polynomial equation,

$$ I^7 + a_6 I^6 + a_5 I^5 + a_4 I^4 + a_3 I^3 + a_2 I^2 + a_1 I + a_0 = 0, \qquad (7) $$

with $a_6 = \frac{3g_l}{g_{nl}}$; $a_5 = \frac{13g_l^2}{4g_{nl}^2}$; $a_4 = \frac{3g_l^3 + \Delta\omega_m g_{nl}}{2g_{nl}^3}$; $a_3 = \frac{g_l(g_l^3 + 3\Delta\omega_m g_{nl})}{4g_{nl}^4}$; $a_2 = \frac{\Delta\omega_m g_l^2}{4g_{nl}^4}$; $a_1 = \frac{\left(\frac{\kappa^2}{4} + \Delta^2\right)\omega_m^2}{(2g_{nl})^4}$; $a_0 = -\frac{\kappa(\alpha^{in}\omega_m)^2}{(2g_{nl})^4}$. These steady state solutions are physically meaningless, unless they are stable. The stability can be given explicitly through Routh-Hurwitz criterion$^{26}$. However, we analyze it here through linear stability theory, and confirm it with parametric instability threshold, since we are in the blue sideband regime. To this end, we start by writing Eq. (5) in the following compact form,

$$ \delta\dot{X} = M\delta X + \varepsilon, \qquad (8) $$

with $\delta X = (\delta\beta, \delta\beta^*, \delta\alpha, \delta\alpha^*)^T$ and $\varepsilon = (\sqrt{\gamma_m}\delta\beta^{in}, \sqrt{\gamma_m}\delta\beta^{in*}, \sqrt{\kappa}\delta\alpha^{in}, \sqrt{\kappa}\delta\alpha^{in*})^T$. The matrix M is given by,

$$ M = \begin{bmatrix} -(i\omega_m + \frac{\gamma_m}{2}) & 0 & i\chi & i\chi \\ 0 & \left(i\omega_m - \frac{\gamma_m}{2}\right) & -i\chi & -i\chi \\ i\chi & i\chi & \left(i\tilde{\Delta} - \frac{\kappa}{2}\right) & i\eta \\ -i\chi & -i\chi & -i\eta & -(i\tilde{\Delta} + \frac{\kappa}{2}) \end{bmatrix}, \qquad (9) $$

with $\tilde{\Delta} = \Delta + \eta$, and where we have chosen the phase reference of the cavity field so that $\alpha_s$ is real.

The system is stable if all the real part of eigenvalues $(\lambda_{i=1..4})$ of the matrix M are negative (Re($\lambda_{i=1..4}$) < 0). This stability depends on steady-state solutions $\alpha_s$ and $\beta_s$, and is shown in Fig. 1a. The blue color in Fig. 1a depicts stable parameters space, while the red area shows the unstable zone. As the self-Kerr term ($g_{nl}$) is increasing, the system becomes unstable and the stability is limited for relatively weak driving strength $\alpha^{in}$. The green dashed line in Fig. 1a holds for the condition $C = 1$, where $C = 4\chi^2/(\gamma_m\kappa) = \gamma_{opt}/\gamma_m$ generalizes the cooperativity at the limit $g_{nl} \approx 0$. This cooperativity, with $\gamma_{opt}$ being the optical damping, depicts the border between stable (linear) and unstable (nonlinear) regimes. Moreover, $C = 1$ defines the threshold of the phonon lasing in the optomechanical blue sideband. It results that the self-Kerr nonlinearity induces low power phonon lasing action, that can be understood by the enhancement of the effective coupling ($\chi$) shown in Fig. 1b. This coupling enhancement is a direct consequence of the fact that $\chi$ highly scales with the phonon number in the presence of $g_{nl}$, leading to a large cooperativity (see inset of Fig. 1b). Hence, as $g_{nl}$ increases, the lasing threshold is shifted towards weak driving strength $\alpha^{in}$, revealing low-power phonon lasing in our proposal.

**Low power phonon lasing.** To evaluate the stimulated emission phonon number, we have simulated the classical equivalent of the nonlinear equation given in Eq. (2). This is valid for large enough photon number in the system ($|\alpha|^2 \gg 1$), so that both fluctuations from intracavity field and mechanical resonator can be neglected. By Fast Fourier Transforming (FFT) the results, and collecting the mechanical peaks at the resonance, we can obtain the phonon number versus the driving field $\alpha^{in}$. Analytical results can be also obtained through an approximated analysis that is detailed in section Methods. Figure 1c displays a log-log scale representation, and it shows that the effect of the nonlinear term $g_{nl}$ is negligible on the phonon number in the linear regime ($\alpha^{in} \lesssim 10^2\sqrt{\omega_m}$). However, this nonlinearity strongly enhances the increase of the phonon number near the instability thresholds. The linear scale representation of the phonon number is displayed in the inset of Fig. 1c, where a comparison is made with the full nonlinear results. Indeed, full and dotted lines show the numerical and analytical results, respectively. This yields good agreement between numerical calculation and analytics, that is detailed in section Methods. The input fields at the lasing threshold ($\alpha_{th}^{in}$) are depicted by the green dots and can be obtained through the condition $\gamma = \gamma_{opt}$ (or $C = 1$). As $g_{nl}$ is increasing, Fig. 1c reveals a low driving strength $\alpha^{in}$ required for phonon lasing threshold. In Fig. 1d, we have represented the phonon number versus $g_{nl}$ for different $\alpha^{in}$. It results that $g_{nl}$ enhances stimulated emission of phonons. Indeed, for $g_{nl} \sim 0$, there is no lasing up to $\alpha^{in} = 3 \times 10^2\sqrt{\omega_m}$ in
<figure><img src="image_2.png" /></figure>

<figure><img src="image_3.png" /><figcaption>**Figure 1.** (a) Stability diagram. Blue (red) color is stable (unstable). The green dashed curve shows the border between stable and unstable regions, and corresponds to the lasing threshold $C = 1$, with the cooperativity $C = 4\chi^2/(\gamma_m\kappa)$. (b) Effective coupling $\chi$ versus $\alpha^{in}$ for different values of $g_{nl}$. Full (dashed) curves are stable (unstable). The inset shows the cooperativity $C$ versus $\alpha^{in}$, for the corresponding values of $g_{nl}$. (c) Phonon number versus $\alpha^{in}$ in log-log scale, showing the fast increase near the thresholds. The inset displays the same figure in the linear scale, where full curves are from numerical simulation of classical version of Eq. (2), while the dotted lines represent the analytical approximation from $N \sim \exp(-\gamma_{eff}t)$. (d) Analytical (approximated) phonon number versus $g_{nl}$. In (c,d), the different lasing thresholds ($C = 1$) are indicated by the green dots, namely $\alpha_{th}^{in}$ and $g_{nl}^{th}$ for the input field and the nonlinear coefficient, respectively. The used parameters are$^{20,23}$, $\gamma_m = 10^{-5}\omega_m$, $\kappa = 2 \times 10^{-2}\omega_m$, $g_l = 2 \times 10^{-5}\omega_m$, and $\Delta_0 = \omega_m$.</figcaption></figure>

Fig. 1d. Therefore, by adding a small amount of nonlinearity to the system, the lasing threshold rises up and is shifted towards small $g_{nl}$ as $\alpha^{in}$ increases. Briefly speaking, the higher is the driving strength, smaller is the amount of nonlinearity ($g_{nl}$) to reach the lasing threshold, and vice-versa. The nonlinear coefficients required to reach the lasing threshold are indicated by the green dots, namely $g_{nl}^{th}$. This can be understood from the dynamics of $b(t)$ in Eq. (2), showing that $g_{nl}$ supplies energy ($\propto g_{nl}|\alpha|^4$) to drive the mechanical resonator. This reveals why self-Kerr nonlinearity studied here, is quite interesting and different from Kerr nonlinearity$^{22,24,26}$, quadratic nonlinearity$^{25,27}$ and Duffing nonlinearities$^{28}$ also studied in optomechanics. Indeed, these nonlinearities either shift the cavity frequency or the mechanical frequency, but none of them is directly driving the mechanical resonator. Free-carrier and thermal nonlinearities can compete with self-Kerr nonlinearity, but their effects are negligible in the weak driving limit considered in our work.

**Low power amplification.** As self-Kerr nonlinearity enhances low-power phonon lasing (Fig. 1c), this also reveals low-power amplification process in the system. The feature of the nondegenerate parametric amplifier here is to convert a pump mode photon into, one photon signal mode and one idler phonon mode. This can lead to the fact that weak incident signals are amplified, with a minimum possible added noise. Such amplification process can be seen from the linearization of the interaction Hamiltonian of our system. Indeed, the interaction of our system is captured by the Hamiltonian,

$$H_{int} = H_{l,int} + H_{nl,int} = -g_{l}a^{\dagger}a(b^{\dagger} + b) - g_{nl}a^{\dagger}aa(b^{\dagger} + b), \qquad (10)$$

which leads to its linearized form,

$$H_{int}^{lin} = -\chi(\delta\alpha^{\dagger}\delta\beta^{\dagger} + \delta\alpha\delta\beta) - \chi(\delta\alpha^{\dagger}\delta\beta + \delta\alpha\delta\beta^{\dagger}), \qquad (11)$$

after having used Eq. (3). Furthermore, Eq. (11) has been obtained by omitting static terms since they are taken into account in the frequency shift $\tilde{\Delta}$, and the higher order fluctuations terms have been neglected as being smaller than $\alpha_s$. The second term on the right-hand side of Eq. (11) stands for the counter rotating terms $H_{CR} = -\chi(\delta\alpha^{\dagger}\delta\beta + \delta\alpha\delta\beta^{\dagger})$ and can be neglected in the rotating wave approximation (RWA)$^1$. However, the first term ($H_R = -\chi(\delta\alpha^{\dagger}\delta\beta^{\dagger} + \delta\alpha\delta\beta)$) describes a nondegenerate parametric amplifier, where a pump mode (pho-
<figure><img src="image_4.png" /><figcaption>Figure 2. (a) Output noise power spectral density (PSD) for different values of g<sub>nl</sub>. The magnitude of S<sub>out</sub> = √⟨δα<sup>out†</sup>δα<sup>out</sup>⟩ reveals the amplitude of the mechanical resonator, detected at the output for a given frequency. (b) Values of the output spectrum (in (a)) at the resonance (ω = −Δ). (c) Amplifier’s gain (see Eq. (16)) versus α<sup>in</sup>, for different values of g<sub>nl</sub>. (d) Cavity frequency shift Δ<sub>nl</sub> = Δ̃ − Δ<sub>0</sub> versus α<sup>in</sup>. The driving strength in (a,b) is α<sup>in</sup> = 2 × 10<sup>2</sup>√ω<sub>m</sub>, with n<sub>th</sub> = n<sub>α</sub> = 1 for (a). The other parameters are the same as in Fig. 1.</figcaption></figure>

ton) is converted into two quanta, one in the signal mode (photon), and the other in the idler (phonon). To characterize this amplifier, we neglect non resonant terms in Eq. (5) and rewrite it in the RWA as,

$$
\begin{equation}
\begin{cases}
\delta\alpha = \left(i\Delta - \frac{\kappa}{2}\right)\delta\alpha + i\chi\delta\beta^{\dagger} + \sqrt{\kappa}\delta\alpha^{in}, \\
\delta\beta^{\dagger} = \left(i\omega_m - \frac{\gamma_m}{2}\right)\delta\beta^{\dagger} - i\chi\delta\alpha + \sqrt{\gamma_m}\delta\beta^{in}.
\end{cases}
\tag{12}
\end{equation}
$$

By solving Eq. (12) in the Fourier space, together with the input-output relation, δα<sup>out</sup> = δα<sup>in</sup> − √κδα<sup>15,18</sup>, we can evaluate the output field δα<sup>out</sup>. This output field is a key element to characterize the gain and added noise of the amplifier. More details on calculations leading to the gain and added noise are reported in section Methods, while specific results are shown in what follows. The input-output relation leads to the output field,

$$
\delta\alpha^{\text{out}} = \left(1 - \sqrt{\kappa} \chi_{\text{eff}}^{c}\right) \delta\alpha^{\text{in}} - i\sqrt{\kappa} \eta_{c} \delta\beta^{\text{in}\dagger}, \quad (13)
$$

where

$$
\eta_c = \frac{\chi_m \chi_c \chi \sqrt{\gamma_m}}{1 - \chi_m \chi_c \chi^2}; \quad \chi_{\text{eff}}^c = \frac{\chi_c \sqrt{\kappa}}{1 - \chi_m \chi_c \chi^2}, \qquad (14)
$$

with the susceptibilities χ_c = [κ/2 − i(ω + Δ)]<sup>−1</sup>, and χ_m = [γ_m/2 − i(ω + ω_m)]<sup>−1</sup>. In Eq. (13), the coefficient in front of the incident signal δα<sup>in</sup> characterizes the amplification gain, while the one in front of the thermal noise informs on the added noise. These characteristics are deduced from the output noise power spectral density (PSD) defined as<sup>17,18</sup>,

$$
S_{\text{out}} = \frac{1}{2} (\langle \delta\alpha^{\text{out}\dagger} \delta\alpha^{\text{out}} + \delta\alpha^{\text{out}} \delta\alpha^{\text{out}\dagger} \rangle). \quad (15)
$$

The output PSD (S<sub>out</sub>) is shown in Fig. 2a, and reveals an amplification process at the resonance ω = −Δ,
induced by the nonlinear term g<sub>nl</sub>. For α<sup>in</sup> = 2 × 10<sup>2</sup>√ω<sub>m</sub> and g<sub>nl</sub> ≤ 1 × 10<sup>−3</sup>g<sub>p</sub>, there is no amplification as
shown in Fig. 1c. However, for g<sub>nl</sub> ≈ 1.85 × 10<sup>−3</sup>g<sub>p</sub>, Fig. 2a clearly shows amplification, meaning that g<sub>nl</sub> brings
the system near the lasing threshold even for a weak driving strength. In order to appreciate this amplification, the
output PSD has been plotted at the resonance $\omega = -\Delta$ (for $\alpha^{in} = 2 \times 10^2 \sqrt{\omega_m}$) versus $g_{nl}$, and shown in Fig. 2b. The red and black curves correspond to ($n_\alpha = 1, n_{th} = 1$) and ($n_\alpha = 1, n_{th} = 0$), respectively. It results that thermal noise has an impact on the amplification, revealing the effect of the added noise on the amplifier. This can be pointed out by evaluating both the power gain $G(\omega)$ and the added noise $N(\omega)$. The amplification of the signal $\delta\alpha^{in}$ is measured through the gain $|1 - \sqrt{\kappa} \chi_{eff}^c|^2$, which can be simplified (see section Methods) at the resonance ($\omega = -\Delta$) as,

$$G = \left| \frac{C + 1}{1 - C} \right|^2. \qquad (16)$$

The noise performance of the amplifier is figured out through the input-referred added noise, defined as $N(\omega) = (S_{out} - S_{in})/G$, where $S_{in} = \frac{1}{2}$ is the vacuum input noise driving the cavity. As shown in section Methods, this expression reduces at the resonance to,

$$N = \frac{4C(n_{th} + \frac{1}{2})}{|C + 1|^2} + n_{\alpha}. \qquad (17)$$

It can be seen from Eq. (16) that the gain is greatly enhanced as the system approaches the lasing threshold at the cooperativity $C \to 1$, revealing the amplification process. As $C$ strongly scales with the self-Kerr term $g_{nl}$, it results an enhancement of amplification induced by $g_{nl}$ as shown in Fig. 2c. Furthermore, we have $N \to (n_{eff} + \frac{1}{2})$ for $C \to 1$, where $n_{eff} = n_{th} + n_\alpha$ is the effective phonon number of the mechanical resonator. This reveals that the amplifier reaches the quantum limit for phase-preserving$^{13,14}$ near the lasing threshold for $n_{eff} = 0$, and if there are no additional loss channels. However, for non zero effective (thermal) phonon number, $N$ linearly increases with $n_{eff}$ degrading the signal amplification purity as depicted by the red curve in Fig. 2b. This is not the case in the amplifier studied in$^{13,17}$ where the gain can take arbitrarily large values, and where any thermal noise contribution is suppressed for a large-gain limit$^{17}$. Despite of this, the practical advantage of our amplifier is its low power consumption due to the presence of $g_{nl}$.

## Discussion

We have carried out investigation on position-modulated Kerr-type nonlinearity in blue sideband optomechanics. This is mainly focused on phonon lasing and amplification process. The amplification, which is based on the recorded output field, is shown to be as a manifestation of phonon lasing inside the cavity. Both phenomena are enhanced by this nonlinear term, through the condition $C \to 1$. Indeed, a cooperativity equal to unity fulfills phonon lasing threshold requirement, large gain amplification, and condition of quantum limit for a phase-preserving amplifier. We have shown that these features happen for low-power strength, which is the figure of merit of the used nonlinear term, compared to those known so far$^{22,24-28}$. This is shown in Fig. 1a, where for $g_{nl} \neq 0$, the system reaches nonlinear regime for weak driving. Owing to such feature, we were able to drop counter rotating terms in our analytical treatment, since they start to manifest in strong driving limit. The green dashed line in Fig. 1a shows the threshold of phonon lasing using the rotating wave approximation, while the border between stable (blue) and unstable (red) regimes shows the same threshold when counter rotating terms are included. It results a good agreement between numerical and analytical calculations, confirming the validity of the rotating wave approximation in our investigation. Moreover, this is highlighted in Fig. 1c where a good matching between numerical (full curves) and analytical (dashed curves) results is induced by $g_{nl}$. For acoustic excitations, using weak driving strength to reach phonon lasing and large gain for the amplified phonon source having minimum added noise could be an interesting achievement. This can be useful for phonon information processing, and the position-modulated Kerr-type nonlinearity can be helpful. At the resonance and for $n_{eff} \to 0$, the amplifier reaches the quantum limit $N \to \frac{1}{2}$ at the threshold $C \to 1$ no matter is the strength of $g_{nl}$. For $n_{eff} \neq 0$ however, the performance of the amplifier is impaired by both thermal and input field fluctuations as shown in Fig. 2b. Close to the quantum limit and for lower sideband regime, this amplifier performs better when the modified Kerr term is involved. This is because the effective quantum population $n_{eff}$ is more reduced for $g_{nl} \neq 0$ as recently shown in$^{20}$. Moreover, a practical benefit of the amplifier discussed here is its low-power consumption to reach high amplification gain (phonon lasing).

## Methods

Figure 2d shows that the cavity frequency shift ($\tilde{\Delta} - \Delta_0$) is weak in the linear regime that the intracavity field remains at the mechanical sideband $\tilde{\Delta} \sim \omega_m$. This allows us to get analytical expression of stimulated phonon number by introducing the slowly varying amplitude variables$^{29}$, $\delta\tilde{\alpha} = \delta\alpha\exp(\frac{\kappa}{2} - i\tilde{\Delta})t$ and $\delta\tilde{\beta} = \delta\beta\exp(i\omega_m)t$. Using these variables in the cavity field equation in Eq. (12) yields,

$$\delta\tilde{\alpha} = \int_{-\infty}^{t} i\chi\delta\tilde{\beta}^{\dagger}e^{\frac{\kappa}{2}\tau}d\tau. \qquad (18)$$

As we are dealing with weak coupling regime in our analysis ($\kappa \gg \chi$, see Fig. 1b), we can adiabatically eliminate $\delta\alpha(t)$, then $\chi$ can be taken out of the integral Eq. (18). Moreover, $\kappa \gg \gamma$ indicates that the evolution of $\delta\tilde{\beta}(t)$ is much slower that $\delta\alpha(t)$, meaning that $\delta\tilde{\beta}(t)$ can be considered as a constant term in Eq. (18). Under these conditions, the integration of Eq. (18) yields,
$$ \delta\alpha = \frac{2i\chi}{\kappa}\delta\beta^\dagger. \qquad (19) $$

By replacing Eq. (19) in Eq. (12), one gets

$$ \delta\beta^\dagger = \left(i\omega_m - \frac{\gamma_{\text{eff}}}{2}\right)\delta\beta^\dagger, \qquad (20) $$

where the effective damping is $\gamma_{\text{eff}} = \gamma_m - \gamma_{\text{opt}}$, with the optical damping $\gamma_{\text{opt}} = \frac{4\chi^2}{\kappa}$. The solution of $C = 1$ (or $\gamma_{\text{eff}} = 0$), with the cooperativity $C = 4\chi^2/(\gamma_m\kappa)$, gives the lasing threshold shown by the green dashed curve in Fig. 1a. The stimulated phonon number depicted in Fig. 1(c,d), is deduced from Eq. (20), by evaluating $N = \langle\delta\beta^\dagger\delta\beta\rangle = N_0e^{-\gamma_{\text{eff}}t}$, with $N_0 = 1$ the phonon number at $t = 0^{\text{th}}$.

The amplification shown in Fig. 2(a,b) are obtained by solving Eq. (12) in the frequency domain, together with the input-output relation$^{15,18}$. In the Fourier space, Eq. (12) leads to,

$$ \left\{ \begin{aligned} -i\omega\delta\alpha &= \left(i\tilde{\Delta} - \frac{\kappa}{2}\right)\delta\alpha + i\chi\delta\beta^\dagger + \sqrt{\kappa}\delta\alpha^{\text{in}}, \\ -i\omega\delta\beta^\dagger &= \left(i\omega_m - \frac{\gamma_m}{2}\right)\delta\beta^\dagger - i\chi\delta\alpha + \sqrt{\gamma_m}\delta\beta^{\text{in}\dagger}. \end{aligned} \right. \qquad (21) $$

After some calculations, one obtains

$$ \left\{ \begin{aligned} \delta\alpha(\omega) &= \chi_{\text{eff}}^c \delta\alpha^{\text{in}} + i\eta_c \delta\beta^{\text{in}\dagger}, \\ \delta\beta^\dagger(\omega) &= \eta_m \delta\beta^{\text{in}\dagger} - i\chi_{\text{eff}}^m \delta\alpha^{\text{in}}, \end{aligned} \right. \qquad (22) $$

where

$$ \left\{ \begin{aligned} \eta_c &= \frac{\chi_m \chi_c \chi \sqrt{\gamma_m}}{1 - \chi_m \chi_c \chi^2}, & \chi_{\text{eff}}^c &= \frac{\chi_c \sqrt{\kappa}}{1 - \chi_m \chi_c \chi^2}, \\ \eta_m &= \frac{\chi_m \sqrt{\gamma_m}}{1 - \chi_m \chi_c \chi^2}, & \chi_{\text{eff}}^m &= \frac{\chi_m \chi_c \chi \sqrt{\kappa}}{1 - \chi_m \chi_c \chi^2}, \end{aligned} \right. \qquad (23) $$

with the susceptibilities $\chi_c = \left[\frac{\kappa}{2} - i(\omega + \tilde{\Delta})\right]^{-1}$, and $\chi_m = \left[\frac{\gamma_m}{2} - i(\omega + \omega_m)\right]^{-1}$. Using the input-output relation, one gets the output field,

$$ \delta\alpha^{\text{out}} = (1 - \sqrt{\kappa} \chi_{\text{eff}}^c) \delta\alpha^{\text{in}} - i\sqrt{\kappa} \eta_c \delta\beta^{\text{in}\dagger}, \qquad (24) $$

which leads to the output PSD,

$$ \begin{aligned} S_{\text{out}} &= \frac{1}{2}((\delta\alpha^{\text{out}\dagger}\delta\alpha^{\text{out}} + \delta\alpha^{\text{out}}\delta\alpha^{\text{out}\dagger})) \\ &= \kappa|\eta_c|^2(n_{th} + \frac{1}{2}) + |1 - \sqrt{\kappa}\chi_{\text{eff}}^c|^2(n_{\alpha} + \frac{1}{2}) \\ &= \kappa|\eta_c|^2(n_{th} + \frac{1}{2}) + \mathcal{G}\left(n_{\alpha} + \frac{1}{2}\right). \end{aligned} \qquad (25) $$

The amplifier is then characterized by the power gain,

$$ \mathcal{G}(\omega) = \left|1 - \sqrt{\kappa} \chi_{\text{eff}}^c\right|^2 = \left| \frac{1 - C_{\omega} - 2\left(1 - \frac{2i}{\gamma_m}(\omega + \omega_m)\right)}{1 - C_{\omega}} \right|^2, \qquad (26) $$

with $C_{\omega} = \frac{4}{\gamma_m \kappa}(\chi^2 + (\omega + \omega_m)(\omega + \tilde{\Delta})) - \frac{2i}{\gamma_m \kappa}(\gamma_m(\omega + \tilde{\Delta}) + \kappa(\omega + \omega_m))$. The input-referred added noise quanta to the amplifier is defined as,

$$ \begin{aligned} \mathcal{N}(\omega) &= (S_{\text{out}} - S_{\text{in}})/\mathcal{G} \\ &= \frac{\kappa|\eta_c|^2}{\mathcal{G}}\left(n_{\text{th}} + \frac{1}{2}\right) + n_{\alpha} \\ &= \frac{4C_{\omega}\left(n_{\text{th}} + \frac{1}{2}\right)}{\left|1 - C_{\omega} - 2\left(1 - \frac{2i}{\gamma_m}(\omega + \omega_m)\right)\right|^2 + n_{\alpha}}, \end{aligned} \qquad (27) $$
where $S_{out}$ is given in Eq. (15), and $S_{in} = \frac{1}{2}$ holds for the vacuum input noise driving the cavity. At the resonance ($\omega = -\Delta$), these amplifier's characteristic reduce to,

$$G = \left| \frac{C + 1}{1 - C} \right|^2, \qquad (28)$$

and

$$N = \frac{4C(n_{th} + \frac{1}{2})}{|C + 1|^2} + n_{\alpha}, \qquad (29)$$

which clearly show: (i) amplification near the lasing threshold ($G \to \infty$ for $C \to 1$), and the fact that (ii) the amplifier reaches the quantum limit for a phase-preserving amplifier ($N \to \frac{1}{2}$) for $C \to 1$ and $n_{eff} \to 0$, with $n_{eff} = n_{th} + n_{\alpha}$.

## References

1. Aspelmeyer, M., Kippenberg, T. J. & Marquardt, F. Cavity optomechanics. *Rev. Mod. Phys* **86**, 1391 (2014).

2. Mahboob, I., Nishiguchi, K., Fujiwara, A. & Yamaguchi, H. Phonon lasing in electromechanical resonator. *Phys. Rev. Lett.* **110**, 127202 (2013).

3. Verlot, P., Tavernarakis, A., Briant, T., Cohadon, P. F. & Heidmann, A. Backaction amplification and quantum limit in optomechanical measurements. *Phys. Rev. Lett.* **104**, 133602 (2010).

4. Khurgin, J. B., Pruessner, M. W., Stievater, T. H. & Rabinovich, W. S. Laser-rate-equation description of optomechanical oscillations. *Phys. Rev. Lett.* **108**, 223904 (2012).

5. Cohen, J. D. et al. Phonon counting and intensity interferometry of a nanomechanical resonator. *Nature* **520**, 522 (2015).

6. Hong, S. et al. Hanbury brown and twiss interferometry of single phonons from optomechanical resonator. *Science* **358**, 203 (2017).

7. Stannigel, K. et al. Optomechanical quantum information processing with photons and phonons. *Phys. Rev. Lett.* **109**, 013603 (2012).

8. Riedinger, R. et al. Remote quantum entanglement between two micromechanical oscillators. *Nature* **556**, 473 (2018).

9. Suh, J. et al. Mechanically detecting and avoiding the quantum fluctuations of a microwave field. *Science* **344**, 1262 (2014).

10. Jing, H. et al. Pt-symmetric phonon laser. *Phys. Rev. Lett.* **113**, 053604 (2014).

11. Wang, B. et al. Phonon laser in the coupled vector cavity optomechanics. *Scientific Reports* **5**, 282 (2018).

12. Bergeal, N. et al. Analog information processing at the quantum limit with a Josephson ring modulator. *Nature Physics* **6**, 296 (2010).

13. Bergeal, N. et al. Phase-preserving amplification near the quantum limit with a Josephson ring modulator. *Nature* **465**, 64 (2010).

14. Caves, C. M. et al. Quantum limits on phase-preserving linear amplifiers. *Phys. Rev. A* **86**, 063802 (2012).

15. Clerk, A. A. et al. Introduction to quantum noise, measurement, and amplification. *Rev. Mod. Phys* **82**, 1155 (2010).

16. Roy, A. & Devoret, M. Introduction to parametric amplification of quantum signals with Josephson circuits. *C. R. Physique* **17**, 740 (2016).

17. Metelmann, A. & Clerk, A. A. Quantum-limited amplification via reservoir engineering. *Phys. Rev. Lett.* **112**, 133904 (2014).

18. Nunnenkamp, A. et al. Quantum-limited amplification and parametric instability in the reversed dissipation regime of cavity optomechanics. *Phys. Rev. Lett.* **113**, 023604 (2014).

19. Liu, Y.-H. et al. Cycling excitation process: An ultra efficient and quiet signal amplification mechanism in semiconductor. *Appl. Phys. Lett.* **107**, 053505 (2015).

20. Mikkelsen, M. et al. Optomechanics with a position-modulated kerr-type nonlinear coupling. *Phys. Rev. A* **96**, 043832 (2017).

21. Bobrovska, N. et al. Interactive optomechanical coupling with nonlinear polaritonic systems. *Phys. Rev. B* **95**, 085309 (2017).

22. Rebic, S. et al. Giant kerr nonlinearities in circuit quantum electrodynamics. *Phys. Rev. Lett.* **103**, 150503 (2009).

23. Teufel, J. D. et al. Sideband cooling of micromechanical motion to the quantum regime. *Nature* **475**, 359 (2011).

24. Gong, Z. R. et al. Effective hamiltonian approach to the kerr nonlinearity in an optomechanical system. *Phys. Rev. A* **80**, 065801 (2009).

25. Xiong, W. et al. Cross-kerr effect on an optomechanical system. *Phys. Rev. A* **93**, 023844 (2016).

26. Aldana, S., Bruder, C. & Nunnenkamp, A. Equivalence between optomechanical system and a kerr medium. *Phys. Rev. A* **88**, 043826 (2013).

27. Thompson, J. D. et al. Strong dispersive coupling of a high-finesse cavity to a micromechanical membrane. *Nature* **452**, 72 (2008).

28. Lu, X.-Y. et al. Steady-state mechanical squeezing in an optomechanical system via duffing nonlinearity. *Phys. Rev. A* **91**, 013834 (2015).

29. Xu, X.-W. et al. Mechanical pt-symmetry in coupled optomechanical systems. *Phys. Rev. A* **92**, 013852 (2015).

30. Bemani, F. et al. Synchronization dynamics of two nanomechanical membranes within a fabry-perot cavity. *Phys. Rev. A* **96**, 023805 (2017).

## Acknowledgements

This work was supported by the European Commission FET OPEN H2020 project PHENOMEN-Grant Agreement No. 713450.

## Author Contributions

P.D. carried out all calculations under the guidance of Y.P. and B.D.R. P.D. and B.D.R. participated in the discussions and interpretation of the results and to the writing of the manuscript.

## Additional Information

**Competing Interests:** The authors declare no competing interests.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
<figure><img src="image_5.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International

License, which permits use, sharing, adaptation, distribution and reproduction in any medium or
format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Cre-
ative Commons license, and indicate if changes were made. The images or other third party material in this
article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the
material. If material is not included in the article's Creative Commons license and your intended use is not per-
mitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the
copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2019