Flip-flop Converter of Dual-bistability Using Cavity and Parametric Amplified Four-Wave Mixing

We study the realization of dual-bistability flip-flop converter in cavity and parametrically amplified four-wave mixing (FWM) process at a four-level cavity atomic system. Using the effect of nonreciprocity optical dual-bistability, we can obtain different output multi-mode states of probe transmission signal and FWM signal. We find the channel equalization ratio and optical contrast between multi-mode states is related to the degree of dual-bistability. Besides, the degree of dual-bistability can be controlled by the input parameters (frequency detuning and powers of the dressing beams). More, using electro-optical modulator and acoustic optical modulator to modulate the powers and frequency detuning, respectively, we can realize the fast conversion between different output states. And the switch speed of this flip-flop converter is about 16 ns. These outcomes may provide foundation for the development of all-optical devices and quantum information processing.

In the last few years, a great deal of optical phenomena based on atomic coherence and quantum interference have attracted intensive attention of many researchers in the multilevel atomic systems. Four-wave mixing (FWM) processes due to atomic coherence have been experimentally reported in all alkali atoms [1][2][3][4][5][6] and various experimental schemes [7][8][9] . With widely studied in applications, FWM plays potential roles in the sub shot noise measurements 1,10,11 , quantum imaging 12,13 , quantum communication 14 , slow light 15 , storage of light 16 and relative intensity squeezing 5 . Besides, with the self-stabilizing function, FWM can cause the energy transfer of different waves [17][18][19] . A spontaneous parametric four wave mixing (SP-FWM) process that generates two weak fields can amplify the seeded signal, the phenomenon is called optical parametrically amplification (OPA) 20,21 . More, optical bistability (OB) behavior based on atomic coherence and quantum interference are practiced in recent decades 22,23 . On the other hand, some schemes for realizing optical stability through multi-wave mixing process in an optical cavity have been studied experimentally and theoretically 15 , where cavity can provide a feedback as an essential factor for the generation of bistability 24 . OB has also been demonstrated without a cavity using degenerate FWM in atomic vapor with two counter-propagating laser beams 25,26 . The OB in multilevel atoms inside optical cavities 27,28 has been the subject of many recent studies because of its broad application prospects in all-optical logic and memories performance 29,30 . Bistability and instability were also observed in cold clouds of cesium atoms inside an optical cavity, where degenerate Zeeman sublevels participate in the dynamic processes 31 .
In this paper, we investigate the nonreciprocity optical dual-bistability (ODB) phenomena of probe transmission signal (PTS) and four-wave mixing (FWM) signal in a composite atom-cavity system both theoretically and experimentally. Based on the relationship of vacuum Ribi splitting and OB 24 , the ODB is obtained by scanning the frequency detuning of beams, which is more sensitive than scanning the powers. The nonreciprocity of frequency (frequency offset, in x direction) and intensity (shape change, in y direction) between the two signals (from the rising and falling edges in one frequency scanning round trip) are attributed to "∞"-shape non-overlapping region 32 , which can reflect the degree of dual-bistability directly. Besides, the frequency offset is caused by the cavity feedback dressing effect and has been amplified by the effect of vacuum induced enhancement in cavity 24 . Moreover, the nonlinear refractive index and the enhancement and suppression of FWM cavity mode have been studied. By changing the frequency detuning and the powers of dressing beam, the feedback intensity will ScIEntIfIc REPORts | (2018) 8:2492 | DOI:10.1038/s41598-018-20962-5 be enhanced or suppressed, which results in many kinds of multi-mode output states of signals. The ODB in FWM process provide a solution for generating the multi-mode output states at a time, which means the signals can convey more information. With the controlling of switch between different output states, the action of this phenomena can be realized as a dual-bistability flip-flop converter. These results may have novel and promising development for generation and potentially applicable in all-optical devices and quantum information processing.

