Introduction

Strong optical nonlinearity gives rise to many important quantum effects, such as photon blockade1,2,3, quantum squeezing4, quantum nondemolition measurements5,6, optical switching with single photon7 and so on8,9,10. These nonlinear optical effects have been demonstrated in cavity QED systems, where the quantum coherence in the atom1,2,3 (or artificial atom11,12,13,14,15,16,17) generates strong effective photon nonlinearities.

Recently, cavity optomechanics has become a rapidly developing research field exploring nonlinear coupling via radiation pressure between the electromagnetic and mechanical systems18,19,20. It has been shown theoretically that strong optical nonlinear effects (and relevant applications, such as generating nonclassical state, photon blockade, multiple sidebands, photon-phonon transistors and optomechanical photon measurement) can be realized in single-mode21,22,23,24,25,26,27,28,29,30,31,32,33 or two-mode optomechanical systems (OMSs)34,35. These phenomena are mainly demonstrated in the single-photon strong-coupling regime, where the optomechanical coupling strength at the single-photon level ga exceeds the cavity decay rate κa (ga > κa). However, in most experiments to date36,37,38, ga is much smaller than κa (gaa ~ 10−3). Only a few new-type optomechanical setups, using ultracold atoms in optical resonators (gaa ~ 10−1)39 or optomechanical crystals (gaa ~ 10−2)40, can one begin to approach the single-photon strong-coupling regime. On the other hand, a strong optical driving field may effectively enhance the optomechanical coupling by a factor , where n is the mean photon number in the cavity41,42,43. But this enhancement comes at the cost of losing the nonlinearity of the interactions. Specifically, under the condition of strong optical driving, the linearized coupling strength between the optical and mechanical modes is largely enhanced, which makes the intrinsic nonlinear optomechanical coupling smaller and negligible.

Given the above, it is highly desirable to find a new method for obtaining strong Kerr nonlinearities in OMSs in the weak-coupling regime, namely the optomechanical coupling strength is much smaller than the optical cavity decay rate (). In this paper, we investigate the Kerr nonlinear effects of the optical field in a hybrid electro-optomechanical system containing a mechanical oscillator coupled to both an optical cavity and a microwave LC resonator (see Fig. 1)44,45,46,47. We find that the eletromechanical subsystem (the mechanical oscillator plus the microwave resonator) displays a quantum criticality. One can drive the electromechanical subsystem close to the quantum critical regime by applying a microwave field with properly chosen frequency and intensity to the microwave resonator. Then the quantum criticality can induce a strong Kerr nonlinearity in the optical cavity, even if the optomechanical systems (the optical cavity and mechanical oscillator) is in the weak-coupling regime. This strong Kerr nonlinearity can be demonstrated by the existences of photon blockade and nonclassical states (e.g., Schrödinger cat states) of the cavity field when the electromechanical subsystem approaches the quantum critical point. Furthermore, the strong Kerr nonlinearity can also be controlled easily by tuning the intensity (or frequency) of the microwave driving field. This provides a promising route for experimentally observing strong Kerr nonlinearities in OMSs in the weak-coupling regime.

Figure 1
figure 1

Schematic diagram of the hybrid electro-optomechanical system.

A mechanical oscillator couples to both an optical cavity and a microwave LC resonator.

Results

Hybrid electro-optomechanical system

In the hybrid electro-optomechanical system of Fig. 1, the mechanical oscillator is parametrically coupled to both the optical cavity and the microwave resonator. The microwave resonator is driven by a strong field with amplitude εc and frequency ωci, where εc is related to the input microwave power P and microwave decay rate κc by . In a frame rotating with frequency ωci, the Hamiltonian for this hybrid systems reads48

where the detuning δc = ωc − ωci and the microwave frequency , ga (gc) denotes the optomechanical (electromechanical) coupling strength at the single-photon level and ( or ) is the annihilation operator of the optical cavity (the mechanical oscillator or the microwave resonator). Under a strong microwave driving field, following the standard linearization procedure49,50,51,52 (shifting and with their steady-state mean values α and β, i.e., , ), the Hamiltonian can be transformed into

