Optimal multi-photon entanglement concentration with the photonic Faraday rotation
Zhou Lan1, 2, Wang Dan-Dan3, Wang Xing-Fu1, Gu Shi-Pu3, Sheng Yu-Bo2, 4, †
College of Mathematics & Physics, Nanjing University of Posts and Telecommunications, Nanjing 210003, China
Key Laboratory of Broadband Wireless Communication and Sensor Network Technology (Ministry of Education), Nanjing University of Posts and Telecommunications, Nanjing 210003, China
College of Electronic Science and Engineering, Nanjing University of Posts and Telecommunications, Nanjing 210003, China
Institute of Signal Processing Transmission, Nanjing University of Posts and Telecommunications, Nanjing 210003, China

 

† Corresponding author. E-mail: shengyb@njupt.edu.cn

Abstract

We put forward an optimal entanglement concentration protocol (ECP) for recovering an arbitrary less-entangled multi-photon Greenberger–Horne–Zeilinger (GHZ) state into the maximally entangled GHZ state based on the photonic Faraday rotation in low-quality (Q) cavity. In the ECP, only one pair of less-entangled multi-photon GHZ state and one auxiliary photon are required, and the concentration task can be realized by local operations. Moreover, our ECP can be used repeatedly to further concentrate the discarded items of conventional ECPs, which can increase its success probability largely. Under the practical imperfect detection condition, our protocol can still work with relatively high success probability. This ECP has application potential in current and future quantum communication.

1. Introduction

In the last few decades, quantum information technology has developed with extremely high speed. Both theory and experimental results demonstrate that the quantum communication and computation are much more efficient than their classical counterparts in many aspects. Entanglement is an extremely important resource in almost all the practical quantum communication and computation tasks, such as quantum cryptography,[1] quantum teleportation,[2,3] quantum dense coding,[4] quantum secure direct communication,[59] quantum repeaters,[10,11] entanglement-based quantum key distribution,[12,13] and some other quantum communication fields.[1418] During the past decade, research groups have successfully generated many different single-particle and multi-particle entanglement states.[1924] Among the large number of entanglement states, the maximally entangled state is one of the most important entanglement states, which is often adopted to construct the quantum entanglement channel in almost all quantum communication tasks. However, as the entangled particles interact with the channel noise during the storage and transmission process, the decoherence is inevitable. The decoherence can make the maximally entangled state degrade to the mixed state or pure less-entangled state, which cannot be used to set up the high quality quantum entanglement channel. Even worse, they may make the quantum communication insecure.[25] Therefore, both the two kinds of decoherence states should be recovered into the maximally entangled states prior to the practical application.

Entanglement concentration is an efficient method for distilling the maximally entangled quantum states from pure less-entangled states.[2663] The concept of entanglement concentration protocol (ECP) was proposed by Bennett et al. in 1996, which is called the Schmidt projection method.[26] In the entanglement concentration protocol, some collective and nondestructive measurements are required, which cannot be realized under current experimental conditions. After that, various research groups proposed many efficient entanglement concentration protocols in the optic and solid systems, successively. For example, in 1999, the entanglement concentration protocol based on entanglement swapping was proposed by the Bose group.[27] Soon after, Shi et al. put forward an improved entanglement concentration protocol.[28] In 2001, Zhao et al. and Yamamoto et al. independently proposed two similar entanglement concentration protocols,[29,30] in which they first adopt the optical polarization beam splitters (PBSs) to complete the parity check measurement for the photons. In 2008, Sheng et al. constructed the quantum non-demolition detector with the help of the cross-Kerr nonlinearity. With the quantum nondemolition detector, they proposed a recyclable entanglement concentration protocol, which can increase the success probability largely.[31] In 2012, Sheng et al. proposed two entanglement concentration protocols assisted with a single photon, which can largely economize the entanglement resource.[32,33] So far, most entanglement concentration protocols cannot deal with the multi-particle entanglement, for its high operational complexity. Only a few entanglement concentration protocols can deal with the decoherence problem in the multi-photon system.[5060]