Experimental setup and Basic theory
A composite atom-cavity system contains a four-level 85 Rb atomic vapor cell in the optical ring cavity. The details of the experimental setup are described in the Methods section. The four relevant energy levels are 5S 1/2 (|0>), 5P 3/2 (|1>), 5D 3/2 , (|2>), and 5D 5/2 , (|3>) as show in Fig. 1(c). In this system, a weak probe beam E 1 (frequency ω 1 , wave vector k 1 , Rabi frequency G 1 , vertically polarized) couples the transition |0> to |1>, while a strong pumping beam E 2 (ω 2 , k 2 , G 2 , horizontally polarized) and another pumping beam E 2' (ω 2' , k 2' , G 2' , vertically polarized) couples upper transition |1> to |2>. Here the detuning Δ i = Ω i − ω i is defined as the difference between the resonant transition frequency Ω i and the laser frequency ω i of E i . Besides, E 1 counter-propagates with E 2 , and E 2' propagates along the optical axis of the cavity having an angle of 2° with E 2 as show in Fig. 1(a). The external-dressing beam E 3 (ω 3 , k 3 , G 3 ) propagates along E 2 direction and drives the other upper transition |1> to |3>. When all incident beams are focused at the center of the rubidium cell by optical lenses, a phase conjugate FWM signal E F (ω F , k F , G F , horizontally polarized) is satisfying the phase-matching condition (k F = k 1 + k 2 − k 2' as show in Fig. 1(b1)) and propagates in the opposite direction of E 2' , which means the signal of E F is mode-matched to the cavity and can form cavity mode. The cavity transmission spectrum of E F leaked from M3 is detected by an avalanche photodiode detector (APD). According to the energy system and Liouville pathways, the generated E F can be obtained by solving the density-matrix equations, the third-order density matrix element as follow: is the Rabi frequency between levels |i>→|j>, and μ ij is the dipole moment; Γ ij = (Γ i + Γ j )/2 is the decoherence rate between |i> and |j>. Then, the generated E F field can oscillate and circulate inside the three-mirrors ring cavity while E 1 , E 2 and E 2' cannot, which is caused by the direction of angle and effect of PBS. Besides, if the atomic transition |0>→|1> is resonant with cavity and the atom cavity coupling is strong enough, the atom-cavity coupling needs to be considered in the system. By solving the master equation of the cavity field evolution and the density matrix operators (the details of the solving process are described in the Methods section), we can get the cavity mode of E F by solving master equation as: ; g is the single-atom-cavity coupling strength and N is the atom number; gN 1/2 represents the atom-cavity coupling which is caused by the combined action between resonant fluorescence of E 1 and vacuum induced enhancement of FWM; the term g 2 N/d 4 represents the cavity dressing effect; the terms |G′ 2 | 2 /d 2 and |G′ 3 | 2 /d 3 are the internal-dressing effect of E 2' and E 3 , respectively. Specifically, by considering FWM process in the atom-cavity system with strong self-Kerr nonlinear effect, there exists an un-neglected cavity feedback effect (also a self-dressing effect) 24 . Clearly, the generated E F field have a self-dressing effect of |G F | 2 , which is derived from the relatively strong feedback effect. This self-dressing effect have similar influence with E 1 , E 2' , E 3 and g 2 N, so Eq. (2) can be rewritten as: Besides, the probe beam E 1 has a sufficiently low power, as well as is far detuned from |0>→|1>. Moreover, with enhanced cavity mode of E F , a SP-FWM process will occur in the system, which can generate two weak fields (Stokes field Es and anti-Stokes field Eas which satisfying the phase-matching condition 2k F = ks + kas). Therefore, the E 1 is naturally injected into the input Stokes or anti-Stokes port of the SP-FWM process as shown in Fig. 1(b2), and the injection will serve as an OPA process (with phase-matching conditions k aS = 2k F − k 1 and k S = 2k F − k 1 ) assisted by the cascaded nonlinear process. The signal of E 1 beam that amplify by OPA-FWM process is named as PTS, and it is detected by the other APD. The photon numbers of the output Stokes and anti-Stokes fields in the amplification process with injection are described in the Methods section.
Further, by turning on probe and coupling fields, the first-order density matrix element with consideration of the dressing effects from E 1 , E 3 and E F is given via Liouville pathways: What's more, similarly to the proposed concept of vacuum induced transparency, the vacuum induced nonreciprocity responses for PTS and FWM are caused by |G F | 2 . As this cavity feedback dressing term |G F | 2 is not equal between the signal curves in the frequency-increasing and frequency-decreasing processes (corresponding to the rising and falling edges in one frequency scanning round trip, respectively), the generated nonreciprocity ODB of folded signals could exist "∞"-shape non-overlapping region for scanning the frequency detuning of probe or cavity, which includes the frequency offset in x direction and the shape change in y direction. At first, there exists frequency offset (Δυ) between the two peaks or dips in the same baseline. Whereas the nonreciprocity can be interpreted by the change of nonlinear refractive index Δn' , which is given as: where term Δσ = Δυn1l/c is the phase delay and Δυ is the frequency difference that can reflect the ODB phenomenon directly; I up (I down ) is the feedback intensities of the signals for the PTS and FWM on the two side ramps generated at the same frequency scan. n 2up (n 2down ) is the nonlinear refractive index coefficient that can be generally expressed as n 2up ≈ n 2down ≈ n 2 =Re[χ (3) /(ε 0 cn 0 )], which is the mainly dominated by the Kerr coefficient of E 2 . The nonlinear susceptibility is: Besides, the shape change of the signals can also advocate the ODB effect in the composite cavity-atom system which can be understood through requirement for the dressing suppression and enhancement. When considering the different feedback dressing on the rising and falling edges, the signal will meet different enhancement (or suppression) conditions. For instance, when we scan Δ 3 , the primary A-T splitting is caused by E 3 whose corresponding eigenvalues are λ ± = [∆ 3 ± (∆ 3 2 + 4|G 3 | 2 ) 1/2 ]/2, so the suppression and enhancement conditions of G 3 are ∆ 1 + ∆ 3 = 0 and ∆ 1 + λ ± = 0, respectively. The secondary A-T splitting is caused by the feedback dressing term G F whose corresponding eigenvalues are λ +± = [∆ 3 ' 2 ± (∆ 3 ' 2 + 4|G F | 2 ) 1/2 ]/2 (∆ 3 ′ = −∆ 1 −λ + ), and the suppression and enhancement conditions are ∆ 1 + ∆′ 3 = 0 and ∆ 1 + λ ± + λ +± = 0, respectively. Further, when we scan Δ ac , the split energy levels are λ +± ± gN 1/2 and λ − , where λ ++± = [∆ ac ′ ± (∆ ac ′ + 4g 2 N) 1/2 ]/2 (∆ ac ′ = −∆ ac − λ +± ), the respective suppression and enhancement conditions are ∆ 1 + ∆ ac ′ = 0 and ∆ 1 + λ ± + λ +± + λ ++± = 0.

