Next Article in Journal
Optical Coherence Tomography Following Panretinal Photocoagulation Demonstrating Choroidal Detachment
Next Article in Special Issue
Integer and Fractional Floquet Resonances in a Driven Three-Well System
Previous Article in Journal
Measuring the pth-Order Correlation Function of Light Field via Two-Level Atoms
Previous Article in Special Issue
Dynamics of Nano-Particles Inside an Optical Cavity in the Quantum Regime
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optical Nonreciprocity in Double Optomechanical Systems with Quadratic Coupling

Department of Physics, Zhejiang Sci-Tech University, Hangzhou 310018, China
*
Author to whom correspondence should be addressed.
Photonics 2022, 9(10), 728; https://doi.org/10.3390/photonics9100728
Submission received: 20 September 2022 / Revised: 1 October 2022 / Accepted: 2 October 2022 / Published: 5 October 2022

Abstract

:
In this paper, the optical nonreciprocal phenomena in double optomechanical systems with quadratic coupling are studied. Our model belongs to an optomechanical system in which three coupling modes coexist, that is, the two cavity fields are coupled with the mechanical oscillator at the same time and the couplings are in the form of quadratic interactions. In addition, there is a linear coupling mode between the two cavity fields. In the entire system, each cavity field is effectively coupled by a control field and a probe field simultaneously. The expression of the transmission coefficient of the probe field is obtained by solving the dynamic evolution equation satisfied by the system. Using numerical analysis, we analyze the change in transmission coefficient of the probe field under the conditions of different physical parameters. The results show that we can realize optical nonreciprocal transmission in this system. Appropriate choices about physical parameters can achieve perfect nonreciprocity. Our theoretical scheme to realize optical nonreciprocal transmission in a double optomechanical system provides a theoretical basis for optical circulators, cyclic amplifiers and directional amplifiers.

1. Introduction

Photons with good characteristics, such as fast response ability, strong parallel processing ability, no charge, strong anti-interference and so on, have become an effective carrier for information transmission and play an extremely important role in optical communication, which has been widely discussed by researchers. As is common knowledge in the field, in the process of optical signal transmission and processing, optical nonreciprocal devices [1,2] are important for information processing. Because they can control information transmission unidirectionally, they play an important role in optical communication. Compared with optical reciprocity, optical nonreciprocity is usually realized by breaking the time-reversal symmetry with the Faraday rotation effect [3,4] or optical nonlinear effect [5,6]. Devices with nonreciprocal properties can greatly simplify the construction of optical networks and are widely used in signal transmission. Therefore, the study of optical nonreciprocity is of great significance.
Early researchers realized optical nonreciprocity by using the magneto-optical effect to break the time-reversal symmetry of system, and then they realized the nonreciprocal transmission of the light field. However, the device obtained in this way has a big size and, thus, it is difficult to meet the requirements of integrated optical communication networks. Therefore, people try to find other ways to realize nonreciprocal transmission. We know that many strange optical phenomena can be obtained by using the nonlinear characteristics of optical systems. For example, recently, people have also successfully achieved optical nonreciprocity by using this nonlinear effect [7]. Shi and his coworkers analytically prove and numerically demonstrate a dynamic reciprocity in nonlinear optical isolators based on Kerr or Kerr-like nonlinearity [8]. Using a method based on optical nonlinearity, Qi’s group [9,10] uses two silicon rings 5 μm in radius; the device is passive yet maintains optical nonreciprocity for a broad range of input power levels. The research of Jing and his coworker [11] shows that in a spinning Kerr resonator, photon blockade can happen when the resonator is driven in one direction but not the other and this effect is called nonreciprocal photon blockade.
In recent years, cavity optomechanical systems [12,13,14,15,16], being a new class of microscale integrable devices, have attracted extensive attention from researchers. We know that cavity optomechanical systems can effectively couple mechanical degrees of freedom and optical degrees of freedom. Based on performing quantum-coherent manipulation on the system, the nonlinear optical effect of the system can be significantly improved, which provides an effective physical carrier for quantum information processing. Based on the cavity optomechanical system, researchers have achieved quantum nonlinear phenomena such as electromagnetically induced transparency [17,18], cooling of mechanical oscillators [19,20], preparation of quantum entanglement [21,22] between macroscopic mechanical oscillators and cavity fields, and photon blockade [23,24,25]. Based on the optomechanical system, some schemes about optical nonreciprocal phenomena have been reported [26,27,28,29,30]. By optically pumping the ring resonator in one direction, the optomechanical coupling is only enhanced in that direction, and consequently, the system exhibits a nonreciprocal response [26]. The optical nonreciprocal behavior of a three-mode optomechanical system is induced by the phase difference between the two optomechanical coupling rates, which breaks the time-reversal symmetry of the three-mode optomechanical system [27]. Chen (one of the authors) and his coworker achieve the nonreciprocal quantum control of photons in a quadratic optomechanical system by using directional nonlinear interactions [28].
In this paper, a scheme about optical nonreciprocal phenomena in double optomechanical systems with quadratic coupling is proposed. We investigate the change in the magnitude and direction of nonreciprocity when the effective phase difference of the system is changed. We consider that the system is in two different situations: red detuning and blue detuning. The conditions for realizing nonreciprocity are discussed. We find that, with appropriate parameter selection, we can achieve perfect nonreciprocity. Our theoretical scheme to realize optical nonreciprocal transmission in a double optomechanical system provides a theoretical basis for optical circulators, cyclic amplifiers and directional amplifiers. The sections of the paper are arranged as follows: Section 2, the physical model and the evolution of the system are introduced. Realization about the nonreciprocal transmission is discussed in Section 3, and some sums are shown in Section 4.

2. Model and Evolution of the System