The cavity is a highly-efficient platform for performing the quantum information tasks, because the interaction between atoms and photons in the cavity can be well controlled.[6478] In the quantum information tasks, we can adopt the spatially separated cavities to be quantum nodes by making the atoms strongly interact with the local cavities. They can also construct a quantum network with the help of the photons, which act as quantum buses.[70,7476] In 2008, the input-output process relevant to optical low-quality (Q) cavities has been achieved.[79] In 2009, the group of An obtained a general expression for the input-output process of the photon-atom system in the low-Q cavity, which can cause a certain angle rotation of the photonic polarization, that is, the photonic Faraday rotation.[80] In 2012, with the help of the photonic Faraday rotation effect, Peng et al. proposed an entanglement concentration protocol for the multi-photon Greenberger–Horne–Zeilinger (GHZ) state.[81] This entanglement concentration protocol is not optimal, for they need to consume two pairs of pure less-entangled multi-photon GHZ states to distill one pair of maximally entangled GHZ states with some probability, and the other one pair of less-entangled multi-photon GHZ states is wasted. In this paper, we put forward an optical entanglement concentration protocol for distilling the maximally entangled photonic GHZ state from one pair of less-entangled photonic GHZ states and an auxiliary single photon based on the photonic Faraday rotation. Our entanglement concentration protocol is an improvement of the entanglement concentration protocol in Ref. [81], for it can save a pair of less-entangled GHZ states and reach the same success probability of Ref. [81]. Moreover, our entanglement concentration protocol is recyclable. By repeating the entanglement concentration protocol, the discarded items of Ref. [81] can also be concentrated, which can largely increase the success probability. Especially, under practical imperfect experimental conditions, our entanglement concentration protocol shows greater application potential than the entanglement concentration protocol proposed in Ref. [81].

The paper is organized as follows. In Section 2, we introduce the photonic Faraday rotation in low-Q cavity. In Section 3, we explain the entanglement concentration for the two-photon state. In Section 4, we introduce the entanglement concentration for the multi-photon GHZ state. In Section 5, we provide a discussion and conclusion.

2. Basic theoretical model of the photonic Faraday rotation

We briefly explain the basic principle of the photonic Faraday rotation. Suppose a three-level atom is trapped in the low-Q cavity. As shown in Fig. 1,[8285] |g L〉 and |g R〉 are two degenerate ground states and |e〉 is an excited state. Suppose a single photon pulse |φin〉 = (1/√2)(|L〉+|R〉) (|L〉: left-circularly polarization, |R〉: right-circularly polarization) enters the low-Q cavity and interacts coherently with the atom. During the interaction process, the atom will achieve the transition |g L〉↔|e〉 (|g R〉↔|e〉) by absorbing or emitting a photon in |L〉 (|R〉), respectively. As shown in Refs. [80] and [86]–[88], the general expression of the reflection coefficient of the atom–cavity system can be written as

Fig. 1. (color online) The basic principle of the interaction between the three-level atom trapped in the low-Q cavity and the single-photon pulse. (a) A three-level atom is trapped in the low-Q cavity. (b) The three-level atom has an excited state |e〉 and two degenerate ground states |g L〉 and |g R〉, respectively. The states |g L〉 and |g R〉 couple with a left (L) polarized and a right (R) polarized photon, respectively.

In Eq. (1), κ and γ are the cavity damping rate and atomic decay rate, respectively. ω p, ω c, and ω 0 are the frequency of the input photon, the cavity, and the atom, respectively. g is the atom–cavity coupling strength. Equation (1) can be rewritten as a pure phase shift r(ω p) ≈ eiϕ . If the atom uncouples to the cavity, which makes g = 0, we can simplify Eq. (1) to

(2)
Equation (2) can be rewritten as r 0(ω p) = e iϕ 0 .

Therefore, under the atom state of |g L〉, the photon state will evolve to

(3)
while under the atom state of |g R〉, the photon state will evolve to
(4)
It can be found that the polarization direction of the output photon rotates an angle or , which is the so-called photonic Faraday rotation.