Results and Discussion
In our experiment, we investigate the multi-mode output states that is modulated by the different input parameters. Firstly, we obtain the nonreciprocity of the probe transmission signal (PTS) and the cavity mode four-wave mixing (FWM) signal by scanning the frequency detuning of external-dressing field E 3 (Δ 3 ) or cavity length (namely, ∆ac). The cavity length is controlled by the PZT (connecting with plate-concave mirror M1) at on left direction is "frequency-rising edge" and at the other direction is "frequency-falling edge" in Fig. 1(a1). That signals come from electromagnetically induced transparency (EIT) window ∆ 1 + ∆ 2 = 0. The curves of PTS perform (where I 0 is the intensity of the probe field without Doppler absorption; ρ (3) 01 and ρ (1) 10 (in Eq. (6)) show as gain peak and EIA dip, respectively). Meanwhile, the intensity of cavity mode FWM signal is related to a FWM in Eq. (3). In Fig. 1(d), with scanning Δ 3 , the left and right curves of the bell shaped represent the frequency-rising edge and frequency-falling edge, respectively, and the vertical curve represents the turning point of the round trip. Correspondingly, there are the signals of F cavity mode FWM, where the left (right) peaks belong to the signals of rising (falling) edge in one frequency scanning round trip. Besides, the frequency detuning between left peak and turning point T is ∆δ 1 , and the frequency detuning between right peak and T is ∆δ 2 . Therefore, when fold the signals on the two edges (rising and falling) from the maxima of the ramp curves point T do not overlap, the frequency offset can express as δ = ∆δ 2 − ∆δ 1 . In our results, we also use this way (δ = ∆δ 2 − ∆δ 1 ) to analyze the frequency offset. Besides, it is obvious that the shape of right peak is different from left peak. Therefore, the degree of nonreciprocity can be demonstrated by non-overlapping region includes frequency offset and shape change which can be approximately viewed with infinite sidebands, so we named this kind of nonreciprocity as "∞"-shape ODB.
Secondary, the action of ODB phenomena can be realized as a dual-bistability flip-flop converter and schematic diagram as show in Fig. 1(e), where E 1 is trigger signal, E F is input signal, FD (frequency detuning) and PO (power) are input parameters, PTS and FWM are output multi-mode states. There also have two states ("zero" state and "trigger" state) of this flip-flop and it mainly works at "trigger" state. FD is controlled by an electro-optical modulator and the speed is 10 ns, while PO is controlled by an acoustic optical modulator and the speed is 12 ns. The switching speed is controlled by the atomic coherence time that is mutable by the phonon effect from microseconds to nanoseconds. The total switching speed (16 ns) of this flip-flop converter is taken to be the quadrature sum of several independent contributions. In following experimental results (Figs 2-4), with changed input parameters, we realize the fast conversion between different output multi-mode states by analyzing "∞-shape non-overlapping region and optical contrast of folded signals.
In Fig. 2, we discuss the output multi-mode states of the PTS and cavity mode FWM signal (E F ) versus the frequency detuning of E 3 (Δ 3 ). By changing frequency detuning (FD) and power (PO) of dressing fields, we realize the switch action between different output multi-mode states. In Fig. 2(a), the signal of E F is measured when versus Δ 3 at discrete Δ 1 , where the dip of signal curve is lower than the baseline that represents suppression effect of E 3 . The curves of FWM signal is related to a FWM in Eq. (3). The dip gets the maximum value in Fig. 2(a4) where satisfies window ∆ 1 + ∆ 2 = 0 and Δ 1 + Δ 3 = 0 together, which means the value of ∆ 2 is equal to Δ 3 . Since the two peaks in Fig. 2(a) have the same baseline, the cavity feedback dressing (also a self-dressing) term |G F | 2 / Γ 00 on the rising and falling edges are not equal in Eq. (3). Hence, the feedback intensity I up is not equal to I down while n 2up ≈ n 2down , which can lead to the occurrence of frequency offset (Δυ) of ODB. From top to bottom, with the increased detuning Δ 1 , the Δυ appeared larger, which is attributed to increased n 2 in Eq. (6), and the Δυ attains maximum value of 6.7 MHz in Fig. 2(a4) where ∆1 = 0 MHz. In this output state, with H-state of PO, FD is changed from B to S, and the output state of FWM switches from "1,0" to "0,1". Succinctly, the optical contrast for switching application can be defined as C = (I 1 − I 0 )/(I 1 + I 0 ), where I 1 is the intensity of "1"-state and I 0 is the intensity of "0"-state. The optical contrast calculated for this case is 100%. Since, it is obvious from the same baseline that the shape of peaks is different in Fig. 2(a), which is caused by the different suppression conditions. Here we consider second-order splitting caused by E 3 and G F that we have been illustrated theoretically. Because of different G F on the two edges, the suppression conditions ∆ 1 + ∆′ 3 = 0 of two peaks are not same while there are different corresponding eigenvalues λ. And the change of shape difference is minute from top to bottom. While Fig. 2(b,c) are the cavity modes of FWM signals by scanning Δ 3 at different external-dressing field power (P 3 ) and probe power (P 1 ), respectively. Comparing with Fig. 2(a), the signals perform gain peak caused by strong enhancement of cavity-dressing effect of gN 1/2 in Fig. 2(b,c). There frequency offset Δυ is also attributed to cavity feedback dressing term |G F | 2 /Γ 00 . From bottom to top, with increased P 3 , both Δυ and shape difference are increasing gradually in Fig. 2(b), which means the power of external-dressing field can alter the state of FWM. Where PO is changed from L to H and the state of FWM switches from "0,0" to "1,1", while the state of FD is B, and the optical contrast calculated for this case is 78%. But in Fig. 2(c), from bottom to top, the change of Δυ increases from minimum to maximum in Fig. 2(c4)-(c2), then it becomes minimum again in Fig. 2(c1), it also attributes to changed P 1 , which is increased first and then is decreased. In this state, with S-state of FD, PO is changed from L to H and then to L, and the state of FWM switches from "1,1" to "0,0" and then to "1,1" again. The optical contrast calculated for this case is 69%.
In Fig. 2(d-f), there are signals of PTS versus Δ 3 corresponding to Fig. 2(a-c), respectively. The signals are only dip visually, since the suppression is much strong. Same with cavity mode FWM signal, the nonreciprocity ODB of PTS also is induced by cavity feedback dressing term (also parametrically amplified effect) |G F | 2 /Γ 00 in Eqs (3) and (4). In Fig. 2(d), with the Δ 1 increased from top to bottom, the n 2 is increased, the signals in Fig. 2(d5) get the maximum value where Δ 3 around 0 MHz. And the state of PTS also switches from "1,0" to "0,1". In Fig. 2(b,c), all signals satisfy the enhancement condition ∆ 1 + λ ± + λ +± = 0. With the decreasing of P 3 and P 1 in Fig. 2(e1-e6) and Fig. 2(f1-f6), the signal intensity of PTS is increased by the dressing term G 3 2 /d 3 in Eq. (2) and G 1 2 /Γ 00 in Eq. (6), respectively. Besides, the change process of "∞"-shape non-overlapping region in Fig. 2(e,f) are same with Fig. 2(c). In Fig. 2(e3), the Δυ gets maximum value at 31.8 MHz with the largest n 2 . For same reason, the Δυ gets maximum value at 20.1 MHz in Fig. 2(f3). After that, with the decreasing of P 3 or P 1 , the effect of cavity feedback dressing becomes weaken. In these two cases, the states of PTS are switched from "0,0" to "1,1" and "1,1" to "1,1", respectively. And the optical contrasts are more than 90%.
In Fig. 3, with beam E 3 blocked, we discuss the output multi-mode states of the cavity mode FWM signal (E F ) versus the frequency detuning of cavity (Δ ac ). By changing frequency detuning of dressing fields, we realize the switch action between different output multi-mode states. Same with Fig. 2, we recognize the state by analyzing the "∞"-shape non-overlapping region of folded signals.
In Fig. 3(a,b), with beam E 3 off, we measured E F versus Δ ac at discrete Δ1, where the E F signal comes from EIT window Δ 1 + Δ 2 = 0. In Fig. 3(a), we change the detuning Δ 1 in a period but in Fig. 3(b) is just half period. When the detuning Δ 1 is adjusted the position of resonant peaks shifts in the direction of , and the peak gets the maximum value in Fig. 3(a3). Here we changed the detuning Δ2 = 0 satisfying the conditions Δ 1 − Δ ac = 0 and Δ 1 − Δ 2 = 0, and the window Δ 1 − Δ ac = 0 comes from dressing effect of gN 1/2 by the term g 2 N/d 4 in Eq. (3). This equation manifests a positive correlation between Δ ac and Δ 1 with a large moving speed satisfying |dΔ ac /dΔ 1 |>2. Same with Fig. 2, the different feedback dressing parts between I up and I down result in the nonreciprocity of signals. In Fig. 3(a1-a3), the detuning Δ 1 increased first then decreased in Fig. 3(a3-a5). In this case, n 2 is changed from increased to decrease accordingly. The corresponding frequency offset Δυ are 49.1 MHz, 113.9 MHz and 47.0 MHz in Fig. 3(a1,a3,a5), respectively. In this state, with L-state of PO, FD of Δ 1 is changed from S state to B state and then to S state again, and the state of FWM switches from "0,0" to "1,1" and then to "0,0" again. The optical contrast calculated for this case is 79%. While in Fig. 3(b), from top to bottom, the decreased detuning Δ 1 result the n 2 to be decreased and corresponding Δυ is changed from 66.8 MHz to 53.8 MHz. Where the transition of state is changed from "0,0" to "1,1", and the optical contrast is 96%. All peaks of signals in Fig. 3(a,b) satisfy enhancement condition ∆ 1 + λ ± + λ +± + λ ++± = 0. In contrast to left peaks, the shape of right peaks is smaller in Fig. 3(a) and is larger in Fig. 3(a), which is attributed to different cavity feedback dressing effect. In Fig. 3(c), Δ ac is scanned at discrete Δ 2 with all beams turned on and setting Δ 3 = 0. When satisfying conditions Δ 1 + Δ 2 = 0, Δ 1 − Δ ac = 0 and Δ 1 + Δ 3 = 0, we can get the maximum value in Fig. 3(c1). From top to bottom, the increased detuning Δ 2 result the decreased n 2 , and one can Figure 3. All the experiment curves are "T" (trigger) state by measuring cavity mode (FWM) with beam E 3 blocked; B: big; S: small; 0,0: output state of signals; 1,1: output state of signals. In a period, (a) versus Δ ac at discrete Δ 1 from "B" state to "S" state then to "B" state; (b) versus Δ ac at discrete Δ 1 in a half period (from "S" state to "B" state); (c) versus Δ ac at discrete Δ 2 from "S" state to "B" state.
witness that the frequency offset Δυ in Fig. 3(c) is approximately 110 MHz. Where the change of ODB phenomena is very small due to saturation effect. Here the transition of state is similar to Fig. 3(b), but the optical contrast is 80%. Same with Fig. 2, in Fig. 4, we discuss the output multi-mode states of the PTS and cavity mode FWM signal (E F ) versus the frequency detuning of E 3 (Δ 3 ). By changing frequency detuning of dressing fields, we realize the switch action between different output multi-mode states. Using the phenomenon of ODB, we investigate this action by analyzing the "∞"-shape non-overlapping region of folded signals. And same with Figs 2 and 3, the nonreciprocity of signals is also caused by different feedback dressing between two edges (rising and falling edges) in Fig. 4.
In Fig. 4(a,b), the curves represent E F signal by scanning the dressing detuning Δ3 with the external dressing effect of E 3 . In Fig. 4(a), when Δ 1 is decreased first in Fig. 4(a1-a3) then is increased in Fig. 4(a3-a5), the measured E F signal is obtained. Compared to Figs 2(a) and 3(e), the difference is that the signals in Fig. 4(a1,a2) have peaks and satisfy the enhancement condition ∆ 1 + λ ± + λ +± = 0. While the condition of Fig. 4(a3-a5) change from semi-suppression and semi-enhancement to suppression (∆ 1 + ∆ 3 = 0). With the changing of Δ 1 , the frequency offset Δυ is bit-changed in Fig. 4(a1-a3) and then is decreased from 90.2 MHz to 35.6 MHz gradually. In this state, with L-state of PO, FD of Δ 1 is changed from B-state to S-state and then to B-state again, and the state of FWM switches from "0,0" to "1,1" and then to "0,0" again. The optical contrast calculated in this case is 65%. With increased cavity detuning Δ ac from top to bottom, the signals of E F as shown in Fig. 4(b). Where the peaks of signals satisfy the enhancement condition ∆ 1 + λ ± + λ +± = 0, the change of degree of nonreciprocity ODB is similar with Fig. 3(b), but the state of FWM is keeping at "1,1" state with H-state of PO. Figures 4(c,d) are signals of PTS corresponding to Fig. 4(a,b), respectively. In Fig. 4(c), the peaks of signals satisfy the enhancement condition ∆ 1 + λ ± + λ +± = 0. From top to bottom, with changed Δ 1 , the enhancement of peaks is weakened first in Fig. 4(c1-c3) then is strengthened in Fig. 4(c3-c5), but the frequency offset Δυ is little-changed which is attributed to the enhancement effect. Here the transition of PTS state is contrary to Fig. 4(a). Compared with Fig. 4(c), the signals have dips in Fig. 4(d1,d2) that satisfy the suppression condition ∆ 1 + ∆ 3 = 0. Then, the condition of signals change from suppression to semi-suppression and semi-enhancement in Fig. 4(d2-d5). Besides, the degree of "∞"-shape non-overlapping region of ODB phenomena is decreased gradually due to changed suppression and enhancement conditions. Where the transition of state is changed from "0,0" to "1,1", and the optical contrast is 96%.
In the Figs 2-4, the corresponding switching ratio of our flip-flop converter is attributed to the degree of non-overlapping region. If the size of non-overlapping region between dips or peaks in same baseline is big enough, it may be benefit for the division of two states. Such results can be exploited for frequency detuning division multiplexing. Here we use the channel equalization ratio = − ∑ − − P S S S 1 ( ( ) / ) 1/2 to measure the de-multiplexing effect, where s is the area of one dip or peak, s i is the area of non-overlapping region. In Fig. 3(a-c), P is near 100%, we shall obtain more balanced and stable spatial channels. In Figs 2(e,f), 3(d-f) and 4(a,c), the channel equalization P can approach 80-90%. We find that the high channel equalization ratio is caused by large of frequency offset Δυ.