where G is the linearized electromechanical coupling strength; Δc and are, respectively, the effective microwave detuning and optical frequency including the radiation-pressure-induced optical resonance shift. Their explicit expressions are given by

Notice that G and Δc can be easily controlled by tuning the power and frequency of the microwave driving field.

Quantum critical property of the electromechanical subsystem

The quantum criticality in the electromechanical subsystem can be shown by diagonalizing the electromechanical subsystem via a Bogoliubov transformation . Here, the canonical operators are and and M is the transformation matrix given by

where the angle θ is defined by

Then, the Hamiltonian becomes

where ω± are the normal mode frequencies of the electromechanical subsystem,

and

are the effective coupling strengths between the optical photon and the normal modes. Equation (5) shows that becomes zero (negative) when

as shown in Fig. 2(a). This corresponds to a critical property53, namely, the normal mode ω will change from a standard harmonic oscillator (G < Gcp) to a free particle and further becomes dynamically unstable (G > Gcp) as G crosses its critical value Gcp.

Figure 2
figure 2

Quantum criticality of the electromechanical subsystem and strong Kerr nonlinearity of the optical field.

(a,b) Quantum criticality of the electromechanical subsystem, characterized by the normal-mode frequency ω±b. As one can see, the normal-mode ω+ continuously passes through the critical point. The quantum criticality is manifested with the normal-mode ω, which becomes purely imaginary after the critical point Gb > 0.5 and Δcb < 1.25 (c,d) Strong Kerr-nonlinearity given by the photon-photon interaction strength η in the optical cavity, as a function of the adjustable parameters G and Δc controlled by the microwave driving field. The pink circles and shaded area in (c,d) correspond, respectively, to the regimes η = κa and η > κa. The black dot-dashed vertical lines indicate the quantum critical points Gcp and Δcp. Other system parameters are taken as: ωb/2π = 10 MHz, gab = gcb = 10−3, κab = 0.1, κcb = 0.127, Δcb = 1.251 (a,c) and Gb = 0.5595 (b,d).

The above critical property can become more transparent with the following canonical relationships:

Here xb, xc are the dimensionless displacements of the mechanical and microwave resonators from their stable points and pb, pc are the corresponding dimensionless momentums. The Hamiltonian of the electromechanical system can then be written in terms of the usual canonical xp variables, He–m = H0 + Hint with

denoting the free Hamiltonian of the microwave and the mechanical resonators and the coupling between them. The potential of the free Hamiltonian (6a) can be further expressed as

It shows that the intrinsic potential of the electro and mechanical resonators is characterized by . Comparing Eq. (7) with the coupling Hamiltonian (6b), one can see that there is an interplay between the intrinsic potential and the coupling interaction between them. This interplay leads to the above critical property. In other words, when G approaches (or exceeds) , the normal mode ω is dragged out from its effective potential and becomes increasingly flat (or inverted) [see the Fig. 3].

Figure 3
figure 3

The critical property of the electromechanical subsystem.

(a) The mechanical and electrical modes couple with each other with the coupling strength G. The black circle indicate the quantum critical point. (b,c,d) The effective potential of the normal mode ω becomes flat and further inverted when increasing the coupling strength G.

Quantum-criticality-induced strong Kerr nonlinearities

As one can see, the last two terms in the Hamiltonian (4) show that optical photons can interact with each other through the exchange of the normal modes , very similar to electrons interacting with each other through the exchange of photons in the QED Hamiltonian. In particular, when the electromechanical subsystem approaches its quantum critical point, the optical cavity shows a strong effective Kerr nonlinearity. This quantum-criticality-induced strong Kerr nonlinearity becomes clear after taking a displacement transformation, , where is a similarity transformation and with . The result is

and η is the photon-photon interaction strength,

Notice that the photon-photon interaction strength η remains unchanged when the system-environment interaction is explicitly included (see the detailed derivation in Methods). On the other hand, it also shows in Figs. 2(c,d) that even in the weak-coupling regime , a strong photon-photon interaction η (η > κa) can still be obtained when G (or Δc) approaches the quantum critical point. In particular, Fig. 2 shows that when the coupling strength G (or the detuning Δc) is close to its quantum critical point, a very small normal mode frequency ω is obtained, which induces a large photon-photon interaction with η 1/ω. The system parameters G and Δc, determined by the power P and the frequency detuning δc of the input microwave driving field, can be directly tuned in experiments. Figure 4 shows explicitly the practical parameter range of P and δc for obtaining the strong Kerr nonlinear parameter η (η > κa).

