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
represents the annihilation operator (frequency) of the
cavity and
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
.
and
represent the momentum operator and position operator of the mechanical resonator, respectively, and define the operator
under the condition of canonical quantization. Here,
is the Planck constant,
is the zero-point fluctuation of the mechanical resonator displacement, and
P and
satisfy the relation of
. In addition, the two cavities are coupled with a linear coupling constant
, and each cavity field is driven by a control field
and a probe field
simultaneously. The model can be implemented in the membrane-in-the-middle optomechanical system. The Hamiltonian can be written as
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 , the Hamiltonian is written as
where
represents the detuning of the
cavity mode relative to the
control field and
represents the detuning of the
probe field relative to the
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:
here,
is the decay rate of the
cavity mode and
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,
, where
represents the following operators:
. 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
It can be seen from Equation (4) that for quadratic coupling, the steady-state value of the position is 0, but the potential energy of the mechanical oscillator is not zero due to the thermal environment and zero-point fluctuations. Therefore, we define and, further, rewrite Equation (3) as follows
where
is the constant introduced by the coupling between the mechanical oscillator and the thermal environment with
being the mean thermal phonon occupation number of the mechanical oscillator at temperature
.
is the Boltzmann constant. Substituting
into Equation (5), the steady-state values of the system can be expressed as
where
is the effective detuning between the
cavity mode and the
control field. The quantum Langevin equations considering only the linear term of operator fluctuation can be written as:
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
and
are defined, and Equation (7) is expressed as the following matrix form:
The coefficient matrix is given by
where
and
are the effective optomechanical coupling constants with the effective phase difference
, and
is the effective mechanical oscillator frequency. From Equation (6), we can change the values of
and
by adjusting the amplitudes and phase of the control field, and then change the values of
and
. Next, for simplicity, we make the value of
a real number through appropriately choosing the amplitudes and phase of the control field, that is,
. At this time, the effective phase difference
.
Supposing the solution of Equation (7) has the form , we set . By comparing the coefficients of and , the solution of and can be given by
The input–output relation of the optomechanical system can be written as
Here,
. Next, assuming that the output field also has the form of
, by comparing the coefficients of
and
, the relations are obtained as follows:
Transmission amplitudes can be expressed as [
27]:
So, by substituting Equation (10) into Equations (12) and (13), transmission amplitudes can be obtained as
3. Results and Discussions about Nonreciprocity
From Equation (14), we can find that transmission amplitudes (
and
) 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
can be transmitted to
through two paths: path 1 (
) and path 2
). In path 2,
,
and
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
and plot the variation of transmission amplitudes (
and
) with standardized detuning
in different effective phase difference
. As shown in
Figure 2, when
, transmission amplitudes
, and the system is reciprocal. When
, transmission amplitudes
, and the system is nonreciprocal. When
and
, the transmission direction of the system is reserved. In the case
, the spectral line is symmetrical about the axis
. In the range of
, transmission amplitudes
, and in the case of
, transmission amplitudes
. 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
or blue detuning
is discussed, respectively. The perfect nonreciprocity means that we can realize the transmission process with
or with
. In the following discussion, only one of the above two processes will be considered. At this time, the signal can be completely transmitted from
to
, but it cannot be completely transmitted from
to
, which can realize unidirectional transmission of the signal, namely
According to Equation (15), we obtain
When and the system is in red detuning , by solving the above equation, we obtain
With
,
. It can be obtained that the perfect nonreciprocity appears at the resonance frequency
. The coupling constant
satisfying the perfect nonreciprocal condition is related to the decay rates
and
of the cavity, and satisfies the condition
. At the same time, it is obvious that the value of
is related to a series of parameters (
), which must satisfy both expressions of
in Equation (17) simultaneously. As
Figure 3 shows, when
and
, with the increase in
, the ratio of the effective mechanical oscillator frequency
to the mechanical oscillator frequency
gradually approaches 1.
Assuming that
is much larger than other parameters, the effective mechanical oscillator frequency
is approximately considered. For the value of
, only the highest power of
and its previous coefficients are retained; thus,
and
are obtained. Substituting the above conditions into Equation (14) shows that although
satisfying perfect nonreciprocity is related to the mean thermal phonon occupation number
, the value of
has no influence on the final transmission amplitude under this condition, so the value of
is arbitrary. In the case of red detuning, compared with the condition
that the double optomechanical systems with linear coupling should achieve perfect nonreciprocity under when
, the value of
of the double optomechanical systems with quadratic coupling is related to
, and the value of
can be controlled by temperature
, so the
can achieve perfect nonreciprocity at different values by controlling temperature
. Under the above conditions, the relationship between transmission amplitudes (
and
) and normalized detuning
at different
is shown in
Figure 4, in which the other parameters are taken as
and the values of
are 0.01, 0.1, 0.5 and 2, in turn. Under the conditions that
and
are fixed, and
is also fixed, as
takes different values, the values of
and
will change, but
is always true. It can be seen from
Figure 4 that when
is relatively small, the spectral line width of transmission amplitude
is proportional to the damping rate of the mechanical oscillator
and perfect nonreciprocity can be achieved even when
is very small. With the increase in
, it leads to the increase in
, and the effective optomechanical coupling constant
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
and blue detuning
. When
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
and the previous coefficient can be retained for the values of
and
, so that
and
can be obtained. Perfect nonreciprocity is also obtained at resonance frequency
. In the case of blue detuning,
is related to the mean thermal phonon occupation number
, which is same as in the case of red detuning, but the value of
has no effect on the final transmission amplitudes, so
can realize the perfect nonreciprocity at different values by controlling the temperature
. By solving the Rouse–Hurwitz stability criterion, we obtain that when
and
, the system is always stable. When
and
, 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 (
and
) and normalized detuning
at different
. The values of
are taken as 0.01, 0.1, 0.8 and 1, in turn. When
and
, it can be seen from
Figure 5 that the system is unstable in the range of
. From
Figure 6c,d, where
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
is multi-peak. As shown in
Figure 6a, when
is small, the spectrum width of transmission amplitude is proportional to the attenuation rate
of the mechanical oscillator. Excluding the extremely small range around
, transmission amplitude
, which means the processes of transmission are reciprocal.
Finally, we make
and let
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
, the maximum value of
gradually decreases and the spectral line width of
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
.
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
and
is much larger than other parameters, there is a group of
and
that meet the conditions of perfect nonreciprocity. When the system is in blue detuning with
and
is much larger than other parameters, there are two groups of
and
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
and
, the perfect nonreciprocity is obtained at the resonance frequency
, whether the system is in red or blue detuning. When other parameters are fixed, the ratio of
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
of the cavity mode and the damping rate
of the mechanical oscillator are fixed, both of them satisfy the same relation for the value of
. However, the value of effective optomechanical coupling constant
is different and is not only related to
and
but also related to
, where the value of
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
and achieve perfect nonreciprocity by keeping the system at different temperatures
. 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
, which means the linear coupling denotes work done by light pressure pushing the oscillator to displacement
. The Hamiltonian of the quadratic coupling, by contrast, is proportional to
, 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
of the oscillator [
27]. However, for the quadratic coupling in our model, the condition that the effective detuning is equal to
will be broken. Here, the realization of optical nonreciprocity requires the effective detuning with
. 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.