The generation of continuous-variable entanglement frequency comb.

Continuous-variable (CV) entanglement frequency comb can be produced by enhanced Raman scattering in an above-threshold optical oscillator cavity in which a hexagonally-poled LiTaO3 crystal resides as a Raman gain medium. The Stokes and anti-Stokes Raman signals are enhanced by a coupled quasi-phase-matching optical parametric process and the entanglement natures among these Raman signals and pump are demonstrated by applying a sufficient inseparability criterion for CV entanglement. Such entanglement frequency comb source with different frequencies and continuously tunable frequency interval may be very significant for the applications in quantum communication and networks.

A valuable feature of quantum optical implementations based on continuous-variable (CV) entanglement, related to their high efficiency, is their unconditionalness1. Multipartite CV entangled beams, produced generally by combining squeezed beams on linear beam splitters, are degenerate2345. However, multipartite CV entangled beams with different frequencies are necessary to connect different physical systems at the nodes of quantum networks6 since they can be separated easily in the systems. In contrast to the degenerate case, multicolor entanglement has to be generated directly by nonlinear optical processes. It was predicted that two-color CV entanglement can be produced by the nondegenerate optical parametric oscillator (NOPO)7, which has been demonstrated experimentally below8, as well as above9, the oscillation threshold. Recently, theory predicts that three-color CV entanglement can be generated by operating an OPO above its threshold10, which has been verified in a recent experiment11. The three-color CV entanglement beams can also be produced by quasi-phase-matching (QPM)1213 third-harmonic generation14 or forth-harmonic generation15, respectively. The bright quadricolor CV entanglement can be generated by four-wave-mixing process16.

Optical frequency comb (FC) can be generated from multimode OPO, which plays an important role in the application of one-way quantum computation17. More recently, a scheme was proposed to produce large CV cluster state using one single-mode vacuum squeezer and a quantum nondemolition gate18, and another scheme was presented to create multioctave FC by stimulated Raman scattering (RS) in a hydrogen-filled hollow-core photonic crystal fiber19. CV cluster states can also produced in the time domain20. The multipartite entanglement of 60 modes of a quantum optical FC was experimental generated in a bimodally pumped OPO21. The quantum frequency comb is easily addressable, robust with respect to decoherence and scalable, which renders it a unique tool for quantum information22.

In the present paper, we propose a new scheme to directly produce multicolor CV entanglement FC by enhanced RS in an above-threshold optical oscillator cavity, in which a two-dimensional hexagonally-poled LiTaO3 (2D-HPLT) crystal serves as a Raman gain medium. This idea is based on the previous experimental result of Ref. 23 and the elicitation derived from the beautiful experiment of Ref. 11. In this case, we verify, for the first time, that all the Stokes and anti-Stokes Raman frequency components in this comb, including the pump, are mutually entangled by applying the sufficient inseparability criterion for CV entanglement24. This scheme of generating multicolor CV entanglement FC comes directly from a single optical superlattice inside the optical cavity, which has been achieved experimentally without the optical cavity23, and exhibited a symmetrical FC of phonon-polariton (P-P) RS enhanced by a coupled QPM optical parametric amplification processes.

Enhanced RS.

We consider a 2D-HPLT crystal as a Raman gain medium located inside an one-sided optical oscillator cavity, as illustrated in Fig.1(a). Initially, a pump with the frequency of ω0 and a signal with the frequency of ω1 are incident onto the crystal23. Then, a beam with the frequency of ω−1 is generated by the difference frequency generation between the pump and signal beams, and the phase mismatching is compensated by the reciprocal G via a process of QPM as given by

Figure 1

(a) Sketch of the one-sided cavity. (b) The sketched momentum geometry of a coupled quasi-phase-matching process. (c) The Raman spectra corresponding P-P frequency ω = 19 cm−1 23.

The intense P-P field supplies the inelastic source and generates the anti-Stokes beams

and Stokes beams

where stands for the order of the RS.