According to Eqs. (1) and (2), it can be seen that under the special condition that ω 0 = ω c, ω p = ω c-κ/2, γ = 0, if g = κ/2, we can obtain r(ω p) = −1 (ϕ = π), while if g = 0, we can obtain r 0(ω p) = i(ϕ 0 = π/2). In this way, the output photon–atom state will evolve to[81]

(5)

3. Efficient concentration for two-photon state

Here, we first explain the concentration for the two-photon state, which is the simplest case. The basic theoretical model is shown in Fig. 2. Alice and Bob share a less-entangled two-photon state generated from the entanglement source S1 with the form of

(6)
where α and β are initial entanglement coefficients, and |α|2+|β|2 = 1.

Fig. 2. (color online) The basic theoretical model of our ECP for the two-photon state. S1 and S2 are the entanglement source and single photon source, respectively. The auxiliary three-level atom is trapped in a low-Q cavity c. QWP represents a quarter-wave plate, and PBS means a polarization beam splitter. D1 and D2 represent two single-photon detectors.

For realizing the entanglement concentration, Alice makes the single photon source S2 generate an auxiliary single photon as

(7)
Meanwhile, she makes a three-level atom with the following form: |ψ〉 = (1/√2)(|g L〉 + |g R〉) encode in a low-Q cavity c. Therefore, the whole photon-atom system can be written as
(8)

Next, Alice makes the photons in a1 and a2 modes enter the low-Q cavity c to interact with the three-level atom, successively. Based on the input–output relation rule in Eq. (5), |φ〉 can evolve to

(9)

After the photons exit the cavity, Alice performs a Hadamard (H) operation on the photon in a4 mode and the three-level atom. The H operation on the auxiliary photon can be easily performed by making the photon pass through a quarter-wave plate (QWP), which makes

(10)
where |H〉 and |V〉 represent the horizontal and vertical polarization, respectively. The H operation on the atom can be performed by driving it with an external classical field (polarized lasers), which will lead to
(11)
After the H operations, the state in Eq. (9) can ultimately evolve to the following equation
(12)

Then, Alice makes the photon in a4 mode enter the polarization beam splitter (PBS), which can transmit the photon in |H〉 and reflect the photon in |V〉, respectively. Finally, she measures the photons in the two output modes by the detectors D1 and D2, and the three-level atom. Based on her measurement results, |φ 2〉 in Eq. (12) may collapse to four possible cases. For example, |φ 2〉 will collapse to

(13)
when the measurement result is , while |φ 2〉 will collapse to
(14)
when the measurement result is . Both and are the maximally entangled two-photon states. If is obtained, Alice can transform to by performing the phase-flip operation. Until now, our concentration process is completed, and we recover the less-entangled two-photon state into the maximally entangled state. The success probability is P 1 = 2|αβ|2, which equals to that of the ECP in Ref. [81].

On the other hand, |φ 2〉 in Eq. (12) will collapse to

(15)
when the measurement result is , equation (12) will collapse to
(16)
when the measurement result is . With the help of the phase-flip operation, Alice can also convert to . Interestingly, it can be found that has the similar form with |ϕ 0〉 in Eq. (6). In this way, is a new less-entangled two-photon state. According to the concentration step described above, Alice can reconcentrate for the next concentration round.

In the second concentration round, S2 generates another auxiliary single photon in a 2 mode with the form

(17)
Alice makes the two photons in her location interact with the three-level atom in the low-Q cavity, successively. Similarly, based on Eq. (5), the whole photon–atom state can evolve to
(18)
After the H operation, |φ 3〉 can evolve to
(19)

Finally, by measuring the output auxiliary photons and the three-level atom, |φ 4〉 in Eq. (19) can also collapse to four possible cases. It will collapse to or , corresponding to the measurement result of or respectively. Therefore, we can obtain the maximally entangled two-photon state in the second concentration round, with the success probability of

(20)
where ‘+’ corresponds to and ‘−’ corresponds to . Both of them can be reconcentrated for the third round.

In this way, it can be found that by providing an auxiliary photon with the form

(21)