We consider a double-cavity optomechanical system [27], in which a common mechanical oscillator is coupled with the two cavities simultaneously as shown in Figure 1, where c i ( ω c i ) represents the annihilation operator (frequency) of the i t h cavity and b ( ω m ) denotes the annihilation operator and vibration frequency of the mechanical resonator. The set-up in Figure 1 can be realized in a Fabry–Pérot optical cavity, which is composed of two parallel mirrors separated by a certain distance. The standing wave can be produced in the optical cavity when we use an incident light field to drive the optical cavity. A partially transparent dielectric membrane is placed in the optical cavity, and it can be viewed as a mechanical oscillator under the interaction of radiation pressure. If the membrane oscillates around the node, the quadratic term of optomechanical coupling strength is much larger compared to the other terms. In addition, the membrane makes the system become a double-cavity system. Each cavity and mechanical resonator are coupled by optomechanical interactions of the form of quadratic coupling with the quadratic coupling constant of g i = 2 ω c i 2 q x 0 2 ( i = 1 ,   2 ) . p and q represent the momentum operator and position operator of the mechanical resonator, respectively, and define the operator P = x 0 p , Q = q x 0 under the condition of canonical quantization. Here,   is the Planck constant, x 0 = m ω m is the zero-point fluctuation of the mechanical resonator displacement, and P and Q satisfy the relation of [ Q , P ] = i . In addition, the two cavities are coupled with a linear coupling constant J , and each cavity field is driven by a control field ϵ L i ( i = 1 ,   2 ) and a probe field ϵ P i ( i = 1 ,   2 ) simultaneously. The model can be implemented in the membrane-in-the-middle optomechanical system. The Hamiltonian can be written as ( = 1 )
= ω c 1 c 1 c 1 + ω c 2 c 2 c 2 + 1 2 ω m ( P 2 + Q 2 )   + J ( c 1 c 2 + c 2 c 1 ) + 2 g 1 c 1 c 1 Q 2 + 2 g 2 c 2 c 2 Q 2 + i ϵ L 1 ( c 1 e i ω L 1 t c 1 e i ω L 1 t ) + i ϵ L 2 ( c 2 e i ω L 2 t c 2 e i ω L 2 t ) + i ( ϵ P 1 c 1 e i ω P 1 t ϵ P 1 * c 1 e i ω P 1 t ) + i ( ϵ P 2 c 2 e i ω P 2 t ϵ P 2 * c 2 e i ω P 2 t ) .
The Hamiltonian of the system is comprised of five terms, namely the free Hamiltonian of the two cavities, the free Hamiltonian of the mechanical oscillator, the interaction between the two cavities, the interaction between the mechanical oscillator and the two cavities, respectively, and the coupling terms between the control (probe) field and each cavity. In the rotating frame with the frequency ω L i , the Hamiltonian is written as
= Δ 1 c 1 c 1 + Δ 2 c 2 c 2 + 1 2 ω m ( P 2 + Q 2 )   + J ( c 1 c 2 + c 2 c 1 ) + 2 g 1 c 1 c 1 Q 2 + 2 g 2 c 2 c 2 Q 2 + i ϵ L 1 ( c 1 c 1 ) + i ϵ L 2 ( c 2 c 2 ) + i ( ϵ P 1 c 1 e i δ 1 t ϵ P 1 * c 1 e i δ 1 t ) + i ( ϵ P 2 c 2 e i δ 2 t ϵ P 2 * c 2 e i δ 2 t ) ,
where Δ i = ω c i ω L i ( i = 1 ,   2 ) represents the detuning of the i t h cavity mode relative to the i t h control field and δ i = ω P i ω L i ( i = 1 ,   2 ) represents the detuning of the i t h probe field relative to the i t h control field. Using the Heisenberg–Langevin equations, the dynamics of the system [12,13,14,15,16] represented by Equation (2) can be described by the following set of equations:
c 1 ˙ = [ k 1 2 + i ( Δ 1 + 2 g 1 Q 2 ) ] c 1 i J c 2 + ϵ L 1 + ϵ P 1 e i δ 1 t , c 2 ˙ = [ k 2 2 + i ( Δ 2 + 2 g 2 Q 2 ) ] c 2 i J c 1 + ϵ L 2 + ϵ P 2 e i δ 2 t , Q ˙ = ω m P , P ˙ = ω m Q 4 ( g 1 c 1 c 1 + g 2 c 2 c 2 ) Q γ m P .
here, k i ( i = 1 ,   2 ) is the decay rate of the i t h cavity mode and γ m is the damping rate of the mechanical oscillator. When the intensity of the control field is far greater than that of the probe field, the above operators can be written as the sum of the steady-state value and an additional fluctuation value, that is, o = o s + δ o , where o represents the following operators: ( c 1 , c 2 , P , Q ) . Substituting this expression into Equation (3), while ignoring the nonlinear terms, we obtain the linearized Heisenberg–Langevin equations and obtain the steady-state values of the system. The steady-state values can be written as
c 1 s = i J ϵ L 2 + [ k 2 2 + i ( Δ 2 + 2 g 2 Q s 2 ) ] ϵ L 1 [ k 2 2 + i ( Δ 2 + 2 g 2 Q s 2 ) ] [ k 1 2 + i ( Δ 1 + 2 g 1 Q s 2 ) ] J 2 , c 2 s = i J ϵ L 1 + [ k 1 2 + i ( Δ 1 + 2 g 1 Q s 2 ) ] ϵ L 2 [ k 2 2 + i ( Δ 2 + 2 g 2 Q s 2 ) ] [ k 1 2 + i ( Δ 1 + 2 g 1 Q s 2 ) ] J 2 , P s = 0 ,     Q s = 0 .
It can be seen from Equation (4) that for quadratic coupling, the steady-state value Q s of the position is 0, but the potential energy 1 2 m ω m Q s 2 of the mechanical oscillator is not zero due to the thermal environment and zero-point fluctuations. Therefore, we define X = Q 2 , Y = P 2 , Z = Q P + P Q and, further, rewrite Equation (3) as follows
c 1 ˙ = [ k 1 2 + i ( Δ 1 + 2 g 1 X ) ] c 1 i J c 2 + ϵ L 1 + ϵ P 1 e i δ 1 t ,   c 2 ˙ = [ k 2 2 + i ( Δ 2 + 2 g 2 X ) ] c 2 i J c 1 + ϵ L 2 + ϵ P 2 e i δ 2 t , X ˙ = ω m Z , Y ˙ = [ ω m + 4 ( g 1 c 1 c 1 + g 2 c 2 c 2 ) ] Z 2 γ m Y + γ m ( 1 + 2 n t h ) ,   Z ˙ = 2 ω m Y 2 [ ω m + 4 ( g 1 c 1 c 1 + g 2 c 2 c 2 ) ] X γ m Z ,  
where γ m ( 1 + 2 n t h ) is the constant introduced by the coupling between the mechanical oscillator and the thermal environment with n t h = ( e ω m k B T 1 ) 1 being the mean thermal phonon occupation number of the mechanical oscillator at temperature T . k B is the Boltzmann constant. Substituting o = o s + δ o ( o = c 1 , c 2 , X , Y , Z ) into Equation (5), the steady-state values of the system can be expressed as
c 1 s = i J ϵ L 2 + ( κ 2 2 + i Δ 2 ) ϵ L 1 ( κ 2 2 + i Δ 2 ) ( κ 1 2 + i Δ 1 ) J 2 , c 2 s = i J ϵ L 1 + ( κ 1 2 + i Δ 1 ) ϵ L 2 ( κ 2 2 + i Δ 2 ) ( κ 1 2 + i Δ 1 ) J 2 , Z s = 0 , X s = ω m Y s ( ω m + 4 g 1 c 1 s c 1 s * + 4 g 2 c 2 s c 2 s * ) , Y s = 1 + 2 n t h 2 ,
where Δ i = Δ i + 2 g i X s ( i = 1 ,   2 ) is the effective detuning between the i t h cavity mode and the i t h control field. The quantum Langevin equations considering only the linear term of operator fluctuation can be written as:
δ c 1 ˙ = ( κ 1 2 + i Δ 1 ) δ c 1 2 i g 1 c 1 s δ X i J δ c 2 + ϵ P 1 e i δ 1 t , δ c 2 ˙ = ( κ 2 2 + i Δ 2 ) δ c 2 2 i g 2 c 2 s δ X i J δ c 1 + ϵ P 2 e i δ 2 t , δ X ˙ = ω m δ Z , δ Y ˙ = [ ω m + 4 ( g 1 c 1 s c 1 s * + g 2 c 2 s c 2 s * ) ] δ Z 2 γ m δ Y ,   δ Z ˙ = 2 ω m δ Y 8 X s [ ( g 1 c 1 s δ c 1 + g 1 c 1 s * δ c 1 ) + ( g 2 c 2 s * δ c 2 + g 2 c 2 s δ c 2 ) ] 2 [ ω m + 4 ( g 1 c 1 s c 1 s * + g 2 c 2 s c 2 s * ) ] δ X γ m δ Z .
The system is not necessarily stable in some parameters, so it is necessary to judge the stability of the system under different conditions by solving the Rouse–Hurwitz stability criterion. The vector V = ( δ R e [ c 1 ] , δ I m [ c 1 ] , δ R e [ c 2 ] , δ I m [ c 2 ] , δ X , δ Y , δ Z ) T and V i n = ( ϵ P 1 cos [ δ 1 t ] , ϵ P 1 sin [ δ 1 t ] , ϵ P 2 cos [ δ 2 t ] , ϵ P 2 sin [ δ 2 t ] , 0 , 0 , 0 ) T are defined, and Equation (7) is expressed as the following matrix form:
V ˙ = R V + V i n .
The coefficient matrix R is given by
κ 1 2 Δ 1 0 J 0 0 0 Δ 1 κ 1 2 J 0 0 2 G 1 0 0 J κ 2 2 Δ 2 0 2 G 2 sin [ θ ] 0 J 0 Δ 2 κ 2 2 0 2 G 2 cos [ θ ] 0 0 0 0 0 2 γ m 0 ω m , eff 0 0 0 0 0 0 ω m 16 G 1 X s 0 16 G 2 X s cos [ θ ] 16 G 2 X s sin [ θ ] 2 ω m 2 ω m , eff γ m ,
where G 1 = g 1 | c 1 s | = g 1 c 1 s e i θ 1 and G 2 = g 2 | c 2 s | = g 2 c 2 s e i θ 2 are the effective optomechanical coupling constants with the effective phase difference θ = θ 2 θ 1 , and ω m , e f f = ω m + 8 ( G 1 2 / g 1 ) + 8 ( G 2 2 / g 2 ) is the effective mechanical oscillator frequency. From Equation (6), we can change the values of c 1 s and c 2 s by adjusting the amplitudes and phase of the control field, and then change the values of G 1 , G 2 and θ . Next, for simplicity, we make the value of c 1 s a real number through appropriately choosing the amplitudes and phase of the control field, that is, θ 1 = 0 . At this time, the effective phase difference θ = θ 2 .
Supposing the solution of Equation (7) has the form o = o + e i δ t + o e i δ t , we set Δ 1 = Δ 2 = Δ , δ 1 = δ 2 = δ . By comparing the coefficients of i δ t and i δ t , the solution of δ c 1 + and δ c 2 + can be given by
δ c 1 + = k 2 ϵ p 2 [ ( i J ) ( A 1 M 1 + H 1 ) A 1 e i θ ] + k 1 ϵ p 1 [ χ 2 ( H 1 + A 1 M 1 ) + A 1 G 1 G 2 ] ( A 1 M 1 + H 1 ) A 2 A 1 H 2 ,   δ c 2 + = k 1 ϵ p 1 [ ( i J ) ( A 1 M 1 + H 1 ) A 1 e i θ ] + k 2 ϵ p 2 [ χ 1 ( H 1 + A 1 M 1 ) + A 1 G 1 G 2 ] ( A 1 M 1 + H 1 ) A 2 A 1 H 2 ,
with
χ 1 = ( k 1 2 + i Δ 1 ) + i δ , χ 11 = ( k 1 2 i Δ 1 ) + i δ ,   χ 2 = ( k 2 2 + i Δ 2 ) + i δ , χ 22 = ( k 2 2 i Δ 2 ) + i δ ,   A 1 = χ 11 χ 22 + J 2 , A 2 = χ 1 χ 2 + J 2 ,   H 1 = χ 11 G 2 G 1 + χ 22 G 1 G 2 2 i J c o s [ θ ] , H 2 = χ 1 G 1 G 2 + χ 2 G 2 G 1 + 2 i J c o s [ θ ] ,   M 1 = ( 2 γ m + i δ ) ( γ m + i δ ) ( i δ ) + 2 ( γ m + i δ ) ω m ω m , e f f 16 i ( 2 γ m + i δ ) G 1 G 2 ω m X s .  
The input–output relation of the optomechanical system can be written as
ϵ p 1 o u t + ϵ p 1 i n e i δ 1 t = κ 1 δ c 1 , ϵ p 2 o u t + ϵ p 2 i n e i δ 2 t = κ 2 δ c 2 .  
Here, ϵ p 1 i n = ϵ p 1 κ 1 , ϵ p 2 i n = ϵ p 2 κ 2   . Next, assuming that the output field also has the form of o = o + e i δ t + o e i δ t , by comparing the coefficients of i δ t and i δ t , the relations are obtained as follows:
ϵ p 1 + o u t = κ 1 δ c 1 + ϵ p 1 κ 1 , ϵ p 1 o u t = 0 , ϵ p 2 + o u t = κ 2 δ c 2 + ϵ p 2 κ 2 , ϵ p 2 o u t = 0 .
Transmission amplitudes can be expressed as [27]:
T 12 = | ϵ p 2 + o u t ϵ p 1 i n | ϵ p 2 i n = 0 , T 21 = | ϵ p 1 + o u t ϵ p 2 i n | ϵ p 1 i n = 0 .  
So, by substituting Equation (10) into Equations (12) and (13), transmission amplitudes can be obtained as
T 12 = | k 1 k 2   [ ( i J ) ( A 1 M 1 + H 1 ) A 1 e i θ ] ( A 1 M 1 + H 1 ) A 2 A 1 H 2 | , T 21 = | k 1 k 2   [ ( i J ) ( A 1 M 1 + H 1 ) A 1 e i θ ] ( A 1 M 1 + H 1 ) A 2 A 1 H 2 | .