Figure 1(b) gives the sketched momentum geometry for the QPM-OPA. The high-order RS beams can be significantly amplified by the pump through the process of coupled QPM-OPA, which is related to the second-order nonlinear interaction, and induces the rapid growth of the RS signals25. In the previous experiment, the eighth-order RS beam was observed and the relative intensities of the RS beams approach to 14%–70%. For example, Fig. 1(c) shows the Raman spectrum which was obtained in the previous experiment23. One can see that the Stokes and anti-Stokes frequency components distribute symmetrically and form a comb-like structure with equal frequency interval.

Output fields

The interaction Hamiltonian for the coupled nonlinear processes can be written as26

where , and are annihilation operators for the cavity modes, respectively. κ0, κ and κ are the effective coupling constants, which are proportional to the nonlinear susceptibility and structure parameters, and taken to be real without loss of generality272829.

The equations of motion for the cavity modes can be obtained by solving the Fokker-Planck equation in the positive-P representation30313233. In this case, for example, we only consider the first- and second-order RS beams in the cavity. Then, one can obtain the equations of motion of the five cavity modes as3031

where α and correspond to the mode operators a and in the positive-P representation, respectively, γ (i = 0, ±1, ±2) stand for the damping rates for the corresponding cavity modes, ε0 and ε1 are the classical pump laser amplitudes, and η are complex Gaussian noise terms with and . In order to simplify the calculation we assume γ1 = γ−1 = γ2 = γ−2 = γ and κ2 = κ−2 = κ−1 = κ1 in the above equations. Also, we consider the pump and injected signal as real fields, i. e. and .

The steady-state solutions can be obtained by letting ∂α/∂t = 0 (i = 0, ±1, ±2) and ignoring the noise terms in Eq. (5). The steady-state value A0 of the pump satisfies the following equation

where Δ0 = γ0A0 − E0, , and . The other steady-state values are related to A0 as , , , and , where A (i = ±1, ±2) are the steady-state values of the cavity modes. Because the Eq. (6) is a fifth-order equation, it is well known, there is no analytical solutions for it. So one can not obtained the threshold of the system although one can solve the Eq. (6) by numerical calculation. In order to discuss expediently, we use ε = γ0γ/κ0, which is the well-known threshold for the OPO, as a terminology for the threshold. Generally, the injected signal ω1 can increase the conversion efficiency and stability, and change the characteristics of the threshold.

In the following, one can decompose the system variables into their steady-state values and small fluctuations around the steady-state values as α = A + δα and use this linearized analysis as a method to rewrite the evolution equations in Eq. (5) as3031

where ; M is the drift matrix with the steady-state values; B contains the steady-state coefficients of the noise terms of Eq. (5), and dW is a vector of Wiener increment3031. The condition for the validity of the linearization process is that the eigenvalues of the drift matrix M have no negative real parts3031. The analytical expressions since the eigenvalues of M can't be obtained here for the fifth equation has no analytical solutions. However, we find the eigenvalues are related to κ1 via numerical calculation. When κ1 ≥ κ1, where , the negative real parts only present in the range of ε0 ≥ ε. With the decrease of κ1, when κ1 < κ1, the positive real parts become present in the range of ε0 ≥ εth. In Fig. 2 We plot the real parts of the eigenvalues versus the pump power parameter ε0/εth with γ0 = 0.03, γ = 0.01, κ0 = 0.05, and /κ0 for κ1 = 10κ1, κ1, 0.1κ1, and 0.01κ1, respectively. In present study, we only focus on the OPO operating above the threshold. From Fig. 2 one can see that there are positive real parts appear in the range of above the threshold only when κ1 < κ1 is satisfied. In addition, there are negative real parts appear at about ε0/ε = 6.8 in Fig. 2(c). With the further decrease of κ1, when κ1 = 0.01κ1, one can see from Fig. 2(d) that the real parts of the eigenvalues are all above zero in the range of ε0 ≤ 10ε. Above analysis indicates that the linearization analysis is valid in the case depicted in Fig. 2(c) and (d).

Figure 2

The real parts of the eigenvalues versus the pump power parameter ε0/ε for different κ1: (a)κ1 = 10κ1; (b)κ1 = κ1; (c)κ1 = 0.1κ1; (d)κ1 = 0.01κ1.