Figure 4
figure 4

Controllability of the system parameters.

Nonlinear parameter η versurs: (a) the coupling strength G and (b) the effective detuning Δc. The inserts present the experimental parameter range (the power P and frequency detuning δc of the input microwave field) for η > κa. The system parameters are the same as those in Fig. 2 except for ωc/2π = 7.5 GHz for the inserts.

Photon blockade

The strong Kerr nonlinearity in the present system can be further demonstrated by the steady-state second-order correlation function of the optical field g(2)(0). g(2)(0) → 0 in the weak-coupling regime signals the photon blockade effect and can be directly detected by a Hanbury-Brown-Twiss Interferometer3. Explicitly, by driving the optical cavity with a weak laser field of frequency ωai and amplitude εa, the Hamiltonian of the system becomes

where all the similarity transformations used before have been taken into account and . The damping effect arising from the coupling of the optical field to the electromagnetic vacuum modes of the environment can also be taken into account and the dissipative dynamics of cavity mode is described by the quantum Langevin equation,

Here κa is the decay rate of cavity mode and is a vacuum noise operator satisfying , .

With a weak optical driving field, the quantum Langevin equation is solved by truncating them to the lowest relevant order in εa. The resulting two-photon correlation function is given by

with

where

is the normalized s-photon probability in the cavity ( in the weak-driving regime) and

Note that is a complex operator including the microwave field and the mechanical mode and is determined by the dynamics of the electromechanical modes Bj, which evolves as

The noise operator , which comes from the environment of the microwave resonator. The environment is initially in the thermal equilibrium state ρth with temperature T and is the initial environment operators of the microwave resonator. Here, we have safely ignored the dissipation of the mechanical oscillator because the mechanical decay rate κb is extremely small, κba, κbc < 10−3. Thus, the effective decay rates κj is determined by the original decay rate of the microwave resonator κc (see the detailed derivation in Methods).

In Fig. 5, we show the dependences of κ± on the system parameters G, Δc and κc. From Fig. 5 one can see that the effective decay rate κ sharply changes from a positive value to a negative value when the system parameter G (or Δc) crosses its quantum critical point Gcp (or Δcp). This result demonstrates that the mode ω will become a gain mode when G > Gcp or Δc < Δcp. Near the quantum critical points Gcp and Δcp, the effective decays κ± almost become constant with G or Δc [see the inserts of Fig. 5(a) and 5(c)]. In Fig. 5(b) κ± is plotted via the microwave field decay rate κc when G (or Δc) is near the quantum critical points. As it is shown, κ± exhibit a linear increase with the decay rate of the microwave field κc.

Figure 5
figure 5

Effective decay rates for the electromechanical normal modes.

Effective decay rates κ± of the normal modes ω± versus: (a) the coupling strength G, (b) the decay rate of microwave field κc and (c) the detuning Δc. The system parameters are Δc/2π = 12.51 MHz, κc/2π = 1 MHz for (a), while G/2π = 5.5924 MHz, Δc/2π = 12.51 MHz for (b) and G/2π = 5.595 MHz, κc/2π = 1 MHz for (c). The purple circles in (b) indicate the parameter regime corresponding to the κ± used in Fig. 6.

When the microwave (mechanical) mode is initially in the coherent state |α〉 (|β〉) and the optical field in the vacuum state, the two-point correlation function exp(−Φ2) and the four-point correlation function exp(−Φ4) are calculated. With numerically integrating Eqs. (13), the dependence of g(2)(0) on κ, G and Δc is shown in Fig. 6. As one see, in the quantum critical regime, the photon antibunching effect g(2)(0) < 1 (even the photon blockade g(2)(0) → 0) occurs because the two-photon transition is largely suppressed in comparison with the single-photon transition when κ/2π > 60 kHz [see the insert in Fig. 6(a)]. Figures 6(b) and (c) further show that the photon blockade [g(2) (0) → 0] occurs when the tunable parameter G (or Δc) approaches its quantum critical value even if the optomechanical coupling ga is very weak.

