Hybrid Three-Mode Correlation and Squeezing in a Pr3+:YSO Crystal

We report the generation of three-mode hybrid intensity-noise correlation and intensity-difference squeezing of spontaneous parametric four-wave mixing (SP-FWM) and fourth-order fluorescence (FL) signals in the heteronuclear-like (three-level Λ-type) molecular structure of a Pr3+:Y2SiO5 (Pr3+:YSO) crystal using the nonlinear cross-Kerr effect under a polarized dressing effect. In the semi-classical view of a Kerr nonlinear medium, the amplitude of two-mode hybrid correlations of this kind is subject to a limit determined by the hybrid maximally entangled state. Whereas the degree of correlation and squeezing is determined by the dressing effects of the input laser fields participating in the SP-FWM process. We also find that the variations in magnitude of three-mode hybrid intensity-noise correlation and intensity-difference squeezing are consistent with nonlinear cross-Kerr processes. Such a three-mode hybrid signal may have potential applications in long-distance communication, dense coding, all-optical communication and quantum storage on photonic chips.

has experimentally demonstrated that increasing the number of modes of entanglement enhances the degree of correlation and squeezing 19 , and this phenomenon also shows potential for application to increase the information rate in dense coding 24 .
Motivated by these achievements, two-mode intensity-noise correlation and intensity-difference squeezing of SP-FWM and fourth-order FL signals have been investigated 25 . In this paper, we focus on a hybrid state with three-mode intensity-noise correlation and intensity-difference squeezing of SP-FWM and hybrid FL signals in a Pr:YSO crystal using the nonlinear cross-Kerr effect under polarized dressing fields as Kerr non linearity is directly proportional to the square of the electric field. We investigated the differences between the three-mode hybrid correlation and squeezing effects induced by the cross-Kerr effect under linearly and circularly polarized dressing fields. We found that the degree of three-mode hybrid correlation and squeezing induced by the cross-Kerr effect decreases when the dressing field is changed from a linear polarization to a circular polarization. We also found that the FL signal interacts with the Stokes/anti-Stokes photon pairs and that the degree of hybrid intensity-noise correlation and intensity-difference squeezing can be controlled by controlling the emission of the hybrid FL signal. Our results obtained in this atom-like medium indicate several obvious advantages, such as an improved integration ability compared to atomic vapours and a long coherence time compared to other traditional solid materials. These achievements also have potential applications in all-optical communication and quantum storage on photonic chips.