Then, from Eq. (7) one can obtain the intracavity spectral matrix in the frequency domain as3031

where I is the identity matrix and ω is the analysis frequency. By applying the well-known input-output relations34 one can obtain the output fields.

CV entanglements in the FC

When q = 1, there are three modes in the cavity, the experiment setup is exactly similar to an above-threshold OPO with an injected signal which has been demonstrated that the three modes, i. e. pump, signal, and idler, are CV entangled with each other35. So when q = 1, one only considers the first order RS, the three modes, i. e. , , and , are CV entangled with each other when the cavity operates above the oscillation threshold. When q = 2, there are five modes in the cavity including two pairs of RS beams and the pump. From the multipartite CV entanglement criteria24, four of the simultaneous satisfactions in the following inequalities

should be sufficient to verify the five-mode CV entanglement. In addition, the first two of the simultaneous satisfactions in above inequalities should be sufficient to demonstrate a0, a1, and a−1 three-mode CV entanglement when five-mode interaction exists. represent the amplitude quadratures of the fields; stand for their phase quadratures; gi are the adjustable parameters chosen to minimize the left of the inequalities.

In Fig. 3(a) and (b) we show that the minima of the five inequalities versus the analysis frequency ω with γ0 = 0.03, γ = 0.01, κ0 = 0.05, ε0 = E0 = 2ε, and for different values of κ1: κ1 = 0.1κ1 and κ1 = 0.01κ1, respectively. One can see in Fig. 3(a) and (b) that the minima of the five inequalities are all below 1 in a wide range of analysis frequency and Eq. (9) is satisfied, which indicates that the two pairs of RS beams and the pump are CV entangled with each other. In addition, the results in Fig. 3(a) and (b) indicate that the entanglement nature of the same-order RS beams is stronger than that between the pump and RS beams as well as among different-order RS beams. In fact, when the optical cavity operates below the threshold value, we don't need to consider the quantum nature of the pump beam since it is undepleted, and only the same-order RS beams are CV entangled36. When ε > ε, the damping of the pump beam can't be ignored, and the quantum nature of the pump should be included. The same-order RS beams are both CV entangled with the pump for the energy conservation and they are amplified through coupled QPM optical parametric processes. The adjacent-order scattering beams should not be entangled because RS is a kind of inelastic scattering process16. However, they are entangled as shown in Fig. 3(a) and (b), the reason for which is that they are both entangled with the pump beam, i. e. the adjacent-order RS beams are CV entangled via the entanglement transfer through the pump beam.

Figure 3

The minima of the inequalities versus the analysis frequency ω for different values of κ1: (a)κ1 = 0.1κ1; (b)κ1 = 0.01κ1.

The minima of five inequalities versus the pump power parameter ε0/ε for different values of κ1: (c)κ1 = 0.1κ1; (d)κ1 = 0.01κ1.

When κ1 = 0, only the first RS exists, the three fields of a0, a1, and a−1 are CV entangled with each other since it is similar to an above-threshold OPO with an injected signal which we have investigated demonstrated that the three fields are CV entangled35. From Figs. 3(a) and (b) one also can see that the first two inequalities are simultaneous satisfied which indicates that the three fields of a0, a1, and a−1 are CV entangled with each other when the five-mode interaction exists. The extended nonlinear interaction do not influence the quantum entanglement nature of the RS beams.

In addition, the degree of entanglement also depends on the intensity difference among the beams. Smaller the intensity difference is, stronger the quantum correlation will be. With the increase of the order of RS, the quantum correlation between the high-order Raman signals and the center of the comb will decrease. Figure 1c), for example, depicts that the Raman beams exhibit a comb-shaped structure with an equal frequency interval which was obtained in the previous experiment for the single-pass nonlinear process23. The relative intensities for the first- to eighth-order Raman signals are obtained to be from 14% to 70%. However, if this setup is placed into an optical cavity, the relative intensity will be significantly enhanced. Then, the quantum correlations between the high-order Raman signals and the center of the comb will increase. In fact, the same order Raman signals are both entangled with the pump beam by the optical parametric amplification process according to Refs. 10, 11, 35, and the adjacent-order Raman signals are also entangled via the entanglement transfer through the pump beam. Therefore, even more optical sidebands are added in the Raman interactions, all the RS signals and the pump are entangled with each other, although the larger intensity difference can decrease the quantum correlation between the high-order RS signals and the pump. Based on above analysis, the demonstrations of three-mode and five-mode entanglement should be enough to verify all the RS beams are entangled with each other. Ultimately, in this scheme a single optical superlattice can directly generate bright multicolor CV entangled FC, which involves multipartite CV entangled beams with different frequencies, having a continuously tunable and equal frequency interval.