Similar to in Section 3, the single photon source S2 generates an auxiliary photon with the form of Eq. (7) in a 2 mode. Alice makes the photons in a 1 and a 2 modes pass through the low-Q cavity c and interact with the three-level atom, successively. After the cavity, the whole photon–atom state will evolve to the following expression

(22)

Then, Alice needs to perform the H operation on the photon in a 4 mode and the atom. After that, |φ N〉 in Eq. (22) will ultimately evolve to

(23)

Finally, the states of the auxiliary single photon and the atom are measured by Alice. If she obtains , will collapse to

(24)
while will collapse to
(25)
corresponding to the measurement result of . Both and are maximally entangled N-photon GHZ states. Alice can also convert to with the help of the phase-flip operation. Therefore, we have successfully recovered the less-entangled N-photon GHZ state into the maximally entangled state. The success probability is still P 1N = 2|αβ|2, which is the same as the success probability P 1 in Section 3.

On the other hand, if Alice obtains the measurement result of , in Eq. (23) will collapse to

(26)
while if she obtains , will collapse to
(27)
which can also be converted to |ϕ N2〉. It is obvious that |ϕ N2〉 is a new less-entangled N-photon GHZ state, which can be reconcentrated for the next round. In this way, Alice only needs to provide a suitable auxiliary single photon with the form in the K-th concentration round; she can repeat the ECP to distill the maximally entangled GHZ state from the discarded items. The success probability in each concentration round is totally the same as that for the concentration of the two-photon state.

5. Discussion and summary

In our ECP, based on the photonic Faraday rotation in the low-Q cavity, we can probabilistically recover the less-entangled multi-photon GHZ state into the maximally entangled state. Our ECP is an improvement for the ECP in Ref. [81] for the following reasons. First, the ECP in Ref. [81] needs to consume two pairs of less-entangled multi-photon GHZ states to distill one pair of maximally entangled states, while our ECP only requires one pair of less-entangled multi-photon GHZ states and an auxiliary photon. It makes our ECP much more economical. Second, in Ref. [81], after performing the ECP, all the N parties should measure the photon states. Then, they can judge whether they obtain the maximally entangled GHZ state by checking all the measurement results of the N parties. It greatly increases the operational complexity, especially under the case of large photon number N. Our ECP only requires local operations by Alice, and she can judge whether the concentration is successful only by her own measurement results. Under the successful concentration case, she tells the other N−1 parties to preserve their photons. In this way, it can simplify the operation largely. Third, if the concentration fails, Alice can repeat it to concentrate discarded items in Ref. [81]. It makes the success probability of our ECP much higher than the ECP in Ref. [81].

In the concentration process, both the atom and photon detections play key roles. Under current experimental conditions, both detections are imperfect: their detection efficiencies are lower than 100%. Therefore, it is necessary to consider the influence of the detection efficiency on the success probability. Here, we suppose that η p and η a are the single photon detection efficiency and the single atom detection efficiency, respectively. In the ECP of Ref. [81], as the parties should measure a pair of N-photon states and the single atom state, its success probability can be written as

(28)
It is obvious that is exponentially reduced with the growth of the photon number N. Fortunately, in each concentration round of our ECP, only the auxiliary single photon and the single atom should be measured. According to the description of Section 3 and Section 4, the success probability in each concentration round can be calculated as
(29)
The total success probability (P total) of our ECP equals to the sum of the success probability in each concentration round. Therefore, we can calculate it as follows:
(30)
It can be found that the photon number N does not influence the P total.

Recently, Heine et al. reported their research result on the single atom detection. They have achieved η a = 66% in the experiment.[89] Moreover, they have shown that with some improvement, the single atom detection efficiency can achieve η a > 95% in theory. The single photon detection has been a big difficulty under current experimental conditions, for the quantum decoherence effect of the photon detector.[90] In the optical range, η p is usually less than 30%.[90,91] Lita et al. reported their experimental result about the near-infrared single-photon detection. They showed the η p at 1556 nm can reach 95± 2%.[92] Based on their research results, we make the numerical simulation on the total success probability (P total) of both ECPs. Figure 4 shows the P total varying with the entanglement coefficient α.

