Number state filtered coherent state
Abstract
Number state filtering in coherent states leads to sub-Poissonian photon statistics. These states are more suitable for phase estimation when compared with the coherent states. Nonclassicality of these states is quantified in terms of the negativity of the Wigner function and the entanglement potential. Filtering of the vacuum from a coherent state is almost like the photon-addition. However, filtering makes the state more resilient against dissipation than photon-addition. Vacuum state filtered coherent states perform better than the photon-added coherent states for a two-way quantum key distribution protocol. A scheme to generate these states in multi-photon atom-field interaction is presented.
I Introduction
Dynamics of the classical electromagnetic field is described by the Maxwell equations. All the features describable within this formalism are classified as the classical aspects of the field. However, this description is inadequate as there are experimental observations such as the anti-bunching of light, quadrature squeezing, sub-Poissonian photon statistics, etc., which cannot be understood within the classical description Gerry and Knight (2005). These observations are explainable in the quantum description of the electromagnetic field. In this formulation, the most fundamental states of the electromagnetic field are the number states, which correspond to states of definite number of photons. It is of interest to know those quantum states of the field which have features similar to the classical field. It turns out that the coherent state, which is a superposition of all the number states with specific superposition coefficients, has that property Glauber (1963a); Sudarshan (1963). This specific choice of the superposition coefficients makes the coherent state very different from these number states. For instance, the coherent states are the eigenstates of the annihilation operator , the corresponding complex eigenvalue is the coherent state amplitude. These states form an over-complete basis set so that any state of the electromagnetic field is expressible as a continuous superposition. As a consequence, the expectation values of an operator-valued function of the annihilation operator is the statistical average of the corresponding function of the coherent state amplitude. The expression for the coherent state in the number state basis is
(1) |
Correlation functions of any order are factorizable for the coherent states Glauber (1963b). For large amplitude coherent state, the photon number fluctuation is much smaller than the mean number of photon. Hence large amplitude coherent states represent the classical field. These states saturate the uncertainties in the field quadratures to their minimum values allowed by quantum mechanics. This, in turn, implies that the phase space distribution of a coherent state is well localized. In the context of a harmonic oscillator, the coherent state wave packet oscillates without dispersion. An alternate definition identifies the coherent state as a displaced vacuum state Gerry and Knight (2005). Interestingly, one of the quasi-probability distributions, namely, the Glauber-Sudarshan P-function for the coherent state is a Dirac delta-function, which is a well defined probability density. For this reason these states are called quasi-classical states. Any other choice of coefficients makes the resultant state non-classical whose P-function is more singular than a delta function. This signals the nonclassical aspect of the quantum state Glauber (1963a); Sudarshan (1963).
Nonclassicality of quantum states is of importance in quantum optics Mandel (1986); Pathak and Ghatak (2018). The notion of nonclassicality is based on the quasi-probability distributions for the quantum states. Negative value of P-function or Wigner function is sufficient to identify a quantum state as nonclassical.
Many interesting nonclassical states have been defined by suitable modifications of the quasi-classical coherent states; for instance, photon-added coherent states Agarwal and Tara (1991), even/odd coherent states Dodonov et al. (1974), truncated coherent states Miranowicz et al. (1994); Sivakumar (2014), etc. The even and odd coherent states are defined as the symmetric and anti-symmetric superposition of two coherent states whose complex amplitudes are of equal magnitude but out of phase by Dodonov et al. (1974). These states exhibit quadrature squeezing and sub-Poissonian statistics for suitable values of the amplitude. Agarwal and Tara introduced a new class of nonclassical states by adding quanta to the coherent state. This states are called photon added coherent state (PACS) Agarwal and Tara (1991). For proper choice of parameters, PACS shows quadrature squeezing and sub-Poissonian photon statistics. Truncated coherent states are obtained by modifying the expression in Eqn. 1 by removing either a finite number of contiguous states starting from the vacuum or by removing all the states beyond a specified number state Miranowicz et al. (1994); Sivakumar (2014). Another important class of nonclassical states is obtained by number state hole burning or filtering in the coherent states Baseia et al. (1998). It is known that any pure state in which the vacuum state is absent is nonclassical Lee (1995); de Lima and Baseia (1996).
A possible generalization is to consider states in which a number state, not necessarily the vacuum, is absent. More generally, states with many of the number states not being present in the superposition have been studied to assess their possible use in practical applications such as optical data storage and optical communication Baseia et al. (1998); Baseia and Malbouisson (2001). It may be noted that schemes to generate such states are knownSanaka et al. (2006); Malbouisson and Baseia (2001); Avelar and Baseia (2004, 2005); Gerry and Benmoussa (2002); Zagury and de Toledo Piza (1994); Ghosh and Gerry (1997); Escher et al. (2004); Cooper et al. (2015); Zou et al. (2016).
In this paper, we study the non classical properties of number state filtered coherent state. The photon counting statistics of these states is discussed in Section II. Sub-Poissonian statistics is shown to improve phase estimation in comparison to the coherent states. In Section III, quantification of nonclassicality of these states is carried out. Entanglement potential and negativity of Wigner function are studied in this context. In Section IV, the possibility of generating these states in the interaction of a 3-level atom with a cavity field is discussed. In Section V, features of vacuum state filtered coherent state are analyzed. Specifically, the robustness of the state against dissipation is exhibited. This feature enables these states to perform better than the single photon added coherent state (SPACS) when it comes to two-way quantum key distribution. Use of these states as a source to generate the even and odd coherent states is discussed in Section VI.
Ii Number state filtering in coherent state
Number state filtered state (NSFS), denoted by , is defined as
(2) |
where is the normalization constant given by
(3) |
This definition implies that the state is obtained if the number state is absent in the superposition defined in Eqn. 1.
The state does not become for any . However, if or vice-versa, then . It is worth noting that the coherent state is a superposition of NSFSs,
(4) |
where is odd and is even.
NSFS is expressible as
(5) |
a superposition of and , where Baseia et al. (1997). The overlap between the states and is
(6) |
which is minimum if . Therefore, maximum nonclassicality of is expected if .
For a given , the state corresponding to different are linearly independent. Hence the set can be used as a basis. Further, for a given
(7) |
which is the completeness relation for the Hilbert space less the number state of a harmonic oscillator.
It is interesting to note that the even and odd coherent states are expressible as
(8) |
and
(9) |
ii.1 Photon Statistics
Photon statistics is an experimentally measurable aspect of the electromagnetic field. Probability of detecting photons in a measurement is , where represents the number state with photons. Using the number state expression given in Eqn. 2,
(10) |
which gives the average number of photons in the state as
(11) |
For or
(12) |
as expected since under these limits. Surprisingly, if ,
(13) |
same as the average number of photons in the coherent state even though overlap between and is minimum. Nevertheless, the nonclassical nature of the photon distribution can be identified by the Mandel Q-parameter
(14) |
If , the state is Poissonian which corresponds to a coherent state. If , the state is sub-Poissonian which is a nonclassical feature. The state is super-Poissonian if . For the state ,
(15) |
where is given in Eqn. 11. The dependence of on is shown in Fig. 1 for and Fig. 1 for . The state is characterized by as well as in different ranges of .
To understand the photon statistics of , consider expressing it as a linear superposition of a lower truncated coherent state (LTCS) and an upper truncated coherent state (UTCS) whose number states expressions are Miranowicz et al. (1994); Sivakumar (2014)
(16) |
and
(17) |
respectively. Here and with the incomplete Gamma function. Note that, the states UTCS and LTCS are orthogonal to each other.
An interesting aspect of these two states is that they are sub-Poissonian for any . As increases to values much larger than unity, the photon number distribution of LTCS becomes Poissonian, that is, tends to zero. In the same limit, the parameter of UTCS becomes -1, that is, sub-Poissonian.
The state is expressible as a linear combination of LTCS and UTCS:
(18) |
where and .
Though the states involved in the superposition are sub-Poissonian, the resultant state is super-Poissonian for a range of values of as shown in Fig. 1. A key observation is that the state exhibits maximum super-Poissonian character when UTCS and LTCS are equally sub-Poissonian. As changes, exhibits sub-Poissonian feature over a small region and becomes nearly Poissonian for larger and smaller values of . This feature emerges here in a way that is analogous to the appearance of super-Poissonian statistics when two number states, which are sub-Poissonian, are superposed. If the superposition coefficients are nearly equal in magnitude, measurement of photon number gives results corresponding to the two number states with nearly equal probability. That is, the measurement results fluctuate and the statistics is super-Poissonian. To relate this with the behaviour of parameter of , and are shown as a function of for different values of in Fig. 2. When one of the coefficients is nearly unity, has negative corresponding to either LTCS or UTCS both of which are sub-Poissonian. When the coefficients are nearly equal, the statistics exhibits fluctuations as in the case of superposition of two number states and the resultant state is super-Poissonian.
The fact that the states exhibit sub-Poissonian character indicates that they could be of use in phase estimation using interferometers. Coherent states provide a phase resolution of Gerry and Knight (2005). There is mild improvement if one uses instead of the coherent states. This feature is shown in Fig. 3, where the minimum detectable phase is shown as a function of . It is clear that sensitivity is maximum when exhibits minimum .
ii.2 Subtraction of photons from NSFS
Photon-addition to a coherent state generates a nonclassical state Agarwal and Tara (1991). Like photon addition, photon subtraction is a process which may introduce nonclassicality in many states. However, photon subtraction leaves the coherent state unaltered as it is an eigenstate of . Action of on other states generates photon-subtracted states Takahashi et al. (2008). One way of realizing photon subtraction experimentally is by using beam splitters. Though is necessarily nonclassical, the effect of photon subtraction is of interest. The resultant state on subtracting a photon from is
(19) |
in which the number state is filtered. Successive action of leads to a sequence of NSFS wherein the successively lower number states are filtered. After photon subtractions, becomes
(20) |
in which the vacuum state is absent. This state is similar to the SPACS Agarwal and Tara (1991). The overlap between these two states is
(21) |
which is plotted as a function of in Fig. 4. It clearly indicates that for small or large , the states are nearly the same as the overlap is nearly unity.
Interestingly, subtraction of one more photon from yields
(22) |
the coherent state. In short, subtraction of photons from generates the coherent state. This has to be compared with the action of on a number state . After successive photon subtractions, becomes the vacuum state and one more subtraction annihilates the state.
It is to be noted that while subtraction of photons from generates a coherent state, addition of photons to a coherent state does not generate , a consequence of the non-commutativity of and .
In -representation, the wavefunction for a coherent state is a Gaussian. Probability density for in is , which is shown in Fig. 5. For comparison, the corresponding probability densities for coherent state of amplitude and the number state with is also shown. From the figure it is clear that, the probability distribution is not a Gaussian. In fact, the distribution is oscillatory.
Considering the time evolution of under , the evolved state is
(24) |
analogous to the evolution of the coherent state .
On filtering the number state from the coherent state , the resultant state has maximum deviation from if , which is the condition for maximum nonclassicality of the NSFS. On the other hand, NSFS exhibits maximum super-Poissonian photon statistics under the same condition. In fact, a state with super-Poissonian photon statistics does not refer to the nonclassicality of the state. In order to clarify the condition for maximum nonclassicality of NSFS, we study the negativity of Wigner function and entanglement potential in the next section.
Iii Nonclassicality of NSFS
Phase space representation of quantum states provides additional insights. For instance, classical aspects of quantum states are readily obtained by such representations. In quantum optics, the most used phase space distributions are the -, - and Wigner-functions. Negative values of Glauber-Sudarshan -function is a sufficient condition for qualifying the state as non-classical. The condition requires the -function to be negative. Hence, the amplitude squeezing is a non-classical effect Gerry and Knight (2005). A state with amplitude squeezing possesses sub-Poissonian photon statistics. The sub-Poissonian character of has been discussed in the Section II.
iii.1 Wigner function
In this section, we discuss the nonclassicality of in terms of its Wigner function. Wigner function, though not singular for any state, can become negative in some regions. The negativity is an indicator of quantumness Wigner (1932). Wigner function for a state is,
(25) |
which is an expression in terms of the displaced number state Gerry and Knight (2005); Bishop and Vourdas (1994).
For the state ,
(26) |
Here
(27) | ||||
(28) |
where is the associated Laguerre polynomial.
Wigner functions for with and are shown in Fig. 6. Both the states are nonclassical as their respective Wigner functions are negative in some regions of phase space. Interestingly, the number of peaks equals . To quantify the nonclassicality, the volume under negative portion of the Wigner function is shown in Fig. 7. It is readily inferred from the figure that the negativity is maximum when . This is consistent with the result that the overlap of with the classical state reaches its minimum when .
iii.2 Beam splitter transformation
Another way to quantify nonclassicality is the entanglement at output of a beam splitter. If one of the input states is a coherent state, it requires a nonclassical state at the second port to generate entanglement at the output Kim et al. (2002). The amount of entanglement is a measure of the nonclassicality of the input state to the second port. The unitary transformation effected by a beam splitter is
(29) |
The input bipartite state transforms to
(30) |
where and . Here is a coherent state.
Entanglement in the output of the beam splitter is quantified by the linear entropy of one of the reduced density operators of the bipartite output. It is defined as Wei et al. (2003),
(31) |
where . Here and are refer to the two output ports. For pure bipartite states, linear entropy is a good measure of entanglement. Non-zero values of imply entanglement in . Fig. 8 shows the entanglement in the output state for various values of with . It is clear from the figure that the state at the output ports of the beam splitter is entangled, which indicates that the input state is nonclassical. For large , and hence the entanglement of the output state is zero. Beam splitter generates maximally entangled state if for a given and independent of . This is consistent with the observation made earlier in Section III that the nonclassicality of is expected to be maximum if .
Iv Generation of NSFS
In this section, a scheme to generate is presented. This process involves driving a system of an atom and a cavity by two consecutive Hamiltonian evolutions. First, a multi-photon interaction drives the atom-cavity system. The resultant state is modified by driving the atomic system with a classical field. Details of the processes are indicated here.
Consider a type 3-level atom interacting with a cavity field of resonance frequency . The states of the 3-level atom are labeled by and in decreasing order of their energies as shown in Fig. 9. Interaction between the field and the atom is given by the multi-photon Hamiltonian Zubairy and Yeh (1980)
(32) |
which connects the states and . The transition from to is achieved by annihilating photons of the cavity field and the transition from to creates photons of the cavity field.
The initial state of the atom-cavity is , where the cavity field is in its vacuum and the atom is in the state . The time evolved state at time is
(33) |
The atomic transition comes in to the picture when the atom is driven by an external laser with Rabi frequency . The corresponding Hamiltonian is
(34) |
where , and . The transition frequency between the states and is . We assume . In the limit of large detuning () and in the rotating frame defined by , the effective Hamiltonian corresponding to is Gerry (2002)
(35) |
where and . If the state of the atom is , then the cavity field is unaffected. If the atomic state is or , the cavity evolution is non-trivial. Cavity driven is conditional depending on the atomic state.
Under the evolution governed by , evolves to
(36) |
where is a coherent state. Using a Ramsey pulse to drive the transitions and , the state transforms to
(37) |
where , and . The resultant state contains the required field state as a component. It is in post-selection of the atomic state, the cavity field is projected onto the required state. If the atom is detected in the state , the state of the field collapses to which is the target state for proper choices of and .
V NSFS for quantum information processing
The elementary unit for quantum computation and information processing is a qubit. Any two orthogonal quantum states of a system form a qubit. In the context of cavities, the vacuum state and the single photon state form a qubit. Various quantum protocols such as quantum teleportation Lombardi et al. (2002); Song et al. (2002), quantum state transfer Meher et al. (2017), entanglement generationTan et al. (1991), etc. have been implemented in this system. SPACS has also been suggested for implementing various quantum information protocols Pinheiro and Ramos (2013); Wang et al. (2014). Performance of any quantum information processing system is affected by external factors such as dissipation, decoherence, etc. It has always been of interest to use qubits which are resilient against such processes.
One way of assessing the robustness of the qubit evolving under a dissipative channel, is to calculate the overlap between its initial and final states. We consider four families of qubits using cavity field. One of the states of the qubit is the vacuum state. The other state can be a single photon state , a coherent state , a SPACS or a NSFS . Of these choices, the one with is not ideal qubit as the states and are not orthogonal. A phenomenological model of damping in a quantum channel is considered to study the effect of dissipation on these families of qubits Gerry and Knight (2005). The fidelity, which is the overlap of evolved state and the initial state, is
(38) |
where . The amplitude damping rate is . The respective fidelities for the states mentioned above are
(39) |
where and . The suffix indicates the state whose fidelity is being studied.
Fidelities for these various quantum states are shown in Fig. 10, assuming . It is to be noted that the state is robust against damping in a quantum channel as its fidelity decays slower than the other states. This makes the state suitable for quantum information processing.
Secure communication protocols such as BB84 are being carried out using the coherent states. The drawback in using the coherent states is that the eavesdropper can extract some photons, which may lead to loss of information. In order to prevent photon loss, Mundarain et. al have shown that SPACS can be used for two-way quantum key distribution Miranda and Mundarain (2017). Their protocol guarantees secure key distribution between two parties A and B even in the presence of a third party E who eavesdrops with a beam splitter. In order to quantify the information gathered by E, the respective probabilities of detecting at least one photon by B and E are calculated. If and SPACS are used, then these probabilities are
(40) |
and
(41) |
where is an angular parameter that characterize the beam splitter. It is clear that for small and , information loss is less.
The state is more robust against dissipative losses in a quantum channel. This may endow as a better choice than SPACS. The respective probabilities of detecting at least one photon by B and E when and are the input to the beam splitter are
(42) |
and
(43) |
In order to compare the advantage of over SPACS, the ratios of detection probabilities are considered. These are defined as
(44) |
and
(45) |
These two ratios are shown in Fig. 11. The ratio is always less than indicates that the diversion of photons by eavesdropper E via beam splitter transformation is less probable if one uses the state instead of SPACS. This shows performs better than SPACS.
Vi Generation of Schrdinger cat state from NSFS
Generation of the even and odd coherent states from is discussed. It is possible to realize the aforementioned states, for instance, from coherent states in cavities Brune et al. (1996), mixing squeezed vacuum and vacuum in beam splitter Wolinsky and Carmichael (1988); Dakna et al. (1997); Glancy and de Vasconcelos (2008), etc. The procedure for the generation of the odd/even coherent states from is similar to that for the generation of these states from the coherent state Brune et al. (1996). The system consists of a type three-level atom as described in Section IV and a cavity with resonance frequency . Consider the initial state
(46) |
A dispersive atom-field coupling is assumed so that the effective Hamiltonian is
(47) |
Time evolved state under is
(48) |
A Ramsey pulse is applied that transforms and . The evolved state at time is
(49) |
It is easy to recognize that if is even and the atom is detected in the state , the field state collapses to the odd coherent state. If is odd and the atom is detected in , the field state is projected on to the even coherent state.
Vii Summary
Coherent states are special as they exhibit minimum uncertainty in the quadratures, Poissonian statistics, factorizability of coherence functions, etc. All these features which make these states the most classical among the quantum states arise due to the specific choice of superposition coefficients. Number state filtering from a coherent state leads to a nonclassical state. These states are realizable in the multi-photon interaction between a three-level atom and a single mode cavity field in a coherent state. The number state that has the maximum weightage in a coherent state corresponds to . This is the condition for maximal nonclassicality for a number state filtered coherent state. The overlap between the coherent state and NSFS is minimum when this condition is satisfied. The negativity of the Wigner function is maximum under these conditions. Also, the entanglement potential of these states is highest when . If deviates from , the resultant state is still nonclassical as shown by the sub-Poissonian. This property enables these states to perform slightly better than the coherent states if used for phase measurements in an interferometric setup. Surprisingly, NSFS is super-Poissonian when . Emergence of this super-Poissonian statistics is understood based on the fact that NSFS is a superposition of suitably truncated coherent states which are sub-Poissonian. The super-Poissonian statistics emerges in an analogous manner as in suitable superposition of two number states which are also sub-Poissonian.
Further, vacuum state filtered coherent states are resilient against the effects of dissipation than the photon-added coherent states. This makes NSFS more suitable for quantum information processing. For instance, a robust qubit can be realized using the vacuum and NSFS as the computational basis states. The vacuum state filtered coherent state performs better than the single photon-added coherent state in two-way quantum key distribution where eavesdropping is carried out using a beam splitter.
Acknowledgement
One of the authors (NM) acknowledges the research fellowship from the Department of Atomic Energy, Government of India.
Footnotes
- preprint: APS/123-QED
References
- C. Gerry and P. Knight, Introductory Quantum Optics, edited by Cambridge (Cambridge (England): Cambridge UP, 2005.).
- R. J. Glauber, Phys. Rev. 131, 2766 (1963a).
- E. C. G. Sudarshan, Phys. Rev. Lett. 10, 277 (1963).
- R. J. Glauber, Phys. Rev. 130, 2529 (1963b).
- L. Mandel, Physica Scripta 1986, 34 (1986).
- A. Pathak and A. Ghatak, Journal of Electromagnetic Waves and Applications 32, 229 (2018), https://doi.org/10.1080/09205071.2017.1398109 .
- G. S. Agarwal and K. Tara, Phys. Rev. A 43, 492 (1991).
- V. Dodonov, I. Malkin, and V. Man’ko, Physica 72, 597 (1974).
- A. Miranowicz, K. Pia¸tek, and R. Tanaś, Phys. Rev. A 50, 3423 (1994).
- S. Sivakumar, International Journal of Theoretical Physics 53, 1697 (2014).
- B. Baseia, M. Moussa, and V. Bagnato, Physics Letters A 240, 277 (1998).
- C. T. Lee, Phys. Rev. A 52, 3374 (1995).
- A. F. de Lima and B. Baseia, Phys. Rev. A 54, 4589 (1996).
- B. Baseia and J. M. C. Malbouisson, Chinese Physics Letters 18, 1467 (2001).
- K. Sanaka, K. J. Resch, and A. Zeilinger, Phys. Rev. Lett. 96, 083601 (2006).
- J. Malbouisson and B. Baseia, Physics Letters A 290, 234 (2001).
- A. Avelar and B. Baseia, Optics Communications 239, 281 (2004).
- A. T. Avelar and B. Baseia, Phys. Rev. A 72, 025801 (2005).
- C. Gerry and A. Benmoussa, Physics Letters A 303, 30 (2002).
- N. Zagury and A. F. R. de Toledo Piza, Phys. Rev. A 50, 2908 (1994).
- H. Ghosh and C. C. Gerry, J. Opt. Soc. Am. B 14, 2782 (1997).
- B. M. Escher, A. T. Avelar, T. M. da Rocha Filho, and B. Baseia, Phys. Rev. A 70, 025801 (2004).
- M. Cooper, E. Slade, M. KarpiÅski, and B. J. Smith, New Journal of Physics 17, 033041 (2015).
- C.-L. Zou, L. Jiang, X.-B. Zou, and G.-C. Guo, Phys. Rev. A 94, 013841 (2016).
- B. Baseia, S. C. G. Granja, and G. C. Marques, Physica Scripta 55, 719 (1997).
- H. Takahashi, K. Wakui, S. Suzuki, M. Takeoka, K. Hayasaka, A. Furusawa, and M. Sasaki, Phys. Rev. Lett. 101, 233605 (2008).
- E. Wigner, Phys. Rev. 40, 749 (1932).
- R. F. Bishop and A. Vourdas, Phys. Rev. A 50, 4488 (1994).
- M. S. Kim, W. Son, V. Bužek, and P. L. Knight, Phys. Rev. A 65, 032323 (2002).
- T.-C. Wei, K. Nemoto, P. M. Goldbart, P. G. Kwiat, W. J. Munro, and F. Verstraete, Phys. Rev. A 67, 022110 (2003).
- M. S. Zubairy and J. J. Yeh, Phys. Rev. A 21, 1624 (1980).
- C. C. Gerry, Phys. Rev. A 65, 063801 (2002).
- E. Lombardi, F. Sciarrino, S. Popescu, and F. De Martini, Phys. Rev. Lett. 88, 070402 (2002).
- K.-H. Song, D.-H. Lu, W.-J. Zhang, and G.-C. Guo, Optics Communications 207, 219 (2002).
- N. Meher, S. Sivakumar, and P. K. Panigrahi, Scientific Reports 7, 9251 (2017).
- S. M. Tan, D. F. Walls, and M. J. Collett, Phys. Rev. Lett. 66, 252 (1991).
- P. V. P. Pinheiro and R. V. Ramos, Quantum Information Processing 12, 537 (2013).
- D. Wang, M. Li, F. Zhu, Z.-Q. Yin, W. Chen, Z.-F. Han, G.-C. Guo, and Q. Wang, Phys. Rev. A 90, 062315 (2014).
- M. Miranda and D. Mundarain, Quantum Information Processing 16, 298 (2017).
- M. Brune, E. Hagley, J. Dreyer, X. Maître, A. Maali, C. Wunderlich, J. M. Raimond, and S. Haroche, Phys. Rev. Lett. 77, 4887 (1996).
- M. Wolinsky and H. J. Carmichael, Phys. Rev. Lett. 60, 1836 (1988).
- M. Dakna, T. Anhut, T. Opatrný, L. Knöll, and D.-G. Welsch, Phys. Rev. A 55, 3184 (1997).
- S. Glancy and H. M. de Vasconcelos, J. Opt. Soc. Am. B 25, 712 (2008).