Figure 6
figure 6

Photon statistics in the hybrid electro-optomechanical system.

Equal-time second-order correlation function g(2)(0) versus: (a) effective decay rate κ, (b) coupling strength G and (c) detuning Δc at T = 0 for the microwave bath. The red circles in (a) indicate the value of κ for g(2)(0) = 1. We have chosen the quantum critical parameters: G/2π = 5595 kHz, Δc/2π = 12521.64 kHz in (a) and the decay rates κ/2π = (500, 250, 50) kHz (corresponding to κc/2π = (1270, 620, 110) kHz) in (b,c). The other parameters are the same as in Fig. 2, except for Δa = η in order to maximize the photon antibunching effect.

Furthermore, we also find that the photon antibunching [g(2)(0) < 1] disappears when κ/2π < 60 kHz [see the inserts in Figs. 6(b) and (c)]. Physically, this is because in the hybrid OMS, a relatively large decay rate κ/2π > 60 kHz) occurs when the electromechanical subsystem approaches the quantum critical point. This decay will significantly suppress the steadystate sideband transition in the electromechanical subsystem. Then, in the quantum critical regime, the hybrid OMS becomes a pure optical nonlinear system and η > κa is the exclusive condition for achieving the photon blockade. Meanwhile, the very small ω → 10 kHz) near the quantum critical point effectively enhances the photon-photon interaction to η > κa because η 1/ω, namely the photon blockade can still be reachable even if the effective electromechanical frequency extends beyond the resolved sideband regime, i.e. ω < κa. Notice that the original mechanical frequency used here is still in the resolved sideband regime () so that there is no problem in cooling the mechanical oscillator at the initial time.

Nonclassical states

As demonstrated in previous studies21,22,23, strong Kerr nonlinearities generally lead to the periodic generation of nonclassical states, (e.g., cat states) of the cavity field. With the help of the Hamiltonian (4), we obtain the time evolution operator in the interaction picture,

where the term corresponding to ζ+ has been omitted due to its negligible effect on the evolution of the cavity mode +b ~ 10−4) near the quantum critical point. If the cavity field is initially in a coherent state , the cavity field at time tn = 2 (n = 1, 2…) will be in the state

The state |Ψa(tn)〉 is a multi-component cat state, depending on the value of η/ω. Figure 7 shows the different multi-component cat states for different values of the tunable parameters G and Δc near the quantum critical point. Figures 7(b,c,d) present the specific realization of two-, three- and four-component cat states, respectively. Here damping effects (given by κa, κc, κb) have been ignored. In principle, this is valid when the cut-off time , 1/κc, 1/κb. The optical field damping is similar to that in a recent cavity-QED experiment54. Moreover, inspired by Ref. 54, the Wigner function can be measured (or reconstructed) by detecting the states of the atoms interacting with the optical field. Nevertheless, the above result indicates that the quantum-criticality-induced strong Kerr nonlinearities in this hybrid OMS can generate nonclassical states by cutting off the optomechanical interaction at the appropriate time, which can be detected via Wigner tomography.

Figure 7
figure 7

Parameter regimes (a) for obtaining the two- (b), three- (c) and four-component (d) Schrödinger cat state.

The quadratures variables are , . The system parameters are the same as in Fig. 2 except for .

Discussion

We have provided a new mechanism for obtaining strong Kerr nonlinear effects in a hybrid OMS in the weak-coupling regime. We found that the electromechanical subsystem displays a critical property when adjusting the intensity (or frequency) of the microwave driving field and a strong controllable photon-photon interaction is induced in the quantum critical regime. Usually, the phonon modulation effect influences the photon statistics in the usual OMSs24 and in general will also weaken the photon-photon interaction effect, except in the single-photon strong-coupling (ga > κa) and the resolved sideband () regime24. The essence of the strong photon-photon interaction presented in this paper can be understood as follows. In the quantum critical regime, the electromechanical normal mode coupled to the optical field is highly softened (or a very small normal-mode frequency ω is obtained). At the same time, the sideband phonon transitions in the electromechanical subsystem are significantly suppressed by the relative large decay rate of the electromechanical normal mode, which makes the hybrid OMS essentially a pure optical nonlinear system. Thus, the quantum-criticality-induced strong self-Kerr nonlinearity is very different from previous investigations in the usual OMSs24,34,35.