Fig. 4. (color online) The success probability (P total) of the ECP in Ref. [81] and our ECP. Here, both the photon state detection and atom state detection efficiencies are controlled to be 80%. (a) P total of our ECP with the recycle number of 5. (b) P total of the ECP in Ref. [81] with N = 5. (c) P total of the ECP in Ref. [81] with N = 10.

In Fig. 4, we assume η p = η a = 80%. In curve a, we repeat our ECP 5 times, while in curve b and curve c, the photon number for the ECP in Ref. [81] is controlled to be N = 5 and N = 10, respectively. It is obvious that P total of the ECP in Ref. [81] reduces largely with the growth of N. When N is large, the . On the other hand, as shown in curve a, by repeating our ECP, the P total can be largely increased. Meanwhile, as the P total of our ECP is not influenced by the photon number N, it can still keep a high value in the practical application with a large value of N.

In summary, we proposed an optimal ECP to distill the maximally entangled multi-photon GHZ state from an arbitrary less-entangled GHZ state based on the input–output relations in the low-Q cavity. Our ECP is an improvement on the ECP in Ref. [81]. On one hand, the ECP in Ref. [81] requires two copies of less-entangled multi-photon GHZ states to be consumed, while our ECP only requires one copy of a less-entangled multi-photon GHZ state and an auxiliary photon. It makes our ECP much more economical. On the other hand, the ECP in Ref. [81] relies on the nonlocal operations from N parties, while our ECP only requires the local operation from Alice. In this way, the experimental operation of our ECP is much easier. Moreover, the distilled photon state can be well preserved and we can repeat the ECP to reuse the discarded items in the ECP of Ref. [81] by repeating our ECP. It makes the success probability of our ECP much higher than that of the ECP in Ref. [81]. Under practical imperfect detection conditions, the success probability of the ECP in Ref. [81] is exponentially reduced with the photon number N, while in our ECP, a relatively high success probability can still be obtained, which is not influenced by N. Based on all the features mentioned above, our ECP has a good application potential in current and future quantum communications.