3. Results and Discussions about Nonreciprocity

From Equation (14), we can find that transmission amplitudes ( T 12 and T 21 ) are related to the effective phase difference θ . The reason why the system produces a nonreciprocal response is that different transmission paths interfere with each other. For example, the input signal from c 1 can be transmitted to c 2 through two paths: path 1 ( c 1 b c 2   ) and path 2 ( c 1 c 2 ). In path 2,   G 1 ,   G 2 and J play the major role in this process. In the following discussions, to study the influence of effective phase difference θ on nonreciprocity, we fix a series of parameters (   k 1 = k 2 = k , Δ = 2 ω m , n t h = 0.5 , k = 2 γ m , ω m = 10 γ m ,   ω m , e f f = 1.01 ω m , J = 1 γ m , G 1 = G 2 = 1 2 γ m ) and plot the variation of transmission amplitudes ( T 12 and T 21 ) with standardized detuning x γ m (   x = δ 2 ω m ) in different effective phase difference θ . As shown in Figure 2, when θ = n π , transmission amplitudes T 12 = T 21 , and the system is reciprocal. When θ n π , transmission amplitudes T 12 T 21 , and the system is nonreciprocal. When θ = π 2 and θ = 3 π 2 , the transmission direction of the system is reserved. In the case θ = π 2 + n π , the spectral line is symmetrical about the axis x γ m = 0 . In the range of 0 < θ < π , transmission amplitudes T 12 > T 21 , and in the case of π < θ < 2 π , transmission amplitudes T 12 < T 21 . Therefore, the phase difference will not only affect the magnitude of nonreciprocal response, but also affect the direction of nonreciprocal response.
Next, the perfect nonreciprocity under the condition with red detuning ( Δ i = 2 ω m ,   x = δ 2 ω m ) or blue detuning ( Δ i = 2 ω m , x = δ + 2 ω m ) is discussed, respectively. The perfect nonreciprocity means that we can realize the transmission process with T 12 = 1 , T 21 = 0 or with T 12 = 0 , T 21 = 1 . In the following discussion, only one of the above two processes will be considered. At this time, the signal can be completely transmitted from c 1 to c 2 , but it cannot be completely transmitted from c 2 to c 1 , which can realize unidirectional transmission of the signal, namely
T 12 = 1 , T 21 = 0 .
According to Equation (15), we obtain
( i J ) ( A 1 M 1 + H 1 ) = A 1 e i θ , | 2 J k 1 k 2 s i n [ θ ] A 2 e i θ + i J H 2 | = 1 .
When θ = π 2 and the system is in red detuning ( Δ = 2 ω m ) , by solving the above equation, we obtain
J = 1 2 k 1 k 2   ,     G 2 2 G 1 2 = k 2 k 1 ( x = 0 )   G 1 = k 1 ω m , e f f ( γ m ( k 1 k 2 32 ω m 2 ) U 1 + 4 ( k 1 + k 2 ) ω m U 2 ) 2 2 2 ( 1 + 2 n t h ) ( k 1 k 2 16 ω m 2 ) ω m ( γ m 2 + ω m 2 ) = k 1 ω m , e f f ( 4 ( k 1 + k 2 ) γ m ω m U 1 ( k 1 k 2 32 ω m 2 ) U 2 ) 2 2 8 ( 1 + 2 n t h ) ( k 1 + k 2 ) ω m 2 ( γ m 2 + ω m 2 ) ,
With U 1 = γ m 2 + ω m ( ω m , e f f + ω m ) ,   U 2 = 2 ω m ( γ m 2 + ω m 2 ) + ω m , e f f ( γ m 2 + 2 ω m 2 ) . It can be obtained that the perfect nonreciprocity appears at the resonance frequency ( x = 0 ) . The coupling constant J satisfying the perfect nonreciprocal condition is related to the decay rates k 1 and k 2 of the cavity, and satisfies the condition J = 1 2 k 1 k 2 . At the same time, it is obvious that the value of G i ( i = 1 ,   2 ) is related to a series of parameters (   k i , γ m , ω m , ω e f f , n t h ), which must satisfy both expressions of G i in Equation (17) simultaneously. As Figure 3 shows, when k 1 = k and k 2 = 10 k , with the increase in ω m , the ratio of the effective mechanical oscillator frequency ω m , e f f to the mechanical oscillator frequency ω m gradually approaches 1.
Assuming that ω m is much larger than other parameters, the effective mechanical oscillator frequency ω m , e f f ω m is approximately considered. For the value of   G i ( i = 1 , 2 ) , only the highest power of ω m and its previous coefficients are retained; thus, J = 1 2 k 1 k 2 , G 1 k 1 γ m 2 1 + 2 n t h and G 2 k 2 γ m 2 1 + 2 n t h are obtained. Substituting the above conditions into Equation (14) shows that although G i satisfying perfect nonreciprocity is related to the mean thermal phonon occupation number n t h , the value of n t h has no influence on the final transmission amplitude under this condition, so the value of n t h is arbitrary. In the case of red detuning, compared with the condition ( J = 1 2 k 1 k 2 , G i = k i γ m 2 ) that the double optomechanical systems with linear coupling should achieve perfect nonreciprocity under when θ = π 2 , the value of G i of the double optomechanical systems with quadratic coupling is related to n t h , and the value of n t h can be controlled by temperature T , so the G i can achieve perfect nonreciprocity at different values by controlling temperature T . Under the above conditions, the relationship between transmission amplitudes ( T 12 and T 21 ) and normalized detuning x k at different γ m k is shown in Figure 4, in which the other parameters are taken as k 1 = k , k 2 = 10 k , ω m = 40 k , ω m , e f f = ω m and the values of γ m k are 0.01, 0.1, 0.5 and 2, in turn. Under the conditions that k 1 and k 2 are fixed, and J is also fixed, as γ m takes different values, the values of G 1 and G 2 will change, but G 2 = 10 G 1 is always true. It can be seen from Figure 4 that when γ m k is relatively small, the spectral line width of transmission amplitude T 21 is proportional to the damping rate of the mechanical oscillator ( ω m γ m ) and perfect nonreciprocity can be achieved even when γ m is very small. With the increase in γ m , it leads to the increase in   G i , and the effective optomechanical coupling constant   G i begins to play a major role in the signal transmission process, and the range of nonreciprocity increases accordingly.
Next, we consider that the system is in the case with θ = π 2 and blue detuning ( Δ i = 2 ω m ) . When ω m is much larger than other parameters, Equation (16) can be solved in a similar way to the previous case of the red detuning, and only the highest power of ω m and the previous coefficient can be retained for the values of G i and J , so that J ± = ( 3 ± 2 2 ) k 1 k 2 2 , G 1 ± 3 ± 2 2 k 1 γ m 2 1 + 2 n t h and G 2 ± 3 ± 2 2 k 2 γ m 2 1 + 2 n t h can be obtained. Perfect nonreciprocity is also obtained at resonance frequency ( x = 0 ) . In the case of blue detuning, G i is related to the mean thermal phonon occupation number n t h , which is same as in the case of red detuning, but the value of n t h has no effect on the final transmission amplitudes, so G i can realize the perfect nonreciprocity at different values by controlling the temperature T . By solving the Rouse–Hurwitz stability criterion, we obtain that when J = J and G = G   , the system is always stable. When J = J + and G = G + , the system is stable only when it meets the conditions shown in Figure 5.
In Figure 6a–d, the values of the parameters are the same as those shown in Figure 5, and we plot the relationship between transmission amplitudes ( T 12 and T 21 ) and normalized detuning x k at different γ m k . The values of γ m k are taken as 0.01, 0.1, 0.8 and 1, in turn. When J = J + and G = G + , it can be seen from Figure 5 that the system is unstable in the range of γ m > 1.16 k . From Figure 6c,d, where γ m k = 0.8 and 1, the system is close to the critical state of stability and instability. In this case, the system displays a gain phenomenon for the probe field. The occurrence of gain is due to nonlinear effects based on the quantum coherence, which means optomechanically induced amplification [31] takes place. At this moment, the peak value of transmission amplitude is far greater than 1. In addition, the spectral line of T 12 is multi-peak. As shown in Figure 6a, when γ m k is small, the spectrum width of transmission amplitude is proportional to the attenuation rate γ m of the mechanical oscillator. Excluding the extremely small range around x = 0 , transmission amplitude T 12 = T 21 , which means the processes of transmission are reciprocal.
Finally, we make J = J , G = G and let γ m k equal 0.01, 0.1, 1 and 2, in turn. The other parameters are the same as those in Figure 6. As shown in Figure 7, with the increase in   γ m k , the maximum value of T 21 gradually decreases and the spectral line width of T 12 gradually increases, which means that the range of nonreciprocal response becomes larger. However, due to the system being stable in the whole range, there will be no transmission amplitudes greater than 1, and multi-peaked spectral lines and perfect nonreciprocity can be achieved at the resonance frequency ( x = 0 ) .