Figure 3(c) and (d) depict the minima of the five inequalities versus the pump power parameter ε0/ε with κ1 = 0.1κ1 and κ1 = 0.01κ1, respectively. One can see from Fig. 3(c) and (d) that the minima of the five inequalities are all below 1 in the wide range of pump power parameter which demonstrates the RS beams and pump are CV entangled each other when the optical cavity operates above its threshold. Moreover, one can see that the minima of the five inequalities decrease initially, then increase with the increase of ε0. Initially, with increasing pump above the threshold, RS beams are enhanced, however, the energy of the pump will be transferred to the RS beams. Then, the quantum nature of the pump emerges and the bright multicolor CV entanglement FC can be obtained. With the further increase of the pump, the gain begins to saturate, so the intensity of the pump becomes larger than that of the RS beams, and then its quantum nature will vanish. In addition, there is an abrupt transition at about ε0/ε = 6.8 in Fig. 3(c) which is good agreement with the case in Fig. 2(c). The negative eigenvalues of the drift matrix lead to the abrupt transition at this position for the abrupt increase of the fluctuations.

In summary

We proposed a new scheme to directly produce bright multicolor CV entanglement FC by enhanced RS in an optical cavity with a 2D-HPLT crystal as the Raman gain medium. We have theoretically verified that the Stokes, anti-Stokes RS beams, and the pump are CV entangled by applying sufficient inseparability criteria for multipartite CV entanglement. The adjacent-order RS beams are CV entangled via the entanglement transfer through the pump beam. Because the single-pass process has been experimentally realized, the multimode squeezed light has been produced for both spatial4041 and frequency42 modes, respectively, and the frequency intervals of the RS beams in this scheme are equal and continuously tunable, so it is experimentally feasible to design an optical cavity, in which all the beams oscillate simultaneously.

This entangled FC source, generated from optical supperlattice and with different frequencies and continuously tunable frequency interval and high beam intensity, may be very useful for the applications in quantum communication and networks. For example, the entangled frequency components could be drastically altered using interspecies quantum teleportation protocols, which would enhance the integrability of disparate nodes in a quantum information network637. Similar to the case in Ref. 17, the entangled state generated in the current scheme also has the equal frequency spacing and forms an optical frequency comb. Certainly, not all types of multipartite states are useful for quantum computing3839. The quantum frequency comb is easily addressable, robust with respect to decoherence and scalable, which renders it a unique tool for quantum information.

In the previous experiment23, the eighth-order RS beam was observed and the relative intensities of the RS beams approach to 14%–70%, since the Stokes and anti-Stokes RS beams are significantly enhanced by coupled QPM optical parametric processes. One can obtain at least 17 modes (including the pump) entanglement frequency comb in present scheme. However, when one take this device into an optical cavity, the RS beams intensity will be further enhanced, more high-order RS beams will be observed and these beams are also entangled with each other by the optical parametric processes. Moreover, the frequency intervals of the RS beams in this scheme are equal and continuously tunable. Therefore, comparing to the previously reported quantum optical frequency comb, this scheme should generate better quantum correlations because the utilization of QPM technique can greatly enhance both the Raman signal and the entanglement.

Author Contributions

Y.B.Y. designed the study. X.M.C. directed the study. Y.B.Y. mainly executed the Project and wrote the manuscript. H.J.W., Z.T.S., J.W.Z., F.M.J., Z.Y. and Y.J.W. contributed to the discussion and edited the manuscript.