Conclusion
In summary, we study the realization of a dual-bistability flip-flop converter in cavity and PA-FWM. The flip-flop action results from the nonreciprocity ODB phenomena which is induced by self-dressing effect (cavity feedback dressing). The degree of this ODB can reflect by the "∞"-shape non-overlapping region of folded signals, which is changed with different input parameters. Therefore, with high optical contrast and stable spatial channels (channel equalization ratio), we can obtain many kinds of output multi-mode states by controlling the input parameters. By switching the states of PO or FD, we also can realize the convert of output multi-mode states. Specifically, the switch speed of this flip-flop converter is about 16 ns. These results are well explained with theoretical model. The observed phenomenon has novel and promising development for generation and potentially applicable in all-optical devices and quantum information processing.

Methods
Experimental setup. The experiment is performed in a composite atom-cavity system which contains a four-level 85 Rb atomic vapor cell in the optical ring cavity as shown in Fig. 1(a). With the length of 38 cm, the cavity is consisted of a plate mirror M3 (with reflectivity of 97.5% at 780 nm) and two plate-concave mirrors (M1 99.9% and M2 97.5% at 780 nm). M1 is mounted on a PZT for adjusting and locking the cavity length. In the setup, a weak probe beam E 1 counter-propagates with a strong pumping beam E 2 . Another pumping beam E 2' propagates along the optical axis of the cavity (indicated by the dashed line) having an angle of 2° with E 2 . Hence a phase conjugate FWM signal E F propagates in the opposite direction of E 2' , which means the signal of E F is modematched to the cavity and can form cavity mode. The external-dressing beam E 3 propagates at E 2 direction. All incident beams are focused at the center of the optical cavity by optical lenses. The temperature of the cell is set to 75 °C in order to have enough atoms in the cavity to enhance the strength of atom-cavity coupling. The cavity transmission spectrum of E F leaked from M3 (FWM) is detected by avalanche photodiode detector (APD), and the absorption of E 1 (PTS) is detected by the other APD.
Evolution of the cavity field and the density matrix operators. Under the weak-cavity field limitation and with all the atoms initially in the ground state |0〉, the evolution of the cavity field and the density matrix operators obey the following linear equations as: where the over-dots represent the first-order derivative with respect to time t; γ is the decay rate for the cavity; Δ ac = ω 10 − ω c is the detuning of the cavity field with cavity resonant frequency ω c ,  ε μρ , g N results from the resonant florescence induced by E 1 and treated as the strength of atom-cavity coupling with single-atom-cavity coupling strength g and atom number N, and the atom-cavity coupling strength has a similar dressing effect with internal-dressing effects of E 1 , E 2 , E 2' and external-dressing effect of E 3 .

The photon numbers of the output Stokes and anti-Stokes fields.
With the amplification of OPA-FWM process, the photon numbers of the output Stokes and anti-Stokes fields with E 1  denotes the intensity of PA-FWM, r is the radius of the probe field, G AF is the PA-FWM amplification factor, and  ε μ ρ ∝ G c N 2/ 2 for E S and E aS , respectively.