1Department of Applied Physics, East China Jiaotong University, Nanchang 330013, China
2Faculty of Information Technology, College of Microelectronics, Beijing University of Technology, Beijing 100124, China
3Center for Emergent Matter Science (CEMS), RIKEN, Wako, Saitama 351-0198, Japan
4Key Laboratory of Low-Dimensional Quantum Structures and Quantum Control of Ministry of Education, Department of Physics and Synergetic Innovation Center for Quantum Effects and Applications, Hunan Normal University, Changsha 410081, China
5Department of Physics, Zhejiang Sci-Tech University, Hangzhou 310018, China
We propose to achieve nonreciprocal quantum control of photons in a quadratic optomechanical (QOM) system based on directional nonlinear interactions. We show that by optically pumping the QOM system in one side, the effective QOM coupling can be enhanced significantly in that side, but not for the other side. This, contrary to the intuitive picture, allows the emergence of a nonreciprocal photon blockade in such optomechanical devices with weak single-photon QOM coupling. Our proposal opens up the prospect of exploring and utilizing quantum nonreciprocal optomechanics, with applications ranging from single-photon nonreciprocal devices to on-chip chiral quantum engineering.
1. INTRODUCTION
Optomechanical systems can be used to test the foundations of physics, and have important applications ranging from gravitational wave detections to quantum information processing (for reviews, see Refs. [1–6]). Recently, optomechanically induced nonreciprocity has been proposed theoretically [7–18] and demonstrated experimentally [19–28]. Nonreciprocal devices [29], such as isolators and circulators, have drawn an immense amount of interest in the past few years, for their irreplaceable role in signal processing and communication. An optomechanical system opens a route for nonreciprocity with nonlinear interactions [7,8,19,20], inherent non-trivial topology [9,21–24], and synthetic magnetism [10–18,25–28].
In this paper, we show that nonreciprocity can be realized in a whispering-gallery-mode (WGM) optomechanical system with quadratic optomechanical (QOM) coupling for quadratic mechanical-position dependent change in the optical frequency. Specifically, by optically pumping the quadratic WGM optomechanical system in one side, the effective QOM coupling can be enhanced significantly in that side, but not for the other side. The directional nonlinear interactions can induce nonreciprocal photon transport. This new possibility, as far as we know, has not been revealed in previous works.
It has been shown that the QOM systems driven by a strong optical driving field can generate strong (second-order) nonlinear photon-phonon interaction, even under weak single-photon optomechanical coupling conditions [30,31]. Nevertheless, we show here that the WGM optomechanical system with QOM coupling can be used to generate directional nonlinear interactions, which opens up the prospect of exploring nonreciprocal photon transport based on nonlinearity without the limitations of dynamic reciprocity [32]. Moreover, we demonstrate that directional nonlinear interactions can not only induce classical nonreciprocity, but also achieve nonreciprocal quantum control of photons by manipulating the statistic of the nonreciprocal transport photons.
Sign up for Photonics Research TOC. Get the latest issue of Photonics Research delivered right to you!Sign up now
Recently, a quantum nonreciprocal effect called nonreciprocal photon blockade (PB) was predicted theoretically [33] in a spinning Kerr resonator [34–37], that PB can emerge when the resonator is driven in one direction but not the other. Different from the nonreciprocal PB induced by the Fizeau–Sagnac drag [33], which relies on the condition of strong single-photon nonlinearity, here we show that the nonreciprocal PB can also be explored even with weak single-photon QOM coupling. Our proposal provides a feasible method to realize nonreciprocal photon blockade.
We note that the WGM optomechanical system has already been used for nonreciprocity with inherent non-trivial topology [9,21–24]. However, the WGM optomechanical systems with linear mechanical-position dependent change in the optical frequency (i.e., linear optomechanical coupling) can only enhance linear photon-phonon interaction, which cannot be used to observe nonreciprocal PB. Our work provides not only a new platform (quadratic WGM optomechanics) but a new approach (i.e., directional nonlinear interaction) to designing nonreciprocal single-photon devices for backaction-immune quantum communication and chiral quantum photonics.
2. MODEL AND HAMILTONIAN
We note that the cavity frequency of a WGM optomechanical system is almost linearly proportional to the mechanical position [38,39], even though the effects of QOM coupling have also been observed in experiments [40,41]. The QOM coupling can facilitate nondestructive measurements of the energy of the mechanical oscillator [42,43], which is instead impossible for a system with linear optomechanical coupling. We present a scheme to generate QOM interactions in the normal optical modes of a WGM optomechanical system by eliminating the linear optomechanical couplings [44–47].
We consider a near-field cavity optomechanical setup consisting of one mechanical resonator optomechanically coupling to two optical resonators () via the optical evanescent field, with each optical resonator supporting two degenerate clockwise (CW) and counter-clockwise (CCW) travelling-wave WGMs, as shown in Figs. 1(a) and 1(b). The model can be described by the optomechanical interaction Hamiltonian where and ( and ) are the annihilation and creation operators of the optical modes with frequency , and and are the dimensionless displacement and momentum operators of the mechanical resonator with frequency . As shown in Figs. 1(a) and 1(a), the mechanical resonator is placed at the middle of the gap of the two WGM resonators, so that the mechanical mode couples to the two WGM optical modes with linear optomechanical coupling of equal magnitude but opposite signs. Due to the tunneling coupling between the optical modes, there is an avoided crossing for the optical resonance frequency, , as shown in Figs. 1(c) and 1(d). We assume that the optical mode is coupled to a waveguide (port ) with strength , and the damping rate of mechanical resonator is . This theoretical model may also be realized with a membrane suspended on WGM microcavities [48] or photonic-crystal microcavities [49].
Figure 1.(a), (b) Schematic diagram for generating QOM coupling, where a mechanical nanostring oscillator is placed between two whispering gallery mode (WGM) resonators. (c), (d) Dispersion of the optical modes as a function of the displacement.
Following the approach in Refs. [44–47] (see Appendix A for more details), where is assumed such that can be treated as a quasi-static variable, the Hamiltonian can be diagonalized quasi-statically as for normal modes and with eigenfrequencies and single-photon QOM coupling .
As already shown in the experiment [47], when the tunneling coupling between the optical modes is strong, i.e., , the transmission spectrum of a laser probe through the coupled optical modes features resonance dips at the normal resonance frequencies , not at the bare optical resonance frequencies . That is to say, the normal modes are coupled to the external waveguides and can be used to describe the input–output characteristic of the coupled optical modes system. Specifically, the normal modes and are coupled to both of the ports with strength , respectively, and the total damping rates are .
We can choose either pair of degenerated quasi-static normal optical modes ( and , or and ), to generate strong nonlinear interaction for a few photons travelling in one direction, but not in the reverse direction. Without loss of generality, the two optical modes are denoted as and , with frequency or , for photons travelling from port 1 to port 2 and the opposite direction, respectively. For example, to enhance the nonlinear optomechanical coupling between the optical mode and the mechanical resonator , the optical mode is pumped by a strong field with amplitude and frequency . In the rotating reference frame with the optical frequency , the system can be described by with detuning . Under the strong driving condition , we perform the displacement transformations: , , , and , where , , , and are the steady state values, and , , , and (on the right side of the arrow) are the quantum fluctuation operators. Adopting the standard linearization technique (see Appendix B), an effective Hamiltonian [30,31] for the quantum fluctuation operators is obtained as where is the effective mechanical frequency, and is the effective (second-order) nonlinear coupling strength between the optical and mechanical modes. Without loss of generality, is assumed to be real in the following.
To investigate the system’s response behavior to a weak probe field, a weak field with amplitude and frequency is input from one of the ports. In the rotating reference frame with the unitary operator , the total Hamiltonian is given by , where with detunings , , and satisfying the condition . describes the driving of the probe field with subscript for weak field input from port 1 and for weak field input from port 2. The effective Hamiltonian describes a directional nonlinear interaction between the optical and mechanical modes. If a probe field is input from port 1, the enhanced nonlinear interaction between and may suppress the excitation of the second photon under a single-photon resonant condition. In contrast, by driving from port 2, as there is no strong nonlinear interaction between and , the photons are free to cross the system.
According to the input–output relations [50], we have and ( and ), and then the transmission coefficient for the weak probe field can be defined by for photon transport from port 1 to port 2, and for photon transport from port 2 to port 1, where and are the mean photon numbers. The isolation for probe field transport from port 1 to port 2 is defined by . Using the input–output relations, the statistic properties of the transmitted photons and can be described by the second-order correlation functions in the steady state () as follows: for photon transport from port 1 to port 2, and for photon transport from port 2 to port 1.
3. NONRECIPROCAL PB
The transmission coefficients and correlation functions can be obtained by numerically solving the master equation for the density matrix of the system [51] where denotes a Lindbland term for an operator , and is the mean thermal phonon number, given by the Bose–Einstein statistics with the Boltzmann constant and the temperature of the reservoir at the thermal equilibrium.
In Fig. 2(a), we show the transmission coefficients for probe field transport from port 1 to port 2, and for probe field transport from port 2 to port 1 versus the detuning . For , there are two peaks at and one dip at . In contrast, there is only one peak at for . Figure 2(b) shows the isolation for probe field transport from port 1 to port 2 as a function of the detuning . The isolation for the direction is more than 17 dB at detuning , and the isolation for the reverse direction , i.e., , is more than 70 dB at detuning .
Figure 2.(a) The transmission coefficients (solid black curve) and (dashed red curve) as a function of the detuning . (b) The isolation as a function of the detuning . (c) The equal-time second-order correlation function () as a function of the detuning . (d) The second-order correlation function as a function of the normalized time delay at detuning . The other parameters are , , , , and .
To explore the statistic properties of the transmitted photons, the equal-time second-order correlation function () is shown as a function of the detuning in Fig. 2(c). The photon transport from port 2 to port 1 is coherent in full frequency, i.e., . The photon transport from port 1 to port 2 exhibits strong antibunching effect, i.e., , around the detunings , and exhibits strong bunching effect, i.e., , around the detunings and . The time duration for nonreciprocal PB at is on the order of , as shown in Fig. 2(d).
The peak for at detuning can be understood by the fact that when the probe field is injected from port 2, only the (linear) optical mode can be excited. Thus the maximum transmission coefficient is reached for the probe field in resonance with the optical mode , i.e., , and the transmitted photons keep the statistic properties of the probe field (a coherent field), i.e., , as there is no nonlinear interaction in the optical path from port 2 to port 1.
In order to understand the origin of the peak for around the detuning and the dip for at , we show energy spectrum (see Appendix C) of the linearized QOM coupling between optical mode and mechanical resonator in Fig. 3. Under weak probe condition, the maximum transmission coefficient is reached for the probe field in resonance with the transition , i.e., , as shown in Fig. 3(a). However, the photon absorbed in the transition blocks the transition for large detuning, so we have around the detuning . The dip for at arises from the quantum interference between the transitions and [Fig. 3(b)], in an equivalent picture as optomechanically induced transparency [52–54], or electromagnetically induced transparency in lambda-type three-level atoms [55,56]. Moreover, when , the transition is suppressed, but the two-photon transition is resonant, which induces two-photon tunneling from port 1 to point 2, i.e., . Similarly, around is induced by the resonant transition .
Figure 3.Schematic energy spectrum of the linearized QOM coupling between optical mode and mechanical resonator , where , , , , , , and represents the Fock state with photons in and phonons in .
We now discuss the effect of thermal phonons on the nonreciprocal PB. Figures 4(a) and 4(b) show the transmission coefficient and the equal-time second-order correlation function versus the detuning with a different mean thermal phonon number (). Thermal phonons have little influence on the transmission coefficient around , but have great effect on the transmission coefficient around . Thermal phonons have great effect on the second-order correlation function around and .
Figure 4.(a) Transmission coefficient . (b) The equal-time second-order correlation function versus the detuning with different mean thermal phonon number (). (c) The isolation . (d) The equal-time second-order correlation function versus the mean thermal phonon number with different detuning (). The other parameters are the same as in Fig. 2.
The isolation and the second-order correlation function versus the mean thermal phonon number are shown in Figs. 4(c) and 4(d) with detuning . The isolation around is robust against the thermal phonons, but the antibunching effect of the transport photons becomes much weaker for greater thermal phonons. Both the isolation and second-order correlation function around are sensitive to the mean thermal phonon number , and this quality may be used in accurate temperature measurement at ultra-low temperature.
More interestingly, two peaks appear around in the transmission coefficient , and the isolation can be improved with a larger thermal phonon number . This abnormal effect is induced by the phonon states, e.g., in Fig. 3(c). As the temperature increases, the population probability in phonon state increases, and the transitions of with resonance frequency become remarkable gradually, which induces the increasing peaks of the transmission coefficient (or the isolation ) around .
4. DISCUSSIONS AND CONCLUSIONS
Let us now discuss the experimental requirements for our proposal. Considering the parameters from a recent near-field WGM optomechanical experiment [39]: , , and for strong coupling, we have QOM coupling , which is too weak to realize . However, many schemes have been proposed to enhance the optomechanical coupling in recent years, such as by using the collective mechanical modes in scatter arrays [57], the squeezed cavity mode [58], or the Josephson effect for electromechanical systems [59]. The optomechanical coupling can be enhanced by several orders of magnitude and the strong coupling regime () can even be approached [57–59]. Furthermore, the frequency range of the SiN nanostring oscillators in the near-field cavity optomechanical setup is about 0.1–10 MHz [38–41], which is too low to observe nonreciprocal PB. One possible solution is replaying the SiN nanostring oscillator with an optomechanical crystal nanobeam, as optomechanical crystal nanobeams with GHz mechanical modes have been designed and demonstrated recently [60,61]. Moreover, the thermal phonons have important influence on the nonreciprocal PB, and the thermal effect can be suppressed for higher frequency mechanical mode.
In addition, there are two driving fields coupling to the system with frequencies and , and the strengths . In order to detect the photon correlation induced by the weak probe field, we should spectrally filter out the strong optical driving field at frequency under the condition . A similar method has already been used to measure correlations of phonons in a recent experiment [62].
In summary, we have shown that the WGM optomechanical system with QOM coupling can be used to obtain directional nonlinear interactions and observe nonreciprocal PB. In addition, the thermal phonons have important influence on the nonreciprocal PB, especially on the statistic properties of the transport photons, and this quality may be used in temperature sensing at ultra-low temperature. Moreover, this work can be extended to study phonon manipulation in double-cavity optomechanics, e.g., nonreciprocal phonon blockades [63], nonreciprocal phonon lasers [64–70], nonreciprocal photon-phonon entanglements and quantum transfers. Directional nonlinear interaction can also be applied in the high-Q microtoroid resonator with phase-matched parametric amplification [71] or superconducting microwave circuit with Josephson parametric converters [72] to achieve nonreciprocal quantum control of photons. Our proposal provides a new routine towards the realization of on-chip quantum nonreciprocal devices, backaction-immune quantum measurement, and chiral quantum physics.
APPENDIX A: QUADRATIC OPTOMECHANICAL COUPLING
To realize quantum nonreciprocality in optomechanics, we consider a near-field cavity optomechanical setup consisting of one mechanical resonator optomechanical coupling to two optical resonators () via the optical evanescent field, with each optical resonator supporting two degenerate clockwise (CW) and counter-clockwise (CCW) travelling-wave whispering-galley-modes (WGMs). The model can be described by the optomechanical interaction Hamiltonian where and ( and ) are the annihilation and creation operators of the optical modes with frequency , is the tunneling amplitude between the optical modes, and are the dimensionless displacement and momentum operators of the mechanical resonator with frequency , and is the linear optomechanical coupling strength between the mechanical resonator and optical modes. We assume that the optical mode is coupled to a waveguide (port ) with strength , and the damping rate of mechanical resonator is .
The Hamiltonian in Eq. (A1) can be diagonalized as in the normal modes basis, and with and eigenfrequencies
Moreover, is assumed [44–47] such that we can Taylor-expand the eigenfrequencies in Eq. (A6) as with frequencies , and the quasi-static Hamiltonian in Eq. (A2) is given approximately by Eq. (2) in the main text with quadratic optomechanical coupling , and quasi-static normal optical modes, and .
APPENDIX B: LINEARIZED OPTOMECHANICAL HAMILTONIAN
We consider a pair of degenerated quasi-static normal optical modes ( and , or and ), to generate strong nonlinear interaction for a few photons travelling in one direction, but not in the reverse direction. Without loss of generality, the two optical modes are denoted as and , with frequency or , for photons travelling from port 1 to port 2 and the opposite direction, respectively. For example, to enhance the nonlinear optomechanical coupling between the optical mode and the mechanical resonator, the optical mode is pumped by a strong field with amplitude and frequency . In the rotating reference frame with the optical frequency , the system can be described by Eq. (3) in the main text with detuning .
Under strong driving condition , we perform the displacement transformations: , , , and , where , , , and are the steady state values, and , , , and (on the right side of the arrow) are the quantum fluctuation operators. The steady state values , , , and can be obtained by the equations of motions yielding and
The operators and for the mechanical resonator can be written in terms of phonon creation and annihilation operators as , , and the effective Hamiltonian for the quantum fluctuation operators reads
We assume that the optical driving field is strong enough so we have , and so that the term in the above equation can be neglected safely. For , the effective Hamiltonian can be further simplified by rotating-wave approximation and neglecting the high frequency terms [30,31], e.g., and . Then, we get Eq. (4) in the main text, where is the effective mechanical frequency, and is the effective nonlinear coupling strength between the optical and mechanical modes.
APPENDIX C: ENERGY SPECTRUM
In this appendix, we will show the energy spectrum of the effective Hamiltonian in the diagonal basis under the condition , and similar results are shown in Refs. [30,31]. In the noncoupling basis, represents the Fock state with photons in and phonons in . It is clear that the vacuum state and the single-phonon state are decoupled from the other states. The degenerate states and are coupled with strengh , so the eigenstates are with energy splitting
The degenerate states and are coupled with strength , so the eigenstates are with energy splitting
The states , , and are degenerate, where and are coupled with strength , and and are coupled with strength , so the eigenstates are with energy splitting
The schematic energy spectrum in the diagonal basis is given in Fig. 3.
[49] H. Wang, Q. Qiao, C. Peng, J. Xia, G. Zhou, Y.-J. Zhao, X.-W. Xu. Two-dimensional optomechanics formed by the graphene sheet and photonic crystal cavity(2018).