Results and Analysis
Experimental setup. A rare-earth-ion-doped (0.05% Pr 3+ ) 3-mm YSO (Y 2 SiO 5 ) crystal was used in the experiment. Coupling can occur between Pr 3+ ions localized at different cation vacancies as a result of induced dipole-dipole interactions; consequently, one can treat two such ions as a heteronuclear-like molecule. Therefore, we can construct a Λ-type three-level system (as shown in Fig. 1(a)) by coupling laser fields to the corresponding atomic transition levels, as shown in Fig. 1(b). Two tunable dye lasers (with a narrow scan of 0.04 cm −1 in linewidth) pumped by an injection-locked single-mode Nd:YAG laser and a Continuum Powerlite DLS 9010 with a 10 Hz repetition rate and a 5 ns pulse width were used to generate the pumping field E 1 (ω 1 , Δ 1 ) and the dressing field E 2 (ω 2 , Δ 2 ) with frequency detunings of ω ∆ = Ω − i m n i , where Ω mn is the corresponding atomic transition frequency between levels |m〉 and |n〉 and ω i (i = 1, 2) is the laser frequency. The pumping field E 1 (for which the power, polarization and detuning were fixed at 5 mW, linear polarization and 0 GHz, respectively) and the dressing field E 2 were coupled into the crystal in opposite directions. A quarter-wave plate (QWP) was used to control the polarization state of the dressing field E 2 , and the power P 2 and detuning Δ 2 of the dressing field E 2 could be modified by changing the parameters of the dye lasers. The Stokes and anti-Stokes photons were produced under the phase-matching condition k 1 + k 2 = k S + k aS and the energy conservation condition ω 1 + ω 2 = ω S + ω aS , where k and ω denote the wave vectors and transition frequencies, respectively, of the pumping field, the dressing field and the generated photon pairs. Meanwhile, the hybrid photons were obtained via a Kerr nonlinear interaction under the influence of the higher photonic flux of the pump beam E 1 (ω 1 , Δ 1 ), which passed through the nonlinear crystal as shown in Fig. 1(b). Two photomultiplier tubes (PMT1 and PMT2) were arranged as shown in Fig. 1(b) to detect the Stokes and anti-Stokes signals reflected by PBSs. The FL signals were simultaneously detected by PMT3 with a fast boxcar gated integrator; as shown in Fig. (2), these signals accompany the SP-FWM process under the cross-Kerr effect of E 1 (pumping field) and E 2 (dressing field). The weak Stokes field (E S ), anti-Stokes field (E aS ) and fluorescence field (FL) constitute the hybrid signal, whose propagation can be described as evolution under the Heisenberg picture 26,27 ; see the Methods section for further details.
Results and discussion. Figure 1(c1) and (c2) show the intensity of the FL+SP-FWM signal in the time domain as the power of E 2 increases from 2 mW (shown in (c1)) to 6 mW (in (c2)) when the polarization state of E 2 remains fixed at 0° (linear polarization), the detuning of E 2 is set to Δ 2 = 0 GHz, and the power of E 1 (P 1 ) is held fixed. Figure 1(d1) and (d2) show the intensity of the pure SP-FWM signal in the time domain (obtained under the phase-matching condition for the SP-FWM process) under the same conditions as in Fig. 1(c). The total intensity of the hybrid SP-FWM+FL signal can be written as ρ ρ , and ρ ρ = FL FL (4) 11( ) (4) . Because of the dressing-induced competition effect between the FL and SP-FWM signals, the delay time between them increases 28 with increasing power, as shown in Fig. 1(c1,c2). The delay times are indicated by the splitting distances (in units of time), labelled as D 1 and D 2 , between the SP-FWM signal (left peak) and the FL signal (right peak) in Fig. 1(c). However, at a sufficiently high power, the FL signal is suppressed and the SP-FWM signal is enhanced because of the dressing-induced competition effect, as shown in Fig. 1(c1)-(c2) and 1(d1)-(d2); these findings agree with those of ref. 29. Thus, the FL signal is more sensitive to dressing than is the SP-FWM signal 30,31 . We can make use of these properties of the FL signal to selectively monitor the input dressing laser's polarization, power and frequency detuning in a Λ-type system to control the degree of hybrid three-mode intensity-noise correlation and intensity-difference squeezing. Figure 3 shows the experimental results obtained for two-mode intensity-noise correlation and intensity-difference squeezing at the high power level of 6 mW represented in Fig. 1(c2), with the polarization of E 2 fixed at 0° and the detuning of E 2 set to Δ 2 = 0 GHz. Under the cross-Kerr effect with the dressing field E 2 in the linearly polarized state, Eqs (5) and (6)  10 , respectively. The two-mode intensity-noise correlation between the Stokes and anti-Stokes signals ( Fig. 3(a1)) can be obtained by substituting the intensity-noise signals detected by PMT1 and PMT2 into Eq. (17). Meanwhile, the intensity-noise correlations between the Stokes and FL signals ( Fig. 3(a2)) and between the anti-Stokes and FL signals ( Fig. 3(a3)) can be obtained by substituting the corresponding detected intensities into Eq. (19). Similarly, by substituting the detected intensity data into Eq. (18), we can also obtain the intensity difference δ −Î I ( ) S a S 2 (blue curve in Fig. 3(b1)) and the noise sum δ +Î I ( ) versus ω between the Stokes and FL signals ( Fig. 3(b2)) and between the anti-Stokes and FL signals ( Fig. 3(b3)) can be obtained by substituting the corresponding detected intensity data into Eq. (20). Figure 3(a1) shows the intensity-noise correlation between the Stokes and anti-Stokes signals from the SP-FWM process. A high correlation is observed between these two coherent signals, indicating that the two signals represent strongly correlated 32, 33 and non-classical states. Furthermore, the degree of correlation (the height of the correlation curve) is approximately 1.5 at τ = 0, which corresponds to τ ≈  G a ( ) 1 2 ; consequently, the output correlation can be treated as indicating a maximally entangled state, as the Stokes and anti-Stokes signals are determined by one another through four-wave mixing. The Stokes and anti-Stokes photons are produced under the phase-matching condition and the energy conservation condition. Hence, the Stokes and anti-Stokes signals are determined by each other, and the degree of correlation in the SP-FWM process is large.
By contrast, Fig. 3(a2) and (a3) show the correlation results for the Stokes/SP-FWM+FL and anti-Stokes/ SP-FWM+FL signals, respectively. Here, the hybrid signal (SP-FWM+FL) is obtained via a cross-Kerr nonlinear interaction as shown in Fig. 2, and its propagation is defined as in Eqs (1) and (2), under a non-phase-matching and non-energy-conserving condition between the pumping field and the FL signal. Consequently, the degrees of correlation based on the cross-Kerr nonlinear processes are approximately 1.25 and 1.1, smaller than that for the pure SP-FWM process. As the SP-FWM+FL signal considered in these correlations is a hybrid signal, one can treat this as a hybrid entangled state ψ = + − a a a ( ) ( 0 1 )/ 2, where |0〉 and |1〉 represent Stokes/ anti-Stokes and hybrid (SP-FMW+FL) photons, respectively. The correlation values obtained in this case are slightly lower than those for the maximally entangled Stokes and anti-Stokes states, which are obtained through a coherent process. This result is attributed to the mixture of the SP-FWM and FL signals under the cross-phase modulation of the Kerr interaction, which results in a larger decay rate of the SP-FWM signal compared to the FL signal. The squeezing between the intensity difference and noise sum of the Stokes/anti-Stokes signals (−15.19 dB) is larger than that of both the Stokes/FL signals (−8.92 dB) and the anti-Stokes/FL signals (−9.98 dB), as shown in Fig. 3(b1-b3), consistent with the correlations shown in Fig. 3(a1-a3).
In other words, because the Stokes and anti-Stokes photons are produced under conditions of phase matching and energy conservation, the Stokes and anti-Stokes signals are determined by each other, so the degree of correlation in the SP-FWM process is large. By contrast, there is no phase matching or energy conservation between the pumping field and the FL signal, so although the hybrid Stokes/SP-FWM+FL and anti-Stokes/SP-FWM+FL signals affect each other, their degrees of correlation via the cross-Kerr processes are smaller than that of the SP-FWM process. Figure 3(a) also depicts the shape of the correlation function. A broad correlation peak is obtained from the coherent process, with a linewidth determined by the atomic coherence time. Thus, the SP-FWM signals from the c o h e r e n t p r o c e s s d e c a y r a p i d l y, a n d t h e c o r r e l a t i o n s h a p e i s d e f i n e d a s , as shown in Fig. 3(a), where Γ = Γ ± ± aS S / 2 2 and the parameter ζ represents the bandwidth of the source laser and is constant. By contrast, the correlation shapes obtained from the coherent and Kerr processes together are sharp and described by  interactions between the hybrid FL and Stokes/anti-Stokes signals, the correlation function shape contains both sharp and broad peaks.
In the following paragraphs, we will discuss the three-mode correlation and squeezing results of the three signals. In the Λ-type system shown in Fig. 1(a), the total intensity of the composite signal is (4) . The intensity fluctuations δ τ I ( ) 1 1 , δ τ I ( ) 2 2 and δ τ I ( ) 3 3 were recorded by three detectors, and the triple-beam (S+aS+FL) intensity-noise correlations are shown in Fig. 4(a1-a4). Figures 5(a1,a2) and 6(a1-a3) were plotted by substituting these intensity-noise data into Eq. (21). Similarly, the corresponding three-mode intensity differences δ − −ˆÎ I I ( ) Here, we first study the differences in hybrid three-mode intensity-noise correlation and intensity-difference squeezing induced by the cross-Kerr effect with the dressing field in a linearly (θ = 0°) or circularly (θ = 45°) polarized state. Figure 4 shows the polarization dependence of the three-mode hybrid intensity-noise correlation when the power and detuning Δ 2 of E 2 are fixed to 6 mW and 0 GHz, respectively. Figure 4(a1-a4) shows the intensity-noise correlation between the Stokes, anti-Stokes and FL signals detected in the setup shown in  (a1-a4) show the power dependence of the three-mode hybrid intensity-noise correlation between the Stokes, anti-Stokes and FL signals; (b1-b4) show the corresponding intensity-difference squeezing as the power of E 2 changes from high (6 mW, left) to low (1.5 mW, right) while P 1 is kept constant. Fig. 1(b). A QWP was inserted to control the polarization from linear to circular. The cross-Kerr effect changes with the polarization of the dressing field as given by Eqs (3)-(6).
Different degrees of correlation are observed for different polarization states, as shown in Fig. 4(a1-a4); this behaviour is determined by the Clebsch-Gordan coefficients and the polarization-induced dressing effects. The dressing effect of the field E 2 , which contains the terms θ θ + c G (cos sin )  (11) and (12); this effect causes the SP-FWM intensity to decrease as θ is changed from 0° ( Fig. 4(a1)) to 15° (Fig. 4(a2)), 30° ( Fig. 4(a3)) and, finally, to 45° ( Fig. 4(a4)). In addition, the self-dressing terms of  (11) and (12) cause the overall SP-FWM signal to decrease as the polarization changes from linear to circular. Similarly, the external dressing effect of the field E 2 , which contains the terms θ θ + c G (cos sin ) 4 , in Eq. (14). Thus, the three-mode hybrid intensity-noise correlation between the SP-FWM signals and the associated FL signals decreases in magnitude because of the dressing effect as shown in Fig. 4(a1-a4). Moreover, each field experiences a non-trivial phase shift as it propagates through the nonlinear medium. When a strong field passes through a Kerr medium, it causes the other output fields to undergo a phase shift. However, in our case, the nonlinear phase characterization ( χ 2 ) has a uniform impact of inducing a non-trivial phase on all outputs; this effect is why the correlation values are positive. Therefore, one can model this three-mode hybrid correlation as a three-mode hybrid entanglement of Stokes, anti-Stokes and FL photons in the form ψ = + + a ( ) ( 001 010 100 )/ 3 . Because the amplitude of the correlation does not approach zero, one cannot treat this scenario as one of unentangled hybrid photons. Similarly, the three-mode squeezing between the intensity difference and the noise sum decreases from −4.59 dB (Fig. 4(b1)) to −1.28 dB (Fig. 4(b4)) and follows the same pattern as the corresponding correlation. These results are consistent with the theoretical formulations. In the treatment of the intensity-noise correlation, the three-mode entanglement criteria that determine the quantum correlation variances, described in Eq. , with I 1 = X 1 , I 2 = X 2 and I 3 = X 3 , where the intensity fluctuations correspond to fluctuations of the amplitude quadratures and the corresponding phases are Y 1 = Y 2 = Y 3 = 0 because of the uniformity of the Kerr effect on all outputs. The results presented in Fig. 4 clearly show the correlation between the two distinct states.
As the polarization of E 2 changes from 0° to 45°, the anti-Stokes signal becomes small, as the Stokes and anti-Stokes signals are produced and determined by each other and the correlation between them is still strong. However, the amount of FL in the hybrid signal that is created by the Stokes and anti-Stokes signals will be reduced; hence, the degree of correlation between the Stokes/anti-Stokes and hybrid FL signals will decrease as the polarization of E 2 changes from 0° to 45°. This effect will induce three-mode correlation changes, as shown in Fig. 5(a1-a4). Because the classical photonic flux part of the SP-FWM+FL signal has an equal role in the three-mode entanglement, the entanglement becomes a three-mode hybrid entanglement. The greater the FL intensity is, the greater is the hybrid entanglement. The shape of the three-mode correlation function is given by Eq. (22). Moreover, the degree of three-mode hybrid correlation and squeezing induced by the cross-Kerr effect decreases as the dressing field changes from being linearly polarized to being circularly polarized.
1 2 versus the delay times (τ 1 , τ 2 ) of the Stokes, anti-Stokes and fourth-order FL signals, and (b1-b3) show the corresponding three-mode intensity-difference squeezing as the detuning Δ 2 is scanned from −200 GHz to +200 GHz for θ = 0°. Figure 5 shows the three-mode hybrid intensity-noise correlation and intensity-difference squeezing of the Stokes, anti-Stokes and hybrid FL composite signals under the linearly polarized cross-Kerr effect as the power decreases from high (6 mW, left) to low (1.5 mW, right) while the polarization and detuning of E 2 remain fixed at 0° and 0 GHz, respectively. Figure 5(a1-a4) show the three-mode intensity-noise correlation of the SP-FWM and hybrid FL signals, which decreases as the power decreases. This result is explained by the decrease in the total intensity of the hybrid signals as power of the dressing field decreases. When the power decreases, the external dressing terms G /d 2 2 10 , G /d 2 2 20 , and Γ G / 2 2 10 and the self-dressing terms for the field E 1 in Eqs (7) and (8) suppress the SP-FWM signal, which is consistent with what we observe in Fig. 1(d). In Fig. 1(d), curves (d1) and (d2) correspond to low and high powers, respectively. However, the gain effect of E 2 , determined by G G  (7) and (8). Thus, at low power, the hybrid FL signal intensity dominates, and at high power, the SP-FWM signal intensity dominates. Therefore, in the correlation function, the magnitude of the broad shape becomes small and the magnitude of the sharp shape becomes large as the power of E 2 decreases. The theoretical formulations agree well with the intensity distributions shown in Fig. 1(c) and (d) for hybrid signals and pure SP-FWM signals, respectively. Because the overall intensity decreases at low power, the correlation peaks also decrease from left to right, as shown in Fig. 5(a1-a4). Figure 5(b1-b4) show that the corresponding squeezing between the intensity difference and the noise sum decreases from −3.63 dB (Fig. 5(b1)) to −0.64 dB (Fig. 5(b4)).
In Fig. 6, we present the degree of correlation and squeezing of the composite signals under the linearly polarized cross-Kerr effect as the frequency detuning Δ 2 changes from −200 GHz to +200 GHz while the power and polarization of E 2 remain fixed at 6 mW and 0°, respectively. The three-mode signals from PMT1, PMT2 and PMT3 ( Fig. 1(b)) were simultaneously recorded. Under off-resonance conditions (large detuning, Δ 2 ), the detuning terms ( = Γ + ∆ − ∆ d i( ) 10 10 1 2 , = Γ + ∆ d i 20 20 2 ) in Eqs (7) and (8)  2 and = Γ − ∆ d i 12 12 2 , enhance the emission of the hybrid FL signal at off-resonance points and suppress it at resonance (zero-detuning) points. This result is consistent with the results of ref. 29, in which the FL signal was found to be more sensitive to dressing than the SP-FWM signal. Thus, the squeezings between the intensity difference and noise sum in Fig. 6(a1) and (a3) are −5.41 dB and −5.68 dB, respectively, and are much larger than that in Fig. 6(a2) (−1.09 dB). Meanwhile, the detuning terms at off-resonance points in Eqs (3) and ρ − FL aS 20( ) , respectively. Therefore, the cross-Kerr effect is enhanced under off-resonance conditions for a linearly polarized dressing field. Consequently, the magnitude of the broad shape in the correlation function is smaller at resonance points and larger at off-resonance points, whereas the magnitude of the sharp shape is larger at resonance points and smaller at off-resonance points. Because of the interactions between the hybrid FL signal and the Stokes/anti-Stokes signals, the correlations between the hybrid FL signal and the Stokes/anti-Stokes signals are lower at resonance points. Thus, the degree of hybrid intensity-noise correlation and intensity-difference squeezing can be controlled by controlling the amount of FL emission via the dressing field.
In conclusion, the two-and three-mode hybrid intensity-noise correlation and intensity-difference squeezing of SP-FWM and FL signals was studied by selectively controlling the polarization, power and frequency detuning of the input dressing laser in a Λ-type three-level atomic system in a Pr 3+ :Y 2 SiO 5 crystal under the polarized cross-Kerr effect. The degree of three-mode hybrid correlation and squeezing degree induced by the cross-Kerr effect was found to decrease as the dressing field was changed from being linearly polarized to being circularly polarized. It was found that the degree of hybrid intensity-noise correlation and intensity-difference squeezing can be controlled by controlling the amount of FL emission from two Kerr processes. In particular, the two-mode correlation obtained from the interaction of the SP-FWM and Kerr processes is lower than that of the maximally entangled Stokes/anti-Stokes state. The variations in magnitude of the two-and three-mode intensity-noise correlation and intensity-difference squeezing are consistent with each other. Such a three-mode hybrid signal may have potential applications in long-distance communication, dense coding, all-optical communication and quantum storage on photonic chips.