Reference
[1] Gisin N Ribordy G Tittel W Zbinden H 2002 Rev. Mod. Phys. 74 145
[2] Bennett C H Brassard G Crepeau C Jozsa R Peres A Wootters W K 1993 Phys. Rev. Lett. 70 1895
[3] Wang M Y Yan F L 2016 Quantum Inform. Process. 15 3383
[4] Bennett C H Wiesner S J 1992 Phys. Rev. Lett. 69 2881
[5] Deng F G Long G L Liu X S 2003 Phys. Rev. A 68 042317
[6] Long G L Liu X S 2002 Phys. Rev. A 65 032302
[7] Hu J Y Yu B Jing M Y Xiao L T Jia S T Qin G Q Long G L 2016 Light Sci. Appl. 5 16144
[8] Zou X F Qiu D W 2014 Sci. China-Phys. Mech. Astron. 57 1696
[9] Zheng C Long G F 2014 Sci. China-Phys. Mech. Astron. 57 1238
[10] Briegel H J Dür W Cirac J I Zoller P 1998 Phys. Rev. Lett. 81 5932
[11] Li T Deng F G 2015 Sci. Rep. 5 15610
[12] Ekert A K 1991 Phys. Rev. Lett. 67 661
[13] Cao D Y Liu B H Wang Z Huang Y F Li C F Guo G C 2015 Sci. Bull 60 1128
[14] Chang Y Zhang S B Yan L L Han G H 2015 Chin. Phys. B 24 080306
[15] Ji Y Q Jin Z Zhu A D Wang H F Zhang S 2014 Chin. Phys. B 23 050306
[16] Wu H Z Yang Z B 2014 Chin. Phys. Lett. 31 024206
[17] Tang S Q Yuan J B Wang X W Kuang L M 2015 Chin. Phys. Lett. 32 040303
[18] Ai Q 2016 Sci. Bull. 61 110
[19] Wu Y L Li S J Ge W Xu Z X Tian L Wang H 2016 Sci. Bull. 61 302
[20] He Y Q Ding D Yan F L Gao T 2015 J. Phys. B: At. Mol. Opt. Phys. 48 055501
[21] He Y Q Ding D Yan F L Gao T 2016 Sci. Rep. 6 19116
[22] He Y Q Ding D Yan F L Gao T 2015 Opt. Express 23 21671
[23] Su X L Jia X J Xie C D Peng K C 2014 Sci. China-Phys. Mech. Astron. 57 1210
[24] Ding D Yan F L Gao T 2014 Sci. China-Phys. Mech. Astron. 57 2098
[25] Duan L M Lukin M D Cirac J I Zoller P 2001 Nature 414 413
[26] Bennett C H Bernstein H J Popescu S Schumacher B 1996 Phys. Rev. A 53 2046
[27] Bose S Vedral V Knight P L 1999 Phys. Rev. A 60 194
[28] Shi B S Jiang Y K Guo G C 2000 Phys. Rev. A 62 054301
[29] Yamamoto T Koashi M Imoto N 2001 Phys. Rev. A 64 012304
[30] Zhao Z Pan J W Zhan M S 2001 Phys. Rev. A 64 014301
[31] Sheng Y B Deng F G Zhou H Y 2008 Phys. Rev. A 77 062325
[32] Sheng Y B Zhou L Zhao S M Zheng B Y 2012 Phys. Rev. A 85 012307
[33] Sheng Y B Zhou L Zhao S M 2012 Phys. Rev. A 85 042302
[34] Deng F G 2012 Phys. Rev. A 85 022311
[35] Ren B C Du F F Deng F G 2013 Phys. Rev. A 88 012302
[36] Sheng Y B Zhou L 2013 Entropy 15 1776
[37] Sheng Y B Liu J Zhao S Y Wang L Zhou L 2014 Chin. Phys. B 23 080305
[38] Wang G Y Li T Deng F G 2015 Quantum Inform. Process. 14 1305
[39] Cao C Ding H Li Y Wang T J Mi S C Zhang R Wang C 2015 Quantum Inform. Process. 14 1265
[40] Wang C Cao C He L Y Zhang C L 2014 Quantum Inform. Process. 13 1025
[41] Sheng Y B Zhou L 2013 Chin. Phys. B 22 110303
[42] Zhou L 2014 Chin. Phys. B 23 050308
[43] Wang C Zhang Y Jin G S 2011 Phys. Rev. A 84 032307
[44] Zhou L Sheng Y B Wang X F 2014 J. Opt. Soc. Am. B 31 503
[45] Liu J Zhao S Y Zhou L Sheng Y B 2014 Chin. Phys. B 23 020313
[46] Zhou L Sheng Y B 2016 Chin. Phys. B 25 020308
[47] Wang C Shen W W Mi S C Zhang Y Wang T J 2015 Sci. Bull. 60 2016
[48] Maimaiti W Li Z Chesi S Wang Y D 2015 Sci. China-Phys. Mech. Astron. 58 050309
[49] Sheng Y B Feng Z F Ou-Yang Y Qu C C Zhou L 2014 Chin. Phys. Lett. 31 050303
[50] Liang B B Hu S Cui W X An C S Xing Y Hu J S Sun G Q Jiang X X Wang H F 2014 Laser Phys. Lett. 11 115202
[51] Shukla C Banerjee A Pathak A 2015 Quantum Inform. Process. 14 2077
[52] Fan L L Xia Y Song J 2014 Quantum Inform. Process. 13 1967
[53] Qu C C Zhou L Sheng Y B 2015 Quantum Inform. Process. 14 4131
[54] Sheng Y B Pan J Guo R Zhou L Wang L 2015 Sci. China-Phys. Mech. Astron. 58 060301
[55] Pan J Zhou L Gu S P Wang X F Sheng Y B Wang Q 2016 Quantum Inform. Process. 15 1669
[56] Zhou L Wang X F Sheng Y B 2014 Int. J. Theor. Phys. 53 1752
[57] Wang M Y Yan F L Xu J Z 2016 J. Phys. B: At. Mol. Opt. Phys. 49 155502
[58] Cao C Wang T J Zhang R Wang C 2015 Laser Phys. Lett. 12 036001
[59] Cao C Wang T J Zhang R Wang C 2015 J. Opt. Soc. Am. B 32 1524
[60] Cao C Chen X Duan Y W Fan L Zhang R Wang T J Wang C 2016 Sci. China-Phys. Mech. Astron. 59 100315
[61] Du F F Deng F G 2015 Sci. China-Phys. Mech. Astron. 58 040303
[62] Li X H Ghose S 2014 Laser Phys. Lett. 11 125201
[63] Li X H Ghose S 2015 Phys. Rev. A 91 062302
[64] Cirac J I Zoller P Kimble H J Mabuchi H 1997 Phys. Rev. Lett. 78 3221
[65] Cirac J I Ekert A K Huelga S F Macchiavello C 1999 Phys. Rev. A 59 4249
[66] Turchette Q A Hood C J Lange W Mabuchi H Kimble H J 1995 Phys. Rev. Lett. 75 4710
[67] Brune M Hagley E Dreyer J Maitre X Maali A Wunderlich C Raimond J M Haroche S 1996 Phys. Rev. Lett. 77 4887
[68] Mattle K Weinfurter H Kwiat P G Zeilinger A 1996 Phys. Rev. Lett. 76 4656
[69] Xue P Xiao Y F 2006 Phys. Rev. Lett. 97 140501
[70] Duan L M Kimble H J 2004 Phys. Rev. Lett. 92 127902
[71] Xiao Y F Lin X M Gao J Yang Y Han Z F Guo G C 2004 Phys. Rev. A 70 042314
[72] Cho J Lee H W 2005 Phys. Rev. Lett. 95 160501
[73] Duan L M Wang B Kimble H J 2005 Phys. Rev. A 72 032333
[74] Lin X M Zhou Z W Ye M Y Xiao Y F Guo G C 2006 Phys. Rev. A 73 012323
[75] Deng Z J Zhang X L Wei H Gao K L Feng M 2007 Phys. Rev. A 76 044305
[76] Wei H Deng Z J Zhang X L Feng M 2007 Phys. Rev. A 76 054304
[77] Zhou L Sheng Y B 2014 Phys. Rev. A 90 024301
[78] Zhou L Sheng Y B 2015 Phys. Rev. A 92 042314
[79] Dayan B Parkins A S Aoki T Ostby E P Vahala K I Kimble H J 2008 Science 319 1062
[80] An J H Feng M Oh C H 2009 Phys. Rev. A 79 032303
[81] Peng Z H Zou J Liu X J Xiao Y J Kuang L M 2012 Phys. Rev. A 86 034305
[82] McCall S L Hahn E L 1967 Phys. Rev. Lett. 18 908
[83] Patel C K N Slusher R E 1967 Phys. Rev. Lett. 19 1019
[84] Courtens E 1968 Phys. Rev. Lett. 21 3
[85] Labeyrie G Miniatura C Kaiserm R 2001 Phys. Rev. A 64 033402
[86] Julsgaard B Kozhekin A Polzik E S 2001 Nature 413 400
[87] Bastos W P Cardoso W B Avelar A T de Almeida N G Baseia B 2012 Quantum Inform. Process. 11 1867
[88] Chen J J An J H Feng M Liu G 2010 J. Phys. B: At. Mol. Opt. Phys. 43 095505
[89] Heine D Rohringer W Fischer D Wilzbach M Raub T Loziczky S Yuan L X Groth S Hessmo B Schmiedmayer J 2010 New J. Phys 12 095005
[90] D'uria V Lee N Amri T Fabre C Laurat J 2011 Phys. Rev. Lett. 107 050504
[91] Henrich D Rehm L Dörner S Hofherr M Ilin K Semenov A Siegel M 2013 IEEE Trans. Appl. Supercond. 23 2200405
[92] Lita A E Miller A J Nam S W 2008 Opt. Express 16 3032