Experimentally, the strong photon-photon interaction achieved in the present hybrid OMS requires driving the electromechanical subsystem into its quantum critical region (shaded area in Fig. 2). Normal-mode splitting in the driven electromechanical system has been observed43. The quantum critical region could be easily reached by increasing the intensity of the microwave driving field. Moreover, as shown in Figs. 2 and 4, the interesting ranges of G and Δc are respectively on the order of 0.1 kHz and 1 kHz for the quantum critical region and this parameter precision is experimentally realizable55. We believe that our proposal will provide a new avenue for experimentally realizing strong optical nonlinearities in the weak-coupling regime and largely enrich the scope of implementing quantum information processing and quantum metrology with cavity OMSs.

Methods

Derivation of the photon-photon interaction with system-environment couplings

The total Hamiltonian of the hybrid OMS plus the environment can be written as

where the system Hamiltonian is given by Eq. (4) and

are the Hamiltonians of the environment and the system-environment interaction, respectively. Notice that the system-environment interaction is invariant to the linearization procedure applied on the electromechanical subsystem. Here , , are the bath operators for , , and Kj(ω) (j = a, b, c) are the corresponding coupling constants. For a slowly-varying bath spectrum, we can simply replace Kj(ω) by the decay rate . Here the last term can be safely neglected because the decay rate κb of the mechanical oscillator is extremely small (κba, κbc < 10−3).

By applying a Bogoliubov transformation to the total Hamiltonian , the hybrid OMS Hamiltonian and the interaction between the system and the environment can be rewritten in terms of the normal-mode canonical operators

while the environment Hamiltonian retains its original form.

To derive the photon-photon interaction, the total Hamiltonian should be further diagonalized in a displaced-oscillator representation, and the result is

where

is the photon-photon interaction strength. This similarity transformation also does not affect the environment Hamiltonian . By comparing with the dissipation Hamiltonian in the original representation [Eq. (18b)], it can bee seen that the last term of Eq. (19b) is induced by the similarity transformation in the displaced-oscillator representation and it may change the photon-photon interaction. However, we will show next that, in the quantum critical regime, this term will not change the photon-photon interaction η and it only induces a negligible pure-dephase of the optical mode.

In the quantum critical regime, the system parameters M13, M14 and ζ+ are negligible compared to the parameters M11, M12 and ζ, due to the relative large frequency ω+ (M13, M14 and ζ+ are smaller than M11, M12 and ζ by about 3 to 4 orders of magnitude). This means that the influence of the normal mode B+ on the dynamics of the optical mode can be safely neglected when the electromechanical subsystem approaches its quantum critical point. By ignoring the normal modes , the dynamics of the bath operator can be determined by the following equation of motion,

Solving Eq. (21), the result is

where is the initial environment operator of the microwave resonator. Substituting the above solution of the bath operator and its hermitian conjugate into the last term of Eq. (19b) and noticing that in the quantum-critical-regime , we have

Therefore, the photon-photon interaction given in Eq. (19a) remains invariant under the interaction with the environment.

Calculation of the effective decay rates for the electromechanical normal modes

Based on the dissipative dynamics of the electro-mechanical subsystem in the original representation, we can find the relationship between the effective decay rates κ± and the original decay rates of the microwave resonator κc.

Considering the thermal environments of the microwave resonator, the Hamiltonian of the electro-mechanical subsystem plus the environment is

Then, the dynamics of the canonical operator is given by

where the coefficient matrix

Here, Γ = diag(κc, κc, 0, 0) denotes the decay rates of the microwave resonator and the mechanical oscillator and are the Langevin forces. Equation (25) shows that the imaginary and real parts of the eigenvalues of D correspond to the eigenfrequencies ω± and the effective decay rates κ± of the normal modes. For the undamped case (κc = 0), the eigenvalues of D are purely imaginary and we obtain the expression Eq. (5) for the normal-mode frequencies. For the general κc, we numerically diagonalized the coefficient matrix D and shown the results in Fig. 5.