4. Conclusions

In sum, we study the optical nonreciprocity of double optomechanical systems with quadratic coupling and find that the magnitude and direction of the nonreciprocity of the system can be changed by adjusting the nonreciprocal phase difference θ . At the same time, it is calculated that when the system is in red detuning with θ = π 2 and ω m is much larger than other parameters, there is a group of G i and J that meet the conditions of perfect nonreciprocity. When the system is in blue detuning with θ = π 2 and ω m is much larger than other parameters, there are two groups of G i and J that can achieve perfect nonreciprocity. Specifically, the system is not always stable under some parameters, one of which causes the maximum value of transmission amplitudes to be far greater than 1 in a certain range. In the case of θ = ± π 2 and Δ i = ± 2 ω m , the perfect nonreciprocity is obtained at the resonance frequency ( x = 0 ) , whether the system is in red or blue detuning. When other parameters are fixed, the ratio of γ m k can affect the magnitude and range of nonreciprocity. In addition, we compare the conditions of realizing perfect nonreciprocity between the double optomechanical systems with linear coupling and the double optomechanical systems with quadratic coupling and find that when the decay rate k i of the cavity mode and the damping rate γ m of the mechanical oscillator are fixed, both of them satisfy the same relation for the value of J . However, the value of effective optomechanical coupling constant G i is different and is not only related to k i and γ m but also related to n t h , where the value of n t h will not affect the transmission amplitudes. Therefore, compared with the double optomechanical systems with linear coupling, the double optomechanical systems with quadratic coupling can change the value of G i and achieve perfect nonreciprocity by keeping the system at different temperatures T . In addition, we discuss the three-mode optomechanical system where the coupling between the cavity field and the mechanical oscillator is in the quadratic case. Comparing with the Ref. [27] in which optical nonreciprocity in three-mode optomechanical systems with linear coupling is investigated, we extend their model to the case with quadratic coupling. The Hamiltonian of the linear coupling is proportional to c 1 c 1 Q , which means the linear coupling denotes work done by light pressure pushing the oscillator to displacement Q . The Hamiltonian of the quadratic coupling, by contrast, is proportional to c 1 c 1 Q 2 , which means the quadratic coupling will change the native frequency of the oscillator. In order to achieve optical nonreciprocity, these two different couplings have different requirements for selections of physical parameters. For the linear coupling, the effective detuning needs to be equal to native frequency ω m of the oscillator [27]. However, for the quadratic coupling in our model, the condition that the effective detuning is equal to ω m will be broken. Here, the realization of optical nonreciprocity requires the effective detuning with 2 ω m . We know that production of these two different kinds of coupling depends on the position of the oscillator. Our studies show that the three-mode optomechanical systems with quadratic coupling can also be used as an alternative scheme for producing optical nonreciprocity.