Methods
Theoretical model. In this experiment, the correlation was calculated from the intensities of the signals obtained from three processes: coherent signals were obtained from the SP-FWM process, and non-coherent (hybrid) signals were obtained from two cross-Kerr nonlinear processes. The SP-FWM process produces Stokes and anti-Stokes signals. The two cross-Kerr nonlinear processes produce two hybrid FL+Stokes and FL+anti-Stokes signals via the interaction of the FL with the Stokes and anti-Stokes photons, respectively.
To generate the hybrid signal, we used a scheme similar to the one explained in ref. 5, where the cross-Kerr nonlinearity is responsible for the interaction of the SP-FWM and FL signals in the nonlinear crystal.
Cross-Kerr nonlinearity is essential for the generation of an entangled state. The Kerr nonlinearity is characterized by a refractive index of n 0 + n 2 |E| 2 , where n 0 is the weak-field linear refractive index term and n 2 = 3χ (3) /2n 0 is a nonlinear refractive index term that is proportional to the field strength |E| 2 . Finally, the nonlinear refractive index can be characterized as n 0 + 3χ (3) /2n 0 |E| 2 .
Because of the coherence of the nonlinear medium, the interaction between the SP-FWM and FL outputs results in a hybrid signal (SP-FWM+FL), whose propagation can be described as follows: Scientific RepoRts | 7: 1743 | DOI:10.1038/s41598-017-01860-8  where θ is the rotation angle of the QWP's axis from the x axis, the c x,y are the anisotropy factors in different directions in the crystal, Γ ij is the transverse decay rate, and |G i | 2 is the Rabi frequency of the field E i. The field intensities of these hybrid signals are given by = aS aS . Equations (1) and (2) describe two cross-Kerr nonlinear processes via the interaction of the FL with the Stokes and anti-Stokes signals, respectively. The first hybrid signal is obtain from the interaction of the Stokes and FL signals (Eq. (1a)), and conversely, the hybrid FL signal may produce the Stokes signal (Eq. (1b)). The second hybrid signal is obtained from the interaction of the anti-Stokes and FL signals (Eq. (2a)); similarly, this hybrid FL signal may produce the anti-Stokes signal (Eq. (2b)). Such Kerr nonlinearity is achieved under a strong magnitude of the pumping flux, which creates an interaction between the coherent SP-FWM and FL signals. The combination of a common energy level, a strong pumping field, a dressing field and propagation through a common medium creates such an interaction between the generated signals and the FL. This effect also allows one to create strongly correlated states 32,33 .
In our scheme, we have demonstrated the correlation and squeezing of a hybrid signal consisting of Stokes, anti-Stokes and FL signals. In our mathematical model, we first adopt pure perturbation theory to obtain the Stokes and anti-Stokes signals under the weak-field approximation. Then, we incorporate dressing perturbation theory (DPT), in which the dressing terms of a strong field are included. In DPT, the density matrix elements of the Stokes and anti-Stokes signals and the associated FL signals are derived using the Liouville equation is the Hamiltonian; ˆ † a a ( ) S S and ˆ † a a ( ) aS aS are creation and annihilation operators acting on the Stokes and anti-Stokes signals, respectively; v is the group velocity; and g is the nonlinear gain. For a Λ-type three-level system with two strong pumping fields E 1 and E 2 , when the self-dressing effect of E 1 and the external-dressing effect of E 2 are considered, the third-order nonlinear density matrix elements of E S and E aS obtained via the perturbation chains ρ ρ ρ ρ (2) 20( ) , respectively, are given by 2 0 2 2 1 , = Γ + ∆ d i 21 21 2 , and = Γ − ∆ d i 12 12 2 . When a QWP is inserted, it polarizes the light beam E 2 into its components along the x and y axes. Thus, considering the different Clebsch-Gordan coefficients for various transitions, the Rabi frequency |G 2 | 2 in Eqs (7) and (8)  Thus, in the dressed state picture, we obtain the third-order density matrices of the corresponding Stokes and anti-Stokes signals by substituting into Eq. (10) the new Rabi frequencies that appear in Eqs (7) and (8), as follows:  respectively.
In the Λ-type three-level system shown in Fig. 1(a), with both the E 1 and E 2 fields open, the fourth-order FL signal is generated as depicted in Fig. 2 and is given as 2 , and = Γ − ∆ d i 12 12 2 . When dressing terms due to polarization changes of the field E 2 are considered, the corresponding equation obtained by modifying Eq. (13) is For two weak multimode fields E S (t, z) and E aS (t, z), under the dressing conditions of the driving fields E 1 and E 2 and the light-matter interaction induced by the finite pulses shown in Fig. 1(a)  . This pair of equations describes the SP-FWM process. The Stokes and anti-Stokes signals are thus created and determined by each other.
These three nonlinear processes, the SP-FWM process and the two cross-Kerr nonlinear processes, create the hybrid FL state; hence, the degree of hybrid intensity-noise correlation and intensity-difference squeezing can be controlled by controlling the amount of FL in the hybrid signal.
The generation of Stokes and anti-Stokes fields via the SP-FWM process is illustrated in Fig. 1(a) , , (3) , the output signals from the FWM process can be modified by applying a dressing field and controlling the power, frequency detuning and polarization of that dressing field. Below, the results for three-mode hybrid intensity-noise correlation and intensity-difference squeezing are analysed.
In the two-mode case, the intensity-noise correlation from the coherent process is given by For this two-mode correlation, the degree of intensity-difference squeezing from the coherent SP-FWM process is given by ref. Similarly, for the two cross-Kerr nonlinear processes illustrated in Fig. 2 (see Eqs (1) and (2)), each involving one signal from the SP-FWM process, the two-mode hybrid correlations are given by is the mean square deviation of the noise sum.
The three-mode correlation of the intensity fluctuations for the three optical beams τ τ τ G ( , , ) The shape of the three-mode correlation function is described as where B and C are parameters that affect the magnitude of the correlation, Ω is the Rabi frequency, L is the length of the medium, and κ is the nonlinear coefficient. The function given in Eq. (22) is related to τ and, as a result, determines the shape of the three-mode correlation function τ τ τ G ( ) . To study the correlations and squeezing of Stokes, anti-Stokes and FL signals transmitted through Pr 3+ :YSO, the dependence of the signal intensities on time, 〈I i 〉 + δI i (τ) (i = aS, S, FL), was recorded using digital oscilloscopes and PMTs. Here, the 〈I i 〉 are the average intensities of the anti-Stokes, Stokes and FL beams, and the δI i (τ) are the corresponding intensity fluctuations. Then, the intensity fluctuations were averaged by a fast gated integrator over 10 pulses. The temporal waveform correlations between the Stokes, anti-Stokes and FL signals were investigated based on Eqs (17), (19) and (21) under a variety of conditions. Additionally, using the intensity fluctuations recorded by the PMTs, the measured intensity fluctuations in the Stokes, anti-Stokes and FL channels were subtracted from and added to each other and were then analysed using a spectrum analyser to investigate the relative noise power between the intensity difference and the noise sum, as described by Eqs (18), (20) and (23).
Furthermore, the three-mode entanglement condition is given by