Author Contributions

Conceptualization and model, S.H. and A.C.; numerical simulation, X.J. and P.P.; writing—original draft preparation, X.J.; writing—review and editing, S.H. and A.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Natural Science Foundation of China through Grants No. 12175199, and No. 12174344, Zhejiang Provincial Natural Science Foundation of China through Grants No. LZ20A040002 and No. LY21A040007, and the Science Foundation of Zhejiang Sci-Tech University through Grants No. 19062408-Y and No. 17062071-Y.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Potton, R.J. Reciprocity in optics. Rep. Prog. Phys. 2004, 67, 717. [Google Scholar] [CrossRef]
  2. Shadrivov, I.V.; Fedotov, V.A.; Powell, D.A.; Kivshar, Y.S.; Zheludev, N.I. Electromagnetic wave analog of an electronic diode. New J. Phys. 2011, 13, 033025. [Google Scholar] [CrossRef]
  3. Zaman, T.R.; Guo, X.; Ram, R.J. Faraday rotation in an InP waveguide. Appl. Phys. Lett. 2007, 90, 023514. [Google Scholar] [CrossRef]
  4. Hadad, Y.; Steinberg, B. Magnetized Spiral Chains of Plasmonic Ellipsoids for One-Way Optical Waveguides. Phys. Rev. Lett. 2010, 105, 233904. [Google Scholar] [CrossRef] [Green Version]
  5. Fan, L.; Varghese, L.T.; Wang, J.; Xuan, Y.; Weiner, A.M.; Qi, M. Silicon optical diode with 40 dB nonreciprocal transmission. Opt. Lett. 2013, 38, 1259–1261. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Lenferink, E.J.; Wei, G.; Stern, N.P. Coherent optical non-reciprocity in axisymmetric resonators. Opt. Express 2014, 22, 16099–16111. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Jalas, D.; Petrov, A.; Eich, M.; Freude, W.; Fan, S.; Yu, Z.; Baets, R.; Popović, M.; Melloni, A.; Joannopoulos, J.D.; et al. What is and what is not: An optical isolator. Nat. Photon. 2013, 7, 579–582. [Google Scholar] [CrossRef]
  8. Shi, Y.; Yu, Z.; Fan, S. Limitations of nonlinear optical isolators due to dynamic reciprocity. Nat. Photon. 2015, 9, 388–392. [Google Scholar] [CrossRef]
  9. Wang, J.; Fan, L.; Varghese, L.T.; Shen, H.; Xuan, Y.; Niu, B.; Qi, M. A theoretical model for an optical diode built with nonlinear silicon microrings. J. Light. Tech. 2013, 31, 313–321. [Google Scholar] [CrossRef]
  10. Fan, L.; Wang, J.; Varghese, L.T.; Shen, H.; Niu, B.; Xuan, Y.; Weiner, A.M.; Qi, M. An All-Silicon Passive Optical Diode. Science 2012, 335, 447–450. [Google Scholar] [CrossRef] [PubMed]
  11. Huang, R.; Miranowicz, A.; Liao, J.-Q.; Nori, F.; Jing, H. Nonreciprocal Photon Blockade. Phys. Rev. Lett. 2018, 121, 153601. [Google Scholar] [CrossRef] [Green Version]
  12. Marquardt, F.; Girvin, S.M. Trend: Optomechanics. Physics 2009, 2, 40. [Google Scholar] [CrossRef] [Green Version]
  13. Meystre, P. A short walk through quantum optomechanics. Ann. Phys. 2012, 525, 215–233. [Google Scholar] [CrossRef] [Green Version]
  14. Aspelmeyer, M.; Meystre, P.; Schwab, K. Quantum optomechanics. Phys. Today 2012, 65, 29–35. [Google Scholar] [CrossRef]
  15. Kippenberg, T.J.; Vahala, K.J. Cavity Optomechanics: Back-Action at the Mesoscale. Science 2008, 321, 1172–1176. [Google Scholar] [CrossRef] [Green Version]
  16. Aspelmeyer, M.; Kippenberg, T.J.; Marquardt, F. Cavity optomechanics. Rev. Mod. Phys. 2014, 86, 1391–1452. [Google Scholar] [CrossRef] [Green Version]
  17. Agarwal, G.S.; Huang, S. Electromagnetically induced transparency in mechanical effects of light. Phys. Rev. A 2010, 81, 041803. [Google Scholar] [CrossRef] [Green Version]
  18. Huang, S.; Agarwal, G.S. Electromagnetically induced transparency with quantized fields in optocavity mechanics. Phys. Rev. A 2011, 83, 043826. [Google Scholar] [CrossRef] [Green Version]
  19. Li, G.; Nie, W.; Li, X.; Chen, A. Dynamics of ground-state cooling and quantum entanglement in a modulated optomechanical system. Phys. Rev. A 2019, 100, 063805. [Google Scholar] [CrossRef]
  20. Huang, S.; Chen, A. Improving the cooling of a mechanical oscillator in a dissipative optomechanical system with an optical parametric amplifier. Phys. Rev. A 2018, 98, 063818. [Google Scholar] [CrossRef]
  21. Abdi, M.; Barzanjeh, S.; Tombesi, P.; Vitali, D. Effect of phase noise on the generation of stationary entanglement in cavity optomechanics. Phys. Rev. A 2011, 84, 032325. [Google Scholar] [CrossRef] [Green Version]
  22. Genes, C.; Vitali, D.; Tombesi, P. Emergence of atom-light-mirror entanglement inside an optical cavity. Phys. Rev. A 2008, 77, 050307. [Google Scholar] [CrossRef] [Green Version]
  23. Deng, H.; Zou, F.; Huang, J.-F.; Liao, J.-Q. Optical normal-mode-induced phonon-sideband splitting in the photon-blockade effect. Phys. Rev. A 2021, 104, 033706. [Google Scholar] [CrossRef]
  24. Wang, D.; Bai, C.H.; Xing, Y.; Liu, S.; Zhang, S.; Wang, H.F. Enhanced photon blockade via driving a trapped Λ-type atom in a hybrid optomechanical system. Phys. Rev. A 2020, 102, 043705. [Google Scholar] [CrossRef]
  25. Zheng, L.-L.; Yin, T.-S.; Bin, Q.; Lü, X.-Y.; Wu, Y. Single-photon-induced phonon blockade in a hybrid spin-optomechanical system. Phys. Rev. A 2019, 99, 013804. [Google Scholar] [CrossRef]
  26. Hafezi, M.; Rabl, P. Optomechanically induced non-reciprocity in microring resonators. Opt. Express 2012, 20, 7672–7684. [Google Scholar] [CrossRef]
  27. Xu, X.W.; Li, Y. Optical nonreciprocity and optomechanical circulator in three-mode optomechanical systems. Phys. Rev. A 2015, 91, 053854. [Google Scholar] [CrossRef] [Green Version]
  28. Xu, X.; Zhao, Y.; Wang, H.; Jing, H.; Chen, A. Quantum nonreciprocality in quadratic optomechanics. Photon. Res. 2020, 8, 143–150. [Google Scholar] [CrossRef]
  29. Li, G.; Xiao, X.; Li, Y.; Wang, X. Tunable optical nonreciprocity and a phonon-photon router in an optomechanical system with coupled mechanical and optical modes. Phys. Rev. A 2018, 97, 023801. [Google Scholar] [CrossRef]
  30. Li, B.; Özdemir, S.K.; Xu, X.W.; Zhang, L.; Kuang, L.M.; Jing, H. Nonreciprocal optical solitons in a spinning Kerr resonator. Phys. Rev. A 2021, 103, 053522. [Google Scholar] [CrossRef]
  31. Yan, X.-B.; Jia, W.Z.; Li, Y.; Wu, J.-H.; Li, X.-L.; Mu, H.-W. Optomechanically induced amplification and perfect transparency in double-cavity optomechanics. Front. Phys. 2015, 10, 351–357. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of the double optomechanical systems: the two optical cavities with optical modes c 1 and c 2 are coupled with a mechanical oscillator with the mechanical mode b . ( p and q are the momentum and position operators of the mechanical oscillator, respectively). The two cavities are coupled with a coupling constant of J , and each cavity is driven by a control field with amplitude ϵ L i ( i = 1 ,   2 ) and a probe field with amplitude ϵ P i ( i = 1 ,   2 ) simultaneously.
Figure 1. Schematic diagram of the double optomechanical systems: the two optical cavities with optical modes c 1 and c 2 are coupled with a mechanical oscillator with the mechanical mode b . ( p and q are the momentum and position operators of the mechanical oscillator, respectively). The two cavities are coupled with a coupling constant of J , and each cavity is driven by a control field with amplitude ϵ L i ( i = 1 ,   2 ) and a probe field with amplitude ϵ P i ( i = 1 ,   2 ) simultaneously.
Photonics 09 00728 g001
Figure 2. Transmission amplitudes T 12 (Red dashed lines) and T 21 (Black solid lines) vary with x γ m (   x = δ 2 ω m ) for different effective phase difference θ : (a) θ = 0 , (b) θ = π 4 , (c) θ = π 2 , (d) θ = 3 π 4 , (e) θ = π , (f) θ = 5 π 4 , (g) θ = 3 π 2 , and (h) θ = 7 π 4 . The other parameters are as follows: k 1 = k 2 = k ,     Δ = 2 ω m ,   n = 0.5 ,   ω m = 10 γ m ,   ω m , e f f = 1.01 ω m ,   J = 1 γ m ,   k = 2 γ m , and G 1 = G 2 = 1 2 γ m .
Figure 2. Transmission amplitudes T 12 (Red dashed lines) and T 21 (Black solid lines) vary with x γ m (   x = δ 2 ω m ) for different effective phase difference θ : (a) θ = 0 , (b) θ = π 4 , (c) θ = π 2 , (d) θ = 3 π 4 , (e) θ = π , (f) θ = 5 π 4 , (g) θ = 3 π 2 , and (h) θ = 7 π 4 . The other parameters are as follows: k 1 = k 2 = k ,     Δ = 2 ω m ,   n = 0.5 ,   ω m = 10 γ m ,   ω m , e f f = 1.01 ω m ,   J = 1 γ m ,   k = 2 γ m , and G 1 = G 2 = 1 2 γ m .
Photonics 09 00728 g002
Figure 3. When θ = π 2 and the system is in red detuning ( Δ = 2 ω m ) , to satisfy the perfect nonreciprocity, the ratio of the effective mechanical oscillator frequency to the mechanical oscillator frequency ω m , e f f ω m varies with ω m k at different γ m k ( 0.01 ,   0.1 ,   0.5 ,   2 ) . Here, γ m k = 0.01 (black), 0.1 (orange), 0.5 (blue) and 2 (red), respectively. The other parameters are k 1 = k , and k 2 = 10 k .
Figure 3. When θ = π 2 and the system is in red detuning ( Δ = 2 ω m ) , to satisfy the perfect nonreciprocity, the ratio of the effective mechanical oscillator frequency to the mechanical oscillator frequency ω m , e f f ω m varies with ω m k at different γ m k ( 0.01 ,   0.1 ,   0.5 ,   2 ) . Here, γ m k = 0.01 (black), 0.1 (orange), 0.5 (blue) and 2 (red), respectively. The other parameters are k 1 = k , and k 2 = 10 k .
Photonics 09 00728 g003
Figure 4. When θ = π 2 and the system is in red detuning ( Δ = 2 ω m ) , transmission amplitudes T 21 (Black solid lines) and T 12 (Red dashed lines) vary with the normalized detuning x k ( x = δ 2 ω m ) for different values of γ m k : (a) γ m k = 0.01 , (b) γ m k = 0.1 , (c) γ m k = 0.5 , and (d) γ m k = 2 . The other parameters are as follows: J = 1 2 k 1 k 2 , G 1 = k 1 γ m 2 1 + 2 n t h , G 2 = k 2 γ m 2 1 + 2 n t h , k 1 = k , k 2 = 10 k , ω m = 40 k , and ω m , e f f = ω m .
Figure 4. When θ = π 2 and the system is in red detuning ( Δ = 2 ω m ) , transmission amplitudes T 21 (Black solid lines) and T 12 (Red dashed lines) vary with the normalized detuning x k ( x = δ 2 ω m ) for different values of γ m k : (a) γ m k = 0.01 , (b) γ m k = 0.1 , (c) γ m k = 0.5 , and (d) γ m k = 2 . The other parameters are as follows: J = 1 2 k 1 k 2 , G 1 = k 1 γ m 2 1 + 2 n t h , G 2 = k 2 γ m 2 1 + 2 n t h , k 1 = k , k 2 = 10 k , ω m = 40 k , and ω m , e f f = ω m .
Photonics 09 00728 g004
Figure 5. A stability diagram of the system versus k and γ m . The used parameters are θ = π 2 , Δ i = 2 ω m ,   k 1 = k 2 = k , J = J + , G = G + , ω m = 40 k , and ω m , e f f = ω m .
Figure 5. A stability diagram of the system versus k and γ m . The used parameters are θ = π 2 , Δ i = 2 ω m ,   k 1 = k 2 = k , J = J + , G = G + , ω m = 40 k , and ω m , e f f = ω m .
Photonics 09 00728 g005
Figure 6. When θ = π 2 and the system is in blue detuning ( Δ = 2 ω m ) , transmission amplitudes T 21 (Black solid lines) and T 12 (Red dashed lines) vary with the normalized detuning x k   ( x = δ + 2 ω m ) for different γ m k . The values of γ m k are as follows: (a) γ m k = 0.01 , (b) γ m k = 0.1 , (c) γ m k = 0.8 , and (d) γ m k = 1 . All other parameters are the same as in Figure 5.
Figure 6. When θ = π 2 and the system is in blue detuning ( Δ = 2 ω m ) , transmission amplitudes T 21 (Black solid lines) and T 12 (Red dashed lines) vary with the normalized detuning x k   ( x = δ + 2 ω m ) for different γ m k . The values of γ m k are as follows: (a) γ m k = 0.01 , (b) γ m k = 0.1 , (c) γ m k = 0.8 , and (d) γ m k = 1 . All other parameters are the same as in Figure 5.
Photonics 09 00728 g006
Figure 7. Transmission amplitudes T 21 (Black solid lines) and T 12 (Red dashed lines) vary with the normalized detuning x k under the parameters of J = J and G = G . The values of γ m k ( x = δ + 2 ω m ) are as follows: (a) γ m k = 0.01 , (b)   γ m k = 0.1 , (c) γ m k = 1 , (d)   γ m k = 2 . All other parameters are the same as those in Figure 6.
Figure 7. Transmission amplitudes T 21 (Black solid lines) and T 12 (Red dashed lines) vary with the normalized detuning x k under the parameters of J = J and G = G . The values of γ m k ( x = δ + 2 ω m ) are as follows: (a) γ m k = 0.01 , (b)   γ m k = 0.1 , (c) γ m k = 1 , (d)   γ m k = 2 . All other parameters are the same as those in Figure 6.
Photonics 09 00728 g007
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ji, X.; Pan, P.; Huang, S.; Chen, A. Optical Nonreciprocity in Double Optomechanical Systems with Quadratic Coupling. Photonics 2022, 9, 728. https://doi.org/10.3390/photonics9100728

AMA Style

Ji X, Pan P, Huang S, Chen A. Optical Nonreciprocity in Double Optomechanical Systems with Quadratic Coupling. Photonics. 2022; 9(10):728. https://doi.org/10.3390/photonics9100728

Chicago/Turabian Style

Ji, Xinyu, Peipei Pan, Sumei Huang, and Aixi Chen. 2022. "Optical Nonreciprocity in Double Optomechanical Systems with Quadratic Coupling" Photonics 9, no. 10: 728. https://doi.org/10.3390/photonics9100728

APA Style

Ji, X., Pan, P., Huang, S., & Chen, A. (2022). Optical Nonreciprocity in Double Optomechanical Systems with Quadratic Coupling. Photonics, 9(10), 728. https://doi.org/10.3390/photonics9100728

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop